▪ Abstract The current status of modeling the evolution and nucleosynthesis of asymptotic giant branch (AGB) stars is reviewed. The principles of AGB evolution have been investigated in recent years ...leading to improved and refined models, for example with regard to hot-bottom burning or the third dredge-up. The postprocessing s-process model yields quantitative results that reproduce many observations. However, these and most other processes in AGB stars are intimately related to the physics of stellar mixing. Mixing in AGB stars is currently not well-enough understood for accurate yield predictions. Several constraints and methods are available to improve the models. Some regimes of AGB evolution have not yet been studied in sufficient detail. These include the super-AGB stars and AGB stars at extremely low or ultra low metallicity.
Full text
Available for:
CMK, DOBA, FMFMET, IZUM, KILJ, NUK, PILJ, PNG, SAZU, UILJ, UKNU, UL, UM, UPUK
Abstract
We update the capabilities of the open-knowledge software instrument Modules for Experiments in Stellar Astrophysics (
MESA
). The new
auto
_
diff
module implements automatic differentiation ...in
MESA
, an enabling capability that alleviates the need for hard-coded analytic expressions or finite-difference approximations. We significantly enhance the treatment of the growth and decay of convection in
MESA
with a new model for time-dependent convection, which is particularly important during late-stage nuclear burning in massive stars and electron-degenerate ignition events. We strengthen
MESA
’s implementation of the equation of state, and we quantify continued improvements to energy accounting and solver accuracy through a discussion of different energy equation features and enhancements. To improve the modeling of stars in
MESA
, we describe key updates to the treatment of stellar atmospheres, molecular opacities, Compton opacities, conductive opacities, element diffusion coefficients, and nuclear reaction rates. We introduce treatments of starspots, an important consideration for low-mass stars, and modifications for superadiabatic convection in radiation-dominated regions. We describe new approaches for increasing the efficiency of calculating monochromatic opacities and radiative levitation, and for increasing the efficiency of evolving the late stages of massive stars with a new operator-split nuclear burning mode. We close by discussing major updates to
MESA
’s software infrastructure that enhance source code development and community engagement.
We present a detailed analysis of the composition and nucleosynthetic origins of the heavy elements in the metal-poor (Fe/H = -1.62 + or - 0.09) star HD 94028. Previous studies revealed that this ...star is mildly enhanced in elements produced by the slow neutron-capture process (s process; e.g., Pb/Fe = +0.79 + or - 0.32) and rapid neutron-capture process (r process; e.g., Eu/Fe = +0.22 + or - 0.12), including unusually large molybdenum (Mo/Fe = +0.97 + or - 0.16) and ruthenium (Ru/Fe = +0.69 + or - 0.17) enhancements. However, this star is not enhanced in carbon (C/Fe = -0.06 + or - 0.19). We analyze an archival near-ultraviolet spectrum of HD 94028, collected using the Space Telescope Imaging Spectrograph on board the Hubble Space Telescope, and other archival optical spectra collected from ground-based telescopes. We report abundances or upper limits derived from 64 species of 56 elements. We compare these observations with s-process yields from low-metallicity AGB evolution and nucleosynthesis models. No combination of s- and r-process patterns can adequately reproduce the observed abundances, including the super-solar As/Ge ratio (+0.99 + or - 0.23) and the enhanced Mo/Fe and Ru/Fe ratios. We can fit these features when including an additional contribution from the intermediate neutron-capture process (i process), which perhaps operated through the ingestion of H in He-burning convective regions in massive stars, super-AGB stars, or low-mass AGB stars. Currently, only the i process appears capable of consistently producing the super-solar As/Ge ratios and ratios among neighboring heavy elements found in HD 94028. Other metal-poor stars also show enhanced As/Ge ratios, hinting that operation of the i process may have been common in the early Galaxy.
We use the OMEGA galactic chemical evolution code to investigate how the assumptions used for the treatment of galactic inflows and outflows impact numerical predictions. The goal is to determine how ...our capacity to reproduce the chemical evolution trends of a galaxy is affected by the choice of implementation used to include those physical processes. In pursuit of this goal, we experiment with three different prescriptions for galactic inflows and outflows and use OMEGA within a Markov Chain Monte Carlo code to recover the set of input parameters that best reproduces the chemical evolution of nine elements in the dwarf spheroidal galaxy Sculptor. This provides a consistent framework for comparing the best-fit solutions generated by our different models. Despite their different degrees of intended physical realism, we found that all three prescriptions can reproduce in an almost identical way the stellar abundance trends observed in Sculptor. This result supports the similar conclusions originally claimed by Romano & Starkenburg for Sculptor. While the three models have the same capacity to fit the data, the best values recovered for the parameters controlling the number of SNe Ia and the strength of galactic outflows, are substantially different and in fact mutually exclusive from one model to another. For the purpose of understanding how a galaxy evolves, we conclude that only reproducing the evolution of a limited number of elements is insufficient and can lead to misleading conclusions. More elements or additional constraints such as the Galaxy's star-formation efficiency and the gas fraction are needed in order to break the degeneracy between the different modeling assumptions. Our results show that the successes and failures of chemical evolution models are predominantly driven by the input stellar yields, rather than by the complexity of the Galaxy model itself. Simple models such as OMEGA are therefore sufficient to test and validate stellar yields. OMEGA is part of the NuGrid chemical evolution package and is publicly available online at http://nugrid.github.io/NuPyCEE.
Based on stellar evolution simulations, we demonstrate that rapidly accreting white dwarfs (WDs) 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. reactions release enough energy to potentially impact convection, and i process is activated through the 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 asymptotic giant branch (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 WD models do not retain any significant amount of the accreted mass, with a He retention efficiency of depending on mass and convective boundary mixing assumptions. This makes the evolutionary path of such systems to supernova Ia explosion highly unlikely.
The Stellar Yields for Galactic Modeling Applications (SYGMA) code is an open-source module that models the chemical ejecta and feedback of simple stellar populations (SSPs). It is intended for use ...in hydrodynamical simulations and semi-analytic models of galactic chemical evolution. The module includes the enrichment from asymptotic giant branch (AGB) stars, massive stars, Type Ia supernovae (SNe Ia), and compact binary mergers. An extensive and extendable stellar yields library includes the NuGrid yields with all elements and many isotopes up to Bi. Stellar feedback from mechanic and frequency-dependent radiative luminosities are computed based on NuGrid stellar models and their synthetic spectra. The module further allows for customizable initial mass functions and SN Ia delay-time distributions to calculate time-dependent ejecta based on stellar yield input. A variety of r-process sites can be included. A comparison of SSP ejecta based on NuGrid yields with those from Portinari et al. and Marigo reveals up to factors of 3.5 and 4.8 less C and N enrichment from AGB stars at low metallicity, a result we attribute to NuGrid's modeling of hot-bottom burning. Different core-collapse supernova explosions and fallback prescriptions may lead to substantial variations for the accumulated ejecta of C, O and Si in the first 107 years at Z = 0.001. An online interface of the open-source SYGMA module enables interactive simulations, analysis, and data extraction of the evolution of all species formed by the evolution of simple stellar populations.
Rapidly accreting white dwarfs (RAWDs) have been proposed as contributors to the chemical evolution of heavy elements in the Galaxy. Here, we test this scenario for the first time and determine the ...contribution of RAWDs to the solar composition of first-peak neutron-capture elements. We add the metallicity-dependent contribution of RAWDs to the one-zone galactic chemical evolution code OMEGA according to RAWD rates from binary stellar population models combined with metallicity-dependent i-process stellar yields calculated following the models of Denissenkov et al. With this approach, we find that the contribution of RAWDs to the evolution of heavy elements in the Galaxy could be responsible for a significant fraction of the solar composition of Kr, Rb, Sr, Y, Zr, Nb, and Mo ranging from 2% to 45% depending on the element, the enrichment history of the Galactic gas, and the total mass ejected per RAWD. This contribution could explain the missing solar Lighter Element Primary Process for some elements (e.g., Sr, Y, and Zr). We do not overproduce any isotope relative to the solar composition, but 96Zr is produced in a similar amount. The i process produces efficiently the Mo stable isotopes 95Mo and 97Mo. When nuclear reaction rate uncertainties are combined with our GCE uncertainties, the upper limits for the predicted RAWD contribution increase by a factor of 1.5-2 for Rb, Sr, Y, and Zr, and by 3.8 and 2.4 for Nb and Mo, respectively. We discuss the implication of the RAWD stellar evolution properties on the single-degenerate SN Ia scenario.
We analyze the production of the element Cr in galactic chemical evolution (GCE) models using the NuGrid nucleosynthesis yields set. We show that the unusually large Cr/Fe abundance at Fe/H 0 ...reported by previous studies using those yields and predicted by our Milky Way model originates from the merging of convective Si-burning and C-burning shells in a 20 model at metallicity Z = 0.01, about an hour before the star explodes. This merger mixes the incomplete burning material in the Si shell, including 51V and 52Cr, out to the edge of the carbon/oxygen (CO) core. The adopted supernova model ejects the outer 2 of the CO core, which includes a significant fraction of the Cr-rich material. When including this 20 M model at Z = 0.01 in the yields interpolation scheme of our GCE model for stars between 15 and 25 , we overestimate Cr/Fe by an order of magnitude at Fe/H 0 relative to observations in the Galactic disk. This raises a number of questions regarding the occurrence of Si-C shell mergers in nature, the accuracy of different simulation approaches, and the impact of such mergers on the presupernova structure and explosion dynamics. According to the conditions in this 1D stellar model, the substantial penetration of C-shell material into the Si shell could launch a convective-reactive global oscillation if a merger does take place. In any case, GCE provides stringent constraints on the outcome of this stellar evolution phase.
ABSTRACT
The abundances of neutron (n)-capture elements in the carbon-enhanced metal-poor (CEMP)-r/s stars agree with predictions of intermediate n-density nucleosynthesis, at Nn ∼ 1013–1015 cm−3, in ...rapidly accreting white dwarfs (RAWDs). We have performed Monte Carlo simulations of this intermediate-process (i-process) nucleosynthesis to determine the impact of (n,γ) reaction rate uncertainties of 164 unstable isotopes, from 131I to 189Hf, on the predicted abundances of 18 elements from Ba to W. The impact study is based on two representative one-zone models with constant values of Nn = 3.16 × 1014 and 3.16 × 1013 cm−3 and on a multizone model based on a realistic stellar evolution simulation of He-shell convection entraining H in a RAWD model with Fe/H = −2.6. For each of the selected elements, we have identified up to two (n,γ) reactions having the strongest correlations between their rate variations constrained by Hauser–Feshbach computations and the predicted abundances, with the Pearson product–moment correlation coefficients |rP| > 0.15. We find that the discrepancies between the predicted and observed abundances of Ba and Pr in the CEMP-i star CS 31062−050 are significantly diminished if the rate of 137Cs(n,γ)138Cs is reduced and the rates of 141Ba(n,γ)142Ba or 141La(n,γ)142La increased. The uncertainties of temperature-dependent β-decay rates of the same unstable isotopes have a negligible effect on the predicted abundances. One-zone Monte Carlo simulations can be used instead of computationally time-consuming multizone Monte Carlo simulations in reaction rate uncertainty studies if they use comparable values of Nn. We discuss the key challenges that RAWD simulations of i process for CEMP-i stars meet by contrasting them with recently published low-Z asymptotic giant branch (AGB) i process.