Radius of the neutron star PSR J0614-3329 from NICER data ✛
Context
Neutron stars compress matter to densities above nuclear density, which cannot be reproduced in laboratories. Therefore, astronomical measurements of the masses and radii of these objects serve as the primary source of knowledge about the equation of state of ultradense matter. In recent years, the NICER telescope has measured the radii of several pulsars, including PSR J0030+0451, PSR J0740+6620, and PSR J0437-4715. This work adds another point to the mass-radius diagram, refining the behavior of matter at densities several times nuclear saturation density.
Methods
The researchers analyzed 1.09 megaseconds of NICER exposure of PSR J0614-3329. They binned X-ray photons by rotation phase and energy, producing a two-dimensional picture of the pulsations. They then fitted a surface model with hot spots, considering configurations with two circular, two oval, and three circular spots. For each photon, the path in curved spacetime was traced using the oblate Schwarzschild approximation. Bayesian analysis with the pocoMC sampler yielded posterior distributions for mass and radius. Additionally, joint fits with XMM-Newton data were performed, but they revealed a systematic deficit of the model relative to observations, so the primary result is based solely on NICER data.
Results
The model with three circular spots was preferred, with the logarithm of the Bayes factor exceeding 10 compared to alternatives. For the equatorial radius, a symmetric 68% confidence interval of 9.88–12.77 km was obtained. The stellar mass is barely constrained by the data and remains near the prior value of 1.44 solar masses. Goodness-of-fit testing gave χ²/dof = 5154.89/5252 (probability 83%), indicating that the model is adequate. Joint analysis with XMM-Newton under the standard assumption of a known background leads to a systematic deficit of predicted counts: for XMM-pn, the deficit reaches 3σ. This points to a possible additional background component not accounted for in the standard processing.
Implications
The new radius measurement of PSR J0614-3329 refines the equation of state of neutron matter. Including this result in the global analysis lowers the median estimate of the maximum mass of a nonrotating neutron star from 2.20 to 2.15 solar masses and reduces the radius of a canonical 1.4-solar-mass star from 12.56 to 12.44 km. This is a small but systematic shift toward a softer equation of state. The work also highlights a problem in combining NICER and XMM-Newton: standard XMM-Newton background processing may lead to an underestimated radius, which is important for future joint analyses.
Future development
Further NICER observations of other pulsars will allow accumulating statistics and narrowing the allowed region of the equation of state. Of particular interest are stars with masses around 2 solar masses, where sensitivity to the stiffness of the equation of state is maximal. Developing methods to account for the XMM-Newton background, for example, by introducing an additional free component, could resolve the discrepancy between instruments and improve the accuracy of joint fits. Also promising is the use of data from future X-ray missions with improved time resolution.
Impact
The results will affect nuclear physics by constraining theoretical models of dense matter, and astrophysics of compact objects by refining the mass-radius relation for neutron stars.
Next steps
Conduct a similar analysis for other pulsars with known masses from radio observations. Develop a more realistic model of the XMM-Newton background for joint fits.
Key open problems
The main unresolved problem in this field is determining the equation of state of matter at densities above nuclear density. This work refines constraints on that equation but does not resolve the question of phase transitions, such as the appearance of hyperons or quark matter. Also open is the issue of systematic errors when combining different X-ray instruments.
The knowledge machine's view
This section is addressed to the author of the paper and their colleagues. Neighbouring works were found by meaning-based search across our archive, not by word overlap. We do not review: we show the field around the work.
What the work does
This work presents a Bayesian analysis of 1.09 Ms of NICER data for the millisecond pulsar PSR J0614-3329, using hot spot modeling (two circular, two oval, three circular) and accounting for spacetime curvature in the oblate Schwarzschild approximation. The inferred radius is 9.88–12.77 km (68% CI), consistent with previous measurements but shifted toward lower values; the three-circular-spot model is preferred (log Bayes factor >10). A joint analysis with XMM-Newton reveals a systematic underestimation by the model, likely due to an incomplete background treatment.
Where its strength is
A strong point is the careful comparison of hot spot configurations and the identification of systematics in the combined NICER and XMM-Newton data, which is important for future joint analyses. The resulting radius constraint adds a valuable point to the mass–radius diagram and refines the equation of state.
Where it could go next
- It might be interesting to compare the obtained radius constraint with results from other NICER measurements, for example for PSR J0030+0451 or PSR J0740+6620, to check consistency and refine the equation of state.based on: Pulsar Mystery: How Hot Spots Revealed the Star's Size arXiv:2602.23743
- A related approach uses generative models to reconstruct the neutron star equation of state; it could be applied to interpret your new radius constraint in terms of microscopic properties of matter.based on: Generative AI Reconstructs Neutron Star Matter arXiv:2608.17457v1
- There is work in the archive discussing constraints on the equation of state from causality and maximum mass; your result could be included in such an analysis to test compatibility with physical limits.based on: The Cosmic Tailor: How Causality Takes Measurements of Neutron Stars arXiv:2606.25186
Worth trying
- Test how the radius estimate changes when using alternative XMM-Newton background models, to reduce systematics.
- Apply Bayesian model comparison with different numbers of spots to other NICER pulsars to unify the analysis.
- Investigate the impact of the metric approximation (e.g., using an exact solution instead of the oblate Schwarzschild) on the radius estimate.
Why this may matter
Refining the neutron star radius directly impacts our understanding of the equation of state of superdense matter, a key challenge in modern astrophysics. The identified systematics in the joint NICER and XMM-Newton data are important for future multi-instrument analyses.
Nearby works in our archive
- Pulsar Mystery: How Hot Spots Revealed the Star's Sizesimilarity 0.71
- How to Weigh a Neutron Star Without Trickssimilarity 0.61
- Cosmic Heavyweight: Neutron Star Sets Mass Recordsimilarity 0.60
- Generative AI Reconstructs Neutron Star Mattersimilarity 0.59
- Stellar Interiors: A Clash of Equationssimilarity 0.58
- The Cosmic Tailor: How Causality Takes Measurements of Neutron Starssimilarity 0.58
The archive is crowded around this work: 1 review nearby, closest at 0.71. This ground is well walked — and there is plenty to compare against.
PSR J0614-3329 spins 318 times per second, meaning a point on its equator moves at about 10% of the speed of light.
Key numbers
- Radius (68% interval): 9.88–12.77 km
- Mass (prior): 1.44 ± 0.07 solar masses
- Rotation period: 3.14 ms (318 Hz)
- NICER exposure: 1.09 Ms
- Log Bayes factor: >10 in favor of three spots