New Insights into Classical Novae

We survey our understanding of classical novae—nonterminal, thermonuclear eruptions on the surfaces of white dwarfs in binary systems. The recent and unexpected discovery of GeV gamma rays from Galactic novae has highlighted the complexity of novae and their value as laboratories for studying shocks and particle acceleration. We review half a century of nova literature through this new lens, and conclude the following: ▪ The basics of the thermonuclear runaway theory of novae are confirmed by observations. The white dwarf sustains surface nuclear burning for some time after runaway, and until recently, it was commonly believed that radiation from this nuclear burning solely determines the nova's bolometric luminosity. ▪ The processes by which novae eject material from the binary system remain poorly understood. Mass loss from novae is complex (sometimes fluctuating in rate, velocity, and morphology) and often prolonged in time over weeks, months, or years. ▪ The complexity of the mass ejection leads to gamma-ray-producing shocks internal to the nova ejecta. When gamma rays are detected (around optical maximum), the shocks are deeply embedded and the surrounding gas is very dense. ▪ Observations of correlated optical and gamma-ray light curves confirm that the shocks are radiative and contribute significantly to the bolometric luminosity of novae. Novae are therefore the closest and most common interaction-powered transients. Expected final online publication date for the Annual Review of Astronomy and Astrophysics, Volume 59 is September 2021. Please see http://www.annualreviews.org/page/journal/pubdates for revised estimates.

Isotopic ratios in the red giant component of the recurrent nova T Coronae Borealis

ABSTRACT We report the determination of abundances and isotopic ratios for C, O, and Si in the photosphere of the red giant (RG) component of the recurrent nova (RN) T Coronae Borealis from new 2.284–2.402 μm and 3.985–4.155 μm spectroscopy. Abundances and isotopic ratios in the photosphere may be affected by (i) processes in the RG interior which are brought to the surface during dredge-up and (ii) contamination of the RG, either during the common envelope phase of the binary evolution or by material synthesized in RN eruptions, or a combination of the two. We find that the abundances of C, O, and Si are reasonably consistent with the expected composition of an RG after first dredge-up, as is the 16O/17O ratio. The 28Si/29Si ratio is found to be 8.6 ± 3.0, and that for 28Si/30Si is 21.5 ± 3.0. The 12C/13C ratio (10 ± 2) is somewhat lower than expected for first dredge-up. The 16O/18O ratio (41 ± 3) is highly inconsistent with that expected either from RG evolution (∼550) or from contamination of the RG by the products of a nova thermonuclear runaway. In particular, the C and O isotopic ratios taken in combination are a puzzle. We urge confirmation of our results using spectroscopy at high resolution. We also encourage a thorough theoretical study of the effects on the secondary star in an RN system of contamination by ejecta having anomalous abundances and isotopic ratios.

Type Ia supernovae from non-accreting progenitors

Type Ia supernovae (SNe Ia) are manifestations of stars that are deficient in hydrogen and helium, and disrupt in a thermonuclear runaway. While explosions of carbon-oxygen white dwarfs are thought to account for the majority of events, part of the observed diversity may be due to varied progenitor channels. We demonstrate that helium stars with masses between ∼1.8 and 2.5 M⊙ may evolve into highly degenerate cores with near-Chandrasekhar mass and helium-free envelopes that subsequently ignite carbon and oxygen explosively at densities of ∼(1.8−5.9) × 109 g cm−3. This occurs either due to core growth from shell burning (when the core has a hybrid CO/NeO composition), or following ignition of residual carbon triggered by exothermic electron captures on 24Mg (for a NeOMg-dominated composition). We argue that the resulting thermonuclear runaway is likely to prevent core collapse, leading to the complete disruption of the star. The available nuclear energy at the onset of explosive oxygen burning suffices to create ejecta with a kinetic energy of ∼1051 erg, as in typical SNe Ia. Conversely, if these runaways result in partial disruptions, the corresponding transients would resemble SN Iax events similar to SN 2002cx. If helium stars in this mass range indeed explode as SNe Ia, then the frequency of events would be comparable to the observed SN Ib/c rates, thereby sufficing to account for the majority of SNe Ia in star-forming galaxies.

Search for 7Be in the outbursts of four recent novae

