Momentum space sampling of neutrinos in $N$-body simulations

Abstract

Including massive neutrinos in $N$-body simulations is a challenging task due to the large thermal velocities of the neutrinos. In particle based codes this leads to problems of shot-noise due to insufficient sampling of the neutrino momentum distribution function. In this paper we investigate the benefits and drawbacks of a scheme first suggested in a paper by Banerjee et al. in which the initial neutrino distribution is symmetrised in momentum space. We confirm that this method reduce shot-noise significantly, but we also find that it generates some spurious power in the neutrino power spectrum at intermediate and small scales. We speculate that this happens because many neutrinos in the simulation sample the same underlying dark matter structures while moving through the simulation. By carefully tuning the number of directions in momentum space of the initial neutrino distribution we show that some improvements can be made over the case where initial neutrino directions are purely random. At redshifts $z\gtrsim 3$ the method works very well, but at smaller redshifts significant improvements are not possible due to the spurious power generation.

14 pages, 8 figures

References in corpus (7)

• The Effect of Thermal Neutrino Motion on the Non-linear Cosmological Matter Power Spectrum
• MassiveNuS: Cosmological Massive Neutrino Simulations
• Simulating nonlinear cosmological structure formation with massive neutrinos
• Relativistic N-body simulations with massive neutrinos
• An Efficient and Accurate Hybrid Method for Simulating Non-Linear Neutrino Structure
• Reducing Noise in Cosmological N-body Simulations with Neutrinos
• $\nu \mathrm{CO}N\mathrm{CEPT}$: Cosmological neutrino simulations from the non-linear Boltzmann hierarchy

Cited by in corpus (8)

• A fast particle-mesh simulation of non-linear cosmological structure formation with massive neutrinos
• Euclid: Modelling massive neutrinos in cosmology -- a code comparison
• Fast simulations of cosmic large-scale structure with massive neutrinos
• The cosmic neutrino background as a collection of fluids in large-scale structure simulations
• Finding closure: approximating Vlasov-Poisson using finitely generated cumulants
• Relativistic Corrections to the Growth of Structure in Modified Gravity
• Hybrid multi-fluid-particle simulations of the cosmic neutrino background
• Improving initialization and evolution accuracy of cosmological neutrino simulations