Pith. Sign in

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

How this review was made: An AI referee (deepseek/deepseek-v4-flash) read the paper on September 11, 2026 and answered a fixed list of questions about it. T0 is the starting level every paper gets. The levels describe how deeply Pith has checked the paper, never the authors.

  1. T0 AI review An AI read the paper and answered a fixed list of questions.
  2. T1 Audit Automatic checks that the references and files the paper relies on exist and match.
  3. T2 Compute The paper's code or calculations re-run with the same results.
  4. T3 Formal statement The claims written out in a formal math language a computer can read.
  5. T4 Proof Those claims proved, and the proof checked by computer.

How Pith checks papers →

· 3 major objections · 7 minor · Ask about this paper · The compiled record · Read the full review →
Challenge this review Comment (0) Read on arXiv ↗
Cite ▾
arXiv ID2609.10651 (version 1)
Pith ID4EK6TQFV, Pith's permanent record of this paper and its review
CitationAlvaro de la Cruz-Dombriz, Savvas Nesseris, Tshepo Mathibela. "Constraining $f(R)$ gravity and evolving dark energy via large-scale structure and phase-space trajectories." arXiv:2609.10651 (2026). https://pith.science/paper/4EK6TQFV
Cite this reviewBibTeX for Pith's review of this paper
Claim paper

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

The authors confront two viable f(R) modified gravity models—the Hu & Sawicki and Starobinsky models—with combined background and growth-of-structure data (BAO, supernovae, CMB distance priors, and measurements of fσ₈, f, and σ₈). Using MCMC sampling, they obtain tight constraints on the deviation parameters: log₁₀ b_HS ≈ −6.3 and b_S ≈ 8×10⁻⁵, both consistent with General Relativity. Model comparison via the Akaike Information Criterion shows a mild preference for the f(R) models over flat ΛCDM (|ΔAIC| ≥ 4), but the Bayesian Information Criterion does not confirm this preference (|ΔBIC| ≤ 2). The paper also introduces a new diagnostic using two-dimensional phase-space diagrams in the (μ, γ) and (μ, Σ) planes, where ΛCDM corresponds to the fixed point (1,1). If future weak-lensing and galaxy surveys find points with μ−1 < 0 and γ−1 > 1 or Σ−1 < 0, the Hu & Sawicki and Starobinsky models could be ruled out.

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

Editorial extensions of the paper, not claims the author makes directly.

  • 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.

Compiled by classifier-1.13.0 via TypeSafe API, question bank abstract-p1-20260923, 2026-09-23: 1 request, 6,140 input tokens.

What the abstract itself states, as probabilities: claims first 0.15 · states limitation 0.85 · human participants 0.01 · controlled comparison 0.72 · randomized assignment 0.01 · states quantitative result 0.99 · compares with prior results 0.90 · mentions code or data availability 0.02.

The abstract, sentence by sentence

arXiv listing abstract v1. Results and the main claim are underlined; .

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

What the classifier reports for each sentence. The number is the probability it gave the label.
#sentencerolethe authors' own workstates the main conclusionstrength
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 kind of paper the classifier reports this is: empirical study 0.68 · evidence from analysis of existing data 1.00.

What was examined, and what was not

sources read
  • arXiv listing abstract v1
not examined
  • full text
  • figures and tables
  • references and citing sentences
  • relations to other papers
levels run
  • P abstract
not run
  • S sections
  • T body sentences
  • C citations
  • R cross-paper relations

Every quoted sentence on this record is the authors'. Every label is from Pith's closed vocabulary, with the probability shown beside it. Every number in the measurements table is the authors', checked by code. Record as JSON

Figures
Figure 1. The evolution of the normalized potentials
Figure 2. The redshift evolution of the growth of structure
Figure 3. The evolution of the light deflection parameter
Figure 4. The evolution of the gravitational light deflection parameter
Figure 5. The constraints for the parameters {Ωm,0, h, σ8} for the standard ΛCDM model. The three contour levels correspond to the 1σ, 2σ, and 3σ interval levels. {27, 3}. For a given set of cosmological parameters θ, the residual vectors for the quantities ∈ {j} are given by, ∆X = fσ8(zi , θ) − (fσ8)obs, (59) ∆Y = f(zi , θ) − (f)obs, (60) ∆Q = σ8(zi , θ) − (σ8)obs, (61) where zi are the i th observed redshifts and (fσ8)obs,(f)obs,(σ8)obs are the observed data for the quantities in {j} as shown in Table IV. V. MCMC RESULTS In this section we provide the cosmological fits for Hu & Sawicki, and…
Figure 6. Constraints for the parameters {Ωm,0, h, σ8, log10bHS} for the standard Hu & Sawicki model. The three contour levels correspond to the 1σ, 2σ, and 3σ interval levels. the parameters {w0, wa} are set to {−1, 0}. Lastly, the initial values for {log10bHS, bS} = {−3, 0.01} for the Hu & Sawicki and Starobinsky models respectively. For all combinations of datasets, we adopt the uniform wide flat priors on the parameters as presented in Table I. Parameter Prior Ωm,0 U[0.1, 0.6] h U[0.4, 1] σ8,0 U[0.2, 0.9] w0 U[−2, 1] wa U[−3, 2] log10bHS U[−8, 0] bS U[−0.5, 0.5] Table I. Priors imposed on…
Figure 8. The constraints for the parameters {Ωm,0, h, σ8} for the w0waCDM model. The three contour levels correspond to the 1σ, 2σ, and 3σ interval levels
Figure 9. The scale and redshift evolution of the LSS functions
Figure 10. The scale and redshift evolution of the LSS functions

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.