Using γ-ray data from the Fermi Large Area Telescope, various groups have identified a clear excess emission in the Inner Galaxy, at energies around a few GeV. This excess resembles remarkably well a signal from dark-matter annihilation. One of the most compelling astrophysical interpretations is that the excess is caused by the combined effect of a previously undetected population of dim γ-ray sources. Because of their spectral similarity, the best candidates are millisecond pulsars. Here, we search for this hypothetical source population, using a novel approach based on wavelet decomposition of the γ-ray sky and the statistics of Gaussian random fields. Using almost seven years of Fermi -LAT data, we detect a clustering of photons as predicted for the hypothetical population of millisecond pulsar, with a statistical significance of 10.0σ. For plausible values of the luminosity function, this population explains 100% of the observed excess emission. We argue that other extragalactic or Galactic sources, a mismodeling of Galactic diffuse emission, or the thick-disk population of pulsars are unlikely to account for this observation.

An anomalous emission component at energies of a few GeV and located towards the inner Galaxy is present in the Fermi -LAT data. It is known as the Fermi -LAT GeV excess. Using almost 8 years of data we reanalyze the characteristics of this excess with SkyFACT, a novel tool that combines image reconstruction with template fitting techniques. We find that an emission profile that traces stellar mass in the boxy and nuclear bulge provides the best description of the excess emission, providing strong circumstantial evidence that the excess is due to a stellar source population in the Galactic bulge. We find a luminosity to stellar mass ratio of (2.1 ± 0.2) × 10 27 erg s −1 M −1 for the boxy bulge, and of (1.4 ± 0.6) × 10 27 erg s −1 M −1 for the nuclear bulge. Stellar mass related templates are preferred over conventional DM profiles with high statistical significance.

Abstract. Over the past decade, extensive studies have been undertaken to search for photon signals from dark matter annihilation or decay for dark matter particle masses above ∼ 1 GeV. However, due to the lacking sensitivity of current experiments at MeV-GeV energies, sometimes dubbed the 'MeV gap', dark matter models with MeV to sub-GeV particle masses have received little attention so far. Various proposed MeV missions (like, e.g., e-ASTROGAM or AMEGO) are aimed at closing this gap in the mid-or long-term future. This, and the absence of clear dark matter signals in the GeV-TeV range, makes it relevant to carefully reconsider the expected experimental instrumental sensitivities in this mass range. The most common two-body annihilation channels for sub-GeV dark matter are to neutrinos, electrons, pions or directly to photons. Among these, only the electron channel has been extensively studied, and almost exclusively in the context of the 511 keV line. In this work, we study the prospects for detecting MeV dark matter annihilation in general in future MeV missions, using e-ASTROGAM as reference, and focusing on dark matter masses in the range 1 MeV-3 GeV. In the case of leptonic annihilation, we emphasise the importance of the often overlooked bremsstrahlung and in-flight annihilation spectral features, which in many cases provide the dominant gamma-ray signal in this regime.

In the cold dark matter paradigm, structures form hierarchically, implying that large structures contain smaller substructures. These subhalos will enhance signatures of dark matter annihilation such as gamma rays. In the literature, typical estimates of this boost factor assume a concentration-mass relation for field halos, to calculate the luminosity of subhalos. However, since subhalos accreted in the gravitational potential of their host lose mass through tidal stripping and dynamical friction, they have a quite characteristic density profile, different from that of the field halos of the same mass. In this work, we quantify the effect of tidal stripping on the boost factor, by developing a semianalytic model that combines the mass-accretion history of both the host and subhalos as well as subhalo accretion rates. We find that when subhalo luminosities are treated consistently the boost factor increases by a factor 2-5, compared to the typical calculation assuming a field-halo concentration. This holds for host halos ranging from subgalaxy to cluster masses and is independent of the subhalo mass function or specific concentration-mass relation. The results are particularly relevant for indirect dark matter searches in the extragalactic gamma-ray sky.