Precise Gamma-Ray Timing and Radio Observations of 17 Fermi Gamma-Ray Pulsars

We present precise phase-connected pulse timing solutions for 16 gamma-ray-selected pulsars recently discovered using the Large Area Telescope (LAT) on the Fermi Gamma-ray Space Telescope plus one very faint radio pulsar (PSR J1124-5916) that is more effectively timed with the LAT. We describe the analysis techniques including a maximum likelihood method for determining pulse times of arrival from unbinned photon data. A major result of this work is improved position determinations, which are crucial for multi-wavelength follow up. For most of the pulsars, we overlay the timing localizations on X-ray images from Swift and describe the status of X-ray counterpart associations. We report glitches measured in PSRs J0007+7303, J1124-5916, and J1813-1246. We analyze a new 20 ks Chandra ACIS observation of PSR J0633+0632 that reveals an arcminute-scale X-ray nebula extending to the south of the pulsar. We were also able to precisely localize the X-ray point source counterpart to the pulsar and find a spectrum that can be described by an absorbed blackbody or neutron star atmosphere with a hard powerlaw component. Another Chandra ACIS image of PSR J1732-3131 reveals a faint X-ray point source at a location consistent with the timing position of the pulsar. Finally, we present a compilation of new and archival searches for radio pulsations from each of the gamma-ray-selected pulsars as well as a new Parkes radio observation of PSR J1124-5916 to establish the gamma-ray to radio phase offset.

💡 Research Summary

The paper presents phase‑connected timing solutions for a set of seventeen pulsars: sixteen gamma‑ray‑selected pulsars discovered with the Fermi Large Area Telescope (LAT) and one very faint radio pulsar, PSR J1124‑5916, whose timing is actually more precise using LAT data. The authors introduce a maximum‑likelihood method for extracting pulse times‑of‑arrival (TOAs) directly from unbinned photon events, overcoming the limitations of traditional binned approaches that suffer when photon counts are low or background is high. By modeling the probability density of each photon’s rotational phase and maximizing the joint likelihood, they obtain TOAs with sub‑millisecond precision even for the faintest sources.

Applying this technique, the authors simultaneously fit rotation frequency, its first and second derivatives, and sky position for each pulsar. The resulting positional uncertainties shrink from arc‑minute scales typical of LAT catalog positions to a few tens of arcseconds, a dramatic improvement that enables direct overlay of timing localizations on X‑ray images from Swift and Chandra.

The high‑precision positions are used to investigate X‑ray counterparts. For PSR J0633+0632, a new 20 ks Chandra ACIS observation reveals an arcminute‑scale nebular structure extending southward, indicative of a pulsar wind nebula (PWN) or shock‑driven outflow. The point‑source counterpart is precisely localized, and its spectrum can be described by an absorbed blackbody (or neutron‑star atmosphere) plus a hard power‑law tail, suggesting thermal surface emission combined with magnetospheric non‑thermal radiation. A separate Chandra ACIS image of PSR J1732‑3131 shows a faint X‑ray point source coincident with the timing position, providing the first direct X‑ray detection of this object.

The timing analysis also uncovers rotational glitches in three pulsars: PSR J0007+7303, PSR J1124‑5916, and PSR J1813‑1246. For each glitch the authors report the fractional change in spin frequency (Δν/ν) and the accompanying change in spin‑down rate (Δ\dotν/\dotν), supplying valuable constraints on neutron‑star interior physics, such as superfluid vortex unpinning and crustal stress release.

In parallel, the paper compiles new and archival radio searches for each gamma‑ray‑selected pulsar, including a fresh Parkes observation of PSR J1124‑5916. Although most of the gamma‑ray‑selected pulsars remain radio‑quiet, the authors succeed in measuring the gamma‑ray‑to‑radio phase offset for PSR J1124‑5916, an essential parameter for modeling emission geometry.

Overall, the study demonstrates that (1) a likelihood‑based TOA extraction from unbinned LAT photons dramatically improves timing precision, (2) refined positions enable multi‑wavelength counterpart identification and detailed X‑ray morphological studies, (3) glitch measurements add to the growing catalog of neutron‑star rotational irregularities, and (4) coordinated radio observations, even when yielding non‑detections, help define the radio‑quiet nature of many LAT pulsars. The methodology and results set a new benchmark for pulsar timing with gamma‑ray data and illustrate the power of combining high‑energy timing with X‑ray imaging and radio searches to probe the physics of rotation‑powered neutron stars.