Context. Monte Carlo methods can be used to evaluate the uncertainty of a reaction rate that arises from many uncertain nuclear inputs. However, until now no attempt has been made to find the effect of correlated energy uncertainties in input resonance parameters. Aims. To investigate the impact of correlated energy uncertainties on reaction rates. Methods. Using a combination of numerical and Monte Carlo variation of resonance energies, the effect of correlations are investigated. Five reactions are considered: two fictional, illustrative cases and three reactions whose rates are of current interest. Results. The effect of correlations in resonance energies depends on the specific reaction cross section and temperatures considered. When several resonances contribute equally to a reaction rate, and are located either side of the Gamow peak, correlations between their energies dilute their effect on reaction rate uncertainties. If they are both located above or below the maximum of the Gamow peak, however, correlations between their resonance energies can increase the reaction rate uncertainties. This effect can be hard to predict for complex reactions with wide and narrow resonances contributing to the reaction rate.

It has been suggested that hydrogen ingestion into the helium shell of massive stars could lead to high $^{13}$C and $^{15}$N excesses when the shock of a core-collapse supernova passes through its helium shell. This prediction questions the origin of extremely high $^{13}$C and $^{15}$N abundances observed in rare presolar SiC grains which is usually attributed to classical novae. In this context $^{13}$N($\alpha$,p)$^{16}$O the reaction plays an important role since it is in competition with $^{13}$N $\beta^+$-decay to $^{13}$C. The $^{13}$N($\alpha$,p)$^{16}$O reaction rate used in stellar evolution calculations comes from the CF88 compilation with very scarce information on the origin of this rate. The goal of this work is to provide a recommended $^{13}$N($\alpha$,p)$^{16}$O reaction rate, based on available experimental data. Unbound nuclear states in the $^{17}$F compound nucleus were studied using the spectroscopic information of the analog states in $^{17}$O nucleus that were measured at the Alto facility using the $^{13}$C($^7$Li,t)$^{17}$O alpha-transfer reaction, and spectroscopic factors were derived using a DWBA analysis. This spectroscopic information was used to calculate a recommended $^{13}$N($\alpha$,p)$^{16}$O reaction rate with meaningful uncertainty using a Monte Carlo approach. The present $^{13}$N($\alpha$,p)$^{16}$O reaction rate is found to be within a factor of two of the previous evaluation, with a typical uncertainty of a factor 2-3. The source of this uncertainty comes from the three resonances at $E_r^{c.m.} = 221$, 741 and 959 keV. This new error estimation translates to an overall uncertainty in the $^{13}$C production of a factor of 50. The main source of uncertainty on the re-evaluated $^{13}$N($\alpha$,p)$^{16}$O reaction rate currently comes from the uncertain alpha-width of relevant $^{17}$F states.

The rate of the $^{25}$Al($p$,$\gamma$)$^{26}$Si reaction is one of the few key remaining nuclear uncertainties required for predicting the production of the cosmic $\gamma$-ray emitter $^{26}$Al in explosive burning in novae. This reaction rate is dominated by three key resonances ($J^{\pi}=0^{+}$, $1^{+}$ and $3^{+}$) in $^{26}$Si. Only the $3^{+}$ resonance strength has been directly constrained by experiment. A high resolution measurement of the $^{25}$Mg($d$,$p$) reaction was used to determine spectroscopic factors for analog states in the mirror nucleus, $^{26}$Mg. A first spectroscopic factor value is reported for the $0^{+}$ state at 6.256 MeV, and a strict upper limit is set on the value for the $1^{+}$ state at 5.691 MeV, that is incompatible with an earlier ($^{4}$He,$^{3}$He) study. These results are used to estimate proton partial widths, and resonance strengths of analog states in $^{26}$Si contributing to the $^{25}$Al($p$,$\gamma$)$^{26}$Si reaction rate in nova burning conditions.

