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.

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.

This is a White Paper in support of the mission concept of the Large Observatory for X-ray Timing (LOFT), proposed as a medium-sized ESA mission. We discuss the potential of LOFT for the study of thermonuclear X-ray bursts on accreting neutron stars. For a summary, we refer to the paper.

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.

Classical novae are expected to contribute to the 1809-keV Galactic $\gamma$-ray emission by producing its precursor $^{26}$Al, but the yield depends on the thermonuclear rate of the unmeasured $^{25}$Al($p,\gamma$)$^{26}$Si reaction. Using the $\beta$ decay of $^{26}$P to populate the key $J^{\pi}=3^+$ resonance in this reaction, we report the first evidence for the observation of its exit channel via a $1741.6 \pm 0.6 (\textrm{stat}) \pm 0.3 (\textrm{syst})$ keV primary $\gamma$ ray, where the uncertainties are statistical and systematic, respectively. By combining the measured $\gamma$-ray energy and intensity with other experimental data on $^{26}$Si, we find the center-of-mass energy and strength of the resonance to be $E_r = 414.9 \pm 0.6(\textrm{stat}) \pm 0.3 (\textrm{syst}) \pm 0.6(\textrm{lit.})$ keV and $\omega\gamma = 23 \pm 6 (\textrm{stat})^{+11}_{-10}(\textrm{lit.})$ meV, respectively, where the last uncertainties are from adopted literature data. We use hydrodynamic nova simulations to model $^{26}$Al production showing that these measurements effectively eliminate the dominant experimental nuclear-physics uncertainty and we estimate that novae may contribute up to 30% of the Galactic $^{26}$Al.

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.

The 33S(p,\gamma)34Cl reaction is important for constraining predictions of certain isotopic abundances in oxygen-neon novae. Models currently predict as much as 150 times the solar abundance of 33S in oxygen-neon nova ejecta. This overproduction factor may, however, vary by orders of magnitude due to uncertainties in the 33S(p,\gamma)34Cl reaction rate at nova peak temperatures. Depending on this rate, 33S could potentially be used as a diagnostic tool for classifying certain types of presolar grains. Better knowledge of the 33S(p,\gamma)34Cl rate would also aid in interpreting nova observations over the S-Ca mass region and contribute to the firm establishment of the maximum endpoint of nova nucleosynthesis. Additionally, the total S elemental abundance which is affected by this reaction has been proposed as a thermometer to study the peak temperatures of novae. Previously, the 33S(p,\gamma)34Cl reaction rate had only been studied directly down to resonance energies of 432 keV. However, for nova peak temperatures of 0.2-0.4 GK there are 7 known states in 34Cl both below the 432 keV resonance and within the Gamow window that could play a dominant role. Direct measurements of the resonance strengths of these states were performed using the DRAGON recoil separator at TRIUMF. Additionally two new states within this energy region are reported. Several hydrodynamic simulations have been performed, using all available experimental information for the 33S(p,\gamma)34Cl rate, to explore the impact of the remaining uncertainty in this rate on nucleosynthesis in nova explosions. These calculations give a range of ~ 20-150 for the expected 33S overproduction factor, and a range of ~ 100-450 for the 32S/33S ratio expected in ONe novae.

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.

The 18F(p,\alpha)15O reaction rate is crucial for constraining model predictions of the \gamma-ray observable radioisotope 18F produced in novae. The determination of this rate is challenging due to particular features of the level scheme of the compound nucleus, 19Ne, which result in interference effects potentially playing a significant role. The dominant uncertainty in this rate arises from interference between J\pi=3/2+ states near the proton threshold (Sp = 6.411 MeV) and a broad J\pi=3/2+ state at 665 keV above threshold. This unknown interference term results in up to a factor of 40 uncertainty in the astrophysical S-factor at nova temperatures. Here we report a new measurement of states in this energy region using the 19F(3He,t)19Ne reaction. In stark contrast with previous assumptions we find at least 3 resonances between the proton threshold and Ecm=50 keV, all with different angular distributions. None of these are consistent with J\pi= 3/2+ angular distributions. We find that the main uncertainty now arises from the unknown proton-width of the 48 keV resonance, not from possible interference effects. Hydrodynamic nova model calculations performed indicate that this unknown width affects 18F production by at least a factor of two in the model considered.

