A Model Independent Method to Study Dark Matter induced Leptons and Gamma rays

By using recent data, we directly determine the dark matter (DM) induced $e^\pm$ spectrum at the source from experimental measurements at the earth, without reference to specific particle physics models. The DM induced gamma rays emitted via inverse Compton scattering are then obtained in a model independent way. However the results depend on the choice of the astrophysical $e^\pm$ background, which is not reliably known. Nevertheless, we calculate, as an illustration, the fluxes of gamma rays from the Fornax cluster in the decaying DM scenario with various astrophysical $e^\pm$ backgrounds. Without any assumptions on details of the DM model, the predictions turn out to be either in disagreement with or only marginally below the upper limits measured recently by the Fermi-LAT Collaboration. In addition, these DM induced ICS gamma rays in the GeV range are shown to be almost independent of choices of cosmic ray propagation model and of DM density profile, when a given astrophysical $e^\pm$ background is assumed. This provides a strong constraint on decaying DM scenario as the gamma rays may be produced in other processes besides inverse Compton scattering, such as the bremsstrahlung and neutral pion decays.

💡 Research Summary

The authors present a novel, model‑independent framework for extracting the dark‑matter‑induced electron/positron source spectrum directly from Earth‑based cosmic‑ray measurements, and then using that spectrum to predict inverse‑Compton‑scattering (ICS) gamma‑ray emission from galaxy clusters. Traditional indirect‑detection studies usually start from a specific particle‑physics model, fixing the dark‑matter mass, decay or annihilation channels, and branching ratios, then compute the source spectrum and propagate it through the Galaxy. In contrast, this work solves the diffusion‑loss equation analytically, recognizing that the solution for the DM‑induced electron density (f_{\rm DM}^e(E,\mathbf r)) can be written as a Volterra integral over the unknown source term (X(E)). By applying a Laplace transform and its inverse, they derive an explicit expression (Eq. 5) that reconstructs (X(E)) solely from the measured flux at the solar position, together with the diffusion parameters (K(E)=K_0(E/{\rm GeV})^\alpha) and the energy‑loss time (\tau_E). The diffusion parameters are taken from three benchmark propagation models (MIN, MED, MAX), which span the range of plausible positron fluxes.

Having obtained (X(E)) without any assumption about the dark‑matter particle, the authors feed it into the standard formula for ICS gamma‑ray production, focusing on the Fornax galaxy cluster as a point‑like source. In Fornax the dominant target photons are the cosmic‑microwave‑background (CMB); starlight and dust contributions are negligible. The gamma‑ray flux is computed by integrating the product of (X(E)), the CMB photon density, and the Klein‑Nishina cross‑section over electron energy. The authors explore several dark‑matter density profiles (NFW, Einasto, Isothermal) and find that, because the relevant electrons have energies above ∼500 GeV and thus travel only short distances before losing energy, the choice of halo profile has essentially no impact on the predicted gamma flux. Propagation model variations affect the reconstructed (X(E)) only below ∼300 GeV, a regime irrelevant for the GeV‑range gamma rays considered.

The key phenomenological result is a comparison with the 18‑month Fermi‑LAT upper limits on gamma‑ray emission from nearby clusters. Using a conventional astrophysical electron/positron background (“model 0”) with a normalization factor (N=0.8), the reconstructed source spectrum yields an ICS gamma‑ray component that either exceeds or lies very close to the Fermi‑LAT limits in the 1–10 GeV band. This already places strong tension on decaying‑dark‑matter explanations of the PAMELA/Fermi positron excesses. If a harder background (normalization (N=1)) is adopted, the low‑energy part of the electron spectrum becomes unphysical (negative flux), but the high‑energy part—and consequently the predicted gamma flux—remains essentially unchanged, leading to the same conclusion.

The analysis highlights that the dominant uncertainty stems from the astrophysical electron/positron background, which is not well constrained below a few hundred GeV. Nevertheless, because the gamma‑ray prediction depends only on the high‑energy tail of the source spectrum, the method provides a robust, model‑independent bound on decaying dark matter. The authors note that the same technique can be applied to annihilating dark‑matter scenarios, and that future, more precise measurements of cluster gamma rays and of the local cosmic‑ray background will sharpen these constraints. In summary, the paper demonstrates that, even without specifying any particle‑physics details, current cosmic‑ray and gamma‑ray data already severely limit the viability of decaying dark‑matter models that aim to explain the observed positron excess.