The hydrogen-rich envelopes accreted by white dwarf stars from their red dwarf companions lead to thermonuclear runaways observed as classical nova eruptions peaking at up to 1 Million solar luminosities. Virtually all nova progenitors are novalike binaries exhibiting high rates of mass transfer to their white dwarfs before and after an eruption. It is a puzzle that binaries indistinguishable from novalikes, but with much lower mass transfer rates, and resulting dwarf nova outbursts, co-exist at the same orbital periods. Nova shells surrounding several dwarf novae demonstrate that at least some novae become dwarf novae between successive nova eruptions, though the mechanisms and timescales governing mass transfer rate variations are poorly understood. Here we report simulations of the multiGyr evolution of novae which self-consistently model every eruption's thermonuclear runaway, mass and angular momentum losses, feedback due to irradiation and variable mass transfer, and orbital size and period changes. The simulations reproduce the observed wide range of mass transfer rates at a given orbital period, with large and cyclic changes in white dwarf-red dwarf binaries emerging on kyr to Myr timescales. They also demonstrate that deep hibernation, (complete stoppage of mass transfer for long periods), occurs only in short-period binaries; that initially very different binaries converge to become nearly identical systems; that while almost all prenovae should be novalike binaries, dwarf novae should also occasionally be observed to give rise to novae; and that the masses of white dwarfs decrease only slightly while their red dwarf companions are consumed.

The supersoft X-ray and optical transient ASASSN-16oh has been interpreted by Maccarone et al. (2019) as having being induced by an accretion event on a massive white dwarf, resembling a dwarf nova super-outburst. These authors argued that the supersoft X-ray spectrum had a different origin than in an atmosphere heated by shell nuclear burning, because no mass was ejected. We find instead that the event's timescale and other characteristics are typical of non-mass ejecting thermonuclear runaways, as already predicted by Shara et al. (1977) and the extensive grid of nova models by Yaron et al. (2005). We suggest that the low X-ray and bolometric luminosity in comparison to the predictions of the models of nuclear burning are due to an optically thick accretion disk, hiding most of the white dwarf surface. If this is the case, we calculated that the optical transient can be explained as a non-ejective thermonuclear event on a WD of $\simeq$1.1M$_\odot$ accreting at the rate of $\simeq3.5{-}5{\times}10^{-7}$M$_\odot$yr$^{-1}$. We make predictions that should prove whether the nature of the transient event was due to thermonuclear burning or to accretion; observational proof should be obtained in the next few years, because a new outburst should occur within $\simeq$10-15 years of the event.

Models have long predicted that the frequency-averaged masses of white dwarfs in Galactic classical novae are twice as large as those of field white dwarfs. Only a handful of dynamically well-determined nova white dwarf masses have been published, leaving the theoretical predictions poorly tested. The recurrence time distributions and mass accretion rate distributions of novae are even more poorly known. To address these deficiencies, we have combined our extensive simulations of nova eruptions with the Strope et al (2010) and Schaefer et al (2010) databases of outburst characteristics of Galactic classical and recurrent novae to determine the masses of 92 white dwarfs in novae. We find that the mean mass (frequency averaged mean mass) of 82 Galactic classical novae is 1.06 (1.13) Msun, while the mean mass of 10 recurrent novae is 1.31 Msun. These masses, and the observed nova outburst amplitude and decline time distributions allow us to determine the long-term mass accretion rate distribution of classical novae. Remarkably, that value is just 1.3 x 10^{-10} Msun/yr, which is an order of magnitude smaller than that of cataclysmic binaries in the decades before and after classical nova eruptions. This predicts that old novae become low mass transfer rate systems, and hence dwarf novae, for most of the time between nova eruptions. We determine the mass accretion rates of each of the 10 known Galactic RN, finding them to be in the range 10^{-7} - 10^{-8} $ Msun/yr. We are able to predict the recurrence time distribution of novae and compare it with the predictions of population synthesis models.

