REVIEW 4 major objections 4 minor 1 cited by
Equation of State Independent Determination on the Radius of a 1.4 $M_{\odot}$ Neutron Star Using Mass-Radius Measurements
T0 review · 4 major / 4 minor · reviewed 2026-08-11 · deepseek-v4-flash
Pith's one-line read This paper claims that the radius of a 1.4-solar-mass neutron star can be inferred directly from NICER mass-radius measurements of two pulsars, without choosing a specific dense-matter equation of state, and that the value is about 12 km.
desk verdict A simple, transparent data-driven recipe for R1.4 from two NICER posteriors, but the 'EOS-independent' title overstates the assumptions. 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 linear interpolation formula $R_{1.4} = R_1 + \frac{R_2-R_1}{M_2-M_1} (1.4\,M_\odot - M_1)$, applied to posterior samples drawn from each pulsar's mass-radius measurement after excluding pairs with equal masses (which would imply twin stars or a phase transition). In the second scenario, the machinery is the product of the two radius posterior distributions, each conditioned on a mass of exactly $1.4\,M_\odot$, so that enforcing a shared radius combines the information from both stars. The interpolation converts two nearby mass-radius measurements into an estimate at the canonical mass, while the Bayesian product tightens the estimate under the no-twin-star assumption.
What would settle it
A concrete test: simulate mass-radius observations from an equation of state with a strong phase transition near $1.4\,M_\odot$, generate two sources bracketing that mass with 5% uncertainties, and apply the interpolation; if the recovered $R_{1.4}$ is biased by more than the reported uncertainty, the smoothness assumption fails.
Extended reading notes
Core claim
The central claim is that $R_{1.4}$ is about 12 km, specifically $11.99^{+1.62}_{-1.57}$ km under the first scenario with the PDT-U hotspot model and $12.09^{+0.81}_{-0.63}$ km under the second scenario with ST+PDT, and that these estimates are largely independent of the equation of state. The paper argues that because both pulsars have masses near $1.4\,M_\odot$, a straight-line interpolation between a lower-mass sample and an upper-mass sample in the mass-radius plane can reach the canonical mass, and in the equal-mass scenario the two independent radius posteriors can be multiplied to sharpen the answer. The author frames this as a data-driven pathway that naturally absorbs EOS model systematic errors into the reported uncertainty, with the residual EOS dependence coming only from assumptions such as the oblate surface shape in pulse-profile modeling and the exclusion of phase transitions.
Load-bearing premise
The load-bearing premise is that the true mass-radius relation between the two pulsars is smooth and monotonic enough that a straight line between any bracketing pair passes through the radius at $1.4\,M_\odot$, which rules out phase transitions and twin stars; a second premise is that the NICER posteriors are faithful despite the equation-of-state-dependent oblate-surface assumption baked into the pulse-profile modeling.
Editorial extensions
If this is right
- The reported $R_{1.4}$ values cluster near 12 km, and the corresponding $ \Lambda_{1.4}$ values are consistent with the GW170817 tidal-deformability measurement.
- Under Scenario 2, assuming both stars are exactly $1.4\,M_\odot$, the radius posterior is narrower than the interpolation-based Scenario 1, with the ST+PDT hotspot configuration giving the tightest 95% credible interval.
- Simulated future observations with about 5% mass-radius precision recover an injected $R_{1.4}$ to within roughly 0.02 km in the two equations of state tested, indicating the method is accurate for ideal two-source data.
- Because the method avoids explicit EOS priors, its uncertainties are argued to absorb EOS model systematics that traditional EOS-based inference treats as a separate step.
Reading between the lines
- Inference: the pairwise interpolation logic generalizes to any set of pulsars whose masses bracket $1.4\,M_\odot$; grouping sources into pairs and multiplying the resulting distributions could tighten $R_{1.4}$ as more near-canonical pulsars are measured.
- Inference: the method's key assumption can be stress-tested by simulating mass-radius posteriors from an equation of state with a first-order phase transition at or near $1.4\,M_\odot$; a linear interpolation should then recover a biased $R_{1.4}$, quantifying the size of the model error.
- Inference: because the NICER posteriors themselves carry some EOS dependence through the oblate-surface assumption, a fully model-free radius will require independent constraints on the stellar shape from future observations.
- Inference: practical application to real data with current uncertainties will likely be dominated by hotspot-model systematics rather than the interpolation itself, since the three J0030 hotspot configurations shift the Scenario 1 central values by about 0.3 km.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper proposes two data-driven methods to infer the radius of a 1.4 solar-mass neutron star (R1.4) from NICER mass-radius posteriors of PSR J0030+0451 (J0030) and PSR J0437-4715 (J0437), without explicitly constructing an equation of state (EOS). Scenario 1 samples pairs of (M,R) from kernel density estimates of the two sources' posteriors, keeps pairs that bracket 1.4 solar masses, and linearly interpolates to obtain an R1.4 distribution. Scenario 2 assumes both stars have mass 1.4 solar masses and the same radius, then combines conditional radius posteriors via a product formula that also includes the probability each star has mass exactly 1.4. The paper reports R1.4 values for several J0030 hotspot models (Table 2), compares with astrophysical and nuclear constraints, and validates Scenario 1 on simulated future X-ray data generated from two relativistic mean-field EOSs (NL3 omega-rho and TM1-2omega-rho).
Significance. If the method is valid, it offers a relatively model-independent way to extract a key neutron-star observable, complementing EOS-inference approaches and providing a check on nuclear-physics constraints. The paper is transparent about several assumptions, and it uses open-source software and public likelihoods, which aids reproducibility. However, the claimed 'equation-of-state independence' is overstated: the input NICER posteriors carry EOS assumptions (oblate shape), and both scenarios explicitly exclude phase transitions and twin-star behavior. The method's novelty is therefore more modest than the title suggests, but the idea of interpolating or conditioning on mass-radius posteriors from nearby-mass sources is useful for future X-ray missions.
major comments (4)
- [Title and Abstract; Section 2.1] The claim of an 'Equation of State Independent' determination is contradicted by the paper's own admissions: Section 2.1 states that NICER inference results are 'not entirely EOS independent' because of the oblate surface assumption, and both scenarios explicitly exclude phase transitions and twin stars (Scenario 1 assumes linear interpolation; Scenario 2 assumes equal mass and equal radius). The title and abstract should be qualified, e.g., 'largely EOS-model-independent' or 'EOS-model-independent given the smoothness assumption,' so that readers are not misled about the strength of the result.
- [Section 2.1, Scenario 1 (Equation 2)] The sampling and interpolation procedure is under-specified: the paper does not report the KDE bandwidth choices, the number (or fraction) of accepted sample pairs that satisfy the bracketing condition M1<1.4<M2 or M2<1.4<M1, or any convergence checks. Because J0437's mass is 1.44±0.07 Msun and J0030's mass posterior is broad and model-dependent, the accepted pairs may be a small subset of the posterior, and the resulting R1.4 distribution in Table 2 could be dominated by the tails of the KDE rather than by the bulk of the mass-radius information. The authors should report the acceptance fraction and demonstrate robustness to KDE choices.
- [Section 2.1, Equation (6)] The joint posterior formula for Scenario 2 is not derived from a clear probabilistic model. The product of the two conditional posteriors P(Ri|Oi,Mi=1.4) is appropriate if the two stars share a common radius, but the extra factors P(Mi=1.4|Oi) are posterior probabilities of a point mass value and require an explicit prior on the mass. Their inclusion is not justified and the resulting expression is not a proper posterior under Bayes' theorem. The authors should either derive Equation (6) from a well-defined hierarchical model or remove the P(Mi|Oi) factors and explain why the conditional posteriors alone suffice.
- [Section 4, Simulated dataset] The simulation study validates Scenario 1 only for two smooth RMF EOSs with no phase transitions, which is exactly the assumption the method is built on. It does not test the failure mode—for example, a first-order phase transition between the masses of the two simulated stars, or a twin-star configuration—so the claim that the method 'is capable of recovering the underlying R1.4' is only demonstrated in favorable, assumption-satisfying cases. The paper should either run a phase-transition EOS in the simulation or explicitly state that the validity is conditional on the smoothness assumption and that no test of the failure mode is provided.
minor comments (4)
- [Section 2.1, paragraph on equal-mass exclusion] The text says 'we exclude pairs where M1 = M2' and then says 'we test the condition M1≠M2 with a precision of 10^-3 Msun.' These two statements are inconsistent; the paper should clarify that exact equality is replaced by a threshold of 10^-3 Msun.
- [Section 3.2, Table 1 and discussion] The phrase 'the PDT-U model most closely represents a 1.4 Msun star' is imprecise; what is meant is that the PDT-U mass posterior places the most probability near 1.4 Msun. Please rephrase to avoid implying the model itself is a physical object.
- [Section 2.2, data description] The ST+PST configuration is described as using the mass-radius measurement from Riley et al. (2019), while later the paper says the updated ST+PST inference from Vinciguerra et al. (2024) is used. Please clarify which posterior samples are actually used for each hotspot model.
- [References] Several citations are to 'in prep' works (Huang 2024; Huang & Chen 2024) and to arXiv preprints. For a journal submission, these should be updated to published versions or clearly marked as unpublished.
Circularity Check
No circularity: R1.4 is a conditional/interpolated summary of the input NICER posteriors, not an input to the fits.
full rationale
The derivation chain in this paper is a statistical re-expression of the NICER mass-radius posteriors, not a circular self-confirmation. In Scenario 1, Eq. (2) linearly interpolates between bracketing (M, R) samples drawn from the KDE-reconstructed posteriors of J0030 and J0437; the target R1.4 is not an input to any fit or likelihood. In Scenario 2, Eqs. (4)-(6) condition the same input posteriors at M = 1.4 M_sun and combine the conditional radius posteriors under the explicit assumption R1 = R2 = R1.4. The quoted values in Table 2 are therefore direct summaries of the input posterior samples, not independent predictions that were fitted to R1.4 itself. The paper also acknowledges the main limitations that prevent the result from being fully EOS-independent: the NICER posterior inputs are 'not entirely EOS independent' because of the oblate-surface assumption based on a preselected set of EOSs (Section 2.1), and Scenario 1 assumes a smooth linear mass-radius relation with no phase transitions. These are correctness/limitation concerns, not circularity. The self-citation to the author's CompactObject package (Huang et al. 2024c) is a software implementation citation and is not load-bearing for the statistical argument; the in-preparation citations (Huang 2024; Huang & Chen 2024) are contextual, not part of the derivation. No circular step can be exhibited in which a prediction reduces by construction to its own input, and no load-bearing uniqueness claim is imported from the author's prior work.
Assumptions & free parameters
free parameters (4)
- KDE bandwidth =
Not specified (Python default)
- Equal-mass exclusion threshold =
10^-3 Msun
- Uniform prior range in Scenario 2 =
6 to 16 km
- Simulated uncertainty level =
5% (68% range)
assumptions (5)
- domain assumption The NICER mass-radius posterior samples for J0030 and J0437 are faithful representations of the true mass-radius measurements.
- ad hoc to paper The true mass-radius relation is smooth and monotonic between the masses of J0030 and J0437, so linear interpolation is valid.
- domain assumption No twin stars exist, i.e., two stars with equal mass have equal radius.
- domain assumption The empirical tidal deformability relation Lambda(R1.4) = 2.88e-6 (R1.4/km)^7.5 holds for the inferred stars.
- standard math The KDE reconstruction accurately represents the posterior, including tails.
Cite this review
Pith. "Pith review of Equation of State Independent Determination on the Radius of a 1.4 $M_{\odot}$ Neutron Star Using Mass-Radius Measurements." pith.science (2026). https://pith.science/paper/ASXIRCSR
@misc{pith2026241210242,
author = {Pith},
title = {Pith review of: Equation of State Independent Determination on the Radius of a 1.4 $M_\odot$ Neutron Star Using Mass-Radius Measurements},
year = {2026},
howpublished = {\url{https://pith.science/paper/ASXIRCSR}},
note = {Machine review of arXiv:2412.10242}
}
abstract
Traditional methods for determining the radius of a 1.4 $M_{\odot}$ neutron star ($R_{1.4}$) rely on specific equations of state (EOS) models that describe various types of dense nuclear matter. This dependence on EOS models can introduce substantial systematic uncertainties, which may exceed the measurement uncertainties when constraining $R_{1.4}$. In this study, we explore a novel approach to constraining $R_{1.4}$ using data from NICER observations of PSR J0030+0451 (J0030) and PSR J0437-4715 (J0437). However, this work presents a more data-driven analysis framework, substantially decreasing the need for EOS assumptions. By analyzing the Mass-Radius measurements of these two neutron stars, we infer $R_{1.4}$ using statistical methods based mostly on observational data. We examine various hotspot configurations for J0030, along with new J0437 observations, and their effects on the inferred radius. Our results are consistent with X-ray timing, gravitational wave, and nuclear physics constraints, while avoiding EOS-related biases. The same method has also been applied to a simulated mass-radius dataset, based on our knowledge of future X-ray telescopes, demonstrating the model's ability to recover the injected $R_{1.4}$ value in certain cases. This method provides a data-driven pathway for extracting neutron star properties and offers a new approach for future observational efforts in neutron star astrophysics.
Figures
Forward citations
Cited by 1 Pith paper
-
Effects of dark matter and magnetic field on neutron star properties in relativistic mean-field theory: A single-fluid approach
Fermionic dark matter and strong central magnetic fields both reduce neutron-star maximum mass and radius and lower tidal deformability in single-fluid RMF models, remaining compatible with GW/NICER constraints over t...
Reference graph
Works this paper leans on
-
[1]
2018, Physical Review Letters, 121, doi: 10.1103/physrevlett.121.161101
Abbott, B., Abbott, R., Abbott, T., et al. 2018, Physical Review Letters, 121, doi: 10.1103/physrevlett.121.161101
-
[2]
Abbott, B. P., Abbott, R., Abbott, T. D., et al. 2017, Physical Review Letters, 119, doi: 10.1103/physrevlett.119.161101 —. 2019, Phys. Rev. X, 9, 011001, doi: 10.1103/PhysRevX.9.011001
-
[4]
Adhikari, D., et al. 2021b, Phys. Rev. Lett., 126, 172502, doi: 10.1103/PhysRevLett.126.172502
-
[5]
Adhikari, D., Albataineh, H., Androic, D., et al. 2022, Phys. Rev. Lett., 129, 042501, doi: 10.1103/PhysRevLett.129.042501
-
[6]
2018, Physical Review Letters, 120, doi: 10.1103/physrevlett.120.172703
Annala, E., Gorda, T., Kurkela, A., & Vuorinen, A. 2018, Physical Review Letters, 120, doi: 10.1103/physrevlett.120.172703
-
[7]
Barret, D., Nandra, K., Barcons, X., et al. 2013, in SF2A-2013: Proceedings of the Annual meeting of the French Society of Astronomy and Astrophysics, ed. L. Cambresy, F. Martins, E. Nuss, & A. Palacios, 447–453, doi: 10.48550/arXiv.1310.3814
-
[8]
Bogdanov, S., Lamb, F. K., Mahmoodifar, S., et al. 2019, ApJL, 887, L26, doi: 10.3847/2041-8213/ab5968
-
[9]
2021, The Journal of Open Source Software, 6, 3001, doi: 10.21105/joss.03001
Buchner, J. 2021, The Journal of Open Source Software, 6, 3001, doi: 10.21105/joss.03001
Show all 45 references
-
[10]
2024, The Astrophysical Journal Letters, 971, L20, doi: 10.3847/2041-8213/ad5a6f
Choudhury, D., Salmi, T., Vinciguerra, S., et al. 2024, The Astrophysical Journal Letters, 971, L20, doi: 10.3847/2041-8213/ad5a6f
2024 doi
-
[11]
J., Miller, M
Dittmann, A. J., Miller, M. C., Lamb, F. K., et al. 2024, ApJ, 974, 295, doi: 10.3847/1538-4357/ad5f1e
2024 doi
-
[12]
M., et al
Drischler, C., Han, S., Lattimer, J. M., et al. 2021, PhRvC, 103, 045808, doi: 10.1103/PhysRevC.103.045808
2021 doi
-
[13]
2021, PhRvL, 127, 192701, doi: 10.1103/PhysRevLett.127.192701
Essick, R., Tews, I., Landry, P., & Schwenk, A. 2021, PhRvL, 127, 192701, doi: 10.1103/PhysRevLett.127.192701
2021 doi
-
[14]
T., Pennucci, T
Fonseca, E., Cromartie, H. T., Pennucci, T. T., et al. 2021, ApJL, 915, L12, doi: 10.3847/2041-8213/ac03b8
2021 doi
-
[15]
S., Roming, P., Ray, P., et al
Froning, C. S., Roming, P., Ray, P., et al. 2024, in Space Telescopes and Instrumentation 2024: Ultraviolet to Gamma Ray, Vol. 13093, SPIE, 729–735
2024
-
[16]
C., Arzoumanian, Z., Adkins, P
Gendreau, K. C., Arzoumanian, Z., Adkins, P. W., et al. 2016, in Space telescopes and instrumentation 2016: Ultraviolet to gamma ray, Vol. 9905, SPIE, 420–435
2016
-
[17]
2024, in prep
Huang, C., & Chen, A. 2024, in prep
2024
-
[18]
L., Tolos, L., & Providˆ encia, C
Huang, C., Raaijmakers, G., Watts, A. L., Tolos, L., & Providˆ encia, C. 2024a, Mon. Not. Roy. Astron. Soc., 529, 4650, doi: 10.1093/mnras/stae844
-
[19]
Huang, C., Tolos, L., Providˆ encia, C., & Watts, A. 2024b. https://arxiv.org/abs/2410.14572
-
[20]
2024c, Journal of Open Source Software
Huang, C., Malik, T., Cartaxo, J., et al. 2024c, Journal of Open Source Software. https://arxiv.org/abs/2411.14615
-
[21]
2023, Phys
Kumar, M., Kumar, S., Thakur, V., et al. 2023, Phys. Rev. C, 107, 055801, doi: 10.1103/PhysRevC.107.055801
2023 doi
-
[22]
M., & Lim, Y
Lattimer, J. M., & Lim, Y. 2013, ApJ, 771, 51, doi: 10.1088/0004-637X/771/1/51
2013 doi
- [23]
-
[24]
2024, PhRvC, 109, 035801, doi: 10.1103/PhysRevC.109.035801
Lim, Y., & Schwenk, A. 2024, PhRvC, 109, 035801, doi: 10.1103/PhysRevC.109.035801
2024 doi
-
[25]
C., Lamb, F
Miller, M. C., Lamb, F. K., Dittmann, A. J., et al. 2019, The Astrophysical Journal Letters, 887, L24, doi: 10.3847/2041-8213/ab50c5
2019 doi
-
[27]
C., Lamb, F
Miller, M. C., Lamb, F. K., Dittmann, A. J., et al. 2021, The Astrophysical Journal Letters, 918, L28, doi: 10.3847/2041-8213/ac089b
2021 doi
-
[28]
2023, Phys
Mondal, C., & Gulminelli, F. 2023, Phys. Rev. C, 107, 015801, doi: 10.1103/PhysRevC.107.015801 Providˆ encia, C. m. c., & Rabhi, A. 2013, Phys. Rev. C, 87, 055801, doi: 10.1103/PhysRevC.87.055801
2023 doi
-
[29]
E., Watts, A
Raaijmakers, G., Riley, T. E., Watts, A. L., et al. 2019, The Astrophysical Journal Letters, 887, L22, doi: 10.3847/2041-8213/ab451a
2019 doi
-
[30]
K., Riley, T
Raaijmakers, G., Greif, S. K., Riley, T. E., et al. 2020, The Astrophysical Journal Letters, 893, L21, doi: 10.3847/2041-8213/ab822f
2020 doi
-
[31]
K., Hebeler, K., et al
Raaijmakers, G., Greif, S. K., Hebeler, K., et al. 2021, The Astrophysical Journal Letters, 918, L29, doi: 10.3847/2041-8213/ac089a
2021 doi
-
[32]
S., Arzoumanian, Z., Ballantyne, D., et al
Ray, P. S., Arzoumanian, Z., Ballantyne, D., et al. 2019, arXiv e-prints, arXiv:1903.03035. https://arxiv.org/abs/1903.03035
2019 arXiv
-
[33]
J., Bailes, M., Shannon, R
Reardon, D. J., Bailes, M., Shannon, R. M., et al. 2024, The Astrophysical Journal Letters, 971, L18, doi: 10.3847/2041-8213/ad614a
2024 doi
-
[34]
T., Fattoyev, F., Horowitz, C., & Piekarewicz, J
Reed, B. T., Fattoyev, F., Horowitz, C., & Piekarewicz, J. 2021, Physical Review Letters, 126, doi: 10.1103/physrevlett.126.172503
2021 doi
-
[35]
2013, Phys
Reinhard, P.-G., Piekarewicz, J., Nazarewicz, W., et al. 2013, Phys. Rev. C, 88, 034325, doi: 10.1103/PhysRevC.88.034325
2013 doi
-
[36]
2021, Phys
Reinhard, P.-G., Roca-Maza, X., & Nazarewicz, W. 2021, Phys. Rev. Lett., 127, 232501, doi: 10.1103/PhysRevLett.127.232501 —. 2022, Phys. Rev. Lett., 129, 232501, doi: 10.1103/PhysRevLett.129.232501
2021 doi
-
[37]
E., Watts, A
Riley, T. E., Watts, A. L., Bogdanov, S., et al. 2019, The Astrophysical Journal Letters, 887, L21, doi: 10.3847/2041-8213/ab481c
2019 doi
-
[38]
E., Watts, A
Riley, T. E., Watts, A. L., Ray, P. S., et al. 2021, ApJL, 918, L27, doi: 10.3847/2041-8213/ac0a81
2021 doi
-
[39]
2024, The Astrophysical Journal Letters, 971, L19, doi: 10.3847/2041-8213/ad5f02
Rutherford, N., Mendes, M., Svensson, I., et al. 2024, The Astrophysical Journal Letters, 971, L19, doi: 10.3847/2041-8213/ad5f02
2024 doi
-
[40]
2022, ApJ, 941, 150, doi: 10.3847/1538-4357/ac983d 11
Salmi, T., Vinciguerra, S., Choudhury, D., et al. 2022, ApJ, 941, 150, doi: 10.3847/1538-4357/ac983d 11
2022 doi
- [41]
-
[42]
2024, The Radius of the High Mass Pulsar PSR J0740+6620 With 3.6 Years of NICER Data
Salmi, T., Choudhury, D., Kini, Y., et al. 2024, The Radius of the High Mass Pulsar PSR J0740+6620 With 3.6 Years of NICER Data. https://arxiv.org/abs/2406.14466
2024 arXiv
-
[43]
2023, Physical Review C, 107, doi: 10.1103/physrevc.107.045806
Scurto, L., Pais, H., & Gulminelli, F. 2023, Physical Review C, 107, doi: 10.1103/physrevc.107.045806
2023 doi
-
[44]
L., et al
Vinciguerra, S., Salmi, T., Watts, A. L., et al. 2024, The Astrophysical Journal, 961, 62, doi: 10.3847/1538-4357/acfb83
2024 doi
-
[45]
L., Andersson, N., Chakrabarty, D., et al
Watts, A. L., Andersson, N., Chakrabarty, D., et al. 2016, Reviews of Modern Physics, 88, 021001, doi: 10.1103/RevModPhys.88.021001
2016 doi
-
[46]
L., Yu, W., Poutanen, J., et al
Watts, A. L., Yu, W., Poutanen, J., et al. 2019, Science China Physics, Mechanics, and Astronomy, 62, 29503, doi: 10.1007/s11433-017-9188-4
2019 doi
-
[47]
2022, Phys
Yue, T.-G., Chen, L.-W., Zhang, Z., & Zhou, Y. 2022, Phys. Rev. Res., 4, L022054, doi: 10.1103/PhysRevResearch.4.L022054
2022 doi
Reviewed August 11, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.