N-body Integrators with Individual Time Steps from Hierarchical Splitting

Abstract

We review the implementation of individual particle time-stepping for N-body dynamics. We present a class of integrators derived from second order Hamiltonian splitting. In contrast to the usual implementation of individual time-stepping, these integrators are momentum conserving and show excellent energy conservation in conjunction with a symmetrized time step criterion. We use an explicit but approximate formula for the time symmetrization that is compatible with the use of individual time steps. No iterative scheme is necessary. We implement these ideas in the HUAYNO (available online at www.amusecode.org) code and present tests of the integrators and show that the presented integration schemes shows good energy conservation, with little or no systematic drift, while conserving momentum and angular momentum to machine precision for long term integrations.

20 pages, 9 figures, accepted for publication in New Astronomy

References in corpus (7)

• A multiphysics and multiscale software environment for modeling astrophysical systems
• N-body simulations of gravitational dynamics
• 6th and 8th Order Hermite Integrator for N-body Simulations
• A Hybrid N-Body Code Incorporating Algorithmic Regularization and Post-Newtonian Forces
• A Time-Symmetric Block Time-Step Algorithm for N-Body Simulations
• MYRIAD: A new N-body code for simulations of Star Clusters
• Grommet: an N-body code for high-resolution simulations of individual galaxies

Cited by in corpus (61)

• Simulating Galaxy Formation with the IllustrisTNG Model
• The Arepo public code release
• Simulating cosmic structure formation with the GADGET-4 code
• Universal structure of dark matter haloes over a mass range of 20 orders of magnitude
• The Astrophysical Multipurpose Software Environment
• The impact of companions on stellar evolution
• 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
• The formation of planets in circumbinary disks
• How Sedna and family were captured in a close encounter with a solar sibling
• Survivability of planetary systems in young and dense star clusters
• Was the nineteenth century giant eruption of Eta Carinae a merger event in a triple system?
• Slowly, slowly in the wind: 3D hydrodynamical simulations of wind mass transfer and angular-momentum loss in AGB binary systems
• Gone with the wind: the impact of wind mass transfer on the orbital evolution of AGB binary systems
• Simulations of common envelope evolution in triple systems: Circumstellar case
• Finding the imprints of stellar encounters in long period comets
• Common envelope evolution of eccentric binaries
• The evolution of triples with a Roche-lobe filling outer star
• Planet-Mediated Precision-Reconstruction of the Evolution of the Cataclysmic Variable HU Aquarius
• Mass transfer between debris discs during close stellar encounters
• The growth of intermediate mass black holes through tidal captures and tidal disruption events
• The evolution of the Sun's birth cluster and the search for the solar siblings with Gaia
• Oort cloud Ecology II: Extra-solar Oort clouds and the origin of asteroidal interlopers
• Capture of exocomets and the erosion of the Oort cloud due to stellar encounters in the Galaxy
• Oort cloud Ecology II: The chronology of the formation of the Oort cloud
• The eccentric behaviour of windy binary stars
• FROST: a momentum-conserving CUDA implementation of a hierarchical fourth-order forward symplectic integrator
• Exploring the role of the Sun's motion in terrestrial comet impacts
• A debris disk under the influence of a wide planetary mass companion: The system of HD106906
• BIFROST: simulating compact subsystems in star clusters using a hierarchical fourth-order forward symplectic integrator code
• The secular evolution of the Kuiper belt after a close stellar encounter
• Mergers of Equal-Mass Binaries with Compact Object Companions from Mass Transfer in Triple Star Systems
• Time-symmetric integration in astrophysics
• A Triple Origin for Twin Blue Stragglers in Close Binaries
• Quest for finding the lost siblings of the Sun
• Towards time symmetric N-body integration
• Non-intrusive hierarchical coupling strategies for multi-scale simulations in gravitational dynamics
• A Keplerian-based Hamiltonian Splitting for Gravitational $N$-body Simulations
• The Origin of Chaos in the Orbit of Comet 1P/Halley
• Changes in orientation and shape of protoplanetary discs moving through an ambient medium
• Disc truncation in embedded star clusters: Dynamical encounters versus face-on accretion
• A connected component-based method for efficiently integrating multiscale $N$-body systems
• Examining the effects of dark matter spikes on eccentric intermediate mass ratio inspirals using N-body simulations
• Dynamical Ejections of Stars due to an Accelerating Gas Filament
• Simulating Hydrodynamics in Cosmology with CRK-HACC
• Fast Multipole Methods for $N$-body Simulations of Collisional Star Systems
• Simulating stellar winds in AMUSE
• Switching integrators reversibly in the astrophysical $N$-body problem
• Evolution of Massive Black Hole Binaries in Collisionally Relaxed Nuclear Star Clusters -- Impact of Mass Segregation
• Reversible time-step adaptation for the integration of few-body systems
• Hybrid multi-fluid-particle simulations of the cosmic neutrino background
• SiMon: Simulation Monitor for Computational Astrophysics
• A direct N-body integrator for modelling the chaotic, tidal dynamics of multi-body extrasolar systems: TIDYMESS
• Expanding shells around young clusters -- S 171/Be 59
• Oort Cloud Ecology. III. The Sun left the parent star cluster shortly after the giant planets formed
• The Signature of a Windy Radio Supernova Progenitor in a Binary System
• The 800pc long tidal tails of the Hyades star cluster: Possible discovery of candidate epicyclic overdensities from an open star cluster
• Venice: a multi-scale operator-splitting algorithm for multi-physics simulations
• Close Encounters of the Interstellar Kind: Examining the Capture of Interstellar Objects in Near Earth Orbit
• Tracing Milky Way substructure with an RR Lyrae hierarchical clustering forest
• The solar siblings in the Gaia era
• Imladris: a detailed and flexible model for galaxy simulations with individual stars