ABSTRACT
GW190521 challenges our understanding of the late-stage evolution of massive stars and the effects of the pair instability in particular. We discuss the possibility that stars at low or zero ...metallicity could retain most of their hydrogen envelope until the pre-supernova stage, avoid the pulsational pair-instability regime, and produce a black hole with a mass in the mass gap by fallback. We present a series of new stellar evolution models at zero and low metallicity computed with the geneva and mesa stellar evolution codes and compare to existing grids of models. Models with a metallicity in the range 0–0.0004 have three properties that favour higher black hole (BH) masses. These are (i) lower mass-loss rates during the post main sequence phase, (ii) a more compact star disfavouring binary interaction, and (iii) possible H–He shell interactions which lower the CO core mass. We conclude that it is possible that GW190521 may be the merger of black holes produced directly by massive stars from the first stellar generations. Our models indicate BH masses up to 70–75 M⊙. Uncertainties related to convective mixing, mass loss, H–He shell interactions, and pair-instability pulsations may increase this limit to ∼85 M⊙.
ABSTRACT
We show that it is not possible to determine the final mass Mfin of a red supergiant (RSG) at the pre-supernova (SN) stage from its luminosity L and effective temperature Teff alone. Using a ...grid of stellar models, we demonstrate that for a given value of L and Teff, an RSG can have a range of Mfin as wide as 3 to 45 M⊙. While the probability distribution within these limits is not flat, any individual determination of Mfin for an RSG will be degenerate. This makes it difficult to determine its evolutionary history and to map Mfin to an initial mass. Single stars produce a narrower range that is difficult to accurately determine without making strong assumptions about mass-loss, convection, and rotation. Binaries would produce a wider range of RSG Mfin. However, the final Helium core mass $M_{\operatorname{He-core}}$ is well determined by the final luminosity and we find $\log (M_{\operatorname{He-core}}/\mathrm{M}_{\odot }) = 0.659 \log (L/\mathrm{L}_{\odot }) -2.630$. Using this relationship, we derive $M_{\operatorname{He-core}}$ for directly imaged SN progenitors and one failed SN candidate. The value of Mfin for stripped star progenitors of SNe IIb is better constrained by L and Teff due to the dependence of Teff on the envelope mass Menv for Menv ≲ 1 M⊙. Given the initial mass function, our results apply to the majority of progenitors of core-collapse SNe, failed SNe, and direct-collapse black holes.
Abstract
We update the capabilities of the open-knowledge software instrument Modules for Experiments in Stellar Astrophysics (
MESA
). The new
auto
_
diff
module implements automatic differentiation ...in
MESA
, an enabling capability that alleviates the need for hard-coded analytic expressions or finite-difference approximations. We significantly enhance the treatment of the growth and decay of convection in
MESA
with a new model for time-dependent convection, which is particularly important during late-stage nuclear burning in massive stars and electron-degenerate ignition events. We strengthen
MESA
’s implementation of the equation of state, and we quantify continued improvements to energy accounting and solver accuracy through a discussion of different energy equation features and enhancements. To improve the modeling of stars in
MESA
, we describe key updates to the treatment of stellar atmospheres, molecular opacities, Compton opacities, conductive opacities, element diffusion coefficients, and nuclear reaction rates. We introduce treatments of starspots, an important consideration for low-mass stars, and modifications for superadiabatic convection in radiation-dominated regions. We describe new approaches for increasing the efficiency of calculating monochromatic opacities and radiative levitation, and for increasing the efficiency of evolving the late stages of massive stars with a new operator-split nuclear burning mode. We close by discussing major updates to
MESA
’s software infrastructure that enhance source code development and community engagement.
ABSTRACT
We investigate the low-luminosity supernova SN 2016bkv and its peculiar early-time interaction. For that, we compute radiative transfer models using the cmfgen code. Because SN 2016bkv shows ...signs of interaction with material expelled by its progenitor, it offers a great opportunity to constrain the uncertain evolutionary channels leading to low-luminosity supernovae. Our models indicate that the progenitor had a mass-loss rate of (6.0 ± 2.0) × 10−4$\mathit {\rm M}_{\odot }\,{\rm yr}^{-1}$(assuming a velocity of 150 ${\rm km\, s^{-1}}$). The surface abundances of the progenitor are consistent with solar contents of He and CNO. If SN 2016bkv’s progenitor evolved as a single star, it was an odd red supergiant that did not undergo the expected dredge up for some reason. We propose that the progenitor more likely evolved through binary interaction. One possibility is that the primary star accreted unprocessed material from a companion and avoided further rotational and convective mixing until the SN explosion. Another possibility is a merger with a lower mass star, with the primary remaining with low N abundance until core collapse. Given the available merger models, we can only put a loose constraint on the pre-explosion mass around 10–20 $\mathit {\rm M}_{\odot }$, with lower values being favoured based on previous observational constraints from the nebular phase.
We investigate the effects of mass loss during the main-sequence (MS) and post-MS phases of massive star evolution on black hole (BH) birth masses. We compute solar metallicity Geneva stellar ...evolution models of an 85 star with mass-loss rate ( ) prescriptions for MS and post-MS phases and analyze under which conditions such models could lead to very massive BHs. Based on the observational constraints for of luminous stars, we discuss two possible scenarios that could produce massive BHs at high metallicity. First, if a massive BH progenitor evolves from the observed population of massive MS stars known as WNh stars, we show that its average post-MS mass-loss rate has to be less than . However, this is lower than the typical observed mass-loss rates of luminous blue variables (LBV). Second, a massive BH progenitor could evolve from a yet undetected population of 80-85 stars with strong surface magnetic fields, which could quench mass loss during the evolution. In this case, the average mass-loss rate during the post-MS LBV phase has to be less than 5 × 10−5 to produce 70 BHs. We suggest that LBVs that explode as SNe have large envelopes and small cores that could be prone to explosion, possibly evolving from binary interaction (either mergers or mass gainers that do not fully mix). Conversely, LBVs that directly collapse to BHs could have evolved from massive single stars or binary-star mergers that fully mix, possessing large cores that would favor BH formation.
ABSTRACT
We investigate a suspected very massive star in one of the most metal-poor dwarf galaxies, PHL 293B. Excitingly, we find the sudden disappearance of the stellar signatures from our 2019 ...spectra, in particular the broad H lines with P Cygni profiles that have been associated with a massive luminous blue variable (LBV) star. Such features are absent from our spectra obtained in 2019 with the Echelle Spectrograph for Rocky Exoplanet- and Stable Spectroscopic Observation and X-shooter instruments of the European Southern Observatory’s Very Large Telescope. We compute radiative transfer models using cmfgen, which fit the observed spectrum of the LBV and are consistent with ground-based and archival Hubble Space Telescope photometry. Our models show that during 2001–2011, the LBV had a luminosity L* = 2.5–3.5 × 106 L⊙, a mass-loss rate $\dot{M} = 0.005{-}0.020 ~{\rm M}_{\odot }$ yr−1, a wind velocity of 1000 km s−1, and effective and stellar temperatures of Teff = 6000–6800 and T* = 9500–15 000 K. These stellar properties indicate an eruptive state. We consider two main hypotheses for the absence of the broad emission components from the spectra obtained since 2011. One possibility is that we are seeing the end of an LBV eruption of a surviving star, with a mild drop in luminosity, a shift to hotter effective temperatures, and some dust obscuration. Alternatively, the LBV could have collapsed to a massive black hole without the production of a bright supernova.
ABSTRACT
We introduce snapshot, a technique to systematically compute stellar structure models in hydrostatic and thermal equilibrium based on three structural properties – core mass Mcore, envelope ...mass Menv, and core composition. This approach allows us to connect these properties of stellar interiors to the luminosity and effective temperature Teff in a more systematic way than with stellar evolution models. We compute core-H burning models with total masses of Mtotal = 8–60 M⊙ and central H mass fractions from 0.70 to 0.05. Using these, we derive an analytical relationship between Mcore, Mtotal, and central H abundance that can be readily used in rapid stellar evolution algorithms. In contrast, core-He burning stars can have a wide range of combinations of Mcore, Menv, and core compositions. We compute core-He burning models with Mcore = 2–9 M⊙, Menv = 0–50 M⊙, and central He mass fractions of 0.50 and 0.01. Models with Mcore/Mtotal from 0.2 to 0.8 have convective envelopes, low Teff and will appear as red supergiants (RSGs). For a given Mcore, they exhibit a small variation in luminosity (0.02 dex) and Teff ($\sim 400\, \mathrm{K}$) over a wide range of Menv ($\sim$2–20 M⊙). This means that it is not possible to derive RSG masses from luminosities and Teff alone. We derive the following relationship between Mcore and the total luminosity of an RSG during core He burning: log Mcore ≃ 0.44log L/L⊙ − 1.38. At Mcore/Mtotal ≈ 0.2, our models exhibit a bistability and jump from an RSG to a BSG structure. Our models with Mcore/Mtotal > 0.8, which correspond to stripped stars produced by mass-loss or binary interaction, show that Teff has a strong dependence on Menv, Mcore, and the core composition. We constrain the mass of one of these stripped stars in a binary system, HD 45166, and find it to be less than its estimated dynamical mass. When a large observational sample of stripped stars becomes available, our results can be used to constrain their Mcore, Menv, mass-loss rates, and the physics of binary interaction.
ABSTRACT
The evolution of massive stars is affected by a variety of physical processes, including convection, rotation, mass-loss, and binary interaction. Because these processes modify the internal ...chemical abundance profiles in multiple ways simultaneously, it can be challenging to determine which properties of the stellar interior are primarily driving the overall evolution. Building on previous work, we develop a new modelling approach called snapshot that allows us to isolate the key features of the internal abundance profile that drive the evolution of massive stars. Using our approach, we compute numerical stellar structure models in thermal equilibrium covering key phases of stellar evolution. For the main sequence, we demonstrate that models with the same mass and very similar surface properties can have different internal distributions of hydrogen and convective core masses. We discuss why massive stars expand after the main sequence and the fundamental reasons for why they become red, blue or yellow supergiants. For the post-main sequence, we demonstrate that small changes in the abundance profile can cause very large effects on the surface properties. We also discuss the effects that produce blue supergiants and the cause of blue loops. Our models show that massive stars with lower metallicity tend to be more compact due to the combined effect of lower CNO abundances in the burning regions and lower opacity in the envelope.
ABSTRACT
The long-term evolution of hydrogen-dominated atmospheres of sub-Neptune-like planets is mostly controlled to by two factors: a slow dissipation of the gravitational energy acquired at the ...formation (known as thermal evolution) and atmospheric mass-loss. Here, we use mesa to self-consistently couple the thermal evolution model of lower atmospheres with a realistic hydrodynamical atmospheric evaporation prescription. To outline the main features of such coupling, we simulate planets with a range of core masses (5–20 M⊕) and initial atmospheric mass fractions (0.5–30 per cent), orbiting a solar-like star at 0.1 au. In addition to our computed evolutionary tracks, we also study the stability of planetary atmospheres, showing that the atmospheres of light planets can be completely removed within 1 Gyr and that compact atmospheres have a better survival rate. From a detailed comparison between our results and the output of the previous-generation models, we show that coupling between thermal evolution and atmospheric evaporation considerably affects the thermal state of atmospheres for low-mass planets and, consequently, changes the relationship between atmospheric mass fraction and planetary parameters. We, therefore, conclude that self-consistent consideration of the thermal evolution and atmospheric evaporation is of crucial importance for evolutionary modelling and a better characterization of planetary atmospheres. From our simulations, we derive an analytical expression between planetary radius and atmospheric mass fraction at different ages. In particular, we find that, for a given observed planetary radius, the predicted atmospheric mass fraction changes as age0.11.
Efforts to understand the origin and growth of massive black holes observed in the early Universe have spurred significant interest in the evolution and fate of rapidly accreting primordial ...(metal-free) stars. Here, we investigate the evolution of such Population III (PopIII) stars under variable accretion rates, focusing on the thermal response and stellar structure, the impact of the luminosity wave encountered early in the pre-main sequence (pre-MS) phase, and the influence of accretion on their subsequent evolution. We employ the Geneva stellar evolution code and simulate ten models with varying accretion histories, covering a final mass range from 491
M
⊙
to 6127
M
⊙
. Our findings indicate that the critical accretion rate delineating the red and blue supergiant regimes during the pre-MS evolution is approximately 2.5 × 10
−2
M
⊙
yr
−1
. Once core hydrogen burning commences, the value of this critical accretion rate drops to 7.0 × 10
−3
M
⊙
yr
−1
. Moreover, we also confirm that the Kelvin–Helmholtz timescale in the outer surface layers is the more relevant timescale for determining the transition between red and blue phases. Regarding the luminosity wave, we find that it affects only the early pre-MS phase of evolution and does not directly influence the transition between red and blue phases, which primarily depends on the accretion rate. Finally, we demonstrate that variable accretion rates significantly impact the lifetimes, surface enrichment, and final mass of the PopIII stars, as well as the time they spend in the red phase. Our study provides a comprehensive understanding of the intricate evolutionary patterns of PopIII stars subjected to variable accretion rates.
Full text
Available for:
FMFMET, NUK, UL, UM, UPUK