The production of about half of the heavy elements found in nature is assigned to a specific astrophysical nucleosynthesis process: the rapid neutron capture process (r-process). Although this idea has been postulated more than six decades ago, the full understanding faces two types of uncertainties/open questions: (a) The nucleosynthesis path in the nuclear chart runs close to the neutron-drip line, where presently only limited experimental information is available, and one has to rely strongly on theoretical predictions for nuclear properties. (b) While for many years the occurrence of the r-process has been associated with supernovae, more recent studies have cast substantial doubts on this environment. Alternative scenarios include the mergers of neutron stars, neutron-star black hole mergers, but possibly also rare classes of supernovae as well as hypernovae/collapsars with polar jet ejecta and also accretion disk outflows related to the collapse of fast rotating massive stars with high magnetic fields. Stellar r-process abundance observations, have provided insights into, and constraints on the frequency of and conditions in the responsible stellar production sites. One of them, neutron star mergers, was just identified and related to the Gravitational Wave event GW170817. High resolution observations, increasingly more precise due to improved experimental atomic data, have been particularly important in defining the heavy element abundance patterns of the old halo stars, and thus determining the extent, and nature, of the earliest nucleosynthesis in our Galaxy. Combining new results and important breakthroughs in the related nuclear, atomic and astronomical fields of science, this review attempts to provide an answer to the question "How Were the Elements from Iron to Uranium Made?" (Abridged)

We present first results and plans for future neutron activation measurements at the Goethe University Frankfurt. The measurements were performed at the Van-de-Graaff accelerator employing the 7Li(p,n) reaction.

Neutron capture cross sections are one of the most important nuclear inputs to models of stellar nucleosynthesis of the elements heavier than iron. The activation technique and the time-of-flight method are mostly used to determine the required data experimentally. Recent developments of experimental techniques allow for new experiments on radioactive isotopes. Monte-Carlo based analysis methods give new insights into the systematic uncertainties of previous measurements. We present an overview over the state-of-the-art experimental techniques, a detailed new evaluation of the $^{197}$Au(n,$\gamma$) cross section in the keV-regime and the corresponding re-evaluation of 63 more isotopes, which have been measured in the past relative to the gold cross section.

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.

X-ray observations of transiently accreting neutron stars during quiescence provide information about the structure of neutron star crusts and the properties of dense matter. Interpretation of the observational data requires an understanding of the nuclear reactions that heat and cool the crust during accretion, and define its nonequilibrium composition. We identify here in detail the typical nuclear reaction sequences down to a depth in the inner crust where the mass density is 2E12 g/cm^3 using a full nuclear reaction network for a range of initial compositions. The reaction sequences differ substantially from previous work. We find a robust reduction of crust impurity at the transition to the inner crust regardless of initial composition, though shell effects can delay the formation of a pure crust somewhat to densities beyond 2E12 g/cm^3. This naturally explains the small inner crust impurity inferred from observations of a broad range of systems. The exception are initial compositions with A >= 102 nuclei, where the inner crust remains impure with an impurity parameter of Qimp~20 due to the N = 82 shell closure. In agreement with previous work we find that nuclear heating is relatively robust and independent of initial composition, while cooling via nuclear Urca cycles in the outer crust depends strongly on initial composition. This work forms a basis for future studies of the sensitivity of crust models to nuclear physics and provides profiles of composition for realistic crust models.

The solar s-process abundances have been analyzed in the framework of a Galactic Chemical Evolution (GCE) model. The aim of this work is to implement the study by Bisterzo et al. (2014), who investigated the effect of one of the major uncertainties of asymptotic giant branch (AGB) yields, the internal structure of the 13C pocket. We present GCE predictions of s-process elements computed with additional tests in the light of the suggestions provided in recent publications. The analysis is extended to different metallicities, by comparing GCE results and updated spectroscopic observations of unevolved field stars. We verify that the GCE predictions obtained with different tests may represent, on average, the evolution of selected neutron-capture elements in the Galaxy. The impact of an additional weak s-process contribution from fast-rotating massive stars is also explored.

The synthesis of heavy, proton rich isotopes in the astrophysical gamma-process proceeds through photodisintegration reactions. For the improved understanding of the process, the rates of the involved nuclear reactions must be known. The reaction 128Ba(g,a)124Xe was found to affect the abundance of the p nucleus 124Xe. Since the stellar rate for this reaction cannot be determined by a measurement directly, the aim of the present work was to measure the cross section of the inverse 124Xe(a,g)128Ba reaction and to compare the results with statistical model predictions. Of great importance is the fact that data below the (a,n) threshold was obtained. Studying simultaneously the 124Xe(a,n)127Ba reaction channel at higher energy allowed to further identify the source of a discrepancy between data and prediction. The 124Xe + alpha cross sections were measured with the activation method using a thin window 124Xe gas cell. The studied energy range was between E = 11 and 15 MeV close above the astrophysically relevant energy range. The obtained cross sections are compared with statistical model calculations. The experimental cross sections are smaller than standard predictions previously used in astrophysical calculations. As dominating source of the difference, the theoretical alpha width was identified. The experimental data suggest an alpha width lower by at least a factor of 0.125 in the astrophysical energy range. An upper limit for the 128Ba(g,a)124Xe stellar rate was inferred from our measurement. The impact of this rate was studied in two different models for core-collapse supernova explosions of 25 solar mass stars. A significant contribution to the 124Xe abundance via this reaction path would only be possible when the rate was increased above the previous standard value. Since the experimental data rule this out, they also demonstrate the closure of this production path.

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.

This paper provides a detailed analysis of the main component of the slow neutron capture process (the s-process), which accounts for the solar abundances of half of the nuclei with 90 <~ A <~ 208. We examine the impact of the uncertainties of the two neutron sources operating in low-mass asymptotic giant branch (AGB) stars: the 13C(alpha, n)16O reaction, which releases neutrons radiatively during interpulse periods (kT ~ 8 keV), and the 22Ne(alpha, n)25Mg reaction, partially activated during the convective thermal pulses (TPs). We focus our attention on the branching points that mainly influence the abundance of s-only isotopes. In our AGB models, the 13C is fully consumed radiatively during interpulse. In this case, we find that the present uncertainty associated to the 13C(alpha, n)16O reaction has marginal effects on s-only nuclei. On the other hand, a reduction of this rate may increase the amount of residual (or unburned) 13C at the end of the interpulse: in this condition, the residual 13C is burned at higher temperature in the convective zone powered by the following TP. The neutron burst produced by the 22Ne(alpha, n)25Mg reaction has major effects on the branches along the s path. The contributions of s-only isotopes with 90 <~ A <= 204 are reproduced within solar and nuclear uncertainties, even if the 22Ne(alpha, n)25Mg rate is varied by a factor of two. Improved beta-decay and neutron capture rates of a few key radioactive nuclides would help to attain a comprehensive understanding of the solar main component.

