Pith. Sign in

REVIEW

A NICER view of PSR J1614$-$2230: a massive and compact millisecond pulsar

High Energy Astrophysical Phenomena · Submitted 31 Aug 2026 · 144 references

Pith's verdict: Solid, well-documented radius for a high-mass MSP; a genuine new data point for the EoS community, though low statistics keep the uncertainties moderate.

T0 · AI review

How this review was made: An AI referee (deepseek/deepseek-v4-flash) read the paper on September 9, 2026 and answered a fixed list of questions about it. T0 is the starting level every paper gets. The levels describe how deeply Pith has checked the paper, never the authors.

  1. T0 AI review An AI read the paper and answered a fixed list of questions.
  2. T1 Audit Automatic checks that the references and files the paper relies on exist and match.
  3. T2 Compute The paper's code or calculations re-run with the same results.
  4. T3 Formal statement The claims written out in a formal math language a computer can read.
  5. T4 Proof Those claims proved, and the proof checked by computer.

How Pith checks papers →

· 3 major objections · 5 minor · Ask about this paper · The compiled record · Read the full review →
Challenge this review Comment (0) Read on arXiv ↗
Cite ▾
arXiv ID2609.00172 (version 1)
Pith IDOJBFJRAN, Pith's permanent record of this paper and its review
CitationAnna L. Watts, Bas Dorsman, Christine Kazantsev, Daniela Huppenkothen, Denis Gonz\'alez-Caniulef, Devarshi Choudhury, Gilles Theureau, Isma\"el Cognard, Jean-Francois Olive, J\'er\^ome Novak, Lami Suleiman, Lucas Guillemot, Lucien Mauviard, Mariska Hoogkamer, Matthew Kerr, Natalie A. Webb, Pierre Stammler, Sebastien Guillot, Tuomo Salmi, Yves Kini. "A NICER view of PSR J1614$-$2230: a massive and compact millisecond pulsar." arXiv:2609.00172 (2026). https://pith.science/paper/OJBFJRAN
Cite this reviewBibTeX for Pith's review of this paper
Claim paper

Why Pith says this

Highlight any text to comment on it.

Referee report
The paper

The paper claims

This paper measures the radius of the massive millisecond pulsar PSR J1614−2230 to be 10.06 km, providing a second high-mass neutron-star radius constraint.

The reading

This paper uses pulse-profile modeling of X-ray data from NICER, XMM-Newton, and Chandra to measure the radius and mass of the millisecond pulsar PSR J1614−2230. The result is an equatorial radius of 10.$06^{{+1.25}}$_{-0.87} km at a mass of 1.$937^{{+0.012}}$_{-0.013} M_sun, making it the second most massive pulsar with a radius measurement. The authors find that the data favor a compact star across all tested hot-spot geometries and atmosphere models, with an upper radius bound of about 12 km in the 90th percentile. If correct, this measurement strengthens constraints on the equation of state at high density, supporting models that produce compact stars near 2 solar masses.

What carries the argument

Pulse-profile modeling using the X-PSI code, which simulates ray tracing from an oblate rotating neutron star surface with relativistic light bending and Doppler effects. The star's oblateness and surface gravity are computed from updated quasi-universal relations (QURs) tailored to the spin frequency and mass range of J1614, derived from a large set of equations of state using the LORENE/nrotstar numerical relativity code. The model is fit to NICER energy-resolved pulse profiles and to XMM-Newton and Chandra phase-averaged spectra via nested sampling with MultiNest, marginalizing over a phase-independent background.

Extended reading notes

Core claim

The paper establishes that PSR J1614−2230 has an equatorial radius of 10.$06^{{+1.25}}$_{-0.87} km and a gravitational mass of 1.$937^{{+0.012}}$_{-0.013} M_sun (68% credible intervals), based on joint pulse-profile modeling of NICER timing data and XMM-Newton/Chandra spectra. Regardless of the geometric complexity of the hot spots or the choice of atmosphere composition (fully ionized hydrogen, partially ionized hydrogen, or fully ionized helium), the radius posterior consistently favors values below 12.08 km at the 90th percentile. The best-fitting geometry features one hot spot near the stellar pole and another near the equator, a configuration that accounts for the single broad pulse and large unpulsed component seen at the high inclination of the system. The mass constraint is driven by the radio-timing prior from Shapiro delay, while the radius is informed by the combination of the pulse shape and the phase-averaged spectral energy distribution.

