Accuracy of power spectra in dissipationless cosmological simulations

Abstract

We exploit a suite of large \emph{N}-body simulations (up to N=$4096^3$) performed with \Abacus, of scale-free models with a range of spectral indices $n$, to better understand and quantify convergence of the matter power spectrum. Using self-similarity to identify converged regions, we show that the maximal wavenumber resolved at a given level of accuracy increases monotonically as a function of time. At the 1\% level it starts at early times from a fraction of $k_{\Lambda }$, the Nyquist wavenumber of the initial grid, and reaches at most, if the force softening is sufficiently small, $\sim 2-3k_{\Lambda }$ at the very latest times we evolve to. At the $5\%$ level, accuracy extends up to wavenumbers of order $5k_{\Lambda }$ at late times. Expressed as a suitable function of the scale-factor, accuracy shows a very simple $n$-dependence, allowing a extrapolation to place conservative bounds on the accuracy of \emph{N}-body simulations of non-scale free models like LCDM. We note that deviations due to discretization in the converged range are not well modelled by shot noise, and subtracting it in fact degrades accuracy. Quantitatively our findings are broadly in line with the conservative assumptions about resolution adopted by recent studies using large cosmological simulations (e.g. Euclid Flagship) aiming to constrain the mildly non-linear regime. On the other hand, we remark that conclusions about small scale clustering (e.g. concerning the validity of stable clustering) obtained using PS data at wavenumbers larger than a few $k_{\Lambda }$ may need revision in light of our convergence analysis.

arXiv comment14 pages, 10 figures, version to match published version. Additional subsection in Shot noise and discreteness effects

References in corpus (38)

• Planck 2013 results. XVI. Cosmological parameters
• LSST: from Science Drivers to Reference Design and Anticipated Data Products
• Stable clustering, the halo model and nonlinear cosmological power spectra
• First results from the IllustrisTNG simulations: matter and galaxy clustering
• The DESI Experiment Part I: Science,Targeting, and Survey Design
• Revising the Halofit Model for the Nonlinear Matter Power Spectrum
• Evolution of structure in cold dark matter universes
• HMcode-2020: Improved modelling of non-linear cosmological power spectra with baryonic feedback
• The Coyote Universe I: Precision Determination of the Nonlinear Matter Power Spectrum
• The Uchuu Simulations: Data Release 1 and Dark Matter Halo Concentrations
• Correcting for the alias effect when measuring the power spectrum using FFT
• AbacusSummit: A Massive Set of High-Accuracy, High-Resolution $N$-Body Simulations
• The BACCO Simulation Project: Exploiting the full power of large-scale structure for cosmology
• Euclid preparation: II. The EuclidEmulator -- A tool to compute the cosmology dependence of the nonlinear matter power spectrum
• Calibrating the Nonlinear Matter Power Spectrum: Requirements for Future Weak Lensing Surveys
• The detectability of baryonic acoustic oscillations in future galaxy surveys
• Matter power spectrum and the challenge of percent accuracy
• The $\text{Abacus}$ Cosmological $N$-body Code
• A High-Fidelity Realization of the Euclid Code Comparison $N$-body Simulation with Abacus
• Fundamental Discreteness Limitations of Cosmological N-Body Clustering Simulations
• Self-similarity and scaling behavior of scale-free gravitational clustering
• Improving Initial Conditions for Cosmological $N$-Body Simulations
• Discreteness Effects in Lambda Cold Dark Matter Simulations: A Wavelet-Statistical View
• General Requirements on Matter Power Spectrum Predictions for Cosmology with Weak Lensing Tomography
• On effects of resolution in dissipationless cosmological simulations
• Self-Similar Evolution of Gravitational Clustering: N-Body Simulations of the n=-2 Spectrum
• Linear perturbative theory of the discrete cosmological N-body problem
• Nonlinear Biasing and Redshift-Space Distortions in Lagrangian Resummation Theory and N-body Simulations
• Towards quantitative control on discreteness error in the non-linear regime of cosmological N body simulations
• Gravitational evolution of a perturbed lattice and its fluid limit
• Power spectrum for the small-scale Universe
• Quantifying resolution in cosmological N-body simulations using self-similarity
• Scaling laws in gravitational clustering for counts-in-cells and mass functions
• The Last Journey. I. An Extreme-Scale Simulation on the Mira Supercomputer
• Testing dark matter halo properties using self-similarity
• On Mitigation of the Uncertainty in Nonlinear Matter Clustering for Cosmic Shear Tomography
• Good and Proper: Self-similarity of N-body Simulations with Proper Force Softening
• Stable clustering and the resolution of dissipationless cosmological N-body simulations

Cited by in corpus (5)

• The MillenniumTNG Project: High-precision predictions for matter clustering and halo statistics
• Convergence of halo statistics: code comparison between Rockstar and CompaSO using scale-free simulations
• Cosmological perturbation theory using generalized Einstein de Sitter cosmologies
• A dynamics-based density profile for dark haloes. I. Algorithm and basic results
• Testing growth rate dependence in cosmological perturbation theory using scale-free models