Neutrons produced by the carbon fusion reaction 12C(12C,n)23Mg play an important role in stellar nucleosynthesis. However, past studies have shown large discrepancies between experimental data and theory, leading to an uncertain cross section extrapolation at astrophysical energies. We present the first direct measurement that extends deep into the astrophysical energy range along with a new and improved extrapolation technique based on experimental data from the mirror reaction 12C(12C,p)23Na. The new reaction rate has been determined with a well-defined uncertainty that exceeds the precision required by astrophysics models. Using our constrained rate, we find that 12C(12C,n)23Mg is crucial to the production of Na and Al in Pop-III Pair Instability Supernovae. It also plays a non-negligible role in the production of weak s-process elements as well as in the production of the important galactic gamma-ray emitter 60Fe.

Baryon density inhomogeneities in the early universe can give rise to a floor of heavy elements (up to $A\approx 270$) produced in a primordial r-process with fission cycling. A parameter study with variation of the global baryon to photon ratio $\eta$, and under inclusion of neutron diffusion effects was performed. New results concerning the dependence of the results on nuclear physics parameters are presented.

We study the s-process abundances (A > 90) at the epoch of the solar-system formation. AGB yields are computed with an updated neutron capture network and updated initial solar abundances. We confirm our previous results obtained with a Galactic Chemical Evolution (GCE) model: (i) as suggested by the s-process spread observed in disk stars and in presolar meteoritic SiC grains, a weighted average of s-process strengths is needed to reproduce the solar s-distribution of isotopes with A > 130; (ii) an additional contribution (of about 25%) is required in order to represent the solar s-process abundances of isotopes from A = 90 to 130. Furthermore, we investigate the effect of different internal structures of the 13C-pocket, which may affect the efficiency of the 13C(a, n)16O reaction, the major neutron source of the s-process. First, keeping the same 13C profile adopted so far, we modify by a factor of two the mass involved in the pocket; second, we assume a flat 13C profile in the pocket, and we test again the effects of the variation of the mass of the pocket. We find that GCE s-predictions at the epoch of the solar-system formation marginally depend on the size and shape of the 13C-pocket once a different weighted range of 13C-pocket strengths is assumed. We ascertain that, independently of the internal structure of the 13C-pocket, the missing solar-system s-process contribution in the range from A = 90 to 130 remains essentially the same.

The temperature in the crust of an accreting neutron star, which comprises its outermost kilometer, is set by heating from nuclear reactions at large densities, neutrino cooling, and heat transport from the interior. The heated crust has been thought to affect observable phenomena at shallower depths, such as thermonuclear bursts in the accreted envelope. Here we report that cycles of electron capture and its inverse, $\beta^-$ decay, involving neutron-rich nuclei at a typical depth of about 150 m, cool the outer neutron star crust by emitting neutrinos while also thermally decoupling the surface layers from the deeper crust. This Urca mechanism has been studied in the context of white dwarfs and Type Ia supernovae, but hitherto was not considered in neutron stars, because previous models computed the crust reactions using a zero-temperature approximation and assumed that only a single nuclear species was present at any given depth. This thermal decoupling means that X-ray bursts and other surface phenomena are largely independent of the strength of deep crustal heating. The unexpectedly short recurrence times, of the order of years, observed for very energetic thermonuclear superbursts are therefore not an indicator of a hot crust, but may point instead to an unknown local heating mechanism near the neutron star surface.

We provide an updated discussion of the sample of CEMP-s and CEMP-s/r stars collected from the literature. Observations are compared with the theoretical nucleosynthesis models of asymptotic giant branch (AGB) stars presented by Bisterzo et al. (2010, 2011, 2012), in the light of the most recent spectroscopic results.

We study the s-process abundances at the epoch of the Solar-system formation as the outcome of nucleosynthesis occurring in AGB stars of various masses and metallicities. The calculations have been performed with the Galactic chemical evolution (GCE) model presented by Travaglio et al. (1999, 2004). With respect to previous works, we used updated solar meteoritic abundances, a neutron capture cross section network that includes the most recent measurements, and we implemented the $s$-process yields with an extended range of AGB initial masses. The new set of AGB yields includes a new evaluation of the 22Ne(alpha, n)25Mg rate, which takes into account the most recent experimental information.

The nucleosynthesis of elements beyond iron is dominated by neutron captures in the s and r processes. However, 32 stable, proton-rich isotopes cannot be formed during those processes, because they are shielded from the s-process flow and r-process beta-decay chains. These nuclei are attributed to the p and rp process. For all those processes, current research in nuclear astrophysics addresses the need for more precise reaction data involving radioactive isotopes. Depending on the particular reaction, direct or inverse kinematics, forward or time-reversed direction are investigated to determine or at least to constrain the desired reaction cross sections. The Facility for Antiproton and Ion Research (FAIR) will offer unique, unprecedented opportunities to investigate many of the important reactions. The high yield of radioactive isotopes, even far away from the valley of stability, allows the investigation of isotopes involved in processes as exotic as the r or rp processes.

About a year after core collapse supernova, dust starts to condense in the ejecta. In meteorites, a fraction of C-rich presolar grains (e.g., silicon carbide (SiC) grains of Type-X and low density graphites) are identified as relics of these events, according to the anomalous isotopic abundances. Several features of these abundances remain unexplained and challenge the understanding of core-collapse supernovae explosions and nucleosynthesis. We show, for the first time, that most of the measured C-rich grain abundances can be accounted for in the C-rich material from explosive He burning in core-collapse supernovae with high shock velocities and consequent high temperatures. The inefficiency of the $^{12}$C($\alpha$,$\gamma$)$^{16}$O reaction relative to the rest of the $\alpha$-capture chain at $T > 3.5\times10^8 \mathrm{K}$ causes the deepest He-shell material to be carbon rich and silicon rich, and depleted in oxygen. The isotopic ratio predictions in part of this material, defined here as the C/Si zone, are in agreement with the grain data. The high-temperature explosive conditions that our models reach at the bottom of the He shell, can also be representative of the nucleosynthesis in hypernovae or in the high-temperature tail of a distribution of conditions in asymmetric supernovae. Finally, our predictions are consistent with the observation of large $^{44}$Ca/$^{40}$Ca observed in the grains. This is due to the production of $^{44}$Ti together with $^{40}$Ca in the C/Si zone, and/or to the strong depletion of $^{40}$Ca by neutron captures.

