2023MNRAS.524.5854S

We report on a likelihood-stacking search for γ-ray pulsars at 362 high-latitude locations that coincide with known radio pulsar positions. We observe a stacked signal conservatively 2.5σ over the background. Stacking their likelihood profiles in spectral parameter space implies a pulsar-like spectral index and a characteristic flux a factor of 2 below the Fermi Large Area Telescope point-source sensitivity, assuming isotropic/unbeamed emission from all sample pulsars. The same procedures performed on empty control fields indicate that the pulsars as a population can be distinguished from the background with a Δ(TS) = 28, where TS refers to test statistic, at the peak location (or 4.8σ), and the stacked spectra of the control fields are distinctly softer than those of the pulsars. This study also probes a unique region of parameter space populated by low $\dot{E}$ pulsars, most of which have no γ-ray ephemeris available, and is sensitive to high duty cycles. We also discuss the possible γ-ray emission mechanism from such pulsars.