Cataclysmic Variables (CVs) and Symbiotic Binaries are close (or not so close) binary star systems which contain both a white dwarf (WD) primary and a larger cooler secondary star that typically fills its Roche Lobe. The cooler star is losing mass through the inner Lagrangian point of the binary and a fraction of this material is accreted by the WD. Here we report on our hydrodynamic studies of the thermonuclear runaway (TNR) in the accreted material that ends in a Classical Nova explosion. We have followed the evolution of the TNRs on both carbon-oxygen (CO) and oxygen-neon (ONe) WDs. We report on 3 studies in this paper. First, simulations in which we accrete only solar matter using NOVA (our 1-D, fully implicit, hydro code). Second, we use MESA for similar studies in which we accrete only Solar matter and compare the results. Third, we accrete solar matter until the TNR is ongoing and then switch the composition in the accreted layers to a mixed composition: either 25% WD and 75% solar or 50% WD and 50% Solar matter. We find that the amount of accreted material is inversely proportional to the initial 12C abundance (as expected). Thus, accreting solar matter results in a larger amount of accreted material to fuel the outburst; much larger than in earlier studies where a mixed composition was assumed from the beginning of the simulation. Our most important result is that all these simulations eject significantly less mass than accreted and, therefore, the WD is growing in mass toward the Chandrasekhar Limit.

Big bang nucleosynthesis provides the earliest probe of standard model physics, at a time when the universe was less than a thousand seconds old. It determines the abundances of the lightest nuclides, which give rise to the subsequent history of the visible matter in the Universe. This work derives new $^7$Be(n,p)$^7$Li thermonuclear reaction rates based on all available experimental information. This reaction sensitively impacts the primordial abundances of $^{7}$Be and $^7$Li during big bang nucleosynthesis. We critically evaluate all available data and disregard experimental results that are questionable. For the nuclear model, we adopt an incoherent sum of single-level, two-channel R-matrix approximation expressions, which are implemented into a hierarchical Bayesian model, to analyze the remaining six data sets we deem most reliable. In the fitting of the data, we consistently model all known sources of uncertainty, including discrepant absolute normalizations of different data sets, and also take the variation of the neutron and proton channel radii into account, hence providing less biased estimates of the $^7$Be(n,p)$^7$Li thermonuclear rates. From the resulting posteriors, we extract R-matrix parameters ($E_r$, $\gamma^2_n$, $\gamma^2_p$) and derive excitation energies, partial and total widths. Our fit is sensitive to the contributions of the first three levels above the neutron threshold. Reaction rates were computed by integrating 10,000 samples of the reduced cross section. Our $^7$Be(n,p)$^7$Li thermonuclear rates have uncertainties between 1.5% and 2.0% at temperatures of $\leq$1 GK. We compare our rates to previous results and find that the $^7$Be(n,p)$^7$Li rates most commonly used in big bang simulations have too optimistic uncertainties.

We report on studies of Classical Nova (CN) explosions where we follow the evolution of thermonuclear runaways (TNRs) on Carbon Oxygen (CO) white dwarfs (WDs). We vary both the mass of the WD (from 0.6 M$_\odot$ to 1.35 M$_\odot$) and the composition of the accreted material. Our simulations are guided by the results of multi-dimensional studies of TNRs in WDs that find sufficient mixing with WD core material occurs after the TNR is well underway, reaching levels of enrichment that agree with observations of CN ejecta abundances. We use NOVA (our 1-dimensional hydrodynamic code) to accrete solar matter until the TNR is ongoing and then switch to a mixed composition (either 25% WD material and 75% solar or 50% WD material and 50% solar). Because the amount of accreted material is inversely proportional to the initial $^{12}$C abundance, by first accreting solar matter the amount of material taking part in the outburst is larger than in those simulations where we assume a mixed composition from the beginning. Our results show large enrichments of $^7$Be in the ejected gases implying that CO CNe may be responsible for a significant fraction ($\sim$ 100 M$_\odot$) of the $^7$Li in the galaxy ($\sim$1000 M$_\odot$). In addition, although the ejected gases are enriched in WD material, the WDs in these simulations eject less material than they accrete. We predict that the WD is growing in mass as a consequence of the accretion-outburst-accretion cycle and CO CNe may be an important channel of Supernova Ia progenitors.

