Code dependencies of pre-supernova evolution and nucleosynthesis in massive stars: Evolution to the end of core helium burning

Massive stars are key sources of radiative, kinetic, and chemical feedback in the universe. Grids of massive star models computed by different groups each using their own codes, input physics choices and numerical approximations, however, lead to inconsistent results for the same stars. We use three of these 1D codes---GENEC, KEPLER and MESA---to compute non-rotating stellar models of $15~mathrm{M}_odot$, $20~mathrm{M}_odot$, and $25~mathrm{M}_odot$ and compare their nucleosynthesis. We follow the evolution from the main sequence until the end of core helium burning. The GENEC and KEPLER models hold physics assumptions used in large grids of published models. The MESA code was set up to use convective core overshooting such that the CO core masses are consistent with those obtained by GENEC. For all models, full nucleosynthesis is computed using the NuGrid post-processing tool MPPNP. We find that the surface abundances predicted by the models are in reasonable agreement. In the helium core, the standard deviation of the elemental overproduction factors for Fe to Mo is less than $30,%$---smaller than the impact of the present nuclear physics uncertainties. For our three initial masses, the three stellar evolution codes yield consistent results. Differences in key properties of the models, e.g., helium and CO core masses and the time spent as a red supergiant, are traced back to the treatment of convection and, to a lesser extent, mass loss. The mixing processes in stars remain the key uncertainty in stellar modelling. Better constrained prescriptions are thus necessary to improve the predictive power of stellar evolution models.

Relics of ancient post-AGB stars in a primitive meteorite

Graphite is one of the many presolar circumstellar condensate species found in primitive meteorites. While the isotopic compositions of low-density graphite grains indicate an origin in core-collapse supernovae, some high-density grains have extreme isotopic anomalies in C, Ca and Ti, which cannot be explained by envelope predictions of asymptotic giant branch (AGB) stars or theoretical supernova models. The Ca and Ti isotopic anomalies, however, match the predictions of He-shell abundances in AGB stars. In this study, we show that the C, Ca, and Ti isotopic anomalies are consistent with nucleosynthesis predictions of the H-ingestion phase during a very late thermal pulse (VLTP) event in post-AGB stars. The low $^{12}$C/$^{13}$C isotopic ratios in these grains are a result of abundant $^{12}$C efficiently capturing the protons that are being ingested during the VLTP. Very high neutron densities of $sim 10^{15}$ cm$^{-3}$, typical of the $i$-process, are achieved during this phase in post-AGB stars. The large $^{42,43,44}$Ca excesses in some graphite grains are indicative of neutron capture nucleosynthesis during VLTP. The comparison of VLTP nucleosynthesis calculations to the graphite data also indicate that apparent anomalies in the Ti isotopic ratios are due to large contributions from $^{46,48}$Ca, which cannot be resolved from the isobars $^{46,48}$Ti during the measurements. We conclude that presolar graphite grains with moderate to extreme Ca and Ti isotopic anomalies originate in post-AGB stars that suffer a very late thermal pulse.

Convective-reactive proton-C12 combustion in Sakurais object (V4334 Sagittarii) and implications for the evolution and yields from the first generations of stars

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 a He-burning zone, for example in convectively unstable shell on top of electron-degenerate cores in AGB 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 zero metal content [...] We focus on the convective-reactive episode in the very-late thermal pulse star Sakurais object (V4334 Sagittarii). Asplund etal. (1999) 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 that are too low to obtain a significant neutron capture nucleosynthesis on the heavy elements. We have carried out 3D hydrodynamic He-shell flash convection [...] 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 [...] we obtain significantly higher neutron densities (~few 10^15 1/cm^3) and reproduce the key observed abundance trends found in Sakurais object. These include an overproduction of Rb, Sr and Y by about 2 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. [...] We determine how our results depend on uncertainties of nuclear reaction rates, for example for the C13(alpha, n)O16 reaction.

Complete nucleosynthesis calculations for low-mass stars from NuGrid

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.