Pulsed Gamma-rays from PSR J2021+3651 with the Fermi Large Area Telescope

We report the detection of pulsed gamma-rays from the young, spin-powered radio pulsar PSR J2021+3651 using data acquired with the Large Area Telescope (LAT) on the Fermi Gamma-ray Space Telescope (formerly GLAST). The light curve consists of two narrow peaks of similar amplitude separated by 0.468 +/- 0.002 in phase. The first peak lags the maximum of the 2 GHz radio pulse by 0.162 +/- 0.004 +/- 0.01 in phase. The integral gamma-ray photon flux above 100 MeV is (56 +/- 3 +/- 11) x 10^{-8} /cm2/s. The photon spectrum is well-described by an exponentially cut-off power law of the form dF/dE = kE^{-\Gamma} e^(-E/E_c) where the energy E is expressed in GeV. The photon index is \Gamma = 1.5 +/- 0.1 +/- 0.1 and the exponential cut-off is E_c = 2.4 +/- 0.3 +/- 0.5 GeV. The first uncertainty is statistical and the second is systematic. The integral photon flux of the bridge is approximately 10% of the pulsed emission, and the upper limit on off-pulse gamma-ray emission from a putative pulsar wind nebula is <10% of the pulsed emission at the 95% confidence level. Radio polarization measurements yield a rotation measure of RM = 524 +/- 4 rad/m^2 but a poorly constrained magnetic geometry. Re-analysis of Chandra data enhanced the significance of the weak X-ray pulsations, and the first peak is roughly phase-aligned with the first gamma-ray peak. We discuss the emission region and beaming geometry based on the shape and spectrum of the gamma-ray light curve combined with radio and X-ray measurements, and the implications for the pulsar distance. Gamma-ray emission from the polar cap region seems unlikely for this pulsar.

💡 Research Summary

The paper reports the first detection of pulsed γ‑ray emission from the young, spin‑down powered radio pulsar PSR J2021+3651 using data from the Large Area Telescope (LAT) aboard the Fermi Gamma‑ray Space Telescope. By folding LAT photons with a contemporaneous radio ephemeris obtained at 2 GHz, the authors construct a high‑precision γ‑ray light curve that exhibits two narrow peaks (P1 and P2). The peaks are separated by 0.468 ± 0.002 in rotational phase, and the first γ‑ray peak lags the radio main pulse by 0.162 ± 0.004 (statistical) ± 0.01 (systematic) in phase. The full‑pulse photon flux above 100 MeV is (56 ± 3 ± 11) × 10⁻⁸ cm⁻² s⁻¹, corresponding to an energy flux of (3.5 ± 0.2 ± 0.5) × 10⁻¹⁰ erg cm⁻² s⁻¹.

Spectral analysis shows that the phase‑averaged spectrum is well described by an exponentially cut‑off power law, dF/dE = k E⁻Γ exp(−E/E_c), with photon index Γ = 1.5 ± 0.1 (stat) ± 0.1 (sys) and cutoff energy E_c = 2.4 ± 0.3 (stat) ± 0.5 (sys) GeV. The bridge emission (the region between the two peaks) contributes roughly 10 % of the total pulsed flux, while the off‑pulse interval yields no significant detection; an upper limit of <10 % of the pulsed flux is placed at the 95 % confidence level, indicating that any nebular or unpulsed component is weak.

Radio polarization measurements give a rotation measure RM = 524 ± 4 rad m⁻², confirming a dense, magnetized line‑of‑sight plasma, but the geometry (magnetic inclination α and viewing angle ζ) remains poorly constrained. Re‑analysis of archival Chandra X‑ray data reveals weak X‑ray pulsations (~3 % pulsed fraction) whose first peak aligns in phase with the γ‑ray P1, suggesting a common emission region for high‑energy photons.

The authors discuss the implications for emission geometry. The observed phase lag of the γ‑ray peak relative to the radio pulse, together with the peak separation, is consistent with outer‑gap or slot‑gap models where γ‑rays are produced at high altitudes in the magnetosphere, well above the radio emission region. The relatively hard photon index and GeV‑scale cutoff are typical of curvature‑radiation‑limited spectra expected from such outer‑magnetospheric zones, whereas polar‑cap models generally predict lower cutoff energies and smaller phase lags, making them unlikely for PSR J2021+3651.

Distance estimates based on the dispersion measure (DM ≈ 369 pc cm⁻³) and the measured RM suggest a distance of roughly 1.8–2.5 kpc. At this distance, the γ‑ray efficiency (γ‑ray luminosity divided by spin‑down power) is on the order of 10 %, comparable to other young, energetic γ‑ray pulsars.

In summary, the detection of two narrow, well‑separated γ‑ray peaks with a GeV cutoff, the phase relationship to radio and X‑ray emission, and the lack of significant off‑pulse flux all point toward high‑altitude outer‑magnetospheric emission mechanisms for PSR J2021+3651. The results strengthen the case that most young, energetic pulsars emit γ‑rays from outer gaps or slot gaps rather than from polar caps, and they provide a valuable multi‑wavelength benchmark for testing pulsar emission models. Future observations with longer LAT exposure and higher‑precision radio polarization will refine the geometry (α, ζ) and may uncover faint nebular emission, further elucidating the particle acceleration processes in this system.