Resonance energies are frequently derived from precisely measured excitation energies and reaction Q-values. The latter quantities are usually calculated from atomic instead of nuclear mass differences. This procedure disregards the energy shift caused by the difference in the total electron binding energies before and after the interaction. Assuming that the interacting nuclei in a stellar plasma are fully ionized, this energy shift can have a significant effect, considering that the resonance energy enters exponentially into the expression for the narrow-resonance thermonuclear reaction rates. As an example, the rate of the $^{36}$Ar(p,$\gamma$)$^{37}$K reaction is discussed, which, at temperatures below 1 GK, depends only on the contributions of a single resonance and direct capture. In this case, disregarding the energy shift caused by the total electron binding energy difference erroneously enhances the rate by $\approx$40\% near temperatures of 70 MK.

We developed a hierarchical Bayesian framework to estimate S-factors and thermonuclear rates for the $^3$He(d,p)$^{4}$He reaction, which impacts the primordial abundances of $^3$He and $^7$Li. The available data are evaluated and all direct measurements are taken into account in our analysis for which we can estimate separate uncertainties for systematic and statistical effects. For the nuclear reaction model, we adopt a single-level, two-channel approximation of R-matrix theory, suitably modified to take the effects of electron screening at lower energies into account. Apart from the usual resonance parameters (resonance location and reduced widths for the incoming and outgoing reaction channel), we include for the first time the channel radii and boundary condition in the fitting process. Our new analysis of the $^3$He(d,p)$^{4}$He S-factor data results in improved estimates for the thermonuclear rates. This work represents the first nuclear rate evaluation using the R-matrix theory embedded into a hierarchical Bayesian framework, properly accounting for all known sources of uncertainty. Therefore, it provides a test bed for future studies of more complex reactions.

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.

The $^3$H(d,n)$^4$He reaction has a large low-energy cross section and will likely be utilized in future commercial fusion reactors. This reaction also takes place during big bang nucleosynthesis. Studies of both scenarios require accurate and precise fusion rates. To this end, we implement a one-level, two-channel R-matrix approximation into a Bayesian model. Our main goals are to predict reliable astrophysical S-factors and to estimate R-matrix parameters using the Bayesian approach. All relevant parameters are sampled in our study, including the channel radii, boundary condition parameters, and data set normalization factors. In addition, we take uncertainties in both measured bombarding energies and S-factors rigorously into account. Thermonuclear rates and reactivities of the $^3$H(d,n)$^4$He reaction are derived by numerically integrating the Bayesian S-factor samples. The present reaction rate uncertainties at temperatures between $1.0$ MK and $1.0$ GK are in the range of 0.2% to 0.6%. Our reaction rates differ from previous results by 2.9% near 1.0 GK. Our reactivities are smaller than previous results, with a maximum deviation of 2.9% near a thermal energy of $4$ keV. The present rate or reactivity uncertainties are more reliable compared to previous studies that did not include the channel radii, boundary condition parameters, and data set normalization factors in the fitting. Finally, we investigate previous claims of electron screening effects in the published $^3$H(d,n)$^4$He data. No such effects are evident and only an upper limit for the electron screening potential can be obtained.

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.

The study of d+d reactions is of major interest since their reaction rates affect the predicted abundances of D, 3He, and 7Li. In particular, recent measurements of primordial D/H ratios call for reduced uncertainties in the theoretical abundances predicted by big bang nucleosynthesis (BBN). Different authors have studied reactions involved in BBN by incorporating new experimental data and a careful treatment of systematic and probabilistic uncertainties. To analyze the experimental data, Coc et al. (2015) used results of ab initio models for the theoretical calculation of the energy dependence of S-factors in conjunction with traditional statistical methods based on Chi-2 minimization. Bayesian methods have now spread to many scientific fields and provide numerous advantages in data analysis. Astrophysical S-factors and reaction rates using Bayesian statistics were calculated by Iliadis et al. (2016). Here we present a similar analysis for two d+d reactions, d(d,n)3He and d(d,p)3H, that has been translated into a total decrease of the predicted D/H value by 0.16%.

Abundance anomalies in globular clusters provide strong evidence for multiple stellar populations within each cluster. These populations are usually interpreted as distinct generations, with the currently observed second-generation stars having formed in part from the ejecta of massive, first-generation "polluter" stars, giving rise to the anomalous abundance patterns. The precise nature of the polluters and their enrichment mechanism are still unclear. Even so, the chemical abundances measured in second-generation stars within the globular cluster NGC 2419 provide insight into this puzzling process. Previous work used Monte Carlo nuclear reaction network calculations to constrain the temperature-density conditions that could reproduce the observed abundances, thereby placing robust limits on the origins of the polluter material. The effect of individual reaction rates on these conditions has not been studied, however. Thus, we perform an exhaustive sensitivity study on the nuclear physics input to determine which reactions have the greatest impact on these predictions. We find that the $^{30}$Si(p,$\gamma$)$^{31}$P, $^{37}$Ar(p,$\gamma$)$^{38}$K, $^{38}$Ar(p,$\gamma$)$^{39}$K, and $^{39}$K(p,$\gamma$)$^{40}$Ca reactions are all critical in determining the temperature-density conditions, and ultimately, the origin of the polluter material. We conclude with recommendations for future experiments.

