PSR J2021+4026: Gamma-ray Pulsar Modes

• PSR J2021+4026 is a radio-quiet gamma-ray pulsar in the Gamma Cygni supernova remnant characterized by quasi-periodic, abrupt reconfigurations of its magnetosphere.
• Multiwavelength observations across gamma-ray, X-ray, and optical bands reveal distinct shifts in flux, pulse profiles, and polar cap alignment during mode changes.
• Long-term Fermi-LAT monitoring indicates secular trends in gamma-ray flux and spin-down rate, offering insights into neutron star evolution and magnetospheric dynamics.

PSR J2021+4026 is a middle-aged (characteristic age ≈ 77 kyr), γ-ray–bright, radio-quiet pulsar in the Gamma Cygni supernova remnant (G78.2+2.1), located at an estimated distance of ~1.5 kpc. It is distinguished among the Fermi-LAT pulsar population as the only isolated γ-ray pulsar exhibiting quasi-periodic, abrupt, and reversible “mode changes,” in which its γ-ray flux and spin-down rate undergo coupled, step-like transitions on timescales of years. Multiwavelength observations (γ-ray, X-ray, and deep optical) reveal that these events correspond to global reconfiguration of the neutron star magnetosphere, likely involving higher-order magnetic multipoles near the stellar surface. PSR J2021+4026 provides a unique laboratory for studying the interplay between magnetospheric structure, particle acceleration, and high-energy pulsar emission.

1. Astrophysical Overview and Multiwavelength Identification

PSR J2021+4026 was discovered as a γ-ray pulsar in blind Fermi-LAT searches (f ≈ 3.7689 Hz, ṡE ≈ 10¹⁰³⁵ erg s⁻¹) and is coincident with the Gamma Cygni SNR shell (collaboration, 2013). Its X-ray counterpart was identified at RA = 20:21:30.733, Dec = +40:26:46.04 (J2000) through Chandra imaging, with a thermal spectrum indicative of hot polar cap emission plus a faint non-thermal tail (Weisskopf et al., 2011). Multi-epoch, deep XMM-Newton observations confirm pulsed, thermal, highly modulated X-ray emission coincident with the neutron star’s position and spin frequency [(Lin et al., 2013); (Rigoselli et al., 2020)].

Despite intensive searches, no radio pulsations have been detected, with stringent upper limits on periodic and single-pulse brightness at S₁.₅ GHz ≲ 0.2 mJy, classifying it as “radio-quiet” [(Shaw et al., 2023); (Trepl et al., 2010)]. Optical observations reach g′ > 26.1 and r′ > 25.3, ruling out any plausible non-degenerate companion or compact binary (Razzano et al., 2023). The X-ray–to–γ-ray and X-ray–to–optical flux ratios confirm the identification as an isolated neutron star.

2. Mode Changes: Timing, Spectral, and Pulse-Profile Phenomenology

PSR J2021+4026 exhibits unique, abrupt mode changes recurring on ~6–7 yr cycles (Fiori et al., 2024, Wang et al., 2023, Takata et al., 2020). The most prominent events occurred at MJD 55850 (2011-10-16), MJD 57000 (2014-12-09), MJD 58150 (2018-02-01), and MJD 59010 (2020-06-10):

| Mode Change | ΔF_γ/F_γ | Δ|ṽ|/|ṽ| | Timescale | Recovery | |-------------|----------------|---------------|------------------|--------------------| | 2011 (A→B) | –18% | +5.8% | ≲1 week | ~3 yr | | 2014 (B→C) | +20% | –5.7% | ~few months | HGF rest. | | 2018 (C→D) | –14% | +2–4% | ≲1 week | ~3 yr | | 2020 (D→E) | +16% | –3.0% | ~few months | HGF rest. |

