skip to main content

Title: Black hole, neutron star, and white dwarf merger rates in AGN discs
ABSTRACT Advanced LIGO and Advanced Virgo are detecting a large number of binary stellar origin black hole (BH) mergers. A promising channel for accelerated BH merger lies in active galactic nucleus (AGN) discs of gas around supermasssive BHs. Here, we investigate the relative number of compact object (CO) mergers in AGN disc models, including BH, neutron stars (NS), and white dwarfs, via Monte Carlo simulations. We find the number of all merger types in the bulk disc grows ∝ t1/3 which is driven by the Hill sphere of the more massive merger component. Median mass ratios of NS–BH mergers in AGN discs are $\tilde{q}=0.07\pm 0.06(0.14\pm 0.07)$ for mass functions (MF) M−1(− 2). If a fraction fAGN of the observed rate of BH–BH mergers (RBH–BH) come from AGN, the rate of NS–BH (NS–NS) mergers in the AGN channel is ${R}_{\mathrm{ BH}\!-\!\mathrm{ NS}} \sim f_{\mathrm{ AGN}}[10,300]\, \rm {Gpc}^{-3}\, \rm {yr}^{-1},({\mathit{ R}}_{NS\!-\!NS} \le \mathit{ f}_{AGN}400\, \rm {Gpc}^{-3}\, \rm {yr}^{-1}$). Given the ratio of NS–NS/BH–BH LIGO search volumes, from preliminary O3 results the AGN channel is not the dominant contribution to observed NS–NS mergers. The number of lower mass gap events expected is a strong function of the nuclear MF and mass segregation efficiency. CO merger ratios more » derived from LIGO can restrict models of MF, mass segregation, and populations embedded in AGN discs. The expected number of electromagnetic (EM) counterparts to NS–BH mergers in AGN discs at z < 1 is $\sim [30,900]\, {\rm {yr}}^{-1}(f_{\mathrm{ AGN}}/0.1)$. EM searches for flaring events in large AGN surveys will complement LIGO constraints on AGN models and the embedded populations that must live in them. « less
; ;
Award ID(s):
Publication Date:
Journal Name:
Monthly Notices of the Royal Astronomical Society
Page Range or eLocation-ID:
4088 to 4094
Sponsoring Org:
National Science Foundation
More Like this
  1. ABSTRACT Accretion discs around supermassive black holes are promising sites for stellar mass black hole mergers detectable with LIGO. Here we present the results of Monte Carlo simulations of black hole mergers within 1-d AGN disc models. For the spin distribution in the disc bulk, key findings are: (1) The distribution of χeff is naturally centred around $\tilde{\chi }_{\rm eff} \approx 0.0$, (2) the width of the χeff distribution is narrow for low natal spins. For the mass distribution in the disc bulk, key findings are: (3) mass ratios $\tilde{q} \sim 0.5\!-\!0.7$, (4) the maximum merger mass in the bulk is $\sim 100\!-\!200\, \mathrm{M}_{\odot }$, (5) $\sim 1{{\ \rm per\ cent}}$ of bulk mergers involve BH $\gt 50\, \mathrm{M}_{\odot }$ with (6) $\simeq 80{{\ \rm per\ cent}}$ of bulk mergers are pairs of first generation BH. Additionally, mergers at a migration trap grow an IMBH with typical merger mass ratios $\tilde{q}\sim 0.1$. Ongoing LIGO non-detections of black holes $\gt 10^{2}\, \mathrm{M}_{\odot }$ puts strong limits on the presence of migration traps in AGN discs (and therefore AGN disc density and structure) as well as median AGN disc lifetime. The highest merger rate occurs for this channel if AGN discs aremore »relatively short-lived (≤1 Myr) so multiple AGN episodes can happen per Galactic nucleus in a Hubble time.« less

    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}$.

  3. ABSTRACT Recent gravitational wave (GW) observations by LIGO/Virgo show evidence for hierarchical mergers, where the merging BHs are the remnants of previous BH merger events. These events may carry important clues about the astrophysical host environments of the GW sources. In this paper, we present the distributions of the effective spin parameter (χeff), the precession spin parameter (χp), and the chirp mass (mchirp) expected in hierarchical mergers. Under a wide range of assumptions, hierarchical mergers produce (i) a monotonic increase of the average of the typical total spin for merging binaries, which we characterize with $\scriptstyle{{\bar{\chi }}_\mathrm{typ}\equiv \overline{(\chi _\mathrm{eff}^2+\chi _\mathrm{p}^2)^{1/2}}}$, up to roughly the maximum mchirp among first-generation (1g) BHs, and (ii) a plateau at ${\bar{\chi }}_\mathrm{typ}\sim 0.6$ at higher mchirp. We suggest that the maximum mass and typical spin magnitudes for 1g BHs can be estimated from ${\bar{\chi }}_\mathrm{typ}$ as a function of mchirp. The GW data observed in LIGO/Virgo O1–O3a prefers an increase in ${\bar{\chi }}_\mathrm{typ}$ at low mchirp, which is consistent with the growth of the BH spin magnitude by hierarchical mergers at ∼2σ confidence. A Bayesian analysis using the χeff, χp, and mchirp distributions suggests that 1g BHs have the maximum mass of ∼15–$30\, {\rm M}_\odot$more »if the majority of mergers are of high-generation BHs (not among 1g–1g BHs), which is consistent with mergers in active galactic nucleus discs and/or nuclear star clusters, while if mergers mainly originate from globular clusters, 1g BHs are favoured to have non-zero spin magnitudes of ∼0.3. We also forecast that signatures for hierarchical mergers in the ${\bar{\chi }}_\mathrm{typ}$ distribution can be confidently recovered once the number of GW events increases to ≳ O(100).« less
  4. ABSTRACT We investigate mass ejection from accretion discs formed in mergers of black holes (BHs) and neutron stars (NSs). The third observing run of the LIGO/Virgo interferometers provided BH–NS candidate events that yielded no electromagnetic (EM) counterparts. The broad range of disc configurations expected from BH–NS mergers motivates a thorough exploration of parameter space to improve EM signal predictions. Here we conduct 27 high-resolution, axisymmetric, long-term hydrodynamic simulations of the viscous evolution of BH accretion discs that include neutrino emission/absorption effects and post-processing with a nuclear reaction network. In the absence of magnetic fields, these simulations provide a lower limit to the fraction of the initial disc mass ejected. We find a nearly linear inverse dependence of this fraction on disc compactness (BH mass over initial disc radius). The dependence is related to the fraction of the disc mass accreted before the ouflow is launched, which depends on the disc position relative to the innermost stable circular orbit. We also characterize a trend of decreasing ejected fraction and decreasing lanthanide/actinide content with increasing disc mass at fixed BH mass. This trend results from a longer time to reach weak freezout and an increasingly dominant role of neutrino absorption at highermore »disc masses. We estimate the radioactive luminosity from the disc outflow alone available to power kilonovae over the range of configurations studied, finding a spread of two orders of magnitude. For most of the BH–NS parameter space, the disc outflow contribution is well below the kilonova mass upper limits for GW190814.« less

    Globular clusters (GCs) are found in all types of galaxies and harbour some of the most extreme stellar systems, including black holes that may dynamically assemble into merging binary black holes (BBHs). Uncertain GC properties, including when they formed, their initial masses and sizes, affect their production rate of BBH mergers. Using the gravitational-wave transient catalogue (GWTC-3), we measure that dynamically assembled BBHs – those that are consistent with isotropic spin directions – make up ${61^{+29}_{-44}\%}$ of the total merger rate, with a local merger rate of ${10.9^{+16.8}_{-9.3}}$ Gpc−3 yr−1 rising to ${58.9^{+149.4}_{-46.0}}$ Gpc−3 yr−1 at z  = 1. We assume that this inferred rate describes the contribution from GCs and compare it against the Cluster Monte Carlo (cmc) simulation catalogue to directly fit for the GC initial mass function, virial radius distribution, and formation history. We find that GC initial masses are consistent with a Schechter function with slope ${\beta _m = -1.9^{+0.8}_{-0.8}}$ . Assuming a mass function slope of βm  = −2 and a mass range between 104–$10^8\, \mathrm{ M}_\odot$ , we infer a GC formation rate at z  = 2 of ${5.0^{+9.4}_{-4.0}}$ Gpc−3 yr−1, or ${2.1^{+3.9}_{-1.7}}\times 10^6\, \mathrm{ M}_\odot$ Gpc−3 yr−1 in terms of mass density.more »We find that the GC formation rate probably rises more steeply than the global star formation rate between z  = 0 and z  = 3 (82 per cent credibility) and implies a local number density that is ${f_\mathrm{ev} = 22.6^{+29.9}_{-16.2}}$ times higher than the observed density of survived GCs. This is consistent with expectations for cluster evaporation, but may suggest that other environments contribute to the rate of BBH mergers with significantly tilted spins.

    « less