Classical novae show a rapid rise in optical brightness over a few hours. Until recently the rise phase, particularly the phenomenon of a pre-maximum halt, was observed sporadically. Solar observation satellites observing Coronal Mass Ejections enable us to observe the pre-maximum phase in unprecedented temporal resolution. We present observations of V5589 Sgr with STEREO HI-1B at a cadence of 40 min, the highest to date. We temporally resolve a pre-maximum halt for the first time, with two examples each rising over 40 min then declining within 80 min. Comparison with a grid of outburst models suggests this double peak, and the overall rise timescale, are consistent with a white dwarf mass, central temperature and accretion rate close to 1.0 solar mass, 5x10^7 K and 10^-10 solar masses per year respectively. The modelling formally predicts mass loss onset at JD 2456038.2391+/-0.0139, 12 hrs before optical maximum. The model assumes a main-sequence donor. Observational evidence is for a subgiant companion; meaning the accretion rate is under-estimated. Post-maximum we see erratic variations commonly associated with much slower novae. Estimating the decline rate difficult, but we place the time to decline two magnitudes as 2.1 < t_2(days) < 3.9 making V5589 Sgr a "very fast" nova. The brightest point defines "day 0" as JD 2456038.8224+/-0.0139, although at this high cadence the meaning of the observed maximum becomes difficult to define. We suggest that such erratic variability normally goes undetected in faster novae due to the low cadence of typical observations; implying erratic behaviour is not necessarily related to the rate of decline.

Can a white dwarf, accreting hydrogen-rich matter from a non-degenerate companion star, ever exceed the Chandrasekhar mass and explode as a type Ia supernova? We explore the range of accretion rates that allow a white dwarf (WD) to secularly grow in mass, and derive limits on the accretion rate and on the initial mass that will allow it to reach $1.4M_\odot$ --- the Chandrasekhar mass. We follow the evolution through a long series of hydrogen flashes, during which a thick helium shell accumulates. This determines the effective helium mass accretion rate for long-term, self-consistent evolutionary runs with helium flashes. We find that net mass accumulation always occurs despite helium flashes. Although the amount of mass lost during the first few helium shell flashes is a significant fraction of that accumulated prior to the flash, that fraction decreases with repeated helium shell flashes. Eventually no mass is ejected at all during subsequent flashes. This unexpected result occurs because of continual heating of the WD interior by the helium shell flashes near its surface. The effect of heating is to lower the electron degeneracy throughout the WD, and especially in the outer layers. This key result yields helium burning that is quasi-steady state, instead of explosive. We thus find a remarkably large parameter space within which long-term, self-consistent simulations show that a WD can grow in mass and reach the Chandrasekhar limit, despite its helium flashes.

A definitive determination of the progenitors of type Ia supernovae (SNIa) has been a conundrum for decades. The single degenerate scenario $-$ a white dwarf (WD) in a semi-detached binary system accreting mass from its secondary $-$ is a plausible path; however, no simulation to date has shown that such an outcome is possible. In this study, we allowed a WD with a near Chandrasekhar mass of $1.4M_\odot$ to evolve over tens of thousands of nova cycles, accumulating mass secularly while undergoing periodic nova eruptions. We present the mass accretion limits within which a SNIa can possibly occur. The results showed, for each parameter combination within the permitted limits, tens of thousands of virtually identical nova cycles where the accreted mass exceeded the ejected mass, i.e. the WD grew slowly but steadily in mass. Finally, the WD became unstable, the maximal temperature rose by nearly two orders of magnitude, heavy element production was enhanced by orders of magnitude and the nuclear and neutrino luminosities became enormous. We also found that this mechanism leading to WD collapse is robust, with WDs in the range $1.0 - 1.38M_\odot$, and an accretion rate of $5*10^{-7}M_\odot/yr$, all growing steadily in mass. These simulations of the onset of a SNIa event make observationally testable predictions about the light curves of pre-SN stars, and about the chemistry of SNIa ejecta.