Type I X-ray bursts are thermonuclear explosions that occur in the envelopes of accreting neutron stars. Detailed observations of these phenomena have prompted numerous studies in theoretical astrophysics and experimental nuclear physics since their discovery over 35 years ago. In this review, we begin by discussing key observational features of these phenomena that may be sensitive to the particular patterns of nucleosynthesis from the associated thermonuclear burning. We then summarize efforts to model type I X-ray bursts, with emphasis on determining the nuclear physics processes involved throughout these bursts. We discuss and evaluate limitations in the models, particularly with regard to key uncertainties in the nuclear physics input. Finally, we examine recent, relevant experimental measurements and outline future prospects to improve our understanding of these unique environments from observational, theoretical and experimental perspectives.

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).

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.

We present a nucleosynthesis study of the merger of a 0.4 solar masses helium white dwarf with a 0.8 solar masses carbon-oxygen white dwarf, coupling the thermodynamic history of Smoothed Particle Hydrodynamics particles with a post-processing code. The resulting chemical abundance pattern, particularly for oxygen and fluorine, is in qualitative agreement with the observed abundances in R Coronae Borealis stars.

Many hydrogen deficient stars are characterised by surface abundance patterns that are hard to reconcile with conventional stellar evolution. Instead, it has been suggested that they may represent the result of a merger episode between a helium and a carbon-oxygen white dwarf. In this Letter, we present a nucleosynthesis study of the merger of a 0.4 M_sol helium white dwarf with a 0.8 M_sol carbon-oxygen white dwarf, by coupling the thermodynamic history of Smoothed Particle Hydrodynamics particles with a post-processing code. The resulting chemical abundance pattern, particularly for oxygen and fluorine, is in qualitative agreement with the observed abundances in R Coronae Borealis stars.

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.

The recurrent nova RS Oph experienced an outburst in 2006, 21 years after its previous explosion in 1985, as expected. It was observed at almost all wavelengths, and important information about its properties is still being extracted. We present theoretical models of the explosion of this fascinating object, which indicate that the mass of the accreting white dwarf should be very close to the Chandrasekhar mass, to allow for such a short recurrence period. In addition, since models suggest that this nova ejects less mass than it accretes, it is an excellent candidate for a thermonuclear supernova explosion, in about $10^5-10^7$ years from now. We also analyze the emission of soft gamma-rays by RS Oph detected with the BAT instrument onboard Swift, and with the PCA onboard RXTE. We rule out that this emission has its origin in radioactive decays in the expanding nova envelope.

The radioactive isotope $^{22}$Na is one of the primary tracer isotopes that may potentially be detectable by gamma-ray spectroscopy. This isotope is predicted to be produced in ONe nova explosions, yet no detection of its 1.275 MeV gamma-ray signature has been reported so far. In particular, 6 years of COMPTEL observations of the galactic bulge, a region that should be particularily rich in nova explosions, have only led to an upper-limit. In this paper we will present the implications of the COMPTEL upper flux limit on the permissible $^{22}$Na yield for ONe novae. We base our analysis on extensive Monte-Carlo simulations of galactic 1.275 MeV emission that were constrained by the most recent results of galactic nova rates and distributions. We demonstrate that the non-detection of the 1.275 MeV line from the galactic bulge implies a solid upper $^{22}$Na yield limit in agreement with current theoretical nucleosynthesis calculations.

Classical novae are important producers of radioactive nuclei, such as be7, n13, f18, na22 and al26. The disintegration of these nuclei produces positrons (except for be7) that through annihilation with electrons produce photons of energies 511 keV and below. Furthermore, be7 and na22 decay producing photons with energies of 478 keV and 1275 keV, respectively, well in the gamma-ray domain. Therefore, novae are potential sources of gamma-ray emission. The properties of gamma-ray spectra and gamma-ray light curves (for the continuum and for the lines at 511, 478 and 1275 keV) have been analyzed, with a special emphasis on the difference between carbon-oxygen and oxygen-neon novae. Predictions of detectability of individual novae by the future SPI spectrometer on board the INTEGRAL satellite are made.

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.