We provide a set of stellar evolution and nucleosynthesis calculations that applies established physics assumptions simultaneously to low- and intermediate-mass and massive star models. Our goal is to provide an internally consistent and comprehensive nuclear production and yield data base for applications in areas such as pre-solar grain studies. Our non-rotating models assume convective boundary mixing where it has been adopted before. We include 8 (12) initial masses for $Z = 0.01$ ($0.02$). Models are followed either until the end of the asymptotic giant branch phase or the end of Si burning, complemented by a simple analytic core-collapse supernova models with two options for fallback and shock velocities. The explosions show which pre-supernova yields will most strongly be effected by the explosive nucleosynthesis. We discuss how these two explosion parameters impacts the light elements and the $s$ and $p$ process. For low- and intermediate-mass models our stellar yields from H to Bi include the effect of convective boundary mixing at the He-intershell boundaries and the stellar evolution feedback of the mixing process that produces the $^{13}$C pocket. All post-processing nucleosynthesis calculations use the same nuclear reaction rate network and nuclear physics input. We provide a discussion of the nuclear production across the entire mass range organized by element group. All our stellar nucleosynthesis profile and time evolution output is available electronically, and tools to explore the data on the NuGrid VOspace hosted by the Canadian Astronomical Data Centre are introduced.

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.

We have measured the neutron capture cross sections of the stable magnesium isotopes Mg-24,Mg-25,Mg-26 in the energy range of interest to the s process using the neutron time-of-flight facility n_TOF at CERN. Capture events from a natural metal sample and from samples enriched in Mg-25 and Mg-26 were recorded using the total energy method based on (C6H6)-H-2 detectors. Neutron resonance parameters were extracted by a simultaneous resonance shape analysis of the present capture data and existing transmission data on a natural isotopic sample. Maxwellian-averaged capture cross sections for the three isotopes were calculated up to thermal energies of 100 keV and their impact on s-process analyses was investigated. At 30 keV the new values of the stellar cross section for Mg-24, Mg-25, and Mg-26 are 3.8 +/- 0.2 mb, 4.1 +/- 0.6 mb, and 0.14 +/- 0.01 mb, respectively.

Over the last 40 years, the 12C +12C fusion reaction has been the subject of considerable experimental efforts to constrain uncertainties at temperatures relevant for stellar nucleosynthesis. Recent studies have indicated that the reaction rate may be higher than that currently used in stellar models. 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 considering five different initial masses, 15, 20, 25, 32 and 60 M⊙, at solar metallicity, were generated using the Geneva Stellar Evolution Code (genec) and were later post-processed with the NuGrid Multi-zone Post-Processing Network tool (mppnp). A dynamic nuclear reaction network of ∼1100 isotopes was used to track the s-process nucleosynthesis. An enhanced 12C +12C reaction 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 (rather than a radiative one). 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 α- and p-exit channels.

Depending on mass and metallicity as well as evolutionary phase, stars occasionally experience convective-reactive nucleosynthesis episodes. We specifically investigate the situation when nucleosynthetically unprocessed, H-rich material is convectively mixed with an He-burning zone, for example in a convectively unstable shell on top of electron-degenerate cores in asymptotic giant branch stars, young white dwarfs, or X-ray bursting neutron stars. Such episodes are frequently encountered in stellar evolution models of stars of extremely low or zerometal content, such as the first stars. We have carried out detailed nucleosynthesis simulations based on stellar evolution models and informed by hydrodynamic simulations. We focus on the convective-reactive episode in the very late thermal pulse star Sakurai's object (V4334 Sagittarii). Asplund et al. determined the abundances of 28 elements, many of which are highly non-solar, ranging from H, He, and Li all the way to Ba and La, plus the C isotopic ratio. Our simulations show that the mixing evolution according to standard, one-dimensional stellar evolution models implies neutron densities in the He intershell (less than or similar to few 10(11) cm(-3)) that are too low to obtain a significant neutron capture nucleosynthesis on the heavy elements. We have carried out three-dimensional hydrodynamic He-shell flash convection simulations in 4 pi geometry to study the entrainment of H-rich material. Guided by these simulations we assume that the ingestion process of H into the He-shell convection zone leads only after some delay time to a sufficient entropy barrier that splits the convection zone into the original one driven by He burning and a new one driven by the rapid burning of ingested H. By making such mixing assumptions that are motivated by our hydrodynamic simulations we obtain significantly higher neutron densities (similar to few 10(15) cm(-3)) and reproduce the key observed abundance trends found in Sakurai's object. These include an overproduction of Rb, Sr, and Y by about two orders of magnitude higher than the overproduction of Ba and La. Such a peculiar nucleosynthesis signature is impossible to obtain with the mixing predictions in our one-dimensional stellar evolution models. The simulated Li abundance and the isotopic ratio C-12/C-13 are as well in agreement with observations. Details of the observed heavy element abundances can be used as a sensitive diagnostic tool for the neutron density, for the neutron exposure and, in general, for the physics of the convective-reactive phases in stellar evolution. For example, the high elemental ratio Sc/Ca and the high Sc production indicate high neutron densities. The diagnostic value of such abundance markers depends on uncertain nuclear physics input. We determine how our results depend on uncertainties of nuclear reaction rates, for example for the C-13(alpha, n)O-16 reaction.

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.

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.

