Improved analytical modeling of the nonlinear power spectrum in modified gravity cosmologies

Reliable analytical modeling of the nonlinear power spectrum (PS) of matter perturbations is among the chief prerequisites for cosmological analyses from the largest sky surveys. This is especially true for the models that extend the standard general-relativity paradigm by adding the fifth force, where numerical simulations can be prohibitively expensive. Here we present a method for building accurate PS models for two modified gravity (MG) variants: namely the Hu-Sawicki f (R ) , and the normal branch of the Dvali-Gabadadze-Porrati (nDGP) braneworld. We start by modifying the standard Halo Model (HM) with respect to the baseline lambda-cold-dark-matter (Λ CDM ) scenario, by using the HM components with specific MG extensions. We find that our P (k )_HM retains 5% accuracy only up to mildly nonlinear scales (k ≲0.3 h /Mpc ) when compared to PS from numerical simulations. At the same time, our HM prescription much more accurately captures the ratio ϒ (k )=P (k )_MG/P (k )_{Λ CDM} up to nonlinear scales. We show that using HM-derived ϒ (k ) together with a viable nonlinear Λ CDM P (k ) prescription (such as HALOFIT), we render a much better and more accurate PS predictions in MG. The new approach yields considerably improved performance, with modeled P (k )M_G being now accurate to within 5% all the way to nonlinear scales of k ≲2.5 - 3 h /Mpc . The magnitude of deviations from GR as fostered by these MG models is typically O (10 %) in these regimes. Therefore reaching 5% PS modeling is enough for forecasting constraints on modern-era cosmological observables.