skip to main content

Title: Relative importance of convective uncertainties in massive stars
ABSTRACT In this work, we investigate the impact of uncertainties due to convective boundary mixing (CBM), commonly called ‘overshoot’, namely the boundary location and the amount of mixing at the convective boundary, on stellar structure and evolution. For this we calculated two grids of stellar evolution models with the MESA code, each with the Ledoux and the Schwarzschild boundary criterion, and vary the amount of CBM. We calculate each grid with the initial masses of 15, 20, and $25\, \rm {M}_\odot$. We present the stellar structure of the models during the hydrogen and helium burning phases. In the latter, we examine the impact on the nucleosynthesis. We find a broadening of the main sequence with more CBM, which is more in agreement with observations. Furthermore, during the core hydrogen burning phase there is a convergence of the convective boundary location due to CBM. The uncertainties of the intermediate convective zone remove this convergence. The behaviour of this convective zone strongly affects the surface evolution of the model, i.e. how fast it evolves redwards. The amount of CBM impacts the size of the convective cores and the nucleosynthesis, e.g. the 12C to 16O ratio and the weak s-process. Lastly, we determine the uncertainty that the range of parameter values investigated introduces and we find differences of up to $70{{\ \rm per\ cent}}$ for the core masses and the total mass of the star.  more » « less
Award ID(s):
1927130 1430152
Author(s) / Creator(s):
; ; ; ; ;
Date Published:
Journal Name:
Monthly Notices of the Royal Astronomical Society
Page Range / eLocation ID:
1967 to 1989
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. ABSTRACT We have modelled the multicycle evolution of rapidly accreting CO white dwarfs (RAWDs) with stable H burning intermittent with strong He-shell flashes on their surfaces for 0.7 ≤ MRAWD/M⊙ ≤ 0.75 and [Fe/H] ranging from 0 to −2.6. We have also computed the i-process nucleosynthesis yields for these models. The i process occurs when convection driven by the He-shell flash ingests protons from the accreted H-rich surface layer, which results in maximum neutron densities Nn, max ≈ 1013–1015 cm−3. The H-ingestion rate and the convective boundary mixing (CBM) parameter ftop adopted in the one-dimensional nucleosynthesis and stellar evolution models are constrained through three-dimensional (3D) hydrodynamic simulations. The mass ingestion rate and, for the first time, the scaling laws for the CBM parameter ftop have been determined from 3D hydrodynamic simulations. We confirm our previous result that the high-metallicity RAWDs have a low mass retention efficiency ($\eta \lesssim 10{{\ \rm per\ cent}}$). A new result is that RAWDs with [Fe/H] $\lesssim -2$ have $\eta \gtrsim 20{{\ \rm per\ cent}}$; therefore, their masses may reach the Chandrasekhar limit and they may eventually explode as SNeIa. This result and the good fits of the i-process yields from the metal-poor RAWDs to the observed chemical composition of the CEMP-r/s stars suggest that some of the present-day CEMP-r/s stars could be former distant members of triple systems, orbiting close binary systems with RAWDs that may have later exploded as SNeIa. 
    more » « less

    Our knowledge of stellar evolution is driven by one-dimensional (1D) simulations. 1D models, however, are severely limited by uncertainties on the exact behaviour of many multidimensional phenomena occurring inside stars, affecting their structure and evolution. Recent advances in computing resources have allowed small sections of a star to be reproduced with multi-D hydrodynamic models, with an unprecedented degree of detail and realism. In this work, we present a set of 3D simulations of a convective neon-burning shell in a 20 M⊙ star run for the first time continuously from its early development through to complete fuel exhaustion, using unaltered input conditions from a 321D-guided 1D stellar model. These simulations help answer some open questions in stellar physics. In particular, they show that convective regions do not grow indefinitely due to entrainment of fresh material, but fuel consumption prevails over entrainment, so when fuel is exhausted convection also starts decaying. Our results show convergence between the multi-D simulations and the new 321D-guided 1D model, concerning the amount of convective boundary mixing to include in stellar models. The size of the convective zones in a star strongly affects its structure and evolution; thus, revising their modelling in 1D will have important implications for the life and fate of stars. This will thus affect theoretical predictions related to nucleosynthesis, supernova explosions, and compact remnants.

    more » « less
  3. null (Ed.)
    The slow neutron-capture process (s-process) efficiency in low-mass AGB stars (1.5 < M/M⊙ < 3) critically depends on how mixing processes in stellar interiors are handled, which is still affected by considerable uncertainties. In this work, we compute the evolution and nucleosynthesis of low-mass AGB stars at low metallicities using the MESA stellar evolution code. The combined data set includes models with initial masses Mini/M⊙=2 and 3 for initial metallicities Z=0.001 and 0.002. The nucleosynthesis was calculated for all relevant isotopes by post-processing with the NuGrid mppnp code. Using these models, we show the impact of the uncertainties affecting the main mixing processes on heavy element nucleosynthesis, such as convection and mixing at convective boundaries. We finally compare our theoretical predictions with observed surface abundances on low-metallicity stars. We find that mixing at the interface between the He-intershell and the CO-core has a critical impact on the s-process at low metallicities, and its importance is comparable to convective boundary mixing processes under the convective envelope, which determine the formation and size of the 13C-pocket. Additionally, our results indicate that models with very low to no mixing below the He-intershell during thermal pulses, and with a 13C-pocket size of at least ∼3 × 10−4 M⊙, are strongly favored in reproducing observations. Online access to complete yield data tables is also provided. 
    more » « less
  4. Abstract

    Observations indicate that the convective cores of stars must ingest a substantial amount of material from the overlying radiative zone, but the extent of this mixing and the mechanism that causes it remain uncertain. Recently, Anders et al. developed a theory of convective penetration and calibrated it with 3D numerical hydrodynamics simulations. Here we employ that theory to predict the extent of convective boundary mixing (CBM) in early-type main-sequence stars. We find that convective penetration produces enough mixing to explain core masses inferred from asteroseismology and eclipsing binary studies, and matches observed trends in mass and age. While there are remaining uncertainties in the theory, this agreement suggests that most CBM in early-type main-sequence stars arises from convective penetration. Finally, we provide a fitting formula for the extent of core convective penetration for main-sequence stars in the mass range from 1.1–60M.

    more » « less
  5. Context: Modeling of the stars in the red clump (RC), that is, core helium-burning stars that have gone through a He flash, is challenging because of the uncertainties associated with the physical processes in their core and during the helium flash. By probing the internal stellar structure, asteroseismology allows us to constrain the core properties of RC stars and eventually, to improve our understanding of this evolutionary phase. Aims: We aim to quantify the impact on the seismic properties of the RC stars of the two main core modeling uncertainties: core boundary mixing, and helium-burning nuclear reaction rates. Methods: Using the MESA stellar evolution code, we computed models with different core boundary mixing as well as different 3α and 12C(alpha, gamma)16O nuclear reaction rates. We investigated the impact of these parameters on the period spacing ΔΠ, which is a probe of the region around the core. Results: We find that different core boundary mixing schemes yield significantly different period spacings, with differences of 30 s between the maximum ΔΠ value computed with semiconvection and maximal overshoot. We show that an increased rate of 12C(alpha, gamma)16O lengthens the core helium-burning phase, which extends the range of ΔΠ covered by the models during their evolution. This results in a difference of 10 s between the models computed with a nominal rate and a rate multiplied by 2, which exceeds the observational uncertainties. The effect of changing the 3α reaction rate is comparatively small. Conclusions. The core boundary mixing is the main source of uncertainty in the seismic modeling of RC stars. Moreover, the effect of the 12C(alpha, gamma)16O is non-negligible, even though it is difficult to distinguish from the effect of the mixing. This degeneracy could be seen more frequently in the future in the new seismic data from the PLATO mission and through theoretical constraints from numerical simulations. 
    more » « less