Despite much effort in the past decades, the C-burning reaction rate is uncertain by several orders of magnitude, and the relative strength between the different channels 12C(12C,alpha)20Ne, 12C(12C,p)23Na and 12C(12C,n)23Mg is poorly determined. Additionally, in C-burning conditions a high 12C+12C rate may lead to lower central C-burning temperatures and to 13C(alpha,n)16O emerging as a more dominant neutron source than 22Ne(alpha,n)25Mg, increasing significantly the s-process production. This is due to the rapid decrease of the 13N(gamma,p)12C with decreasing temperature, causing the 13C production via 13N(beta+)13C. Presented here is the impact of the 12C+12C reaction uncertainties on the s-process and on explosive p-process nucleosynthesis in massive stars, including also fast rotating massive stars at low metallicity. Using various 12C+12C rates, in particular an upper and lower rate limit of ~ 50000 higher and ~ 20 lower than the standard rate at 5*10^8 K, five 25 Msun stellar models are calculated. The enhanced s-process signature due to 13C(alpha,n)16O activation is considered, taking into account the impact of the uncertainty of all three C-burning reaction branches. Consequently, we show that the p-process abundances have an average production factor increased up to about a factor of 8 compared to the standard case, efficiently producing the elusive Mo and Ru proton-rich isotopes. We also show that an s-process being driven by 13C(alpha,n)16O is a secondary process, even though the abundance of 13C does not depend on the initial metal content. Finally, implications for the Sr-peak elements inventory in the Solar System and at low metallicity are discussed.

Asymptotic giant branch (AGB) stars with low initial mass (1 - 3 Msun) are responsible for the production of neutron-capture elements through the main s-process (main slow neutron capture process). The major neutron source is 13C(alpha, n)16O, which burns radiatively during the interpulse periods at about 8 keV and produces a rather low neutron density (10^7 n/cm^3). The second neutron source 22Ne(alpha, n)25Mg, partially activated during the convective thermal pulses when the energy reaches about 23 keV, gives rise to a small neutron exposure but a peaked neutron density (Nn(peak) > 10^11 n/cm^3). At metallicities close to solar, it does not substantially change the final s-process abundances, but mainly affects the isotopic ratios near s-path branchings sensitive to the neutron density. We examine the effect of the present uncertainties of the two neutron sources operating in AGB stars, as well as the competition with the 22Ne(alpha, gamma)26Mg reaction. The analysis is carried out on AGB the main-s process component (reproduced by an average between M(AGB; ini) = 1.5 and 3 Msun at half solar metallicity, see Arlandini et al. 1999), using a set of updated nucleosynthesis models. Major effects are seen close to the branching points. In particular, 13C(alpha, n)16O mainly affects 86Kr and 87Rb owing to the branching at 85Kr, while small variations are shown for heavy isotopes by decreasing or increasing our adopted rate by a factor of 2 - 3. By changing our 22Ne(alpha, n)25Mg rate within a factor of 2, a plausible reproduction of solar s-only isotopes is still obtained. We provide a general overview of the major consequences of these variations on the s-path. A complete description of each branching will be presented in Bisterzo et al., in preparation.

We consider the astrophysical reaction rates for radiative neutron capture reactions ($n,\gamma$) in the crust of a neutron star. The presence of degenerate neutrons at high densities (mainly in the inner crust) can drastically affect the reaction rates. Standard rates assuming a Maxwell-Boltzmann distribution for neutrons can underestimate the rates by several orders of magnitude. We derive simple analytical expressions for reaction rates at a variety of conditions with account for neutron degeneracy. We also discuss the plasma effects on the outgoing radiative transition channel in neutron radiative capture reactions and show that these effects can also increase the reaction rates by a few orders of magnitude. In addition, using detailed balance, we analyze the effects of neutron degeneracy and plasma physics on reverse ($\gamma,n$) photodisintegration. We discuss the dependence of the reaction rates on temperature and neutron chemical potential and outline the efficiency of these reactions in the neutron star crust.

Numerous nuclear reactions in the crust of accreting neutron stars are strongly affected by dense plasma environment. Simulations of superbursts, deep crustal heating and other nuclear burning phenomena in neutron stars require astrophysical S-factors for these reactions (as a function of center-of-mass energy E of colliding nuclei). A large database of S-factors is created for about 5000 non-resonant fusion reactions involving stable and unstable isotopes of Be, B, C, N, O, F, Ne, Na, Mg, and Si. It extends the previous database of about 1000 reactions involving isotopes of C, O, Ne, and Mg. The calculations are performed using the Sao Paulo potential and the barrier penetration formalism. All calculated S-data are parameterized by an analytic model for S(E) proposed before [Phys. Rev. C 82, 044609 (2010)] and further elaborated here. For a given reaction, the present S(E)-model contains three parameters. These parameters are easily interpolated along reactions involving isotopes of the same elements with only seven input parameters, giving an ultracompact, accurate, simple, and uniform database. The S(E) approximation can also be used to estimate theoretical uncertainties of S(E) and nuclear reaction rates in dense matter, as illustrated for the case of the 34Ne+34Ne reaction in the inner crust of an accreting neutron star.

[Shortened] The 12C + 12C fusion reaction has been the subject of considerable experimental efforts to constrain uncertainties at temperatures relevant for stellar nucleosynthesis. In order to investigate the effect of an enhanced carbon burning rate on massive star structure and nucleosynthesis, new stellar evolution models and their yields are presented exploring the impact of three different 12C + 12C reaction rates. Non-rotating stellar models were generated using the Geneva Stellar Evolution Code and were later post-processed with the NuGrid Multi-zone Post-Processing Network tool. The enhanced rate causes core carbon burning to be ignited more promptly and at lower temperature. This reduces the neutrino losses, which increases the core carbon burning lifetime. An increased carbon burning rate also increases the upper initial mass limit for which a star exhibits a convective carbon core. Carbon shell burning is also affected, with fewer convective-shell episodes and convection zones that tend to be larger in mass. Consequently, the chance of an overlap between the ashes of carbon core burning and the following carbon shell convection zones is increased, which can cause a portion of the ashes of carbon core burning to be included in the carbon shell. Therefore, during the supernova explosion, the ejecta will be enriched by s-process nuclides synthesized from the carbon core s process. The yields were used to estimate the weak s-process component in order to compare with the solar system abundance distribution. The enhanced rate models were found to produce a significant proportion of Kr, Sr, Y, Zr, Mo, Ru, Pd and Cd in the weak component, which is primarily the signature of the carbon-core s process. Consequently, it is shown that the production of isotopes in the Kr-Sr region can be used to constrain the 12C + 12C rate using the current branching ratio for a- and p-exit channels.

