skip to main content

Title: Schwarzschild and Ledoux are Equivalent on Evolutionary Timescales

Stellar evolution models calculate convective boundaries using either the Schwarzschild or Ledoux criterion, but confusion remains regarding which criterion to use. Here we present a 3D hydrodynamical simulation of a convection zone and adjacent radiative zone, including both thermal and compositional buoyancy forces. As expected, regions that are unstable according to the Ledoux criterion are convective. Initially, the radiative zone adjacent to the convection zone is Schwarzschild unstable but Ledoux stable due to a composition gradient. Over many convective overturn timescales, the convection zone grows via entrainment. The convection zone saturates at the size originally predicted by the Schwarzschild criterion, although in this final state the Schwarzschild and Ledoux criteria agree. Therefore, the Schwarzschild criterion should be used to determine the size of stellar convection zones, except possibly during short-lived evolutionary stages in which entrainment persists.

; ; ; ; ; ;
Award ID(s):
Publication Date:
Journal Name:
The Astrophysical Journal Letters
Page Range or eLocation-ID:
Article No. L10
DOI PREFIX: 10.3847
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    We study the properties of oscillatory double-diffusive convection (ODDC) in the presence of a uniform vertical background magnetic field. ODDC takes place in stellar regions that are unstable according to the Schwarzschild criterion and stable according to the Ledoux criterion (sometimes called semiconvective regions), which are often predicted to reside just outside the core of intermediate-mass main-sequence stars. Previous hydrodynamic studies of ODDC have shown that the basic instability saturates into a state of weak wave-like convection, but that a secondary instability can sometimes transform it into a state of layered convection, where layers then rapidly merge and grow until the entire region is fully convective. We find that magnetized ODDC has very similar properties overall, with some important quantitative differences. A linear stability analysis reveals that the fastest-growing modes are unaffected by the field, but that other modes are. Numerically, the magnetic field is seen to influence the saturation of the basic instability, overall reducing the turbulent fluxes of temperature and composition. This in turn affects layer formation, usually delaying it, and occasionally suppressing it entirely for sufficiently strong fields. Further work will be needed, however, to determine the field strength above which layer formation is actually suppressedmore »in stars. Potential observational implications are briefly discussed.

    « less
  2. 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 determinemore »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
  3. 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
  4. Abstract

    Subsurface convection zones are ubiquitous in early-type stars. Driven by narrow opacity peaks, these thin convective regions transport little heat but play an important role in setting the magnetic properties and surface variability of stars. Here we demonstrate that these convection zones arenotpresent in as wide a range of stars as previously believed. In particular, there are regions which 1D stellar evolution models report to be convectively unstable but which fall below the critical Rayleigh number for onset of convection. For sub-solar metallicity this opens up astability windowin which there are no subsurface convection zones. For Large Magellanic Cloud metallicity this surface stability region extends roughly between 8 and 16M, increasing to 8–35Mfor Small Magellanic Cloud metallicity. Such windows are then an excellent target for probing the relative influence of subsurface convection and other sources of photometric variability in massive stars.

  5. Our ability to predict the structure and evolution of stars is in part limited by complex, 3D hydrodynamic processes such as convective boundary mixing. Hydrodynamic simulations help us understand the dynamics of stellar convection and convective boundaries. However, the codes used to compute such simulations are usually tested on extremely simple problems and the reliability and reproducibility of their predictions for turbulent flows is unclear. We define a test problem involving turbulent convection in a plane-parallel box, which leads to mass entrainment from, and internal-wave generation in, a stably stratified layer. We compare the outputs from the codes FLASH , MUSIC , PPMSTAR , PROMPI , and SLH , which have been widely employed to study hydrodynamic problems in stellar interiors. The convection is dominated by the largest scales that fit into the simulation box. All time-averaged profiles of velocity components, fluctuation amplitudes, and fluxes of enthalpy and kinetic energy are within ≲3 σ of the mean of all simulations on a given grid (128 3 and 256 3 grid cells), where σ describes the statistical variation due to the flow’s time dependence. They also agree well with a 512 3 reference run. The 128 3 and 256 3 simulationsmore »agree within 9% and 4%, respectively, on the total mass entrained into the convective layer. The entrainment rate appears to be set by the amount of energy that can be converted to work in our setup and details of the small-scale flows in the boundary layer seem to be largely irrelevant. Our results lend credence to hydrodynamic simulations of flows in stellar interiors. We provide in electronic form all outputs of our simulations as well as all information needed to reproduce or extend our study.« less