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 more » 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. « less
; ; ; ; ;
Award ID(s):
1927130 1430152
Publication Date:
Journal Name:
Monthly Notices of the Royal Astronomical Society
Page Range or eLocation-ID:
1967 to 1989
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 chemicalmore »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.« less
  2. 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 leastmore »∼3 × 10−4 M⊙, are strongly favored in reproducing observations. Online access to complete yield data tables is also provided.« less
  3. 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.


    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.

  5. Context. A realistic parametrization of convection and convective boundary mixing in conventional stellar evolution codes is still the subject of ongoing research. To improve the current situation, multidimensional hydrodynamic simulations are used to study convection in stellar interiors. Such simulations are numerically challenging, especially for flows at low Mach numbers which are typical for convection during early evolutionary stages. Aims. We explore the benefits of using a low-Mach hydrodynamic flux solver and demonstrate its usability for simulations in the astrophysical context. Simulations of convection for a realistic stellar profile are analyzed regarding the properties of convective boundary mixing. Methods. The time-implicit Seven-League Hydro (SLH) code was used to perform multidimensional simulations of convective helium shell burning based on a 25  M ⊙ star model. The results obtained with the low-Mach AUSM + -up solver were compared to results when using its non low-Mach variant AUSM B + -up. We applied well-balancing of the gravitational source term to maintain the initial hydrostatic background stratification. The computational grids have resolutions ranging from 180 × 90 2 to 810 × 540 2 cells and the nuclear energy release was boosted by factors of 3 × 10 3 , 1 × 10 4 , and 3 × 10 4 tomore »study the dependence of the results on these parameters. Results. The boosted energy input results in convection at Mach numbers in the range of 10 −3 –10 −2 . Standard mixing-length theory predicts convective velocities of about 1.6 × 10 −4 if no boosting is applied. The simulations with AUSM + -up show a Kolmogorov-like inertial range in the kinetic energy spectrum that extends further toward smaller scales compared with its non low-Mach variant. The kinetic energy dissipation of the AUSM + -up solver already converges at a lower resolution compared to AUSM B + -up. The extracted entrainment rates at the boundaries of the convection zone are well represented by the bulk Richardson entrainment law and the corresponding fitting parameters are in agreement with published results for carbon shell burning. However, our study needs to be validated by simulations at higher resolution. Further, we find that a general increase in the entropy in the convection zone may significantly contribute to the measured entrainment of the top boundary. Conclusion. This study demonstrates the successful application of the AUSM + -up solver to a realistic astrophysical setup. Compressible simulations of convection in early phases at nominal stellar luminosity will benefit from its low-Mach capabilities. Similar to other studies, our extrapolated entrainment rate for the helium-burning shell would lead to an unrealistic growth of the convection zone if it is applied over the lifetime of the zone. Studies at nominal stellar luminosities and different phases of the same convection zone are needed to detect a possible evolution of the entrainment rate and the impact of radiation on convective boundary mixing.« less