skip to main content


Title: Is GW190521 the merger of black holes from the first stellar generations?
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
Award ID(s):
1927130
NSF-PAR ID:
10193320
Author(s) / Creator(s):
; ; ; ; ; ; ;
Date Published:
Journal Name:
ArXivorg
ISSN:
2331-8422
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. null (Ed.)
    Our understanding of the formation and evolution of binary black holes (BBHs) is significantly impacted by the recent discoveries made by the LIGO/Virgo collaboration. Of utmost importance is the detection of the most massive BBH system, GW190521. Here we investigate what it takes for field massive stellar binaries to account for the formation of such massive BBHs. Whether the high mass end of the BH mass function is populated by remnants of massive stars that either formed at extremely low metallicities and avoid the pair-instability mass gap or increase their birth mass beyond the pair-instability mass gap through the accretion of gas from the surrounding medium. We show that assuming that massive stars at very low metallicities can form massive BHs by avoiding pair-instability supernova, coupled with a correspondingly high formation efficiency for BBHs, can explain the observed BH mass function. To this end, one requires a relation between the initial and final mass of the progenitor stars at low metallicities that is shallower than what is expected from wind mass loss alone. On the other hand, assuming pair-instability operates at all metallicities, one can account for the observed BH mass function if at least about 10% of the BHs born at very low metallicities double their mass before they merge because of accretion of ambient gas. Such BBHs will have to spend about a Gyr within a parsec length-scale of their parent atomic cooling halos or a shorter timescale if they reside in the inner sub-parsecs of their host dark matter halos. Future stellar evolution calculations of massive stars at very low metallicity and hydrodynamical simulations of gas accretion onto BBHs born in atomic cooling halos can shed light on this debate. 
    more » « less
  2. 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
  3. Abstract The recent detection of GW190521 stimulated ideas on how to populate the predicted black hole (BH) pair-instability (PI) mass gap. One proposal is the dynamical merger of two stars below the PI regime forming a star with a small core and an oversized envelope. We outline the main challenges this scenario faces to form one BH in the gap. In particular, the core needs to avoid growing during the merger, and the merger product needs to retain enough mass, including in the subsequent evolution, and at core collapse (CC). We explore this scenario with detailed stellar evolution calculations, starting with ad hoc initial conditions enforcing no core growth during the merger. We find that these massive merger products are likely to be helium-rich and spend most of their remaining lifetime within regions of instabilities in the Herzsprung–Russell diagram, such as luminous blue variable eruptions. An energetic estimate of the amount of mass loss neglecting the back reaction of the star suggests that the total amount of mass that can be removed at low metallicity is ≲1 M ⊙ . This is small enough that at CC our models are retaining sufficient mass to form BHs in the PI gap similar to the recent ones detected by LIGO/Virgo. However, mass loss at the time of merger, the resulting core structure, and the mass loss at CC still need to be quantified for these models to confirm the viability of this scenario. 
    more » « less
  4. 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
  5. ABSTRACT

    Stellar-mass black holes (BHs) can be retained in globular clusters (GCs) until the present. Simulations of GC evolution find that the relaxation driven mass-loss rate is elevated if BHs are present, especially near dissolution. We capture this behaviour in a parametrized mass-loss rate, bench marked by results from N-body simulations, and use it to evolve an initial GC mass function (GCMF), similar to that of young massive clusters in the Local Universe, to an age of 12 Gyr. Low-metallicity GCs ([Fe/H] ≲ −1.5) have the highest mass-loss rates, because of their relatively high BH masses, which combined with their more radial orbits and stronger tidal field in the past explains the high turnover mass of the GCMF ($\sim 10^5\, {\rm M}_\odot$ ) at large Galactic radii ($\gtrsim 10\, {\rm kpc}$ ). The turnover mass at smaller Galactic radii is similar because of the upper mass truncation of the initial GCMF and the lower mass-loss rate due to the higher metallicities. The density profile in the Galaxy of mass lost from massive GCs ($\gtrsim 10^{5}\, {\rm M}_\odot$ ) resembles that of nitrogen-rich stars in the halo, confirming that these stars originated from GCs. We conclude that two-body relaxation is the dominant effect in shaping the GCMF from a universal initial GCMF, because including the effect of BHs reduces the need for additional disruption mechanisms.

     
    more » « less