Higher-order initial conditions with massive neutrinos

Abstract

The discovery that neutrinos have mass has important consequences for cosmology. The main effect of massive neutrinos is to suppress the growth of cosmic structure on small scales. Such growth can be accurately modelled using cosmological $N$-body simulations, but doing so requires accurate initial conditions (ICs). There is a trade-off, especially with first-order ICs, between truncation errors for late starts and discreteness and relativistic errors for early starts. Errors can be minimized by starting simulations at late times using higher-order ICs. In this paper, we show that neutrino effects can be absorbed into scale-independent coefficients in higher-order Lagrangian perturbation theory (LPT). This clears the way for the use of higher-order ICs for massive neutrino simulations. We demonstrate that going to higher order substantially improves the accuracy of simulations. To match the sensitivity of surveys like DESI and Euclid, errors in the matter power spectrum should be well below 1%. However, we find that first-order Zel'dovich ICs lead to much larger errors, even when starting as early as $z=127$, exceeding 1% at $z=0$ for $k>0.5{\text{ Mpc}}^{-1}$ for the power spectrum and $k>0.1{\text{ Mpc}}^{-1}$ for the equilateral bispectrum in our simulations. Ratios of power spectra with different neutrino masses are more robust than absolute statistics, but still depend on the choice of ICs. For all statistics considered, we obtain 1% agreement between 2LPT and 3LPT at $z=0$.

15 pages, 9 figures, with link to software, minor update of version submitted to MNRAS

References in corpus (21)

• Massive neutrinos and cosmology
• Transients from Initial Conditions in Cosmological Simulations
• Large-scale dark matter simulations
• The Effect of Thermal Neutrino Motion on the Non-linear Cosmological Matter Power Spectrum
• Impact of massive neutrinos on nonlinear matter power spectrum
• On the most constraining cosmological neutrino mass bounds
• Accurate initial conditions for cosmological N-body simulations: Minimizing truncation and discreteness errors
• Hunting for Galaxies and Halos in simulations with VELOCIraptor
• Massive Neutrinos in Cosmology: Analytic Solutions and Fluid Approximation
• Structure formation with massive neutrinos: going beyond linear theory
• Relativistic N-body simulations with massive neutrinos
• An optimal nonlinear method for simulating relic neutrinos
• Higher-order initial conditions for mixed baryon-CDM simulations
• The halo model in a massive neutrino cosmology
• Linear perturbative theory of the discrete cosmological N-body problem
• A Lagrangian Perturbation Theory in the presence of massive neutrinos
• How smooth are particle trajectories in a $\Lambda$CDM Universe?
• The DESI $N$-body Simulation Project I: Testing the Robustness of Simulations for the DESI Dark Time Survey
• Cosmological $N$-body simulations including radiation perturbations
• Cosmological perturbations for two cold fluids in $\Lambda$CDM
• The cosmic neutrino background as a collection of fluids in large-scale structure simulations

Cited by in corpus (12)

• The FLAMINGO project: cosmological hydrodynamical simulations for large-scale structure and galaxy cluster surveys
• Fast computation of non-linear power spectrum in cosmologies with massive neutrinos
• Aemulus $\nu$: Precise Predictions for Matter and Biased Tracer Power Spectra in the Presence of Neutrinos
• Precision Redshift-Space Galaxy Power Spectra using Zel'dovich Control Variates
• Cosmological constraints from the cross-correlation of DESI Luminous Red Galaxies with CMB lensing from Planck PR4 and ACT DR6
• The FLAMINGO project: Baryon effects on the matter power spectrum
• Geodesic motion and phase-space evolution of massive neutrinos
• The FLAMINGO project: Baryonic impact on weak gravitational lensing convergence peak counts
• Cosmological feedback from a halo assembly perspective
• Improving initialization and evolution accuracy of cosmological neutrino simulations
• BullFrog: Multi-step perturbation theory as a time integrator for cosmological simulations
• The FLAMINGO simulations data release