A NICER view of the millisecond pulsar PSR J2124$-$3358: evidence for a helium atmosphere

Abstract

Pulse profile modeling has proven to be a powerful technique for determining the mass and radius of neutron stars. To date, this method has been applied to a handful of millisecond pulsars observed by the Neutron Star Interior Composition Explorer (NICER). However, analyses of more millisecond pulsars are necessary to determine tight constraints on the equation of state of superdense matter. In this study, we present an analysis of the isolated, rotation-powered millisecond pulsar PSR J2124$-$3358 using the X-ray Pulse Simulation and Inference (X-PSI) package, a publicly available state-of-the-art code for neutron-star relativistic ray tracing and Bayesian parameter inference. We use NICER and Chandra observations of this pulsar, exploring different neutron star atmospheric compositions and different configurations of the hot polar caps responsible for the pulsed X-ray emission. Our analyses favor a helium atmospheric composition, plausibly originating from accretion and subsequent evaporation of a former hydrogen-depleted binary companion. For this composition, and given the faint nature of the source and the low signal-to-noise of the data sets, we obtain broad posterior distributions yielding a mass $M = 1.8\pm0.5\,M_\odot$ and an equatorial radius $R_{\mathrm{eq}} = 11.7^{+2.6}_{-3.0}$ km (medians and $68\%$ credible intervals), and infer a configuration consisting of two slightly non-antipodal hot spots. By contrast, when using a hydrogen atmosphere model, the mass and radius decrease by $\sim 0.5\,M_\odot$ and $\sim 1$ km, respectively. Future multiwavelength studies, particularly those incorporating radio and gamma-ray pulse-emission, may provide tighter constraints on the geometry and physical properties of this source.