Motivated by recent reports concerning the observation of limited enrichment in He but excess K in stars of globular clusters, we revisit the H-burning conditions that lead to the chemical properties of multiple stellar populations in these systems. In particular, we are interested in correlations of He and K with other elements, like O, Na, Al, Mg and Si, reported in stars of NGC 2808. We perform calculations of nucleosynthesis at constant temperature and density, exploring the temperature range of 25 to 200 10$^6$ K (25 to 200 MK), using a detailed nuclear reaction network and the most up-to-date nuclear reaction rates. We find that Mg is the most sensitive "thermometer" of hydrostatic H-burning conditions, pointing to a temperature range of 70-80 MK for NGC~2808, while He is a lesser - but not negligible - constraint. K can be produced at the levels reported for NGC 2808 at temperatures $>$180 MK and Si at T$>$80 MK. However, in the former temperature range Al and Na are totally destroyed and no correlation can be obtained, in contrast to the reported observations. None of the putative polluter sources proposed so far seem to satisfy the ensemble of nucleosynthesis constraints.

The problem of estimating non-resonant astrophysical S-factors and thermonuclear reaction rates, based on measured nuclear cross sections, is of major interest for nuclear energy generation, neutrino physics, and element synthesis. Many different methods have been applied in the past to this problem, almost all of them based on traditional statistics. Bayesian methods, on the other hand, are now in widespread use in the physical sciences. In astronomy, for example, Bayesian statistics is applied to the observation of extra-solar planets, gravitational waves, and type Ia supernovae. However, nuclear physics, in particular, has been slow to adopt Bayesian methods. We present astrophysical S-factors and reaction rates based on Bayesian statistics. We develop a framework that incorporates robust parameter estimation, systematic effects, and non-Gaussian uncertainties in a consistent manner. The method is applied to the d(p,$\gamma$)$^3$He, $^3$He($^3$He,2p)$^4$He, and $^3$He($\alpha$,$\gamma$)$^7$Be reactions, important for deuterium burning, solar neutrinos, and big bang nucleosynthesis.

Nova explosions occur on the white dwarf component of a Cataclysmic Variable binary stellar system that is accreting matter lost by its companion. When sufficient material has been accreted by the white dwarf, a thermonuclear runaway occurs and ejects material in what is observed as a Classical Nova explosion. We describe both the recent advances in our understanding of the progress of the outburst and outline some of the puzzles that are still outstanding. We report on the effects of improving both the nuclear reaction rate library and including a modern nuclear reaction network in our one-dimensional, fully implicit, hydrodynamic computer code. In addition, there has been progress in observational studies of Supernovae Ia with implications about the progenitors and we discuss that in this review.

This white paper informs the nuclear astrophysics community and funding agencies about the scientific directions and priorities of the field and provides input from this community for the 2015 Nuclear Science Long Range Plan. It summarizes the outcome of the nuclear astrophysics town meeting that was held on August 21-23, 2014 in College Station at the campus of Texas A&M University in preparation of the NSAC Nuclear Science Long Range Plan. It also reflects the outcome of an earlier town meeting of the nuclear astrophysics community organized by the Joint Institute for Nuclear Astrophysics (JINA) on October 9- 10, 2012 Detroit, Michigan, with the purpose of developing a vision for nuclear astrophysics in light of the recent NRC decadal surveys in nuclear physics (NP2010) and astronomy (ASTRO2010). The white paper is furthermore informed by the town meeting of the Association of Research at University Nuclear Accelerators (ARUNA) that took place at the University of Notre Dame on June 12-13, 2014. In summary we find that nuclear astrophysics is a modern and vibrant field addressing fundamental science questions at the intersection of nuclear physics and astrophysics. These questions relate to the origin of the elements, the nuclear engines that drive life and death of stars, and the properties of dense matter. A broad range of nuclear accelerator facilities, astronomical observatories, theory efforts, and computational capabilities are needed. With the developments outlined in this white paper, answers to long standing key questions are well within reach in the coming decade.