We explore the nucleosynthesis trends from two mechanisms during freeze-out expansions in core-collapse supernovae. The first mechanism is related to the convection and instabilities within homogeneous stellar progenitor matter that is accreted through the supernova shock. The second mechanism is related to the impact of the supersonic wind termination shock (reverse shock) within the tumultuous inner regions of the ejecta above the proto-neutron star. Our results suggest that isotopes in the mass range 12<=A<=122 that are produced during the freeze-out expansions may be classified in two families. The isotopes of the first family manifest a common mass fraction evolutionary profile, whose specific shape per isotope depends on the characteristic transition between two equilibrium states (equilibrium state transition) during each type of freeze-out expansion. The first family includes the majority of isotopes in this mass range. The second family is limited to magic nuclei and isotopes in their locality, which do not sustain any transition, become nuclear flow hubs, and dominate the final composition. We use exponential and power-law adiabatic profiles to identify dynamic large-scale and small-scale equilibrium patterns among nuclear reactions. (truncated abstract) In addition, we introduce non-monotonic parameterized profiles to probe the impact of the reverse shock and multi-dimensional explosion asymmetries on nucleosynthesis. Cases are shown in which the non-monotonic profiles favor the production of radioactivities. Non-monotonic freeze-out profiles involve longer non-equilibrium nucleosynthesis intervals compared with the exponential and power-law profiles, resulting in mass fraction trends and yield distributions that may not be achieved by the monotonic freeze-out profiles.

The location of electron capture heat sources in the crust of accreting neutron stars depends on the masses of extremely neutron-rich nuclei. We present first results from a new implementation of the time-of-flight technique to measure nuclear masses of rare isotopes at the National Superconducting Cyclotron Laboratory. The masses of 16 neutron-rich nuclei in the scandium -- nickel range were determined simultaneously, improving the accuracy compared to previous data in 12 cases. The masses of $^{61}${V}, $^{63}${Cr}, $^{66}${Mn}, and $^{74}${Ni} were measured for the first time with mass excesses of $-30.510(890)$ MeV, $-35.280(650)$ MeV, $-36.900(790)$ MeV, and $-49.210(990)$ MeV, respectively. With the measurement of the $^{66}$Mn mass, the locations of the two dominant electron capture heat sources in the outer crust of accreting neutron stars that exhibit superbursts are now experimentally constrained. We find that the location of the $^{66}$Fe$\rightarrow^{66}$Mn electron capture transition occurs significantly closer to the surface than previously assumed because our new experimental Q-value is 2.1 MeV (2.6$\sigma$) smaller than predicted by the FRDM mass model.

The contribution by massive stars (M > 9 solar masses) to the weak s-process component of the solar system abundances is primarily due to the 22Ne neutron source, which is activated near the end of helium-core burning. The residual 22Ne left over from helium-core burning is then reignited during carbon burning, initiating further s-processing that modifies the isotopic distribution. This modification is sensitive to the stellar structure and the carbon burning reaction rate. Recent work on the 12C + 12C reaction suggests that resonances located within the Gamow peak may exist, causing a strong increase in the astrophysical S-factor and consequently the reaction rate. To investigate the effect of an increased rate, 25 solar mass stellar models with three different carbon burning rates, at solar metallicity, were generated using the Geneva Stellar Evolution Code (GENEC) with nucleosynthesis post-processing calculated using the NuGrid Multi-zone Post-Processing Network code (MPPNP). The strongest rate caused carbon burning to occur in a large convective core rather than a radiative one. The presence of this large convective core leads to an overlap with the subsequent convective carbon-shell, significantly altering the initial composition of the carbon-shell. In addition, an enhanced rate causes carbon-shell burning episodes to ignite earlier in the evolution of the star, igniting the 22Ne source at lower temperatures and reducing the neutron density.

Nuclear reaction rates determine the abundances of isotopes in stellar burning processes. A multitude of reactions determine the reaction flow pattern which is described in terms of reaction network simulations. The reaction rates are determined by laboratory experiments supplemented by nuclear reaction and structure theory. We will discuss the experimental approach as well as the theoretical tools for obtaining the stellar reaction rates. A detailed analysis of a reaction is only possible for a few selected cases which will be highlighted in this section. The bulk of nuclear reaction processes is however described in terms of a statistical model approach, which relies on global nuclear structure and reaction parameters such as level density and mass and barrier penetration, respectively. We will discuss a variety of experimental facilities and techniques used in the field, this includes low energy stable beam experiments, measurements at radioactive beam accelerators, and neutron beam facilities.

We propose a physically transparent analytic model of astrophysical S-factors as a function of a center-of-mass energy E of colliding nuclei (below and above the Coulomb barrier) for non-resonant fusion reactions. For any given reaction, the S(E)-model contains four parameters [two of which approximate the barrier potential, U(r)]. They are easily interpolated along many reactions involving isotopes of the same elements; they give accurate practical expressions for S(E) with only several input parameters for many reactions. The model reproduces the suppression of S(E) at low energies (of astrophysical importance) due to the shape of the low-r wing of U(r). The model can be used to reconstruct U(r) from computed or measured S(E). For illustration, we parameterize our recent calculations of S(E) (using the Sao Paulo potential and the barrier penetration formalism) for 946 reactions involving stable and unstable isotopes of C, O, Ne, and Mg (with 9 parameters for all reactions involving many isotopes of the same elements, e.g., C+O). In addition, we analyze astrophysically important 12C+12C reaction, compare theoretical models with experimental data, and discuss the problem of interpolating reliably known S(E) values to low energies (E <= 2-3 MeV).

We summarize and critically evaluate the available data on nuclear fusion cross sections important to energy generation in the Sun and other hydrogen-burning stars and to solar neutrino production. Recommended values and uncertainties are provided for key cross sections, and a recommended spectrum is given for 8B solar neutrinos. We also discuss opportunities for further increasing the precision of key rates, including new facilities, new experimental techniques, and improvements in theory. This review, which summarizes the conclusions of a workshop held at the Institute for Nuclear Theory, Seattle, in January 2009, is intended as a 10-year update and supplement to Reviews of Modern Physics 70 (1998) 1265.