Load-bearing premise

The paper assumes that the pulsar's spin axis is perfectly aligned with the binary orbital axis, so that the radio-measured orbital inclination (sin i = 0.9999019) can be taken as the spin inclination, and that two circular hot spots with uniform temperature describe the entire X-ray emission from the surface.

Editorial extensions

If this is right

  • If this radius is correct, it rules out equations of state that produce radii larger than about 12 km at 1.9 M_sun, favoring stiff high-density matter.
  • Together with the radius measurement of PSR J0740+6620, it provides a second anchor at high mass for Bayesian inference of the neutron-star equation of state.
  • The absence of a detectable power-law component in the X-ray emission suggests that the observed flux is purely thermal, simplifying future modeling of this and similar sources.
  • The updated QURs for oblateness and surface gravity, which are now specific to the spin frequency of J1614, can be applied to other rapid rotators to improve radius accuracy.
  • The geometry with one polar and one equatorial hot spot, if generic among edge-on millisecond pulsars, may help constrain magnetic field configurations.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • The strong reliance on the radio-timing mass and inclination priors means that the radius constraint is only as secure as the assumption that spin and orbital axes are aligned. A future detection of misalignment would require a reanalysis and could widen the radius uncertainty.
  • The fact that the CST (ring) geometry, though unphysical, is favored by Bayesian evidence over simpler models when data are sparse underscores the need for physics-informed priors on hot-spot shaping, especially for faint sources.
  • The authors' inference that radius leans toward low values across all models suggests that deeper X-ray exposures are unlikely to shift the central value above 12 km but could substantially reduce the error bars.
Original abstract

Using pulse profile modeling, we obtain the mass-radius measurement of a millisecond pulsar (MSP) with data from the Neutron Star Interior Composition ExploreR, XMM-Newton and the Chandra X-ray Observatory. We report here the radius of PSR J1614$-$2230, the second most massive MSP confirmed by radio timing. All of the data sets are well described by a simple model composed of two circular hot spots. The final result yields an equatorial radius of $R_{\rm eq}=10.06 ^{+1.25}_{-0.87}\,$km, and a gravitational mass of $M=1.937^{+0.012}_{-0.013}\,M_{\odot}$ (equally tailed 68% credible intervals). Although a non-thermal component was previously reported at higher energies, we find no sign of it in either our phase-averaged or phase-resolved spectral analyses. Using new relations linking the compactness to oblateness or surface gravity, and tailored to the spin frequency of PSR J1614$-$2230, we infer a configuration with one hot spot near the pole, and another near the equator. The tight mass posterior is essentially informed by radio timing, while the radius constraint is not as tight due to the low source signal (8.5$\sigma$ X-ray pulse significance). However, over all geometries and atmosphere models tested, the radius posterior tends toward low values ($\lesssim12.08\,$km, 90th percentile in all cases).

The compiled record

The paper's own sentences, arranged, kept beside the review and separate from it. Each label is drawn from a closed vocabulary; the number beside it is the probability the classifier gave that label. Numbers in the tables are the authors', checked by code. No sentence in this record was written by a model.

Compiled by classifier-1.13.0 via TypeSafe API, question bank abstract-p1-20260923, 2026-09-23: 1 request, 6,586 input tokens.

What the abstract itself states, as probabilities: claims first 0.04 · states limitation 0.96 · human participants 0.01 · controlled comparison 0.04 · randomized assignment 0.01 · states quantitative result 0.99 · compares with prior results 0.81 · mentions code or data availability 0.02.

The abstract, sentence by sentence

arXiv listing abstract v1. Results and the main claim are underlined; .

