No Arabic abstract
We analytically model the non-linear effects induced by massive neutrinos on the total matter power spectrum using the halo model reaction framework of Cataneo et al. 2019. In this approach the halo model is used to determine the relative change to the matter power spectrum caused by new physics beyond the concordance cosmology. Using standard fitting functions for the halo abundance and the halo mass-concentration relation, the total matter power spectrum in the presence of massive neutrinos is predicted to percent-level accuracy, out to $k=10 , h , {rm Mpc}^{-1}$. We find that refining the prescriptions for the halo properties using $N$-body simulations improves the recovered accuracy to better than 1%. This paper serves as another demonstration for how the halo model reaction framework, in combination with a single suite of standard $Lambda$CDM simulations, can recover percent-level accurate predictions for beyond-$Lambda$CDM matter power spectra, well into the non-linear regime.
In the context of forthcoming galaxy surveys, to ensure unbiased constraints on cosmology and gravity when using non-linear structure information, percent-level accuracy is required when modelling the power spectrum. This calls for frameworks that can accurately capture the relevant physical effects, while allowing for deviations from $Lambda$CDM. Massive neutrino and baryonic physics are two of the most relevant such effects. We present an integration of the halo model reaction frameworks for massive neutrinos and beyond-$Lambda$CDM cosmologies. The integrated halo model reaction, combined with a pseudo power spectrum modelled by HMCode2020 is then compared against $N$-body simulations that include both massive neutrinos and an $f(R)$ modification to gravity. We find that the framework is 5% accurate down to at least $kapprox 3 , h/{rm Mpc}$ for a modification to gravity of $|f_{rm R0}|leq 10^{-5}$ and for the total neutrino mass $M_ u equiv sum m_ u leq 0.15$ eV. We also find that the framework is 4(1)% consistent with the Bacco (EuclidEmulator2) emulator for $ u w$CDM cosmologies down to at least $k approx 3 , h$/Mpc. Finally, we compare against hydrodynamical simulations employing HMCode2020s baryonic feedback modelling on top of the halo model reaction. For $ u Lambda$CDM cosmologies we find 2% accuracy for $M_ u leq 0.48$eV down to at least $kapprox 5h$/Mpc. Similar accuracy is found when comparing to $ u w$CDM hydrodynamical simulations with $M_ u = 0.06$eV. This offers the first non-linear and theoretically general means of accurately including massive neutrinos for beyond-$Lambda$CDM cosmologies, and further suggests that baryonic effects can be reliably modelled independently of massive neutrino and dark energy physics. These extensions have been integrated into the publicly available ReACT code.
We introduce an emulator approach to predict the non-linear matter power spectrum for broad classes of beyond-$Lambda$CDM cosmologies, using only a suite of $Lambda$CDM $N$-body simulations. By including a range of suitably modified initial conditions in the simulations, and rescaling the resulting emulator predictions with analytical `halo model reactions, accurate non-linear matter power spectra for general extensions to the standard $Lambda$CDM model can be calculated. We optimise the emulator design by substituting the simulation suite with non-linear predictions from the standard {sc halofit} tool. We review the performance of the emulator for artificially generated departures from the standard cosmology as well as for theoretically motivated models, such as $f (R)$ gravity and massive neutrinos. For the majority of cosmologies we have tested, the emulator can reproduce the matter power spectrum with errors $lesssim 1%$ deep into the highly non-linear regime. This work demonstrates that with a well-designed suite of $Lambda$CDM simulations, extensions to the standard cosmological model can be tested in the non-linear regime without any reliance on expensive beyond-$Lambda$CDM simulations.
We present a general method to compute the nonlinear matter power spectrum for dark energy and modified gravity scenarios with percent-level accuracy. By adopting the halo model and nonlinear perturbation theory, we predict the reaction of a $Lambda$CDM matter power spectrum to the physics of an extended cosmological parameter space. By comparing our predictions to $N$-body simulations we demonstrate that with no-free parameters we can recover the nonlinear matter power spectrum for a wide range of different $w_0$-$w_a$ dark energy models to better than 1% accuracy out to $k approx 1 , h , {rm Mpc}^{-1}$. We obtain a similar performance for both DGP and $f(R)$ gravity, with the nonlinear matter power spectrum predicted to better than 3% accuracy over the same range of scales. When including direct measurements of the halo mass function from the simulations, this accuracy improves to 1%. With a single suite of standard $Lambda$CDM $N$-body simulations, our methodology provides a direct route to constrain a wide range of non-standard extensions to the concordance cosmology in the high signal-to-noise nonlinear regime.
To effectively exploit large-scale structure surveys, we depend on accurate and reliable predictions of non-linear cosmological structure formation. Tools for efficient and comprehensive computational modelling are therefore essential to perform cosmological parameter inference analyses. We present the public software package ReACT, demonstrating its capability for the fast and accurate calculation of non-linear power spectra from non-standard physics. We showcase ReACT through a series of analyses on the DGP and $f(R)$ gravity models, adopting LSST-like cosmic shear power spectra. Accurate non-linear modelling with ReACT has the potential to more than double LSSTs constraining power on the $f(R)$ parameter, in contrast to an analysis that is limited to the quasi-linear regime. We find that ReACT is sufficiently robust for the inference of consistent constraints on theories beyond $Lambda$CDM for current and ongoing surveys. With further improvement, particularly in terms of the accuracy of the non-linear $Lambda$CDM power spectrum, ReACT can, in principle, meet the accuracy requirements for future surveys such as Euclid and LSST.
Numerical simulations of massive neutrino cosmologies consistently find a spoon-like feature in the non-linear matter power spectrum ratios of cosmological models that differ only in the neutrino mass fraction f_N. Typically, the ratio approaches unity at low wave numbers k, decreases by ~ 10 f_N at k ~ 1 h/Mpc, and turns up again at large k. Using the halo model of large-scale structure, we show that this spoon feature originates in the transition from the two-halo power spectrum to the one-halo power spectrum. The formers sensitivity to f_N rises with k, while that of the latter decreases with k. The presence of this spoon feature is robust with respect to different choices of the halo mass function and the halo density profile, and does not require any parameter tuning within the halo model. We demonstrate that a standard halo model calculation is already able to predict the depth, width, and position of this spoon as well as its evolution with redshift z with remarkable accuracy. Predictions at z >= 1 can be further improved using non-linear perturbative inputs.