Secular Dynamics around a Supermassive Black Hole via Multipole Expansion

Abstract

In galactic nuclei, the gravitational potential is dominated by the central supermassive black hole, so stars follow quasi-Keplerian orbits. These orbits are distorted by gravitational forces from other stars, leading to long-term orbital relaxation. The direct numerical study of these processes is challenging because the fast orbital motion imposed by the central black hole requires very small timesteps. An alternative approach, pioneered by Gauss, is to use the secular approximation of smearing out the $N$ stars over their Keplerian orbits, using $K$ nodes along each orbit. In this study we propose three novel improvements to this method. First, we re-formulate the discretisation of the rates of change of the variables describing the orbital states to ensure that all conservation laws are exactly satisfied. Second, we replace the pairwise sum over nodes by a multipole expansion up to order ${\ell }_{\max }$, reducing the overall computational costs from $O(N^2{K}^2)$ to $O(NK{\ell }_{\max }^2)$. Finally, we show that the averaged dynamical system is equivalent to $2N$ interacting unit spin vectors and provide two time integrators: a second-order symplectic scheme and a fourth-order Lie-group Runge-Kutta method, both of which are straightforward to generalize to higher order. These new simulations recover the diffusion coefficients of stellar eccentricities obtained through analytical calculations of the secular dynamics.

arXiv comment20 pages, 11 figures, submitted to AAS

References in corpus (28)

• Coevolution (Or Not) of Supermassive Black Holes and Host Galaxies
• The Co-Evolution of Galaxies and Supermassive Black Holes: Insights from Surveys of the Contemporary Universe
• REBOUND: An open-source multi-purpose N-body code for collisional dynamics
• Cosmological SPH simulations: The entropy equation
• IAS15: A fast, adaptive, high-order integrator for gravitational dynamics, accurate to machine precision over a billion orbits
• An Update on Monitoring Stellar Orbits in the Galactic Center
• Resonant Relaxation in Stellar Systems
• Satellite dynamics on the Laplace surface
• Stellar Dynamics and Stellar Phenomena Near A Massive Black Hole
• Stellar Dynamics of Extreme-Mass-Ratio Inspirals
• Steady state relativistic stellar dynamics around a massive black hole
• Secular dynamics of hierarchical multiple systems composed of nested binaries, with an arbitrary number of bodies and arbitrary hierarchical structure. First applications to multiplanet and multistar systems
• Gauss's method for secular dynamics, softened
• Secular Stellar Dynamics near a Massive Black Hole
• A numerical study of vector resonant relaxation
• Black Hole Disks in Galactic Nuclei
• Direct N-body simulation of the Galactic centre
• Resonant relaxation near a massive black hole: the dependence on eccentricity
• Stellar Energy Relaxation around A Massive Black Hole
• Scalar Resonant Relaxation of Stars Around a Massive Black Hole
• Relativistic dynamics of stars near a supermassive black hole
• Stellar Dynamics around a Massive Black Hole II: Resonant Relaxation
• Vector Resonant Relaxation of Stars around a Massive Black Hole
• Long-lived lop-sided modes of annular disks orbiting a central mass
• Expansion Techniques for Collisionless Stellar Dynamical Simulations
• Symplectic integrators for spin systems
• The Doubling of Stellar Black Hole Nuclei
• Resonant Relaxation near the Massive Black Hole in the Galactic Center

Cited by in corpus (5)

• Stellar Collisions in Galactic Nuclei: Impact on Destructive Events Near a Supermassive Black Hole
• Dynamical Data Mining Captures Disc-Halo Couplings that Structure Galaxies
• Kinetic theory of one-dimensional inhomogeneous long-range interacting $N$-body systems at order $1/N^2$ without collective effects
• A Lopsided Outer Solar System
• Constraining Axion Dark Matter with Galactic-Centre Resonant Dynamics