arxiv: 2512.19568 · v2 · submitted 2025-12-22 · 🌀 gr-qc

Recognition: 2 theorem links

· Lean Theorem

On the Metric f(R) gravity Viability in Accounting for the Binned Supernovae Data

A. Valletta , G. Montani , M. G. Dainotti , E. Fazzari

Authors on Pith no claims yet

Pith reviewed 2026-05-16 20:29 UTC · model grok-4.3

classification 🌀 gr-qc

keywords metric f(R) gravitysupernovae dataHubble constant runningCauchy problemscalar field massJordan framecosmic expansion

0 comments

The pith

An extra condition on the modified Friedmann equation restores consistent dynamics and finite positive scalar mass in metric f(R) gravity while fitting binned supernova data.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper investigates metric f(R) gravity models in the Jordan frame to describe late-time cosmic expansion and explain the effective running of the Hubble constant seen in binned Type Ia supernovae from the Pantheon and Master samples. It defines an effective running Hubble constant as a diagnostic ratio and approximates the gravitational Lagrangian via a low-redshift Taylor expansion of a redshift-dependent function. This general model matches the data but produces an unphysical scalar field mass due to an implicit restriction on the Cauchy problem. Introducing an additional condition on the modified Friedmann equation resolves the issue by restoring a consistent Cauchy problem, yielding a finite positive scalar mass, and preserving the data agreement. The condition is shown to arise from the need for proper dynamical formulation rather than being imposed arbitrarily.

Core claim

The central claim is that the additional condition on the modified Friedmann equation has a precise dynamical origin, being necessary to restore a consistent Cauchy problem and to ensure a finite, positive scalar field mass in metric f(R) gravity models that account for the observed running of the Hubble parameter in binned supernova data.

What carries the argument

The additional condition imposed on the modified Friedmann equation, which enables a fully consistent reformulation of the scalar field dynamics in the general metric f(R) model.

If this is right

The resulting models preserve agreement with the binned Supernova Ia data.
The framework supplies a dynamical justification for the condition used in earlier late-time cosmological analyses.
The scalar degree of freedom acquires a finite and positive mass.
The effective running Hubble constant is reproduced by the low-redshift Taylor approximation.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

Higher-redshift supernova samples could test whether the second-order Taylor expansion remains adequate or demands higher-order corrections in the effective function.
Similar consistency conditions may be required in other scalar-tensor models to avoid unphysical masses when applied to cosmic expansion data.
Precise measurements of the scalar mass from future surveys could provide an independent check on the viability of these f(R) descriptions.

Load-bearing premise

The binned Pantheon and Master samples genuinely indicate a running Hubble parameter that can be captured by a low-redshift Taylor expansion of an effective f(R) function without higher-order terms or selection effects altering the fit.

What would settle it

A calculation or observation showing that the scalar field mass remains negative or divergent in an f(R) model without the additional condition, or new binned supernova data that cannot be fit once the condition is enforced.

Figures

Figures reproduced from arXiv: 2512.19568 by A. Valletta, E. Fazzari, G. Montani, M. G. Dainotti.

**Figure 1.** Figure 1: Reconstruction of H(z) from 1 compared with the binned Pantheon data. at low redshift. In particular, the model should be able to satisfy the luminosityredshift relation, according to Riess et al. (2016); Efstathiou (2021); Fazzari et al. (2025b). Moreover, the parameter U¯ 0 = 0 is fixed without loss of generality, since the fitting procedure depends only weakly on it. Once these cosmographic conditions… view at source ↗

**Figure 2.** Figure 2: Squared mass of the scalar field for the parameters shown in Fig. 1. [PITH_FULL_IMAGE:figures/full_fig_p014_2.png] view at source ↗

**Figure 3.** Figure 3: Reconstruction of H(z) from Tab. 3 compared with the binned Pantheon data. Parameter Value H0 [km s−1 Mpc −1 ] 69.830 ± 0.060 Ωm0 0.3185 ± 0.0048 a −0.025 ± 0.019 b 0.0360 ± 0.0060 w0 −1.0122 ± 0.0093 wa −0.014 ± 0.037 [PITH_FULL_IMAGE:figures/full_fig_p016_3.png] view at source ↗

