More Like this

1. Evolution of Wolf–Rayet stars as black hole progenitors

   https://doi.org/10.1093/mnras/stab1548  

   Higgins, E R ; Sander, A A ; Vink, J S ; Hirschi, R ( June 2021 , Monthly Notices of the Royal Astronomical Society)

   null (Ed.)

   ABSTRACT Evolved Wolf–Rayet stars form a key aspect of massive star evolution, and their strong outflows determine their final fates. In this study, we calculate grids of stellar models for a wide range of initial masses at five metallicities (ranging from solar down to just 2 per cent solar). We compare a recent hydrodynamically consistent wind prescription with two earlier frequently used wind recipes in stellar evolution and population synthesis modelling, and we present the ranges of maximum final masses at core He-exhaustion for each wind prescription and metallicity Z. Our model grids reveal qualitative differences in mass-loss behaviour of the wind prescriptions in terms of ‘convergence’. Using the prescription from Nugis & Lamers the maximum stellar black hole is found to converge to a value of 20–30 M⊙, independent of host metallicity; however, when utilizing the new physically motivated prescription from Sander & Vink there is no convergence to a maximum black hole mass value. The final mass is simply larger for larger initial He-star mass, which implies that the upper black hole limit for He-stars below the pair-instability gap is set by prior evolution with mass loss, or the pair instability itself. Quantitatively, we find the critical Z for pair-instability (ZPI) to be as high as 50 per cent Z⊙, corresponding to the host metallicity of the Large Magellanic Cloud. Moreover, while the Nugis & Lamers prescription would not predict any black holes above the approx 130 M⊙ pair-instability limit, with Sander & Vink winds included, we demonstrate a potential channel for very massive helium stars to form such massive black holes at ∼2 per cent Z⊙ or below. 

   more » « less
2. Is GW190521 the merger of black holes from the first stellar generations?

   Farrell, E.J. ; Groh, J.H. ; Hirschi, R. ; Murphy, L. ; Kaiser, E. ; Ekström, S. ; Georgy, C. ; Meynet, G. ( September 2020 , ArXivorg)

   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 which favour higher BH masses as compared to higher metallicity models. These are (i) lower mass-loss rates during the post-MS 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 Msun. Uncertainties related to convective mixing, mass loss, H-He shell interactions and pair-instability pulsations may increase this limit to ~85 Msun. 

   more » « less
3. Grids of stellar models with rotation VII: models from 0.8 to 300 M⊙ at supersolar metallicity ( Z  = 0.020)

   https://doi.org/10.1093/mnras/stac230  

   Yusof, Norhasliza ; Hirschi, Raphael ; Eggenberger, Patrick ; Ekström, Sylvia ; Georgy, Cyril ; Sibony, Yves ; Crowther, Paul A. ; Meynet, Georges ; Kassim, Hasan Abu ; Harun, Wan Aishah Wan ; et al ( January 2022 , Monthly Notices of the Royal Astronomical Society)

   We present a grid of stellar models at supersolar metallicity (Z = 0.020) extending the previous grids of Geneva models at solar and sub-solar metallicities. A metallicity of Z = 0.020 was chosen to match that of the inner Galactic disc. A modest increase of 43 per cent (= 0.02/0.014) in metallicity compared to solar models means that the models evolve similarly to solar models but with slightly larger mass-loss. Mass-loss limits the final total masses of the supersolar models to 35 M⊙ even for stars with initial masses much larger than 100 M⊙. Mass-loss is strong enough in stars above 20 M⊙ for rotating stars (25 M⊙ for non-rotating stars) to remove the entire hydrogen-rich envelope. Our models thus predict SNII below 20 M⊙ for rotating stars (25 M⊙ for non-rotating stars) and SNIb (possibly SNIc) above that. We computed both isochrones and synthetic clusters to compare our supersolar models to the Westerlund 1 (Wd1) massive young cluster. A synthetic cluster combining rotating and non-rotating models with an age spread between log10(age/yr) = 6.7 and 7.0 is able to reproduce qualitatively the observed populations of WR, RSG, and YSG stars in Wd1, in particular their simultaneous presence at $\log _{10}(L/\mathit {\mathrm{ L}}_{\odot })$ = 5–5.5. The quantitative agreement is imperfect and we discuss the likely causes: synthetic cluster parameters, binary interactions, mass-loss and their related uncertainties. In particular, mass-loss in the cool part of the HRD plays a key role.

    

   more » « less
