skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: The Binary Black Hole Merger Rate Deviates from the Cosmic Star Formation Rate: A Tug of War between Metallicity and Delay Times
Abstract Gravitational-wave detectors are now making it possible to investigate how the merger rate of binary black holes (BBHs) evolves with redshift. In this study, we examine whether the BBH merger rate of isolated binaries deviates from a scaled star formation rate density (SFRD)—a frequently used model in state-of-the-art research. To address this question, we conduct population synthesis simulations using COMPAS with a grid of stellar evolution models, calculate their cosmological merger rates, and compare them to a scaled SFRD. We find that our simulated rates deviate by factors up to 3.5 atz∼ 0 and 5 atz∼ 9 due to two main phenomena: (i) the formation efficiency of BBHs is an order of magnitude higher at low metallicities than at solar metallicity, and (ii) BBHs experience a wide range of delays (from a few megayears to many gigayears) between formation and merger. The deviations are similar when comparing to a delayed SFRD, and even larger (up to ∼10×) when comparing to SFRD-based models scaled to the local merger rate. Interestingly, our simulations find that the BBH delay time distribution is redshift dependent, increasing the complexity of the redshift distribution of mergers. We find similar results for simulated merger rates of black hole–neutron stars (BHNSs) and binary neutron stars (BNSs). We conclude that the rate of BBH, BHNS, and BNS mergers from the isolated channel can significantly deviate from a scaled SFRD, and that future measurements of the merger rate will provide insights into the formation pathways of gravitational-wave sources.  more » « less
Award ID(s):
2206110
PAR ID:
10653903
Author(s) / Creator(s):
; ;
Publisher / Repository:
The Astrophysical Journal
Date Published:
Journal Name:
The Astrophysical Journal
Volume:
976
Issue:
1
ISSN:
0004-637X
Page Range / eLocation ID:
23
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract Gravitational-wave (GW) detectors are observing compact object mergers from increasingly far distances, revealing the redshift evolution of the binary black hole (BBH)—and soon the black hole–neutron star (BHNS) and binary neutron star (BNS)—merger rate. To help interpret these observations, we investigate the expected redshift evolution of the compact object merger rate from the isolated binary evolution channel. We present a publicly available catalog of compact object mergers and their accompanying cosmological merger rates from population synthesis simulations conducted with the COMPAS software. To explore the impact of uncertainties in stellar and binary evolution, our simulations use two-parameter grids of binary evolution models that vary the common-envelope efficiency with mass transfer accretion efficiency and supernova (SN) remnant mass prescription with SN natal kick velocity, respectively. We quantify the redshift evolution of our simulated merger rates using the local (z∼ 0) rate, the redshift at which the merger rate peaks, and the normalized differential rates (as a proxy for slope). We find that although the local rates span a range of ∼103across our model variations, their redshift evolutions are remarkably similar for BBHs, BHNSs, and BNSs, with differentials typically within a factor 3 and peaks ofz≈ 1.2–2.4 across models. Furthermore, several trends in our simulated rates are correlated with the model parameters we explore. We conclude that future observations of the redshift evolution of the compact object merger rate can help constrain binary models for stellar evolution and GW formation channels. 
    more » « less
  2. 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. 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. 
    more » « less
  3. Context.Despite the nearly hundred gravitational-wave detections reported by the LIGO-Virgo-KAGRA Collaboration, the question of the cosmological origin of merging binary black holes (BBHs) remains open. The two main formation channels generally considered are from isolated field binaries or via dynamical assembly in dense star clusters. Aims.Here we focus on understanding the dynamical formation of merging BBHs within massive clusters in galaxies of different masses. Methods.To this end, we applied a new framework to consistently model the formation and evolution of massive star clusters in zoom-in cosmological simulations of galaxies. Each simulation, taken from the FIRE project, provides a realistic star formation environment, with a unique star formation history, that hosts realistic giant molecular clouds that constitute the birthplace of star clusters. Combined with the code for star cluster evolutionCMC, we are able to produce populations of dynamically formed merging BBHs across cosmic time in different environments. Results.As the most massive star clusters preferentially form in dense massive clouds of gas, we find that, despite their low metallicities favouring the creation of black holes, low-mass galaxies contain few massive clusters and therefore make a limited contribution to the global production of dynamically formed merging BBHs. Furthermore, we find that massive clusters can host hierarchical BBH mergers with clear, identifiable physical properties. Looking at the evolution of the BBH merger rate in different galaxies, we find strong correlations between BBH mergers and the most extreme episodes of star formation. Finally, we discuss the implications for future LIGO-Virgo-KAGRA gravitational wave observations. 
    more » « less
  4. ABSTRACT Galactic nuclei are promising sites for stellar origin black hole (BH) mergers, as part of merger hierarchies in deep potential wells. We show that binary black hole (BBH) merger rates in active galactic nuclei (AGNs) should always exceed merger rates in quiescent galactic nuclei (nuclear star clusters, NSCs) around supermassive black holes (SMBHs) without accretion discs. This is primarily due to average binary lifetimes in AGNs that are significantly shorter than those in NSCs. The lifetime difference comes from rapid hardening of BBHs in AGNs, such that their semimajor axes are smaller than the hard–soft boundary of their parent NSC; this contrasts with the large average lifetime to merger for BBHs in NSCs around SMBHs, due to binary ionization mechanisms. Secondarily, merger rates in AGNs are enhanced by gas-driven binary formation mechanisms. Formation of new BHs in AGN discs is a minor contributor to the rate differences. With the gravitational wave detection of several BBHs with at least one progenitor in the upper mass gap, and signatures of dynamical formation channels in the χeff distribution, we argue that AGNs could contribute $$\sim 25{\!-\!}80{{\ \rm per\ cent}}$$ of the LIGO–Virgo measured rate of $$\sim 24\, \rm {Gpc}^{-3} \rm {yr}^{-1}$$. 
    more » « less
  5. Abstract In their most recent observing run, the LIGO-Virgo-KAGRA Collaboration observed gravitational waves from compact binary mergers with highly asymmetric mass ratios, including both binary black holes (BBHs) and neutron star-black holes (NSBHs). It appears that NSBHs with mass ratiosq≃ 0.2 are more common than equally asymmetric BBHs, but the reason for this remains unclear. We use the binary population synthesis codecosmicto investigate the evolutionary pathways leading to the formation and merger of asymmetric compact binaries. We find that within the context of isolated binary stellar evolution, most asymmetric mergers start off as asymmetric stellar binaries. Because of the initial asymmetry, these systems tend to first undergo a dynamically unstable mass transfer phase. However, after the first star collapses into a compact object, the mass ratio is close to unity and the second phase of mass transfer is usually stable. According to our simulations, this stable mass transfer fails to shrink the orbit enough on its own for the system to merge. Instead, the natal kick received by the second-born compact object during its collapse is key in determining how many of these systems can merge. For the most asymmetric systems with mass ratios ofq≤ 0.1, the merging systems in our models receive an average kick magnitude of 255 km s−1during the second collapse, while the average kick for non-merging systems is 59 km s−1. Because lower mass compact objects, like neutron stars, are expected to receive larger natal kicks than higher mass BHs, this may explain why asymmetric NSBH systems merge more frequently than asymmetric BBH systems. 
    more » « less