2022ApJ...934...30K

The Fermi Large Area Telescope receives ≪1 photon per rotation from any γ-ray pulsar. However, out of the billions of monitored rotations of the bright pulsars Vela (PSR J0835-4510) and Geminga (PSR J0633+1746), a few thousand have ≥2 pulsed photons. These rare pairs encode information about the variability of pulse amplitude and shape. We have cataloged such pairs and find the observed number to be in good agreement with simple Poisson statistics, limiting any amplitude variations to <19% (Vela) and <22% (Geminga) at 2σ confidence. Using an array of basis functions to model pulse-shape variability, the observed pulse phase distribution of the pairs limits the scale of pulse-shape variations of Vela to <13%, while for Geminga we find a hint of ∼20% single-pulse-shape variability most associated with the pulse peaks. If variations last longer than a single rotation, more pairs can be collected, and we have calculated upper limits on amplitude and shape variations for assumed coherence times up to 100 rotations, finding limits of ∼1% (amplitude) and ∼3% (shape) for both pulsars. Because a large volume of the pulsar magnetosphere contributes to γ-ray pulse production, we conclude that the magnetospheres of these two energetic pulsars are stable over one rotation and very stable on longer timescales. All other γ-ray pulsars are too faint for similar analyses. These results provide useful constraints on rapidly improving simulations of pulsar magnetospheres, which have revealed a variety of large-scale instabilities in the thin equatorial current sheets where the bulk of GeV γ-ray emission is thought to originate.