Flux transitions between “high γ-ray flux / low |ṽ|” (HGF/LSD) and “low γ-ray flux / high |ṽ|” (LGF/HSD) states are accompanied by substantial pulse-profile changes: suppression and reappearance of the bridge (“BR”) emission, narrowing and amplitude reduction of the first γ-ray peak (P1), and phase-resolved softening of the spectrum [(collaboration, 2013); (Zhao et al., 2017); (Fiori et al., 2024)]. Phase-averaged spectral fits yield photon index transitions (ΔΓ ~ +0.1) and a cutoff-energy reduction (ΔE₍c₎ ~ 200–300 MeV) in low-flux states [(Takata et al., 2020); (collaboration, 2013)]. The timescales for transitions are ≪ orbital periods (instantaneous on Fermi monitoring timescales), with full recovery to the original state over ~months.

3. Interpretive Frameworks: Magnetospheric Reconfiguration Mechanisms

Detailed multiwavelength timing shows that mode changes are not associated with conventional spin glitches or crust-core vortex unpinning typical of radio-loud glitching pulsars (Zhao et al., 2017, Cozzolongo et al., 7 Jan 2026). Instead, the events represent abrupt, global rearrangements of the magnetospheric structure:

• Multipolar Magnetosphere Model: Surface magnetic topology is modeled as a superposition of a dipolar field (dominant at large radii) and a quadrupolar field (controlling near-surface geometry) (Razzano et al., 2023, Fiori et al., 2024). In quiescent states, particular quadrupole poles are mapped via open dipole lines to the outer-magnetosphere current sheet, leading to focused polar cap heating and a single broad X-ray pulse (Rigoselli et al., 2020). Mode changes can correspond to shifts in the quadrupole–current sheet connectivity, leading to reorganization of the current system, and predictable phase relocation of the thermal X-ray hot spot by Δφ ≈ 0.21 (Razzano et al., 2023).
• Force-free and Dissipative Magnetospheres: The global current and torque are governed by the inclination angle, α, and global magnetospheric conductivity, σ. The “outer-gap” accelerator model links the γ-ray power (L_γ ≈ f₍gap₎³L₍sd₎) and the return current to the size and pair-creation activity in the gap region (Ng et al., 2016). Small changes Δα or σ—whether by crustal plate movement (plate tectonics) or Hall-driven multipolar evolution—directly modulate both the emitted γ-ray flux and spin-down rate.
• Precession Scenario: Long-term, anti-correlated variations in |ṽ| (spin-down) and F_γ (γ-ray flux) can be modeled as damped free precession, with the inclination angle α modulated on ~6 yr timescales and the amplitude decaying due to internal friction, matching the observed cyclicity and amplitude evolution (Tong et al., 27 Jan 2025).
• Crustal Activity as Trigger: The mode transitions are plausibly triggered by small-scale crustal events—e.g., starquakes or plate motion—which perturb the local magnetic field at the polar cap, reconfiguring accelerator geometry and modifying global current flow (Ng et al., 2016, Takata et al., 2020). The fractional displacement and energetics are consistent with observed γ-ray and torque changes.

4. Multiwavelength Pulse Alignment and Polar-Cap Geometry

A key observation is the abrupt shift in the alignment between the X-ray pulse and the principal γ-ray pulse by Δφ = 0.21 ± 0.02 in phase, coincident with mode changes (Δφ moves from ~0.15 to ~0.36 after the first state transition) (Razzano et al., 2023). Maximum-likelihood fits to the unbinned phaseograms from XMM-Newton and Fermi-LAT confirm this shift, whereas traditional cross-correlation methods yield substantially greater uncertainties.

The broad single X-ray pulse is best explained by heated polar-cap emission from only one magnetic pole, requiring a highly asymmetric magnetospheric geometry (Wang et al., 2018, Rigoselli et al., 2020). The polar cap inferred from blackbody modeling has a radius of ≈ 340 m, while magnetized hydrogen atmosphere modeling yields an extended (~5–6 km) spot with temperature T ≈ 1 MK. This spot size is much larger than canonical dipole polar caps, suggesting heating over an extended zone—likely regulated by multipolar field components. The return-current luminosity and pulsed fraction are consistent with outer-gap accelerator predictions.

