skip to main content

Title: Monte Carlo simulations of black hole mergers in AGN discs: Low χeff mergers and predictions for LIGO
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 are more » relatively short-lived (≤1 Myr) so multiple AGN episodes can happen per Galactic nucleus in a Hubble time. « less
; ; ;
Award ID(s):
Publication Date:
Journal Name:
Monthly Notices of the Royal Astronomical Society
Page Range or eLocation-ID:
1203 to 1216
Sponsoring Org:
National Science Foundation
More Like this
  1. 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 aremore »$\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 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
  2. 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}}}$,more »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$ 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
  3. ABSTRACT Direct collapse black holes (BHs) are promising candidates for producing massive z ≳ 6 quasars, but their formation requires fine-tuned conditions. In this work, we use cosmological zoom simulations to study systematically the impact of requiring: (1) low gas angular momentum (spin), and (2) a minimum incident Lyman–Werner (LW) flux in order to form BH seeds. We probe the formation of seeds (with initial masses of $M_{\rm seed} \sim 10^4\!-\!10^6\, \mathrm{M}_{\odot }\, h^{-1})$ in haloes with a total mass >3000 × Mseed and a dense, metal-poor gas mass >5 × Mseed. Within this framework, we find that the seed-forming haloes have a prior history ofmore »star formation and metal enrichment, but they also contain pockets of dense, metal-poor gas. When seeding is further restricted to haloes with low gas spins, the number of seeds formed is suppressed by factors of ∼6 compared to the baseline model, regardless of the seed mass. Seed formation is much more strongly impacted if the dense, metal-poor gas is required to have a critical LW flux (Jcrit). Even for Jcrit values as low as 50J21, no $8\times 10^{5}~\mathrm{M}_{\odot }\, h^{-1}$ seeds are formed. While lower mass ($1.25\times 10^{4},1\times 10^{5}~\mathrm{M}_{\odot }\, h^{-1}$) seeds do form, they are strongly suppressed (by factors of ∼10–100) compared to the baseline model at gas mass resolutions of $\sim 10^4~\mathrm{M}_{\odot }\, h^{-1}$ (with even stronger suppression at higher resolutions). As a result, BH merger rates are also similarly suppressed. Since early BH growth is dominated by mergers in our models, none of the seeds are able to grow to the supermassive regime ($\gtrsim 10^6~\mathrm{M}_{\odot }\, h^{-1}$) by z = 7. Our results hint that producing the bulk of the z ≳ 6 supermassive BH population may require alternate seeding scenarios that do not depend on the LW flux, early BH growth dominated by rapid or super-Eddington accretion, or a combination of these possibilities.« less

    We use the simba cosmological galaxy formation simulation to investigate the relationship between major mergers ($\lesssim$4:1), starbursts, and galaxy quenching. Mergers are identified via sudden jumps in stellar mass M* well above that expected from in situ star formation, while quenching is defined as going from specific star formation rate (sSFR) $\gt t_{\rm H}^{-1}$ to $\lt 0.2t_{\rm H}^{-1}$, where tH is the Hubble time. At z ≈ 0–3, mergers show ∼2–3× higher SFR than a mass-matched sample of star-forming galaxies, but globally represent $\lesssim 1{{\ \rm per\ cent}}$ of the cosmic SF budget. At low masses, the increase in SFRmore »in mergers is mostly attributed to an increase in the H2 content, but for $M_*\gtrsim 10^{10.5} \,\mathrm{ M}_{\odot }$ mergers also show an elevated star formation efficiency suggesting denser gas within merging galaxies. The merger rate for star-forming galaxies shows a rapid increase with redshift, ∝(1 + z)3.5, but the quenching rate evolves much more slowly, ∝(1 + z)0.9; there are insufficient mergers to explain the quenching rate at $z\lesssim 1.5$. simba first quenches galaxies at $z\gtrsim 3$, with a number density in good agreement with observations. The quenching time-scales τq are strongly bimodal, with ‘slow’ quenchings (τq ∼ 0.1tH) dominating overall, but ‘fast’ quenchings (τq ∼ 0.01tH) dominating in M* ∼ 1010–1010.5 M$\odot$ galaxies, likely induced by simba’s jet-mode black hole feedback. The delay time distribution between mergers and quenching events suggests no physical connection to either fast or slow quenching. Hence, simba predicts that major mergers induce starbursts, but are unrelated to quenching in either fast or slow mode.

    « less
  5. ABSTRACT In this work, we establish and test methods for implementing dynamical friction (DF) for massive black hole pairs that form in large volume cosmological hydrodynamical simulations that include galaxy formation and black hole growth. We verify our models and parameters both for individual black hole dynamics and for the black hole population in cosmological volumes. Using our model of DF from collisionless particles, black holes can effectively sink close to the galaxy centre, provided that the black hole’s dynamical mass is at least twice that of the lowest mass resolution particles in the simulation. Gas drag also plays amore »role in assisting the black holes’ orbital decay, but it is typically less effective than that from collisionless particles, especially after the first billion years of the black hole’s evolution. DF from gas becomes less than $1{{\ \rm per\ cent}}$ of DF from collisionless particles for BH masses >107 M⊙. Using our best DF model, we calculate the merger rate down to z = 1.1 using an Lbox = 35 Mpc h−1 simulation box. We predict ∼2 mergers per year for z > 1.1 peaking at z ∼ 2. These merger rates are within the range obtained in previous work using similar resolution hydrodynamical simulations. We show that the rate is enhanced by factor of ∼2 when DF is taken into account in the simulations compared to the no-DF run. This is due to ${\gt}40{{\ \rm per\ cent}}$ more black holes reaching the centre of their host halo when DF is added.« less