Quantifying resolution in cosmological N-body simulations using self-similarity

Abstract

We demonstrate that testing for self-similarity in scale-free simulations provides an excellent tool to quantify the resolution at small scales of cosmological N-body simulations. Analysing two-point correlation functions measured in simulations using ABACUS, we show how observed deviations from self-similarity reveal the range of time and distance scales in which convergence is obtained. While the well-converged scales show accuracy below 1 percent, our results show that, with a small force softening length, the spatial resolution is essentially determined by the mass resolution. At later times the lower cut-off scale on convergence evolves in comoving units as $a^{-1/2}$ ($a$ being the scale factor), consistent with a hypothesis that it is set by two-body collisionality. A corollary of our results is that N-body simulations, particularly at high red-shift, contain a significant spatial range in which clustering appears converged with respect to the time-stepping and force softening but has not actually converged to the physical continuum result. The method developed can be applied to determine the resolution of any clustering statistic and extended to infer resolution limits for non-scale-free simulations.

15 pages, 9 figures, version published in MNRAS, small changes (more discussion of time-stepping)

References in corpus (11)

• A universal model for halo concentrations
• Disruption of Dark Matter Substructure: Fact or Fiction?
• Dark Matter Substructure in Numerical Simulations: A Tale of Discreteness Noise, Runaway Instabilities, and Artificial Disruption
• Corrfunc --- A Suite of Blazing Fast Correlation Functions on the CPU
• Discreteness Effects in Lambda Cold Dark Matter Simulations: A Wavelet-Statistical View
• The tidal evolution of dark matter substructure -- I. Subhalo density profiles
• Linear perturbative theory of the discrete cosmological N-body problem
• Towards quantitative control on discreteness error in the non-linear regime of cosmological N body simulations
• Einasto Profiles and the Dark Matter Power Spectrum
• Dark Matter Halo Profiles in Scale-Free Cosmologies
• Quantification of discreteness effects in cosmological N-body simulations: II. Evolution up to shell crossing

Cited by in corpus (19)

• AbacusSummit: A Massive Set of High-Accuracy, High-Resolution $N$-Body Simulations
• Large-scale dark matter simulations
• The $\text{Abacus}$ Cosmological $N$-body Code
• The cosmology dependence of galaxy clustering and lensing from a hybrid $N$-body-perturbation theory model
• How Biased Are Halo Properties in Cosmological Simulations?
• Flybys, orbits, splashback: subhalos and the importance of the halo boundary
• Halo Formation from Yukawa Forces in the Very Early Universe
• Perturbation theory with dispersion and higher cumulants: non-linear regime
• Good and Proper: Self-similarity of N-body Simulations with Proper Force Softening
• Simulating Hydrodynamics in Cosmology with CRK-HACC
• Novel Conservative Methods for Adaptive Force Softening in Collisionless and Multi-Species N-Body Simulations
• Accuracy of power spectra in dissipationless cosmological simulations
• Self-Similarity of $k$-Nearest Neighbor Distributions in Scale-Free Simulations
• Accurate halo mass functions from the simplest excursion set theory
• Convergence of halo statistics: code comparison between Rockstar and CompaSO using scale-free simulations
• Cosmological perturbation theory using generalized Einstein de Sitter cosmologies
• Stochastic Super-resolution of Cosmological Simulations with Denoising Diffusion Models
• Testing growth rate dependence in cosmological perturbation theory using scale-free models
• A dynamics-based density profile for dark haloes. I. Algorithm and basic results