Using pulse profile modeling, we obtain the mass-radius measurement of a millisecond pulsar (MSP) with data from the Neutron Star Interior Composition ExploreR, XMM-Newton and the Chandra X-ray Observatory. We report here the radius of PSR J1614$-$2230, the second most massive MSP confirmed by radio timing. All of the data sets are well described by a simple model composed of two circular hot spots. The final result yields an equatorial radius of $R_{\rm eq}=10.06 ^{+1.25}_{-0.87}\,$km, and a gravitational mass of $M=1.937^{+0.012}_{-0.013}\,M_{\odot}$ (equally tailed 68% credible intervals). Although a non-thermal component was previously reported at higher energies, we find no sign of it in either our phase-averaged or phase-resolved spectral analyses. Using new relations linking the compactness to oblateness or surface gravity, and tailored to the spin frequency of PSR J1614$-$2230, we infer a configuration with one hot spot near the pole, and another near the equator. The tight mass posterior is essentially informed by radio timing, while the radius constraint is not as tight due to the low source signal (8.5$\sigma$ X-ray pulse significance). However, over all geometries and atmosphere models tested, the radius posterior tends toward low values ($\lesssim12.08\,$km, 90th percentile in all cases).

background contribution method result interpretation main claim

What the classifier reports for each sentence. The number is the probability it gave the label.
#sentencerolethe authors' own workstates the main conclusionstrength
s0 Using pulse profile modeling, we obtain the mass-radius measurement of a millisecond pulsar (MSP) with data from the Neutron Star Interior Composition ExploreR, XMM-Newton and the Chandra X-ray Observatory. method 0.81 0.96 0.06 asserted 1.00
s1 We report here the radius of PSR J1614$-$2230, the second most massive MSP confirmed by radio timing. contribution 0.78 0.72 0.18 asserted 0.99
s2 All of the data sets are well described by a simple model composed of two circular hot spots. method 0.56 0.85 0.04 asserted 1.00
s3 The final result yields an equatorial radius of $R_{\rm eq}=10.06 ^{+1.25}_{-0.87}\,$km, and a gravitational mass of $M=1.937^{+0.012}_{-0.013}\,M_{\odot}$ (equally tailed 68% credible intervals). result 1.00 0.97 0.79 asserted 1.00
s4 Although a non-thermal component was previously reported at higher energies, we find no sign of it in either our phase-averaged or phase-resolved spectral analyses. result 1.00 0.95 0.04 asserted 1.00
s5 Using new relations linking the compactness to oblateness or surface gravity, and tailored to the spin frequency of PSR J1614$-$2230, we infer a configuration with one hot spot near the pole, and another near the equator. result 0.48 0.94 0.08 asserted 0.87
s6 The tight mass posterior is essentially informed by radio timing, while the radius constraint is not as tight due to the low source signal (8.5$\sigma$ X-ray pulse significance). limitation 0.88 0.84 0.05 asserted 1.00
s7 However, over all geometries and atmosphere models tested, the radius posterior tends toward low values ($\lesssim12.08\,$km, 90th percentile in all cases). result 0.90 0.93 0.12 asserted 0.67

What kind of paper the classifier reports this is: empirical study 1.00 · evidence from analysis of existing data 0.75.

What was examined, and what was not

sources read
  • arXiv listing abstract v1
not examined
  • full text
  • figures and tables
  • references and citing sentences
  • relations to other papers
levels run
  • P abstract
not run
  • S sections
  • T body sentences
  • C citations
  • R cross-paper relations

Every quoted sentence on this record is the authors'. Every label is from Pith's closed vocabulary, with the probability shown beside it. Every number in the measurements table is the authors', checked by code. Record as JSON

