2015ApJ...804...86M

Due to the low gamma-ray fluxes from pulsars above 50 GeV and the small collecting area of space-based telescopes, the gamma-ray emission discovered by the Fermi Large Area Telescope (LAT) in ∼150 pulsars is largely unexplored at these energies. In this regime, the uncertainties on the spectral data points and/or the constraints from upper limits are not sufficient to provide robust tests of competing emission models in individual pulsars. The discovery of power-law-type emission from the Crab pulsar at energies exceeding 100 GeV provides a compelling justification for exploration of other pulsars at these energies. We applied the method of aperture photometry to measure pulsar emission spectra from Fermi-LAT data and present a stacked analysis of 115 pulsars selected from the Second Fermi-LAT catalog of gamma-ray pulsars. This analysis, which uses an average of ∼4.2 yr of data per pulsar, aggregates low-level emission which cannot be resolved in individual objects but can be detected in an ensemble. We find no significant stacked excess at energies above 50 GeV. An upper limit of 30% of the Crab pulsar level is found for the average flux from 115 pulsars in the 100-177 GeV energy range at the 95% confidence level. Stacked searches exclusive to the young pulsar sample, the millisecond pulsar sample, and several other promising sub-samples also return no significant excesses above 50 GeV.