We compare the yields of Ti44 and Ni56 produced from post-processing the thermodynamic trajectories from three different core-collapse models -- a Cassiopeia A progenitor, a double shock hypernova progenitor, and a rotating 2D explosion -- with the yields from exponential and power-law trajectories. The peak temperatures and densities achieved in these core-collapse models span several of the distinct nucleosynthesis regions we identify, resulting in different trends in the Ti44 and Ni56 yields for different mass elements. The Ti44 and Ni56 mass fraction profiles from the exponential and power-law profiles generally explain the tendencies of the post-processed yields, depending on which regions are traversed by the model. We find integrated yields of Ti44 and Ni56 from the exponential and power-law trajectories are generally within a factor 2 or less of the post-process yields. We also analyze the influence of specific nuclear reactions on the Ti44 and Ni56 abundance evolution. Reactions that affect all yields globally are the 3a, p(e-,nu)n and n(e+,nubar)p. The rest of the reactions are ranked according to their degree of impact on the synthesis of Ti44. The primary ones include Ti44(a,p)V47, Ca40(a,g)Ti44, V45(p,g)Cr46, Ca40(a,p)Sc43, F17(a,p)Ne20, Na21(a,p)Mg24, Sc41(p,g)Ti42, Sc43(p,g)Ti44, Ti44(p,g)V45, and Ni57(p,g)Cu58, along with numerous weak reactions. Our analysis suggests that not all Ti44 need be produced in an a-rich freeze-out in core-collapse events, and that reaction rate equilibria in combination with timescale effects for the expansion profile may account for the paucity of Ti44 observed in supernovae remnants.

The slow neutron capture process in massive stars (the weak s-process) produces most of the s-only isotopes in the mass region 60 < A < 90. The nuclear reaction rates used in simulations of this process have a profound effect on the final s-process yields. We generated 1D stellar models of a 25 solar mass star varying the 12C + 12C rate by a factor of 10 and calculated full nucleosynthesis using the post-processing code PPN. Increasing or decreasing the rate by a factor of 10 affects the convective history and nucleosynthesis, and consequently the final yields.

Using the Sao Paulo potential and the barrier penetration formalism we have calculated the astrophysical factor S(E) for 946 fusion reactions involving stable and neutron-rich isotopes of C, O, Ne, and Mg for center-of-mass energies E varying from 2 MeV to 18-30 MeV (covering the range below and above the Coulomb barrier). We have parameterized the energy dependence S(E) by an accurate universal 9-parameter analytic expression and present tables of fit parameters for all the reactions. We also discuss the reduced 3-parameter version of our fit which is highly accurate at energies below the Coulomb barrier, and outline the procedure for calculating the reaction rates. The results can be easily converted to thermonuclear or pycnonuclear reaction rates to simulate various nuclear burning phenomena, in particular, stellar burning at high temperatures and nucleosynthesis in high density environments.

The neutron capture cross section of 14C is of relevance for several nucleosynthesis scenarios such as inhomogeneous Big Bang models, neutron induced CNO cycles, and neutrino driven wind models for the r process. The 14C(n,g) reaction is also important for the validation of the Coulomb dissociation method, where the (n,g) cross section can be indirectly obtained via the time-reversed process. So far, the example of 14C is the only case with neutrons where both, direct measurement and indirect Coulomb dissociation, have been applied. Unfortunately, the interpretation is obscured by discrepancies between several experiments and theory. Therefore, we report on new direct measurements of the 14C(n,g) reaction with neutron energies ranging from 20 to 800 keV.

This white paper, directed to the Stars and Stellar Evolution panel, has three objectives: 1) to provide the Astro2010 Decadal Survey with a vista into the goals of the nuclear physics and nuclear astrophysics community; 2) to alert the astronomical community of joint opportunities for discoveries at the interface between nuclear physics and astronomy; and 3) to delineate efforts in nuclear physics and describe the observational and theoretical advances in astrophysics necessary to make progress towards answering the following questions in the Nuclear Science 2007 Long Range Plan: 1) What is the origin and distribution of the elements? 2) What are the nuclear reactions that power stars and stellar explosions? 3) What is the nature of dense matter? The scope of this white paper concerns the specific area of "low energy" nuclear astrophysics. We define this as the area of overlap between astrophysics and the study of nuclear structure and reactions. Of the questions listed above, two -- What is the origin of the elements? and What is the nature of dense matter? -- were specifically listed in the National Academies Study, "Connecting Quarks with the Cosmos".

The recent experimental evaluation of the 18F(a,p)21Ne reaction rate, when considering its associated uncertainties, presented significant differences compared to the theoretical Hauser-Feshbach rate. This was most apparent at the low temperatures relevant for He-shell burning in asymptotic giant branch (AGB) stars. Investigations into the effect on AGB nucleosynthesis revealed that the upper limit resulted in an enhanced production of 19F and 21Ne in carbon-rich AGB models, but the recommended and lower limits presented no differences from using the theoretical rate. This was the case for models spanning a range in metallicity from solar to [Fe/H] ~ -2.3. The results of this study are relevant for observations of F and C-enriched AGB stars in the Galaxy, and to the Ne composition of mainstream silicon carbide grains, that supposedly formed in the outflows of cool, carbon-rich giant stars. We discuss the mechanism that produces the extra F and summarize our main findings.

We investigate the physical conditions where 44Ti and 56Ni are created in core-collapse supernovae. In this preliminary work we use a series of post-processing network calculations with parametrized expansion profiles that are representative of the wide range of temperatures, densities and electron-to-baryon ratios found in 3D supernova simulations. Critical flows that affect the final yields of 44Ti and 56Ni are assessed.

The goal of this paper is to analyze the impact of a primary neutron source on the s-process nucleosynthesis in massive stars at halo metallicity. Recent stellar models including rotation at very low metallicity predict a strong production of primary N14. Part of the nitrogen produced in the H-burning shell diffuses by rotational mixing into the He core where it is converted to Ne22 providing additional neutrons for the s process. We present nucleosynthesis calculations for a 25 Msun star at [Fe/H] = -3, -4, where in the convective core He-burning about 0.8 % in mass is made of primary Ne22. The usual weak s-process shape is changed by the additional neutron source with a peak between Sr and Ba, where the s-process yields increase by orders of magnitude with respect to the yields obtained without rotation. Iron seeds are fully consumed and the maximum production of Sr, Y and Zr is reached. On the other hand, the s-process efficiency beyond Sr and the ratio Sr/Ba are strongly affected by the amount of Ne22 and by nuclear uncertainties, first of all by the Ne22(alpha,n)Mg25 reaction. Finally, assuming that Ne22 is primary in the considered metallicity range, the s-process efficiency decreases with metallicity due to the effect of the major neutron poisons Mg25 and Ne22. This work represents a first step towards the study of primary neutron source effect in fast rotating massive stars, and its implications are discussed in the light of spectroscopic observations of heavy elements at halo metallicity.

