An optimum time-stepping scheme for N-body simulations

Abstract

We present a new time-stepping criterion for N-body simulations that is based on the true dynamical time of a particle. This allows us to follow the orbits of particles correctly in all environments since it has better adaptivity than previous time-stepping criteria used in N-body simulations. Furthermore, it requires far fewer force evaluations in low density regions of the simulation and has no dependence on artificial parameters such as, for example, the softening length. This can be orders of magnitude faster than conventional ad-hoc methods that employ combinations of acceleration and softening and is ideally suited for hard problems, such as obtaining the correct dynamics in the very central regions of dark matter haloes. We also derive an eccentricity correction for a general leapfrog integration scheme that can follow gravitational scattering events for orbits with eccentricity e -> 1 with high precision. These new approaches allow us to study a range of problems in collisionless and collisional dynamics from few body problems to cosmological structure formation. We present tests of the time-stepping scheme in N-body simulations of 2-body orbits with eccentricity e -> 1 (elliptic and hyperbolic), equilibrium haloes and a hierarchical cosmological structure formation run.

15 pages, 10 figures, replaced with version that matches published version

References in corpus (5)

• The cosmological simulation code GADGET-2
• The Inner Structure of LambdaCDM Halos I: A Numerical Convergence Study
• Dark Matter Halos within Clusters
• A simple expression for the terms in the Baker-Campbell-Hausdorff series
• A Time-Symmetric Block Time-Step Algorithm for N-Body Simulations

Cited by in corpus (26)

• Clumps and streams in the local dark matter distribution
• Quantifying the heart of darkness with GHALO - a multi-billion particle simulation of our galactic halo
• The AGORA High-Resolution Galaxy Simulations Comparison Project
• The Dark Matter Annihilation Signal from Galactic Substructure: Predictions for GLAST
• N-body simulations of gravitational dynamics
• Determining orbits for the Milky Way's dwarfs
• The Graininess of Dark Matter Haloes
• An Off-center Density Peak in the Milky Way's Dark Matter Halo?
• Cosmic ray lepton puzzle in the light of cosmological N-body simulations
• Cold Dark Matter Substructures in Early-Type Galaxy Halos
• Multimass spherical structure models for N-body simulations
• The Via Lactea INCITE Simulation: Galactic Dark Matter Substructure at High Resolution
• Dynamics of barred galaxies: effects of disk height
• Halo Retention and Evolution of Coalescing Compact Binaries in Cosmological Simulations of Structure Formation: Implications for Short Gamma-Ray Bursts
• The behaviour of shape and velocity anisotropy in dark matter haloes
• Sharpening the dark matter signature in gravitational waveforms II: Numerical simulations with the NbodyIMRI code
• Novel Adaptive softening for collisionless N-body simulations: Eliminating spurious halos
• Perturbations induced by a molecular cloud on the young stellar disc in the Galactic Centre
• Towards time symmetric N-body integration
• A general-purpose timestep criterion for simulations with gravity
• The supermassive black hole merger driven evolution of high-redshift red nuggets into present-day cored early-type galaxies
• The influence of dense gas rings on the dynamics of a stellar disk in the Galactic center
• Evolution of eccentric stellar disks around supermassive black holes: the complex disk disruption dynamics and the milliparsec stars
• Grommet: an N-body code for high-resolution simulations of individual galaxies
• The Structure of Cold Dark Matter Halos: Recent Insights from High Resolution Simulations
• Simulating gravitational motion, gas dynamics, and structure in the cosmos