Geodesic motion and phase-space evolution of massive neutrinos

Abstract

The non-trivial phase-space distribution of relic neutrinos is responsible for the erasure of primordial density perturbations on small scales, which is one of the main cosmological signatures of neutrino mass. In this paper, we present a new code, FastDF, for generating 1%-accurate particle realisations of the neutrino phase-space distribution using relativistic perturbation theory. We use the geodesic equation to derive equations of motion for massive particles moving in a weakly perturbed spacetime and integrate particles accordingly. We demonstrate how to combine geodesic-based initial conditions with the $\delta f$ method to minimise shot noise and clarify the definition of the neutrino momentum, finding that large errors result if the wrong parametrisation is used. Compared to standard Lagrangian methods with ad-hoc thermal motions, FastDF achieves substantial improvements in accuracy. We outline the approximation schemes used to speed up the code and to ensure symplectic integration that preserves phase-space density. Finally, we discuss implications for neutrino particles in cosmological N-body simulations. In particular, we argue that particle methods can accurately describe the neutrino distribution from $z=10^9$, when neutrinos are linear and ultra-relativistic, down to $z=0$, when they are nonlinear and non-relativistic. FastDF can be used to set up accurate initial conditions (ICs) for N-body simulations and has been integrated into the higher-order IC code monofonIC.

arXiv comment17 pages, 6 figures, contains link to code

References in corpus (47)

• Cosmological Perturbation Theory in the Synchronous vs. Conformal Newtonian Gauge
• Cosmological Perturbation Theory in the Synchronous and Conformal Newtonian Gauges
• Dark Energy Survey Year 3 Results: Cosmological Constraints from Galaxy Clustering and Weak Lensing
• Measurement of Day and Night Neutrino Energy Spectra at SNO and Constraints on Neutrino Mixing Parameters
• The Cosmic Linear Anisotropy Solving System (CLASS) I: Overview
• Weighing Neutrinos with Galaxy Surveys
• The Cosmic Linear Anisotropy Solving System (CLASS) IV: Efficient implementation of non-cold relics
• Measurement of the solar neutrino energy spectrum using neutrino-electron scattering
• The Quijote simulations
• The effect of neutrinos on the matter distribution as probed by the Intergalactic Medium
• Cosmological $N$-body simulations with suppressed variance
• DEMNUni: The clustering of large-scale structures in the presence of massive neutrinos
• Structure Formation with Cold + Hot Dark Matter
• An efficient implementation of massive neutrinos in non-linear structure formation simulations
• Updated results on neutrino mass and mass hierarchy from cosmology with Planck 2018 likelihoods
• The Effect of Thermal Neutrino Motion on the Non-linear Cosmological Matter Power Spectrum
• On the most constraining cosmological neutrino mass bounds
• Simulating nonlinear cosmological structure formation with massive neutrinos
• Grid Based Linear Neutrino Perturbations in Cosmological N-body Simulations
• Accurate initial conditions for cosmological N-body simulations: Minimizing truncation and discreteness errors
• Relativistic N-body simulations with massive neutrinos
• General relativistic corrections to $N$-body simulations and the Zel'dovich approximation
• An Efficient and Accurate Hybrid Method for Simulating Non-Linear Neutrino Structure
• An optimal nonlinear method for simulating relic neutrinos
• Cosmological neutrino simulations at extreme scale
• Precision reconstruction of the dark matter-neutrino relative velocity from N-body simulations
• Reducing Noise in Cosmological N-body Simulations with Neutrinos
• $\nu \mathrm{CO}N\mathrm{CEPT}$: Cosmological neutrino simulations from the non-linear Boltzmann hierarchy
• COLA with massive neutrinos
• How to add massive neutrinos to your $\Lambda$CDM simulation -- extending cosmology rescaling algorithms
• A Lagrangian Perturbation Theory in the presence of massive neutrinos
• Fully relativistic treatment of light neutrinos in $N$-body simulations
• General relativistic weak-field limit and Newtonian N-body simulations
• A Calculation of the Full Neutrino Phase Space in Cold+Hot Dark Matter Models
• Higher-order initial conditions with massive neutrinos
• A new approach to cosmological structure formation with massive neutrinos
• One line to run them all: SuperEasy massive neutrino linear response in $N$-body simulations
• Extending Lagrangian perturbation theory to a fluid with velocity dispersion
• Fast simulations of cosmic large-scale structure with massive neutrinos
• Cosmological Vlasov-Poisson Simulations of Structure Formation with Relic Neutrinos: Nonlinear Clustering and the Neutrino Mass
• Efficient calculation of cosmological neutrino clustering with both linear and non-linear gravity
• The cosmic neutrino background as a collection of fluids in large-scale structure simulations
• Neutrinos in Non-linear Structure Formation - a Simple SPH Approach
• Perturbation theory in Lagrangian hydrodynamics for a cosmological fluid with velocity dispersion
• The Cosmic Neutrino Background Anisotropy - Linear Theory
• A Minimal Model for Massive Neutrinos in Newtonian N-body Simulations
• Multi-Messenger Astrophysics with the Cosmic Neutrino Background

Cited by in corpus (5)

• The FLAMINGO project: cosmological hydrodynamical simulations for large-scale structure and galaxy cluster surveys
• SWIFT: A modern highly-parallel gravity and smoothed particle hydrodynamics solver for astrophysical and cosmological applications
• Euclid: Modelling massive neutrinos in cosmology -- a code comparison
• Hybrid multi-fluid-particle simulations of the cosmic neutrino background
• Euclid preparation. L. Calibration of the linear halo bias in $\Lambda (\nu )$CDM cosmologies