We describe a new code to simulate the stellar evolution of a close interacting binary system. It is then used to calculate the evolution of a classical nova system composed of a 1.25 Msun Main-Sequence (MS) star and a 1.0 Msun white dwarf (WD) companion. The system begins as a well separated non-interacting binary system. Initially, the two stars evolve independently of each other. However, Roche lobe overflow begins as the MS star expands on its way to become a Red Giant. We follow the mass accreted onto the WD and the ensuing nuclear runaways for several thousand flashes. The main finding is that the Roche-Lobe mass transfer rate is modulated by oscillations in the MS star, with a period that is somewhat shorter than the thermal time scale of the star. This periodically modulates the rate of thermonuclear flashes on the WD, between once every 12000 yrs, such that the WD can cool, to once every 300 yrs, such that it cannot. The system is further complicated by the secular drift in the secondary modulation. Such secondary modulation could explain systems like T Pyxidis. Last, we find that the overall process of mass gain by the WD has an efficiency of roughly 9%, thus requiring a donor with an initial mass of larger than about 5 Msun MS for an initial 1 Msun WD, if the WD is to reach the Chandrasekhar mass.

We study C/O white dwarfs with masses of 1.0 to 1.4 Msun accreting solar-composition material at very high accretion rates. We address the secular changes in the WDs, and in particular, the question whether accretion and the thermonuclear runaways result is net accretion or erosion. The present calculation is unique in that it follows a large number of cycles, thus revealing the secular evolution of the WD system. We find that counter to previous studies, accretion does not give rise to steady state burning. Instead, it produces cyclic thermonuclear runaways of two types. During most of the evolution, many small cycles of hydrogen ignition and burning build a helium layer over the surface of the white dwarf. This He layer gradually thickens and progressively becomes more degenerate. Once a sufficient amount of He has accumulated, several very large helium burning flashes take place and expel the accreted envelope, leaving no net mass accumulation. The results imply that such a system will not undergo an accretion induced collapse, nor will it lead to a SN Type Ia, unless a major new physical process is found.

Radiation from accretion discs in cataclysmic variable stars (CVs) provides fundamental information about the properties of these close binary systems and about the physics of accretion in general. The detailed diagnostics of accretion disc structure can be achieved by including in its description all the relevant heating and cooling physical mechanism, in particular the convective energy transport that, although dominant at temperatures less than about 10 000 K, is usually not taken into account when calculating spectra of accretion discs. We constructed a radiative transfer code coupled with a code determining the disc's hydrostatic vertical structure. We have obtained for the first time model spectra of cold, convective accretion discs. As expected, these spectra are mostly flat in the optical wavelengths with no contribution from the UV, which in quiescence must be emitted by the white dwarf. The disc structures obtained with our radiative-transfer code compare well with the solutions of equations used to describe the dwarf-nova outburst cycle according to the thermal-viscous disc instability model thus allowing the two to be combined. Our code allows calculating the spectral evolution of dwarf nova stars through their whole outburst cycle, providing a new tool for testing models of accretion discs in cataclysmic variables. We show that convection plays an important role in determining the vertical disc structure and substantially affects emitted spectra when, as often the case, it is effective at optical depths tau ~ 1. The emergent spectrum is independent of the parameters of the convection model.(Abstract shortened)

Binary evolution predicts a population of helium core (M < 0.5 Msol) white dwarfs (WDs) that are slowly accreting hydrogen-rich material from low mass main sequence or brown dwarf donors with orbital periods less than four hours. Four binaries are presently known in the Milky Way that will reach such a mass-transferring state in a few Gyr. Despite these predictions and observations of progenitor binaries, there are still no secure cases of helium core WDs among the mass-transferring cataclysmic variables (CVs). This led us to calculate the fate of He WDs once accretion begins at a rate Mdot < 1e-10 Msol/yr set by angular momentum losses. We show here that the cold He core temperatures (T_c < 1e7 K) and low Mdot result in ~ 1e-3 Msol of accumulated H-rich material at the onset of the thermonuclear runaway. Shara and collaborators noted that these large accumulated masses may lead to exceptionally long classical nova (CN) events. For a typical donor star of 0.2 Msol, such binaries will only yield a few hundred CNe, making these events rare amongst all CNe. We calculate the reheating of the accreting WD, allowing a comparison to the measured WD effective temperatures in quiescent dwarf novae and raising the possibility that WD seismology may be the best way to confirm the presence of a He WD. We also find that a very long (> 1000 yr) stable burning phase occurs after the CN outburst, potentially explaining enigmatic short orbital period supersoft sources like RX J0537-7034 (P_orb = 3.5 hr) and 1E 0035.4-7230 (P_orb = 4.1 hr).

We present spectra of accretion discs around white-dwarfs calculated with an improved and updated version of the Shaviv & Wehrse (1991) model. The new version includes line opacities and convective energy transport and can be used to calculate spectra of hot discs in bright systems (nova--like variables or dwarf novae in outburst) as well as spectra of cold accretion discs in quiescent dwarf novae.

We use the Shaviv and Wehrse (1991) code to model the vertical structure and the emission properties of quiescent dwarf nova discs. We find that in the case of HT Cas the quiescent disc must be optically thin, in contradiction with the requirements of the standard disc instability model. We find a viscosity parameter alpha > 1. Although this is much less than values (~ 100) obtained in isothermal slab models it is not consistent with the accretion disc model assumptions.