Search

The transition zone between the white dwarf (WD) envelope and a circumstellar accretion disk in classical novae, the boundary layer, is a region of strong dissipation and intense vorticity. In this strongly sheared layer, the hydrogen-rich accreted gas is expected to mix with the underlying WD outermost layers so the conditions for the onset of the thermonuclear runaway (TNR) in classical nova will be different from the the standard treatment of the onset and subsequent mixing. We applied the critical layer instability (CLI) to the boundary between a disk-accreted H/He zone and the C/O - or O/Ne - rich outer layers of a mass-accreting WD in a cataclysmic binary and then used the resulting structure as input to one-dimensional nuclear-hydrodynamic simulations of the nova outburst. We simulated the subsonic mixing process in two dimensions for conditions appropriate for the inner disk and a CO 0.8 solar mass and CO and ONe 1.25 solar mass WDs using the compressible hydrodynamics code PLUTO. The resulting compositional profile was then imported into the one-dimensional nuclear-hydrodynamics code SHIVA to simulate the triggering and growth rate for the TNR and subsequent envelope ejection. We find that the deep shear-driven mixing changes the triggering and development of the TNR. In particular, the time to reach peak temperature is significantly shorter, and the ejected mass and maximum velocity of the ejecta substantially greater, than the current treatment. The 7Li yield is reduced by about an order of magnitude relative to the current treatments., Comment: 12 pages, 6 figures

X-Ray bursts (XRB) are powerful thermonuclear events on the surface of accreting neutron stars (NS), where nucleosynthesis of intermediate-mass elements occurs. Their predicted and observed luminosities sometimes exceed Eddington's value, thus some of the material may escape by means of a stellar wind. This work seeks to determine the mass-loss and chemical composition of the material ejected through radiation-driven winds and its significance for Galactic abundances. It also reports on the evolution of pysical quantities during the wind phase that could help constrain the mass-radius relation in neutron stars. A non-relativistic radiative wind model was implemented and linked, through a new technique, to a series of XRB hydrodynamic simulations, that include over 300 isotopes. This allows us to construct a quasi-stationary time evolution of the wind during the XRB. The simulations resulted in the first realistic quantification of mass-loss for each isotope synthesized in the XRB. The total mass ejected by the wind was about $6\times10^{19}g$, the average ejected mass per unit time represents 2.6% of the accretion rate, with 0.1% of the envelope mass ejected per burst and ~90% of the ejecta composed by $^{60}$Ni, $^{64}$Zn, $^{68}$Ge and $^{58}$Ni. The ejected material also contained a small fraction ($10^{-4}-10^{-5}$) of some light p-nuclei, but not enough to account for their Galactic abundances. Additionally, the observable magnitudes during the wind phase showed remarkable correlations, some of which involve wind parameters like energy and mass outflows, that are determined by the conditions at the base of the wind envelope. These correlations could be used to link observable magnitudes to the physics of the innermost parts of the envelope, close to its interface with the NS crust. This is a promising result regarding the issue of NS radii determination.

Presolar stardust grains found in primitive meteorites are believed to retain the isotopic composition of stellar outflows at the time of grain condensation. Therefore, laboratory measurements of their isotopic ratios represent sensitive probes for investigating open questions related to stellar evolution, stellar explosions, nucleosynthesis, mixing mechanisms, dust formation, and galactic chemical evolution. For a few selected presolar grains, classical novae have been discussed as a potential source. For SiC, silicate, and graphite presolar grains, the association is based on the observation of small $N(^{12}$C)/$N(^{13}$C) and $N(^{14}$N)/$N(^{15}$N) number abundance ratios compared to solar values, and abundance excesses in $^{30}$Si relative to $^{29}$Si, as previously predicted by models of classical novae. We report on a direct measurement of the $^{29}$Si(p,$\gamma$)$^{30}$P reaction, which strongly impacts simulated $\delta ^{29}$Si values from classical novae. Our new experimental $^{29}$Si(p,$\gamma$)$^{30}$P thermonuclear reaction rate differs from previous results by up to 50\% in the classical nova temperature range ($T$ $=$ $100$ $-$ $400$~MK), while the rate uncertainty is reduced by up to a factor of $3$. Using our new reaction rate in Monte Carlo reaction network and hydrodynamic simulations of classical novae, we estimate $\delta ^{29}$Si values with much reduced uncertainties. Our results establish $\delta ^{29}$Si values measured in presolar grains as a sensitive probe for assessing their classical nova paternity. We also demonstrate that $\delta ^{30}$Si values from nova simulations are presently not a useful diagnostic tool unless the large uncertainty of the $^{30}$P(p,$\gamma$)$^{31}$S reaction rate can be significantly reduced., Comment: 5 figures

Photospheric radius expansion during X-ray bursts can be used to measure neutron star radii and help constrain the equation of state of neutron star matter. Understanding the stellar wind dynamics is important for interpreting observations. Stellar wind models, though studied in past decades, have thus regained interest and need to be revisited with updated data and methods. In this work we study the radiative wind model in the context of XRBs, with modern techniques and physics input. We focus on characterization of the solutions and study of observable magnitudes as a function of free model parameters. We implement a spherically-symmetric non-relativistic wind model in a stationary regime, with updated opacity tables and modern numerical techniques. Total mass and energy outflows $(\dot M,\dot E)$ are treated as free parameters. A high resolution parameter space exploration was performed to allow better characterization of observable magnitudes. High correlation was found between different photospheric magnitudes and free parameters. For instance, the photospheric ratio of gravitational energy outflow to radiative luminosity is in direct proportion to the photospheric wind velocity. The correlations found could help determine the physical conditions of the inner layers, where nuclear reactions take place, by means of observable photospheric values. Further studies are needed to determine the range of physical conditions in which the correlations are valid., Comment: Accepted in Astronomy and Astrophysics. New version with corrected typo in eq. (36) and some minor changes after language edition

High-resolution spectroscopy has revealed large concentrations of CNO and sometimes other intermediate-mass elements 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 the thermonuclear runaway. 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. This paper explores a new methodology that combines 1-D and 3-D simulations. The early stages of the explosion (i.e., mass-accretion and initiation of the runaway) have been computed with the 1-D hydrodynamic code SHIVA. When convection extends throughout the entire envelope, the structures for each model were mapped into 3-D Cartesian grids and were subsequently followed with the multidimensional code FLASH. Two key physical quantities were extracted from the 3-D simulations and 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. More massive envelopes than those reported from previous models with pre-enrichment have been found. This results in more violent outbursts, characterized by higher peak temperatures and greater ejected masses, with metallicity enhancements in agreement with observations., Comment: 8 pages, accepted for publication in Astronomy & Astrophysics

