Sub-GeV Dark Matter particles upscattered by cosmic rays gain enough kinetic energy to pass the thresholds of large volume detectors on Earth. We then use public Super-Kamiokande and MiniBooNE data to derive a novel limit on the scattering cross section of Dark Matter with electrons that extends down to sub-keV masses, closing a previously allowed wide region of parameter space. We finally discuss search strategies and prospects at existing and planned neutrino facilities.

We point out that current and planned telescopes have the potential of probing annihilating Dark Matter (DM) with a mass of O(100) TeV and beyond. As a target for such searches, we propose models where DM annihilates into lighter mediators, themselves decaying into Standard Model (SM) particles. These models allow to reliably compute the energy spectra of the SM final states, and to naturally evade the unitarity bound on the DM mass. Indeed, long-lived mediators may cause an early matter-dominated phase in the evolution of the Universe and, upon decaying, dilute the density of preexisting relics thus allowing for very large DM masses. We compute this dilution in detail and provide results in a ready-to-use form. Considering for concreteness a model of dark U(1) DM, we then study both dilution and the signals at various high energy telescopes observing gamma rays, neutrinos and charged cosmic rays. This study enriches the physics case of these experiments, and opens a new observational window on heavy new physics sectors.

Dwarf spheroidal galaxies are among the most promising targets for detecting signals of Dark Matter (DM) annihilations. The H.E.S.S. experiment has observed five of these systems for a total of about 130 hours. The data are re-analyzed here, and, in the absence of any detected signals, are interpreted in terms of limits on the DM annihilation cross section. Two scenarios are considered: i) DM annihilation into mono-energetic gamma-rays and ii) DM in the form of pure WIMP multiplets that, annihilating into all electroweak bosons, produce a distinctive gamma-ray spectral shape with a high-energy peak at the DM mass and a lower-energy continuum. For case i), upper limits at 95\% confidence level of about $\langle \sigma v \rangle \lesssim 3 \times 10^{-25}$ cm$^3$ s$^{-1}$ are obtained in the mass range of 400 GeV to 1 TeV. For case ii), the full spectral shape of the models is used and several excluded regions are identified, but the thermal masses of the candidates are not robustly ruled out.

High energy ${\gamma}$-rays are powerful probes in the search for annihilations of dark matter (DM) par- ticles in dense environments. In several DM particle models their annihilation produces characteristic features such as lines, bumps or cut-offs in their energy spectrum. The High Energy Stereoscopic System (H.E.S.S.) of imaging atmospheric Cherenkov telescopes is perfectly suited to search for such features from multi-TeV mass DM particles. The Dwarf Spheroidal Galaxies (dSphs) of the Local Group are the most common satellites of the Milky Way and assumed to be gravitationally bound dominantly by DM, with up to O(10 3 ) times more mass in DM than in visible matter. Over the past decade, several observational campaigns on dwarf satellite galaxies were launched by H.E.S.S. amounting to more than 140 hours of exposure in total. The observations are reviewed here. In the absence of clear signals, the expected spectral and spatial morphologies of signal and background are used to derive constraints on the DM particle annihilation cross- section for particle models producing line-like signals. The combination of the data of all the dwarf galaxies allows a significant improvement in the HESS sensitivity.

Dark matter (DM) charged under a dark U(1) force appears in many extensions of the Standard Model, and has been invoked to explain anomalies in cosmic-ray data, as well as a self-interacting DM candidate. In this paper, we perform a comprehensive phenomenological analysis of such a model, assuming that the DM abundance arises from the thermal freeze-out of the dark interactions. We include, for the first time, bound-state effects both in the DM production and in the indirect detection signals, and quantify their importance for Fermi, AMS, and CMB experiments. We find that DM in the mass range 1 GeV to 100 TeV, annihilating into dark photons of MeV to GeV mass, is in conflict with observations. Instead, DM annihilation into heavier dark photons is viable. We point out that the late decays of multi-GeV dark photons can produce significant entropy and thus dilute the DM density. This can lower considerably the dark coupling needed to obtain the DM abundance, and in turn relax the existing constraints.

We provide CTA sensitivities to Dark Matter (DM) annihilation in $\gamma$-ray lines, from the observation of the Galactic Center (GC) as well as, for the first time, of dwarf Spheroidal galaxies (dSphs). We compare the GC reach with that of dSphs as a function of a putative core radius of the DM distribution, which is itself poorly known. We find that the currently best dSph candidates constitute a more promising target than the GC, for core radii of one to a few kpc. We use the most recent instrument response functions and background estimations by CTA, on top of which we add the diffuse photon component. Our analysis is of particular interest for TeV-scale electroweak multiplets as DM candidates, such as the supersymmetric Wino and the Minimal Dark Matter fiveplet, whose predictions we compare with our projected sensitivities.

We reconsider the model of Minimal Dark Matter (a fermionic, hypercharge-less quintuplet of the EW interactions) and compute its gamma ray signatures. We compare them with a number of gamma ray probes: the galactic halo diffuse measurements, the galactic center line searches and recent dwarf galaxies observations. We find that the original minimal model, whose mass is fixed at 9.4 TeV by the relic abundance requirement, is constrained by the line searches from the Galactic Center: it is ruled out if the Milky Way possesses a cuspy profile such as NFW but it is still allowed if it has a cored one. Observations of dwarf spheroidal galaxies are also relevant (in particular searches for lines), and ongoing astrophysical progresses on these systems have the potential to eventually rule out the model. We also explore a wider mass range, which applies to the case in which the relic abundance requirement is relaxed. Most of our results can be safely extended to the larger class of multi-TeV WIMP DM annihilating into massive gauge bosons.

We provide ingredients and recipes for computing signals of TeV-scale Dark Matter annihilations and decays in the Galaxy and beyond. For each DM channel, we present the energy spectra of electrons and positrons, antiprotons, antideuterons, gamma rays, neutrinos and antineutrinos e, mu, tau at production, computed by high-statistics simulations. We estimate the Monte Carlo uncertainty by comparing the results yielded by the Pythia and Herwig event generators. We then provide the propagation functions for charged particles in the Galaxy, for several DM distribution profiles and sets of propagation parameters. Propagation of electrons and positrons is performed with an improved semi-analytic method that takes into account position-dependent energy losses in the Milky Way. Using such propagation functions, we compute the energy spectra of electrons and positrons, antiprotons and antideuterons at the location of the Earth. We then present the gamma ray fluxes, both from prompt emission and from Inverse Compton scattering in the galactic halo. Finally, we provide the spectra of extragalactic gamma rays. All results are available in numerical form and ready to be consumed.

The computation of the energy spectra of Standard Model particles originated from the annihilation/decay of dark matter particles is of primary importance in indirect searches of dark matter. We compute how the inclusion of electroweak corrections significantly alter such spectra when the mass M of dark matter particles is larger than the electroweak scale: soft electroweak gauge bosons are copiously radiated opening new channels in the final states which otherwise would be forbidden if such corrections are neglected. All stable particles are therefore present in the final spectrum, independently of the primary channel of dark matter annihilation/decay. Such corrections are model independent.