REVIEW
Constraining $f(R)$ gravity and evolving dark energy via large-scale structure and phase-space trajectories
General Physics · Submitted 9 Sep 2026 · 58 references
Pith's verdict: Solid MCMC analysis of f(R) models, but the missing wavenumber specification for growth predictions is a load-bearing flaw that undermines the claimed constraints.
T0 · AI review
Why Pith says this
Highlight any text to comment on it.
Referee report
The paper
The paper claims
This paper uses large-scale structure data to constrain two popular f(R) gravity models and constructs a null-test based on phase-space trajectories of gravitational parameters that could rule them out with future observations.
The reading
What carries the argument
The dynamical system approach for the background evolution (variables x, y, χ, Ω, h) and the linear perturbation equations in the sub-horizon approximation, which produce the effective gravitational coupling μ, the slip parameter γ, and the light-deflection parameter Σ for f(R) gravity. The chi-squared likelihood construction for multiple datasets and the MCMC sampling (using wide flat priors) are the statistical machinery that delivers parameter constraints. The phase-space null-test then uses the predicted trajectories of (μ,γ) and (μ,Σ) across redshifts and scales to define the region where these models live.
Extended reading notes
Core claim
The central claim is that the Hu & Sawicki and Starobinsky f(R) models, with their deviation parameters set to best-fit values, are not excluded by current cosmological probes, but the data prefer ΛCDM when the Bayesian penalty for extra parameters is accounted for. The paper further claims that the allowed region in the (μ−1, γ−1) and (μ−1, Σ−1) phase spaces occupies the second and first quadrants respectively, providing a crisp null-test: any future data point falling in the opposite quadrants would falsify these models.
Load-bearing premise
The analysis assumes that the Hu & Sawicki and Starobinsky models with n=1 are the correct description of dark energy, that the fifth force is unscreened on cosmological scales (i.e., the chameleon mechanism does not suppress it at the scales probed), and that the sub-horizon approximation (k ≫ aH) holds for the perturbation equations down to k ≈ 30 H₀.
Editorial extensions
If this is right
- If the AIC-based mild preference for f(R) models holds, future data with smaller errors could either strengthen the preference or rule it out, pushing the deviation parameters further toward the GR limit.
- The phase-space null-test provides a model-specific but observationally accessible way to discard the Hu & Sawicki and Starobinsky models using upcoming lensing and clustering surveys (e.g., Euclid, Rubin) without requiring explicit parameter fitting.
- The constraints on the DESI dark-energy equation of state parameters w₀ and wₐ are also reported, allowing comparison with w₀wₐCDM models that show a stronger AIC preference over ΛCDM.
- The fact that the BIC does not support f(R) models over ΛCDM suggests that the additional degree of freedom is not warranted by the current data, reinforcing the standard cosmological model.
- The analytical expressions for μ, γ, and Σ in the sub-horizon limit can be used in future full-shape analyses of galaxy clustering and weak lensing to perform a joint test of f(R) gravity.
Reading between the lines
- The phase-space diagnostic introduced here could be extended to other f(R) models beyond Hu & Sawicki and Starobinsky; if those models also occupy the same quadrants, the test becomes a broader null-test for a class of modified gravity theories.
- Because the sub-horizon approximation may break down at very large scales (k ~ 30 H₀), future surveys with wide angular coverage might see deviations from the predicted trajectories that are not due to the f(R) models but to the approximation itself; a full calculation without the sub-horizon assumption would strengthen the diagnostic.
- The double-parameter comparison (AIC vs. BIC) highlights the well-known tension between model selection criteria; a reader interested in Bayesian model comparison might want to see the Bayes factor explicitly, which the paper does not compute.
- The chameleon screening mechanism is invoked to pass solar-system tests but is not modelled in the growth equations; the derived constraints on b_HS and b_S assume the fifth force is active on cosmological scales, so a failure of screening on those scales would require reinterpreting the bounds.
Original abstract
We present a joint observational analysis confronting viable $f(R)$ modified gravity theories, specifically the Hu \& Sawicki and Starobinsky models, with background and large-scale structure (LSS) data. Utilizing Monte-Carlo Markov chain (MCMC) sampling across datasets including baryon acoustic oscillations (BAO), type Ia supernovae (SNeIa), cosmic microwave background (CMB) distance priors, and linear growth measurements ($f\sigma_8$, $f$, $\sigma_8$), we place tight constraints on the model parameters governing deviations from General Relativity. For the full dataset combination, we obtain $\log_{10} b_\mathrm{HS} = -6.325_{-1.138}^{+1.216}$ for the Hu \& Sawicki model and $b_\mathrm{S} = (0.8\pm61.0)\times10^{-4}$ for the Starobinsky model. Model comparison based on the Akaike Information Criterion indicates that these $f(R)$ extensions are statistically favored over flat $\Lambda\text{CDM}$ ($|\Delta\text{AIC}| \ge 3.99$) for the combined data. However, when considering the Bayesian Information Criterion, the evidence for support is significantly reduced. Furthermore, we construct two-dimensional phase-space diagrams in the $(\mu, \gamma)$ and $(\mu, \Sigma)$ planes across several redshifts, establishing a novel diagnostic null-test allowing us to probe for deviations from $\Lambda\text{CDM}$, corresponding to the fixed point $(1,1)$ in both planes, using LSS observables. Should future weak-lensing and galaxy surveys provide data points with $\mu-1<0$ and $\gamma-1>1$ or $\Sigma-1 < 0 $, then the aforementioned models could be directly ruled out.
The compiled record
The paper's own sentences, arranged, kept beside the review and separate from it. Each label is drawn from a closed vocabulary; the number beside it is the probability the classifier gave that label. Numbers in the tables are the authors', checked by code. No sentence in this record was written by a model.
The abstract, sentence by sentence
We present a joint observational analysis confronting viable $f(R)$ modified gravity theories, specifically the Hu \& Sawicki and Starobinsky models, with background and large-scale structure (LSS) data. Utilizing Monte-Carlo Markov chain (MCMC) sampling across datasets including baryon acoustic oscillations (BAO), type Ia supernovae (SNeIa), cosmic microwave background (CMB) distance priors, and linear growth measurements ($f\sigma_8$, $f$, $\sigma_8$), we place tight constraints on the model parameters governing deviations from General Relativity. For the full dataset combination, we obtain $\log_{10} b_\mathrm{HS} = -6.325_{-1.138}^{+1.216}$ for the Hu \& Sawicki model and $b_\mathrm{S} = (0.8\pm61.0)\times10^{-4}$ for the Starobinsky model. Model comparison based on the Akaike Information Criterion indicates that these $f(R)$ extensions are statistically favored over flat $\Lambda\text{CDM}$ ($|\Delta\text{AIC}| \ge 3.99$) for the combined data. However, when considering the Bayesian Information Criterion, the evidence for support is significantly reduced. Furthermore, we construct two-dimensional phase-space diagrams in the $(\mu, \gamma)$ and $(\mu, \Sigma)$ planes across several redshifts, establishing a novel diagnostic null-test allowing us to probe for deviations from $\Lambda\text{CDM}$, corresponding to the fixed point $(1,1)$ in both planes, using LSS observables. Should future weak-lensing and galaxy surveys provide data points with $\mu-1<0$ and $\gamma-1>1$ or $\Sigma-1 < 0 $, then the aforementioned models could be directly ruled out.
background contribution method result interpretation main claim
| # | sentence | role | the authors' own work | states the main conclusion | strength |
|---|---|---|---|---|---|
| s0 | We present a joint observational analysis confronting viable $f(R)$ modified gravity theories, specifically the Hu \& Sawicki and Starobinsky models, with background and large-scale structure (LSS) data. | contribution 1.00 | 0.93 | 0.05 | asserted 0.98 |
| s1 | Utilizing Monte-Carlo Markov chain (MCMC) sampling across datasets including baryon acoustic oscillations (BAO), type Ia supernovae (SNeIa), cosmic microwave background (CMB) distance priors, and linear growth measurements ($f\sigma_8$, $f$, $\sigma_8$), we place tight constraints on the model parameters governing deviations from General Relativity. | method 1.00 | 0.95 | 0.07 | asserted 1.00 |
| s2 | For the full dataset combination, we obtain $\log_{10} b_\mathrm{HS} = -6.325_{-1.138}^{+1.216}$ for the Hu \& Sawicki model and $b_\mathrm{S} = (0.8\pm61.0)\times10^{-4}$ for the Starobinsky model. | result 1.00 | 0.96 | 0.12 | asserted 1.00 |
| s3 | Model comparison based on the Akaike Information Criterion indicates that these $f(R)$ extensions are statistically favored over flat $\Lambda\text{CDM}$ ($|\Delta\text{AIC}| \ge 3.99$) for the combined data. | result 0.91 | 0.94 | 0.33 | hedged possibility 0.51 |
| s4 | However, when considering the Bayesian Information Criterion, the evidence for support is significantly reduced. | result 0.53 | 0.84 | 0.09 | asserted 0.98 |
| s5 | Furthermore, we construct two-dimensional phase-space diagrams in the $(\mu, \gamma)$ and $(\mu, \Sigma)$ planes across several redshifts, establishing a novel diagnostic null-test allowing us to probe for deviations from $\Lambda\text{CDM}$, corresponding to the fixed point $(1,1)$ in both planes, using LSS observables. | contribution 0.88 | 0.96 | 0.12 | asserted 1.00 |
| s6 | Should future weak-lensing and galaxy surveys provide data points with $\mu-1<0$ and $\gamma-1>1$ or $\Sigma-1 < 0 $, then the aforementioned models could be directly ruled out. | future work 0.59 | 0.19 | 0.09 | hedged possibility 1.00 |
What was examined, and what was not
- arXiv listing abstract v1
- full text
- figures and tables
- references and citing sentences
- relations to other papers
- P abstract
- S sections
- T body sentences
- C citations
- R cross-paper relations
Figures
Showing 9 of the paper's 10 figures. See all figures in the PDF ↗ Figures reproduced from arXiv:2609.10651 by the authors.
References
Discussion
Disagree with the review, or know something it missed? Leave a comment, highlight any passage above to comment on it, or challenge this review formally.