Figures
Figure 1. Data for J1614 obtained after reduction and phase-folding of the NICER events. The top panel displays the 0.3 − 3.0 keV pulse profile, displayed with the associated Poisson error in gray. The bottom panel shows the ener￾gy-resolved pulse profile of J1614. Phases are divided into 32 equal-size bins per cycle and the energies are split by PI channels. The plotted data are the exact same data that are used for inference. The data are duplicated over two rota￾tional phases for visualization purposes. observations of the target (see Section 3.4). This pro￾cedure resulted in a final…
Figure 2. Count spectra of the EPIC and Chandra instruments used in the present work. From left to right, PN, MOS1, MOS2 and ACIS. The source spectra are in blue, and the rescaled background spectra in green. The green shaded region represents the 3σ uncertainties on the background spectra. For visualization purposes only, the binned source and background spectra are provided with the associated 1σ error bars in violet and orange, respectively. For ACIS, the background is barely visible. and 28′′, centered on the source as well. This choice of background region avoids all the nearby sources…
Figure 3. Scatter plot of the dimensionless radius ratio R¯/Ω¯ 2 with respect to the compactness, computed using the adapted LORENE/nrotstar tool for 1000 EoS with 17 lati￾tude points. For each EoS, we compute R¯ for 10 central enthalpy values, chosen so that the masses are all within ±5σ of the mass range of J1614. The QUR from M. Al￾Gendy & S. M. Morsink (2014) is shown in black, while the reevaluated functional for J1614 is in purple. The statistical distribution (68% and 95%) of the points is represented in blue. while the previously used formula was valid across a large mass and frequency
Figure 4. Posterior distributions of the radius, compact￾ness, and mass for the ST-U model, using the previous and updated fit to the QUR, different data sets, as well as the headline result. Dashed–dotted black lines on the diagonal plots show the marginalized prior distributions, shared across all runs beside the previous QUR fit one. However, if present, the prior of the previous QUR run would be hardly distin￾guishable. The shaded vertical bands show the 68.3% credi￾ble intervals (CI) of the headline result. The contours in the 2D posteriors show the 68.3%, 95.4%, and 99.7% credible…
Figure 5. The CST maximum likelihood geometry as seen from Earth, without gravitational effects for clarity. The CST forms a ring going from pole to pole, whose width is thinner than the lines on this plot. We deem this non-head￾line model unrealistic, although it underlines the need for better physics-informed priors on the hot spot geometry. For most of these, we retrieve improved evidence over the ST-U model (see
Figure 6. Representation of the best-fit geometries, as seen from Earth, for the two modes present in the headline ST-U run. The pulsar is represented without gravitational or relativistic effects at phase zero of the data. The hotter spot is represented in blue, while the cooler one is in orange. These two modes are North/South mirror images due to degeneracy in the model. The secondary hot spots are located behind the star in both plots. Both geometries have similar likelihood. These representations do not encompass the full posterior distributions of the geometric parameters. 5.2. Geometry…
Figure 7. Posterior predictive distribution and maximum likelihood of the NICER observables for the headline ST-U geometry: bolometric pulse on the left panel and phase-averaged spectrum on the right panel. The posterior predictive distribution represents the posterior probability distribution of the observables. The contributions of the individual model components are shown in blue for the primary ST (polar hot spot), orange for the secondary ST (equatorial hot spot) and purple for the inferred background (BKG). The combined model components are in green for ST-U, and in light blue for the…
Figure 8. Mass-radius constraints obtained from MSPs with X-PSI. The J1614 posterior is the headline ST-U model from this work, while other posteriors were extracted from T. Salmi et al. (2024a) for PSR J0740+6620, D. Choud￾hury et al. (2024) for PSR J0437−4715, L. Mauviard et al. (2025) for PSR J0614−3329, and Y. Kini et al. (2026) for PSR J0030+0451. The contours show the 68.3% and 95.4% credible regions for each source. The causality limit is plotted here for visualization purposes only, as the limit used during inference is slightly different and is applied to the polar ra￾dius (the…
Figure 9. Model comparison plots using the Leave-One-Out (LOO) metric. Models are ordered by preference, with the better model on top. The gray horizontal lines provide the uncertainty of the difference of the LOO metric to the preferred topmost model. The dashed vertical line shows the LOO value for the preferred model. Left panel shows this comparison using the XMM and Chandra data curated in Section 2, while the right panel shows the results using the data from B. Pancrazi et al. (2012). Moreover, the normalization of the power law component in the NS+PL model with our dataset was on the…
Figure 10. Mock pl data alongside best fit model and residuals. Left panel: The spectrum was simulated and fitted with the non-truncated matrix. Right panel: The spectrum was simulated with the non-truncated matrix, but fitted with the truncated matrix with the default specextract option, showing a count deficit at low energies
Figure 11. Posterior distribution of the radius, hot spot geometric parameters, temperatures, and hydrogen column density for the headline ST-U run. Individual mode distributions (purple and orange) are plotted along with the full distribution (green). Mode 1 corresponds to the left plot from

Figures reproduced from arXiv:2609.00172 by the authors.

References

Discussion

Disagree with the review, or know something it missed? Leave a comment, highlight any passage above to comment on it, or challenge this review formally.