91̽��

Abstract:

The black hole (BH) binaries in active galactic nuclei (AGN) are expected to form mainly through scattering encounters in the ambient gaseous medium. Recent simulations, including our own, have confirmed this formation pathway is highly efficient. We perform 3D smoothed particle hydrodynamics (SPH) simulations of BH scattering encounters in AGN discs. Using a range of impact parameters, we probe the necessary conditions for binary capture and how different orbital trajectories affect the dissipative effects from the gas. We identify a single range of impact parameters, typically of width ∼0.86−1.59 binary Hill radii depending on AGN disc density, that reliably leads to binary formation. The periapsis of the first encounter is the primary variable that determines the outcome of the initial scattering. We find an associated power law between the energy dissipated and the periapsis depth to be ΔE ∝ r^−b with b = 0.42 ± 0.16, where deeper encounters dissipate more energy. Excluding accretion physics does not significantly alter these results. We identify the region of parameter space in initial energy versus impact parameter where a scattering leads to binary formation. Based on our findings, we provide a ready-to-use analytic criterion that utilizes these two pre-encounter parameters to determine the outcome of an encounter, with a reliability rate of >90 per cent. As the criterion is based directly on our simulations, it provides a reliable and highly physically motivated criterion for predicting binary scattering outcomes which can be used in population studies of BH binaries and mergers around AGN.

Abstract:

Galactic nuclei, the densest stellar environments in the Universe, exhibit a complex geometrical structure. The stars orbiting the central supermassive black hole follow a mass segregated distribution both in the radial distance from the center and in the inclination angle of the orbital planes. The latter distribution may represent the equilibrium state of vector resonant relaxation. In this paper, we build simple models to understand the equilibrium distribution found previously in numerical simulations. Using the method of maximizing the total entropy and the quadrupole mean-field approximation, we determine the equilibrium distribution of axisymmetric two-component gravitating systems with two distinct masses, semimajor axes, and eccentricities. We also examine the limiting case when one of the components dominates over the total energy and angular momentum, approximately acting as a heat bath, which may represent the surrounding astrophysical environment such as the tidal perturbation from the galaxy, a massive perturber, a gas torus, or a nearby stellar system. Remarkably, the bodies above a critical mass in the subdominant component condense into a disk in a ubiquitous way. We identify the system parameters where the transition is smooth and where it is discontinuous. The latter cases exhibit a phase transition between an ordered disklike state and a disordered nearly spherical distribution both in the canonical and in the microcanonical ensembles for these long-range interacting systems.