Secular dynamics around a supermassive black hole via multipole expansion

• Verfasst von: Fouvry, Jean-Baptiste [VerfasserIn]
• Verfasst von: Dehnen, Walter [VerfasserIn]
• Verfasst von: Tremaine, Scott [VerfasserIn]
• Verfasst von: Bar-Or, Ben [VerfasserIn]
• Titel: Secular dynamics around a supermassive black hole via multipole expansion
• Verf.angabe: Jean-Baptiste Fouvry, Walter Dehnen, Scott Tremaine, and Ben Bar-Or
• E-Jahr: 2022
• Jahr: 2022 May 18
• Umfang: 17 S.
• Fussnoten: Gesehen am 28.06.2022
• Titel Quelle: Enthalten in: The astrophysical journal / 1
• Ort Quelle: London : Institute of Physics Publ., 1995
• Jahr Quelle: 2022
• Band/Heft Quelle: 931(2022), 1, Artikel-ID 8, Seite 1-17
• ISSN Quelle: 1538-4357
• DOI: doi:10.3847/1538-4357/ac602e
• Datenträger: Online-Ressource
• Sprache: eng
• K10plus-PPN: 1808027116

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 Gauß, is to use the secular approximation of smearing out 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 reformulate the discretization 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 , reducing the overall computational cost from O(N 2 K 2) to . 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.