Quantification of discreteness effects in cosmological N-body simulations: II. Evolution up to shell crossing

Abstract

We apply a recently developed perturbative formalism which describes the evolution under their self-gravity of particles displaced from a perfect lattice to quantify precisely, up to shell crossing, the effects of discreteness in dissipationless cosmological N-body simulations. We give simple expressions, explicitly dependent on the particle density, for the evolution of power in each mode as a function of red-shift. For typical starting red-shifts the effect of finite particle number is to {\it slow down} slightly the growth of power compared to that in the fluid limit (e.g. by about ten percent at half the Nyquist frequency), and to induce also dispersion in the growth as a function of direction at a comparable level. In the limit that the initial red-shift tends to infinity, at fixed particle density, the evolution in fact diverges from that in the fluid limit (described by the Zeldovich approximation). Contrary to widely held belief, this means that a simulation started at a red-shift much higher than the red-shift of shell crossing actually gives a worse, rather than a better, result. We also study how these effects are modified when there is a small-scale regularization of the gravitational force. We show that such a smoothing may reduce the anisotropy of the discreteness effects, but it then {\it increases} their average effect. This behaviour illustrates the fact that the discreteness effects described here are distinct from those usually considered in this context, due to two-body collisions. Indeed the characteristic time for divergence from the collisionless limit is proportional to $N^{2/3}$, rather than $N/\log N$ in the latter case.

18 pages, 10 figures, sequel to astro-ph/0410451 (on initial conditions), final version with minor changes, to appear in Phys. Rev. D

References in corpus (6)

• Discreteness effects in simulations of Hot/Warm dark matter
• Linear perturbative theory of the discrete cosmological N-body problem
• Gravitational Dynamics of an Infinite Shuffled Lattice of Particles
• Gravitational Dynamics of an Infinite Shuffled Lattice: Particle Coarse-grainings, Non-linear Clustering and the Continuum Limit
• Transients from initial conditions based on Lagrangian perturbation theory in N-body simulations
• Particle linear theory on a self-gravitating perturbed cubic Bravais lattice

Cited by in corpus (15)

• Large-scale dark matter simulations
• Accurate estimators of power spectra in N-body simulations
• Evaluating backreaction with the peak model of structure formation
• Discreteness Effects in Lambda Cold Dark Matter Simulations: A Wavelet-Statistical View
• Towards quantitative control on discreteness error in the non-linear regime of cosmological N body simulations
• Perturbation theory with dispersion and higher cumulants: non-linear regime
• 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
• Particle number dependence in the non-linear evolution of N-body self-gravitating systems
• Numerical convergence of pre-initial conditions on dark matter halo properties
• Convergence of halo statistics: code comparison between Rockstar and CompaSO using scale-free simulations
• On the origin of transient features in cosmological N-Body Simulations
• The Hidden Role of Anisotropies in Shaping Structure Formation in Cosmological N-Body Simulations
• Particle loads for cosmological simulations with equal-mass dark matter and baryonic particles
• Tidal adaptive softening and artificial fragmentation in cosmological simulations