Impact of massive binary star and cosmic evolution on gravitational wave observations I: black hole–neutron star mergers

• Published in: Monthly notices of the Royal Astronomical Society Vol. 508; no. 4; pp. 5028 - 5063
• Main Authors: Broekgaarden, Floor S, Berger, Edo, Neijssel, Coenraad J, Vigna-Gómez, Alejandro, Chattopadhyay, Debatri, Stevenson, Simon, Chruslinska, Martyna, Justham, Stephen, de Mink, Selma E, Mandel, Ilya
• Format: Journal Article
• Language: English
• Published: Oxford University Press 01.12.2021
• Subjects: stars: evolution; gravitational waves; (transients:) black hole–neutron star mergers
• ISSN: 0035-8711; 1365-2966
• DOI: 10.1093/mnras/stab2716

ABSTRACT Mergers of black hole–neutron star (BHNS) binaries have now been observed by gravitational wave (GW) detectors with the recent announcement of GW200105 and GW200115. Such observations not only provide confirmation that these systems exist but will also give unique insights into the death of massive stars, the evolution of binary systems and their possible association with gamma-ray bursts, r-process enrichment, and kilonovae. Here, we perform binary population synthesis of isolated BHNS systems in order to present their merger rate and characteristics for ground-based GW observatories. We present the results for 420 different model permutations that explore key uncertainties in our assumptions about massive binary star evolution (e.g. mass transfer, common-envelope evolution, supernovae), and the metallicity-specific star formation rate density, and characterize their relative impacts on our predictions. We find intrinsic local BHNS merger rates spanning $\mathcal {R}_{\rm {m}}^0 \approx$ 4–830 $\, \rm {Gpc}^{-3}$$\, \rm {yr}^{-1}$ for our full range of assumptions. This encompasses the rate inferred from recent BHNS GW detections and would yield detection rates of $\mathcal {R}_{\rm {det}} \approx 1$–180$\, \rm {yr}^{-1}$ for a GW network consisting of LIGO, Virgo, and KAGRA at design sensitivity. We find that the binary evolution and metallicity-specific star formation rate density each impacts the predicted merger rates by order $\mathcal {O}(10)$. We also present predictions for the GW-detected BHNS merger properties and find that all 420 model variations predict that $\lesssim 5{{\ \rm per\ cent}}$ of the BHNS mergers have BH masses $m_{\rm {BH}} \gtrsim 18\, \rm {M}_{\odot }$, total masses $m_{\rm {tot}} \gtrsim 20\, \rm {M}_{\odot }$, chirp masses ${\mathcal {M}}_{\rm {c}} \gtrsim 5.5\, \rm {M}_{\odot }$, and mass ratios qf ≳ 12 or qf ≲ 2. Moreover, we find that massive NSs with $m_{\rm {NS}} \gt 2\, \rm {M}_{\odot }$ are expected to be commonly detected in BHNS mergers in almost all our model variations. Finally, a wide range of $\sim 0{{\ \rm per\ cent}}$ to $70{{\ \rm per\ cent}}$ of the BHNS mergers are predicted to eject mass during the merger. Our results highlight the importance of considering variations in binary evolution and cosmological models when predicting, and eventually evaluating, populations of BHNS mergers.