**Figure 4.** Figure 4: Scalar field mass squared u 2 (z) obtained from the MCMC analysis of the Master Sample. On the contrary, the values obtained for the ΛCDM are χ 2 red, ΛCDM = 2.28 and AICΛCDM = 47.60, BICΛCDM = 48.59. (34) Hence ∆AIC = −2.44, ∆BIC = −4.80. (35) Indicating a weak preference for the AIC and a moderate one for the BIC in favor of the ΛCDM, according to Tab. 2. The scalar field mass squared reconstructed from … view at source ↗

**Figure 5.** Figure 5: Solution of Eq. (43) for α = 1, expressing the required relation between d f dz and d 2 f dz2 at z = 0 for a finite scalar mass. From the initial conditions of the problem, ϕ ′ (0) = 0, see Eq. (15). The analyses performed on both the Pantheon and the Master samples show that the model introduced in Sec. 4.1 is statistically able to reproduce the observed behaviour of the effective H0 at low redshift. Howe… view at source ↗

**Figure 6.** Figure 6: Reconstruction of H(z) from Tab. 4 compared with the binned Pantheon data. The resulting mean values of the fit parameters are reported in Tab. 4, while the same results are shown in [PITH_FULL_IMAGE:figures/full_fig_p023_6.png] view at source ↗

**Figure 7.** Figure 7: Squared mass of the scalar field for this model. [PITH_FULL_IMAGE:figures/full_fig_p024_7.png] view at source ↗

**Figure 8.** Figure 8: Reconstruction of H(z) from Tab. 5 compared with the binned Master data sample. The mean values of the fit parameters values for the model parameters are collected in [PITH_FULL_IMAGE:figures/full_fig_p025_8.png] view at source ↗

**Figure 9.** Figure 9: Squared mass of the scalar field for the Master Sample fit. [PITH_FULL_IMAGE:figures/full_fig_p026_9.png] view at source ↗

**Figure 10.** Figure 10: Comparison between the numerical solution of [PITH_FULL_IMAGE:figures/full_fig_p030_10.png] view at source ↗

**Figure 11.** Figure 11: Comparison between the numerical potential [PITH_FULL_IMAGE:figures/full_fig_p031_11.png] view at source ↗

read the original abstract

In this work, two models of metric $f(R)$ gravity in the Jordan frame are investigated as a dynamical description of the late-time cosmic expansion using binned Type Ia Supernovae data. The aim is to provide an explanation for the effective running of the Hubble constant observed in both the binned Pantheon Sample and the Master Sample. To this end, the effective running Hubble constant $\mathcal{H}(z)$ is defined as the ratio between the modified Hubble parameter and that of the $\Lambda$CDM, multiplied by $H_0$. $\mathcal{H}(z)$ serves as a diagnostic tool to capture deviations from the $\Lambda$CDM model. The first model used is a general representation of metric $f(R)$ gravity in which the gravitational Lagrangian is encoded in an effective redshift-dependent function that mimics the evolution of the Hubble parameter. This function can be approximated by a second-order Taylor expansion at low redshift due to the limited redshift range covered by the Supernovae data. While this general formulation yields a phenomenological fit compatible with that of the $\Lambda$CDM model for the binned Pantheon Sample, the model generically leads to the emergence of an unphysical mass of the scalar field. This issue originates from an implicit restriction imposed on the Cauchy problem for the scalar field. To address this limitation, following previous studies, an additional condition on the modified Friedmann equation is introduced, enabling a fully consistent reformulation of the dynamics. It is clarified that this additional condition has a precise dynamical origin, being necessary to restore a consistent Cauchy problem and to ensure a finite, positive scalar field mass. The resulting framework not only preserves the agreement with binned Supernova Ia data, but also provides a physical justification for the additional condition adopted in earlier analyses of late-time cosmological dynamics.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Desk Editor's Note private letter to a colleague

