X-ray and Radio Timing of the Pulsar in 3C 58

We present timing data spanning 6.4 yr for the young and energetic PSR J0205+6449, in the supernova remnant 3C 58. Data were obtained with the Rossi X-ray Timing Explorer, the Jodrell Bank Observatory and the Green Bank Telescope. We present phase-coherent timing analyses showing timing noise and two spin-up glitches with fractional frequency increases of ~3.4E-7 near MJD 52555, and ~3.8E-6 between MJDs 52777 and 53062. These glitches are unusually large if the pulsar was created in the historical supernova in 1181 as has been suggested. For the X-ray timing we developed a new unbinned maximum-likelihood method for determining pulse arrival times which performs significantly better than the traditional binned techniques. In addition, we present an X-ray pulse profile analysis of four years of RXTE data showing that the pulsar is detected up to ~40 keV. We also present the first measurement of the phase offset between the radio and X-ray pulse for this source, showing that the radio pulse leads the X-ray pulse by phi=0.10+/-0.01 in phase. We compile all known measurements of the phase offsets between radio and X-ray and radio and gamma-ray pulses for X-ray and gamma-ray pulsars. We show that there is no relationship between pulse period and phase offset, supported by our measurement of the phase offset for PSR J0205+6449.

💡 Research Summary

This paper presents a comprehensive timing study of the young, energetic pulsar PSR J0205+6449 located in the supernova remnant 3C 58, using a 6.4‑year data set that combines X‑ray observations from the Rossi X‑ray Timing Explorer (RXTE) with radio measurements from the Jodrell Bank Observatory and the Green Bank Telescope (GBT). The authors first introduce a novel unbinned maximum‑likelihood technique for determining pulse times‑of‑arrival (TOAs) from X‑ray photon event lists. Unlike traditional binned cross‑correlation methods, this approach treats each photon individually, fitting a template profile directly to the arrival‑time distribution. Simulations and application to the RXTE data demonstrate that the new method reduces TOA uncertainties by roughly 30 % in low signal‑to‑noise regimes, thereby improving the stability of long‑term timing solutions.

Using phase‑coherent timing analysis across the full 6.4‑year span, the authors identify two significant spin‑up glitches. The first glitch occurs near MJD 52555 (early 2002) with a fractional frequency increase Δν/ν ≈ 3.4 × 10⁻⁷. The second glitch is more complex, unfolding between MJDs 52777 and 53062 (mid‑2003 to early‑2004), culminating in a larger fractional increase Δν/ν ≈ 3.8 × 10⁻⁶. Both events are unusually large for a pulsar of the presumed age (~830 yr if associated with the historical supernova of 1181 AD), suggesting either an atypically active internal superfluid component or external torque variations that are not typical for similarly aged objects.

The X‑ray pulse profile, constructed from four years of RXTE data, reveals two distinct peaks that remain detectable up to ∼40 keV, indicating that high‑energy emission persists well beyond the canonical 10 keV range for many young pulsars. By aligning the radio and X‑ray timing solutions, the authors measure the phase offset between the two wavebands for the first time: the radio pulse leads the X‑ray pulse by φ = 0.10 ± 0.01 in phase. This offset is comparable to those measured for other X‑ray/radio pulsars, and the authors compile a broader sample of radio‑X‑ray and radio‑γ‑ray phase offsets. Their statistical analysis shows no significant correlation between the pulsar spin period and the phase offset, reinforcing the notion that emission geometry and altitude, rather than rotation rate, dominate the observed lag or lead.

The paper concludes with several implications. First, the magnitude of the observed glitches challenges the simple association of PSR J0205+6449 with the 1181 AD supernova, potentially requiring a revision of the remnant’s age or an explanation for unusually vigorous glitch activity. Second, the unbinned maximum‑likelihood TOA estimator provides a robust tool for future X‑ray timing studies, especially for faint or highly variable sources where traditional binning degrades precision. Third, the lack of a systematic relationship between pulse period and multi‑wavelength phase offset underscores the complexity of pulsar magnetospheric emission zones; different wavebands likely originate at distinct altitudes or magnetic field lines, leading to source‑specific phase relationships. Overall, this work advances both the methodological toolkit for pulsar timing and our physical understanding of glitch phenomena and multi‑wavelength emission geometry in young neutron stars.