English

Multiple timestep reversible $N$-body integrators for close encounters in planetary systems

Earth and Planetary Astrophysics 2024-04-09 v2 Astrophysics of Galaxies Instrumentation and Methods for Astrophysics Chaotic Dynamics Computational Physics

Abstract

We present new almost time-reversible integrators for solution of planetary systems consisting of "planets" and a dominant mass ("star"). The algorithms can be considered adaptive generalizations of the Wisdom--Holman method, in which all pairs of planets can be assigned timesteps. These timesteps, along with the global timestep, can be adapted time-reversibly, often at no appreciable additional compute cost, without sacrificing any of the long-term error benefits of the Wisdom--Holman method. The method can also be considered a simpler and more flexible version of the \texttt{SYMBA} symplectic code. We perform tests on several challenging problems with close encounters and find the reversible algorithms are up to 2.62.6 times faster than a code based on \texttt{SYMBA}. The codes presented here are available on Github. We also find adapting a global timestep reversibly and discretely must be done in block-synchronized manner or similar.

Keywords

Cite

@article{arxiv.2401.07113,
  title  = {Multiple timestep reversible $N$-body integrators for close encounters in planetary systems},
  author = {David M. Hernandez and Walter Dehnen},
  journal= {arXiv preprint arXiv:2401.07113},
  year   = {2024}
}

Comments

11 pages, 7 Figures. Matches accepted MNRAS version

R2 v1 2026-06-28T14:16:02.648Z