TRACE: a code for Time-Reversible Astrophysical Close Encounters

Abstract

We present TRACE, an almost time-reversible hybrid integrator for the planetary N-body problem. Like hybrid symplectic integrators, TRACE can resolve close encounters between particles while retaining many of the accuracy and speed advantages of a fixed time-step symplectic method such the Wisdom-Holman map. TRACE switches methods time-reversibly during close encounters following the prescription of Hernandez & Dehnen. In this paper we describe the derivation and implementation of TRACE and study its performance for a variety of astrophysical systems. In all our test cases, TRACE is at least as accurate and fast as the hybrid symplectic integrator MERCURIUS. In many cases, TRACE's performance is vastly superior to that of MERCURIUS. In test cases with planet-planet close encounters, TRACE is as accurate as MECURIUS with a 12x speed-up. If close encounters with the central star are considered, TRACE achieves good error performance while MERCURIUS fails to give qualitatively correct results. In ensemble tests of violent scattering systems, TRACE matches the high-accuracy IAS15 while providing a 15x speed-up. In large N systems simulating lunar accretion, TRACE qualitatively gives the same results as IAS15 but at a 41x speed-up. We also discuss some cases such as von Zeipel-Lidov-Kozai cycles where hybrid integrators perform poorly and provide some guidance on which integrator to use for which system. TRACE is freely available within the REBOUND package.

Published in MNRAS. Minor changes

References in corpus (27)

• REBOUND: An open-source multi-purpose N-body code for collisional dynamics
• La2010: A new orbital solution for the long term motion of the Earth
• Dynamical Outcomes of Planet-Planet Scattering
• The Eccentric Kozai-Lidov Effect and Its Applications
• Architecture and Dynamics of Kepler's Candidate Multiple Transiting Planet Systems
• IAS15: A fast, adaptive, high-order integrator for gravitational dynamics, accurate to machine precision over a billion orbits
• Origins of Eccentric Extrasolar Planets: Testing the Planet-Planet Scattering Model
• WHFast: A fast and unbiased implementation of a symplectic Wisdom-Holman integrator for long term gravitational simulations
• High-resolution simulations of the final assembly of Earth-like planets 1: terrestrial accretion and dynamics
• REBOUNDx: A Library for Adding Conservative and Dissipative Forces to Otherwise Symplectic N-body Integrations
• Hybrid Symplectic Integrators for Planetary Dynamics
• Halo and Galaxy Formation Histories from the Millennium Simulation: Public release of a VO-oriented and SQL-queryable database for studying the evolution of galaxies in the LambdaCDM cosmogony
• High order symplectic integrators for planetary dynamics and their implementation in REBOUND
• Self-Consistent Spin, Tidal and Dynamical Equations of Motion in the REBOUNDx Framework
• A study of symplectic integrators for planetary system problems: error analysis and comparisons
• Symplectic fourth-order maps for the collisional N-body problem
• Symplectic integration for the collisional gravitational $N$-body problem
• Fast and reliable symplectic integration for planetary system $N$-body problems
• Time-symmetric integration in astrophysics
• Towards time symmetric N-body integration
• Stepsize errors in the $N$-body problem: discerning Mercury's true possible long-term orbits
• A new timestep criterion for N-body simulations
• Switching integrators reversibly in the astrophysical $N$-body problem
• Reversible time-step adaptation for the integration of few-body systems
• Improving the accuracy of simulated chaotic $N$-body orbits using smoothness
• WHFast512: A symplectic N-body integrator for planetary systems optimized with AVX512 instructions
• Planet-Planet Scattering and ZLK Migration -- The Dynamical History of HAT-P-11

Cited by in corpus (7)

• Breaking long-period resonance chains with stellar flybys
• Self-consistent Dynamical and Chaotic Tides in the REBOUNDx framework
• Measuring the Mean Plane of the Distant Kuiper Belt
• The Perturbation Theory Approach to Stability in the Scattered Disk
• Close-in compact super-Earth systems emerging from resonant chains: slow destabilization by unseen remnants of formation
• Origins of Mercury's Big Heart of Iron: Exploring Pathways to Form High Core Mass Fraction (CMF) Planets via N-body Simulations
• A Dormant Captured Oort Cloud Comet Awakens: (18916) 2000 OG44