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
NSF-PAR ID:
10189185
Author(s) / Creator(s):
; ; ; ; ;
Date Published:
Journal Name:
Monthly Notices of the Royal Astronomical Society
Volume:
496
Issue:
2
ISSN:
0035-8711
Page Range / eLocation ID:
1967 to 1989
Format(s):
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
  2. ABSTRACT

    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. ABSTRACT

    Our understanding of stellar structure and evolution coming from one-dimensional (1D) stellar models is limited by uncertainties related to multidimensional processes taking place in stellar interiors. 1D models, however, can now be tested and improved with the help of detailed three-dimensional (3D) hydrodynamics models, which can reproduce complex multidimensional processes over short time-scales, thanks to the recent advances in computing resources. Among these processes, turbulent entrainment leading to mixing across convective boundaries is one of the least understood and most impactful. Here, we present the results from a set of hydrodynamics simulations of the neon-burning shell in a massive star, and interpret them in the framework of the turbulent entrainment law from geophysics. Our simulations differ from previous studies in their unprecedented degree of realism in reproducing the stellar environment. Importantly, the strong entrainment found in the simulations highlights the major flaws of the current implementation of convective boundary mixing in 1D stellar models. This study therefore calls for major revisions of how convective boundaries are modelled in 1D, and in particular the implementation of entrainment in these models. This will have important implications for supernova theory, nucleosynthesis, neutron stars, and black holes physics.

     
    more » « less