We present detailed models of low and intermediate-mass asymptotic giant branch (AGB) stars with and without the 18F(a,p)21Ne reaction included in the nuclear network, where the rate for this reaction has been recently experimentally evaluated for the first time. The lower and recommended measured rates for this reaction produce negligible changes to the stellar yields, whereas the upper limit of the rate affects the production of 19F and 21Ne. The stellar yields increase by ~50% to up to a factor of 4.5 for 19F, and by factors of ~2 to 9.6 for 21Ne. While the 18}F(a,p)21Ne reaction competes with 18O production, the extra protons released are captured by 18O to facilitate the 18O(p,a)15N(a,g)19F chain. The higher abundances of 19F obtained using the upper limit of the rate helps to match the [F/O] ratios observed in AGB stars, but only for large C/O ratios. Extra-mixing processes are proposed to help to solve this problem. Some evidence that the 18F(a,p)21Ne rate might be closer to its upper limit is provided by the fact that the higher calculated 21Ne/22Ne ratios in the He intershell provide an explanation for the Ne isotopic composition of silicon-carbide grains from AGB stars. This needs to be confirmed by future experiments of the 18F(a,p)21Ne reaction rate. The availability of accurate fluorine yields from AGB stars will be fundamental for interpreting observations of this element in carbon-enhanced metal-poor stars.

The breakout reaction $^{15}$O($\alpha,\gamma$)$^{19}$Ne, which regulates the flow between the hot CNO cycle and the rp-process, is critical for the explanation of the burst amplitude and periodicity of X-ray bursters. We report on the first successful measurement of the critical $\alpha$-decay branching ratios of relevant states in $^{19}$Ne populated via $^{19}$F($^3$He,t)$^{19}$Ne. Based on the experimental results and our previous lifetime measurements of these states, we derive the first experimental rate of $^{15}$O($\alpha,\gamma$)$^{19}$Ne. The impact of our experimental results on the burst pattern and periodicity for a range of accretion rates is analyzed.

Neutron stars in close binary star systems often accrete matter from their companion stars. Thermonuclear ignition of the accreted material in the atmosphere of the neutron star leads to a thermonuclear explosion which is observed as an X-ray burst occurring periodically between hours and days depending on the accretion rate. The ignition conditions are characterized by a sensitive interplay between the accretion rate of the fuel supply and its depletion rate by nuclear burning in the hot CNO cycle and the rp-process. For accretion rates close to stable burning the burst ignition therefore depends critically on the hot CNO breakout reaction, O15(alpha,gamma)Ne19, that regulates the flow between the hot CNO cycle and the rapid proton capture process. Until recently, the O15(alpha,gamma)Ne19-reaction rate was not known experimentally and the theoretical estimates carried significant uncertainties. In this paper we perform a parameter study of the uncertainty of this reaction rate and determine the astrophysical consequences of the first measurement of this reaction rate. Our results corroborate earlier predictions and show that theoretically burning remains unstable up to accretion rates near the Eddington limit, in contrast to astronomical observations.

We analyze thermonuclear and pycnonuclear fusion reactions in dense matter containing atomic nuclei of different types. We extend a phenomenological expression for the reaction rate, proposed recently by Gasques et al. (2005) for the one-component plasma of nuclei, to the multi-component plasma. The expression contains several fit parameters which we adjust to reproduce the best microscopic calculations available in the literature. Furthermore, we show that pycnonuclear burning is drastically affected by an (unknown) structure of the multi-component matter (a regular lattice, a uniform mix, etc.). We apply the results to study nuclear burning in a carbon_12-oxygen_16 mixture. In this context we present new calculations of the astrophysical S-factors for carbon-oxygen and oxygen-oxygen fusion reactions. We show that the presence of a CO lattice can strongly suppress carbon ignition in white dwarf cores and neutron star crusts at densities > 3e9 g cm^{-3} and temperatures T<1e8 K.

The astrophysical p process, which is responsible for the origin of the proton rich stable nuclei heavier than iron, was investigated using a full nuclear reaction network for a type II supernova explosion when the shock front passes through the O/Ne layer. Calculations were performed with a multi-layer model adopting the seed of a pre-explosion evolution of a 25 solar mass star. The reaction flux was calculated to determine the main reaction path and branching points responsible for synthesizing the proton rich nuclei. In order to investigate the impact of nuclear reaction rates on the predicted p-process abundances, extensive simulations with different sets of collectively and individually modified neutron-, proton-, alpha-capture and photodisintegration rates have been performed. These results are not only relevant to explore the nuclear physics related uncertainties in p-process calculations but are also important for identifying the strategy and planning of future experiments.

One of the two breakout reactions from the hot CNOcycle is 15O(a,g)19Ne, which at low temperatures depends strongly on the resonance strength of the 4.033 MeV state in 19Ne. An experimental upper limit has been placed on its strength, but the lower limit on the resonance strength and thereby the astrophysical reaction rate is unconstrained experimentally. However, this breakout reaction is crucial to the thermonuclear runaway which causes type I X-ray bursts on accreting neutron stars. In this paper we exploit astronomical observations in an attempt to constrain the relevant nuclear physics and deduce a lower limit on the reaction rate. Our sensitivity study implies that if the rate were sufficiently small, accreting material would burn stably without bursts. The existence of type I X-ray bursts and superbursts consequently suggests a lower limit on the 15O(a,g)19Ne reaction rate at low temperatures.

We present new rates for the 22Ne(alpha, n)25Mg and 22Ne(alpha,gamma)26Mg reactions, with uncertainties that have been considerably reduced compared to previous estimates, and we study how these new rates affect the production of the heavy magnesium isotopes in models of intermediate mass Asymptotic Giant Branch (AGB) stars of different initial compositions. All the models have deep third dredge-up, hot bottom burning and mass loss. Calculations have been performed using the two most commonly used estimates of the 22Ne + alpha rates as well as the new recommended rates, and with combinations of their upper and lower limits. The main result of the present study is that with the new rates, uncertainties on the production of isotopes from Mg to P coming from the 22Ne + alpha-capture rates have been considerably reduced. We have therefore removed one of the important sources of uncertainty to effect models of AGB stars. We have studied the effects of varying the mass-loss rate on nucleosynthesis and discuss other uncertainties related to the physics employed in the computation of stellar structure, such as the modeling of convection, the inclusion of a partial mixing zone and the definition of convective borders. These uncertainties are found to be much larger than those coming from 22Ne + alpha-capture rates, when using our new estimates. Much effort is needed to improve the situation for AGB models.

