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: Investigating the Cosmological Rate of Compact Object Mergers from Isolated Massive Binary Stars
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
Award ID(s):
2206110
PAR ID:
10653902
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:
24
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. 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
  2. Abstract Recent gravitational wave (GW) observations include possible detections of black hole—neutron star binary mergers. As with binary black hole mergers, numerical simulations help characterize the sources. For binary systems with neutron star components, the simulations help to predict the imprint of tidal deformations and disruptions on the GW signals. In a previous study, we investigated how the mass of the black hole has an impact on the disruption of the neutron star and, as a consequence, on the shape of the GWs emitted. We extend these results to study the effects of varying the compactness of the neutron star. We consider neutron star compactness in the 0.113–0.2 range for binaries with mass ratios of 3 and 5. As the compactness and the mass ratio increase, the binary system behaves during the late inspiral and merger more like a black hole binary. For the cases with the least compact neutron star, the GWs emitted, in terms of mismatches, are the most distinguishable from those by a binary black hole. The disruption of the star significantly suppresses the kicks on the final black hole. The disruption also affects, although not dramatically, the spin of the final black hole. Lastly, for neutron stars with low compactness, the quasi-normal ringing of the black hole after the merger does not show a clean quasi-normal ringing because of the late accretion of debris from the neutron star. 
    more » « less
  3. Abstract The delay time distribution of neutron star mergers provides critical insights into binary evolution processes and the merger rate evolution of compact object binaries. However, current observational constraints on this delay time distribution rely on the small sample of Galactic double neutron stars (with uncertain selection effects), a single multimessenger gravitational wave event, and indirect evidence of neutron star mergers based on r -process enrichment. We use a sample of 68 host galaxies of short gamma-ray bursts to place novel constraints on the delay time distribution and leverage this result to infer the merger rate evolution of compact object binaries containing neutron stars. We recover a power-law slope of α = − 1.83 − 0.39 + 0.35 (median and 90% credible interval) with α < −1.31 at 99% credibility, a minimum delay time of t min = 184 − 79 + 67 Myr with t min > 72 Myr at 99% credibility, and a maximum delay time constrained to t max > 7.95 Gyr at 99% credibility. We find these constraints to be broadly consistent with theoretical expectations, although our recovered power-law slope is substantially steeper than the conventional value of α = −1, and our minimum delay time is larger than the typically assumed value of 10 Myr. Pairing this cosmological probe of the fate of compact object binary systems with the Galactic population of double neutron stars will be crucial for understanding the unique selection effects governing both of these populations. In addition to probing a significantly larger redshift regime of neutron star mergers than possible with current gravitational wave detectors, complementing our results with future multimessenger gravitational wave events will also help determine if short gamma-ray bursts ubiquitously result from compact object binary mergers. 
    more » « less
  4. null (Ed.)
    Neutron stars (NSs) are extraordinary not only because they are the densest form of matter in the visible Universe but also because they can generate magnetic fields ten orders of magnitude larger than those currently constructed on earth. The combination of extreme gravity with the enormous electromagnetic (EM) fields gives rise to spectacular phenomena like those observed on August 2017 with the merger of a binary neutron star system, an event that generated a gravitational wave (GW) signal, a short γ-ray burst (sGRB), and a kilonova. This event serves as the highlight so far of the era of multimessenger astronomy. In this review, we present the current state of our theoretical understanding of compact binary mergers containing NSs as gleaned from the latest general relativistic magnetohydrodynamic simulations. Such mergers can lead to events like the one on August 2017, GW170817, and its EM counterparts, GRB 170817 and AT 2017gfo. In addition to exploring the GW emission from binary black hole-neutron star and neutron star-neutron star mergers, we also focus on their counterpart EM signals. In particular, we are interested in identifying the conditions under which a relativistic jet can be launched following these mergers. Such a jet is an essential feature of most sGRB models and provides the main conduit of energy from the central object to the outer radiation regions. Jet properties, including their lifetimes and Poynting luminosities, the effects of the initial magnetic field geometries and spins of the coalescing NSs, as well as their governing equation of state, are discussed. Lastly, we present our current understanding of how the Blandford-Znajek mechanism arises from merger remnants as the trigger for launching jets, if, when and how a horizon is necessary for this mechanism, and the possibility that it can turn on in magnetized neutron ergostars, which contain ergoregions, but no horizons. 
    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