Discreteness effects in $N$-body simulations of warm dark matter

Abstract

In cosmological $N$-body simulations of warm dark matter, thermal velocities of dark-matter particles are sometimes taken into account by adding random initial velocities to the particles of simulation. However, a particle in the $N$-body system represents a huge collection of dark-matter particles, whose average thermal velocity is very close to zero. We consider justification of the procedure of adding thermal velocities in $N$-body simulations and build a simple model of their influence on the power spectrum. Our model captures the physical effect of suppression of the power spectrum at small wave numbers and also explains its artificial enhancement at large wave numbers, observed in numerical simulations with added thermal velocities. The cause of this enhancement is the disturbance of the growth rate of the density profile introduced when adding random initial thermal velocities. Specifically, the model predicts a turnover in the behavior of the simulated power spectrum at a certain wave number $k_{\ast }$, beyond which it grows as $P(k)\propto k^2$. Our treatment is generalized to a system consisting of several matter components with different thermal velocity dispersion. We also estimate the effects of discreteness related to the bulk velocity field and establish the conditions under which these effects dominate over those of thermal velocities.

arXiv comment34 pages, 5 figures, text improved, one figure added, made comparison with existing numerical simulations, matches the published version

References in corpus (15)

• Sterile Neutrinos as Dark Matter
• Constraining Warm Dark Matter candidates including sterile neutrinos and light gravitinos with WMAP and the Lyman-alpha forest
• Massive sterile neutrinos as warm Dark Matter
• Lyman-alpha constraints on warm and on warm-plus-cold dark matter models
• Sterile Neutrino Dark Matter
• Discreteness effects in simulations of Hot/Warm dark matter
• The effect of neutrinos on the matter distribution as probed by the Intergalactic Medium
• Structure Formation with Cold + Hot Dark Matter
• The Effect of Thermal Neutrino Motion on the Non-linear Cosmological Matter Power Spectrum
• On the Structure of Dark Matter Halos at the Damping Scale of the Power Spectrum with and without Relict Velocities
• An adaptively refined phase-space element method for cosmological simulations and collisionless dynamics
• Discreteness Effects in Lambda Cold Dark Matter Simulations: A Wavelet-Statistical View
• ColDICE: a parallel Vlasov-Poisson solver using moving adaptive simplicial tessellation
• Simulating the Complexity of the Dark Matter Sheet I: Numerical Algorithms
• The Effect of Thermal Velocities on Structure Formation in N-body Simulations of Warm Dark Matter

Cited by in corpus (1)

• The CosmoVerse White Paper: Addressing observational tensions in cosmology with systematics and fundamental physics