Isomeric states in the nuclei along the rapid proton capture process path are studied by the projected shell model. Emphasis is given to two waiting point nuclei Se-68 and Kr-72 that are characterized by shape coexistence. Energy surface calculations indicate that the ground state of these nuclei corresponds to an oblate-deformed minimum, while the lowest state at the prolate-deformed minimum can be considered as a shape isomer. Due to occupation of the orbitals with large K-components, states built upon two-quasiparticle excitations at the oblate-deformed minimum may form high $K$-isomers. The impact of the isomer states on isotopic abundance in X-ray bursts is studied in a multi-mass-zone X-ray burst model by assuming an upper-lower limit approach.

In this paper we analyze the nuclear fusion rate between equal nuclei for all five different nuclear burning regimes in dense matter (two thermonuclear regimes, two pycnonuclear ones, and the intermediate regime). The rate is determined by Coulomb barrier penetration in dense environments and by the astrophysical S-factor at low energies. We evaluate previous studies of the Coulomb barrier problem and propose a simple phenomenological formula for the reaction rate which covers all cases. The parameters of this formula can be varied, taking into account current theoretical uncertainties in the reaction rate. The results are illustrated for the example of the ^{12}C+^{12}C fusion reaction. This reaction is very important for the understanding of nuclear burning in evolved stars, in exploding white dwarfs producing type Ia supernovae, and in accreting neutron stars. The S-factor at stellar energies depends on a reliable fit and extrapolation of the experimental data. We calculate the energy dependence of the S-factor using a recently developed parameter-free model for the nuclear interaction, taking into account the effects of the Pauli nonlocality. For illustration, we analyze the efficiency of carbon burning in a wide range of densities and temperatures of stellar matter with the emphasis on carbon ignition at densities rho > 10^9 g/cc.

In the light of recent recalculations of the \el{19}{F}(alpha,p)\el{22}{Ne} reaction rate we present results of the expected yield of \el{19}{F} from Wolf-Rayet (WR) stars. We have computed models using the upper and lower limits for the rate in addition to the recommended rate and hence we constrain the uncertainty in the yield with respect to this reaction. We find a yield of 3.1x10^{-4}M_sol of \el{19}{F} with our recommended rate and a difference of a factor of two between the yields computed with the upper and lower limits. In comparison with previous work we find a difference in the yield of approximately a factor of 4, connected to a different choice of mass loss. Model uncertainties must be carefully evaluated in order to obtain a reliable estimate of the yield of fluorine from WR stars together with its uncertainties.

Isomeric states in the nuclei along the rapid proton capture process path are studied by the projected shell model. Emphasis is given to two waiting point nuclei 68Se and 72Kr that are characterized by shape coexistence. Energy surface calculations indicate that the ground state of these nuclei corresponds to an oblate-deformed minimum, while the lowest state at the prolate-deformed minimum can be considered as a shape isomer. The impact of these isomer states on isotopic abundance in x-ray bursts is studied in a multi-mass-zone x-ray burst model by assuming an upper-lower limit approach.

This paper reports on the detailed rp-process reaction flow on an accreting neutron star and the resulting ashes of a type I X-ray burst. It is obtained by coupling a 298 isotope reaction network to a self-consistent one-dimensional model calculation with a constant accretion rate of dM/dt=1.0e17g/s (0.09 Eddington).

The astrophysical p-process is responsible for the origin of the proton rich nuclei,which are heavier than iron. A huge network involving thousands of reaction rates is necessary to calculate the final p-abundances. But not all rates included in the network have a strong influence on the p-nuclei abundances. The p-process was investigated using a full nuclear reaction network for a type II supernovae explosion when the shock front passes through the O/Ne layer. Calculations were done with a multi-layer model adopting the seed of a pre-explosion evolution of a 25 mass star. In extensive simulations we investigated the impact of single reaction rates on the final p-abundances. The results are important for the strategy of future experiments in this field.

We present nucleosynthesis calculations and the resulting 19F stellar yields for a large set of models with different masses and metallicity. We find that the production of fluorine depends on the temperature of the convective pulses, the amount of primary 12C mixed into the envelope by third dredge up and the extent of the partial mixing zone. Then we perform a detailed analysis of the reaction rates involved in the production of 19F and the effects of their uncertainties. We find that the major uncertainties are associated with the 14C(alpha,gamma)18O and the 19F(alpha,p)22Ne reaction rates. For these two reactions we present new estimates of the rates and their uncertainties. The importance of the partial mixing zone is reduced when using our estimate for the 14C(alpha,gamma)18O rate. Taking into account both the uncertainties related to the partial mixing zone and those related to nuclear reactions, the highest values of 19F enhancements observed in AGB stars are not matched by the models. This is a problem that will have to be revised by providing a better understanding of the formation and nucleosynthesis in the partial mixing zone, also in relation to reducing the uncertainties of the 14C(alpha,gamma)18O reaction rate. At the same time the possible effect of Cool Bottom Processing at the base of the convective envelope should be included in the computation of AGB nucleosynthesis. This process could in principle help matching the highest 19F abundances observed by decreasing the C/O ratio at the surface of the star, while leaving the 19F abundance unchanged.

Multi-zone models of Type I X-ray bursts are presented that use an adaptive nuclear reaction network of unprecedented size, up to 1300 isotopes. Sequences of up to 15 bursts are followed for two choices of accretion rate and metallicity. At 0.1 Eddington (and 0.02 Eddington for low metallicity), combined hydrogen-helium flashes occur. The rise times, shapes, and tails of these light curves are sensitive to the efficiency of nuclear burning at various waiting points along the rp-process path and these sensitivities are explored. The bursts show "compositional inertia", in that their properties depend on the fact that accretion occurs onto the ashes of previous bursts which contain left-over hydrogen, helium and CNO nuclei. This acts to reduce the sensitivity of burst properties to metallicity. For the accretion rates studied, only the first anomalous burst in one model produces nuclei as heavy as A=100, other bursts make chiefly nuclei with A~64. The amount of carbon remaining after hydrogen-helium bursts is typically <1% by mass, and decreases further as the ashes are periodically heated by subsequent bursts. At the lower accretion rate of 0.02 Eddington and solar metallicity, the bursts ignite in a hydrogen-free helium layer. At the base of this layer, up to 90% of the helium has already burned to carbon prior to the unstable ignition. These helium-ignited bursts have briefer, brighter light curves with shorter tails, very rapid rise times (<0.1 s), and ashes lighter than the iron group.

