Novel Conservative Methods for Adaptive Force Softening in Collisionless and Multi-Species N-Body Simulations

Abstract

Modeling self-gravity of collisionless fluids (e.g. ensembles of dark matter, stars, black holes, dust, planetary bodies) in simulations is challenging and requires some force softening. It is often desirable to allow softenings to evolve adaptively, in any high-dynamic range simulation, but this poses unique challenges of consistency, conservation, and accuracy, especially in multi-physics simulations where species with different softening laws may interact. We therefore derive a generalized form of the energy-and-momentum conserving gravitational equations of motion, applicable to arbitrary rules used to determine the force softening, together with consistent associated timestep criteria, interaction terms between species with different softening laws, and arbitrary maximum/minimum softenings. We also derive new methods to maintain better accuracy and conservation when symmetrizing forces between particles. We review and extend previously-discussed adaptive softening schemes based on the local neighbor particle density, and present several new schemes for scaling the softening with properties of the gravitational field, i.e. the potential or acceleration or tidal tensor. We show that the tidal softening scheme not only represents a physically-motivated, translation and Galilean invariant and equivalence-principle respecting (and therefore conservative) method, but imposes negligible timestep or other computational penalties, ensures that pairwise two-body scattering is small compared to smooth background forces, and can resolve outstanding challenges in properly capturing tidal disruption of substructures (minimizing artificial destruction) while also avoiding excessive N-body heating. We make all of this public in the GIZMO code.

25 pages, 13 figures. Updated to match MNRAS accepted version. Six appendices with additional details and tests expanded

References in corpus (27)

• Smoothed Particle Hydrodynamics and Magnetohydrodynamics
• Modeling Collapse and Accretion in Turbulent Gas Clouds: Implementation and Comparison of Sink Particles in AMR and SPH
• Disruption of Dark Matter Substructure: Fact or Fiction?
• Dark Matter Substructure in Numerical Simulations: A Tale of Discreteness Noise, Runaway Instabilities, and Artificial Disruption
• The dynamics and outcome of star formation with jets, radiation, winds, and supernovae in concert
• N-body simulations of gravitational dynamics
• FIRE-3: Updated Stellar Evolution Models, Yields, & Microphysics and Fitting Functions for Applications in Galaxy Simulations
• The Radical Consequences of Realistic Satellite Orbits for the Heating and Implied Merger Histories of Galactic Disks
• Energy equipartition between stellar and dark matter particles in cosmological simulations results in spurious growth of galaxy sizes
• An implementation of radiative transfer in the cosmological simulation code GADGET
• The tidal evolution of dark matter substructure -- II. The impact of artificial disruption on subhalo mass functions and radial profiles
• Spurious heating of stellar motions in simulated galactic disks by dark matter halo particles
• Symphony: Cosmological Zoom-in Simulation Suites over Four Decades of Host Halo Mass
• Cluster assembly and the origin of mass segregation in the STARFORGE simulations
• GANDALF - Graphical Astrophysics code for N-body Dynamics And Lagrangian Fluids
• The impact of spurious collisional heating on the morphological evolution of simulated galactic discs
• Quantification of discreteness effects in cosmological N-body simulations: II. Evolution up to shell crossing
• LCDM halo substructure properties revealed with high resolution and large volume cosmological simulations
• Resolving the Complex Evolution of a Supermassive Black Hole Triplet in a Cosmological Simulation
• Effects of the environment on the multiplicity properties of stars in the STARFORGE simulations
• Tidal Tracks and Artificial Disruption of Cold Dark Matter Halos
• Collisionless Stellar Hydrodynamics as an Efficient Alternative to N-body Methods
• A New Discrete Dynamical Friction Estimator Based on $N$-body Simulations
• BIFROST: simulating compact subsystems in star clusters using a hierarchical fourth-order forward symplectic integrator code
• Modelling the accretion and feedback of supermassive black hole binaries in gas-rich galaxy mergers
• Accuracy and precision of triaxial orbit models I: SMBH mass, stellar mass and dark-matter halo
• Integer Lattice Dynamics for Vlasov-Poisson

Cited by in corpus (6)

• SWIFT: A modern highly-parallel gravity and smoothed particle hydrodynamics solver for astrophysical and cosmological applications
• Numerical challenges for energy conservation in N-body simulations of collapsing self-interacting dark matter halos
• From Seeds to Supermassive Black Holes: Capture, Growth, Migration, and Pairing in Dense Proto-Bulge Environments
• New Cold Dark Matter Crisis Revealed by Multiscale Cluster Lensing
• Early Stages of Dusty Tori: The First Infrared Spectra from a Highly Multiscale Quasar Simulation
• Tidal adaptive softening and artificial fragmentation in cosmological simulations