On the origin of transient features in cosmological N-Body Simulations

Abstract

We study the effect of gravitational clustering at small scales on larger scales by studying mode coupling between virialised halos. We build on the calculation by Peebles (1974) where it was shown that a virialised halo does not contribute any mode coupling terms at small wave numbers $k$. Using a perturbative expansion in wave number, we show that this effect is small and arises from the deviation of halo shapes from spherical and also on tidal interactions between halos. We connect this with the impact of finite mass resolution of cosmological N-Body simulations on the evolution of perturbations at early times. This difference between the expected evolution and the evolution obtained in cosmological N-Body simulations can be quantified using such an estimate. We also explore the impact of a finite shortest scale up to which the desired power spectrum is realised in simulations. Several simulation studies have shown that this effect is small in comparison with the effect of perturbations at large scales on smaller scales. It is nevertheless important to study these effects and develop a general approach for estimating their magnitude. This is especially relevant in the present era of precision cosmology. We provide basic estimates of the magnitude of these effects and their power spectrum dependence. We find that the impact of small scale cutoff in the initial power spectrum and discreteness increases with $(n+3)$, with $n$ being the index of the power spectrum. In general, we recommend that cosmological simulation data should be used only if the scale of non-linearity, defined as the scale where the linearly extrapolated {\it rms} amplitude of fluctuations is unity, is larger than the average inter-particle separation.

13 pages, 5 figures

References in corpus (21)

• Simulating cosmic structure formation with the GADGET-4 code
• Large-scale dark matter simulations
• Accurate initial conditions for cosmological N-body simulations: Minimizing truncation and discreteness errors
• Discreteness Effects in Lambda Cold Dark Matter Simulations: A Wavelet-Statistical View
• Linear perturbative theory of the discrete cosmological N-body problem
• How Biased Are Halo Properties in Cosmological Simulations?
• Comments on the size of the simulation box in cosmological N-Body simulations
• Simulations of Baryon Acoustic Oscillations I: Growth of Large-Scale Density Fluctuations
• Quantification of discreteness effects in cosmological N-body simulations: II. Evolution up to shell crossing
• Gravitational Dynamics of an Infinite Shuffled Lattice of Particles
• Effects of the size of cosmological N-Body simulations on physical quantities - III: Skewness
• On the Onset of Stochasticity in $\Lambda$CDM Cosmological Simulations
• Perturbation theory with dispersion and higher cumulants: framework and linear theory
• Gravitational Dynamics of an Infinite Shuffled Lattice: Particle Coarse-grainings, Non-linear Clustering and the Continuum Limit
• Force distribution in a randomly perturbed lattice of identical particles with $1/r^2$ pair interaction
• Mock halo catalogs: assigning unresolved halo properties using correlations with local halo environment
• Gravitational collapse in an expanding background and the role of substructure I: Planar collapse
• Renormalization group and UV completion of cosmological perturbations: Gravitational collapse as a critical phenomenon
• Numerical convergence of pre-initial conditions on dark matter halo properties
• Halo mass function in scale invariant models
• Kinetic Field Theory: Perturbation theory beyond first order

Cited by in corpus (2)

• The Effects of Linear Matter Power Spectrum Enhancement on Dark Matter Substructure
• Particle loads for cosmological simulations with equal-mass dark matter and baryonic particles