5. Long-Term Evolution and Secular Trends

Seventeen years of Fermi-LAT monitoring reveal that, superimposed on the discrete state transitions, PSR J2021+4026 exhibits additional secular evolution in its γ-ray flux and torque (Liu et al., 20 Jan 2026):

• Secular Evolution Phases: The jump-corrected energy flux δF_γ(t) undergoes a three-phase piecewise-linear evolution: (1) ~10 yr of rise (+2.0% yr⁻¹); (2) ~6 yr of decline (–3.7% yr⁻¹); (3) a rapid recent rise (+15% yr⁻¹).
• Long-Term Convergence: The mean LGF mode flux baseline is gradually approaching the HGF level at +0.7% yr⁻¹, suggesting a dissipative relaxation toward a long-term equilibrium.
• Enhanced Torque Noise: Substantial stochastic variability in |ṽ| is observed within states, with enhanced amplitude during flux-decline episodes. The secular variations in δ|ṽ| within each state do not correlate linearly with δF_γ, contrasting with their anti-correlation across state transitions.

These findings imply that both fast magnetospheric switches and slow dissipative evolution contribute to the long-term emission behavior. The secular trend is qualitatively consistent with damping of precessional wobble, as modeled in (Tong et al., 27 Jan 2025).

6. Open Questions, Uniqueness, and Comparative Context

PSR J2021+4026 is unique among >300 Fermi-LAT pulsars in showing tight, repeatable coupling of γ-ray emission to spin-down changes via abrupt, global magnetospheric reconfigurations (Cozzolongo et al., 7 Jan 2026, Fiori et al., 2024). Analogous mode-switching in radio pulsars is observed only at much lower energies, with different phenomenology. No direct evidence connects classic large-amplitude neutron-star glitches to γ-ray variability in other young pulsars or to long-term emission shifts (e.g., Vela shows no such changes across glitches) (Wang et al., 2023).

No radio emission—neither periodic nor RRAT-like—is detected in any state despite deep searches (Shaw et al., 2023), possibly due to unfavorable beaming or intrinsic faintness. The neutron star remains consistently radio-quiet, strengthening its classification as a Geminga-like object (Trepl et al., 2010).

A plausible implication is that only PSR J2021+4026 occupies a parameter space (field configuration, age, surface composition, internal structure) that allows for such frequent, non-destructive, multipolar field reconfigurations or that hosts a combination of secular precession and magnetospheric switching. The detailed mechanism—whether controlled by Hall-driven crustal multipole evolution, quasiperiodic crust-tectonics, or slow precession—remains under active investigation.

7. Future Prospects and Observational Tests

Continued γ-ray monitoring with Fermi-LAT (and successor instruments) is crucial to resolve the long-term convergence of emission states, to identify potential damping or acceleration of the secular modulation, and to constrain the timescale and amplitude of future mode changes (Liu et al., 20 Jan 2026). Next-generation X-ray observatories (e.g., Einstein Probe, eXTP) can test the predicted phase shifts, spot geometry, and pulse-profile evolution with higher precision, especially across future transitions (Razzano et al., 2023, Wang et al., 2018).

Modelling efforts using time-dependent global dissipative MHD and PIC simulations, including realistic multipolar fields and crust-magnetosphere coupling, are essential to verify whether the observed phenomenology can be reproduced from first principles (Fiori et al., 2024).

The unique phenomenology of PSR J2021+4026 offers the most direct constraint to date on how neutron star magnetospheres couple internal and external physics, how higher-order multipoles shape emission, and how global state transitions modulate the observables across the electromagnetic spectrum. Its continuing study is therefore central to understanding mode-changing behavior and the broader dynamics of rotation-powered pulsars.