Do cosmological data rule out $f(\mathcal{R})$ with $w\neq-1$?

TL;DR

This paper uses the Equation of State approach to analyze $f( ext{R})$ gravity models of dark energy, constraining their parameters with current cosmological data and showing they are tightly limited when $w eq -1$.

Contribution

It provides a numerically stable EoS framework for $f( ext{R})$ models and derives observational constraints on their parameters using Planck and BAO data.

Findings

$B_0<0.006$ for $w=-1$ models

$B_0<0.0045$ and $|w+1|<0.002$ for $w eq -1$ models

$f( ext{R})$ models are more tightly constrained than $w$CDM models

Abstract

We review the Equation of State (EoS) approach to dark sector perturbations and apply it to $f(\mathcal{R})$ gravity models of dark energy. We show that the EoS approach is numerically stable and use it to set observational constraints on designer models. Within the EoS approach we build an analytical understanding of the dynamics of cosmological perturbations for the designer class of $f(\mathcal{R})$ gravity models, characterised by the parameter $B_0$ and the background equation of state of dark energy $w$. When we use the Planck Cosmic Microwave Background (CMB) temperature anisotropy, polarisation and lensing data as well as the Baryonic Acoustic Oscillation (BAO) data from SDSS and WiggleZ, we find $B_0<0.006$ (95\%CL) for the designer models with $w=-1$. Furthermore, we find $B_0<0.0045$ and $|w+1|<0.002$ (95\%CL) for the designer models with $w\neq -1$. Previous analyses found similar results for designer and Hu-Sawicki $f(\mathcal{R})$ gravity models using the Effective Field Theory (EFT) approach (Raveri et al. 2014, Hu et al. 2016), therefore this hints for the fact that generic $f(\mathcal{R})$ models with $w\neq-1$ can be tightly constrained by current cosmological data, complementary to solar system tests (Brax et al. 2008, Faulkner et al. 2007). When compared to a $w$CDM fluid with the same sound speed, we find that the equation of state for $f(\mathcal{R})$ models is better constrained to be close to -1 by about an order of magnitude, due to the strong dependence of the perturbations on $w$.