Presolar grains originating in oxygen-neon novae may be identified by their sulfur isotopic ratios compared with theoretical estimates. These ratios depend on reliable $^{33}$S($p, \gamma$)$^{34}$Cl and $^{34}$S($p, \gamma$)$^{35}$Cl reaction rates. The latter rate has recently been computed based on experimental input, and many new excited states in $^{35}$Cl were discovered above the proton threshold. The experimental $^{34}$S($p, \gamma$)$^{35}$Cl rate was found to be 2 - 5 times smaller than the theoretical one, and the simulated $^{34}$S/$^{32}$S isotopic ratio for nova presolar grains was thus predicted to be smaller than that of type II supernova grains by up to a factor of 3.7. The present study was performed to confirm the existence of these new resonances, and to improve the remaining uncertainties in the $^{34}$S($p, \gamma$)$^{35}$Cl reaction rate. Energies and spin-parities of the $^{35}$Cl excited levels were investigated with an Enge split-pole spectrograph using the $^{32}$S($\alpha, p$)$^{35}$Cl reaction. Differential cross sections of the outgoing protons were measured at $E_{\alpha}$ = 21 MeV. The existence of the newly discovered states are largely confirmed, although a few states were not observed in this study. The spins and parities of several $^{35}$Cl states were assigned tentatively for the first time. The present $^{34}$S($p, \gamma$)$^{35}$Cl experimental thermonuclear reaction rate is consistent within 1$\sigma$ with the previous evaluation. However, our rate uncertainty is larger due to a more realistic treatment of the experimental uncertainties. The uncertainty in the present rate is up to a factor of 3.5 at nova temperatures. We recommend future work to focus on the unknown properties of four excited states of $^{35}$Cl at 6643 keV, 6761 keV, 6780 keV, and 6800 keV.

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.

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.

The magnesium-potassium anti-correlation observed in globular cluster NGC2419 can be explained by nuclear burning of hydrogen in hot environments. The exact site of this nuclear burning is, as yet, unknown. In order to constrain the sites responsible for this anti-correlation, the nuclear reactions involved must be well understood. The $^{39}$K+p reactions are one such pair of reactions. Here, we report a new evaluation of the $^{39}$K(p,$\gamma$)$^{40}$Ca reaction rate by taking into account ambiguities and measurement uncertainties in the nuclear data. The uncertainty in the $^{39}$K(p,$\gamma$)$^{40}$Ca reaction rate is larger than previously assumed, and its influence on nucleosynthesis models is demonstrated. We find the $^{39}$K(p,$\gamma$)$^{40}$Ca reaction cross section should be the focus of future experimental study to help constrain models aimed at explaining the magnesium-potassium anti-correlation in globular clusters.

Context. Monte Carlo methods have enabled nuclear reaction rates from uncertain inputs to be presented in a statistically meaningful manner. However, these uncertainties are currently computed assuming no correlations between the physical quantities that enter those calculations. This is not always an appropriate assumption. Astrophysically important reactions are often dominated by resonances, whose properties are usually normalized to a well-known reference resonance. This insight provides a basis from which to develop a flexible framework for including correlations in Monte Carlo reaction rate calculations. Aims. The aim of this work is to develop and test a method for including correlations in Monte Carlo reaction rate calculations when the input has been normalized to a common reference. Methods. A mathematical framework is developed for including correlations between input parameters in Monte Carlo reaction rate calculations. The magnitude of those correlations is calculated from the uncertainties typically reported in experimental papers, where full correlation information is not available. The method is applied to four illustrative examples: a fictional 3-resonance reaction, $^{27}$Al(p,$\gamma$)$^{28}$Si, $^{23}$Na(p,$\alpha$)$^{20}$Ne, and $^{23}$Na($\alpha$,p)$^{26}$Al. Results. Reaction rates at low temperatures that are dominated by a few isolated resonances are found to minimally impacted by correlation effects. However, reaction rates determined from many overlapping resonances can be significantly affected. Uncertainties in the $^{23}$Na($\alpha$,p)$^{26}$Al reaction, for example, increase by up to a factor of 5. This highlights the need to take correlation effects into account in reaction rate calculations, and provide insight into which cases are expected to be most affected by them. The impact of correlation effects on nucleosynthesis is also investigated.

