Modified Hermite Integrators of Arbitrary Order
arXiv:2006.02501 · doi:10.1093/mnras/staa1631
Abstract
We present a family of modified Hermite integrators of arbitrary order possessing superior behaviour for the integration of Keplerian and near-Keplerian orbits. After recounting the derivation of Hermite N-body integrators of arbitrary order, we derive a corrector expression that minimises integrated errors in the argument of periapsis for any such integrator. In addition to providing an alternate derivation of the modified corrector for the 4th-order Hermite integrator, we focus on improved correctors for the 6th- and 8th-order integrators previously featured in the literature. We present a set of numerical examples and find that the higher-order schemes improve performance, even when considering their slightly higher cost in floating point operations. The algorithms presented herein hold promise for systems dominated by central potentials, such as planetary systems and the centres of galaxies. Existing Hermite integrators of any order can be modified to use the expressions presented here with minimal effort. Accordingly the schemes presented herein can be easily implemented on massively parallel architectures.
8 pages, 5 figures, 1 table, accepted for publication in MNRAS
References in corpus (14)
- Stellar disk in the galactic center -- a remnant of a dense accretion disk?
- WASP-47: A Hot Jupiter System with Two Additional Planets Discovered by K2
- The DRAGON simulations: globular cluster evolution with a million stars
- Accelerating NBODY6 with Graphics Processing Units
- GRAPE-6: The massively-parallel special-purpose computer for astrophysical particle simulation
- Monte-Carlo simulations of black hole mergers in AGN disks: Low mergers and predictions for LIGO
- Formation of Terrestrial Planets from Protoplanets under a Realistic Accretion Condition
- The Role of the Kozai-Lidov Mechanism in Black Hole Binary Mergers in Galactic Centers
- 6th and 8th Order Hermite Integrator for N-body Simulations
- RUN DMC: An efficient, parallel code for analyzing Radial Velocity Observations using N-body Integrations and Differential Evolution Markov chain Monte Carlo
- Million-Body Star Cluster Simulations: Comparisons between Monte Carlo and Direct -body
- Star Formation in Accretion Disks and SMBH Growth
- A Time-Symmetric Block Time-Step Algorithm for N-Body Simulations
- Formation of Protoplanets from Massive Planetesimals in Binary Systems