Novel Conservative Methods for Adaptive Force Softening in Collisionless and Multi-Species N-Body Simulations
arXiv:2212.06851 · doi:10.1093/mnras/stad2548
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.
arXiv comment25 pages, 13 figures. Updated to match MNRAS accepted version. Six appendices with additional details and tests expanded
References in corpus (89)
- The cosmological simulation code GADGET-2
- E pur si muove: Galiliean-invariant cosmological hydrodynamical simulations on a moving mesh
- The Aquarius Project: the subhalos of galactic halos
- Toward a halo mass function for precision cosmology: the limits of universality
- Galaxies on FIRE (Feedback In Realistic Environments): Stellar Feedback Explains Cosmologically Inefficient Star Formation
- The Rockstar Phase-Space Temporal Halo Finder and the Velocity Offsets of Cluster Cores
- GIZMO: A New Class of Accurate, Mesh-Free Hydrodynamic Simulation Methods
- FIRE-2 Simulations: Physics versus Numerics in Galaxy Formation
- The Inner Structure of LambdaCDM Halos I: A Numerical Convergence Study
- Cosmological SPH simulations: The entropy equation
- Improving convergence in smoothed particle hydrodynamics simulations without pairing instability
- Smoothed Particle Hydrodynamics and Magnetohydrodynamics
- Adaptive Refinement Tree - a new high-resolution N-body code for cosmological simulations
- Fundamental differences between SPH and grid methods
- A General Class of Lagrangian Smoothed Particle Hydrodynamics Methods and Implications for Fluid Mixing Problems
- Modeling Collapse and Accretion in Turbulent Gas Clouds: Implementation and Comparison of Sink Particles in AMR and SPH
- Quantifying the heart of darkness with GHALO - a multi-billion particle simulation of our galactic halo
- 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
- Convergence and scatter of cluster density profiles
- NBODY6++GPU: Ready for the gravitational million-body problem
- Moving mesh cosmology: numerical techniques and global statistics
- A Necessary Condition for Individual Time-Steps in SPH Simulations
- Towards optimal softening in 3D N-body codes: I. Minimizing the force error
- STARFORGE: Toward a comprehensive numerical model of star cluster formation and feedback
- Tracing the Dark Matter Sheet in Phase Space
- Accurate, Meshless Methods for Magneto-Hydrodynamics
- Sweating the small stuff: simulating dwarf galaxies, ultra-faint dwarf galaxies, and their own tiny satellites
- Resolving mixing in Smoothed Particle Hydrodynamics
- 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
- BSE versus StarTrack: implementations of new wind, remnant-formation, and natal-kick schemes in NBODY7 and their astrophysical consequences
- Anisotropic Diffusion in Mesh-Free Numerical Magnetohydrodynamics
- A Constrained-Gradient Method to Control Divergence Errors in Numerical MHD
- The Radical Consequences of Realistic Satellite Orbits for the Heating and Implied Merger Histories of Galactic Disks
- The Abacus Cosmos: A Suite of Cosmological N-body Simulations
- Two--body heating in numerical galaxy formation experiments
- A new approach to simulating collisionless dark matter fluids
- Energy equipartition between stellar and dark matter particles in cosmological simulations results in spurious growth of galaxy sizes
- Numerical convergence of simulations of galaxy formation: the abundance and internal structure of cold dark matter haloes
- Optimal softening for force calculations in collisionless N-body simulations
- An implementation of radiative transfer in the cosmological simulation code GADGET
- The Bar--Halo Interaction--I. From Fundamental Dynamics to Revised N-body Requirements
- Sphenix: Smoothed Particle Hydrodynamics for the next generation of galaxy formation simulations
- SatGen: a semi-analytical satellite galaxy generator -- I. The model and its application to Local-Group satellite statistics
- STARFORGE: The effects of protostellar outflows on the IMF
- Post-Newtonian dynamical modeling of supermassive black holes in galactic-scale simulations
- Optimal Smoothing for N-Body Codes
- 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
- Semi-implicit anisotropic cosmic ray transport on an unstructured moving mesh
- Seeds Don't Sink: Even Massive Black Hole "Seeds" Cannot Migrate to Galaxy Centers Efficiently
- Statistics of Dark Matter Substructure: III. Halo-to-Halo Variance
- rpSPH: a novel Smoothed Particle Hydrodynamics Algorithm
- Stellar mass spectrum within massive collapsing clumps II. Thermodynamics and tidal forces of the first Larson core
- Demonstrating Discreteness and Collision Error in Cosmological N-body Simulations of Dark Matter Gravitational Clustering
- Toward an accurate mass function for precision cosmology
- ColDICE: a parallel Vlasov-Poisson solver using moving adaptive simplicial tessellation
- Symphony: Cosmological Zoom-in Simulation Suites over Four Decades of Host Halo Mass
- Two-Body Relaxation in Cosmological Simulations
- Cluster assembly and the origin of mass segregation in the STARFORGE simulations
- How Biased Are Halo Properties in Cosmological Simulations?
- GANDALF - Graphical Astrophysics code for N-body Dynamics And Lagrangian Fluids
- Direct Integration of the Collisionless Boltzmann Equation in Six-dimensional Phase Space: Self-gravitating Systems
- What is the role of stellar radiative feedback in setting the stellar mass spectrum?
- THE OPTIMAL N-BODY METHOD FOR STABILITY STUDIES OF GALAXIES
- Isothermal Fragmentation: Is there a low-mass cut-off?
- The impact of spurious collisional heating on the morphological evolution of simulated galactic discs
- Gravitational softening as a smoothing operation
- 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
- Adaptive gravitational softening in GADGET
- Tidal Tracks and Artificial Disruption of Cold Dark Matter Halos
- Multidimensional VlasovPoisson Simulations with High-order Monotonicity- and Positivity-preserving Schemes
- Collisionless Stellar Hydrodynamics as an Efficient Alternative to N-body Methods
- Optimal Choice of the Softening Length and Time-Step in N-body Simulations
- Quantifying resolution in cosmological N-body simulations using self-similarity
- BIFROST: simulating compact subsystems in star clusters using a hierarchical fourth-order forward symplectic integrator code
- A New Discrete Dynamical Friction Estimator Based on -body Simulations
- Modelling the accretion and feedback of supermassive black hole binaries in gas-rich galaxy mergers
- A general-purpose timestep criterion for simulations with gravity
- Accuracy and precision of triaxial orbit models I: SMBH mass, stellar mass and dark-matter halo
- The effect of softening on dynamical simulations of galaxies
- Accelerating self-gravitating hydrodynamics simulations with adaptive force updates
- 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