The "collapsar" engine for gamma-ray bursts invokes as its energy source the failure of a normal supernova and the formation of a black hole.Here we present the results of the first threedimensional simulation of the collapse of a massive star down to a black hole, including the subsequent accretion and explosion.The explosion differs significantly from the axisymmetric scenario obtained in two-dimensional simulations; this has important consequences for the nucleosynthetic yields.We compare the nucleosynthetic yields to those of hypernovae.Calculating yields from three-dimensional explosions requires new strategies in post-process nucleosynthesis; we discuss NuGrid's plan for three-dimensional yields.

We present preliminary results from recent high-resolution double-degenerate merger simulations with the Smooth Particle Hydrodynamics (SPH) technique.We put particular emphasis on verification and validation in our effort and show the importance of details in the initial condition setup for the final outcome of the simulation.We also stress the dynamical importance of including shocks in the simulations.These results represent a first step toward a suite of simulations that will shed light on the question whether double-degenerate mergers are a viable path toward type 1a supernovae.In future simulations, we will make use of the capabilities of the NuGrid collaboration in post-processing SPH particle trajectories with a complete nuclear network to follow the detailed nuclear reactions during the dynamic merger phase.

Many nucleosynthesis and mixing processes of low-mass stars as they evolve from the Main Sequence to the thermal-pulse Asymptotic Giant Branch phase (TP-AGB) are well understood (although of course important physics components, e.g. rotation, magnetic fields, gravity wave mixing, remain poorly known). Nevertheless, in the last years presolar grain measurements with high resolution have presented new puzzling problems and strong constraints on nucleosynthesis processes in stars. The goal of the NuGrid collaboration is to present uniform yields for a large range of masses and metallicities, including low-mass stars and massive stars and their explosions. Here we present the first calculations of stellar evolution and high-resolution, post-processing simulations of an AGB star with an initial mass of 2 M_sun and solar-like metallicity (Z=0.01), based on the post-processing code PPN. In particular, we analyze the formation and evolution of the radiative 13C-pocket between the 17th TP and the 18th TP. The s-process nucleosynthesis profile of a sample of heavy isotopes is also discussed, before the next convective TP occurrence.

We examine observed heavy element abundances in the Cassiopeia A supernova remnant as a constraint on the nature of the Cas A supernova. We compare bulk abundances from 1D and 3D explosion models and spatial distribution of elements in 3D models with those derived from X-ray observations. We also examine the cospatial production of 26Al with other species. We find that the most reliable indicator of the presence of 26Al in unmixed ejecta is a very low S/Si ratio ( 0.05). Production of N in O/S/Si-rich regions is also indicative. The biologically important element P is produced at its highest abundance in the same regions. Proxies should be detectable in supernova ejecta with high spatial resolution multiwavelength observations.

We investigate the physical conditions where 44 Ti and 56 Ni are created in core-collapse supernovae.In this preliminary work we use a series of post-processing network calculations with parameterized 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 44 Ti and 56 Ni are assessed.

