The detection of the binary black hole merger GW190521, with primary black hole mass $85^{+21}_{-14} {\rm M}_{\odot }$, proved the existence of black holes in the theoretically predicted pair-instability gap ($\sim 60-120 \, {\rm M}_{\odot }$) of their mass spectrum. Some recent studies suggest that such massive black holes could be produced by the collision of an evolved star with a carbon–oxygen core and a main sequence star. Such a post-coalescence star could end its life avoiding the pair-instability regime and with a direct collapse of its very massive envelope. It is still not clear, however, how the collision shapes the structure of the newly produced star and how much mass is actually lost in the impact. We investigated this issue by means of hydrodynamical simulations with the smoothed particle hydrodynamics code StarSmasher, finding that a head-on collision can remove up to 12 per cent of the initial mass of the colliding stars. This is a non-negligible percentage of the initial mass and could affect the further evolution of the stellar remnant, particularly in terms of the final mass of a possibly forming black hole. We also found that the main sequence star can plunge down to the outer boundary of the core of the primary, changing the inner chemical composition of the remnant. The collision expels the outer layers of the primary, leaving a remnant with an helium-enriched envelope (reaching He fractions of about 0.4 at the surface). These more complex abundance profiles can be directly used in stellar evolution simulations of the collision product.
Formation of black holes in the pair-instability mass gap: Hydrodynamical simulation of a head-on massive star collision
The detection of the binary black hole merger GW190521, with primary mass 85+21−14 M⊙ , proved the existence of black holes in the theoretically predicted pair-instability gap ( ∼60−120M⊙ ) of their mass spectrum. Some recent studies suggest that such massive black holes could be produced by the collision of an evolved star with a carbon-oxygen core and a main sequence star. Such a post-coalescence star could end its life avoiding the pair-instability regime and with a direct collapse of its very massive envelope. It is still not clear, however, how the collision shapes the structure of the newly produced star and how much mass is actually lost in the impact. We investigated this issue by means of hydrodynamical simulations with the smoothed particle hydrodynamics code StarSmasher, finding that a head-on collision can remove up to 12% of the initial mass of the colliding stars. This is a non-negligible percentage of the initial mass and could affect the further evolution of the stellar remnant, particularly in terms of the final mass of a possibly forming black hole. We also found that the main sequence star can plunge down to the outer boundary of the carbon-oxygen core of the primary, changing the inner chemical composition of the remnant. The collision expels the outer layers of the primary, leaving a remnant with an helium-enriched envelope (reaching He fractions of about 0.4 at the surface). These more complex abundance profiles can be directly used in stellar evolution simulations of the collision product.
more »
« less
- Award ID(s):
- 1909203
- NSF-PAR ID:
- 10358797
- Date Published:
- Journal Name:
- ArXivorg
- ISSN:
- 2331-8422
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
more » « less
ABSTRACT -
null (Ed.)Context. Grids of stellar models, computed with the same physical ingredients, allow one to study the impact of a given physics on a broad range of initial conditions and they are a key ingredient for modeling the evolution of galaxies. Aims. We present here a grid of single star models for masses between 0.8 and 120 M ⊙ , with and without rotation for a mass fraction of heavy element Z = 0.006, representative of the Large Magellanic Cloud (LMC). Methods. We used the GENeva stellar Evolution Code. The evolution was computed until the end of the central carbon-burning phase, the early asymptotic giant branch phase, or the core helium-flash for massive, intermediate, and low mass stars, respectively. Results. The outputs of the present stellar models are well framed by the outputs of the two grids obtained by our group for metallicities above and below the one considered here. The models of the present work provide a good fit to the nitrogen surface enrichments observed during the main sequence for stars in the LMC with initial masses around 15 M ⊙ . They also reproduce the slope of the luminosity function of red supergiants of the LMC well, which is a feature that is sensitive to the time-averaged mass loss rate over the red supergiant phase. The most massive black hole that can be formed from the present models at Z = 0.006 is around 55 M ⊙ . No model in the range of mass considered will enter into the pair-instability supernova regime, while the minimal mass to enter the region of pair pulsation instability is around 60 M ⊙ for the rotating models and 85 M ⊙ for the nonrotating ones. Conclusions. The present models are of particular interest for comparisons with observations in the LMC and also in the outer regions of the Milky Way. We provide public access to numerical tables that can be used for computing interpolated tracks and for population synthesis studies.more » « less
-
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
-
Abstract We study close encounters of a 1 M ⊙ middle-age main-sequence star (modeled using MESA) with massive black holes through hydrodynamic simulations, and explore in particular the dependence of the outcomes on the black hole mass. We consider here black holes in the intermediate-mass range, M BH = 100–10 4 M ⊙ . Possible outcomes vary from a small tidal perturbation for weak encounters all the way to partial or full disruption for stronger encounters. We find that stronger encounters lead to increased mass loss at the first pericenter passage, in many cases ejecting the partially disrupted star on an unbound orbit. For encounters that initially produce a bound system, with only partial stripping of the star, the fraction of mass stripped from the star increases with each subsequent pericenter passage and a stellar remnant of finite mass is ultimately ejected in all cases. The critical penetration depth that separates bound and unbound remnants has a dependence on the black hole mass when M BH ≲ 10 3 M ⊙ . We also find that the number of successive close passages before ejection decreases as we go from the stellar-mass black hole to the intermediate-mass black hole regime. For instance, after an initial encounter right at the classical tidal disruption limit, a 1 M ⊙ star undergoes 16 (5) pericenter passages before ejection from a 10 M ⊙ (100 M ⊙ ) black hole. Observations of periodic flares from these repeated close passages could in principle indicate signatures of a partial tidal disruption event.more » « less
-
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
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- The DOI is not currently available.