Gamma-ray astronomy explores the most energetic photons in nature to address some of the most pressing puzzles in contemporary astrophysics. It encompasses a wide range of objects and phenomena: stars, supernovae, novae, neutron stars, stellar-mass black holes, nucleosynthesis, the interstellar medium, cosmic rays and relativistic-particle acceleration, and the evolution of galaxies. MeV gamma-rays provide a unique probe of nuclear processes in astronomy, directly measuring radioactive decay, nuclear de-excitation, and positron annihilation. The substantial information carried by gamma-ray photons allows us to see deeper into these objects, the bulk of the power is often emitted at gamma-ray energies, and radioactivity provides a natural physical clock that adds unique information. New science will be driven by time-domain population studies at gamma-ray energies. This science is enabled by next-generation gamma-ray instruments with one to two orders of magnitude better sensitivity, larger sky coverage, and faster cadence than all previous gamma-ray instruments. This transformative capability permits: (a) the accurate identification of the gamma-ray emitting objects and correlations with observations taken at other wavelengths and with other messengers; (b) construction of new gamma-ray maps of the Milky Way and other nearby galaxies where extended regions are distinguished from point sources; and (c) considerable serendipitous science of scarce events -- nearby neutron star mergers, for example. Advances in technology push the performance of new gamma-ray instruments to address a wide set of astrophysical questions., Comment: 14 pages including 3 figures

In this White Paper we present the potential of the enhanced X-ray Timing and Polarimetry (eXTP) mission for studies related to Observatory Science targets. These include flaring stars, supernova remnants, accreting white dwarfs, low and high mass X-ray binaries, radio quiet and radio loud active galactic nuclei, tidal disruption events, and gamma-ray bursts. eXTP will be excellently suited to study one common aspect of these objects: their often transient nature. Developed by an international Consortium led by the Institute of High Energy Physics of the Chinese Academy of Science, the eXTP mission is expected to be launched in the mid 2020s., Comment: Accepted for publication on Sci. China Phys. Mech. Astron. (2019)

