Return of the templates: Revisiting the Galactic Center excess with multimessenger observations

Abstract

The Galactic center excess (GCE) remains one of the most intriguing discoveries from the Fermi Large Area Telescope (LAT) observations. We revisit the characteristics of the GCE by first producing a new set of high-resolution galactic diffuse gamma-ray emission templates. This diffuse emission, which accounts for the bulk of the observed gamma rays, is ultimately due to cosmic-ray interactions with the interstellar medium. Using recent high-precision cosmic-ray observations, in addition to the continuing Fermi-LAT observations and observations from lower energy photons, we constrain the properties of the galactic diffuse emission. We describe a large set of diffuse gamma-ray emission templates which account for a very wide range of initial assumptions on the physical conditions in the inner galaxy. The broad properties of the GCE that we find in this work are qualitatively unchanged despite the introduction of this new set of templates, though its quantitative features appear mildly different than those obtained in previous analyses. In particular, we find a high-energy tail at higher significance than previously reported. This tail is very prominent in the northern hemisphere, and less so in the southern hemisphere. This strongly affects one prominent interpretation of the excess: known millisecond pulsars are incapable of producing this high-energy emission, even in the relatively softer southern hemisphere, and are therefore disfavored as the sole explanation of the GCE. The annihilation of dark matter particles of mass ${40}_{\ensuremath{-}7}^{+10}\text{ }\text{ }\mathrm{GeV}$ (95% CL) to $b$ quarks with a cross-section of $⟨{\ensuremath{\sigma}}_{A}v⟩=1.{4}_{\ensuremath{-}0.3}^{+0.6}\ifmmode\times\else\texttimes\fi{}{10}^{\ensuremath{-}26}\text{ }\text{ }{\mathrm{cm}}^{3}\text{ }{\mathrm{s}}^{\ensuremath{-}1}$ provides a good fit to the excess especially in the relatively cleaner southern sky. Dark matter of the same mass range annihilating to $b$ quarks or heavier dark matter particles annihilating to heavier Standard Model bosons can combine with millisecond pulsars to provide a good fit to the southern hemisphere emission as well, as can a broken power-law spectrum which would be related to recent cosmic-ray burst activity. As part of this paper, we make publicly available all of our templates and the data covariance matrix we have generated to account for systematic uncertainties.