Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userThe fates of massive stars: exploring uncertainties in stellar evolution with METISSE
169
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
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.
rate research
Read More
We assess the systematic uncertainties in stellar evolutionary calculations for low- to intermediate-mass, main-sequence stars. We compare published stellar tracks from several different evolution codes with our own tracks computed using the stellar codes STARS and MESA. In particular, we focus on tracks of 1 and 3 solar masses at solar metallicity. We find that the spread in the available 1 solar mass tracks (computed before the recent solar composition revision by Asplund et al.) can be covered by tracks between 0.97-1.01 solar masses computed with the STARS code. We assess some possible causes of the origin of this uncertainty, including how the choice of input physics and the solar constraints used to perform the solar calibration affect the tracks. We find that for a 1 solar mass track, uncertainties of around 10% in the initial hydrogen abundance and initial metallicity produce around a 2% error in mass. For the 3 solar mass tracks, there is very little difference between the tracks from the various different stellar codes. The main difference comes in the extent of the main sequence, which we believe results from the different choices of the implementation of convective overshooting in the core. Uncertainties in the initial abundances lead to a 1-2% error in the mass determination. These uncertainties cover only part of the total error budget, which should also include uncertainties in the input physics (e.g., reaction rates, opacities, convective models) and any missing physics (e.g., radiative levitation, rotation, magnetic fields). Uncertainties in stellar surface properties such as luminosity and effective temperature will further reduce the accuracy of any potential mass determinations.
The formation of supermassive stars (SMSs) via rapid mass accretion and their direct collapse into black holes (BHs) is a promising pathway for sowing seeds of supermassive BHs in the early universe. We calculate the evolution of rapidly accreting SMSs by solving the stellar structure equations including nuclear burning as well as general relativistic (GR) effects up to the onset of the collapse. We find that such SMSs have less concentrated structure than fully-convective counterpart, which is often postulated for non-accreting ones. This effect stabilizes the stars against GR instability even above the classical upper mass limit $gtrsim 10^5~M_odot$ derived for the fully-convective stars. The accreting SMS begins to collapse at the higher mass with the higher accretion rate. The collapse occurs when the nuclear fuel is exhausted only for cases with $dot M lesssim 0.1~M_odot~{rm yr}^{-1}$. With $dot{M} simeq 0.3 - 1~M_odot~{rm yr}^{-1}$, the star becomes GR-unstable during the helium-burning stage at $M simeq 2 - 3.5~times 10^5~M_odot$. In an extreme case with $10~M_odot~{rm yr}^{-1}$, the star does not collapse until the mass reaches $simeq 8.0times 10^5~M_odot$, where it is still in the hydrogen-burning stage. We expect that BHs with roughly the same mass will be left behind after the collapse in all the cases.
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.
Numerous physical aspects of stellar physics have been presented in Ses- sion 2 and the underlying uncertainties have been tentatively assessed. We try here to highlight some specific points raised after the talks and during the general discus- sion at the end of the session and eventually at the end of the workshop. A table of model uncertainties is then drawn with the help of the participants in order to give the state of the art in stellar modeling uncertainties as of July 2013.
The rotational evolution of cool dwarfs is poorly constrained after around 1-2 Gyr due to a lack of precise ages and rotation periods for old main-sequence stars. In this work we use velocity dispersion as an age proxy to reveal the temperature-dependent rotational evolution of low-mass Kepler dwarfs, and demonstrate that kinematic ages could be a useful tool for calibrating gyrochronology in the future. We find that a linear gyrochronology model, calibrated to fit the period-Teff relationship of the Praesepe cluster, does not apply to stars older than around 1 Gyr. Although late-K dwarfs spin more slowly than early-K dwarfs when they are young, at old ages we find that late-K dwarfs rotate at the same rate or faster than early-K dwarfs of the same age. This result agrees qualitatively with semi-empirical models that vary the rate of surface-to-core angular momentum transport as a function of time and mass. It also aligns with recent observations of stars in the NGC 6811 cluster, which indicate that the surface rotation rates of K dwarfs go through an epoch of inhibited evolution. We find that the oldest Kepler stars with measured rotation periods are late-K and early-M dwarfs, indicating that these stars maintain spotted surfaces and stay magnetically active longer than more massive stars. Finally, based on their kinematics, we confirm that many rapidly rotating GKM dwarfs are likely to be synchronized binaries.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