Primordial or big bang nucleosynthesis (BBN) is one of the three historical strong evidences for the big bang model. Standard BBN is now a parameter free theory, since the baryonic density of the Universe has been deduced with an unprecedented precision from observations of the anisotropies of the cosmic microwave background (CMB) radiation. There is a good agreement between the primordial abundances of 4He, D, 3He and 7Li deduced from observations and from primordial nucleosynthesis calculations. However, the 7Li calculated abundance is significantly higher than the one deduced from spectroscopic observations and remains an open problem. In addition, recent deuterium observations have drastically reduced the uncertainty on D/H, to reach a value of 1.6%. It needs to be matched by BBN predictions whose precision is now limited by thermonuclear reaction rate uncertainties. This is especially important as many attempts to reconcile Li observations with models lead to an increased D prediction. Here, we re-evaluates the D(p,g)3He, D(d,n)3He and D(d,p)3H reaction rates that govern deuterium destruction, incorporating new experimental data and carefully accounting for systematic uncertainties. Contrary to previous evaluations, we use theoretical ab initio models for the energy dependence of the S-factors. As a result, these rates increase at BBN temperatures, leading to a reduced value of D/H = (2.45$\pm0.10)\times10^{-5}$ (2$\sigma$), in agreement with observations.

The $^{18}$Ne($\alpha$,$p$)$^{21}$Na reaction impacts the break-out from the hot CNO-cycles to the $rp$-process in type I X-ray bursts. We present a revised thermonuclear reaction rate, which is based on the latest experimental data. The new rate is derived from Monte-Carlo calculations, taking into account the uncertainties of all nuclear physics input quantities. In addition, we present the reaction rate uncertainty and probability density versus temperature. Our results are also consistent with estimates obtained using different indirect approaches.

Rigorous statistical methods for estimating thermonuclear reaction rates and nucleosynthesis are becoming increasingly established in nuclear astrophysics. The main challenge being faced is that experimental reaction rates are highly complex quantities derived from a multitude of different measured nuclear parameters (e.g., astrophysical S-factors, resonance energies and strengths, particle and gamma-ray partial widths). We discuss the application of the Monte Carlo method to two distinct, but related, questions. First, given a set of measured nuclear parameters, how can one best estimate the resulting thermonuclear reaction rates and associated uncertainties? Second, given a set of appropriate reaction rates, how can one best estimate the abundances from nucleosynthesis (i.e., reaction network) calculations? The techniques described here provide probability density functions that can be used to derive statistically meaningful reaction rates and final abundances for any desired coverage probability. Examples are given for applications to s-process neutron sources, core-collapse supernovae, classical novae, and big bang nucleosynthesis.

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.

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.

Context: Recent reaction rate evaluations include reaction rate uncertainties that have been determined in a statistically meaningful manner. Furthermore, reaction rate probability density distributions have been determined and published in the form of lognormal parameters with the specific goal of pursuing Monte Carlo nucleosynthesis studies. Aims: To test and assess different methods of randomly sampling over reaction rate probability densities and to determine the most accurate method for estimating elemental abundance uncertainties. Methods: Experimental Monte Carlo reaction rates are first computed for the 22Ne+alpha, 20Ne(p,g)21Na, 25Mg(p,g)26Al, and 18F(p,alpha)15O reactions, which are used to calculate reference nucleosynthesis yields for 16 nuclei affected by nucleosynthesis in massive stars and classical novae. Five different methods of randomly sampling over these reaction rate probability distributions are then developed, tested, and compared with the reference nucleosynthesis yields. Results: Given that the reaction rate probability density distributions can be described accurately with a lognormal distribution, Monte Carlo nucleosynthesis variations arising from the parametrised estimates for the reaction rate variations agree remarkably well with those obtained from the true rate samples. Most significantly, the most simple parametrisation agrees within just a few percent, meaning that Monte Carlo nucleosynthesis studies can be performed reliably using lognormal parametrisations of reaction rate probability density functions.

