2026ApJ..1000L..48M

Neutron Star Interior Composition Explorer (NICER) data have been used to estimate the masses and radii of the rotation-powered millisecond pulsars PSR J0030+0451, PSR J0740+6620, PSR J0437-4715, PSR J1231-1411, and PSR J0614-3329, sometimes in joint analyses with XMM-Newton data. These measurements provide invaluable information about the properties of cold, catalyzed matter beyond nuclear saturation density. Here, we present the results of our modeling of NICER data on PSR J0437-4715 using several different models of hot thermal X-ray-emitting spots on the stellar surface. For this pulsar, previous Nuclear Spectroscopic Telescope Array observations established that there is also a modulated nonthermal component to the emission, but the previously published analysis of NICER data did not model this component. We find that the Bayesian evidence is significantly higher when the modulated nonthermal component is included and that omission of this component leads to poorer fits to the bolometric NICER data and thus risks bias in the resulting radius estimates. Our models, which we pursue to inferential convergence, therefore have modulated nonthermal emission, and our headline model has in addition three uniform-temperature thermally emitting circular spots. Using this model, the symmetric 68% credible range in the radius is 11.9-15.5 km, which at the independently measured mass of M = 1.418 ± 0.044 M_⊙ is consistent with previous reports of the radius of the ∼1.4 M_⊙ pulsar PSR J0030+0451. We discuss the implications of this measurement for the equation of state of dense matter.