With the exception of the Big Bang, responsible for 1,2H, 3,4He, and 7Li, stars act as sources for the composition of the interstellar medium. Cosmic rays are related to the latter and very probably due to acceleration of the mixed interstellar medium by shock waves from supernova remnants. Thus, the understanding of the abundance evolution in the interstellar medium and especially the enrichment of heavy elements, as a function of space and time, is essential. It reflects the history of star formation and the lifetimes of the diverse contributing stellar objects. Therefore, the understanding of the endpoints of stellar evolution is essential as well. These are mainly planetary nebulae and type II/Ib/Ic supernovae as evolutionary endpoints of single stars, but also events in binary systems can contribute, like e.g. supernovae of type Ia, novae and possibly X-ray bursts and neutron star or neutron star - black hole mergers. Despite many efforts, a full and self-consistent understanding of supernovae (the main contributors to nucleosynthesis in galaxies) is not existing, yet. Their fingerprints, however, seen either in spectra, lightcurves, radioactivities/decay gamma-rays or in galactic evolution, can help to constrain the composition of their ejecta and related model uncertainties.

The structure of collective states and quasi-particle excitations in $^{68}$Se is analyzed by using the projected shell model. It is found that the $g_{9/2}$ orbitals play a decisive role for the structure of this nucleus. At the prolate minimum, alignment of the low-K states causes the observed first backbending in moment of inertia and we predict a second backbending that has not yet been reported. At the oblate minimum, quasi-particles having high-K can form long-lived isomeric states. The structure of $^{68}$Se is discussed with respect to the isomers and their possible impact on the determination of its nuclear mass.

We calculate the rapid proton (rp) capture process of hydrogen burning on the surface of an accreting neutron star with an updated reaction network that extends up to Xe, far beyond previous work. In both steady-state nuclear burning appropriate for rapidly accreting neutron stars (such as the magnetic polar caps of accreting X-ray pulsars) and unstable burning of Type I X-ray bursts, we find that the rp process ends in a closed SnSbTe cycle. This prevents the synthesis of elements heavier than Te and has important consequences for X-ray burst profiles, the composition of accreting neutron stars, and potentially galactic nucleosynthesis of light p nuclei.

Except for 1H, 2H, 3He, 4He, and 7Li, originating from the Big Bang, all heavier elements are made in stellar evolution and stellar explosions. Nuclear physics, and in many cases nuclear structure far from stability, enters in a crucial way. Therefore, we examine in this review the role of nuclear physics in astrophysics in general and in particular how it affects stellar events and the resulting nucleosynthesis. Stellar modeling addresses four major aspects: 1. energy generation and nucleosynthesis, 2. energy transport via conduction, radiation or possibly convection, 3. hydrodynamics/hydrostatics, and finally 4. thermodynamic properties of the matter involved. Nuclear Physics enters via nuclear reaction cross sections and nuclear structure (affecting the composition changes and nuclear energy generation), neutrino-nucleon and neutrino-nucleus cross sections (affecting neutrino opacities and transport), and e.g. the equation of state at and beyond nuclear densities which creates a relation between the nuclear many-body problem and the hydrodynamic response like pressure and entropy. In the following we review these four topics by highlighting the role and impact of nuclear physics in each of these aspects of stellar modeling. The main emphasis is put on the connection to element synthesis.

The hydrogen-rich accreted envelopes of neutron stars in binary systems are the site for the rapid proton capture nucleosynthethic process (rp-process), which involves nuclei close to the proton dripline. An overview of the relevant reactions and nuclear properties for rp-process studies is given, along with motivations for further experimental nuclear physics studies.

The reaction rates for alpha capture processes on self-conjugate nuclei in the mass range A=20-40 have been investigated. The rates were calculated using the statistical model code NON-SMOKER taking into account isospin suppression rules. These theoretical predictions are compared with rates derived from the available experimental data about the alpha capture reactions taking also into account additional experimental information from different reaction channels populating the alpha unbound states of the self-conjugate compound nuclei.

The reaction rates of neutron capture reactions on light nuclei are important for reliably simulating nucleosynthesis in a variety of stellar scenarios. Neutron capture reaction rates on neutron-rich C-, N-, and O-isotopes are calculated in the framework of a hybrid compound and direct capture model. The results are tabulated and compared with the results of previous calculations as well as with experimental results.

We calculate the nucleosynthesis during stable nuclear burning on an accreting neutron star. This is appropriate for weakly magnetic neutron stars accreting at near-Eddington rates in low mass X-ray binaries, and for most accreting X-ray pulsars. We show that the nuclear burning proceeds via the rapid proton capture process (rp process), and makes nuclei far beyond the iron group. The final mixture of nuclei consists of elements with a range of masses between approximately A=60 and A=100. The average nuclear mass of the ashes is set by the extent of helium burning via (alpha,p) reactions, and depends on the local accretion rate. Our results imply that the crust of these accreting neutron stars is made from a complex mixture of heavy nuclei, with important implications for its thermal, electrical and structural properties. A crustal lattice as impure as our results suggest will have a conductivity set mostly by impurity scattering, allowing more rapid Ohmic diffusion of magnetic fields than previously estimated.

Hot (explosive) hydrogen burning or the Rapid Proton Capture Process (rp-process) occurs in a number of astrophysical environments. Novae and X-ray bursts are the most prominent ones, but accretion disks around black holes and other sites are candidates as well. The expensive and often multidimensional hydro calculations for such events require an accurate prediction of the thermonuclear energy generation, while avoiding full nucleosynthesis network calculations. In the present investigation we present an approximation scheme applicable in a temperature range which covers the whole range of all presently known astrophysical sites. It is based on the concept of slowly varying hydrogen and helium abundances and assumes a kind of local steady flow by requiring that all reactions entering and leaving a nucleus add up to a zero flux. This scheme can adapt itself automatically and covers situations at low temperatures, characterized by a steady flow of reactions, as well as high temperature regimes where a $(p,\gamma)-(\gamma,p)$-equilibrium is established. In addition to a gain of a factor of 15 in computational speed over a full network calculation, and an energy generation accurate to more than 15 %, this scheme also allows to predict correctly individual isotopic abundances. Thus, it delivers all features of a full network at a highly reduced cost and can easily be implemented in hydro calculations.