Recognition: unknown
Combining the Mass--Radius Posteriors of J0030+0451 Allowing for Unknown Model Systematics
Pith reviewed 2026-05-10 15:11 UTC · model grok-4.3
The pith
A Bayesian framework combines eight published mass-radius posteriors for PSR J0030+0451 while allowing for unknown model systematics.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
We adapt a robust Bayesian combination framework to the published M–R posteriors of PSR J0030+0451 while allowing for unknown systematic uncertainties that might have led to the apparently divergent results. Using this technique, we combine eight existing M–R posteriors into a single conservative and reproducible posterior that incorporates unknown model systematics across the currently available analyses and is suitable for direct use in EoS studies. The resulting constraint is M = 1.46^{+0.09}_{-0.08} M_⊙, R = 12.69^{+0.64}_{-0.55} km, and C = 0.172^{+0.006}_{-0.007} (68% credible interval).
What carries the argument
Bayesian combination framework that introduces a set of unknown systematic parameters whose prior is chosen without reference to any specific equation of state to capture discrepancies between published posteriors.
If this is right
- The combined posterior can be used directly in equation-of-state inference without choosing one analysis over the others.
- Joint analysis with PSR J0437-4715 and GW170817 yields R_{1.4} = 11.98^{+0.58}_{-0.68} km and Lambda_{1.4} = 320^{+216}_{-138}.
- The method supplies a reproducible way to fold cross-model uncertainty into future neutron-star equation-of-state pipelines.
- The reported mass, radius, and compactness values at 68% credible interval are the direct output of the combination.
Where Pith is reading between the lines
- The same combination procedure could be applied to other NICER pulsars that have multiple independent analyses.
- If additional hotspot models are published later, they can be folded into an updated combined posterior using the identical framework.
- The approach separates analysis-specific modeling choices from the underlying observational data, which may reduce the impact of any single modeling decision on EoS constraints.
Load-bearing premise
The discrepancies between the eight published posteriors can be adequately captured by a single set of unknown systematic parameters whose prior does not depend on the target equation of state.
What would settle it
Apply the same combination procedure to simulated posteriors generated from a known true mass-radius value with injected model differences; if the output posterior fails to contain the true value at the stated credible-interval rate, the method is miscalibrated.
Figures
read the original abstract
The NASA Neutron star Interior Composition Explorer (\emph{NICER}) mission measures the X-ray pulse profiles of select millisecond pulsars and uses sophisticated pulse profile modeling (PPM) techniques to constrain their masses ($M$) and radii ($R$), in order to probe the state of matter in their interiors. One of the most studied pulsars, PSR J0030+0451, has been analyzed by multiple groups using different choices of hotspot models. The different choices of hotspot prescriptions to fit the same observational data led to different $M$--$R$ posteriors that do not completely agree with one another, resulting in a practical bottleneck for dense-matter equation-of-state (EoS) inference. In this paper, we adapt a robust Bayesian combination framework to the published $M$--$R$ posteriors of PSR J0030+0451 while allowing for unknown systematic uncertainties that might have led to the apparently divergent results. Using this technique, we combine eight existing $M$--$R$ posteriors into a single conservative and reproducible posterior that incorporates unknown model systematics across the currently available analyses and is suitable for direct use in EoS studies. The resulting constraint is $M = 1.46^{+0.09}_{-0.08}\,M_\odot$, $R = 12.69^{+0.64}_{-0.55}\,\mathrm{km}$, and compactness $C = 0.172^{+0.006}_{-0.007}$ (68\% credible interval). Incorporating this combined J0030+0451 constraint in an EoS-agnostic joint analysis with PSR~J0437--4715 and GW170817 yields $R_{1.4} = 11.98^{+0.58}_{-0.68}\,\mathrm{km}$ and $\Lambda_{1.4} = 320^{+216}_{-138}$. Our results provide a combined $M$--$R$ constraint for J0030+0451 and a practical framework for incorporating cross-model uncertainty into neutron star EoS inference pipelines.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript adapts a Bayesian combination framework to merge eight published M–R posteriors for PSR J0030+0451 (from independent NICER PPM analyses using different hotspot models), explicitly marginalizing over unknown model systematics via additional parameters with an EoS-independent prior. The resulting combined posterior is reported as M = 1.46^{+0.09}_{-0.08} M_⊙, R = 12.69^{+0.64}_{-0.55} km, C = 0.172^{+0.006}_{-0.007} (68% CI) and is then used in an EoS-agnostic joint analysis with J0437–4715 and GW170817 to obtain R_{1.4} = 11.98^{+0.58}_{-0.68} km and Λ_{1.4} = 320^{+216}_{-138}. The central claim is that this yields a conservative, reproducible constraint suitable for direct EoS inference.
Significance. If validated, the work supplies a practical, reproducible way to handle discrepant M–R posteriors arising from discrete modeling choices in NICER analyses, which is a recurring bottleneck for neutron-star EoS studies. Explicitly parameterizing and marginalizing unknown systematics (rather than ad-hoc selection or simple averaging) is a constructive contribution, and releasing the combined posterior enhances immediate usability by the community.
major comments (2)
- [Methods (Bayesian combination)] Methods section on the Bayesian combination framework: no simulation-based calibration is presented in which synthetic pulse-profile data are generated from a known true M–R, refitted with the same eight hotspot models, and the recovered combined posterior checked for coverage, bias, or over-/under-conservatism. This test is load-bearing for the claim that the output posterior is conservative and properly calibrated when the true systematics are unknown.
- [Methods (systematic parameterization)] Description of the systematic-parameter prior (likely §3): the framework assumes that discrepancies among the eight input posteriors can be captured by a single set of unknown systematic parameters whose prior is chosen independently of the target EoS. If the actual modeling differences are correlated with the data in ways not spanned by this parameterization, the combined posterior will be miscalibrated; the manuscript provides no robustness checks or alternative prior specifications to quantify this risk.
minor comments (1)
- [Abstract] The abstract states the numerical results but does not explicitly note that eight posteriors are being combined until the body; adding this detail to the abstract would improve immediate clarity for readers.
Simulated Author's Rebuttal
We thank the referee for their thoughtful and constructive comments. We respond to each major comment below with the strongest honest defense of the manuscript, indicating where revisions have been made.
read point-by-point responses
-
Referee: [Methods (Bayesian combination)] Methods section on the Bayesian combination framework: no simulation-based calibration is presented in which synthetic pulse-profile data are generated from a known true M–R, refitted with the same eight hotspot models, and the recovered combined posterior checked for coverage, bias, or over-/under-conservatism. This test is load-bearing for the claim that the output posterior is conservative and properly calibrated when the true systematics are unknown.
Authors: We agree that an end-to-end simulation-based calibration would be the most direct way to verify coverage and calibration properties. However, performing this test requires generating synthetic NICER pulse-profile data under each of the eight distinct hotspot models and re-executing the full PPM analyses, which is computationally prohibitive and depends on proprietary details of the original pipelines and data reduction steps that are not accessible to us. The framework we adapt produces conservative posteriors by design through explicit marginalization over unknown systematics with a deliberately broad, EoS-independent prior; this construction prevents the combined posterior from being narrower than warranted by the input discrepancies. In the revised manuscript we have added a dedicated paragraph in the Methods section that explains this rationale, cites the conservative properties of the prior, and notes the practical barriers to full synthetic calibration as a limitation for future work. revision: partial
-
Referee: [Methods (systematic parameterization)] Description of the systematic-parameter prior (likely §3): the framework assumes that discrepancies among the eight input posteriors can be captured by a single set of unknown systematic parameters whose prior is chosen independently of the target EoS. If the actual modeling differences are correlated with the data in ways not spanned by this parameterization, the combined posterior will be miscalibrated; the manuscript provides no robustness checks or alternative prior specifications to quantify this risk.
Authors: The EoS-independent prior on the systematic parameters was chosen precisely to avoid injecting model-dependent bias into the combination. In the revised manuscript we have added an appendix containing explicit robustness checks: we re-ran the combination under alternative prior widths, different hyperparameter values, and a hierarchical variant of the systematic model. These tests show that the resulting M, R, and C posteriors for J0030+0451 shift by amounts well within the quoted 68% credible intervals. We have also expanded the main text to discuss the assumption that the chosen parameterization spans the dominant modeling differences and to note that any residual unmodeled correlations would be buffered by the breadth of the prior. These additions directly quantify the sensitivity the referee correctly identifies. revision: yes
Circularity Check
No significant circularity in the combination framework
full rationale
The paper takes eight independent, previously published M-R posteriors for PSR J0030+0451 as direct inputs and applies an adapted Bayesian combination procedure that introduces additional systematic parameters whose prior is chosen independently of the target EoS. The output posterior (M = 1.46^{+0.09}_{-0.08} M_⊙, R = 12.69^{+0.64}_{-0.55} km) is produced by marginalizing over those systematics; it is not equivalent to any input by definition, nor is any fitted parameter inside the paper renamed as a prediction. No self-definitional equations, fitted-input predictions, or load-bearing self-citations that reduce the central claim to an unverified prior result are present in the derivation chain. The method is therefore self-contained against external benchmarks.
Axiom & Free-Parameter Ledger
Reference graph
Works this paper leans on
-
[1]
Abbott, B. P., et al. 2017, PhRvL, 119, 161101, doi: 10.1103/PhysRevLett.119.161101
-
[2]
Abbott, B. P., Abbott, R., Abbott, T. D., et al. 2018, PhRvL, 121, 161101, doi: 10.1103/PhysRevLett.121.161101
-
[3]
Generic conditions for stable hybrid stars.Phys
Alford, M. G., Han, S., & Prakash, M. 2013, Physical Review D, 88, doi: 10.1103/physrevd.88.083013
-
[4]
AlGendy, M., & Morsink, S. M. 2014, ApJ, 791, 78, doi: 10.1088/0004-637X/791/2/78
-
[5]
Beloborodov, A. M. 2002, ApJL, 566, L85, doi: 10.1086/339511
-
[6]
L., & Peacock, J
Bernal, J. L., & Peacock, J. A. 2018, Journal of Cosmology and Astroparticle Physics, 2018, 002
2018
-
[7]
Bilous, A. V., Watts, A. L., Harding, A. K., et al. 2019, The Astrophysical Journal Letters, 887, L23, doi: 10.3847/2041-8213/ab53e7
-
[8]
Bogdanov, S., & Grindlay, J. E. 2009, ApJ, 703, 1557, doi: 10.1088/0004-637X/703/2/1557
-
[9]
2021, The Journal of Open Source Software, 6, 3001, doi: 10.21105/joss.03001 22Rackham et al
Buchner, J. 2021, The Journal of Open Source Software, 6, 3001, doi: 10.21105/joss.03001
-
[10]
Cadeau, C., Morsink, S. M., Leahy, D., & Campbell, S. S. 2007, ApJ, 654, 458, doi: 10.1086/509103
-
[11]
Y., Yuan, Y., & Vasilopoulos, G
Chen, A. Y., Yuan, Y., & Vasilopoulos, G. 2020, The Astrophysical Journal Letters, 893, L38, doi: 10.3847/2041-8213/ab85c5
-
[12]
Choudhury, D., Salmi, T., Vinciguerra, S., et al. 2024, arXiv e-prints, arXiv:2407.06789, doi: 10.48550/arXiv.2407.06789
-
[13]
Choudhury, D., Watts, A. L., Dittmann, A. J., et al. 2024, The Astrophysical Journal, 975, 202, doi: 10.3847/1538-4357/ad7255
-
[14]
Choudhury, D., Salmi, T., Vinciguerra, S., et al. 2024, Reproduction package for: ‘A NICER View of the Nearest and Brightest Millisecond Pulsar: PSR J0437–4715’, Zenodo, doi: 10.5281/zenodo.13766753
-
[15]
Cromartie, H. T., Fonseca, E., Ransom, S. M., et al. 2019, Nature Astronomy, 4, 72–76, doi: 10.1038/s41550-019-0880-2
-
[16]
Demorest, P. B., Pennucci, T., Ransom, S. M., Roberts, M. S. E., & Hessels, J. W. T. 2010, Nature, 467, 1081–1083, doi: 10.1038/nature09466
-
[17]
Dittmann, A. J., Miller, M. C., Lamb, F. K., et al. 2024, arXiv e-prints, arXiv:2406.14467, doi: 10.48550/arXiv.2406.14467
-
[18]
2017, Nature Astronomy, 1, 895, doi: 10.1038/s41550-017-0301-3
Gendreau, K., & Arzoumanian, Z. 2017, Nature Astronomy, 1, 895, doi: 10.1038/s41550-017-0301-3
-
[19]
Space Telescopes and Instrumentation 2016: Ultraviolet to Gamma Ray , year = 2016, editor =
Gendreau, K. C., Arzoumanian, Z., Adkins, P. W., et al. 2016, in Society of Photo-Optical Instrumentation Engineers (SPIE) Conference Series, Vol. 9905, Space Telescopes and Instrumentation 2016: Ultraviolet to Gamma Ray, ed. J.-W. A. den Herder, T. Takahashi, & M. Bautz, 99051H, doi: 10.1117/12.2231304
-
[20]
Ho, W. C. G., & Lai, D. 2001, MNRAS, 327, 1081, doi: 10.1046/j.1365-8711.2001.04801.x —. 2003, MNRAS, 338, 233, doi: 10.1046/j.1365-8711.2003.06047.x
-
[21]
2002, Monthly Notices of the Royal Astronomical Society, 335, 377
Hobson, M., Bridle, S., & Lahav, O. 2002, Monthly Notices of the Royal Astronomical Society, 335, 377
2002
-
[22]
Huang, C. 2025a, Astrophys. J. Lett., 978, L14, doi: 10.3847/2041-8213/ad9f3c —. 2025b, Astrophys. J., 985, 216, doi: 10.3847/1538-4357/add5ef
-
[23]
Huang, C., & Chen, A. Y. 2025, Astrophys. J., 991, 90, doi: 10.3847/1538-4357/adf747
-
[24]
Huang, C., & Sourav, S. 2025, Astrophys. J., 983, 17, doi: 10.3847/1538-4357/adbb67
-
[25]
Huang, C., Malik, T., Cartaxo, J., et al. 2024, CompactObject: An open-source Python package for full-scope neutron star equation of state inference. https://arxiv.org/abs/2411.14615 14
-
[26]
2021, The Astrophysical Journal, 907, 63, doi: 10.3847/1538-4357/abcec0
Kazanas, D. 2021, The Astrophysical Journal, 907, 63, doi: 10.3847/1538-4357/abcec0
-
[27]
2026, arXiv preprint arXiv:2602.23743
Kini, Y., Mauviard, L., Salmi, T., et al. 2026, arXiv preprint arXiv:2602.23743
-
[28]
Kini, Y., Mauviard, L., Salmi, T., et al. 2026,A NICER View of PSR J0030+0451: Updated Constraints from Six Years of NICER Observations, Zenodo, doi: 10.5281/zenodo.18741941
-
[29]
2000, Monthly Notices of the Royal Astronomical Society, 315, L45
Lahav, O., Bridle, S., Hobson, M., Lasenby, A., & Sodre Jr, L. 2000, Monthly Notices of the Royal Astronomical Society, 315, L45
2000
-
[30]
Lattimer, J. M. 2021, Annual Review of Nuclear and Particle Science, 71, 433, doi: 10.1146/annurev-nucl-102419-124827
-
[31]
Lattimer, J. M., & Prakash, M. 2001, ApJ, 550, 426, doi: 10.1086/319702 —. 2007, PhR, 442, 109, doi: 10.1016/j.physrep.2007.02.003
-
[32]
Lattimer, J. M., & Prakash, M. 2016, Physics Reports, 621, 127–164, doi: 10.1016/j.physrep.2015.12.005
-
[33]
Lo, K. H., Miller, M. C., Bhattacharyya, S., & Lamb, F. K. 2013, ApJ, 776, 19, doi: 10.1088/0004-637X/776/1/19
-
[34]
Lockhart, W., Gralla, S. E., ¨Ozel, F., & Psaltis, D. 2019, MNRAS, 490, 1774, doi: 10.1093/mnras/stz2524
-
[35]
Lommen, A. N., Kipphorn, R. A., Nice, D. J., et al. 2006, ApJ, 642, 1012, doi: 10.1086/501067
-
[36]
Miller, M. C., & Lamb, F. K. 2015, ApJ, 808, 31, doi: 10.1088/0004-637X/808/1/31 —. 2016, European Physical Journal A, 52, 63, doi: 10.1140/epja/i2016-16063-8
-
[37]
C., Lamb, F
Miller, M. C., Lamb, F. K., Dittmann, A., et al. 2019, The Astrophysical Journal Letters, 887, L24
2019
-
[38]
Miller, M. C., Lamb, F. K., Dittmann, A. J., et al. 2019, NICER PSR J0030+0451 Illinois-Maryland MCMC Samples, 1.0.0, Zenodo, doi: 10.5281/zenodo.3473466 —. 2021, The Astrophysical Journal Letters, 918, L28, doi: 10.3847/2041-8213/ac089b Monta˜ na, G., Tol´ os, L., Hanauske, M., & Rezzolla, L. 2019, Physical Review D, 99, doi: 10.1103/physrevd.99.103009
-
[39]
Morsink, S. M., Leahy, D. A., Cadeau, C., & Braga, J. 2007, ApJ, 663, 1244, doi: 10.1086/518648 N¨ attil¨ a, J., & Pihajoki, P. 2018, A&A, 615, A50, doi: 10.1051/0004-6361/201630261
-
[40]
Oertel, M., Hempel, M., Kl¨ ahn, T., & Typel, S. 2017, Reviews of Modern Physics, 89, 015007, doi: 10.1103/RevModPhys.89.015007 ¨Ozel, F., & Freire, P. 2016, ARA&A, 54, 401, doi: 10.1146/annurev-astro-081915-023322
-
[41]
Pechenick, K. R., Ftaclas, C., & Cohen, J. M. 1983, ApJ, 274, 846, doi: 10.1086/161498
-
[42]
2025, Monthly Notices of the Royal Astronomical Society, 539, 2608
Phillips, A., & Kochanek, C. 2025, Monthly Notices of the Royal Astronomical Society, 539, 2608
2025
-
[43]
Accurate analytic formula for light bending in Schwarzschild metric.Astron
Poutanen, J. 2020, A&A, 640, A24, doi: 10.1051/0004-6361/202037471
-
[44]
Poutanen, J., & Beloborodov, A. M. 2006, MNRAS, 373, 836, doi: 10.1111/j.1365-2966.2006.11088.x
-
[45]
Press, W. H. 1996, arXiv preprint astro-ph/9604126
work page internal anchor Pith review arXiv 1996
-
[46]
Prigozhin, G., Gendreau, K., Doty, J. P., et al. 2016, in Society of Photo-Optical Instrumentation Engineers (SPIE) Conference Series, Vol. 9905, Space Telescopes and Instrumentation 2016: Ultraviolet to Gamma Ray, ed. J.-W. A. den Herder, T. Takahashi, & M. Bautz, 99051I, doi: 10.1117/12.2231718
-
[47]
K., Hebeler, K., Hinderer, T., et al
Raaijmakers, G., Greif, S. K., Hebeler, K., Hinderer, T., et al. 2021, The Astrophysical Journal Letters, 918, L29, doi: 10.3847/2041-8213/ac089a
-
[48]
Raaijmakers, G., Riley, T. E., Watts, A. L., et al. 2019, The Astrophysical Journal Letters, 887, L22, doi: 10.3847/2041-8213/ab451a
-
[49]
Reardon, D. J., Bailes, M., Shannon, R. M., et al. 2024, The Astrophysical Journal Letters, 971, L18, doi: 10.3847/2041-8213/ad614a
-
[50]
E., Watts, A
Riley, T. E., Watts, A. L., Bogdanov, S., et al. 2019, The Astrophysical Journal Letters, 887, L21
2019
-
[51]
Riley, T. E., Watts, A. L., Bogdanov, S., et al. 2019,A NICER View of PSR J0030+0451: Nested Samples for Millisecond Pulsar Parameter Estimation, Zenodo, doi: 10.5281/zenodo.3386449
-
[52]
Riley, T. E., Watts, A. L., Ray, P. S., et al. 2021, The Astrophysical Journal Letters, 918, L27, doi: 10.3847/2041-8213/ac0a81
-
[53]
2024, ApJL, 971, L19, doi: 10.3847/2041-8213/ad5f02
Rutherford, N., Mendes, M., Svensson, I., et al. 2024, ApJL, 971, L19, doi: 10.3847/2041-8213/ad5f02
-
[54]
Salmi, T., Choudhury, D., Kini, Y., et al. 2024, arXiv e-prints, arXiv:2406.14466, doi: 10.48550/arXiv.2406.14466
-
[55]
1986, Density Estimation for Statistics and Data Analysis, Chapman & Hall/CRC Monographs on Statistics & Applied Probability (Taylor & Francis)
Silverman, B. 1986, Density Estimation for Statistics and Data Analysis, Chapman & Hall/CRC Monographs on Statistics & Applied Probability (Taylor & Francis). https://books.google.com/books?id=e-xsrjsL7WkC
1986
-
[56]
2023, Physical Review C, 107, doi: 10.1103/physrevc.107.025801
Somasundaram, R., Tews, I., & Margueron, J. 2023, Physical Review C, 107, doi: 10.1103/physrevc.107.025801
-
[57]
Stairs, I. H. 2003, Living Reviews in Relativity, 6, doi: 10.12942/lrr-2003-5
-
[58]
Vinciguerra, S., Salmi, T., Watts, A. L., et al. 2023, An updated mass-radius analysis of the 2017-2018 NICER data set of PSR J0030+0451, Zenodo, doi: 10.5281/zenodo.8239000 15
-
[59]
Vinciguerra, S., Salmi, T., Watts, A. L., et al. 2024, ApJ, 961, 62, doi: 10.3847/1538-4357/acfb83
-
[60]
Virtanen, P., Gommers, R., Oliphant, T. E., et al. 2020, Nature Methods, 17, 261, doi: 10.1038/s41592-019-0686-2
-
[61]
Waskom, M. L. 2021, Journal of Open Source Software, 6, 3021, doi: 10.21105/joss.03021
-
[62]
Watts, A. L. 2019, in American Institute of Physics Conference Series, Vol. 2127, Xiamen-CUSTIPEN Workshop on the Equation of State of Dense Neutron-Rich Matter in the Era of Gravitational Wave Astronomy (AIP), 020008, doi: 10.1063/1.5117798
-
[63]
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
-
[64]
2025, arXiv e-prints, arXiv:2510.07764, doi: 10.48550/arXiv.2510.07764
Zhou, T., & Huang, C. 2025, arXiv e-prints, arXiv:2510.07764, doi: 10.48550/arXiv.2510.07764
work page internal anchor Pith review doi:10.48550/arxiv.2510.07764 2025
-
[65]
Zhou, T., & Huang, C. 2025, arXiv e-prints. https://arxiv.org/abs/2504.08662
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.