2025ApJ...995...60M

Four neutron star radius measurements have already been obtained by modeling the X-ray pulses of rotation-powered millisecond pulsars observed by the Neutron Star Interior Composition ExploreR (NICER). We report here the radius measurement of PSR J0614-3329 employing the same method with NICER and XMM-Newton data using Bayesian inference. For all different models tested, including one with unrestricted inclination prior, we retrieve very similar nonantipodal hot region geometries and radii. For the preferred model, we infer an equatorial radius of ${R}_{{\rm{eq}}}=10.2{9}_{-0.86}^{+1.01}\,$km for a mass of $M=1.4{4}_{-0.07}^{+0.06}\,{M}_{\odot }$ (median values with equal-tailed 68% credible interval), the latter being essentially constrained from radio timing priors obtained by MeerKAT. A more complex model, fitting the data equally well, resulted in a consistent inferred radius. We find that, for all different models, the pulse emission originates from two hot regions, one at the pole and the other at the equator. The resulting radius constraint is consistent with previous X-ray and gravitational wave measurements of neutron stars in the same mass range. Equation of state inferences, including previous NICER and gravitational wave results, slightly soften the equation of state with PSR J0614-3329 included and shift the allowed mass-radius region toward lower radii by ∼300 m, which is compatible with previous analyses to within less than one standard deviation.