Optimal softening for force calculations in collisionless N-body simulations

Abstract

In N-body simulations the force calculated between particles representing a given mass distribution is usually softened, to diminish the effect of graininess. In this paper we study the effect of such a smoothing, with the aim of finding an optimal value of the softening parameter. As already shown by Merritt (1996), for too small a softening the estimates of the forces will be too noisy, while for too large a softening the force estimates are systematically misrepresented. In between there is an optimal softening, for which the forces in the configuration approach best the true forces. The value of this optimal softening depends both on the mass distribution and on the number of particles used to represent it. For higher number of particles the optimal softening is smaller. More concentrated mass distributions necessitate smaller softening, but the softened forces are never as good an approximation of the true forces as for not centrally concentrated configurations. We give good estimates of the optimal softening for homogeneous spheres, Plummer spheres, and Dehnen spheres. We also give a rough estimate of this quantity for other mass distributions, based on the harmonic mean distance to the $k$th neighbour ($k$ = 1, .., 12), the mean being taken over all particles in the configuration. Comparing homogeneous Ferrers ellipsoids of different shapes we show that the axial ratios do not influence the value of the optimal softening. Finally we compare two different types of softening, a spline softening (Hernquist & Katz 1989) and a generalisation of the standard Plummer softening to higher values of the exponent. We find that the spline softening fares roughly as well as the higher powers of the power-law softening and both give a better representation of the forces than the standard Plummer softening.

16 pages Latex, 19 figures, accepted for publication in MNRAS, corrected typos, minor changes mainly in sec. 5

Cited by in corpus (63)

• GADGET: A code for collisionless and gasdynamical cosmological simulations
• FIRE-2 Simulations: Physics versus Numerics in Galaxy Formation
• Phantom: A smoothed particle hydrodynamics and magnetohydrodynamics code for astrophysics
• An energy-conserving formalism for adaptive gravitational force softening in SPH and N-body codes
• Disruption of Dark Matter Substructure: Fact or Fiction?
• Dark Matter Substructure in Numerical Simulations: A Tale of Discreteness Noise, Runaway Instabilities, and Artificial Disruption
• A Hierarchical O(N) Force Calculation Algorithm
• Towards optimal softening in 3D N-body codes: I. Minimizing the force error
• Initial conditions to cosmological N-body simulations, or how to run an ensemble of simulations
• And yet it moves: The dangers of artificially fixing the Milky Way center of mass in the presence of a massive Large Magellanic Cloud
• Can bars be destroyed by a central mass concentration? I. Simulations
• Galaxy merging, the Fundamental Plane of elliptical galaxies, and the Mbh-sigma relation
• Do dwarf galaxies form in tidal tails?
• Numerical requirements for simulations of self gravitating and non-self gravitating disks
• The fine-grained phase-space structure of Cold Dark Matter halos
• Dynamical Friction on Star Clusters near the Galactic Center
• Formation and evolution of galactic disks with a multiphase numerical model
• Density profiles of ultracompact minihalos: Implications for constraining the primordial power spectrum
• An Iterative Method for Constructing Equilibrium Phase Models of Stellar Systems
• VINE -- A numerical code for simulating astrophysical systems using particles I: Description of the physics and the numerical methods
• Tidal breakup of binary stars at the Galactic Center. II. Hydrodynamic simulations
• A new code for orbit analysis and Schwarzschild modelling of triaxial stellar systems
• Minimum Velocity Dispersion in Stable Stellar Disks. Numerical Simulations
• Radial orbital anisotropy and the Fundamental Plane of elliptical galaxies
• The Vertical Stellar Kinematics in Face-On Barred Galaxies: Estimating the Ages of Bars
• The Effect of Many Minor Mergers on the Size Growth of Compact Quiescent Galaxies
• The Effects of Boxy/Peanut Bulges on Galaxy Models
• Gravitational softening as a smoothing operation
• EvoL: The new Padova T-SPH parallel code for cosmological simulations - I. Basic code: gravity and hydrodynamics
• On the isolated dwarf galaxies: from cuspy to flat dark matter density profiles and metallicity gradients
• Adaptive gravitational softening in GADGET
• On the equilibrium morphology of systems drawn from spherical collapse experiments
• The optimal gravitational softening length for cosmological N-body simulations
• Twin Binaries: Studies of Stability, Mass Transfer, and Coalescence
• Expansion Techniques for Collisionless Stellar Dynamical Simulations
• The Effect of Tidal Stripping on Composite Stellar Populations in Dwarf Spheroidal Galaxies
• A New Discrete Dynamical Friction Estimator Based on $N$-body Simulations
• Figure rotation of dark halos in CDM simulations
• Sharpening the dark matter signature in gravitational waveforms II: Numerical simulations with the NbodyIMRI code
• Dark influences II: gas and star formation in minor mergers of dwarf galaxies with dark satellites
• Dark influences: imprints of dark satellites on dwarf galaxies
• The Adaptive TreePM: An Adaptive Resolution Code for Cosmological N-body Simulations
• The influence of numerical parameters on tidally triggered bar formation
• Optimal Softening for N-Body Halo Simulations
• Effect of Population III Multiplicity on Dark Star Formation
• Description of the evolution of inhomogeneities on a Dark Matter halo with the Vlasov equation
• Regular and chaotic motion in softened gravitational systems
• A Natural Symmetrization for the Plummer Potential
• Novel Conservative Methods for Adaptive Force Softening in Collisionless and Multi-Species N-Body Simulations
• Relaxation times calculated from angular deflections
• Characterization and dynamics of the peculiar stream Jhelum, A tentative role for the Sagittarius dwarf galaxy
• Performance analysis of parallel gravitational $N$-body codes on large GPU cluster
• The effect of softening on dynamical simulations of galaxies
• Resolution criteria to avoid artificial clumping in Lagrangian hydrodynamic simulations with a multi-phase interstellar medium
• Premature supermassive black hole mergers in cosmological simulations of structure formation
• Collapsed and Extended Cold Dark Matter Haloes in Softened N-Body Gravity
• Long-Time Evolution of Gas-Free Disc Galaxies in Binary Systems
• Comparing Extended System Interactions with Motions in Softened Potentials
• The status of isochrony in the formation and evolution of self-gravitating systems
• A Simulation of the Dependence of Tidal Interaction on Galaxy Type in Compact Groups
• Numerical Study of Bar Suppression in Galaxy Models Due to Disc Heating
• The efficiency of the spiral-in of a black hole to the Galactic centre
• Gravitational Waves from Hyperbolic Encounters of Primordial Black Holes in Dwarf Galaxies