Optimal Choice of the Softening Length and Time-Step in N-body Simulations

Abstract

A criterion for the choice of optimal softening length $\epsilon$ and time-step $dt$ for $N$-body simulations of a collisionless stellar system is analyzed. Plummer and Hernquist spheres are used as models to follow how changes in various parameters of an initially equilibrium stable model depend on $\epsilon$ and $dt$. These dependences are used to derive a criterion for choosing $\epsilon$ and $dt$. The resulting criterion is compared to Merritt's criterion for choosing the softening length, which is based on minimizing the mean irregular force acting on a particle with unit mass. Our criterion for choosing $\epsilon$ and $dt$ indicate that $\epsilon$ must be a factor of 1.5-2 smaller than the mean distance between particles in the densest regions to be resolved. The time-step must always be adjusted to the chosen $\epsilon$ (the particle must, on average, travel a distance smaller than $0.5\epsilon$ during one time-step). An algorithm for solving N-body problems with adaptive variations of the softening length is discussed in connection with the task of choosing $\epsilon$, but is found not to be promising.

arXiv commentTo be published in Astonomy Reports, 11 pages, 5 figures

References in corpus (4)

• Towards optimal softening in 3D N-body codes: I. Minimizing the force error
• Optimal Smoothing for N-Body Codes
• Performance and accuracy of a GRAPE-3 system for collisionless N-body simulations
• Relaxation times calculated from angular deflections

Cited by in corpus (18)

• FIRE-2 Simulations: Physics versus Numerics in Galaxy Formation
• An energy-conserving formalism for adaptive gravitational force softening in SPH and N-body codes
• The fine-grained phase-space structure of Cold Dark Matter halos
• An Iterative Method for Constructing Equilibrium Phase Models of Stellar Systems
• A parallel-GPU code for asteroid aggregation problems with angular particles
• A New Discrete Dynamical Friction Estimator Based on $N$-body Simulations
• Sharpening the dark matter signature in gravitational waveforms II: Numerical simulations with the NbodyIMRI code
• GalactICS with Gas
• Novel Conservative Methods for Adaptive Force Softening in Collisionless and Multi-Species N-Body Simulations
• Phase models of the Milky Way stellar disc
• Discreteness effects, $N-$body chaos and the onset of radial-orbit instability
• Resolution criteria to avoid artificial clumping in Lagrangian hydrodynamic simulations with a multi-phase interstellar medium
• The effect of softening on dynamical simulations of galaxies
• Forming disc galaxies in major mergers II. The central mass concentration problem and a comparison of GADGET3 with GIZMO
• Metal distribution in the ICM - a comprehensive numerical study of twelve galaxy clusters
• Long-Time Evolution of Gas-Free Disc Galaxies in Binary Systems
• Numerical Study of Bar Suppression in Galaxy Models Due to Disc Heating
• A Simulation of the Dependence of Tidal Interaction on Galaxy Type in Compact Groups