The paper derives the needed extra condition on the modified Friedmann equation from the Cauchy problem itself so that f(R) models yield a finite positive scalar mass while still fitting the binned supernova data.

read the letter

The main point is that the authors start from a general metric f(R) setup encoded in an effective redshift-dependent function, approximate it with a second-order Taylor expansion at low z, and then show that the usual inconsistency with the scalar-field mass is fixed by an additional condition on the modified Friedmann equation. They trace this condition to the requirement of a well-posed Cauchy problem rather than adding it by hand, which is the clearest new step relative to the papers they cite. The resulting model keeps the effective running Hubble parameter compatible with the binned Pantheon and Master samples, and the scalar mass comes out finite and positive as m² ~ 1/(3f''(R)). That internal consistency argument is solid on its own terms and removes one obvious objection to this class of models. The fit itself remains phenomenological, which the paper states plainly. The soft spot is the dependence on binned data to signal genuine running of H(z). Binning can induce apparent trends even in flat LambdaCDM, and low-z supernova calibration or selection effects are known to produce similar signals, so the data interpretation carries more uncertainty than the consistency fix. The low-order truncation also means higher terms, which feed directly into the mass expression, are assumed negligible; dropping the approximation could reintroduce problems. No new data or machine-checked derivations are supplied, so the strength stays in the dynamical justification. This is useful reading for anyone working on f(R) viability or late-time modified gravity alternatives to dark energy. A colleague focused on the Hubble tension or scalar-tensor constraints would get a concrete check on the Cauchy-problem route. I would send it to peer review because the central claim about the dynamical origin of the condition is worth referee scrutiny on the data handling and expansion assumptions.

Referee Report

2 major / 2 minor

Summary. The manuscript investigates two metric f(R) gravity models in the Jordan frame to explain the effective running of the Hubble constant observed in binned Type Ia Supernovae data from the Pantheon and Master samples. It defines an effective running Hubble parameter H(z) and approximates the gravitational Lagrangian with a second-order Taylor expansion at low redshift. While the general model fits the binned data similarly to LambdaCDM, it leads to unphysical scalar field masses; an additional condition on the modified Friedmann equation is introduced to restore a consistent Cauchy problem and ensure a finite positive scalar mass m^2 ~ 1/(3f''(R)), preserving the data agreement.

Significance. If the result holds, it provides a dynamical justification for the additional condition used in prior f(R) analyses, linking it to the well-posedness of the initial-value problem rather than ad hoc fitting. This strengthens the viability of metric f(R) for late-time cosmology by ensuring physical scalar degrees of freedom, though the reliance on binned data limits broader implications.

major comments (2)

[Derivation of the additional condition] The section deriving the fixing condition: the claim that the additional condition on the modified Friedmann equation has a precise dynamical origin (restoring consistent Cauchy evolution and finite positive scalar mass) requires an explicit step-by-step derivation from the scalar field equation of motion to the mass expression m^2 ~ 1/(3f''(R)); without this, the link to the IVP remains asserted rather than demonstrated.
[§4] §4 (data analysis): the compatibility of the Taylor-expanded f(R) with binned Pantheon/Master samples is shown via H(z) plots, but no reduced chi-squared, degrees of freedom, or test against simulated LambdaCDM binned data is reported; this is load-bearing because binning artifacts could mimic the running without genuine deviation, undermining the viability interpretation.

minor comments (2)

[Abstract and §2] The definition of the effective running Hubble constant H(z) as a ratio multiplied by H0 should include an explicit equation number for clarity.
[Throughout] Notation for the Taylor coefficients of the effective function should be standardized across text and equations to avoid ambiguity in the truncation.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the careful reading of our manuscript and the constructive comments. We address each major comment point by point below. We agree that both points require strengthening and will revise the manuscript accordingly.

read point-by-point responses

Referee: [Derivation of the additional condition] The section deriving the fixing condition: the claim that the additional condition on the modified Friedmann equation has a precise dynamical origin (restoring consistent Cauchy evolution and finite positive scalar mass) requires an explicit step-by-step derivation from the scalar field equation of motion to the mass expression m^2 ~ 1/(3f''(R)); without this, the link to the IVP remains asserted rather than demonstrated.