4. The Aarhus red giants challenge: I. Stellar structures in the red giant branch phase

   https://doi.org/10.1051/0004-6361/201935843  

   Silva Aguirre, V. ; Christensen-Dalsgaard, J. ; Cassisi, S. ; Miller Bertolami, M. ; Serenelli, A. ; Stello, D. ; Weiss, A. ; Angelou, G. ; Jiang, C. ; Lebreton, Y. ; et al ( March 2020 , Astronomy & Astrophysics)

   Context. With the advent of space-based asteroseismology, determining accurate properties of red-giant stars using their observed oscillations has become the focus of many investigations due to their implications in a variety of fields in astrophysics. Stellar models are fundamental in predicting quantities such as stellar age, and their reliability critically depends on the numerical implementation of the physics at play in this evolutionary phase. Aims. We introduce the Aarhus red giants challenge, a series of detailed comparisons between widely used stellar evolution and oscillation codes that aim to establish the minimum level of uncertainties in properties of red giants arising solely from numerical implementations. We present the first set of results focusing on stellar evolution tracks and structures in the red-giant-branch (RGB) phase. Methods. Using nine state-of-the-art stellar evolution codes, we defined a set of input physics and physical constants for our calculations and calibrated the convective efficiency to a specific point on the main sequence. We produced evolutionary tracks and stellar structure models at a fixed radius along the red-giant branch for masses of 1.0  M ⊙ , 1.5  M ⊙ , 2.0  M ⊙ , and 2.5  M ⊙ , and compared the predicted stellar properties. Results. Once models have been calibrated on the main sequence, we find a residual spread in the predicted effective temperatures across all codes of ∼20 K at solar radius and ∼30–40 K in the RGB regardless of the considered stellar mass. The predicted ages show variations of 2–5% (increasing with stellar mass), which we attribute to differences in the numerical implementation of energy generation. The luminosity of the RGB-bump shows a spread of about 10% for the considered codes, which translates into magnitude differences of ∼0.1 mag in the optical V -band. We also compare the predicted [C/N] abundance ratio and find a spread of 0.1 dex or more for all considered masses. Conclusions. Our comparisons show that differences at the level of a few percent still remain in evolutionary calculations of red giants branch stars despite the use of the same input physics. These are mostly due to differences in the energy generation routines and interpolation across opacities, and they call for further investigation on these matters in the context of using properties of red giants as benchmarks for astrophysical studies. 

   more » « less
5. Evolutionary roads leading to low effective spins, high black hole masses, and O1/O2 rates for LIGO/Virgo binary black holes

   https://doi.org/10.1051/0004-6361/201936528  

   Belczynski, K. ; Klencki, J. ; Fields, C. E. ; Olejak, A. ; Berti, E. ; Meynet, G. ; Fryer, C. L. ; Holz, D. E. ; O’Shaughnessy, R. ; Brown, D. A. ; et al ( April 2020 , Astronomy & Astrophysics)

   All ten LIGO/Virgo binary black hole (BH-BH) coalescences reported following the O1/O2 runs have near-zero effective spins. There are only three potential explanations for this. If the BH spin magnitudes are large, then: (i) either both BH spin vectors must be nearly in the orbital plane or (ii) the spin angular momenta of the BHs must be oppositely directed and similar in magnitude. Then there is also the possibility that (iii) the BH spin magnitudes are small. We consider the third hypothesis within the framework of the classical isolated binary evolution scenario of the BH-BH merger formation. We test three models of angular momentum transport in massive stars: a mildly efficient transport by meridional currents (as employed in the Geneva code), an efficient transport by the Tayler-Spruit magnetic dynamo (as implemented in the MESA code), and a very-efficient transport (as proposed by Fuller et al.) to calculate natal BH spins. We allow for binary evolution to increase the BH spins through accretion and account for the potential spin-up of stars through tidal interactions. Additionally, we update the calculations of the stellar-origin BH masses, including revisions to the history of star formation and to the chemical evolution across cosmic time. We find that we can simultaneously match the observed BH-BH merger rate density and BH masses and BH-BH effective spins. Models with efficient angular momentum transport are favored. The updated stellar-mass weighted gas-phase metallicity evolution now used in our models appears to be key for obtaining an improved reproduction of the LIGO/Virgo merger rate estimate. Mass losses during the pair-instability pulsation supernova phase are likely to be overestimated if the merger GW170729 hosts a BH more massive than 50  M ⊙ . We also estimate rates of black hole-neutron star (BH-NS) mergers from recent LIGO/Virgo observations. If, in fact. angular momentum transport in massive stars is efficient, then any (electromagnetic or gravitational wave) observation of a rapidly spinning BH would indicate either a very effective tidal spin up of the progenitor star (homogeneous evolution, high-mass X-ray binary formation through case A mass transfer, or a spin- up of a Wolf-Rayet star in a close binary by a close companion), significant mass accretion by the hole, or a BH formation through the merger of two or more BHs (in a dense stellar cluster). 

   more » « less