ABSTRACT Cosmology is poised to measure the neutrino mass sum Mν and has identified several smaller-scale observables sensitive to neutrinos, necessitating accurate predictions of neutrino clustering over a wide range of length scales. The FlowsForTheMasses non-linear perturbation theory for the the massive neutrino power spectrum, $\Delta ^2_\nu (k)$, agrees with its companion N-body simulation at the $10~{{\ \rm per\ cent}}-15~{{\ \rm per\ cent}}$ level for k ≤ 1 h Mpc−1. Building upon the Mira-Titan IV emulator for the cold matter, we use FlowsForTheMasses to construct an emulator for $\Delta ^2_\nu (k)$, Cosmic-Eν, which covers a large range of cosmological parameters and neutrino fractions Ων, 0h2 ≤ 0.01 (Mν ≤ 0.93 eV). Consistent with FlowsForTheMasses at the 3.5 per cent level, it returns a power spectrum in milliseconds. Ranking the neutrinos by initial momenta, we also emulate the power spectra of momentum deciles, providing information about their perturbed distribution function. Comparing a Mν = 0.15 eV model to a wide range of N-body simulation methods, we find agreement to 3 per cent for k ≤ 3kFS = 0.17 h Mpc−1 and to 19 per cent for k ≤ 0.4 h Mpc−1. We find that the enhancement factor, the ratio of $\Delta ^2_\nu (k)$ to its linear-response equivalent, is most strongly correlated with Ων, 0h2, and also with the clustering amplitude σ8. Furthermore, non-linearities enhance the free-streaming-limit scaling $\partial \log (\Delta ^2_\nu /\Delta ^2_{\rm m}) / \partial \log (M_\nu)$ beyond its linear value of 4, increasing the Mν-sensitivity of the small-scale neutrino density.
Read full abstract