No Arabic 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 $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 red-wards. The amount of CBM impacts the size of the convective cores and the nucleosynthesis, e.g. the $^{12}$C to $^{16}$O ratio and the weak s-process. Lastly, we determine the uncertainty that the range of parameter values investigated introduce and we find differences of up to $70%$ for the core masses and the total mass of the star.
Convection plays a key role in the evolution of stars due to energy transport and mixing of composition. Despite its importance, this process is still not well understood. One longstanding conundrum in all 1D stellar evolution codes is the treatment of convective boundaries. In this study we compare two convective uncertainties, the boundary location (Ledoux versus Schwarzschild) and the amount of extra mixing, and their impact on the early evolution of massive stars. With increasing convective boundary mixing (CBM), we find a convergence of the two different boundary locations, a decreasing blue to red super giant ratio and a reduced importance of semiconvection.
We find significant fluctuations of angular momentum within the convective helium shell of a pre-collapse massive star - a core-collapse supernova progenitor - which may facilitate the formation of accretion disks and jets that can explode the star. The convective flow in our model of an evolved M_ZAMS=15Msun star, computed with the sub-sonic hydrodynamic solver MAESTRO, contains entire shells with net angular momentum in different directions. This phenomenon may have important implications for the late evolutionary stages of massive stars, and for the dynamics of core-collapse.
The convection that takes place in the innermost shells of massive stars plays an important role in the formation of core-collapse supernova explosions. Upon encountering the supernova shock, additional turbulence is generated, amplifying the explosion. In this work, we study how the convective perturbations evolve during the stellar collapse. Our main aim is to establish their physical properties right before they reach the supernova shock. To this end, we solve the linearized hydrodynamics equations perturbed on a stationary background flow. The latter is approximated by the spherical transonic Bondi accretion, while the convective perturbations are modeled as a combination of entropy and vorticity waves. We follow their evolution from large radii, where convective shells are initially located, down to small radii, where they are expected to encounter the accretion shock above the proto-neutron star. Considering typical vorticity perturbations with a Mach number $sim 0.1$ and entropy perturbations with magnitude $sim 0.05 k_mathrm{b}/mathrm{baryon}$, we find that the advection of these perturbations down to the shock generates acoustic waves with a relative amplitude $delta p/gamma p lesssim 10%$, in agreement with published numerical simulations. The velocity perturbations consist of contributions from acoustic and vorticity waves with values reaching $sim 10%$ of the sound speed ahead of the shock. The perturbation amplitudes decrease with increasing $ell$ and initial radii of the convective shells.
In the era of advanced electromagnetic and gravitational wave detectors, it has become increasingly important to effectively combine and study the impact of stellar evolution on binaries and dynamical systems of stars. Systematic studies dedicated to exploring uncertain parameters in stellar evolution are required to account for the recent observations of the stellar populations. We present a new approach to the commonly used Single-Star Evolution (SSE) fitting formulae, one that is more adaptable: Method of Interpolation for Single Star Evolution (METISSE). It makes use of interpolation between sets of pre-computed stellar tracks to approximate evolution parameters for a population of stars. We have used METISSE with detailed stellar tracks computed by the Modules for Experiments in Stellar Astrophysics (MESA), Bonn Evolutionary Code (BEC) and Cambridge STARS code. METISSE better reproduces stellar tracks computed using the STARS code compared to SSE, and is on average three times faster. Using stellar tracks computed with MESA and BEC, we apply METISSE to explore the differences in the remnant masses, the maximum radial expansion, and the main-sequence lifetime of massive stars. We find that different physical ingredients used in the evolution of stars, such as the treatment of radiation dominated envelopes, can impact their evolutionary outcome. For stars in the mass range 9 to 100 M$_odot$, the predictions of remnant masses can vary by up to 20 M$_odot$, while the maximum radial expansion achieved by a star can differ by an order of magnitude between different stellar models.
Spectroscopic studies of Galactic O and B stars show that many stars with masses above 8 M$_{odot}$ are observed in the HR diagram just beyond the Main-Sequence (MS) band predicted by stellar models computed with a moderate overshooting. This may be an indication that the convective core sizes in stars in the upper part of the HR diagram are larger than predicted by these models. Combining stellar evolution models and spectroscopic parameters derived for a large sample of Galactic O and B stars, including brand new information about their projected rotational velocities, we reexamine the question of the convective core size in MS massive stars. We confirm that for stars more massive than about 8 M$_{odot}$, the convective core size at the end of the MS phase increases more rapidly with the mass than in models computed with a constant step overshoot chosen to reproduce the main sequence width in the low mass range (around 2 M$_{odot}$). This conclusion is valid for both the cases of non-rotating models and rotating models either with a moderate or a strong angular momentum transport. The increase of the convective core mass with the mass obtained from the TAMS position is, however, larger than the one deduced from the surface velocity drop for masses above about 15 M$_{odot}$. Although observations available at the moment cannot decide what is the best choice between the core sizes given by the TAMS and the velocity drop, we discuss different methods to get out of this dilemma. At the moment, comparisons with eclipsing binaries seem to favor the solution given by the velocity drop. While we confirm the need for larger convective cores at higher masses, we find tensions in-between different methods for stars more massive than 15 M$_{odot}$. The use of single-aged stellar populations (non-interacting binaries or stellar clusters) would be a great asset to resolve this tension.