Stardust grains that originated in ancient stars and supernovae are recovered from meteorites and carry the detailed composition of their astronomical sites of origin. We present evidence that the majority of large ($\mu$m-sized) meteoritic silicon carbide (SiC) grains formed in C-rich asymptotic giant branch (AGB) stars that were more metal-rich than the Sun. In the framework of the slow neutron-captures (the s process) that occurs in AGB stars the lower-than-solar 88Sr/86Sr isotopic ratios measured in the large SiC grains can only be accompanied by Ce/Y elemental ratios that are also lower than solar, and predominately observed in metal-rich barium stars - the binary companions of AGB stars. Such an origin suggests that these large grains represent the material from high-metallicity AGB stars needed to explain the s-process nucleosynthesis variations observed in bulk meteorites (Ek et al. 2020). In the outflows of metal-rich, C-rich AGB stars SiC grains are predicted to be small ($\simeq$ 0.2 $\mu$m-sized); large ($\simeq$ $\mu$m-sized) SiC grains can grow if the number of dust seeds is two to three orders of magnitude lower than the standard value of $10^{-13}$ times the number of H atoms. We therefore predict that with increasing metallicity the number of dust seeds might decrease, resulting in the production of larger SiC grains.

The competing $^{22}$Ne($\alpha,\gamma$)$^{26}$Mg and $^{22}$Ne($\alpha,n$)$^{25}$Mg reactions control the production of neutrons for the weak $s$-process in massive and AGB stars. In both systems, the ratio between the corresponding reaction rates strongly impacts the total neutron budget and strongly influences the final nucleosynthesis. The $^{22}$Ne($\alpha,\gamma$)$^{26}$Mg and $^{22}$Ne($\alpha,n$)$^{25}$Mg reaction rates was re-evaluated by using newly available information on $^{26}$Mg given by various recent experimental studies. Evaluations of The evaluated $^{22}$Ne($\alpha,\gamma$)$^{26}$Mg reaction rate remains substantially similar to that of Longland {\it et al.} but, including recent results from Texas A\&M, the $^{22}$Ne($\alpha,n$)$^{25}$Mg reaction rate is lower at a range of astrophysically important temperatures. Stellar models computed with NEWTON and MESA predict decreased production of the weak branch $s$-process due to the decreased efficiency of $^{22}$Ne as a neutron source. Using the new reaction rates in the MESA model results in $^{96}$Zr/$^{94}$Zr and $^{135}$Ba/$^{136}$Ba ratios in much better agreement with the measured ratios from presolar SiC grains.

Aims. We investigate the s-process during the AGB phase of stellar models whose cores are enforced to rotate at rates consistent with asteroseismology observations of their progenitors and successors. Methods. We calculated new 2M$_{\odot}$, Z=0.01 models, rotating at 0, 125, and 250 km/s at the start of main sequence. An artificial, additional viscosity was added to enhance the transport of angular momentum in order to reduce the core rotation rates to be in agreement with asteroseismology observations. We compared rotation rates of our models with observed rotation rates during the MS up to the end of core He burning, and the white dwarf phase. Results. We present nucleosynthesis calculations for these rotating AGB models that were enforced to match the asteroseismic constraints on rotation rates of MS, RGB, He-burning, and WD stars. In particular, we calculated one model that matches the upper limit of observed rotation rates of core He-burning stars and we also included a model that rotates one order of magnitude faster than the upper limit of the observations. The s-process production in both of these models is comparable to that of non-rotating models. Conclusions. Slowing down the core rotation rate in stars to match the above mentioned asteroseismic constraints reduces the rotationally induced mixing processes to the point that they have no effect on the s-process nucleosynthesis. This result is independent of the initial rotation rate of the stellar evolution model. However, there are uncertainties remaining in the treatment of rotation in stellar evolution, which need to be reduced in order to confirm our conclusions, including the physical nature of our approach to reduce the core rotation rates of our models, and magnetic processes.

The production of the neutron-capture isotopes beyond iron that we observe today in the solar system is the result of the combined contribution of the r-process, the s- process and possibly the i-process. Low-mass AGB (2 < M/Msun < 3) and massive (M >10 Msun ) stars have been identified as the sites of the s-process. In this work we consider the evolution and nucleosynthesis of low-mass AGB stars. We provide an update of the NuGrid Set models, adopting the same general physics assumptions but using an updated convective-boundary mixing model accounting for the contribution from internal gravity waves. The combined data set includes the initial masses Mzams/Msun = 2, 3 for Z = 0.03, 0.02, 0.01. These models are computed with the MESA stellar code and the evolution is followed up to the end of the AGB phase. The nucleosynthesis was calculated for all isotopes in post-processing with the NuGrid mppnp code. The convective boundary mixing model leads to the formation of a 13C-pocket three times wider compared to the one obtained in the previous set of models, bringing the simulation results now in closer agreement with observations. We also discuss the potential impact of other processes inducing mixing, like rotation, adopting parametric models compatible with theory and observations. Complete yield data tables, derived data products and online analytic data access are provided.