Authors: We agree that the link to the initial-value problem should be demonstrated explicitly rather than asserted. In the revised manuscript we will insert a dedicated subsection deriving the scalar-field mass from the equation of motion. Starting from the trace of the metric field equations in the Jordan frame, we obtain the second-order equation for the scalar degree of freedom φ = f'(R). Linearizing around a background solution yields the effective mass m² = 1/(3f''(R)) (with the sign convention ensuring m² > 0 for stability). We then show that the additional condition imposed on the modified Friedmann equation is precisely the requirement that the coefficient of the highest derivative remains non-vanishing, thereby restoring a well-posed Cauchy problem and guaranteeing a finite positive mass. This derivation will be self-contained and will cite the relevant literature only for context. revision: yes
Referee: [§4] §4 (data analysis): the compatibility of the Taylor-expanded f(R) with binned Pantheon/Master samples is shown via H(z) plots, but no reduced chi-squared, degrees of freedom, or test against simulated LambdaCDM binned data is reported; this is load-bearing because binning artifacts could mimic the running without genuine deviation, undermining the viability interpretation.

Authors: We acknowledge that visual agreement in the H(z) plots is insufficient to rule out binning artifacts. In the revised §4 we will report the reduced χ² values and the number of degrees of freedom for both the f(R) and ΛCDM fits to the binned Pantheon and Master samples. In addition, we will generate mock binned supernova catalogs from a fiducial ΛCDM cosmology using the same redshift binning and error structure as the real data, and we will repeat the entire fitting procedure on these mocks. The resulting χ² distribution will be used to quantify the probability that the observed running arises from binning alone. These additions will be presented in a new table and accompanying text. revision: yes

Circularity Check

0 steps flagged

Derivation self-contained; additional condition from Cauchy consistency, not data fit

full rationale

The paper defines an effective redshift-dependent f(R) via a low-z Taylor expansion chosen to match binned SN Hubble evolution, then imposes an extra Friedmann condition to ensure a well-posed Cauchy problem and positive finite scalar mass m² ~ 1/(3f''(R)). This condition is justified by dynamical requirements (restoring consistent initial-value evolution) rather than by fitting the supernova data; the fit itself remains an independent phenomenological step. No equation reduces to its input by construction, no self-citation is load-bearing for the central viability claim, and the mass positivity follows from the imposed condition without tautological renaming of the data trend. The derivation is therefore self-contained against external benchmarks.

Axiom & Free-Parameter Ledger

1 free parameters · 2 axioms · 0 invented entities

The central claim rests on the validity of encoding f(R) as a redshift-dependent effective function, the sufficiency of a second-order Taylor expansion over the supernova redshift range, and the interpretation of binned data as evidence for running Hubble parameter.

free parameters (1)

Taylor expansion coefficients
The second-order polynomial coefficients in the effective redshift-dependent function are adjusted to reproduce the binned supernova trend.

axioms (2)

domain assumption Metric f(R) gravity in the Jordan frame admits an effective redshift-dependent Lagrangian function
This is the general representation adopted for the first model.
domain assumption Binned supernova data reflect a genuine effective running of the Hubble parameter
The diagnostic H(z) is defined under this interpretation.

pith-pipeline@v0.9.0 · 5649 in / 1485 out tokens · 31105 ms · 2026-05-16T20:29:21.518346+00:00 · methodology

discussion (0)

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

IndisputableMonolith/Cost/FunctionalEquation.lean washburn_uniqueness_aczel unclear

?

unclear
Relation between the paper passage and the cited Recognition theorem.

the additional condition on the modified Friedmann equation ... necessary to restore a consistent Cauchy problem and to ensure a finite, positive scalar field mass
IndisputableMonolith/Foundation/DimensionForcing.lean reality_from_one_distinction unclear

?

unclear
Relation between the paper passage and the cited Recognition theorem.

μ² = 1/3 (f_R / f_RR − R) > 0