We review some of the uncertainties in calculating nucleosynthetic yields, focusing on the explosion mechanism.Current yield calculations tend to either use a piston, energy injection, or enhancement of neutrino opacities to drive an explosion.We show that the energy injection, or more accurately, an entropy injection mechanism is best-suited to mimic our current understanding of the convection-enhanced supernova engine.The enhanced neutrino-opacity technique is in qualitative disagreement with simulations of core-collapse supernovae and will likely produce errors in the yields.But piston-driven explosions are the most discrepant.Piston-driven explosion severely underestimate the amount of fallback, leading to order-of-magnitude errors in the yields of heavy elements.To obtain yields accurate to the factor of a few level, we must use entropy or energy injection and this has become the NuGrid collaboration approach.

The nucleosynthetic yield from a supernova explosion depends upon a variety of effects: progenitor evolution, explosion process, details of the nuclear network, and nuclear rates.Especially in studies of integrated stellar yields, simplifications reduce these uncertainties.But nature is much more complex, and to actually study nuclear rates, we will have to understand the full, complex set of processes involved in nucleosynthesis.Here we discuss a few of these complexities and detail how the NuGrid collaboration will address them.

Simulations of nucleosynthesis in astrophysical environments are at the intersection of nuclear physics reaction rate research and astrophysical applications, for example in the area of galactic chemical evolution or near-field cosmology. Unfortunately, at present the available yields for such applications are based on heterogeneous assumptions between the various contributing nuclear production sites, both in terms of modeling the thermodynamic environment itself as well as the choice of specifc nuclear reaction rates and compilations. On the other side, new nuclear reaction rate determinations are often taking a long time to be included in astrophysical applications. The NuGrid project addresses these issues by providing a set of codes and a framework in which these codes interact. In this contribution we describe the motivation, goals and first results of the NuGrid project. At the core is a new and evolving post-processing nuclesoynthesis code (PPN) that can follow quiescent and explosive nucleosynthesis following multi-zone 1D-stellar evolution as well as multi-zone hydrodynamic input, including explosions. First results are available in the areas of AGB and massive stars.

R. Hirschi abc, U. Frischknechtd, M. Pignatariabe, F.-K. Thielemannd, M. E. Bennettab, S. Diehla f g, C. L. Fryera f , F. Herwigabh, A. Hungerfordag, G. Magkotsiosaei, G. Rockefellerag, F. X. Timmesai, and P. Youngai aThe NuGrid Collaboration bAstrophysics Group, Keele University, ST5 5BG, UK cIPMU, University of Tokyo, Kashiwa, Chiba 277-8582, Japan dTheoretical Astrophysics Group, University of Basel, Basel, 4056, Switzerland eJoint Institute for Nuclear Astrophysics, University of Notre Dame, IN, 46556, USA f Theoretical Astrophysics Group (T-6), Los Alamos National Laboratory, Los Alamos, NM, 87544, USA gComputational Methods (CCS-2), Los Alamos National Laboratory, Los Alamos, NM, 87544, USA hDept. of Physics & Astronomy, Victoria, BC, V8W 3P6, Canada iSchool of Earth and Space Exploration, Arizona State University, Tempe, AZ 85287, USA E-mail:r.hirschi@epsam.keele.ac.uk

We present preliminary results from recent high-resolution double-degenerate merger simulations with the Smooth Particle Hydrodynamics (SPH) technique. We put particular emphasis on verification and validation in our effort and show the importance of details in the initial condition setup for the final outcome of the simulation. We also stress the dynamical importance of including shocks in the simulations. These results represent a first step toward a suite of simulations that will shed light on the question whether double-degenerate mergers are a viable path toward type 1a supernovae. In future simulations, we will make use of the capabilities of the NuGrid collaboration in post-processing SPH particle trajectories with a complete nuclear network to follow the detailed nuclear reactions during the dynamic merger phase.

We review some of the uncertainties in calculating nucleosynthetic yields, focusing on the explosion mechanism. Current yield calculations tend to either use a piston, energy injection, or enhancement of neutrino opacities to drive an explosion. We show that the energy injection, or more accurately, an entropy injection mechanism is best-suited to mimic our current understanding of the convection-enhanced supernova engine. The enhanced neutrino-opacity technique is in qualitative disagreement with simulations of core-collapse supernovae and will likely produce errors in the yields. But piston-driven explosions are the most discrepant. Piston-driven explosion severely underestimate the amount of fallback, leading to order-of-magnitude errors in the yields of heavy elements. To obtain yields accurate to the factor of a few level, we must use entropy or energy injection and this has become the NuGrid collaboration approach.