We investigate properties of carbon-oxygen white dwarfs with respect to the composite uncertainties in the reaction rates using the stellar evolution toolkit, Modules for Experiments in Stellar Astrophysics (MESA) and the probability density functions in the reaction rate library STARLIB. These are the first Monte Carlo stellar evolution studies that use complete stellar models. Focusing on 3 M$_{\odot}$ models evolved from the pre main-sequence to the first thermal pulse, we survey the remnant core mass, composition, and structure properties as a function of 26 STARLIB reaction rates covering hydrogen and helium burning using a Principal Component Analysis and Spearman Rank-Order Correlation. Relative to the arithmetic mean value, we find the width of the 95\% confidence interval to be $\Delta M_{{\rm 1TP}}$ $\approx$ 0.019 M$_{\odot}$ for the core mass at the first thermal pulse, $\Delta$$t_{\rm{1TP}}$ $\approx$ 12.50 Myr for the age, $\Delta \log(T_{{\rm c}}/{\rm K}) \approx$ 0.013 for the central temperature, $\Delta \log(\rho_{{\rm c}}/{\rm g \ cm}^{-3}) \approx$ 0.060 for the central density, $\Delta Y_{\rm{e,c}} \approx$ 2.6$\times$10$^{-5}$ for the central electron fraction, $\Delta X_{\rm c}(^{22}\rm{Ne}) \approx$ 5.8$\times$10$^{-4}$, $\Delta X_{\rm c}(^{12}\rm{C}) \approx$ 0.392, and $\Delta X_{\rm c}(^{16}\rm{O}) \approx$ 0.392. Uncertainties in the experimental $^{12}$C($\alpha,\gamma)^{16}\rm{O}$, triple-$\alpha$, and $^{14}$N($p,\gamma)^{15}\rm{O}$ reaction rates dominate these variations. We also consider a grid of 1 to 6 M$_{\odot}$ models evolved from the pre main-sequence to the final white dwarf to probe the sensitivity of the initial-final mass relation to experimental uncertainties in the hydrogen and helium reaction rates.

Globular clusters are of paramount importance for testing theories of stellar evolution and early galaxy formation. Strong evidence for multiple populations of stars in globular clusters derives from observed abundance anomalies. A puzzling example is the recently detected Mg-K anticorrelation in NGC 2419. We perform Monte Carlo nuclear reaction network calculations to constrain the temperature-density conditions that gave rise to the elemental abundances observed in this elusive cluster. We find a correlation between stellar temperature and density values that provide a satisfactory match between simulated and observed abundances in NGC 2419 for all relevant elements (Mg, Si, K, Ca, Sc, Ti, and V). Except at the highest densities ($\rho \gtrsim 10^8$~g/cm$^3$), the acceptable conditions range from $\approx$ $100$~MK at $\approx$ $10^8$~g/cm$^3$ to $\approx$ $200$~MK at $\approx$ $10^{-4}$~g/cm$^3$. This result accounts for uncertainties in nuclear reaction rates and variations in the assumed initial composition. We review hydrogen burning sites and find that low-mass stars, AGB stars, massive stars, or supermassive stars cannot account for the observed abundance anomalies in NGC 2419. Super-AGB stars could be viable candidates for the polluter stars if stellar model parameters can be fine-tuned to produce higher temperatures. Novae, either involving CO or ONe white dwarfs, could be interesting polluter candidates, but a current lack of low-metallicity nova models precludes firmer conclusions. We also discuss if additional constraints for the first-generation polluters can be obtained by future measurements of oxygen, or by evolving models of second-generation low-mass stars with a non-canonical initial composition.

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.

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.

