Euclid preparation LXXI. Simulations and nonlinearities beyond $\mathsf{\Lambda}$CDM. 3. Constraints on $f(R)$ models from the photometric primary probes

TL;DR

This paper investigates how different nonlinear matter power spectrum models impact constraints on $f(R)$ gravity from Euclid photometric probes, highlighting biases and the importance of accounting for theoretical uncertainties.

Contribution

It compares multiple nonlinear modeling approaches for $f(R)$ gravity and assesses their effects on parameter constraints and biases in Euclid-like weak lensing and galaxy clustering analyses.

Findings

Bias in $f(R)$ parameter constraints depends on the nonlinear model used.

Baryonic physics can bias $f(R)$ constraints if not properly modeled.

Including theoretical errors reduces bias but weakens constraints.

Abstract

We study the constraint on $f(R)$ gravity that can be obtained by photometric primary probes of the Euclid mission. Our focus is the dependence of the constraint on the theoretical modelling of the nonlinear matter power spectrum. In the Hu-Sawicki $f(R)$ gravity model, we consider four different predictions for the ratio between the power spectrum in $f(R)$ and that in $\Lambda$CDM: a fitting formula, the halo model reaction approach, ReACT and two emulators based on dark matter only $N$-body simulations, FORGE and e-Mantis. These predictions are added to the MontePython implementation to predict the angular power spectra for weak lensing (WL), photometric galaxy clustering and their cross-correlation. By running Markov Chain Monte Carlo, we compare constraints on parameters and investigate the bias of the recovered $f(R)$ parameter if the data are created by a different model. For the pessimistic setting of WL, one dimensional bias for the $f(R)$ parameter, $\log_{10}|f_{R0}|$, is found to be $0.5 \sigma$ when FORGE is used to create the synthetic data with $\log_{10}|f_{R0}| =-5.301$ and fitted by e-Mantis. The impact of baryonic physics on WL is studied by using a baryonification emulator BCemu. For the optimistic setting, the $f(R)$ parameter and two main baryon parameters are well constrained despite the degeneracies among these parameters. However, the difference in the nonlinear dark matter prediction can be compensated by the adjustment of baryon parameters, and the one-dimensional marginalised constraint on $\log_{10}|f_{R0}|$ is biased. This bias can be avoided in the pessimistic setting at the expense of weaker constraints. For the pessimistic setting, using the $\Lambda$CDM synthetic data for WL, we obtain the prior-independent upper limit of $\log_{10}|f_{R0}|< -5.6$. Finally, we implement a method to include theoretical errors to avoid the bias.