Multidimensional nucleosynthesis studies with hundreds of nuclei linked through thousands of nuclear processes are still computationally prohibitive. To date, most nucleosynthesis studies rely either on hydrostatic/hydrodynamic simulations in spherical symmetry, or on post-processing simulations using temperature and density versus time profiles directly linked to huge nuclear reaction networks. Parallel computing has been regarded as the main permitting factor of computationally intensive simulations. This paper explores the different pros and cons in the parallelization of stellar codes, providing recommendations on when and how parallelization may help in improving the performance of a code for astrophysical applications. We report on different parallelization strategies succesfully applied to the spherically symmetric, Lagrangian, implicit hydrodynamic code SHIVA, extensively used in the modeling of classical novae and type I X-ray bursts. Speed-up factors of 26 and 35 have been obtained with a parallelized version of SHIVA, in a 200-shell simulation of a type I X-ray burst carried out with two nuclear reaction networks: a reduced one, consisting of 324 isotopes and 1392 reactions, and a more extended network with 606 nuclides and 3551 nuclear interactions. Maximum speed-ups of 41 (324-isotope network) and 85 (606-isotope network), are also predicted for 200 cores, stressing that the number of shells of the computational domain constitutes an effective upper limit for the maximum number of cores that could be used in a parallel application.

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 a white dwarf (either CO or ONe). The material accumulates for 10+4 - 10+5 yr until ignition under degenerate conditions, resulting in a thermonuclear runaway. The nuclear energy released produces peak temperatures of about 0.1 - 0.4 GK. During these events, 10-7 - 10-3 Msun 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. Here we investigate the dredge-up, driven by Kelvin-Helmholtz instabilities, for white dwarf masses in the range 0.8-1.25 Msun and different core compositions (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 on the explosion characteristics. Results. We observe greater mixing (about 30 per cent higher when measured in the same conditions), at the time we stop the simulations, and more energetic outbursts for ONe-rich substrates than for CO-rich substrates and for more massive white dwarfs., Comment: Accepted for Publication in Astronomy and Astrophysics

INTEGRAL observed the nova V5668 Sgr around the time of its optical maximum on March 21, 2015. Studies at UV wavelengths showed spectral lines of freshly produced Be-7. This could be measurable also in gamma-rays at 478 keV from the decay to Li-7. Novae are also expected to synthesise Na-22 which decays to Ne-22, emitting a 1275 keV photon. About one week before the optical maximum, a strong gamma-ray flash on time-scales of hours is expected from short-lived radioactive nuclei, such as N-13 and F-18. These beta-plus-unstable nuclei should yield emission up to 511 keV, but which has never been observed. The spectrometer SPI aboard INTEGRAL pointed towards V5668 by chance. We use these observations to search for possible gamma-ray emission of decaying Be-7, and to directly measure the synthesised mass during explosive burning. We also aim to constrain possible burst-like emission days to weeks before the optical maximum using the SPI anticoincidence shield (ACS). We extract spectral and temporal information to determine the fluxes of gamma-ray lines at 478 keV, 511 keV, and 1275 keV. A measured flux value directly converts into abundances produced by the nova. The SPI-ACS rates are analysed for burst-like emission using a nova model light-curve. For the obtained nova flash candidate events, we discuss possible origins. No significant excess for the expected gamma-ray lines is found. Our upper limits on the synthesised Be-7 and Na-22 mass depend on the uncertainties of the distance to the nova: The Be-7 mass is constrained to less than $4.8\times10^{-9}\,(d/kpc)^2$, and Na-22 to less than $2.4\times10^{-8}\,(d/kpc)^2$ solar masses. For the Be-7 mass estimate from UV studies, the distance to V5668 Sgr must be larger than 1.2 kpc. During three weeks before the optical maximum, we find 23 burst-like events in the ACS rate, of which six could possibly be associated with V5668 Sgr., Comment: 11 pages, 7 figures, 1 table, accepted by A&A

About 30% to 40% of classical novae produce dust 20-100 days after the outburst, but no presolar stardust grains from classical novae have been unambiguously identified yet. Although several studies claimed a nova paternity for certain grains, the measured and simulated isotopic ratios could only be reconciled assuming that the grains condensed after the nova ejecta mixed with a much larger amount of close-to-solar matter. However, the source and mechanism of this potential post-explosion dilution of the ejecta remains a mystery. A major problem with previous studies is the small number of simulations performed and the implied poor exploration of the large nova parameter space. We report the results of a different strategy, based on a Monte Carlo technique, that involves the random sampling over the most important nova model parameters: the white dwarf composition; the mixing of the outer white dwarf layers with the accreted material before the explosion; the peak temperature and density; the explosion time scales; and the possible dilution of the ejecta after the outburst. We discuss and take into account the systematic uncertainties for both the presolar grain measurements and the simulation results. Only those simulations that are consistent with all measured isotopic ratios of a given grain are accepted for further analysis. We also present the numerical results of the model parameters. We identify 18 presolar grains with measured isotopic signatures consistent with a CO nova origin, without assuming any dilution of the ejecta. Among these, the grains G270 2, M11-334-2, G278, M11-347-4, M11-151-4, and Ag2 6 have the highest probability of a CO nova paternity., Comment: 8 figures

Context. Classical novae are thermonuclear explosions hosted by accreting white dwarfs in stellar binary systems. Material piles up on top of the white dwarf star under mildly degenerate conditions, driving a thermonuclear runaway. The energy released by the suite of nuclear processes operating at the envelope (mostly proton-capture reactions and beta-decays) heats the material up to peak temperatures ranging from 100 to 400 MK. In these events, about 10-3 - 10-7 Msun, enriched in CNO and, sometimes, other intermediate-mass elements (e.g., Ne, Na, Mg, Al), are ejected into the interstellar medium. Aims. To date, most of the efforts undertaken in the modeling of classical nova outbursts have focused on the early stages of the explosion and ejection, ignoring the interaction of the ejecta, first with the accretion disk orbiting the white dwarf, and ultimately with the secondary star. Methods. A suite of three-dimensional, SPH simulations of the interaction between the nova ejecta, the accretion disk, and the stellar companion have been performed to fill this gap, aimed at testing the influence of the different parameters (i.e., mass and velocity of the ejecta, mass and geometry of the accretion disk) on the dynamical and chemical properties of the system. Results. We discuss the conditions that lead to the disruption of the accretion disk and to mass loss from the binary system. In addition, we discuss the likelihood of chemical contamination of the stellar secondary induced by the impact with the nova ejecta and its potential effect on the next nova cycle., Comment: Accepted for Publication in Astronomy and Astrophysics

Classical novae are thermonuclear explosions that take place in the envelopes of accreting white dwarfs in binary systems. The material piles up under degenerate conditions, driving a thermonuclear runaway. The energy released by the suite of nuclear processes operating at the envelope heats the material up to peak temperatures about 100 - 400 MK. During these events, about 10-3 - 10-7 Msun, enriched in CNO and, sometimes, other intermediate-mass elements (e.g., Ne, Na, Mg, Al) are ejected into the interstellar medium. To account for the gross observational properties of classical novae (in particular, the large concentrations of metals spectroscopically inferred in the ejecta), models require mixing between the (solar-like) material transferred from the secondary and the outermost layers (CO- or ONe-rich) of the underlying white dwarf. Recent multidimensional simulations have demonstrated that Kelvin-Helmholtz instabilities can naturally produce self-enrichment of the accreted envelope with material from the underlying white dwarf at levels that agree with observations. However, the feasibility of this mechanism has been explored in the framework of CO white dwarfs, while mixing with different substrates still needs to be properly addressed. Three-dimensional simulations of mixing at the core-envelope interface during nova outbursts have been performed with the multidimensional code FLASH, for two types of substrates: CO- and ONe-rich. We show that the presence of an ONe-rich substrate, as in "neon novae", yields larger metallicity enhancements in the ejecta, compared to CO,rich substrates (i.e., non-neon novae). A number of requirements and constraints for such 3-D simulations (e.g., minimum resolution, size of the computational domain) are also outlined., Comment: Accepted for publication in Astronomy & Astrophysics

Presolar grains constitute remnants of stars that existed before the formation of the solar system. In addition to providing direct information on the materials from which the solar system formed, these grains provide ground-truth information for models of stellar evolution and nucleosynthesis. Here we report the in-situ identification of two unique presolar graphite grains from the primitive meteorite LaPaz Icefield 031117. Based on these two graphite grains, we estimate a bulk presolar graphite abundance of 5(-3)(+7) ppm in this meteorite. One of the grains (LAP-141) is characterized by an enrichment in 12C and depletions in 33,34S, and contains a small iron sulfide subgrain, representing the first unambiguous identification of presolar iron sulfide. The other grain (LAP-149) is extremely 13C-rich and 15N-poor, with one of the lowest 12C/13C ratios observed among presolar grains. Comparison of its isotopic compositions with new stellar nucleosynthesis and dust condensation models indicates an origin in the ejecta of a low-mass CO nova. Grain LAP-149 is the first putative nova grain that quantitatively best matches nova model predictions, providing the first strong evidence for graphite condensation in nova ejecta. Our discovery confirms that CO nova graphite and presolar iron sulfide contributed to the original building blocks of the solar system., Comment: Accepted for publication in The Astrophysical Journal

Context. Classical novae are thermonuclear explosions that take place in the envelopes of accreting white dwarfs in stellar binary systems. The material transferred onto the white dwarf piles up under degenerate conditions, driving a thermonuclear runaway. In those outbursts, about 10-7 - 10-3 Msun, enriched in CNO and, sometimes, other intermediate-mass elements (e.g., Ne, Na, Mg, or Al, for ONe novae) are ejected into the interstellar medium. The large concentrations of metals spectroscopically inferred in the nova ejecta reveal that the (solar-like) material transferred from the secondary mixes with the outermost layers of the underlying white dwarf. Aims. Most theoretical models of nova outbursts reported to date yield, on average, outflows characterized by O > C, from which only oxidized condensates (e.g, O-rich grains) would be expected, in principle. Methods. To specifically address whether CO novae can actually produce C-rich dust, six different hydrodynamic nova models have been evolved, from accretion to the expansion and ejection stages, with different choices for the composition of the substrate with which the solar-like accreted material mixes. Updated chemical profiles inside the H-exhausted core have been used, based on stellar evolution calculations for a progenitor of 8 Msun through H and He-burning phases. Results. We show that these profiles lead to C-rich ejecta after the nova outburst. This extends the possible contribution of novae to the inventory of presolar grains identified in meteorites, particularly in a number of carbonaceous phases (i.e., nanodiamonds, silicon carbides and graphites)., Comment: 5 pages, accepted for publication in Astronomy & Astrophysics

Extreme excesses of $^{13}C$ ($^{12}C$/$^{13}C$<10) and $^{15}N$ ($^{14}N$/$^{15}N$<20) in rare presolar SiC grains have been considered diagnostic of an origin in classical novae, though an origin in core collapse supernovae (CCSNe) has also been proposed. We report C, N, and Si isotope data for 14 submicron- to micron-sized $^{13}C$- and $^{15}N$-enriched presolar SiC grains ($^{12}C$/$^{13}C$<16 and $^{14}N$/$^{15}N$<~100) from Murchison, and their correlated Mg-Al, S, and Ca-Ti isotope data when available. These grains are enriched in $^{13}C$ and $^{15}N$, but with quite diverse Si isotopic signatures. Four grains with $^{29,30}Si$ excesses similar to those of type C SiC grains likely came from CCSNe, which experienced explosive H burning occurred during explosions. The independent coexistence of proton- and neutron-capture isotopic signatures in these grains strongly supports heterogeneous H ingestion into the He shell in pre-supernovae. Two of the seven putative nova grains with $^{30}Si$ excesses and $^{29}Si$ depletions show lower-than-solar $^{34}S$/$^{32}S$ ratios that cannot be explained by classical nova nucleosynthetic models. We discuss these signatures within the CCSN scenario. For the remaining five putative nova grains, both nova and supernova origins are viable because explosive H burning in the two stellar sites could result in quite similar proton-capture isotopic signatures. Three of the grains are sub-type AB grains that are also $^{13}C$ enriched, but have a range of higher $^{14}N$/$^{15}N$. We found that $^{15}N$-enriched AB grains (~50<$^{14}N$/$^{15}N$<~100) have distinctive isotopic signatures compared to putative nova grains, such as higher $^{14}N$/$^{15}N$, lower $^{26}Al$/$^{27}Al$, and lack of $^{30}Si$ excess, indicating weaker proton-capture nucleosynthetic environments., Comment: fix typo in one of the authors' names

Aims: The aim of this work is to compare the performance of three reaction network integration methods used in stellar nucleosynthesis calculations. These are the Gear's backward differentiation method, Wagoner's method (a 2nd-order Runge-Kutta method), and the Bader-Deuflehard semi-implicit multi-step method. Methods: To investigate the efficiency of each of the integration methods considered here, a test suite of temperature and density versus time profiles is used. This suite provides a range of situations ranging from constant temperature and density to the dramatically varying conditions present in white dwarf mergers, novae, and x-ray bursts. Some of these profiles are obtained separately from full hydrodynamic calculations. The integration efficiencies are investigated with respect to input parameters that constrain the desired accuracy and precision. Results: Gear's backward differentiation method is found to improve accuracy, performance, and stability in integrating nuclear reaction networks. For temperature-density profiles that vary strongly with time, it is found to outperform the Bader-Deuflehard method (although that method is very powerful for more smoothly varying profiles). Wagoner's method, while relatively fast for many scenarios, exhibits hard-to-predict inaccuracies for some choices of integration parameters owing to its lack of error estimations., Comment: 13 pages, 12 figures, accepted to Astronomy and Astrophysics (section 15) - corrected units in Figs. 6-11

Classical novae are caused by mass transfer episodes from a main sequence star onto a white dwarf via Roche lobe overflow. This material forms an accretion disk around the white dwarf. Ultimately, a fraction of this material spirals in and piles up on the white dwarf surface under electron-degenerate conditions. The subsequently occurring thermonuclear runaway reaches hundreds of megakelvin and explosively ejects matter into the interstellar medium. The exact peak temperature strongly depends on the underlying white dwarf mass, the accreted mass and metallicity, and the initial white dwarf luminosity. Observations of elemental abundance enrichments in these classical nova events imply that the ejected matter consists not only of processed solar material from the main sequence partner but also of material from the outer layers of the underlying white dwarf. This indicates that white dwarf and accreted matter mix prior to the thermonuclear runaway. The processes by which this mixing occurs require further investigation to be understood. In this work, we analyze elemental abundances ejected from hydrodynamic nova models in search of elemental abundance ratios that are useful indicators of the total amount of mixing. We identify the abundance ratios $\Sigma$CNO/H, Ne/H, Mg/H, Al/H, and Si/H as useful mixing meters in ONe novae. The impact of thermonuclear reaction rate uncertainties on the mixing meters is investigated using Monte Carlo post-processing network calculations with temperature-density evolutions of all mass zones computed by the hydrodynamic models. We find that the current uncertainties in the $^{30}$P($p$,$\gamma$)$^{31}$S rate influence the Si/H abundance ratio, but overall the mixing meters found here are robust against nuclear physics uncertainties. A comparison of our results with observations of ONe novae provides strong constraints for classical nova models.

In old, dense stellar systems collisions of white dwarfs are a rather frequent phenomenon. Here we present the results of a comprehensive set of Smoothed Particle Hydrodynamics simulations of close encounters of white dwarfs aimed to explore the outcome of the interaction and the nature of the final remnants for different initial conditions. Depending on the initial conditions and the white dwarf masses, three different outcomes are possible. Specifically, the outcome of the interaction can be either a direct or a lateral collision or the interaction can result in the formation of an eccentric binary system. In those cases in which a collision occurs, the infalling material is compressed and heated such that the physical conditions for a detonation may be reached during the most violent phases of the merger. While we find that detonations occur in a significant number of our simulations, in some of them the temperature increase in the shocked region rapidly lifts degeneracy, leading to the quenching of the burning. We thus characterize under which circumstances a detonation is likely to occur as a result of the impact of the disrupted star on the surface of the more massive white dwarf. Finally, we also study which interactions result in bound systems, and in which ones the more massive white dwarf is also disrupted as a consequence of the dynamical interaction. The sizable number of simulations performed in this work allows to find how the outcome of the interaction depends on the distance at closest approach, and on the masses of the colliding white dwarfs, and which is the chemical pattern of the nuclearly processed material. Finally, we also discuss the influence of the masses and core chemical compositions of the interacting white dwarfs and the different kinds of impact in the properties of the remnants., Comment: 18 pages, 6 figures. Accepted for publication in MNRAS

The impact of nuclear physics uncertainties on nucleosynthesis in thermonuclear supernovae has not been fully explored using comprehensive and systematic studies with multiple models. To better constrain predictions of yields from these phenomena, we have performed a sensitivity study by post-processing thermodynamic histories from two different hydrodynamic, Chandrasekhar-mass explosion models. We have individually varied all input reaction and, for the first time, weak interaction rates by a factor of ten and compared the yields in each case to yields using standard rates. Of the 2305 nuclear reactions in our network, we find that the rates of only 53 reactions affect the yield of any species with an abundance of at least 10^-8 M_sun by at least a factor of two, in either model. The rates of the 12C(a,g), 12C+12C, 20Ne(a,p), 20Ne(a,g) and 30Si(p,g) reactions are among those that modify the most yields when varied by a factor of ten. From the individual variation of 658 weak interaction rates in our network by a factor of ten, only the stellar 28Si(b+)28Al, 32S(b+)32P and 36Ar(b+)36Cl rates significantly affect the yields of species in a model. Additional tests reveal that reaction rate changes over temperatures T > 1.5 GK have the greatest impact, and that ratios of radionuclides that may be used as explosion diagnostics change by a factor of less than two from the variation of individual rates by a factor of 10. Nucleosynthesis in the two adopted models is relatively robust to variations in individual nuclear reaction and weak interaction rates. Laboratory measurements of a limited number of reactions would help to further constrain predictions. As well, we confirm the need for a consistent treatment for relevant stellar weak interaction rates since simultaneous variation of these rates (as opposed to individual variation) has a significant effect on yields in our models., Comment: accepted by A&A, 14 pages, 5 figures, 2 tables

Classical novae are stellar explosions occurring in binary systems, consisting of a white dwarf and a main sequence companion. Thermonuclear runaways on the surface of massive white dwarfs, consisting of oxygen and neon, are believed to reach peak temperatures of several hundred million kelvin. These temperatures are strongly correlated with the underlying white dwarf mass. The observational counterparts of such models are likely associated with outbursts that show strong spectral lines of neon in their shells (neon novae). The goals of this work are to investigate how useful elemental abundances are for constraining the peak temperatures achieved during these outbursts and determine how robust "nova thermometers" are with respect to uncertain nuclear physics input. We present updated observed abundances in neon novae and perform a series of hydrodynamic simulations for several white dwarf masses. We find that the most useful thermometers, N/O, N/Al, O/S, S/Al, O/Na, Na/Al, O/P, and P/Al, are those with the steepest monotonic dependence on peak temperature. The sensitivity of these thermometers to thermonuclear reaction rate variations is explored using post-processing nucleosynthesis simulations. The ratios N/O, N/Al, O/Na, and Na/Al are robust, meaning they are minimally affected by uncertain rates. However, their dependence on peak temperature is relatively weak. The ratios O/S, S/Al, O/P, and P/Al reveal strong dependences on temperature and the poorly known 30P(p,g)31S rate. We compare our model predictions to neon nova observations and obtain the following estimates for the underlying white dwarf masses: 1.34-1.35 solar masses (V838 Her), 1.18-1.21 solar masses (V382 Vel), <1.3 solar masses (V693 CrA), <1.2 solar masses (LMC 1990#1), and <1.2 solar masses (QU Vul)., Comment: 12 pages, 7 figures, accepted to ApJ

The origin of R Coronae Borealis stars has been elusive for over 200 years. Currently, two theories for their formation have been presented. These are the Final Flash scenario, in which a dying asymptotic giant branch (AGB) star throws off its atmosphere to reveal the hydrogen poor, heavily processed material underneath, and the double degenerate scenario, in which two white dwarfs merge to produce a new star with renewed vigour. Some theories predict that the temperatures reached during the latter scenario would destroy any lithium originally present in the white dwarfs. The observed lithium content of some R Coronae Borealis stars, therefore, is often interpreted as an indication that the Final Flash scenario best describes their formation. In this paper, it is shown that lithium production can, indeed, occur in the merging of a helium white dwarf with a carbon-oxygen white dwarf if their chemical composition, particularly that of 3He, is fully considered. The production mechanism is described in detail, and the sensitivity of lithium production to the merger environment is investigated. Nucleosynthesis post-processing calculations of smoothed-particle hydrodynamics (SPH) tracer particles are performed to show that any lithium produced in these environments will be concentrated towards the cloud of material surrounding the R CrB star. Measurements of the lithium content of these clouds would, therefore, provide a valuable insight into the formation mechanism of R CrB stars., Comment: 6 Pages, 2 Figures. To be published in Astronomy and Astrophysics

The Rapid Burster (MXB 1730-335) is a unique object, showing both type I and type II X-ray bursts. A type I burst of the Rapid Burster was observed with Swift/XRT on 2009 March 5, showing photospheric radius expansion for the first time in this source. We report here on the mass and radius determination from this photospheric radius expansion burst using a Bayesian approach. After marginalization over the likely distance of the system (5.8-10 kpc) we obtain M=1.1+/-0.3 M_sun and R=9.6+/-1.5 km (1-sigma uncertainties) for the compact object, ruling out the stiffest equations of state for the neutron star. We study the sensitivity of the results to the distance, the color correction factor, and the hydrogen mass fraction in the envelope. We find that only the distance plays a crucial role., Comment: 7 pages, 8 figures, Accepted for publication in The Astrophysical Journal

Half a century has passed since the foundation of nuclear astrophysics. Since then, this discipline has reached its maturity. Today, nuclear astrophysics constitutes a multidisciplinary crucible of knowledge that combines the achievements in theoretical astrophysics, observational astronomy, cosmochemistry and nuclear physics. New tools and developments have revolutionized our understanding of the origin of the elements: supercomputers have provided astrophysicists with the required computational capabilities to study the evolution of stars in a multidimensional framework; the emergence of high-energy astrophysics with space-borne observatories has opened new windows to observe the Universe, from a novel panchromatic perspective; cosmochemists have isolated tiny pieces of stardust embedded in primitive meteorites, giving clues on the processes operating in stars as well as on the way matter condenses to form solids; and nuclear physicists have measured reactions near stellar energies, through the combined efforts using stable and radioactive ion beam facilities. This review provides comprehensive insight into the nuclear history of the Universe and related topics: starting from the Big Bang, when the ashes from the primordial explosion were transformed to hydrogen, helium, and few trace elements, to the rich variety of nucleosynthesis mechanisms and sites in the Universe. Particular attention is paid to the hydrostatic processes governing the evolution of low-mass stars, red giants and asymptotic giant-branch stars, as well as to the explosive nucleosynthesis occurring in core-collapse and thermonuclear supernovae, gamma-ray bursts, classical novae, X-ray bursts, superbursts, and stellar mergers., Comment: Invited Review. Accepted for publication in "Reports on Progress in Physics" (version with low-resolution figures)

We present a Monte Carlo calculation of the astrophysical rate of the 15O(alpha,gamma)19Ne reaction based on an evaluation of published experimental data. By considering the likelihood distributions of individual resonance parameters derived from measurements, estimates of upper and lower limits on the reaction rate at the 99.73% confidence level are derived in addition to the recommended, median value. These three reaction rates are used as input for three separate calculations of Type I x-ray bursts using spherically symmetric, hydrodynamic simulations of an accreting neutron star. In this way the influence of the 15O(alpha,gamma)19Ne reaction rate on the peak luminosity, recurrence time, and associated nucleosynthesis in models of Type I x-ray bursts is studied. Contrary to previous findings, no substantial effect on any of these quantities is observed in a sequence of four bursts when varying the reaction rate between its lower and upper limits. Rather, the differences in these quantities are comparable to the burst-to-burst variations with a fixed reaction rate, indicating that uncertainties in the 15O(alpha,gamma)19Ne reaction rate do not strongly affect the predictions of this Type I x-ray burst model., Comment: Accepted for publication by the Astrophysical Journal

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 a white dwarf. The material piles up under degenerate conditions and a thermonuclear runaway ensues. The energy released by the suite of nuclear processes operating at the envelope heats the material up to peak temperatures of ~ (1 - 4) \times 108 K. During these events, about 10-4 - 10-5 M\odot, enriched in CNO and other intermediate-mass elements, are ejected into the interstellar medium. To account for the gross observational properties of classical novae (in particular, a metallicity enhancement in the ejecta above solar values), numerical models assume mixing between the (solar-like) material transferred from the companion and the outermost layers (CO- or ONe-rich) of the underlying white dwarf. The nature of the mixing mechanism that operates at the core-envelope interface has puzzled stellar modelers for about 40 years. Here we investigate the role of Kelvin-Helmholtz instabilities as a natural mechanism for self-enrichment of the accreted envelope with core material. The feasibility of this mechanism is studied by means of the multidimensional code FLASH. Here, we present a series of 9 numerical simulations perfomed in two dimensions aimed at testing the possible influence of the initial perturbation (duration, strength, location, and size), the resolution adopted, or the size of the computational domain on the results. We show that results do not depend substantially on the specific choice of these parameters, demonstrating that Kelvin- Helmholtz instabilities can naturally lead to self-enrichment of the accreted envelope with core material, at levels that agree with observations., Comment: Accepted for publication in A&A

Type I X-ray bursts are thermonuclear stellar explosions driven by charged-particle reactions. In the regime for combined H/He-ignition, the main nuclear flow is dominated by the rp-process (rapid proton-captures and beta+ decays), the 3 alpha-reaction, and the alpha-p-process (a suite of (alpha,p) and (p,gamma) reactions). The main flow is expected to proceed away from the valley of stability, eventually reaching the proton drip-line beyond A = 38. Detailed analysis of the relevant reactions along the main path has only been scarcely addressed, mainly in the context of parameterized one-zone models. In this paper, we present a detailed study of the nucleosynthesis and nuclear processes powering type I X-ray bursts. The reported 11 bursts have been computed by means of a spherically symmetric (1D), Lagrangian, hydrodynamic code, linked to a nuclear reaction network that contains 325 isotopes (from 1H to 107Te), and 1392 nuclear processes. These evolutionary sequences, followed from the onset of accretion up to the explosion and expansion stages, have been performed for 2 different metallicities to explore the dependence between the extension of the main nuclear flow and the initial metal content. We carefully analyze the dominant reactions and the products of nucleosynthesis, together with the the physical parameters that determine the light curve (including recurrence times, ratios between persistent and burst luminosities, or the extent of the envelope expansion). Results are in qualitative agreement with the observed properties of some well-studied bursting sources. Leakage from the predicted SbSnTe-cycle cannot be discarded in some of our models. Production of 12C (and implications for the mechanism that powers superbursts), light p-nuclei, and the amount of H left over after the bursting episodes will also be discussed., Comment: 78 pages (pdf), including 34 figures. Accepted for publication in The Astrophysical Journal Suppl. Series

Classical novae are powered by thermonuclear runaways that occur on the white dwarf component of close binary systems. During these violent stellar events, whose energy release is only exceeded by gamma-ray bursts and supernova explosions, about 10-4 10-5 Msun of material is ejected into the interstellar medium. Because of the high peak temperatures attained during the explosion, Tpeak ~ (1-4)x10+8 K, the ejecta are enriched in nuclear-processed material relative to solar abundances, containing significant amounts of 13C, 15N, and 17O and traces of other isotopes. The origin of these metal enhancements observed in the ejecta is not wellknown and has puzzled theoreticians for about 40 years. In this paper, we present new 2-D simulations of mixing at the core-envelope interface. We show that Kelvin-Helmholtz instabilities can naturally lead to self-enrichment of the solar-like accreted envelopes with material from the outermost layers of the underlying white dwarf core, at levels that agree with observations.

We examine the sensitivity of nucleosynthesis in Type I X-ray bursts to variations in nuclear rates. As a large number of nuclear processes are involved in these phenomena -with the vast majority of reaction rates only determined theoretically due to the lack of any experimental information- our results can provide a means for determining which rates play significant roles in the thermonuclear runaway. These results may then motivate new experiments. For our studies, we have performed a comprehensive series of one-zone post-processing calculations in conjunction with various representative X-ray burst thermodynamic histories. We present those reactions whose rate variations have the largest effects on yields in our studies., Comment: 8 pages, accepted for publication in New Astronomy Reviews, Special Issue on "Astronomy with Radioactivities VI" workshop, Ringberg Castle, Germany, Jan. 2008

Type I X-ray bursts are violent stellar events that take place on the H/He-rich envelopes of accreting neutron stars. We have investigated the role played by uncertainties in nuclear processes on the nucleosynthesis accompanying these explosive phenomena. Two different approaches have been adopted, in the framework of post-processing calculations. In the first one, nuclear rates are varied individually within uncertainties. Ten different models, covering the characteristic parameter space for these stellar events, have been considered. The second, somewhat complementary approach involves a Monte Carlo code in which all nuclear rates are randomly varied within uncertainty limits simultaneously. All in all, about 50,000 post-processing calculations, with a network containing 606 nuclides (H to 113Xe) and more than 3500 nuclear processes, have been performed in this work. A brief comparison between both procedures is outlined together with an overall account of the key nuclear reactions whose uncertainties have the largest impact in our X-ray burst nucleosynthesis studies., Comment: 91 pages, submitted to "The Astrophysical Journal Suppl. Series"

Classical novae are fascinating stellar events, at the crossroads of astrophysics, nuclear physics and cosmochemistry. In this review, we outline the history of nova modeling with special emphasis on recent advances and perspectives in multidimensional simulations. Among the topics that are covered, we analyze the interplay between nova outbursts and the Galactic chemical abundances, the synthesis of radioactive nuclei of interest for gamma-ray astronomy, such as 7Li, 22Na or 26Al, and the recent discovery of presolar meteoritic grains, likely condensed in nova shells., Comment: Invited review presented at the "Nuclei in the Cosmos VIII", Vancouver, july 19-23, 2004

Infrared and ultraviolet observations of nova light curves have confirmed grain formation in their expanding shells that are ejected into the interstellar medium by a thermonuclear runaway. In this paper, we present isotopic ratios of intermediate-mass elements up to silicon for the ejecta of CO and ONe novae, based on 20 hydrodynamic models of nova explosions. These theoretical estimates will help to properly identify nova grains in primitive meteorites. In addition, equilibrium condensation calculations are used to predict the types of grains that can be expected in the nova ejecta, providing some hints on the puzzling formation of C-rich dust in O>C environments. These results show that SiC grains can condense in ONe novae, in concert with an inferred (ONe) nova origin for several presolar SiC grains., Comment: 42 pages. Accepted for publication in The Astrophysical Journal

We analyse the impact of the initial abundances of the underlying white dwarf core on the nucleosynthesis accompanying classical nova outbursts, in the framework of hydrodynamic models of the explosion. Specifically, we take into account the chemical stratification of the white dwarf. It turns out that the presence of a thick CO-buffer on top of the ONe-rich core, as results from detailed calculations of previous evolution of the progenitor star, may lead to significant amounts of both 7Li and 26Al, after an outburst that, due to the lack of neon isotopes in the ejecta, would be misclassified as a non-neon nova (i.e., CO nova)., Comment: Accepted for publication in The Astrophysical Journal Letters

In the present paper, we study the Poincare map associated to a periodic perturbation, both in space and time, of a linear Hamiltonian system. The dynamical system embodies the essential physics of stellar pulsations and provides a global and qualitative explanation of the chaotic oscillations observed in some stars. We show that this map is an area preserving one with an oscillating rotation number function. The nonmonotonic property of the rotation number function induced by the triplication of the elliptic fixed point is superposed on the nonmonotonic character due to the oscillating perturbation. This superposition leads to the co-manifestation of generic phenomena such as reconnection and meandering, with the nongeneric scenario of creation of vortices. The nonmonotonic property due to the triplication bifurcation is shown to be different from that exhibited by the cubic Henon map, which can be considered as the prototype of area preserving maps which undergo a triplication followed by the twistless bifurcation. Our study exploits the reversibility property of the initial system, which induces the time-reversal symmetry of the Poincare map., Comment: To appear in Communications in Nonlinear Science and Numerical Simulation Vol.8 (2003); Figures 11,14,15,16 of better resolution are available from the authors

There is growing observational evidence that the irregular changes in the light curves of certain variable stars might be due to deterministic chaos. Supporting these conclusions, several simple models of non-linear oscillators have been shown to be capable of reproducing the observed complex behaviour. In this work, we introduce a non-linear, non-adiabatic one-zone model intended to reveal the factors leading to irregular luminosity variations in some pulsating stars. We have studied and characterized the dynamical behaviour of the oscillator as the input parameters are varied. The parametric study implied values corresponding to stellar models in the family of Long Period Variables and in particular of Mira-type stars. We draw the attention on certain solutions that reproduce with reasonable accuracy the observed behaviour of some peculiar Mira variables., Comment: 8 pages, 6 figures; accepted for publication in MNRAS

One of the observational evidences in support of the "thermonuclear runaway model" for the classical nova outburst relies on the accompanying nucleosynthesis. In this paper, we stress the relevant role played by nucleosynthesis in our understanding of the nova phenomenon by constraining models through a comparison with both the atomic abundance determinations from the ejecta and the isotopic ratios measured in presolar grains of a likely nova origin. Furthermore, the endpoint of nova nucleosynthesis provides hints for the understanding of the mixing process responsible for the enhanced metallicities found in the ejecta, and reveals also information on the properties of the underlying white dwarf (mass, luminosity...). We discuss first the interplay between nova outbursts and the Galactic chemical abundances: Classical nova outbursts are expected to be the major source of 13C, 15N and 17O in the Galaxy, and to contribute to the abundances of other species with A < 40, such as 7Li or 26Al. We describe the main nuclear path during the course of the explosion, with special emphasis on the synthesis of radioactive species, of particular interest for the gamma-ray output predicted from novae (7Li, 18F, 22Na, 26Al). An overview of the recent discovery of presolar nova candidate grains, as well as a discussion of the role played by nuclear uncertainties associated with key reactions of the NeNa-MgAl and Si-Ca regions, are also given., Comment: 10 pages. To appear in "Classical Nova Explosions", M. Hernanz and J. Jose (eds.), AIP

When intermediate mass stars reach their last stages of evolution they show pronounced oscillations. This phenomenon happens when these stars reach the so-called Asymptotic Giant Branch (AGB), which is a region of the Hertzsprung-Russell diagram located at about the same region of effective temperatures but at larger luminosities than those of regular giant stars. The period of these oscillations depends on the mass of the star. There is growing evidence that these oscillations are highly correlated with mass loss and that, as the mass loss increases, the pulsations become more chaotic. In this paper we study a simple oscillator which accounts for the observed properties of this kind of stars. This oscillator was first proposed and studied by Icke et al. [Astron.Astrophys. 258, 341 (1992)] and we extend their study to the region of more massive and luminous stars - the region of Super-AGB stars. The oscillator consists of a periodic nonlinear perturbation of a linear Hamiltonian system. The formalism of dynamical systems theory has been used to explore the associated Poincare map for the range of parameters typical of those stars. We have studied and characterized the dynamical behaviour of the oscillator as the parameters of the model are varied, leading us to explore a sequence of local and global bifurcations. Among these, a tripling bifurcation is remarkable, which allows us to show that the Poincare map is a nontwist area preserving map. Meandering curves, hierarchical-islands traps and sticky orbits also show up. We discuss the implications of the stickiness phenomenon in the evolution and stability of the Super-AGB stars., Comment: 13 pages, 9 figures

We investigate the effects of thermonuclear reaction rate uncertainties on nova nucleosynthesis. One-zone nucleosynthesis calculations have been performed by adopting temperature-density-time profiles of the hottest hydrogen-burning zone (i.e., the region in which most of the nucleosynthesis takes place). We obtain our profiles from 7 different, recently published, hydrodynamic nova simulations covering peak temperatures in the range from Tpeak=0.145-0.418 GK. For each of these profiles, we individually varied the rates of 175 reactions within their associated errors and analyzed the resulting abundance changes of 142 isotopes in the mass range below A=40. In total, we performed 7350 nuclear reaction network calculations. We use the most recent thermonuclear reaction rate evaluations for the mass ranges A=1-20 and A=20-40. For the theoretical astrophysicist, our results indicate the extent to which nova nucleosynthesis calculations depend on presently uncertain nuclear physics input, while for the experimental nuclear physicist our results represent at least a qualitative guide for future measurements at stable and radioactive ion beam facilities. We find that present reaction rate estimates are reliable for predictions of Li, Be, C and N abundances in nova nucleosynthesis. However, rate uncertainties of several reactions have to be reduced significantly in order to predict more reliable O, F, Ne, Na, Mg, Al, Si, S, Cl and Ar abundances. Results are presented in tabular form for each adopted nova simulation., Comment: Accepted for publication in The Astrophysical Journal Suppl. Series

He-accreting white dwarfs with sub-Chandrasekhar mass are revisited. The impact of the use of an extended reaction network on the predicted energy production and characteristics of the detonating layers is studied. It is shown that the considered scenario can be the site of an $\alpha$p-process combined with a p-process and with a variant of the rp-process we refer to as the pn-process. We define the conditions under which the derived distribution of the abundances of the p-nuclides in the ejecta, including the puzzling light Mo and Ru isotopes, mimics the solar-system one., Comment: 4 pages; A&A Letters in press

Thermonuclear runaways driven by accretion into degenerate white dwarf cores are the source that power classical nova outbursts. In this paper, we identify the dominant nuclear paths involved in the synthesis of intermediate-mass elements, from Si to Ca, during such violent events. New evolutionary sequences of 1.35 Mo ONe novae have been computed, using updated nuclear reaction rates. The main nuclear activity in this region is powered by the leakage from the NeNa-MgAl region, where the activity is confined during the early stages of the explosion. We discuss the critical role played by 30P(p,gamma) in the synthesis of nuclear species beyond sulfur and point out the large uncertainty that affects its rate, which has dramatic consequences for studies of nova nucleosynthesis in the Si-Ca mass region., Comment: Accepted for publication in The Astrophysical Journal. 23 pages (aastex), 6 Postscript figures

Gamma-ray emission from classical novae is dominated, during the first hours, by positron annihilation resulting from the beta decay of radioactive nuclei. The main contribution comes from the decay of 18F and hence is directly related to 18F formation during the outburst. A good knowledge of the nuclear reaction rates of production and destruction of 18F is required to study 18F synthesis in novae and the resulting gamma-ray emission. The rates relevant for the main mode of 18F destruction (i.e, through proton captures) have been the object of many recent experiments. However, subsequent analyses were focused on providing rates for X-ray burst nucleosynthesis not valid at nova temperatures (lower than 3.5 10^8 K). Accordingly, it is crucial to propose and discuss new reaction rates, incorporating all new experimental results, down to the domain of nova nucleosynthesis. We show that in this temperature regime, the 18F(p,gamma) and (p,alpha) reaction rates remain uncertain and deserve further experimental and theoretical efforts. Our hydrodynamic calculations including the new nuclear rates demonstrate that their impact on 18F synthesis in nova explosions is quite large and, consequently, the early gamma-ray emission from classical novae is also affected., Comment: To appear in Astron. Astrophys., 14 pages with 9 figures

The mechanism of classical novae explosions is explained, together with some of their observational properties. The scarce but not null impact of novae in the chemical evolution of the Milky Way is analyzed, as well as their relevance for the radioactivity in the Galaxy. A special emphasis is given to the predicted gamma-ray emission from novae and its relationship with the thermonuclear model itself and its related nucleosynthesis., Comment: 10 pages, 2 figures, LaTex (aipproc.sty), review to appear in "Cosmic Explosions", proceedings. of the 10th Annual October Maryland Astrophysics Conference, eds. S. Holt and W. W. Zhang, AIP

Classical novae eject significant amounts of nuclear processed material into the interstellar medium. Among the isotopes synthesized during such explosions, two radioactive nuclei deserve a particular attention: 22Na and 26Al. In this paper, we investigate the nuclear paths leading to 22Na and 26Al production during nova outbursts by means of an implicit, hydrodynamic code that follows the course of the thermonuclear runaway from the onset of accretion up to the ejection stage. New evolutionary sequences of ONe novae have been computed, using updated nuclear reaction rates relevant to 22Na and 26Al production. Special attention is focused on the role played by nuclear uncertainties within the NeNa and MgAl cycles in the synthesis of such radioactive species. From the series of hydrodynamic models, which assume upper, recommended or lower estimates of the reaction rates, we derive limits on the production of both 22Na and 26Al. We outline a list of nuclear reactions which deserve new experimental investigations in order to reduce the wide dispersion introduced by nuclear uncertainties in the 22Na and 26Al yields., Comment: 46 pages, 4 figures. Accepted for publication in The Astrophysical Journal

Classical novae eject significant amounts of matter into the interstellar medium, as a result of thermonuclear runaways. Nucleosynthesis associated with nova outbursts includes products from explosive H burning, such as 17O, 15N and 13C, and also radioactive species like 7Be, 22Na and 26Al. In this paper we report on new hydrodynamic calculations of nova outbursts, from the onset of accretion to mass ejection. We stress the role played by the nuclear uncertainties associated with key reactions of the NeNa-MgAl cycles on the synthesis of 22Na and 26Al., Comment: 4 pages (to appear in "Nuclei in the Cosmos V" Conference Proceedings)

We investigate the evolution of cooling helium atmosphere white dwarfs using a full evolutionary code, specifically developed for following the effects of element diffusion and gravitational settling on white dwarf cooling. The major difference between this work and previous work is that we use more recent opacity data from the OPAL project. Since, in general, these opacities are higher than those available ten years ago, at a given effective temperature, convection zones go deeper than in models with older opacity data. Thus convective dredge-up of observationally detectable carbon in helium atmosphere white dwarfs can occur for thicker helium layers than found by Pelletier et al (1986). We find that the range of observed C to He ratios in different DQ white dwarfs of similar effective temperature is well explained by a range of initial helium layer mass between $10^{-3}$ and $10^{-2} M_{\odot}$, in good agreement with stellar evolution theory, assuming a typical white dwarf mass of $0.6 M_{\odot}$. We also predict that oxygen will be present in DQ white dwarf atmospheres in detectable amounts if the helium layer mass is near the lower limit compatible with stellar evolution theory. Determination of the oxygen abundance has the potential of providing information on the profile of oxygen in the core and hence on the important $^{12}$C$(\alpha,\gamma)^{16}$O reaction rate., Comment: Accepted for publication in MNRAS, 8 pages, 6 figures

Detailed nucleosynthesis in the ejecta of classical novae has been determined for a grid of hydrodynamic nova models. The reported 14 evolutionary sequences, followed from the onset of accretion up to the explosion and ejection stages, span a range of CO and ONe white dwarf masses (0.8-1.35 $M_o$) and mixing levels between the accreted envelope and the underlying white dwarf core (25-75%). The synthesis of each isotope, from ${}^1$H to ${}^{40}$Ca, is discussed along with its sensitivity to model parameters. Special emphasis is given to isotopes such as ${}^{13}$C, ${}^{15}$N, and ${}^{17}$O, which may account for a significant fraction of their Galactic content. Production of the radioactive isotopes ${}^7$Be, ${}^{22}$Na, and ${}^{26}$Al is also analyzed, since they may provide a direct test of the thermonuclear runaway model through their $\gamma$-ray emission. The resulting elemental yields reproduce fairly well the spectroscopic abundance determinations of several well studied classical novae., Comment: 38 pages, Latex file, with 3 ps figures. Accepted for publication in The Astrophysical Journal

The production of 26Al by explosions of classical novae has been computed by means of a hydrodynamic code that follows both the accretion and the explosion stages. A special emphasis has been put on the analysis of the influence of the initial abundances of the accreted envelope, as well as on the nuclear reaction rates involved. With the most recent values of chemical composition and reaction rates available, 26Al production is lowered with respect to previous computations. According to our results, the final contribution of novae to the galactic 26Al is at most 0.4 M_sun, which is a small part of the estimated 26Al in the Galaxy derived from COMPTEL observations of the 1809 keV emission., Comment: Accepted for publication in ApJ Letters (16 pages, 1 figure, aasms4)