We have carried out photodisintegration cross-section measurements on 86Kr using monoenergetic photon beams ranging from the neutron separation energy, S_n = 9.86 MeV, to 13 MeV. We combine our experimental 86Kr(g,n)85Kr cross section with results from our recent 86Kr(g,g') measurement below the neutron separation energy to obtain the complete nuclear dipole response of 86Kr. The new experimental information is used to predict the neutron capture cross section of 85Kr, an important branching point nucleus on the abundance flow path during s-process nucleosynthesis. Our new and more precise 85Kr(n,g)86Kr cross section allows to produce more precise predictions of the 86Kr abundance from s-process models. In particular, we find that the models of the s-process in asymptotic giant branch stars of mass < 1.5 Msun, where the 13C neutron source burns convectively rather than radiatively, represent a possible solution for the highest 86Kr/82Kr ratios observed in meteoritic stardust SiC grains.

STARLIB is a next-generation, all-purpose nuclear reaction-rate library. For the first time, this library provides the rate probability density at all temperature grid points for convenient implementation in models of stellar phenomena. The recommended rate and its associated uncertainties are also included. Currently, uncertainties are absent from all other rate libraries, and, although estimates have been attempted in previous evaluations and compilations, these are generally not based on rigorous statistical definitions. A common standard for deriving uncertainties is clearly warranted. STARLIB represents a first step in addressing this deficiency by providing a tabular, up-to-date database that supplies not only the rate and its uncertainty but also its distribution. Because a majority of rates are lognormally distributed, this allows the construction of rate probability densities from the columns of STARLIB. This structure is based on a recently suggested Monte Carlo method to calculate reaction rates, where uncertainties are rigorously defined. In STARLIB, experimental rates are supplemented with: (i) theoretical TALYS rates for reactions for which no experimental input is available, and (ii) laboratory and theoretical weak rates. STARLIB includes all types of reactions of astrophysical interest to Z = 83, such as (p,g), (p,a), (a,n), and corresponding reverse rates. Strong rates account for thermal target excitations. Here, we summarize our Monte Carlo formalism, introduce the library, compare methods of correcting rates for stellar environments, and discuss how to implement our library in Monte Carlo nucleosynthesis studies. We also present a method for accessing STARLIB on the Internet and outline updated Monte Carlo-based rates.

For stars between 0.8-8.0 solar masses, nucleosynthesis enters its final phase during the asymptotic giant branch (AGB) stage. During this evolutionary period, grain condensation occurs in the stellar atmosphere, and the star experiences significant mass loss. The production of presolar grains can often be attributed to this unique stellar environment. A subset of presolar oxide grains features dramatic 18O depletion that cannot be explained by the standard AGB star burning stages and dredge-up models. An extra mixing process, referred to as "cool bottom processing" (CBP), was proposed for low-mass AGB stars. The 18O depletion observed within certain stellar environments and within presolar grain samples may result from the 18O+p processes during CBP. We report here on a study of the 18O(p,gamma)19F reaction at low energies. Based on our new results, we found that the resonance at Er = 95 keV (lab) has a negligible affect on the reaction rate at the temperatures associated with CBP. We also determined that the direct capture S-factor is almost a factor of 2 lower than the previously recommended value at low energies. An improved thermonuclear reaction rate for 18O(p,gamma)19F is presented.

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.

Symbiotic binaries are systems containing white dwarfs (WDs) and red giants. Symbiotic novae are those systems in which thermonuclear eruptions occur on the WD components. These are to be distinguished from events driven by accretion disk instabilities analogous to dwarf novae eruptions in cataclysmic variable outbursts. Another class of symbiotic systems are those in which the WD is extremely luminous and it seems likely that quiescent nuclear burning is ongoing on the accreting WD. A fundamental question is the secular evolution of the WD. Do the repeated outbursts or quiescent burning in these accreting systems cause the WD to gain or lose mass? If it is gaining mass, can it eventually reach the Chandrasekhar Limit and become a supernova (a SN Ia if it can hide the hydrogen and helium in the system)? In order to better understand these systems, we have begun a new study of the evolution of Thermonuclear Runaways (TNRs) in the accreted envelopes of WDs using a variety of initial WD masses, luminosities and mass accretion rates. We use our 1-D hydro code, NOVA, which includes the new convective algorithm of Arnett, Meakin and Young, the Hix and Thielemann nuclear reaction solver, the Iliadis reaction rate library, the Timmes equation of state, and the OPAL opacities. We assume a solar composition (Lodders abundance distribution) and do not allow any mixing of accreted material with core material. This assumption strongly influences our results. We report here (1) that the WD grows in mass for all simulations so that canonical `steady burning' does not occur, and (2) that only a small fraction of the accreted matter is ejected in some (but not all) simulations. We also find that the accreting systems, before thermonuclear runaway, are too cool to be seen in X-ray searches for SN Ia progenitors.

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

We present a brief summary of the Single Degenerate Scenario for the progenitors of Type Ia Supernovae in which it is assumed that a low mass carbon-oxygen white dwarf is growing in mass as a result of accretion from a secondary star in a close binary system. Recent hydrodynamic simulations of accretion of solar material onto white dwarfs without mixing always produce a thermonuclear runaway and steady burning does not occur. For a broad range in WD mass (0.4 Solar masses to 1.35 Solar Masses), the maximum ejected material occurs for the 1.25 Solar Mass sequences and then decreases as the white dwarf mass decreases. Therefore, the white dwarfs are growing in mass as a consequence of the accretion of solar material and as long as there is no mixing of accreted material with core material. In contrast, a thermonuclear runaway in the accreted hydrogen-rich layers on the low luminosity WDs in close binary systems where mixing of core matter with accreted material has occurred is the outburst mechanism for Classical, Recurrent, and Symbiotic novae. The differences in characteristics of these systems is likely the WD mass and mass accretion rate. The high levels of enrichment of CN ejecta in elements ranging from carbon to sulfur confirm that there is dredge-up of matter from the core of the WD and enable them to contribute to the chemical enrichment of the interstellar medium. Therefore, studies of CNe can lead to an improved understanding of Galactic nucleosynthesis, some sources of pre-solar grains, and the Extragalactic distance scale. The characteristics of the outburst depend on the white dwarf mass, luminosity, mass accretion rate, and the chemical composition of both the accreting material and WD material. The properties of the outburst also depends on when, how, and if the accreted layers are mixed with the WD core and the mixing mechanism is still unknown.

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.

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.

We investigate the effects of thermonuclear reaction rate variations on 26Al production in massive stars. The dominant production sites in such events were recently investigated by using stellar model calculations: explosive neon-carbon burning, convective shell carbon burning, and convective core hydrogen burning. Post-processing nucleosynthesis calculations are performed for each of these sites by adopting temperature-density-time profiles from recent stellar evolution models. For each profile, we individually multiplied the rates of all relevant reactions by factors of 10, 2, 0.5 and 0.1, and analyzed the resulting abundance changes of 26Al. Our simulations are based on a next-generation nuclear physics library, called STARLIB, which contains a recent evaluation of Monte Carlo reaction rates. Particular attention is paid to quantifying the rate uncertainties of those reactions that most sensitively influence 26Al production. For stellar modelers our results indicate to what degree predictions of 26Al nucleosynthesis depend on currently uncertain nuclear physics input, while for nuclear experimentalists our results represent a guide for future measurements. We tabulate the results of our reaction rate sensitivity study for each of the three distinct massive star sites referred to above. It is found that several current reaction rate uncertainties influence the production of 26Al. Particularly important reactions are 26Al(n,p)26Mg, 25Mg(alpha,n)28Si, 24Mg(n,gamma)25Mg and 23Na(alpha,p)26Mg. These reactions should be prime targets for future measurements. Overall, we estimate that the nuclear physics uncertainty of the 26Al yield predicted by the massive star models explored here amounts to about a factor of 3.

The reactions that destroy 26Al in massive stars have significance in a number of astrophysical contexts. We evaluate the reaction rates of 26Al(n,p)26Mg and 26Al(n,a)23Na using cross sections obtained from the codes EMPIRE and TALYS. These have been compared to the published rates obtained from the non-smoker code and to some experimental data. We show that the results obtained from EMPIRE and TALYS are comparable to those from non-smoker. We also show how the theoretical results vary with respect to changes in the input parameters. Finally, we present recommended rates for these reactions using the available experimental data and our new theoretical results.

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.

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.

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

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.

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.

In the first lecture of this volume, we will present the basic fundamental ideas regarding nuclear processes occurring in stars. We start from stellar observations, will then elaborate on some important quantum-mechanical phenomena governing nuclear reactions, continue with how nuclear reactions proceed in a hot stellar plasma and, finally, we will provide an overview of stellar burning stages. At the end, the current knowledge regarding the origin of the elements is briefly summarized. This lecture is directed towards the student of nuclear astrophysics. Our intention is to present seemingly unrelated phenomena of nuclear physics and astrophysics in a coherent framework.

Nucleosynthesis in Type I X-ray bursts may involve up to several thousand nuclear processes. The majority of these processes have only been determined theoretically due to the lack of sufficient experimental information. Accurate reaction Q-values are essential for reliable theoretical estimates of reaction-rates. Those reactions with small Q-values (< 1 MeV) are of particular interest in these environments as they may represent waiting points for a continuous abundance flow towards heavier-mass nuclei. To explore the nature of these waiting points, we have performed a comprehensive series of post-processing calculations which examine the sensitivity of nucleosynthesis in Type I X-ray bursts to uncertainties in reaction Q-values. We discuss and list the relatively few critical masses for which measurements could better constrain the results of our studies. In particular, we stress the importance of measuring the mass of 65As to obtain an experimental Q-value for the 64Ge(p,gamma)65As reaction.

We have continued our studies of the Classical Nova outburst by evolving TNRs on 1.25Msun and 1.35Msun WDs (ONeMg composition) under conditions which produce mass ejection and a rapid increase in the emitted light, by examining the effects of changes in the nuclear reaction rates on both the observable features and the nucleosynthesis during the outburst. In order to improve our calculations over previous work, we have incorporated a modern nuclear reaction network into our hydrodynamic computer code. We find that the updates in the nuclear reaction rate libraries change the amount of ejected mass, peak luminosity, and the resulting nucleosynthesis. In addition, as a result of our improvements, we discovered that the pep reaction was not included in our previous studies of CN explosions. Although the energy production from this reaction is not important in the Sun, the densities in WD envelopes can exceed $10^4$ gm cm$^{-3}$ and the presence of this reaction increases the energy generation during the time that the p-p chain is operating. The effect of the increased energy generation is to reduce the evolution time to the peak of the TNR and, thereby, the accreted mass as compared to the evolutionary sequences done without this reaction included. As expected from our previous work, the reduction in accreted mass has important consequences on the characteristics of the resulting TNR and is discussed in this paper.

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.

The 16O(p,gamma)17F reaction rate is revisited with special emphasis on the stellar temperature range of T=60-100 MK important for hot bottom burning in asymptotic giant branch (AGB) stars. We evaluate existing cross section data that were obtained since 1958 and, if appropriate, correct published data for systematic errors that were not noticed previously, including the effects of coincidence summing and updated effective stopping powers. The data are interpreted by using two different models of nuclear reactions, that is, a potential model and R-matrix theory. A new astrophysical S-factor and recommended thermonuclear reaction rates are presented. As a result of our work, the 16O(p,gamma)17F reaction has now the most precisely known rate involving any target nucleus in the mass A >= 12 range, with reaction rate errors of about 7% over the entire temperature region of astrophysical interest (T=0.01-2.5 GK). The impact of the present improved reaction rate with its significantly reduced uncertainties on the hot bottom burning in AGB stars is discussed. In contrast to earlier results we find now that there is not clear evidence to date for any stellar grain origin from massive AGB stars.

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.

Stardust is a class of presolar grains each of which presents an ideally uncontaminated stellar sample. Mainstream silicon carbide (SiC) stardust formed in the extended envelopes of carbon-rich asymptotic giant branch (AGB) stars and incorporated the radioactive nucleus 26Al as a trace element. The aim of this paper is to analyse in detail the effect of nuclear uncertainties, in particular the large uncertainties of up to four orders of magnitude related to the 26Al_g+(p,gamma)27Si reaction rate, on the production of 26Al in AGB stars and compare model predictions to data obtained from laboratory analysis of SiC stardust grains. Stellar uncertainties are also briefly discussed. We use a detailed nucleosynthesis postprocessing code to calculate the 26Al/27Al ratios at the surface of AGB stars of different masses (M = 1.75, 3, and 5 M_sun) and metallicities (Z = 0.02, 0.012, and 0.008). For the lower limit and recommended value of the 26Al_g(p,gamma)27Si reaction rate, the predicted 26Al/27Al ratios replicate the upper values of the range of the 26Al/27Al ratios measured in SiC grains. For the upper limit of the 26Al_g(p,gamma)27Si reaction rate, instead, the predicted 26Al/27Al ratios are approximately 100 times lower and lie below the range observed in SiC grains. When considering models of different masses and metallicities, the spread of more than an order of magnitude in the 26Al/27Al ratios measured in stellar SiC grains is not reproduced. We propose two scenarios to explain the spread of the 26Al/27Al ratios observed in mainstream SiC, depending on the choice of the 26Al_g+p reaction rate. One involves different times of stardust formation, the other involves extra-mixing processes. Stronger conclusions will be possible after more information is available from future nuclear experiments on the 26Al_g+p reaction.

Classical novae participate in the cycle of Galactic chemical evolution in which grains and metal enriched gas in their ejecta, supplementing those of supernovae, AGB stars, and Wolf-Rayet stars, are a source of heavy elements for the ISM. Once in the diffuse gas, this material is mixed with the existing gases and then incorporated into young stars and planetary systems during star formation. Infrared observations have confirmed the presence of carbon, SiC, hydrocarbons, and oxygen-rich silicate grains in nova ejecta, suggesting that some fraction of the pre-solar grains identified in meteoritic material come from novae. The mean mass returned by a nova outburst to the ISM probably exceeds ~2 x 10^{-4} Solar Masses. Using the observed nova rate of 35 per year in our Galaxy, it follows that novae introduce more than ~7 x 10^{-3} Solar Masses per year of processed matter into the ISM. Novae are expected to be the major source of 15N and 17O in the Galaxy and to contribute to the abundances of other isotopes in this atomic mass range. Here, we report on how changes in the nuclear reaction rates affect the properties of the outburst and alter the predictions of the contributions of novae to Galactic chemical evolution. We also discuss the necessity of including the pep reaction in studies of thermonuclear runaways in material accreted onto white dwarfs.

Hydrogen-burning is the root cause of the star-to-star abundance variations of light nuclei in Galactic globular clusters (GC). In the present work we constrain the physical conditions that gave rise to the observed abundance patterns of Li, C, N, O, Na, Mg, Al, as well as Mg isotopes in the typical case of NGC6752. We perform nucleosynthesis calculations at constant temperature, adopting realistic initial abundances for the proto-cluster gas. We use a detailed nuclear reaction network and state-of-the-art nuclear reaction rates. Although simplistic, our analysis provides original results and new constraints on the self-enrichment scenario for GCs. Our parametrized calculations allow us to determine a narrow range of temperature where the observed extreme abundances of all light elements and isotopes in NGC6752 are nicely reproduced simultaneously. This agreement is obtained after mixing of the H-processed material with 30 % of unprocessed gas. The observed C-N, O-Na, Mg-Al, Li-Na and F-Na anticorrelations, as well as the behaviour of the Mg isotopes can be recovered by assuming mixing with even larger dilution factors. Li production by the stars that build up the other abundance anomalies is not mandatory in the case of NGC 6752. Observations of O, Na, Mg and Al constrain the temperature range for H-burning; such temperatures are encountered in the two main candidate ``polluters'' proposed for GCs, namely massive AGBs and the most massive main-sequence stars. (Abridged).

We test the effect of proton-capture reaction rate uncertainties on the abundances of the Ne, Na, Mg and Al isotopes processed by the NeNa and MgAl chains during hot bottom burning (HBB) in asymptotic giant branch (AGB) stars of intermediate mass between 4 and 6 solar masses and metallicities between Z=0.0001 and 0.02. We provide uncertainty ranges for the AGB stellar yields, for inclusion in galactic chemical evolution models, and indicate which reaction rates are most important and should be better determined. We use a fast synthetic algorithm based on detailed AGB models. We run a large number of stellar models, varying one reaction per time for a very fine grid of values, as well as all reactions simultaneously. We show that there are uncertainties in the yields of all the Ne, Na, Mg and Al isotopes due to uncertain proton-capture reaction rates. The most uncertain yields are those of 26Al and 23Na (variations of two orders of magnitude), 24Mg and 27Al (variations of more than one order of magnitude), 20Ne and 22Ne (variations between factors 2 and 7). In order to obtain more reliable Ne, Na, Mg and Al yields from IM-AGB stars the rates that require more accurate determination are: 22Ne(p,g)23Na, 23Na(p,g)24Mg, 25Mg(p,g)26Al, 26Mg(p,g)27Al and 26Al(p,g)27Si. Detailed galactic chemical evolution models should be constructed to address the impact of our uncertainty ranges on the observational constraints related to HBB nucleosynthesis, such as globular cluster chemical anomalies.

Microscopic presolar grains extracted from primitive meteorites have extremely anomalous isotopic compositions revealing the stellar origin of these grains. The composition of presolar spinel grain OC2 is different from that of all other presolar spinel grains. Large excesses of the heavy Mg isotopes are present and thus an origin from an intermediate-mass (IM) asymptotic giant branch (AGB) star was previously proposed for this grain. We discuss the isotopic compositions of presolar spinel grain OC2 and compare them to theoretical predictions. We show that the isotopic composition of O, Mg and Al in OC2 could be the signature of an AGB star of IM and metallicity close to solar experiencing hot bottom burning, or of an AGB star of low mass (LM) and low metallicity suffering very efficient cool bottom processing. Large measurement uncertainty in the Fe isotopic composition prevents us from discriminating which model better represents the parent star of OC2. However, the Cr isotopic composition of the grain favors an origin in an IM-AGB star of metallicity close to solar. Our IM-AGB models produce a self-consistent solution to match the composition of OC2 within the uncertainties related to reaction rates. Within this solution we predict that the 16O(p,g)17F and the 17O(p,a)14N reaction rates should be close to their lower and upper limits, respectively. By finding more grains like OC2 and by precisely measuring their Fe and Cr isotopic compositions, it may be possible in the future to derive constraints on massive AGB models from the study of presolar grains.

We study the effect of uncertainties in the proton-capture reaction rates of the NeNa and MgAl chains on nucleosynthesis due to the operation of hot bottom burning (HBB) in intermediate-mass asymptotic giant branch (AGB) stars. HBB nucleosynthesis is associated with the production of sodium, radioactive Al26 and the heavy magnesium isotopes, and it is possibly responsible for the O, Na, Mg and Al abundance anomalies observed in globular cluster stars. We model HBB with an analytic code based on full stellar evolution models so we can quickly cover a large parameter space. The reaction rates are varied first individually, then all together. This creates a knock-on effect, where an increase of one reaction rate affects production of an isotope further down the reaction chain. We find the yields of Ne22, Na23 and Al26 to be the most susceptible to current nuclear reaction rate uncertainties.

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.