An in-depth study of Gamma rays from the Starburst Galaxy M 82 with VERITAS

Assuming Galactic cosmic rays originate in supernovae and the winds of massive stars, starburst galaxies should produce very-high-energy (VHE; E$>$100 GeV) gamma-ray emission via the interaction of their copious quantities of cosmic rays with the large reservoirs of dense gas within the galaxies. Such VHE emission was detected by VERITAS from the starburst galaxy M 82 in 2008-09. An extensive, multi-year campaign followed these initial observations, yielding a total of 254 h of good quality VERITAS data on M 82. Leveraging modern analysis techniques and the larger exposure, these VERITAS data show a more statistically significant VHE signal ($\sim$6.5 standard deviations ($\sigma$)). The corresponding photon spectrum is well fit by a power law ($\Gamma = 2.3 \pm 0.3_{stat} \pm0.2_{sys}$) and the observed integral flux is F($>$450 GeV) = $(3.2 \pm0.6_{stat} \pm 0.6_{sys}) \times 10^{-13}~\mathrm{cm^{-2}~s}^{-1}$, or $\sim$0.4\% of the Crab Nebula flux above the same energy threshold. The improved VERITAS measurements, when combined with various multi-wavelength data, enable modeling of the underlying emission and transport processes. A purely leptonic scenario is found to be a poor representation of the gamma-ray spectral energy distribution (SED). A lepto-hadronic scenario with cosmic rays following a power-law spectrum in momentum (index $s\simeq 2.25$), and with significant bremsstrahlung below $1$~GeV, provides a good match to the observed SED. The synchrotron emission from the secondary electrons indicates that efficient non-radiative losses of cosmic-ray electrons may be related to advective escape from the starburst core.