REVIEW 4 major objections 5 minor 1 cited by
Bayesian Inference of Neutron Star Properties in $f(Q)$ Gravity Using NICER Observations
T0 review · 4 major / 5 minor · reviewed 2026-08-03 · deepseek-v4-flash
Pith's one-line read Bayesian analysis of NICER pulsar data favors the exponential f(Q) gravity model and predicts neutron-star maximum masses near 3 solar masses, extending into the lower mass gap.
desk verdict First Bayesian NICER constraint on f(Q) neutron-star parameters, but the preferred nonlinear models violate the paper's own Eq. (28), so the headline claims are unsupported as written. 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 central object is the nonmetricity scalar Q(r) = −2 e^{−λ}(ν′ + 1/r)/r built from the static spherical metric, along with the modified TOV system (Eqs. 25–34) that generalizes hydrostatic equilibrium to f(Q) gravity. A Bayesian sampler maps the f(Q) parameters (α, β, Q₀) to neutron-star observables through these equations, and Bayes factors rank the three functional forms.
What would settle it
For the best-fit exponential model parameters, compute the left-hand side of Eq. (28), (cotθ/2) Q′ f_QQ, along the stellar profile; if it does not vanish, the reported mass–radius solutions are not solutions of the field equations, and the statistical ranking would need to be recomputed with the full teleparallel connection.
Extended reading notes
Core claim
Confronting f(Q)-gravity neutron-star models with NICER data does more than constrain free parameters—it yields a sharp prediction. The exponential model, f(Q) = −Q + αQ₀[1 − exp(−β√(Q/Q₀))], is favored over the linear and logarithmic models by Bayes factors, with its scale parameter Q₀ tightly constrained. All three models show 95% confidence bands of the mass–radius relation extending to nearly 3 M⊙, crossing into the lower mass gap (2.5–5 M⊙). The paper interprets this as evidence that nonmetricity-based gravity enhances high-density support, making mass-gap neutron stars a natural, testable outcome.
Load-bearing premise
The analysis assumes that the coincident-gauge, static spherically symmetric metric ansatz with a purely radial nonmetricity Q(r) satisfies the full f(Q) field equations; for the logarithmic and exponential models this is not the case, because Eq. (28) requires either f_QQ = 0 or Q′ = 0, and those models have neither.
Editorial extensions
If this is right
- Exponential f(Q) is the statistically preferred model among the three, with well-constrained Q₀ and near-unimodal posteriors.
- All three f(Q) models predict 95% confidence regions for the maximum mass that reach about 2.98 M⊙, placing stable configurations in the lower mass gap.
- For the exponential model, the radius and tidal deformability at 1.4 M⊙ (R₁.₄ = 11.27 km, Λ₁.₄ = 156.95) are consistent with NICER and gravitational-wave constraints.
- The mass-gap prediction is a clear, testable signature: precise mass measurements of neutron stars above 2.5 M⊙ would support f(Q) gravity.
- The logarithmic model predicts the highest median maximum mass (2.36 M⊙), yet is less preferred statistically than the exponential model.
Reading between the lines
- If the mass-gap prediction holds, it offers a rare direct probe of symmetric teleparallel gravity in the strong-field regime: a population of 2.5–3 M⊙ neutron stars would be difficult to explain in general relativity with conventional equations of state.
- The paper fixes the DDME2 equation of state to isolate gravity effects; a joint inference over both EoS and f(Q) parameters, or the inclusion of GW170817 tidal constraints as an additional likelihood, could tighten or overturn the model ranking.
- The Bayesian methodology could be extended to other f(Q) functional forms or to torsion-based teleparallel gravity, enabling a systematic observational ranking of teleparallel gravity models.
- The reported preference for the exponential model might be sensitive to the chosen priors; a prior-sensitivity analysis would clarify whether the conclusion is robust.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript performs a Bayesian inference for three f(Q) gravity models (linear, logarithmic, and exponential), with the DDME2 equation of state fixed, using NICER mass–radius posteriors for PSR J0030+0451, PSR J0740+6620, PSR J0437+4715, and PSR J0614+3329. It derives modified TOV equations, computes mass–radius relations, tidal deformabilities, and maximum masses, and reports that the exponential model is statistically preferred while all models allow maximum masses extending toward ~2.98 M_sun and into the lower mass gap. The central technical problem is that the nonlinear logarithmic and exponential solutions do not satisfy the full field equations as written, because Eq. (28) requires either f_QQ=0 or Q'=0, neither of which is true for these models. The reported posterior distributions, mass–radius bands, and Bayes factors therefore are not supported by the stated theory.
Significance. If the results were valid, this would be a useful first systematic Bayesian constraint on f(Q) gravity from neutron-star observations, with a falsifiable mass-gap prediction. The paper addresses a currently underexplored comparison between symmetric teleparallel gravity and NICER data, and the Bayesian setup is standard. However, the main conclusions rest on nonlinear f(Q) solutions that are not admitted by the coincident-gauge ansatz used in the derivation. The exponential-model preference and the mass-gap prediction are precisely the parts of the analysis that are invalidated by this inconsistency. The paper also contains internal inconsistencies in the reported Bayes factors. These are not merely presentational issues; they affect the central claims.
major comments (4)
- [Sec. III, Eq. (28)] Equation (28), an independent component of the field equations for the metric (23), states 0 = (cotθ/2) Q' f_QQ. For the logarithmic model (36), f_QQ = -α/Q^2, and for the exponential model (37), f_QQ is nonzero in the posterior region. The numerical Q(r) from Eq. (29) is not constant, so Q' ≠ 0. Thus the nonlinear solutions do not satisfy the full field equations. The paper integrates the reduced system (25)-(27)/(33)-(34) and uses the resulting M-R curves in the NICER likelihood and Bayes factors, but these curves are not solutions of the stated theory. The linear model (35) is consistent because f_QQ=0, but the headline claims concern the logarithmic and exponential models. The standard remedy, a non-trivial teleparallel connection as in Ref. [23], changes the modified TOV system; this is not implemented or discussed. This is a load-bearing gap.
- [Table IV and Sec. VIII] The Bayes-factor evidence is internally inconsistent. Table IV lists ln(BF_12) = -1.81 for linear versus exponential and -0.43 for logarithmic versus exponential, with exponential as H2. The concluding section states 'values of −0.35 and −1.28' and describes the evidence as strong or decisive. On the paper's own scale in Sec. VI B 2, ln(BF) = -0.43 is only weak-to-moderate evidence for the exponential model, not 'definitely the best.' Please correct the numbers and re-evaluate the model-ranking claim.
- [Table I and Sec. IV B-C] The prior for Q0 is given as uniform on [-1,1], which includes Q0=0. Both the logarithmic model, f(Q) = -Q + α ln(β Q/Q0), and the exponential model, f(Q) = -Q + α Q0 [1 - exp(-β sqrt(Q/Q0))], are singular or undefined at Q0=0. No exclusion, regularization, or reparameterization is described. As written, the Bayesian models are not well defined on the stated prior support. If the implementation de facto excludes Q0=0, that restriction must be stated explicitly and motivated.
- [Sec. III, Eq. (33)] Equation (33) is printed as λ' = -kr(ρ+P)e^λ - ν'/f(Q). In the units used throughout the paper, f(Q) has dimension length^{-2}, so the term ν'/f(Q) has dimension length, while λ', r(ρ+P), and the left-hand side have dimension length^{-1}. This equation is therefore dimensionally inconsistent as it stands, and the GR limit f(Q) = -Q does not reproduce the standard TOV equation. If this is a typographical error, the corrected form of the modified TOV system must be given explicitly, since this is the system that is numerically integrated to produce all the paper's results.
minor comments (5)
- [Abstract] The sentence 'and exhibits well-constrained' is incomplete; the quantity that is well constrained should be stated.
- [Sec. III, Eq. (30)] Equation (30) contains two equal signs and is not a well-formed equation as printed. Please correct the displayed expression so that the pressure-gradient equation is unambiguous.
- [Sec. VII B and Fig. 4] The statement that the M-R curves 'pass through the centers of the observational bands' is an in-sample consistency check, because the f(Q) parameters were fitted to exactly those NICER posteriors. It is not by itself evidence of model success. A posterior predictive check or out-of-sample test would be more informative.
- [Sec. VII B / Table III] The abstract and Sec. VII B refer to 95% confidence regions, while Table III reports 68%(90%) intervals. Please clarify the coverage used in the quoted maximum-mass values and in Fig. 4.
- [Sec. V] The boundary condition λ(0)=0 is non-standard for the coordinate r; in Schwarzschild-like coordinates one usually has e^{-λ(0)}=1 but λ(0) itself can be set to zero by redefinition. The shooting condition e^{ν(R)}=e^{-λ(R)} should be stated in terms of the asymptotic Schwarzschild mass, which is not explicitly defined.
Circularity Check
No significant circularity; the central inference chain is self-contained, though it has a separate field-equation consistency issue.
full rationale
The paper's derivation chain is not circular. The f(Q) model parameters are fitted to NICER mass-radius posteriors through the likelihood in Eq. (39), and the reported R_1.4, Lambda_1.4, M-R curves, and M_max are posterior outputs of that same Bayesian procedure. The statement that M-R curves "pass through the centers of the observational bands" is an in-sample fit description, not an independent prediction, and the paper does not disguise it as a forecast. The headline forward-looking claim, that "all three constrained models predict maximum neutron star masses reaching M_max ~ 2.98 M_sun" and that this "mass-gap prediction emerges naturally from the Bayesian-constrained parameter space," is a posterior extrapolation: no mass-gap object was used as a fitting target, so the result does not reduce to the inputs by construction. Self-citations such as refs. [25], [26], and [64] are used only to motivate functional forms and cosmological viability; they are not load-bearing for the TOV derivation or the Bayesian ranking. The main substantive weakness noted by the reader—Eq. (28), 0 = (cotθ/2) Q' f_QQ, requiring either f_QQ = 0 or Q' = 0 while the logarithmic and exponential models have f_QQ ≠ 0 and r-dependent Q—is a potential internal inconsistency or omitted consistency condition, not a circularity; it belongs to correctness assessment rather than to this circularity pass. Therefore, no circular step can be exhibited under the required standard.
Assumptions & free parameters
free parameters (3)
- α =
linear 0.28, logarithmic 1.55, exponential −5.18 (posterior medians)
- β =
linear −1.03, logarithmic 1.11, exponential 2.64 (posterior medians)
- Q0 =
linear 0.02, logarithmic −0.26, exponential −0.35 (posterior medians)
assumptions (6)
- domain assumption The coincident-gauge static spherically symmetric ansatz (metric Eq. 23, Q from Eq. 29) is a valid description of non-linear f(Q) neutron stars.
- domain assumption The exterior vacuum is the Schwarzschild metric and the matching e^ν(R) = e^{−λ}(R) is sufficient.
- domain assumption The DDME2 EoS captures dense-matter physics adequately.
- domain assumption The hypermomentum satisfies ∇_μ ∇_ν H^{μν}_α = 0.
- domain assumption Published NICER mass–radius posteriors can be used directly as likelihoods.
- ad hoc to paper The prior ranges in Table I are broad enough and well defined on the support of the models.
Cite this review
Pith. "Pith review of Bayesian Inference of Neutron Star Properties in $f(Q)$ Gravity Using NICER Observations." pith.science (2026). https://pith.science/paper/NYMCDWYD
@misc{pith2026260116227,
author = {Pith},
title = {Pith review of: Bayesian Inference of Neutron Star Properties in $f(Q)$ Gravity Using NICER Observations},
year = {2026},
howpublished = {\url{https://pith.science/paper/NYMCDWYD}},
note = {Machine review of arXiv:2601.16227}
}
abstract
In this work, we investigate neutron stars (NSs) in the strong field regime within the framework of symmetric teleparallel $f(Q)$ gravity, considering three representative models: linear, logarithmic, and exponential. While Bayesian studies of NS observations are well established in general relativity and curvature based modified gravity theories, such analyses in $f(Q)$ gravity remain largely unexplored. For the first time we perform a Bayesian inference analysis by confronting theoretical NS mass-radius predictions with NICER observations of PSR J0030+0451, PSR J0740+6620, PSR J0437+4715, and PSR J0614+3329 in the background of nonmetricity based gravity. The dense matter equation of state is fixed to DDME2 in order to isolate the effects of modified gravity on NS structure. Our results show that the exponential $f(Q)$ model is statistically preferred over the linear and logarithmic cases, as confirmed by Bayes factor comparisons, and exhibits well-constrained. For this model, we obtain a radius and tidal deformability at $1.4\,M_\odot$ of $R_{1.4} = 11.27^{+0.53}_{-0.36}\,\mathrm{km}$ and $\Lambda_{1.4} = 156.95^{+84.02}_{-41.73}$, respectively, consistent with current observational constraints. Remarkably, all three constrained models predict maximum neutron star masses reaching $M_{\max} \simeq 2.98\,M_{\odot}$, with the $95\%$ confidence regions extending into the lower mass gap ($\sim 2.5$--$5\,M_{\odot}$). This mass-gap prediction emerges naturally from the Bayesian-constrained parameter space. These results highlight the potential of NSs as powerful probes of symmetric teleparallel gravity in the strong field regime.
Figures
Figures from the paper (3 more)
Forward citations
Cited by 1 Pith paper
-
Astrophysical Objects in Modified Theories of Gravity
The work constructs analytical compact-star solutions in torsion- and nonmetricity-based gravity and shows that modified-gravity parameters can alter maximum mass and stability while remaining compatible with NICER ob...
Reference graph
Works this paper leans on
-
[23]
J. B. Jim´ enez, L. Heisenberg, and T. Koivisto, Phys. Rev. D98, 044048 (2018)
2018
-
[1]
Log of the Bayes Factor Taking thelogarithmof the Bayes factor provides a more practical scale for interpretation and calculation [77]. The log Bayes factor is: ln(BF12) = ln P(D|H 1) P(D|H 2) .(42) This transformation provides a more practical scale for quantifying the strength of evidence that favors one hypothesis. TABLE I: Uniform priors on thef(Q) mo...
-
[2]
•0<ln(BF)<1: Weak to moderate evidence for hypothesisH 1
Interpretation of the Bayes Factor To assess the strength of evidence, the values of the Bayes factor are often assessed using the following conventional scale [78, 79]: •ln(BF)>1: Strong evidence for hypothesisH 1. •0<ln(BF)<1: Weak to moderate evidence for hypothesisH 1. •ln(BF) = 0: No evidence for either hypothesis. •−1<ln(BF)<0: Weak to moderate evid...
2025
-
[3]
Einstein, Annalen der Physik49, 769 (1916)
A. Einstein, Annalen der Physik49, 769 (1916)
1916
-
[4]
R. M. Wald,General Relativity(Chicago Univ. Pr., Chicago, USA, 1984)
1984
-
[5]
C. M. Will, Living Reviews in Relativity17(2014), 10.12942/lrr-2014-4
-
[6]
T. Clifton, P. G. Ferreira, A. Padilla, and C. Skordis, Phys. Rep.513, 1 (2012), arXiv:1106.2476 [astro-ph.CO]
arXiv 2012
-
[7]
S. e. a. Perlmutter, Astrophysical Journal517, 565 (1999)
1999
Show all 81 references
-
[8]
A. G. e. a. Riess, Astron. J.116, 1009 (1998), arXiv:astro-ph/9805201
1998 arXiv
-
[9]
D. N. e. a. Spergel, Astrophysical Journal Supplement Series170, 377 (2007)
2007
-
[10]
M. e. a. Tegmark, Phys. Rev. D71, 103515 (2005)
2005
-
[11]
Sahni and A
V. Sahni and A. A. Starobinsky, International Journal of Modern Physics D9, 373 (2000)
2000
-
[12]
Padmanabhan, Physics Reports380, 235 (2003)
T. Padmanabhan, Physics Reports380, 235 (2003)
2003
-
[13]
E. V. Linder, General Relativity and Gravitation40, 329 (2008)
2008
-
[14]
Frieman, M
J. Frieman, M. Turner, and D. Huterer, Annual Review of Astronomy and Astrophysics46, 385 (2008)
2008
-
[15]
R. R. Caldwell and M. Kamionkowski, Annual Review of Nuclear and Particle Science59, 397 (2009)
2009
-
[16]
De Felice and S
A. De Felice and S. Tsujikawa, Living Reviews in Relativity13, 3 (2010)
2010
-
[17]
T. P. Sotiriou and V. Faraoni, Reviews of Modern Physics 82, 451 (2010)
2010
-
[18]
Capozziello and M
S. Capozziello and M. De Laurentis, Physics Reports 509, 167 (2011)
2011
-
[19]
Nojiri and S
S. Nojiri and S. D. Odintsov, Physics Reports505, 59 (2011)
2011
-
[20]
Harada, Phys
J. Harada, Phys. Rev. D101, 024053 (2020)
2020
-
[21]
Teleparallel gravity: An introduction,
R. Aldrovandi and J. G. Pereira, “Teleparallel gravity: An introduction,” (2013)
2013
-
[22]
J. W. Maluf, Annalen der Physik525, 339 (2013)
2013
-
[24]
Zhao, Eur
D. Zhao, Eur. Phys. J. C82, 303 (2022)
2022
-
[25]
W. Wang, H. Chen, and T. Katsuragawa, Phys. Rev. D 105, 024060 (2022)
2022
-
[26]
J. A. Najera, C. A. Alvarado, and C. Escamilla-Rivera, Mon. Not. Roy. Astron. Soc.524, 5280 (2023)
2023
-
[27]
G. N. Gadbail and P. K. Sahoo, Chin. J. Phys.89, 1754 (2024)
2024
-
[28]
Sokoliuk, S
O. Sokoliuk, S. Arora, S. Praharaj, A. Baransky, and P. K. Sahoo, Mon. Not. Roy. Astron. Soc.522, 252 (2023)
2023
-
[29]
F. K. Anagnostopoulos, V. Gakis, E. N. Saridakis, and S. Basilakos, Eur. Phys. J. C83, 58 (2023)
2023
-
[30]
Hinderer, Astrophys
T. Hinderer, Astrophys. J.677, 1216 (2008)
2008
-
[31]
Hinderer, B
T. Hinderer, B. D. Lackey, R. N. Lang, and J. S. Read, Phys. Rev. D81, 123016 (2010)
2010
-
[32]
Andersson and K
N. Andersson and K. D. Kokkotas, Mon. Not. R. Astron. Soc.299, 1059 (1998)
1998
-
[33]
Sotani, K
H. Sotani, K. D. Kokkotas, and N. Stergioulas, Mon. Not. R. Astron. Soc.375, 261 (2007)
2007
-
[34]
J. M. Lattimer and M. Prakash, Astrophys. J.550, 426 (2001)
2001
-
[35]
¨Ozel and P
F. ¨Ozel and P. Freire, Annu. Rev. Astron. Astrophys.54, 401 (2016)
2016
-
[36]
P. B. Demorest, T. Pennucci, S. M. Ransom, M. S. E. Roberts, and J. W. T. Hessels, Nature467, 1081 (2010)
2010
-
[37]
J. e. a. Antoniadis, Science340, 1233232 (2013)
2013
-
[38]
H. T. e. a. Cromartie, Nature Astronomy4, 72 (2020)
2020
-
[39]
E. e. a. Fonseca, Astrophys. Jour. Lett.915, L12 (2021)
2021
-
[40]
K. C. Gendreau and Z. Arzoumanian, Nature Astronomy 1, 895 (2017)
2017
-
[41]
M. C. e. a. Miller, Astrophys. Jour. Lett.918, L28 (2021)
2021
-
[42]
T. E. e. a. Riley, Astrophys. Jour. Lett.918, L27 (2021)
2021
-
[43]
Karan, A
A. Karan, A. Kumar, M. Sinha, and R. Mallick, Monthly Notices of the Royal Astronomical Society543, 83–94 (2025)
2025
-
[44]
B. P. e. a. Abbott, Phys. Rev. Lett.121, 161101 (2018)
2018
-
[45]
E. e. a. Berti, Class. and Quant. Grav.32, 243001 (2015)
2015
-
[46]
S. e. a. Capozziello, Phys. Rev. D93, 023501 (2016)
2016
-
[47]
K. V. e. a. Staykov, Phys. Rev. D90, 023005 (2014)
2014
-
[48]
P. N. Moreno, A. Wojnar, and F. J. Llanes-Estrada, Journal of Cosmology and Astroparticle Physics2025, 015 (2025)
2025
-
[49]
Cheoun, C
M.-K. Cheoun, C. Deliduman, C. G¨ ung¨ or, V. Kele¸ s, C. Ryu, T. Kajino, and G. J. Mathews, Journal of Cosmology and Astroparticle Physics2013, 021–021 (2013)
2013
-
[50]
A. V. Astashenok, S. Capozziello, and S. D. Odintsov, Journal of Cosmology and Astroparticle Physics2013, 040–040 (2013)
2013
-
[51]
Hendi, G
S. Hendi, G. Bordbar, B. E. Panah, and S. Panahiyan, Journal of Cosmology and Astroparticle Physics2016, 013–013 (2016). 16
2016
-
[52]
Vilhena, S
S. Vilhena, S. Duarte, M. Dutra, and P. Pompeia, Journal of Cosmology and Astroparticle Physics2023, 044 (2023)
2023
-
[53]
G. e. a. Raaijmakers, Astrophys. Jour. Lett.918, L29 (2021)
2021
-
[54]
R. e. a. Essick, Phys. Rev. D102, 043003 (2020)
2020
-
[55]
Jiang, C
J.-L. Jiang, C. Ecker, and L. Rezzolla, The Astrophys. Jour.949, 11 (2023)
2023
-
[56]
Lim and J
Y. Lim and J. W. Holt, Phys. Rev. Lett.121, 062701 (2018)
2018
-
[57]
D. Roy, T. Malik, S. Bhattacharya, and S. Banik, The Astrophys. Jour.968, 124 (2024)
2024
-
[58]
Nobleson, S
K. Nobleson, S. Banik, and T. Malik, Phys. Rev. D107, 124045 (2023), arXiv:2306.01054 [gr-qc]
2023 arXiv
-
[59]
Biswas, E
B. Biswas, E. Smyrniotis, I. Liodis, and N. Stergioulas, Phys. Rev. D109, 064048 (2024)
2024
-
[60]
e. a. C. J. Xia, Commun. in Theor. Phys.74, 095303 (2022)
2022
-
[61]
G. A. Lalazissis, T. Niksic, D. Vretenar, and P. Ring, Phys. Rev. C71, 024312 (2005)
2005
-
[62]
J. B. Jimenez, L. Heisenberg, and T. S. Koivisto, J. Cosmol. Astropart. Phys.2018, 039 (2018)
2018
-
[63]
J. B. Jim´ enez, L. Heisenberg, and T. S. Koivisto, Phys. Rev. D98, 044019 (2018)
2018
-
[64]
J. B. Jimenez, L. Heisenberg, T. Koivisto, and S. Pekar, Phys. Rev. D101, 103507 (2020)
2020
-
[65]
J. F. Donoghue, Phys. Rev. D50, 3874 (1994)
1994
-
[66]
Khyllepet al., Phys
W. Khyllepet al., Phys. Rev. D107, 044022 (2023)
2023
-
[67]
T. E. Rileyet al., Astrophys. J. Lett.887, L21 (2019)
2019
-
[68]
https://zenodo.org/records/8239000
-
[69]
M. C. Milleret al., Astrophys. J. Lett.887, L24 (2019)
2019
-
[70]
https://zenodo.org/record/3473466
-
[71]
https://zenodo.org/records/4697625
-
[72]
https://zenodo.org/records/4670689
-
[73]
Choudhuryet al., Astrophys
D. Choudhuryet al., Astrophys. J. Lett.971, L20 (2024)
2024
-
[74]
https://zenodo.org/records/13766753
-
[75]
Mauviardet al., Astrophys
L. Mauviardet al., Astrophys. J.995, 60 (2025)
2025
-
[76]
https://zenodo.org/records/15603406
-
[77]
J. R. Oppenheimer and G. M. Volkoff, Phys. Rev.55, 374 (1939)
1939
-
[78]
R. C. Tolman, Phys. Rev.55, 364 (1939)
1939
-
[79]
S. M. A. Imam and N. K. Patra, Phys. Rev. D112, 103018 (2025), arXiv:2509.07109 [nucl-th]
2025
-
[80]
R. E. Kass and A. E. Raftery, J. Am. Stat. Assoc.90, 773 (1995)
1995
-
[81]
Jarosz and J
A. Jarosz and J. Wiley, The Journal of Problem Solving 7(2014), 10.7771/1932-6246.1167
2014
Reviewed August 3, 2026 · model on record in the stance chip above.
Discussion (0). Sign in to comment.