Evolution of Supermassive Black Hole Pairs on Inclined Orbits in Post-Merger Galaxies

Theoretical models of the evolution of supermassive black hole (SMBH) pairs in post-merger remnant galaxies are necessary to motivate observational searches for dual active galactic nuclei (AGN) and gravitational wave sources. Studies have explored the dynamical evolution of SMBH pairs under the influence of dynamical friction to calculate pairing times and predict the expected population of dual-AGNs at various redshifts. We formulate a three-dimensional dynamical model of SMBH pairs in the innermost kiloparsec of a post-merger galaxy to investigate the impact of orbital inclination with respect to the galactic disk on pairing times. The SMBH pairs are evolved in 81 different galaxy configurations initialized using a Gauss-Seidel Poisson solver. The dynamics are calculated for 12 distinct initial inclinations ranging from 0 to 75 degrees in each of the galaxies to gauge the impact of inclination on pairing time. Orbits characterized by initial inclinations greater than 20 degrees frequently require longer pairing times when compared to uninclined orbits. Pairing times for orbits with inclinations $\gtrsim 45$ degrees often exceed 14 Gyr. Galaxies with higher mass SMBH pairs and faster rotating disks generally shorten pairing times relative to galaxies with less massive or slower rotating disks when the inclination is $\lesssim 45$ degrees. The model suggests that SMBH pairs that form from mergers at inclinations $\lesssim 20$ degrees are likely progenitors of dual-AGN and gravitational wave sources.

💡 Research Summary

This paper presents a three‑dimensional dynamical study of supermassive black‑hole (SMBH) pairs in the innermost kiloparsec of post‑merger galaxies, focusing on how the inclination of the secondary’s orbit relative to the galactic disk influences the time required for the pair to reach a separation of 10 pc (the “pairing time”). The authors construct axisymmetric galaxy models composed of a stellar disk, a gaseous disk, and a stellar bulge, deliberately omitting a dark‑matter halo because its contribution to the potential within ~1 kpc is negligible. The density profiles follow double‑exponential forms for the stellar disk, an exponential law for the gas disk, and a modified Binney & Tremaine model for the bulge. Key galaxy parameters—central SMBH mass (M₁ = 10⁶, 10⁷, 10⁸ M⊙), central gas number density (n_gd = 100, 200, 300 cm⁻³), gas‑to‑stellar disk mass fraction (f_g = 0.3, 0.5, 0.7), and disk rotation speed as a fraction of the local circular velocity (v_g = 0.3, 0.5, 0.7 v_c)—are varied systematically. The secondary SMBH mass is fixed at M₂ = M₁/9, a ratio motivated by cosmological merger statistics and chosen to keep the primary effectively stationary.

The gravitational potential is computed once for each galaxy configuration using a Gauss‑Seidel relaxation Poisson solver on a 200 × 200 cylindrical grid. To resolve the steep vertical density gradients near the mid‑plane, a refined mesh covering –0.05 R_g < z < 0.05 R_g is over‑laid, with cell heights as small as 1–3 pc depending on the galaxy size. Convergence tolerances are set to 2 × 10⁻⁶ Φ_ε for the outer region and 1 × 10⁻⁷ Φ_ε for the inner region, ensuring that variations in grid resolution alter pairing times by less than 10 %.

The orbital evolution of the secondary SMBH is integrated with a fourth‑order Runge‑Kutta scheme, employing an adaptive timestep limited to 1 % of the instantaneous orbital period to preserve energy and angular‑momentum conservation within 1 %. The equations of motion are solved in cylindrical coordinates (r, φ, z). Dynamical friction (DF) from stars and gas is included. Stellar DF follows the vector formulation of Antonini & Merritt (2012), allowing a straightforward computation of the vertical component from the vertical velocity. Gas DF is more complex because the secondary traverses a stratified medium; the authors adopt the Ostriker (1999) drag law, which depends on the Mach number of the motion relative to the local sound speed. The local sound speed is estimated from a modified Toomre stability criterion, yielding values around 10 km s⁻¹ (≈10⁴ K). The vertical Mach number M_z = |ẋz|/c_s determines the strength of “vertical DF,” and the drag force is expressed as F_g = –4π(GM₂)² ρ_g v_z² I(M_z), where I(M_z) is the piecewise function from Ostriker (1999).

A total of 81 distinct galaxy configurations are combined with 12 initial inclination angles (i₀ = 0°, 2.5°, 5°, 7.5°, 10°, 15°, 20°, 25°, 30°, 45°, 60°, 75°), yielding 972 simulations. The pairing time is defined as the moment when the secondary reaches a radial distance of 10 pc from the primary. Results show a clear inclination dependence: for i₀ > 20°, pairing times increase dramatically; for i₀ ≈ 45° or higher, most models exceed the Hubble time (~14 Gyr), implying that such inclined pairs are unlikely to merge within a cosmic age. Conversely, low‑inclination orbits (i₀ ≤ 20°) typically pair within ≤ 1 Gyr, especially when the primary SMBH is massive (10⁸ M⊙) and the disk rotates rapidly (v_g ≈ 0.7 v_c). Higher gas densities and larger gas fractions also shorten pairing times by enhancing DF, but the benefit diminishes for highly inclined orbits because the secondary spends most of its orbit outside the dense mid‑plane.

The study concludes that SMBH pairs formed with small spin‑orbit misalignments (≤ 20°) are the most promising progenitors of observable dual AGN (dAGN) and future gravitational‑wave sources detectable by space‑based interferometers such as LISA. Massive, fast‑rotating disks further accelerate the pairing process for low‑inclination systems. The authors acknowledge several limitations: the neglect of dark matter, the assumption of a fixed primary, simplified gas thermodynamics, and the finite spatial resolution of the Poisson solver. Future work is suggested to incorporate feedback processes, star formation, and fully cosmological simulations to assess the statistical distribution of inclination angles and their impact on SMBH merger rates.