REVIEW 3 major objections 5 minor 61 references
Inferences for f(R) Models from Late-Time Megamaser Observational Data
T0 review · 3 major / 5 minor · reviewed 2026-08-01 · deepseek-v4-flash
Pith's one-line read Using megamaser angular-diameter distances, the paper finds three f(R) gravity models statistically indistinguishable from ΛCDM, with the deviation parameter b near zero and H0 near 73 km/s/Mpc.
desk verdict A legitimate but overstated application of P20 megamaser data to three f(R) models: the H0 anchor is fine, but the b posteriors do not support the claim that f(R) mimics ΛCDM. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The load-bearing object is the deviation parameter b that enters the Hubble-rate expressions H(z; H0, b, Ωm) for the three f(R) models, each written as ΛCDM plus polynomial terms in b. Because b → 0 reduces each model to ΛCDM, the posterior on b measures how strongly data demand modified-gravity corrections. The analysis also leans on a combined likelihood that adds a peculiar-velocity uncertainty of 250 km/s in quadrature to each recession-velocity measurement, and on MCMC sampling to marginalize over the six galaxy velocities alongside H0, b, and Ωm.
What would settle it
Re-fit the six megamaser objects with the exact, untruncated H(z) expressions for each f(R) model; if the marginalized b then shifts away from zero by more than the quoted 1σ uncertainty, the paper's central conclusion is falsified. Alternatively, a future sample of ~30 megamaser hosts with 3% distance precision that rules out b = 0 at more than 2σ would overturn it.
Extended reading notes
Core claim
The central claim is that the megamaser sample — six galaxies with geometric angular-diameter distances and recession velocities — yields marginalized estimates of the f(R) deviation parameter b close to zero for all three models (b ≈ 0.011 ± 1.00 for Hu-Sawicki, 0.023 ± 1.04 for Starobinsky, −0.005 ± 1.01 for ArcTanh), while H0 is constrained near 73 km/s/Mpc in every model. The paper interprets b ≈ 0 as the signature that these f(R) models reduce to the ΛCDM expansion history at low redshift. Since the 1σ uncertainties on b are of order unity, the claim is not that b is tightly measured, only that the data do not prefer any departure from ΛCDM. The model-comparison statistics (ΔAIC, ΔBIC ≤
Load-bearing premise
The weakest point is the assumption that the H(z) formulas for the three f(R) models, imported from earlier work and truncated at low order in b, are the correct Hubble rates for those theories over the full prior range of b; if that mapping is wrong, the b ≈ 0 conclusion is an artifact.
Editorial extensions
If this is right
- If correct, the megamaser geometry independently supports H0 ≈ 73 km/s/Mpc, reinforcing the disagreement with early-universe CMB-based estimates.
- The b ≈ 0 result implies that, for background expansion, these f(R) models are not distinguishable from ΛCDM with current data; any discriminating power must come from structure growth or higher redshift.
- The weak Ωm constraint (≈0.5 with σ≈0.4) demonstrates that low-redshift distance data alone cannot determine the matter density, so combined analyses with SNe Ia, BAO, or H(z) would be needed to sharpen all parameters.
- AIC/BIC differences of order 2 mean that, by standard model-selection rules, the f(R) models are statistically equivalent to ΛCDM; the extra parameter b is not penalized into preference either way.
Reading between the lines
- Editorial inference: Because the b posterior is only constrained to |b| ≲ 1, the polynomial truncation in the H(z) expressions could matter; recomputing with higher-order terms would test whether b ≈ 0 is robust.
- Editorial inference: The paper works at the background level; one natural extension is to add growth-rate data (e.g., fσ8 measurements) to see whether the same f(R) models, with b consistent with zero from background, remain viable when perturbations are included.
- Editorial inference: The fixed 250 km/s peculiar-velocity uncertainty is conservative; a full marginalization over per-object peculiar velocities with realistic priors could tighten or widen the quoted errors, though the paper's stochastic check suggests little change.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper uses six megamaser angular-diameter distances and recession velocities from the Megamaser Cosmology Project (Pesce et al. 2020) to constrain H0, Omega_m, six recession-velocity nuisance parameters, and, for the three f(R) models, a deviation parameter b. The likelihood combines velocity and distance terms with a peculiar-velocity uncertainty of 250 km/s. MCMC posteriors give H0 around 73.5–73.7 km/s/Mpc for all four models, Omega_m near 0.5 with errors around 0.35–0.4, and b near 0 with errors around 1.0. AIC/BIC comparisons are used to claim that the f(R) models are statistically indistinguishable from LambdaCDM, and the paper concludes that this late-time dataset predicts that f(R) models mimic LambdaCDM.
Significance. The H0 constraint, if correct, is an interesting independent late-time geometric probe favoring the local value near 73 km/s/Mpc, and the paper demonstrates a clean likelihood treatment of the P20 sample. However, the advertised f(R) result is not supported by the reported errors: the b posteriors have 1-sigma widths comparable to the prior, so the data do not constrain b. The model-comparison section also contains an internal inconsistency with its own table. The paper is therefore more valuable for its H0 result than for its f(R) discrimination, which is effectively non-constraining at the current sample size.
major comments (3)
- [Abstract; Table 3] The claim that the P20 dataset 'predicts that f(R) models mimic LambdaCDM' is not supported by the quoted constraints. Table 3 gives b = 0.011 +/- 1.00, 0.023 +/- 1.04, and -0.005 +/- 1.01, while the prior is U(-1.5, 1.5). These 1-sigma intervals cover most of the prior; the data do not constrain b away from the prior. The abstract and Section 4 should be softened to state that the data are consistent with b=0 but provide no significant evidence either way.
- [Sec. 2.3, Eqs. (7)-(9)] The Hubble rates for the three f(R) models are imported from Sultana et al. (2022) with only a parenthetical reference. The paper does not define the f(R) actions, the parameter b, or the expansion used. Eqs. (7)-(9) are truncated series in b (up to order b^2 or b^4), yet the prior allows |b| up to 1.5 and the posteriors have errors of order 1, so the validity of the truncation is not established. The authors must state the model actions and b definition, justify the truncation, or use exact expressions; otherwise the b constraints are not interpretable.
- [Sec. 3.1, Table 5] The text states that '0 <= |Delta X| <= 2' and uses this to conclude that the f(R) models are statistically indistinguishable from LambdaCDM. Table 5 lists Delta AIC = 2.15-2.16 and Delta BIC = 2.64-2.65, which contradict the stated inequality. The comparison should be rephrased using a defined evidence scale, and the table values must be reconciled with the text.
minor comments (5)
- [Abstract] Typo: 'predict' should be 'predicts'.
- [Introduction] Typo: 'Active Galctic Nucleus' should be 'Active Galactic Nucleus'.
- [Sec. 4] The notation 'b- -> 0' and 'b- -> infinity' is unclear; presumably this should be the single parameter b.
- [References] Several references have minor formatting issues, e.g., 'V ol.' in Schinckel et al. and the Hogg (1999) entry is incomplete.
- [Sec. 2.2] The robustness check with sigma_pec drawn from U(150,250) is only described qualitatively. Please report the resulting parameter shifts or remove the claim.
Circularity Check
The claim that P20 predicts f(R) mimics ΛCDM rests on a fit of the deviation parameter b whose posterior equals the prior and whose zero value is ΛCDM by construction; the 'prediction' is therefore inherited from the ansatz.
-
fitted input called prediction
[Abstract; Section 2.3, Eqs. (7)–(9); Table 2; Table 3; Section 4]
"The marginalized estimates of the deviation parameter, b, for the three f(R) gravity models lie close to zero. This late-time dataset, thus predicts that f(R) models mimic ΛCDM. ... In the limit b−→0, these models are constructed to reproduce the ΛCDM to ensure consistency with current observational constraints."
Each of Eqs. (7)–(9) is written as the ΛCDM expression 1−Ωm+(1+z)^3 Ωm plus polynomial terms multiplied by powers of b, so b=0 returns ΛCDM by construction. Table 3 reports b=0.011±1.00, 0.023±1.04, −0.005±1.01, while Table 2 sets the prior b∼U(−1.5,1.5); the 1σ widths are essentially the full prior width, so the P20 data do not constrain b. The central value near zero is the mean of a prior-dominated, nearly flat posterior. Calling this a 'prediction' that f(R) mimics ΛCDM is thus a restatement of the model definition and prior, not an independent, data-driven result.
full rationale
The P20 angular-diameter-distance and velocity data are external to the f(R) constructions, and the H0 constraint (~73 km/s/Mpc) is largely model-independent, so that part of the analysis is not circular. The paper's central 'f(R) mimics ΛCDM' claim, however, is not independently established: the Hubble rates in Eqs. (7)–(9) are imported without derivation from Sultana et al. (2022), and b is constructed so that b=0 is exactly ΛCDM. Because the fitted b posteriors (0.011±1.00, 0.023±1.04, −0.005±1.01) are as wide as the prior U(−1.5,1.5), the data carry no information about b, and the abstract's 'predicts' overstates what is a prior/ansatz-derived central value. The model-comparison section is also internally inconsistent (Section 3.1 says '0≤|ΔX|≤2' while Table 5 lists ΔAIC≈2.15 and ΔBIC≈2.64), which is a correctness issue rather than circularity. Overall, the core inference about b and ΛCDM mimicry reduces, at least in part, to the construction of b and the choice of prior; hence a moderate circularity score is warranted.
Assumptions & free parameters
free parameters (4)
- H0 =
73.583±2.99 (ΛCDM); 73.726±4.79 (Hu-Sawicki); 73.519±3.46 (Starobinsky); 73.458±4.96 (ArcTanh) km/s/Mpc
- Ωm =
0.511±0.39 (ΛCDM); 0.514±0.38 (Hu-Sawicki); 0.515±0.39 (Starobinsky); 0.505±0.34 (ArcTanh)
- b (f(R) deviation parameter) =
0.011±1.00 (Hu-Sawicki); 0.023±1.04 (Starobinsky); -0.005±1.01 (ArcTanh)
- v_i (six galaxy recession velocities, i=1..6) =
Table 4: v1≈3514, v2≈10190, v3≈7813, v4≈8535, v5≈7088, v6≈559 km/s
assumptions (5)
- standard math The FLRW angular-diameter distance formula (Eq. 5) applies for all models.
- domain assumption Equations (7)-(9) correctly give the Hubble rates of the Hu-Sawicki, Starobinsky, and ArcTanh f(R) models.
- domain assumption The low-redshift relation z_i = v_i/c (Eq. 1) is accurate for the megamaser sample.
- domain assumption A single peculiar-velocity uncertainty σ_pec=250 km/s applies to all six galaxies.
- domain assumption The combined Gaussian likelihood (Eq. 4) with flat priors is a valid statistical model for the data.
Cite this review
Pith. "Pith review of Inferences for f(R) Models from Late-Time Megamaser Observational Data." pith.science (2026). https://pith.science/paper/LRO6A7GL
@misc{pith2026260721129,
author = {Pith},
title = {Pith review of: Inferences for f(R) Models from Late-Time Megamaser Observational Data},
year = {2026},
howpublished = {\url{https://pith.science/paper/LRO6A7GL}},
note = {Machine review of arXiv:2607.21129}
}
abstract
In this work, we study three widely used models of f(R) gravity, namely HuSawicki, Starobinsky and ArcTanh along with the standard cosmology model ($\Lambda$CDM). For this, we employ the megamaser angular diameter distance and velocity measurements from the Megamaser Cosmology Project, which provide a purely geometric determination of the Hubble constant. We constrain the parameters using the Markov Chain Monte Carlo method. Our results show that values of the Hubble Constant, $H_{0}$, obtained for all four models are in concordance with its value obtained from other late-time observational data such as SNe Ia. The constraints on $H_{0}$ in all the models under study are restrictive and the marginalized estimates lie close to 73 $\mathrm {km s^{-1} Mpc^{-1}}$. The marginalized estimates of the deviation parameter, b, for the three f(R) gravity models lie close to zero. This late-time dataset, thus predicts that f(R) models mimic $\Lambda$CDM. However, the matter density, $\Omega_m$, remains weakly constrained for all the models with its marginalized estimate close to 0.5. Further, comparison of the four models (f(R) models and $\Lambda$CDM) using information criteria such as Akaike Information Criterion and Bayesian Information Criterion shows that within current uncertainties, the dataset finds f(R) models statistically indistinguishable from $\Lambda$CDM. This is consistent with the fact that the favoured value of b for each of the f(R) models lies close to zero.
Figures
Figures from the paper (5 more)
Reference graph
Works this paper leans on
-
[1]
2007, Phys
Amendola, L., Gannouji, R., Polarski, D., & Tsujikawa, S. 2007, Phys. Rev. D, 75, 083504 2
2007
-
[2]
2010, Dark energy: theory and observations (Cambridge University Press) 2
Amendola, L., & Tsujikawa, S. 2010, Dark energy: theory and observations (Cambridge University Press) 2
2010
-
[3]
Armendariz-Picon, C., Mukhanov, V ., & Steinhardt, P. J. 2001, Phys. Rev. D, 63, 103510 2
2001
-
[4]
2025, ApJ, 992, 205 3
Bai, J., Xia, J.-Q., & Zhao, G.-B. 2025, ApJ, 992, 205 3
2025
-
[5]
Bamba, K., Capozziello, S., Nojiri, S., & Odintsov, S. D. 2012, Ap&SS, 342, 155 2
2012
-
[6]
M., & Alcaniz, J
Barboza, E. M., & Alcaniz, J. S. 2008, Physics Letters B, 666, 415 2
2008
-
[7]
2025, Ap&SS, 370, 62 3, 6, 7
Barua, S., Ramakrishnan, V ., & Desai, S. 2025, Ap&SS, 370, 62 3, 6, 7
2025
-
[8]
2013, Phys
Basilakos, S., Nesseris, S., & Perivolaropoulos, L. 2013, Phys. Rev. D, 87, 123529 2
2013
Show all 61 references
-
[9]
J., Cappellari, M., et al
Birrer, S., Buckley-Geer, E. J., Cappellari, M., et al. 2025, Astronomy & Astrophysics, 704, A63 12
2025
-
[10]
2022, ApJ, 938, 110 12
Brout, D., et al. 2022, ApJ, 938, 110 12
2022
-
[11]
Caldwell, R. R. 2002, Physics Letters B, 545, 23 2
2002
-
[12]
Campista, M., Santos, B., Santos, J., & Alcaniz, J. S. 2011, Physics Letters B, 699, 320 2 Inferences forf(R)Models 17
2011
-
[13]
M., Duvvuri, V ., Trodden, M., & Turner, M
Carroll, S. M., Duvvuri, V ., Trodden, M., & Turner, M. S. 2004, Phys. Rev. D, 70, 043528 2
2004
-
[14]
M., Hoffman, M., & Trodden, M
Carroll, S. M., Hoffman, M., & Trodden, M. 2003, Phys. Rev. D, 68, 023509 2
2003
-
[15]
2020, Research in Astronomy and Astrophysics, 20, 074 15
Chen, R.-R., Zhang, H.-Y ., Jin, C.-J., et al. 2020, Research in Astronomy and Astrophysics, 20, 074 15
2020
-
[16]
2001, International Journal of Modern Physics D, 10, 213 2
Chevallier, M., & Polarski, D. 2001, International Journal of Modern Physics D, 10, 213 2
2001
-
[17]
2003, Physics Letters B, 575, 1 2
Chiba, T. 2003, Physics Letters B, 575, 1 2
2003
-
[18]
G., Padilla, A., & Skordis, C
Clifton, T., Ferreira, P. G., Padilla, A., & Skordis, C. 2012, Phys. Rep., 513, 1 2 D’Agostino, R., & Nunes, R. C. 2019, Phys. Rev. D, 100, 044041 3 D’Agostino, R., & Nunes, R. C. 2020, Phys. Rev. D, 101, 103505 3 Dark Energy Survey and Kilo-Degree Survey Collaboration, Abbott...
2012
-
[19]
2008, Phys
Dev, A., Jain, D., Jhingan, S., et al. 2008, Phys. Rev. D, 78, 083515 2
2008
-
[20]
D., & Kawasaki, M
Dolgov, A. D., & Kawasaki, M. 2003, Physics Letters B, 573, 1 2
2003
-
[21]
Feng, C.-J., Shen, X.-Y ., Li, P., & Li, X.-Z. 2012, J. Cosmol. Astropart. Phys., 2012, 023 2
2012
-
[22]
Freedman, W. L. 2021, ApJ, 919, 16 1
2021
-
[23]
Gelman, A., & Rubin, D. B. 1992, Statistical science, 7, 457 6
1992
-
[24]
R., Moran, J
Herrnstein, J. R., Moran, J. M., Greenhill, L. J., et al. 1999, Nature, 400, 539 3
1999
-
[25]
Hogg, D. W. 1999, arXiv e-prints, astro 14
1999
-
[26]
2007, Phys
Hu, W., & Sawicki, I. 2007, Phys. Rev. D, 76, 064004 5
2007
-
[27]
Hui, L., & Greene, P. B. 2006, Phys. Rev. D, 73, 123526 15
2006
-
[28]
c., & Noh, H
Hwang, J. c., & Noh, H. 2001, Physics Letters B, 506, 13 2
2001
-
[29]
K., Bagla, J
Jassal, H. K., Bagla, J. S., & Padmanabhan, T. 2005, MNRAS, 356, L11 2
2005
-
[30]
K., Ray, S., & Zhang, F
Kumar, D., Dhankar, P. K., Ray, S., & Zhang, F. 2025, Physics of the Dark Universe, 49, 101989 3, 6, 12, 14, 15, 16
2025
-
[31]
C., Pan, S., & Yadav, P
Kumar, S., Nunes, R. C., Pan, S., & Yadav, P. 2023, Physics of the Dark Universe, 42, 101281 12, 14, 15, 16
2023
-
[32]
Leizerovich, M., Kraiselburd, L., Landau, S., & Sc´occola, C. G. 2022, Phys. Rev. D, 105, 103526 3
2022
-
[33]
Lewis, A. 2025, J. Cosmol. Astropart. Phys., 2025, 025 6
2025
-
[34]
Linder, E. V . 2003, Phys. Rev. Lett., 90, 091301 2
2003
-
[35]
Ludwick, K. J. 2017, Modern Physics Letters A, 32, 1730025 2
2017
-
[36]
2023, Fortschritte der Physik, 71, 2300133 2
Majumder, B., Ray, S., & Manna, G. 2023, Fortschritte der Physik, 71, 2300133 2
2023
-
[37]
Nojiri, S., & Odintsov, S. D. 2003, Phys. Rev. D, 68, 123512 2
2003
-
[38]
Nojiri, S., & Odintsov, S. D. 2007, International Journal of Geometric Methods in Modern Physics, 04, 115 2
2007
-
[39]
Nojiri, S., & Odintsov, S. D. 2011, Phys. Rep., 505, 59 2
2011
-
[40]
D., & Oikonomou, V
Nojiri, S., Odintsov, S. D., & Oikonomou, V . K. 2017, Phys. Rep., 692, 1 2
2017
-
[41]
C., Pan, S., Saridakis, E
Nunes, R. C., Pan, S., Saridakis, E. N., & Abreu, E. M. 2017, Journal of Cosmology and Astroparticle Physics, 2017, 005 12, 14, 15, 16
2017
-
[42]
C., Pan, S., Saridakis, E
Nunes, R. C., Pan, S., Saridakis, E. N., & Abreu, E. M. C. 2017, J. Cosmol. Astropart. Phys., 2017, 005 2
2017
-
[43]
D., G´omez, D
Odintsov, S. D., G´omez, D. S.-C., & Sharov, G. S. 2019, Physical Review D, 99, 024003 3
2019
-
[44]
D., S´aez-Chill´on G´omez, D., & Sharov, G
Odintsov, S. D., S´aez-Chill´on G´omez, D., & Sharov, G. S. 2021, Nuclear Physics B, 966, 115377 3 O’Sullivan, N. K., & Aigrain, S. 2024, MNRAS, 531, 4181 12, 13 Pascual-S´anchez, J.-F. 1999, Modern Physics Letters A, 14, 1539 2
2021
-
[45]
J., & Ratra, B
Peebles, P. J., & Ratra, B. 2003, Reviews of Modern Physics, 75, 559 1 P´erez-Romero, J., & Nesseris, S. 2018, Phys. Rev. D, 97, 023525 2, 6
2003
-
[46]
W., Braatz, J
Pesce, D. W., Braatz, J. A., Reid, M. J., et al. 2020, ApJ, 891, L1 3, 4, 15 Planck Collaboration, Aghanim, N., Akrami, Y ., et al. 2020, A&A, 641, A6 14, 16
2020
-
[47]
2026, Classical and Quantum Gravity, 43, 085012 3
Ravi, K. 2026, Classical and Quantum Gravity, 43, 085012 3
2026
-
[48]
J., Braatz, J
Reid, M. J., Braatz, J. A., Condon, J. J., et al. 2009, ApJ, 695, 287 3
2009
-
[49]
J., Pesce, D
Reid, M. J., Pesce, D. W., & Riess, A. G. 2019, ApJ, 886, L27 3
2019
-
[50]
Santos, B., Campista, M., Santos, J., & Alcaniz, J. S. 2012, A&A, 548, A31 2 18 Umang et al
2012
-
[51]
E., Bunton, J
Schinckel, A. E., Bunton, J. D., Cornwell, T. J., Feain, I., & Hay, S. G. 2012, in Society of Photo- Optical Instrumentation Engineers (SPIE) Conference Series, V ol. 8444, Ground-based and Airborne Telescopes IV , ed. L. M. Stepp, R. Gilmozzi, & H. J. Hall, 84442A 15
2012
-
[52]
2007, Phys
Song, Y .-S., Hu, W., & Sawicki, I. 2007, Phys. Rev. D, 75, 044004 2
2007
-
[53]
P., & Faraoni, V
Sotiriou, T. P., & Faraoni, V . 2010, Reviews of Modern Physics, 82, 451 2
2010
-
[54]
1980, Physics Letters B, 91, 99 2
Starobinsky, A. 1980, Physics Letters B, 91, 99 2
1980
-
[55]
Starobinsky, A. A. 2007, Soviet Journal of Experimental and Theoretical Physics Letters, 86, 157 5, 14
2007
-
[56]
K., Melia, F., & Kazanas, D
Sultana, J., Yennapureddy, M. K., Melia, F., & Kazanas, D. 2022, MNRAS, 514, 5827 3, 4, 13
2022
-
[57]
2016, arXiv e-prints, arXiv:1601.05256 2
Sundell, P., & Vilja, I. 2016, arXiv e-prints, arXiv:1601.05256 2
2016 arXiv
-
[58]
2010, 800, 99 14
Tsujikawa, S. 2010, 800, 99 14
2010
-
[59]
2013, Classical and Quantum Gravity, 30, 214003 2
Tsujikawa, S. 2013, Classical and Quantum Gravity, 30, 214003 2
2013
-
[60]
2021, Statistical Science, 36, 518 6
Vats, D., & Knudson, C. 2021, Statistical Science, 36, 518 6
2021
-
[61]
2000, arXiv e-prints, astro 1
Weinberg, S. 2000, arXiv e-prints, astro 1
2000
Reviewed August 1, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.