Christine Kazantsev

Accreting millisecond pulsars (AMSPs) and transitional millisecond pulsars (tMSPs) are neutron star low-mass X-ray binaries that can evolve into recycled radio millisecond pulsars. In both types of systems, X-ray pulsations have been detected during phases of X-ray activity when matter accretion through a disc was turned on. On the other hand, when accretion stops, and these systems enter the quiescent low-luminosity X-ray state, and only tMSPs become visible as radio pulsars. Despite several attempts, radio pulsations have never been detected in quiescent AMSPs, except for IGR J18245−2452. We present the results of two observational campaigns performed on the AMSPs IGR J00291+5934 and MAXI J1957+032 with the Five-hundred-meter Aperture Spherical Telescope ( FAST ) in L band (1–1.5 GHz). The two sources have most likely been observed in quiescence, as suggested by the upper limits on their X-ray and optical flux obtained with Swift and the Las Cumbres Observatory, respectively. We performed a deep search for coherent periodicities in radio, but found no significant candidate signal at the known spin frequency of these sources or at other frequencies. Assuming a pulse duty cycle of 10%, we derive upper limits on the pulsed radio flux density of 3.3 μJy and 5.6 μJy for IGR J00291+5934 and MAXI J1957+032, respectively, which are the most stringent limits for any known persistent AMSP so far.

Pulse-profile modeling of rotation-powered millisecond pulsars targeted by NICER has enabled mass–radius constraints of several neutron star sources, with implications for the dense-matter equation of state. For the bright isolated pulsar PSR J0030+0451, the inferred mass and radius were previously found to depend strongly on the assumed hot-spot model. These hot-spot models yielded different mass–radius constraints, with the statistically preferred model exhibiting some mild tension with results inferred for PSR J0437−4715, PSR J0614−3329, and GW170817. We present an updated pulse-profile analysis of PSR J0030+0451 using new NICER observations obtained between 2017 July to 2023 January, increasing the number of X-ray counts by about 50% compared to previous analyses. We jointly analyze the NICER data with archival XMM-Newton observations to better constrain the source spectrum and background. The new analysis, using more conservative nested-sampling settings together with the new NICER observations, significantly reduces the discrepancy between the hot-spot models. The inferred mass and radius are M = 1.4 3 − 0.17 + 0.20 M ⊙ and R eq = 12.6 8 − 1.04 + 1.31 km (68% credible intervals), reducing the tension with the results from other sources. In addition, the inferred hot-spot configurations suggest the presence of intraspot temperature gradients.

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 eq = 10.2 9 − 0.86 + 1.01 km for a mass of M = 1.4 4 − 0.07 + 0.06 M ⊙ (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.