Sensitivity to sub-GeV dark matter from cosmic-ray scattering with very-high-energy gamma-ray observatories

Huge efforts have been deployed to detect dark matter (DM) in the GeV-TeV mass range involving various detection techniques, and led to strong constraints in the available parameter space. We compute here the sensitivity to sub-GeV DM that can be probed from the inevitable cosmic-ray scattering onto DM particles populating the Milky Way halo. Inelastic scattering of energetic cosmic rays off DM would produce high-energy gamma rays in the final state, providing a new avenue to probe the poorly-constrained so far sub-GeV dark matter mass range. In this work we derive sensitivity forecasts for the inelastic cosmic-ray proton - DM cross section for current and future very-high-energy gamma-ray observatories such as H.E.S.S., LHAASO, CTA and SWGO in the 100 eV to 100 MeV mass range. These inelastic cross section constraints are converted to the elastic proton - DM cross section to highlight further complementarity with cosmological, collider and direct detection searches. The sensitivity computed at 95% confidence level on the elastic cross section reaches ~2× 10⁻³² cm² for a 100 keV DM mass for H.E.S.S.-like and ~7× 10⁻³⁴ cm² for a ~1 keV DM mass for LHAASO. The sensitivity prospects for CTA and a strawman SWGO model reach ~6× 10⁻³⁴ cm² and ~4× 10⁻³⁵ cm², for DM masses of 10 keV and 1 keV, respectively. The sensitivity reach of the gamma-ray observatories considered here enables to probe an uncharted region of the DM mass - cross section parameter space.