More Like this

1. Monte Carlo simulations of black hole mergers in AGN discs: Low χeff mergers and predictions for LIGO

   https://doi.org/10.1093/mnras/staa740  

   McKernan, B ; Ford, K E ; O’Shaugnessy, R ; Wysocki, D ( May 2020 , Monthly Notices of the Royal Astronomical Society)

   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
2. Binary black hole merger rates in AGN discs versus nuclear star clusters: loud beats quiet

   https://doi.org/10.1093/mnras/stac2861  

   Ford, K. E. Saavik ; McKernan, Barry ( October 2022 , Monthly Notices of the Royal Astronomical Society)

   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. Signatures of hierarchical mergers in black hole spin and mass distribution

   https://doi.org/10.1093/mnras/stab2315  

   Tagawa, Hiromichi ; Haiman, Zoltán ; Bartos, Imre ; Kocsis, Bence ; Omukai, Kazuyuki ( September 2021 , Monthly Notices of the Royal Astronomical Society)

   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. The landscape of disc outflows from black hole–neutron star mergers

   https://doi.org/10.1093/mnras/staa2209  

   Fernández, Rodrigo ; Foucart, Francois ; Lippuner, Jonas ( September 2020 , Monthly Notices of the Royal Astronomical Society)

   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
5. Globular cluster formation histories, masses, and radii inferred from gravitational waves

   https://doi.org/10.1093/mnras/stad1364  

   Fishbach, Maya ; Fragione, Giacomo ( May 2023 , Monthly Notices of the Royal Astronomical Society)

   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