Bence Kocsis

Abstract:

Abstract Merging black hole (BH) binaries in AGN discs formed through two-body scatterings via the “gas-capture” process may explain a significant fraction of BH mergers in AGN and a non-negligible contribution to the observed rate from LIGO-VIRGO-KAGRA. We perform Monte Carlo simulations of binary BH formation, evolution and mergers across the observed AGN mass function using a novel physically motivated treatment for the gas-capture process derived from hydrodynamical simulations of BH-BH encounters in AGN. Our models suggest that gas-captured binaries could result in merger rates of 0.73 − 7.1Gpc−3yr−1. Mergers from AGN are dominated by AGN with supermassive BH masses of ∼107M⊙, with 90 % of mergers occurring in the range ∼106M⊙ − 108M⊙. The merging mass distribution is flatter than the initial BH mass power law by a factor Δξ = 1.1 − 1.2, as larger BHs align with the disc and form binaries more efficiently. Similarly, the merging mass ratio distribution is flatter, therefore the AGN channel could explain high mass and unequal mass ratio detections such as GW190521 and GW190814. Using a simpler dynamical friction treatment for the binary formation process, the results are similar, where the primary bottleneck is the alignment time with the disc. The most influential parameters are the anticipated number of BHs and their mass function. Given the many uncertainties that remain in the AGN channel, we expect the true uncertainty extends beyond our predicted rates. Nonetheless, we conclude that AGN remain an important channel for consideration, particularly for gravitational wave detections involving one or two high mass BHs.

Abstract:

ABSTRACT We investigate the evolution of black holes on orbits with small inclinations ($i < 2^\circ$) to the gaseous discs of active galactic nuclei (AGNs). We perform 3D adiabatic hydrodynamic simulations within a shearing frame, studying the damping of inclination by black hole-gas gravitation. We find that for objects with $i< 3H_0R_0^{-1}$, where $H_0R_0^{-1}$ is the disc aspect ratio, the inclination lost per mid-plane crossing is proportional to the inclination preceding the crossing, resulting in a net exponential decay in inclination. For objects with $i>3H_0R_0^{-1}$, damping efficiency decreases for higher inclinations. We consider a variety of different AGN environments, finding that damping is stronger for systems with a higher ambient Hill mass: the initial gas mass within the BH sphere of influence. We provide a fitting formula for the inclination changes as a function of Hill mass. We find reasonable agreement between the damping driven by gas gravity in the simulations and the damping driven by accretion under a Hill-limited Bondi–Hoyle–Lyttleton prescription. We find that gas dynamical friction consistently overestimates the strength of damping, especially for lower inclination systems, by at least an order of magnitude. For regions in the AGN disc where coplanar binary black hole formation by gas dissipation is efficient, we find that the simulated damping time-scales are especially short with $\tau _d < 10P_\mathrm{SMBH}$. We conclude that as the time-scales for inclination damping are shorter than the expected interaction time between isolated black holes, the vast majority of binaries formed from gas capture should form from components with negligible inclination to the AGN disc.

Abstract:

The frequency of compact object interactions in AGN discs is naturally tied to the number of objects embedded within it. We investigate the evolution of black holes in the nuclear stellar cluster on inclined orbits to the AGN disc by performing adiabatic hydrodynamical simulations of isolated black hole disc crossings over a range of disc densities and inclinations . We find radiation dominates the pressure in the wake that forms around the BH across the full inclination and disc density range. We identify no well defined steady state wake morphology due to the thin geometry of the disc and the vertical exponential density drop off, where the wake morphology depends on the vertical depth of the transit within the disc. The inclination damping relative the pre-transit inclination behaves as a power law in and the ambient Hill mass as . The drag on the BH is dominated by the gravity of the wake for the majority of our inclination range until accretion effects become comparable at , where is the disc aspect ratio. At low inclinations () the wake morphology becomes more spherical, leading to a regime change in the inclination damping behaviour. Our results suggest that the inclination damping time-scale is shorter than expected from only episodic Bondi–Hoyle–Lyttelton accretion events during each transit, implying inclined objects may be captured by the AGN disc earlier in its lifetime than previously thought.

Abstract:

We investigate close encounters between initially unbound black holes (BHs) in the gaseous discs of active galactic nuclei (AGNs), performing the first 3D non-isothermal hydrodynamical simulations of gas-assisted binary BH formation. We discuss a suite of 135 simulations, considering nine AGN disc environments and 15 BH impact parameters. We find that the gas distribution within the Hill sphere about an isolated embedded BH is akin to a spherically symmetric star with a low-mass convective envelope and a BH core, with large convective currents driving strong outflows away from the mid-plane. We find that Coriolis force acting on the outflow results in winds, analogous to cyclones, that counter-rotate with respect to the mid-plane flow within the Hill sphere. We confirm the existence of strong thermal blasts due to minidisc collisions during BH close encounters, as predicted in our previous 2D studies. We document binary formation across a wide range of environments, finding formation likelihood is increased when the gas mass in the Hill sphere is large, allowing for easier binary formation in the outer AGN disc. We provide a comprehensive overview of the supermassive black hole’s role in binary formation, investigating how binary formation in intermediate density environments is biased towards certain binary orientations. We offer two models for predicting dissipation by gas during close encounters, as a function of the ambient Hill mass alone, or with the periapsis depth. We use these models to motivate a prescription for binary formation likelihood that can be readily applied to Monte Carlo simulations of AGN evolution.