The 22Ne({\alpha},n)25Mg reaction is an important source of neutrons for the s-process. In massive stars responsible for the weak component of the s-process, 22Ne({\alpha},n)25Mg is the dominant source of neutrons, both during core helium burning and in shell carbon burning. For the main s-process component produced in Asymptotic Giant Branch (AGB) stars, the 13C({\alpha},n)16O reaction is the dominant source of neutrons operating during the interpulse period, with the 22Ne+{\alpha} source affecting mainly the s-process branchings during a thermal pulse. Rate uncertainties in the competing 22Ne({\alpha},n)25Mg and 22Ne({\alpha},{\gamma})26Mg reactions result in large variations of s-process nucleosynthesis. Here, we present up-to-date and statistically rigorous 22Ne+{\alpha} reaction rates using recent experimental results and Monte Carlo sampling. Our new rates are used in post-processing nucleosynthesis calculations both for massive stars and AGB stars. We demonstrate that the nucleosynthesis uncertainties arising from the new rates are dramatically reduced in comparison to previously published results, but several ambiguities in the present data must still be addressed. Recommendations for further study to resolve these issues are provided.

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.

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.

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.

A method based on Monte Carlo techniques is presented for evaluating thermonuclear reaction rates. We begin by reviewing commonly applied procedures and point out that reaction rates that have been reported up to now in the literature have no rigorous statistical meaning. Subsequently, we associate each nuclear physics quantity entering in the calculation of reaction rates with a specific probability density function, including Gaussian, lognormal and chi-squared distributions. Based on these probability density functions the total reaction rate is randomly sampled many times until the required statistical precision is achieved. This procedure results in a median (Monte Carlo) rate which agrees under certain conditions with the commonly reported recommended "classical" rate. In addition, we present at each temperature a low rate and a high rate, corresponding to the 0.16 and 0.84 quantiles of the cumulative reaction rate distribution. These quantities are in general different from the statistically meaningless "minimum" (or "lower limit") and "maximum" (or "upper limit") reaction rates which are commonly reported. Furthermore, we approximate the output reaction rate probability density function by a lognormal distribution and present, at each temperature, the lognormal parameters miu and sigma. The values of these quantities will be crucial for future Monte Carlo nucleosynthesis studies. Our new reaction rates, appropriate for bare nuclei in the laboratory, are tabulated in the second paper of this series (Paper II). The nuclear physics input used to derive our reaction rates is presented in the third paper of this series (Paper III). In the fourth paper of this series (Paper IV) we compare our new reaction rates to previous results.

Numerical values of charged-particle thermonuclear reaction rates for nuclei in the A=14 to 40 region are tabulated. The results are obtained using a method, based on Monte Carlo techniques, that has been described in the preceding paper of this series (Paper I). We present a low rate, median rate and high rate which correspond to the 0.16, 0.50 and 0.84 quantiles, respectively, of the cumulative reaction rate distribution. The meaning of these quantities is in general different from the commonly reported, but statistically meaningless expressions, "lower limit", "nominal value" and "upper limit" of the total reaction rate. In addition, we approximate the Monte Carlo probability density function of the total reaction rate by a lognormal distribution and tabulate the lognormal parameters {\mu} and {\sigma} at each temperature. We also provide a quantitative measure (Anderson-Darling test statistic) for the reliability of the lognormal approximation. The user can implement the approximate lognormal reaction rate probability density functions directly in a stellar model code for studies of stellar energy generation and nucleosynthesis. For each reaction, the Monte Carlo reaction rate probability density functions, together with their lognormal approximations, are displayed graphically for selected temperatures in order to provide a visual impression. Our new reaction rates are appropriate for bare nuclei in the laboratory. The nuclear physics input used to derive our reaction rates is presented in the subsequent paper of this series (Paper III). In the fourth paper of this series (Paper IV) we compare our new reaction rates to previous results.

The nuclear physics input used to compute the Monte Carlo reaction rates and probability density functions that are tabulated in the second paper of this series (Paper II) is presented. Specifically, we publish the input files to the Monte Carlo reaction rate code RatesMC, which is based on the formalism presented in the first paper of this series (Paper I). This data base contains overwhelmingly experimental nuclear physics information. The survey of literature for this review was concluded in November 2009.

We compare our Monte Carlo reaction rates (see Paper II of this series) to previous results that were obtained by using the classical method of computing thermonuclear reaction rates. For each reaction, the comparison is presented using two types of graphs: the first shows the change in reaction rate uncertainties, while the second displays our new results normalized to the previously recommended reaction rate. We find that the rates have changed significantly for almost all reactions considered here. The changes are caused by (i) our new Monte Carlo method of computing reaction rates (see Paper I of this series), and (ii) newly available nuclear physics information (see Paper III of this series).