Based on stellar evolution simulations, we demonstrate that rapidly accreting white dwarfs in close binary systems are an astrophysical site for the intermediate neutron-capture process. During recurrent and very strong He-shell flashes in the stable H-burning accretion regime H-rich material enters the He-shell flash convection zone. $^{12}$C(p,$\gamma)^{13}$N reactions release enough energy to potentially impact convection, and i process is activated through the $^{13}$C($\alpha$,n)$^{16}$O reaction. The H-ingestion flash may not cause a split of the convection zone as it was seen in simulations of He-shell flashes in post-AGB and low-Z AGB stars. We estimate that for the production of first-peak heavy elements this site can be of similar importance for galactic chemical evolution as the s-process production by low-mass AGB stars. The He-shell flashes result in the expansion and, ultimately, ejection of the accreted and then i-process enriched material, via super-Eddington luminosity winds or Roche-lobe overflow. The white dwarf models do not retain any significant amount of the accreted mass, with a He retention efficiency of $\leq 10\%$ depending on mass and convective boundary mixing assumptions. This makes the evolutionary path of such systems to supernova Ia explosion highly unlikely.

The $s$-process nucleosynthesis in Asymptotic Giant Branch (AGB) stars depends on the modeling of convective boundaries. We present models and s-process simulations that adopt a treatment of convective boundaries based on the results of hydrodynamic simulations and on the theory of mixing due to gravity waves in the vicinity of convective boundaries. Hydrodynamics simulations suggest the presence of convective boundary mixing (CBM) at the bottom of the thermal pulse-driven convective zone. Similarly, convection-induced mixing processes are proposed for the mixing below the convective envelope during third dredge-up where the 13C pocket for the s process in AGB stars forms. In this work we apply a CBM model motivated by simulations and theory to models with initial mass $M = 2$ and $M = 3M_\odot$, and with initial metal content Z = 0.01 and Z = 0.02. As reported previously, the He-intershell abundance of 12C and 16O are increased by CBM at the bottom of pulse-driven convection zone. This mixing is affecting the $^{22}Ne(\alpha,n)^{25}Mg$ activation and the s-process effciency in the 13C-pocket. In our model CBM at the bottom of the convective envelope during the third dredgeup represents gravity wave mixing. We take further into account that hydrodynamic simulations indicate a declining mixing efficiency already about a pressure scale height from the convective boundaries, compared to mixing-length theory. We obtain the formation of the 13C-pocket with a mass of $\approx 10^{-4}M_\odot$. The final $s$-process abundances are characterized by 0.36 < [s=Fe] < 0.78 and the heavy-to-light s-process ratio is 0.23 < [hs=ls] < 0.45. Finally, we compare our results with stellar observations, pre-solar grain measurements and previous work.

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.

Classical novae are the result of thermonuclear flashes of hydrogen accreted by CO or ONe white dwarfs, leading eventually to the dynamic ejection of the surface layers. These are observationally known to be enriched in heavy elements, such as C, O and Ne that must originate in layers below the H-flash convection zone. Building on our previous work, we now present stellar evolution simulations of ONe novae and provide a comprehensive comparison of our models with published ones. Some of our models include exponential convective boundary mixing to account for the observed enrichment of the nova ejecta even when accreted material has a solar abundance distribution. Our models produce maximum temperature evolution profiles and nucleosynthesis yields in good agreement with models that generate enriched ejecta by assuming that the accreted material was pre-mixed. We confirm for ONe novae the result we reported previously, i.e.\ we found that $^3$He could be produced {\it in situ} in solar-composition envelopes accreted with slow rates ($\dot{M} < 10^{-10}\,M_\odot/\mbox{yr}$) by cold ($T_{\rm WD} < 10^7$ K) CO WDs, and that convection was triggered by $^3$He burning before the nova outburst in that case. In addition, we now find that the interplay between the $^3$He production and destruction in the solar-composition envelope accreted with an intermediate rate, e.g.\ $\dot{M} = 10^{-10}\,M_\odot/\mbox{yr}$, by the $1.15\,M_\odot$ ONe WD with a relatively high initial central temperature, e.g.\ $T_{\rm WD} = 15\times 10^6$ K, leads to the formation of a thick radiative buffer zone that separates the bottom of the convective envelope from the WD surface. (Abridged)