No Arabic abstract
MeV dark matter (DM) particles annihilating or decaying to electron-positron pairs cannot, in principle, be observed via local cosmic-ray (CR) measurements because of the shielding solar magnetic field. In this letter, we take advantage of spacecraft Voyager 1s capacity for detecting interstellar CRs since it crossed the heliopause in 2012. This opens up a new avenue to probe DM in the sub-GeV energy/mass range that we exploit here for the first time. From a complete description of the transport of electrons and positrons at low energy, we derive predictions for both the secondary astrophysical background and the pair production mechanisms relevant to DM annihilation or decay down to the MeV mass range. Interestingly, we show that reacceleration may push positrons up to energies larger than the DM particle mass. We combine the constraints from the Voyager and AMS-02 data to get novel limits covering a very extended DM particle mass range, from MeV to TeV. In the MeV mass range, our limits reach annihilation cross sections of order $langle sigma vrangle sim 10^{-28}{rm cm^3/s}$. An interesting aspect is that these limits barely depend on the details of cosmic-ray propagation in the weak reacceleration case, a configuration which seems to be favored by the most recent boron-to-carbon ($B/C$) data. Though extracted from a completely different and new probe, these bounds have a strength similar to those obtained with the cosmic microwave background --- they are even more stringent for $p$-wave annihilation.
We recently proposed a method to constrain $s$-wave annihilating MeV dark matter from a combination of the Voyager 1 and the AMS-02 data on cosmic-ray electrons and positrons. Voyager 1 actually provides an unprecedented probe of dark matter annihilation to cosmic rays down to $sim 10$ MeV in an energy range where the signal is mostly immune to uncertainties in cosmic-ray propagation. In this article, we derive for the first time new constraints on $p$-wave annihilation down to the MeV mass range using cosmic-ray data. To proceed, we derive a self-consistent velocity distribution for the dark matter across the Milky Way by means of the Eddington inversion technique and its extension to anisotropic systems. As inputs, we consider state-of-the-art Galactic mass models including baryons and constrained on recent kinematic data. They allow for both a cored or a cuspy halo. We then calculate the flux of cosmic-ray electrons and positrons induced by $p$-wave annihilating dark matter and obtain very stringent limits in the MeV mass range, robustly excluding cross sections greater than $sim 10^{-22}{rm cm^3/s}$ (including theoretical uncertainties), about 5 orders of magnitude better than current CMB constraints. This limit assumes that dark matter annihilation is the sole source of cosmic rays and could therefore be made even more stringent when reliable models of astrophysical backgrounds are included.
Some direct detection experiments have recently collected excess events that could be interpreted as a dark matter (DM) signal, pointing to particles in the $sim$10 GeV mass range. We show that scenarios in which DM can self-annihilate with significant couplings to quarks are likely excluded by the cosmic-ray (CR) antiproton data, provided the annihilation is S-wave dominated when DM decouples in the early universe. These limits apply to most of supersymmetric candidates, eg in the minimal supersymmetric standard model (MSSM) and in the next-to-MSSM (NMSSM), and more generally to any thermal DM particle with hadronizing annihilation final states.
Galactic Cosmic-ray (CR) transport parameters are usually constrained by the boron-to-carbon ratio. This procedure is generically plagued with degeneracies between the diffusion coefficient and the vertical extent of the Galactic magnetic halo. The latter is of paramount importance for indirect dark matter (DM) searches, because it fixes the amount of DM annihilation or decay that contributes to the local antimatter CR flux. These degeneracies could be broken by using secondary radioactive species, but the current data still have large error bars, and this method is extremely sensitive to the very local interstellar medium (ISM) properties. Here, we propose to use the low-energy CR positrons in the GeV range as another direct constraint on diffusion models. We show that the PAMELA data disfavor small diffusion halo ($Llesssim 3$ kpc) and large diffusion slope models, and exclude the minimal ({em min}) configuration (Maurin et al. 2001, Donato et al. 2004) widely used in the literature to bracket the uncertainties in the DM signal predictions. This is complementary to indirect constraints (diffuse radio and gamma-ray emissions) and has strong impact on DM searches. Indeed this makes the antiproton constraints more robust while enhancing the discovery/exclusion potential of current and future experiments, like AMS-02 and GAPS, especially in the antiproton and antideuteron channels.
The Galactic positrons, as observed by their annihilation gamma-ray line at 0.511 MeV, are difficult to account for with astrophysical sources. It has been proposed that they are produced instead by dark matter annihilation or decay in the inner Galactic halo. To avoid other constraints, these processes are required to occur invisibly, such that the eventual positron annihilation is the only detectable signal. However, electromagnetic radiative corrections to these processes inevitably produce real gamma rays (``internal bremsstrahlung); this emission violates COMPTEL and EGRET constraints unless the dark matter mass is less than about 20 MeV.
Dark matter annihilation into charged particles is necessarily accompanied by gamma rays, produced via radiative corrections. Internal bremsstrahlung from the final state particles can produce hard gamma rays up to the dark matter mass, with an approximately model-independent spectrum. Focusing on annihilation into electrons, we compute robust upper bounds on the dark matter self annihilation cross section $<sigma_A v >_{e^+e^-}$ using gamma ray data from the Milky Way spanning a wide range of energies, $sim10^{-3} - 10^4$ GeV. We also compute corresponding bounds for the other charged leptons. We make conservative assumptions about the astrophysical inputs, and demonstrate how our derived bounds would be strengthened if stronger assumptions about these inputs are adopted. The fraction of hard gamma rays near the endpoint accompanying annihilation to $e^+e^-$ is only a factor of $alt 10^2$ lower than for annihilation directly to monoenergetic gamma rays. The bound on $<sigma_A v >_{e^+e^-}$ is thus weaker than that for $<sigma_A v >_{gammagamma}$ by this same factor. The upper bounds on the annihilation cross sections to charged leptons are compared with an upper bound on the {it total} annihilation cross section defined by neutrinos.