Evidence of TeV Halos Around Millisecond Pulsars
  Using data from the HAWC gamma-ray Telescope, we have studied a sample of 37
millisecond pulsars (MSPs), selected for their spindown power and proximity.
From among these MSP, we have identified four which favor the presence of very
high-energy gamma-ray emission at a level of $(2\Delta \ln \mathcal{L})^{1/2}
\ge 2.5$. Adopting a correlation between the spindown power and gamma-ray
luminosity of each pulsar, we performed a stacked likelihood analysis of these
37 MSPs, finding that the data supports the conclusion that these sources emit
very high-energy gamma-rays at a level of $(2\Delta \ln \mathcal{L})^{1/2} =
4.24$. Among sets of randomly selected sky locations within HAWC's
field-of-view, less than 1\% of such realizations yielded such high statistical
significance. Our analysis suggests that MSPs produce very high-energy
gamma-ray emission with a similar efficiency to that observed from the Geminga
TeV-halo, $\eta_{\rm MSP} = (0.39-1.08) \times \eta_{\rm Geminga}$. This
conclusion poses a significant challenge for pulsar interpretations of the
Galactic Center gamma-ray excess, as it suggests that any population of MSPs
potentially capable of producing the GeV excess would also produce TeV-scale
emission in excess of that observed by HESS from this region. Future
observations by CTA will be able to substantially clarify this situation.
