Tidal adaptive softening and artificial fragmentation in cosmological simulations
arXiv:2507.16930 · doi:10.1093/mnras/staf1229
Abstract
Traditional N-body methods introduce localised perturbations in the gravitational forces governing their evolution. These perturbations lead to an artificial fragmentation in the filamentary network of the Large Scale Structure, often referred to as "beads-on-a-string." This issue is particularly apparent in cosmologies with a suppression of the matter power spectrum at small spatial scales, such as warm dark matter models, where the perturbations induced by the N-body discretisation dominate the cosmological power at the suppressed scales. Initial conditions based on third-order Lagrangian perturbation theory, which allow for a late-starting redshift, have been shown to minimise numerical errors contributing to such artefacts. In this work, we investigate whether the additional use of a spatially adaptive softening for dark matter particles, based on the gravitational tidal field, can reduce the severity of artificial fragmentation. Tidal adaptive softening significantly improves force accuracy in idealised filamentary collapse simulations over a fixed softening approach. However, it does not substantially reduce spurious haloes in cosmological simulations when paired with such optimised initial conditions. Nevertheless, tidal adaptive softening induces a shift in halo formation times in warm dark matter simulations compared to a fixed softening counterpart, an effect not seen in cold dark matter simulations. Furthermore, initialising the initial conditions at an earlier redshift generally results in z=0 haloes forming from Lagrangian volumes with lower average sphericity. This sphericity difference could impact post-processing algorithms identifying spurious objects based on Lagrangian volume morphology. We propose potential strategies for reducing spurious haloes without abandoning current N-body methods.
13 pages, 10 figures; accepted for publication in MNRAS
References in corpus (59)
- SciPy 1.0--Fundamental Algorithms for Scientific Computing in Python
- Array Programming with NumPy
- The cosmological simulation code GADGET-2
- Dark Matter Substructure in Galactic Halos
- Where are the missing galactic satellites?
- Baryonic Features in the Matter Transfer Function
- Ultralight scalars as cosmological dark matter
- Small-Scale Challenges to the CDM Paradigm
- GIZMO: A New Class of Accurate, Mesh-Free Hydrodynamic Simulation Methods
- FIRE-2 Simulations: Physics versus Numerics in Galaxy Formation
- Halo Formation in Warm Dark Matter Models
- Dark Matter Self-interactions and Small Scale Structure
- The Inner Structure of LambdaCDM Halos I: A Numerical Convergence Study
- Constraining Warm Dark Matter candidates including sterile neutrinos and light gravitinos with WMAP and the Lyman-alpha forest
- Ahf: Amiga's Halo Finder
- Multi-scale initial conditions for cosmological simulations
- The large-scale structure of the Universe
- Subhaloes in Self-Interacting Galactic Dark Matter Haloes
- The properties of warm dark matter haloes
- The evolution substructure I: a new identification method
- An energy-conserving formalism for adaptive gravitational force softening in SPH and N-body codes
- Milky Way Satellite Census. III. Constraints on Dark Matter Properties from Observations of Milky Way Satellite Galaxies
- Nonlinear Evolution of Cosmological Structures in Warm Dark Matter Models
- Discreteness effects in simulations of Hot/Warm dark matter
- Halo Mass Function and the Free Streaming Scale
- The Warm DM halo mass function below the cut-off scale
- The COpernicus COmplexio: Statistical Properties of Warm Dark Matter Haloes
- Sussing Merger Trees: The Merger Trees Comparison Project
- A new approach to simulating collisionless dark matter fluids
- Accurate initial conditions for cosmological N-body simulations: Minimizing truncation and discreteness errors
- Numerical convergence of simulations of galaxy formation: the abundance and internal structure of cold dark matter haloes
- An adaptively refined phase-space element method for cosmological simulations and collisionless dynamics
- Improving Initial Conditions for Cosmological -Body Simulations
- Discreteness Effects in Lambda Cold Dark Matter Simulations: A Wavelet-Statistical View
- On effects of resolution in dissipationless cosmological simulations
- Discreteness effects in cosmological N-body simulations
- Warm Dark Matter Constraints Using Milky-Way Satellite Observations and Subhalo Evolution Modeling
- Toward an accurate mass function for precision cosmology
- ColDICE: a parallel Vlasov-Poisson solver using moving adaptive simplicial tessellation
- Addressing the too big to fail problem with baryon physics and sterile neutrino dark matter
- Linear perturbative theory of the discrete cosmological N-body problem
- How Biased Are Halo Properties in Cosmological Simulations?
- Quantification of discreteness effects in cosmological N-body simulations: I. Initial Conditions
- Quantification of discreteness effects in cosmological N-body simulations: II. Evolution up to shell crossing
- Forecasts for Galaxy Formation and Dark Matter Constraints from Dwarf Galaxy Surveys
- Simulating the Complexity of the Dark Matter Sheet I: Numerical Algorithms
- The optimal gravitational softening length for cosmological N-body simulations
- Spurious Small-Scale Structure & Discreteness-Driven Relaxation in Cosmological Simulations
- Gravitational Instability in Collisionless Cosmological Pancakes
- Sussing Merger Trees: Stability and Convergence
- Novel Adaptive softening for collisionless N-body simulations: Eliminating spurious halos
- Seeking Observable Imprints of Small-Scale Structure on the Properties of Dark Matter Haloes
- Optimal Softening for N-Body Halo Simulations
- Structural properties of artificial halos in non-standard dark matter simulations
- Degeneracies between baryons and dark matter: the challenge of constraining the nature of dark matter with JWST
- Novel Conservative Methods for Adaptive Force Softening in Collisionless and Multi-Species N-Body Simulations
- Particle linear theory on a self-gravitating perturbed cubic Bravais lattice
- Breaking the baryon-dark matter degeneracy in a model-independent way through the Sunyaev-Zeldovich effect
- Using dust to constrain dark matter models