What do these tags mean?

matches: The paper's claim is directly supported by a theorem in the formal canon.
supports: The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends: The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses: The paper appears to rely on the theorem as machinery.
contradicts: The paper's claim conflicts with a theorem or certificate in the canon.
unclear: Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Forward citations

Cited by 3 Pith papers

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score.

On the Physical Nature of the Scalar Mode Mass in the Jordan frame of a Metric $f(R)$ gravity
gr-qc 2026-03 unverdicted novelty 6.0

The scalar mode mass in Jordan-frame metric f(R) gravity must exceed the Hubble scale by several orders of magnitude.
Revisiting the Matter Creation Process: Observational Constraints on Gravitationally Induced Dark Energy and the Hubble Tension
astro-ph.CO 2026-01 conditional novelty 5.0

Gravitationally induced particle creation models fit cosmological data as well as ΛCDM and reduce the Hubble tension from 4.3σ to 2.4–3σ.
Mapping the redshift drift at various redshifts through cosmography
astro-ph.CO 2026-04 unverdicted novelty 4.0

Cosmographic Taylor and Padé models fitted to Pantheon+SH0ES+GRB+DESI BAO data yield redshift drift predictions compatible with ΛCDM and ω0ω1CDM at 1-2σ, with mock drift data tightening q0 and j0 bounds.

Reference graph

Works this paper leans on

5 extracted references · 5 canonical work pages · cited by 3 Pith papers · 3 internal anchors

[1]

Gravitational Stability and Bulk Cosmology

doi:10.3847/1538-4357/ac8e04. Capozziello, S., Cardone, V .F., Elizalde, E., Nojiri, S., Odintsov, S.D., 2006. Observational constraints on dark energy with generalized equations of state. Physical Review D—Particles, Fields, Gravitation, and Cosmology 73, 043512. Capozziello, S., De Laurentis, M., 2011. Extended theories of gravity. Physics Reports 509, ...

work page internal anchor Pith review Pith/arXiv arXiv doi:10.3847/1538-4357/ac8e04 2006
[2]

arXiv preprint arXiv:2510.19931

Determining the hubble constant through cross-correlation of galaxies and gravitational waves. arXiv preprint arXiv:2510.19931 . Perlmutter, S., et al., 1999. Measurements ofωandλfrom 42 high-redshift su- pernovae. The Astrophysical Journal 517, 565–586. doi:10.1086/307221. Postnikov, S., Dainotti, M.G., Hernandez, X., Capozziello, S., 2014. Nonpara- metr...

work page internal anchor Pith review doi:10.1086/307221 1999
[3]

Silva, E., Sabogal, M.A., Scherer, M., Nunes, R.C., Di Valentino, E., Kumar, S.,

doi:10.3847/1538-4357/aab9bb. Silva, E., Sabogal, M.A., Scherer, M., Nunes, R.C., Di Valentino, E., Kumar, S.,

work page doi:10.3847/1538-4357/aab9bb
[4]

f(R) Theories Of Gravity

New constraints on interacting dark energy from desi dr2 bao observa- tions. Physical Review D 111, 123511. Sotiriou, T.P., Faraoni, V ., 2010. f(r) theories of gravity. Reviews of Modern Physics 82, 451–497. doi:10.1103/RevModPhys.82.451, arXiv:0805.1726. Su, P., Gong, Y ., Xiong, Q., Hu, D., Lin, H., Deng, F., Chen, X., 2025. Exploring joint observation...

work page internal anchor Pith review Pith/arXiv arXiv doi:10.1103/revmodphys.82.451 2010
[5]

arXiv preprint arXiv:2503.15659

A low-redshift preference for an interacting dark energy model. arXiv preprint arXiv:2503.15659 . Zhan, Y ., Wang, D., Yi, S.X., Wang, F.Y ., 2025. Precise hubble constant measure- ment from quasi-periodic eruptions as electromagnetic counterparts to extreme mass ratio inspirals. arXiv preprint arXiv:2506.14150 . 43

work page arXiv 2025