No Arabic abstract
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.
Recurrent novae (RNe) are cataclysmic variables with two or more nova eruptions within a century. Classical novae (CNe) are similar systems with only one such eruption. Many of the so-called CNe are actually RNe for which only one eruption has been discovered. Since RNe are candidate Type Ia supernova progenitors, it is important to know whether there are enough in our galaxy to provide the supernova rate, and therefore to know how many RNe are masquerading as CNe. To quantify this, we collected all available information on the light curves and spectra of a Galactic, time-limited sample of 237 CNe and the 10 known RNe, as well as exhaustive discovery efficiency records. We recognize RNe as having (a) outburst amplitude smaller than 14.5 - 4.5 * log(t_3), (b) orbital period >0.6 days, (c) infrared colors of J-H > 0.7 mag and H-K > 0.1 mag, (d) FWHM of H-alpha > 2000 km/s, (e) high excitation lines, such as Fe X or He II near peak, (f) eruption light curves with a plateau, and (g) white dwarf mass greater than 1.2 M_solar. Using these criteria, we identify V1721 Aql, DE Cir, CP Cru, KT Eri, V838 Her, V2672 Oph, V4160 Sgr, V4643 Sgr, V4739 Sgr, and V477 Sct as strong RN candidates. We evaluate the RN fraction amongst the known CNe using three methods to get 24% +/- 4%, 12% +/- 3%, and 35% +/- 3%. With roughly a quarter of the 394 known Galactic novae actually being RNe, there should be approximately a hundred such systems masquerading as CNe.
Novae are some of the most commonly detected optical transients and have the potential to provide valuable information about binary evolution. Binary population synthesis codes have emerged as the most effective tool for modelling populations of binary systems, but such codes have traditionally employed greatly simplified nova physics, precluding detailed study. In this work, we implement a model treating H and He novae as individual events into the binary population synthesis code binaryc. This treatment of novae represents a significant improvement on the `averaging treatment currently employed in modern population synthesis codes. We discuss the evolutionary pathways leading to these phenomena and present nova event rates and distributions of several important physical parameters. Most novae are produced on massive white dwarfs, with approximately 70 and 55 per cent of nova events occurring on O/Ne white dwarfs for H and He novae respectively. Only 15 per cent of H-nova systems undergo a common-envelope phase, but these systems are responsible for the majority of H nova events. All He-accreting He-nova systems are considered post-common-envelope systems, and almost all will merge with their donor star in a gravitational-wave driven inspiral. We estimate the current annual rate of novae in M31 (Andromeda) to be approximately $41 pm 4$ for H novae, underpredicting the current observational estimate of $65^{+15}_{-16}$, and $0.14pm0.015$ for He novae. When varying common-envelope parameters, the H nova rate varies between 20 and 80 events per year.
Cataclysmic Variables (CVs) and Symbiotic Binaries are close (or not so close) binary star systems which contain both a white dwarf (WD) primary and a larger cooler secondary star that typically fills its Roche Lobe. The cooler star is losing mass through the inner Lagrangian point of the binary and a fraction of this material is accreted by the WD. Here we report on our hydrodynamic studies of the thermonuclear runaway (TNR) in the accreted material that ends in a Classical Nova explosion. We have followed the evolution of the TNRs on both carbon-oxygen (CO) and oxygen-neon (ONe) WDs. We report on 3 studies in this paper. First, simulations in which we accrete only solar matter using NOVA (our 1-D, fully implicit, hydro code). Second, we use MESA for similar studies in which we accrete only Solar matter and compare the results. Third, we accrete solar matter until the TNR is ongoing and then switch the composition in the accreted layers to a mixed composition: either 25% WD and 75% solar or 50% WD and 50% Solar matter. We find that the amount of accreted material is inversely proportional to the initial 12C abundance (as expected). Thus, accreting solar matter results in a larger amount of accreted material to fuel the outburst; much larger than in earlier studies where a mixed composition was assumed from the beginning of the simulation. Our most important result is that all these simulations eject significantly less mass than accreted and, therefore, the WD is growing in mass toward the Chandrasekhar Limit.
Recent observations of a large number of DA and DB white dwarfs show evidence of debris disks, which are the remnants of old planetary systems. The infrared excess detected with emph{Spitzer} and the lines of heavy elements observed in their atmospheres with high-resolution spectroscopy converge on the idea that planetary material accretes onto these stars. Accretion rates have been derived by several authors with the assumption of a steady state between accretion and gravitational settling. The results are unrealistically different for DA and DB white dwarfs. When heavy matter is accreted onto stars, it induces an inverse $mu$-gradient that leads to fingering (thermohaline) convection. The aim of this letter is to study the impact of this specific process on the derived accretion rates in white dwarfs and on the difference between DA and DB. We solve the diffusion equation for the accreted heavy elements with a time-dependent method. The models we use have been obtained both with the IRAP code, which computes static models, and the La Plata code, which computes evolutionary sequences. Computations with pure gravitational settling are compared with computations that include fingering convection. The most important result is that fingering convection has very important effects on DAs but is inefficient in DBs. When only gravitational settling is taken into account, the time-dependent computations lead to a steady state, as postulated by previous authors. When fingering convection is added, this steady state occurs much later. The surprising difference found in the past for the accretion rates derived for DA and DB white dwarfs disappears. The derived accretion rates for DAs are increased when fingering convection is taken into account, whereas those for DBs are not modified. More precise and developed results will be given in a forthcoming paper.
We have examined the optical/X-ray light curves of seven well-observed recurrent novae, V745 Sco, M31N 2008-12a, LMC N 1968, U Sco, RS Oph, LMC N 2009a, T Pyx, and one recurrent nova candidate LMC N 2012a. Six novae out of the eight show a simple relation that the duration of supersoft X-ray source (SSS) phase is 0.70 times the total duration of the outburst ($=$ X-ray turnoff time), i.e., $t_{rm SSS}=0.70 t_{rm off}$, the total duration of which ranges from 10 days to 260 days. These six recurrent novae show a broad rectangular X-ray light curve shape, first half a period of which is highly variable in the X-ray count rate. The SSS phase corresponds also to an optical plateau phase that indicates a large accretion disk irradiated by a hydrogen-burning WD. The other two recurrent novae, T Pyx and V745 Sco, show a narrow triangular shape of X-ray light curve without an optical plateau phase. Their relations between $t_{rm SSS}$ and $t_{rm off}$ are rather different from the above six recurrent novae. We also present theoretical SSS durations for recurrent novae with various WD masses and stellar metallicities ($Z=$0.004, 0.01, 0.02, and 0.05) and compare with observed durations of these recurrent novae. We show that the SSS duration is a good indicator of the WD mass in the recurrent novae with a broad rectangular X-ray light curve shape.