Spectral Chebyshev Approximation of Cosmic Expansion in f(R) Gravity

Abstract

We present a numerical framework to study the cosmological background evolution in f(R) gravity by employing a spectral Chebyshev collocation approach. Unlike standard integration methods such as Runge--Kutta that often encounter stiffness and accuracy issues, this formulation expands the normalized Hubble function E(z) = H(z)/H0 as a finite Chebyshev series. The modified Friedmann equation is then enforced at selected Chebyshev--Gauss--Lobatto points, converting the original nonlinear differential equation into a system of algebraic relations for the series coefficients. This transformation yields exponentially convergent and numerically stable solutions over the entire redshift domain, 0<z<zmax, eliminating the need for adaptive step-size control. We apply the method to two widely studied f(R) models, Hu--Sawicki and Starobinsky, and perform a combined analysis using cosmic chronometer H(z) data and the Union~3.0 supernova compilation. The reconstructed expansion histories match observations to within 2σ over 0 < z < 2, producing best-fit parameters of approximately (m0, H0, eff) (0.29, 68, 1.2--2.5\,H02). These results indicates that both models reproduce the observed late-time acceleration while permitting small geometric corrections to . Overall, the spectral Chebyshev method provides a precise and computationally efficient framework for probing modified-gravity cosmologies in the precision-data era.

0

Turn this paper into a full lesson

ArcXiv compiles a staged curriculum from this paper: 8-12 lessons across beginner → advanced, synthesised section guides, visuals, flashcards, a quiz, exercises, and on-demand deep dives per section. Grounded in the abstract, never invented.

Discussion (0)

Sign in to join the discussion.

Loading comments…