skip to main content

This content will become publicly available on July 1, 2023

Title: Constraining High-redshift Stellar-mass Primordial Black Holes with Next-generation Ground-based Gravitational-wave Detectors
Abstract The possible existence of primordial black holes in the stellar-mass window has received considerable attention because their mergers may contribute to current and future gravitational-wave detections. Primordial black hole mergers, together with mergers of black holes originating from Population III stars, are expected to dominate at high redshifts ( z ≳ 10). However, the primordial black hole merger rate density is expected to rise monotonically with redshift, while Population III mergers can only occur after the birth of the first stars. Next-generation gravitational-wave detectors such as the Cosmic Explorer (CE) and Einstein Telescope (ET) can access this distinctive feature in the merger rates as functions of redshift, allowing for direct measurement of the abundance of the two populations and hence for robust constraints on the abundance of primordial black holes. We simulate four months’ worth of data observed by a CE-ET detector network and perform hierarchical Bayesian analysis to recover the merger rate densities. We find that if the universe has no primordial black holes with masses of  ( 10 M ⊙ ) , the projected upper limit on their abundance f PBH as a fraction of dark matter energy density may be as low as f PBH more » ∼  ( 10 − 5 ) , about two orders of magnitude lower than the current upper limits in this mass range. If instead f PBH ≳ 10 −4 , future gravitational-wave observations would exclude f PBH = 0 at the 95% credible interval. « less
; ; ; ; ;
Award ID(s):
Publication Date:
Journal Name:
The Astrophysical Journal Letters
Page Range or eLocation-ID:
Sponsoring Org:
National Science Foundation
More Like this
  1. The assembly of massive black holes in the early universe remains a poorly constrained open question in astrophysics. The merger and accretion of light seeds (remnants of Population III stars with mass below ∼ 1000 M ) or heavy seeds (in the mass range 104−106 M ) could both explain the formation of massive black holes, but the abundance of seeds and their merging mechanism are highly uncertain. In the next decades, the gravitational-wave observatories coming online are expected to observe very highredshift mergers, shedding light on the seeding of the first black holes. In this Letter we explore the potential and limitations for LISA, Cosmic Explorer and Einstein Telescope to constrain the mixture ratio of light and heavy seeds as well as the probability that central black holes in merging galaxies merge as well. Since the third generation ground-based gravitational-wave detectors will only observe light seed mergers, we demonstrate two scenarios in which the inference of the seed mixture ratio and merging probability can be limited. The synergy of multi-band gravitational-wave observations and electromagnetic observations will likely be necessary in order to fully characterize the process of high-redshift black hole formation.
  2. Long-duration gamma-ray bursts are thought to be associated with the core-collapse of massive, rapidly spinning stars and the formation of black holes. However, efficient angular momentum transport in stellar interiors, currently supported by asteroseismic and gravitational-wave constraints, leads to predominantly slowly-spinning stellar cores. Here, we report on binary stellar evolution and population synthesis calculations, showing that tidal interactions in close binaries not only can explain the observed subpopulation of spinning, merging binary black holes but also lead to long gamma-ray bursts at the time of black-hole formation. Given our model calibration against the distribution of isotropic-equivalent energies of luminous long gamma-ray bursts, we find that ≈10% of the GWTC-2 reported binary black holes had a luminous long gamma-ray burst associated with their formation, with GW190517 and GW190719 having a probability of ≈85% and ≈60%, respectively, being among them. Moreover, given an assumption about their average beaming fraction, our model predicts the rate density of long gamma-ray bursts, as a function of redshift, originating from this channel. For a constant beaming fraction f B  ∼ 0.05 our model predicts a rate density comparable to the observed one, throughout the redshift range, while, at redshift z  ∈ [0, 2.5], a tentative comparison with the metallicitymore »distribution of observed LGRB host galaxies implies that between 20% to 85% of the observed long gamma-ray bursts may originate from progenitors of merging binary black holes. The proposed link between a potentially significant fraction of observed, luminous long gamma-ray bursts and the progenitors of spinning binary black-hole mergers allows us to probe the latter well outside the horizon of current-generation gravitational wave observatories, and out to cosmological distances.« less
  3. 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. Wemore »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).« less
  4. Abstract

    Gravitational-wave detectors are starting to reveal the redshift evolution of the binary black hole (BBH) merger rate,RBBH(z). We make predictions forRBBH(z) as a function of black hole mass for systems originating from isolated binaries. To this end, we investigate correlations between the delay time and black hole mass by means of the suite of binary population synthesis simulations,COMPAS. We distinguish two channels: the common envelope (CE), and the stable Roche-lobe overflow (RLOF) channel, characterized by whether the system has experienced a common envelope or not. We find that the CE channel preferentially produces BHs with masses below about 30Mand short delay times (tdelay≲ 1 Gyr), while the stable RLOF channel primarily forms systems with BH masses above 30Mand long delay times (tdelay≳ 1 Gyr). We provide a new fit for the metallicity-dependent specific star formation rate density based on the Illustris TNG simulations, and use this to convert the delay time distributions into a prediction ofRBBH(z). This leads to a distinct redshift evolution ofRBBH(z) for high and low primary BH masses. We furthermore find that, at high redshift,RBBH(z) is dominated by the CE channel, while at low redshift, it contains a large contribution (∼40%) from the stable RLOF channel.more »Our results predict that, for increasing redshifts, BBHs with component masses above 30Mwill become increasingly scarce relative to less massive BBH systems. Evidence of this distinct evolution ofRBBH(z) for different BH masses can be tested with future detectors.

    « less
  5. Abstract

    We describe the public release of the Cluster Monte Carlo (CMC) code, a parallel, star-by-starN-body code for modeling dense star clusters.CMCtreats collisional stellar dynamics using Hénon’s method, where the cumulative effect of many two-body encounters is statistically reproduced as a single effective encounter between nearest-neighbor particles on a relaxation timescale. The star-by-star approach allows for the inclusion of additional physics, including strong gravitational three- and four-body encounters, two-body tidal and gravitational-wave captures, mass loss in arbitrary galactic tidal fields, and stellar evolution for both single and binary stars. The public release ofCMCis pinned directly to theCOSMICpopulation synthesis code, allowing dynamical star cluster simulations and population synthesis studies to be performed using identical assumptions about the stellar physics and initial conditions. As a demonstration, we present two examples of star cluster modeling: first, we perform the largest (N= 108) star-by-starN-body simulation of a Plummer sphere evolving to core collapse, reproducing the expected self-similar density profile over more than 15 orders of magnitude; second, we generate realistic models for typical globular clusters, and we show that their dynamical evolution can produce significant numbers of black hole mergers with masses greater than those produced from isolated binary evolution (such as GW190521, amore »recently reported merger with component masses in the pulsational pair-instability mass gap).

    « less