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.

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.

Ten weeks of daily imaging of the giant elliptical galaxy M87 with the Hubble Space Telescope (HST) has yielded 41 nova light curves of unprecedented quality for extragalactic cataclysmic variables. We have recently used these light curves to demonstrate that the observational scatter in the so-called Maximum-Magnitude Rate of Decline (MMRD) relation for classical novae is so large as to render the nova-MMRD useless as a standard candle. Here we demonstrate that a modified Buscombe - de Vaucouleurs hypothesis, namely that novae with decline times t2 > 10 days converge to nearly the same absolute magnitude about two weeks after maximum light in a giant elliptical galaxy, is supported by our M87 nova data. For 13 novae with daily-sampled light curves, well determined times of maximum light in both the F606W and F814W filters, and decline times $t2 > 10 days we find that M87 novae display M(606W,15) = -6.37 +/- 0.46 and M(814W,15) = -6.11 +/- 0.43. If very fast novae with decline times t2 < 10 days are excluded, the distances to novae in elliptical galaxies with stellar binary populations similar to those of M87 should be determinable with 1 sigma accuracies of +/-20% with the above calibrations.

Horizontal branch (HB) morphology is a complex multiple-parameter problem. Besides the metallicity, two other leading parameters are the mass loss rate (MLR) and the initial He abundance of the HB progenitors. Using the STAREV stellar evolution code, we produce a wide array of Extreme Horizontal Branch (EHB) stars and also examine their post-HB evolution. EHB stars are produced in our calculations by the so called `delayed (late) hot core flash' scenario. The MLR is increased on the red giant branch (RGB) to the extent that, prior to reaching core flash conditions, only a very thin H-rich envelope remains and helium ignition takes place at hotter positions on the HRD. We perform an extensive, self-consistent parameter study, covering populations I and II (Z=0.0001-0.03), for both normal initial helium abundances and He-enriched models (up to Y=0.40). For completeness of the study and in order to point to complete trends, we chose NOT to cut out several combinations (or results of) that may extend beyond realistic limits. We present results and complete evolutionary tracks for the covered parameter space, showing in particular that: a) Increased He abundance ON ITS OWN -- without having a significant-enough MLR on the RGB -- DOES NOT lead to the production of EHB stars; however, b) The bluest (hottest) HB positions do result from the COMBINED effect of He-enhancement and increased MLR; c) The general trend is that the effective temperature on the HB increases with decreasing metallicity, but there is an indication for a halt, or even a reversal of this trend, as Z further decreases below $10^{-3}$; d) EHB stars can serve as major contributors to the UV flux emanating from their host system. Thus, the present comprehensive study both complements and lends support to previous, more restricted studies of the HB phase, and adds results for unexplored regions of the parameter space.

The extensive grid of numerical simulations of nova eruptions of Yaron et al.(2005) first predicted that some classical novae might deviate significantly from the Maximum Magnitude - Rate of Decline (MMRD) relation, which purports to characterise novae as standard candles. Kasliwal et al. (2011) have announced the observational detection of an apparently new class of faint, fast classical novae in the Andromeda galaxy. These objects deviate strongly from the MMRD relationship, exactly as predicted by Yaron et al. (2005). Shara et al. (2016) recently reported the first detections of faint, fast novae in M87. These previously overlooked objects are as common in the giant elliptical galaxy M87 as they are in the giant spiral M31; they comprise about 40% of all classical nova eruptions and greatly increase the observational scatter in the MMRD relation. We use the extensive grid of nova simulations of Yaron et al. (2005) to identify the underlying causes of the existence of faint, fast novae. These are systems which have accreted, and can thus eject, only very low mass envelopes, of order 10^-7 - 10^-8 Msun, on massive white dwarfs. Such binaries include, but are not limited to, the recurrent novae. These same models predict the existence of ultrafast novae which display decline times t2 as short as five hours. We outline a strategy for their future detection.

The internal structure of gas giant planets may be more complex than the commonly assumed core-envelope structure with an adiabatic temperature profile. Different primordial internal structures as well as various physical processes can lead to non-homogenous compositional distributions. A non-homogenous internal structure has a significant impact on the thermal evolution and final structure of the planets. In this paper, we present alternative structure and evolution models for Jupiter and Saturn allowing for non-adiabatic primordial structures and the mixing of heavy elements by convection as these planets evolve. We present the evolution of the planets accounting for various initial composition gradients, and in the case of Saturn, include the formation of a helium-rich region as a result of helium rain. We investigate the stability of regions with composition gradients against convection, and find that the helium shell in Saturn remains stable and does not mix with the rest of the envelope. In other cases, convection mixes the planetary interior despite the existence of compositional gradients, leading to the enrichment of the envelope with heavy elements. We show that non-adiabatic structures (and cooling histories) for both Jupiter and Saturn are feasible. The interior temperatures in that case are much higher that for standard adiabatic models. We conclude that the internal structure is directly linked to the formation and evolution history of the planet. These alternative internal structures of Jupiter and Saturn should be considered when interpreting the upcoming Juno and Cassini data.

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.

The primordial internal structures of gas giant planets are unknown. Often giant planets are modeled under the assumption that they are adiabatic, convective, and homogeneously mixed, but this is not necessarily correct. In this work, we present the first self-consistent calculation of convective transport of both heat and material as the planets evolve. We examine how planetary evolution depends on the initial composition and its distribution, whether the internal structure changes with time, and if so, how it affects the evolution. We consider various primordial distributions, different compositions, and different mixing efficiencies and follow the distribution of heavy elements in a Jupiter-mass planet as it evolves. We show that a heavy-element core cannot be eroded by convection if there is a sharp compositional change at the core-envelope boundary. If the heavy elements are initially distributed within the planet according to some compositional gradient, mixing occurs in the outer regions resulting in a compositionally homogeneous outer envelope. Mixing of heavy materials that are injected in a convective gaseous envelope are found to mix efficiently. Our work demonstrates that the primordial internal structure of a giant planet plays a substantial role in determining its long-term evolution and that giant planets can have non-adiabatic interiors. These results emphasize the importance of coupling formation, evolution, and internal structure models of giant planets self-consistently.

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.

We model the evolution of planets with various masses and compositions. We investigate the effects of the composition and its depth dependence on the long-term evolution of the planets. The effects of opacity and stellar irradiation are also considered. It is shown that the change in radius due to various compositions can be significantly smaller than the change in radius caused by the opacity. Irradiation also affects the planetary contraction but is found to be less important than the opacity effects. We suggest that the mass-radius relationship used for characterization of observed extrasolar planets should be taken with great caution since different physical conditions can result in very different mass-radius relationships.

Extremely luminous, red eruptive variables like RV in M31 are being suggested as exemplars of a new class of astrophysical object. Our greatly extended series of nova simulations shows that classical nova models can produce very red, luminous eruptions. In a poorly studied corner of 3-D nova parameter space (very cold, low-mass white dwarfs, accreting at very low rates) we find bona fide classical novae that are very luminous and red because they eject very slowly moving, massive envelopes. A crucial prediction of these nova models - in contrast to the predictions of merging star ("mergeburst") models - is that a hot remnant, the underlying white dwarf, will emerge after the massive ejected envelope has expanded enough to become optically thin. This blue remnant must fade on a timescale of decades - much faster than a "mergeburst", which must fade on timescales of millennia or longer. Furthermore, the cooling nova white dwarf and its expanding ejecta must become redder in the years after eruption, while a contracting mergeburst must become hotter and bluer. We predict that red novae will always brighten to L~1000 Lsun for about 1 year before rising to maximum luminosity at L~10^6 - 10^7 Lsun. The maximum luminosity attainable by a nova is likely to be L~10^7 Lsun, corresponding to M-12. In an accompanying paper we describe a fading, luminous blue candidate for the remnant of M31-RV; it is observed with HST to be compatible only with the nova model.

M31-RV was an extraordinarily luminous (~10^6 Lsun) eruptive variable, displaying very cool temperatures (roughly 1000 Kelvins) as it faded. The photometric behavior of M31-RV (and several other very red novae, i.e. luminous eruptive red variables) has led to several models of this apparently new class of astrophysical object. One of the most detailed models is that of "mergebursts": hypothetical mergers of close binary stars. These are predicted to rival or exceed the brightest classical novae in luminosity, but to be much cooler and redder than classical novae, and to become slowly hotter and bluer as they age. This prediction suggests two stringent and definitive tests of the mergeburst hypothesis. First, there should always be a cool red remnant, and NOT a hot blue remnant at the site of such an outburst. Second, the inflated envelope of a mergeburst event should be slowly contracting, hence it must display a slowly rising effective temperature. We have located a luminous, UV-bright object within 0.4 arcsec (1.5 sigma of the astrometric position) of M31-RV in archival WFPC2 images taken 10 years after the outburst: it resembles an old nova. Twenty years after the outburst, the object remains much too hot to be a mergeburst. Its behavior remains consistent with that of theoretical nova models which erupt on a low mass white dwarf. Future Hubble UV and visible images could determine if the M31-RV analogs (in M85 and in M99) are also behaving like old novae.

The total masses ejected during classical nova eruptions are needed to answer two questions with broad astrophysical implications: Can accreting white dwarfs be pushed towards the Chandrasekhar mass limit to yield type Ia supernovae? Are Ultra-luminous red variables a new kind of astrophysical phenomenon, or merely extreme classical novae? We review the methods used to determine nova ejecta masses. Except for the unique case of BT Mon (nova 1939), all nova ejecta mass determinations depend on untested assumptions and multi-parameter modeling. The remarkably simple assumption of equipartition between kinetic and radiated energy (E_kin and E_rad, respectively) in nova ejecta has been invoked as a way around this conundrum for the ultra-luminous red variable in M31. The deduced mass is far larger than that produced by any classical nova model. Our nova eruption simulations show that radiation and kinetic energy in nova ejecta are very far from being in energy equipartition, with variations of four orders of magnitude in the ratio E_kin/E_rad being commonplace. The assumption of equipartition must not be used to deduce nova ejecta masses; any such "determinations" can be overestimates by a factor of up to 10,000. We data-mined our extensive series of nova simulations to search for correlations that could yield nova ejecta masses. Remarkably, the mass ejected during a nova eruption is dependent only on (and is directly proportional to) E_rad. If we measure the distance to an erupting nova and its bolometric light curve then E_rad and hence the mass ejected can be directly measured.

We present a new stellar evolution code and a set of results, demonstrating its capability at calculating full evolutionary tracks for a wide range of masses and metallicities. The code is fast and efficient, and is capable of following through all evolutionary phases, without interruption or human intervention. It is meant to be used also in the context of modeling the evolution of dense stellar systems, for performing live calculations for both normal star models and merger-products. The code is based on a fully implicit, adaptive-grid numerical scheme that solves simultaneously for structure, mesh and chemical composition. Full details are given for the treatment of convection, equation of state, opacity, nuclear reactions and mass loss. Results of evolutionary calculations are shown for a solar model that matches the characteristics of the present sun to an accuracy of better than 1%; a 1 Msun model for a wide range of metallicities; a series of models of stellar populations I and II, for the mass range 0.25 to 64 Msun, followed from pre-main-sequence to a cool white dwarf or core collapse. An initial final-mass relationship is derived and compared with previous studies. Finally, we briefly address the evolution of non-canonical configurations, merger-products of low-mass main-sequence parents.

We present a calculation of the sedimentation of grains in a giant gaseous protoplanet such as that resulting from a disk instability of the type envisioned by Boss (1998). Boss (1998) has suggested that such protoplanets would form cores through the settling of small grains. We have tested this suggestion by following the sedimentation of small silicate grains as the protoplanet contracts and evolves. We find that during the course of the initial contraction of the protoplanet, which lasts some $4\times 10^5$ years, even very small (> 1 micron) silicate grains can sediment to create a core both for convective and non-convective envelopes, although the sedimentation time is substantially longer if the envelope is convective, and grains are allowed to be carried back up into the envelope by convection. Grains composed of organic material will mostly be evaporated before they get to the core region, while water ice grains will be completely evaporated. These results suggest that if giant planets are formed via the gravitational instability mechanism, a small heavy element core can be formed due to sedimentation of grains, but it will be composed almost entirely of refractory material. Including planetesimal capture, we find core masses between 1 and 10 M$_{\oplus}$, and a total high-Z enhancement of ~40 M$_{\oplus}$. The refractories in the envelope will be mostly water vapor and organic residuals.

Multicycle nova evolution models have been calculated over the past twenty years, the number being limited by numerical constraints. Here we present a long-term evolution code that enables a continuous calculation through an unlimited number of nova cycles for an unlimited evolution time, even up to (or exceeding) a Hubble time. Starting with two sets of the three independent nova parameters -- the white dwarf mass, the temperature of its isothermal core, and the rate of mass transfer on to it -- we have followed the evolution of two models, with initial masses of 1 and 0.65 solar masses, accretion rates (constant throughout each calculation) of 1e-11 and 1e-9 solar-masses/yr, and relatively high initial temperatures (as they are likely to be at the onset of the outburst phase), through over 1000 and over 3000 cycles, respectively. The results show that although on the short-term consecutive outbursts are almost identical, on the long-term scale the characteristics change. This is mainly due to the changing core temperature, which decreases very similarly to that of a cooling white dwarf for a time, but at a slower rate thereafter. As the white dwarf's mass continually decreases, since both models lose more mass than they accrete, the central pressure decreases accordingly. The outbursts on the massive white dwarf change gradually from fast to moderately fast, and the other characteristics (velocity, abundance ratios, isotopic ratios) change, too. Very slowly, a steady state is reached, where all characteristics, both in quiescence and in outburst, remain almost constant. For the less massive white dwarf accreting at a high rate, outbursts are similar throughout the evolution.

We follow the contraction and evolution of a typical Jupiter-mass clump created by the disk instability mechanism, and compute the rate of planetesimal capture during this evolution. We show that such a clump has a slow contraction phase lasting ~3x10^5 years. By following the trajectories of planetesimals as they pass through the envelope of the protoplanet, we compute the cross-section for planetesimal capture at all stages of the protoplanet's evolution. We show that the protoplanet can capture a large fraction of the solid material in its feeding zone, which will lead to an enrichment of the protoplanet in heavy elements. The exact amount of this enrichment depends upon, but is not very sensitive to the size and random speed of the planetesimals.

This paper is a sequel to an earlier paper devoted to multiple, multicycle nova evolution models (Prialnik & Kovetz 1995, first paper of the series), which showed that the different characteristics of nova outbursts can be reproduced by varying the values of three basic and independent parameters: the white dwarf mass-M_{WD}, the temperature of its isothermal core-T_{WD} and the mass transfer rate-Mdot. Apart from being the largest computational classical novae parameter-space survey, we show here that the parameter space is constrained by several analytical considerations and find its limiting surfaces. Consequently, we extend the grid of multicycle nova evolution models presented in Paper I almost to its limits, adding multicycle nova outburst calculations for a considerable number of new parameter combinations. In particular, the extended parameter pace that produces nova eruptions includes low mass transfer rates down to 5e-13 Msun/yr, and more models for low T_{WD}. Resulting characteristics of these runs are added to the former parameter combination results, to provide a full grid spanning the entire parameter space for Carbon-Oxygen white dwarfs. The full grid covers the entire range of observed nova characteristics, even those of peculiar objects, which have not been numerically reproduced until now. Most remarkably, runs for very low Mdot lead to very high values for some characteristics, such as outburst amplitude A>~20, high super-Eddington luminosities at maximum, heavy element abundance of the ejecta Z_{ej}~0.63 and high ejected masses m_{ej}~7e-4 Msun.