ABSTRACT Following the recent detection of 7Be ii in the outburst spectra of classical novae, we report the search for this isotope in the outbursts of four recent bright novae by means of high-resolution Ultraviolet and Visual Echelle Spectrograph (UVES) observations. The 7Be ii λλ313.0583, 313.1228 nm doublet resonance lines are detected in the high-velocity components of Nova Mus 2018 and ASASSN-18fv during outbursts. However, 7Be ii is not detected in ASASSN-17hx and possibly not in Nova Cir 2018, which shows that 7Be is not always ejected in the thermonuclear runaway. Taking into account the 7Be decay, we find X(7Be)/X(H) ≈ 1.5 × 10−5 and 2.2 × 10−5 in Nova Mus 2018 and ASASSN-18fv, respectively. A value of 7Be/H ≈ 2 × 10−5 is found in five out of the seven extant measurements, and it can be considered as a typical 7Be yield for novae. However, this value is almost one order of magnitude larger than predicted by current theoretical models. We argue that the variety of high 7Be/H abundances could be the result of a higher than solar content of 3He in the donor star. The cases with 7Be not detected might be related to the small mass of the white dwarf (WD) or to relatively little mixing with the core material of the WD. The 7Be/H, or 7Li/H, abundance is ≈ 4 dex above meteoritic abundance, thus confirming the novae as the main sources of 7Li in the Milky Way.

123–321 models of classical novae

Context. High-resolution spectroscopy has revealed large concentrations of CNO and sometimes other intermediate-mass elements (e.g., Ne, Na, Mg, or Al, for ONe novae) in the shells ejected during nova outbursts, suggesting that the solar composition material transferred from the secondary mixes with the outermost layers of the underlying white dwarf during thermonuclear runaway. Aims. Multidimensional simulations have shown that Kelvin-Helmholtz instabilities provide self-enrichment of the accreted envelope with material from the outermost layers of the white dwarf, at levels that agree with observations. However, the Eulerian and time-explicit nature of most multidimensional codes used to date and the overwhelming computational load have limited their applicability, and no multidimensional simulation has been conducted for a full nova cycle. Methods. This paper explores a new methodology that combines 1D and 3D simulations. The early stages of the explosion (i.e., mass-accretion and initiation of the runaway) were computed with the 1D hydrodynamic code SHIVA. When convection extended throughout the entire envelope, the structures for each model were mapped into 3D Cartesian grids and were subsequently followed with the multidimensional code FLASH. Two key physical quantities were extracted from the 3D simulations and were subsequently implemented into SHIVA, which was used to complete the simulation through the late expansion and ejection stages: the time-dependent amount of mass dredged-up from the outer white dwarf layers, and the time-dependent convective velocity profile throughout the envelope. Results. This work explores for the first time the effect of the inverse energy cascade that characterizes turbulent convection in nova outbursts. More massive envelopes have been found that are those reported from previous models with pre-enrichment. These result in more violent outbursts, characterized by higher peak temperatures and greater ejected masses, with metallicity enhancements in agreement with observations.

Two-dimensional simulations of mixing in classical novae: The effect of white dwarf composition and mass

Context. Classical novae are explosive phenomena that take place in stellar binary systems. They are powered by mass transfer from a low-mass main sequence star onto either a CO or ONe white dwarf. The material accumulates for 104–105 yr until ignition under degenerate conditions, resulting in a thermonuclear runaway. The nuclear energy released produces peak temperatures of ∼0.1–0.4 GK. During these events, 10−7−10−3 M⊙ enriched in intermediate-mass elements, with respect to solar abundances, are ejected into the interstellar medium. However, the origin of the large metallicity enhancements and the inhomogeneous distribution of chemical species observed in high-resolution spectra of ejected nova shells is not fully understood. Aims. Recent multidimensional simulations have demonstrated that Kelvin-Helmholtz instabilities that operate at the core-envelope interface can naturally produce self-enrichment of the accreted envelope with material from the underlying white dwarf at levels that agree with observations. However, such multidimensional simulations have been performed for a small number of cases and much of the parameter space remains unexplored. Methods. We investigated the dredge-up, driven by Kelvin-Helmholtz instabilities, for white dwarf masses in the range 0.8–1.25 M⊙ and different core compositions, that is, CO-rich and ONe-rich substrates. We present a set of five numerical simulations performed in two dimensions aimed at analyzing the possible impact of the white dwarf mass, and composition, on the metallicity enhancement and explosion characteristics. Results. At the time we stop the simulations, we observe greater mixing (∼30% higher when measured in the same conditions) and more energetic outbursts for ONe-rich substrates than for CO-rich substrates and more massive white dwarfs.