Kinematic model-independent reconstruction of Palatini $f(R)$ cosmology

TL;DR

This paper develops a model-independent method using rational approximations to reconstruct Palatini $f(R)$ cosmology from luminosity distance data, revealing slight deviations from General Relativity and providing new insights into cosmic expansion.

Contribution

It introduces a kinematic, model-independent approach with rational approximations to reconstruct $f(R)$ gravity in the Palatini formalism using observational data.

Findings

Reconstructed $f(R)$ models closely match numerical solutions.

Identified specific functional forms for $f(R)$ with fitted parameters.

Compared results with $ ext{Lambda}$CDM and discussed deviations from GR.

Abstract

A kinematic treatment to trace out the form of $f(R)$ cosmology, within the Palatini formalism, is discussed by only postulating the universe homogeneity and isotropy. To figure this out we build model-independent approximations of the luminosity distance through rational expansions. These approximants extend the Taylor convergence radii computed for usual cosmographic series. We thus consider both Pad\'e and the rational Chebyshev polynomials. They can be used to accurately describe the universe late-time expansion history, providing further information on the thermal properties of all effective cosmic fluids entering the energy momentum tensor of Palatini's gravity. To perform our numerical analysis, we relate the Palatini's Ricci scalar with the Hubble parameter $H$ and thus we write down a single differential equation in terms of the redshift $z$. Therefore, to bound $f(R)$, we make use of the most recent outcomes over the cosmographic parameters obtained from combined data surveys. In particular our clue is to select two scenarios, i.e. $(2,2)$ Pad\'e and $(2,1)$ Chebyshev approximations, since they well approximate the luminosity distance at the lowest possible order. We find that best analytical matches to the numerical solutions lead to $f(R)=a+bR^n$ with free parameters given by the set $(a, b, n)=(-1.627, 0.866, 1.074)$ for $(2,2)$ Pad\'e approximation, whereas $f(R)=\alpha+\beta R^m$ with $(\alpha, \beta, m)=(-1.332, 0.749, 1.124)$ for $(2,1)$ rational Chebyshev approximation. Finally, our results are compared with the $\Lambda$CDM predictions and with previous studies in the literature. Slight departures from General Relativity are also discussed.