# Loopy Constraints on Leptophilic Dark Matter and Internal Bremsstrahlung

###### Abstract

A sharp and spatially extended peak in an astrophysical gamma ray spectrum would provide very strong evidence for the existence of dark matter (DM), given that there are no known astrophysical processes that could mimic such a signal. From the particle physics perspective, perhaps the simplest explanation for a gamma ray peak is internal bremsstrahlung in DM annihilation through a charged -channel mediator close in mass to the DM particle . Since DM annihilation to quarks is already tightly constrained in this scenario, we focus here on the leptophilic case. We compute the electromagnetic anapole and dipole moments that DM acquires at 1-loop, and we find an interesting enhancement of these moments if the DM particle and the mediator are close in mass. We constrain the DM anapole and dipole moments using direct detection data, and then translate these limits into bounds on the DM annihilation cross section. Our bounds are highly competitive with those from astrophysical gamma ray searches. In the second part of the paper, we derive complementary constraints on internal bremsstrahlung in DM annihilation using LEP mono-photon data, measurements of the anomalous magnetic moments of the electron and the muon, and searches for lepton flavor violation. We also comment on the impact of the internal bremsstrahlung scenario on the hyperfine splitting of true muonium.

## I Introduction

One of the cleanest signatures in indirect dark matter (DM) searches are peaks
in the cosmic gamma ray spectrum from the Galactic Center or other regions of
high DM density. On the one hand, there are no known astrophysical sources
that could mimic such a signal.^{1}^{1}1The authors of
Ref. Aharonian et al. (2012) show that a particular composition of a pulsar
wind nebula could generate a peaked gamma ray signal, but an observation of a
peak at the same energy in different regions of the galaxy would rule out this
possibility. On the other hand, gamma ray observatories are making tremendous
progress in terms of statistics, resolution and control of systematic
uncertainties.

From the particle physics point of view, peaks in the gamma ray spectrum can originate from DM annihilation or decay to two photons, a photon and a boson, or a photon and a Higgs boson. However, since DM is electrically neutral, these processes can only happen at the 1-loop level, making it likely that DM is first discovered in other annihilation or decay channels. There is, however, a class of models where the first experimental hint for DM is a gamma ray peak. Namely, this can happen in models where DM annihilates via a charged -channel mediator, so that a photon can be emitted from the mediator, see Fig. 1. This process is called virtual internal bremsstrahlung (VIB) Bringmann et al. (2008); Bell et al. (2011); Bringmann et al. (2012). If the mediator mass and the DM mass are close to each other, the resulting photon energy is strongly peaked (see Fig. 2) and can yield a line-like gamma ray signal if the width of the peak is below the detector resolution.

Of particular interest in this context are models in which the DM couples preferentially to leptons. These leptophilic dark matter (LDM) models are motivated by the fact that DM couplings to quarks are strongly constrained by gamma ray emission from dwarf galaxies Geringer-Sameth and Koushiappas (2011); Ackermann et al. (2013), by direct detection bounds Hisano et al. (2011); Garny et al. (2012, 2013), and by LHC searches Goodman et al. (2011); Bai et al. (2010); Fox et al. (2012); Lin et al. (2013); The CMS collaboration (2013); The ATLAS collaboration (2012); Aad et al. (2013); Boyd (2013); Richman (2013). Additional motivation could be provided by various cosmic ray anomalies. For instance, attempts to explain the cosmic ray positron excess observed by PAMELA Adriani et al. (2009, 2013), Fermi-LAT Ackermann et al. (2012a) and AMS-02 Aguilar et al. (2013) in terms of DM annihilation typically require a leptophilic DM model Cirelli et al. (2009); Donato et al. (2009); Nardi et al. (2009); Bertone et al. (2009); Fox and Poppitz (2009); Evoli et al. (2011); Kopp (2013); Cholis and Hooper (2013); Bergstrom et al. (2013); Ibarra et al. (2013) in order not to exceed the measured antiproton flux Adriani et al. (2010); Bartoli et al. (2012). Finally, it is intriguing that the possible anomalies in the gamma ray signal from the Galactic Center Hooper and Goodenough (2011); Hooper and Linden (2011); Hooper (2012), in the gamma ray emission from the Fermi Bubbles Su et al. (2010); Hooper and Slatyer (2013); Huang et al. (2013a, b), and in radio signals from filamentary structures in the inner galaxy Linden et al. (2011) could be explained in leptophilic DM models. (Note, however, that some of them can also be understood if dark matter annihilates to final states.) Direct detection constraints on leptophilic DM have been studied in Kopp et al. (2009); Essig et al. (2012a, b).

In the present paper, we derive new constraints on leptophilic DM, and we translate these constraints into bounds on the cross section for internal bremsstrahlung. We also discuss the prospects for probing the parameter space of leptophilic DM even further with future experiments. We work in a simplified model which augments the Standard Model (SM) by a fermionic DM candidate and a charged scalar mediator , with a coupling of the form , where is a charged lepton field. This effective scenario can be realized in supersymmetry (SUSY) (see for instance Bringmann et al. (2012)), where could be identified with the lightest neutralino, and would be a slepton. It also applies to certain radiative neutrino mass models, whose direct detection phenomenology has been discussed in Schmidt et al. (2012). A simplified framework of the form used here has been employed, for instance, to explain an anomalous line-like feature at GeV in the Fermi-LAT gamma ray data Bringmann et al. (2012); Ackermann et al. (2012b); The Fermi-LAT collaboration (2013). Even though the statistical significance of this feature is not yet convincing The Fermi-LAT collaboration (2013), and there are (inconclusive) indications that poorly understood systematic effects may play a role Boyarsky et al. (2012); Whiteson (2012); Hektor et al. (2012); Finkbeiner et al. (2012); Whiteson (2013), it demonstrates the relevance of internal bremsstrahlung signatures as considered here if anomalous peaked features are found in future gamma ray observations.

Our starting point is the observation that even in leptophilic models, loop processes endow the DM with nonzero electromagnetic moments, which in turn allow it to interact in direct detection experiments. If DM is a Majorana fermion, only an anapole moment is generated Radescu (1985); Kayser and Goldhaber (1983), while for Dirac fermions, also a magnetic dipole moment can exist. DM with anapole interactions has been studied previously in Ho and Scherrer (2012); Gresham and Zurek (2013); Gao et al. (2013); Del Nobile et al. (2014) using an effective field theory framework, and DM with magnetic dipole moments has been investigated in Heo (2010); Masso et al. (2009); Schmidt et al. (2012); Del Nobile et al. (2012); Weiner and Yavin (2012a); Barger et al. (2012a); Gresham and Zurek (2013). The importance of loop processes even for hadrophilic DM has been studied in the context of LHC searches in Haisch and Kahlhoefer (2013).

Loop processes involving DM particles can also modify electromagnetic properties of leptons, in particular their anomalous magnetic moments and the energy levels of dilepton systems such as positronium and muonium. In the most general case, also lepton flavor violation could be induced by DM loops. Finally, if DM couples to electrons, it can be directly produced at LEP or at a future linear collider, allowing us to derive constraints from searches for mono-photons plus missing energy.

The paper is organized as follows. In Sec. II, we introduce the simplified model framework which we will use in the rest of the paper, we establish its connection to supersymmetric scenarios, and we review the expected indirect detection (internal bremsstrahlung) signals from DM annihilation in this model. In Sec. III, we compute the electromagnetic form factors of DM and the resulting direct detection cross sections. We compare these to LUX Akerib et al. (2013) and XENON100 Aprile et al. (2012) data, and to the expected future sensitivity of XENON1T and LUX-ZEPLIN to derive constraints. We translate these constraints into limits on the intensity of possible internal bremsstrahlung signals. To illustrate the strength of direct detection limits, we show that for flavor-universal DM couplings to leptons, the explanation of the aforementioned 135 GeV feature in the Fermi-LAT data Bringmann et al. (2012); Ackermann et al. (2012b); The Fermi-LAT collaboration (2013) in terms of internal bremsstrahlung is severely constrained. We then investigate in Sec. IV the complementary constraints from collider data, and in Sec. V the bounds from measurements of the anomalous magnetic moment of the muon and the electron, from searches for lepton flavor violation, and from possible future experiments on true muonium spectroscopy. We summarize our findings and conclude in Sec. VI.

## Ii Internal bremsstrahlung in simplified models

### ii.1 The simplest model

The simplest theoretical models that feature internal bremsstrahlung in DM annihilation extend the Standard Model by a neutral DM candidate and a charged mediator Bringmann et al. (2012). can be either a Majorana fermion (as in most supersymmetric theories) or a Dirac fermion (as for example in supersymmetric theories with preserved -symmetry Chun et al. (2010); Buckley et al. (2013)). As explained above, we are mostly interested in leptophilic models because DM couplings to quarks are already tightly constrained. In the simplest case, we thus start with the interaction Lagrangian

(1) |

where denotes the photon field, is the fermionic DM candidate, is the unit electric charge, is a SM lepton field, is the right-handed chiral projection operator, and is the Yukawa coupling constant of the DM–lepton interaction. Unless indicated otherwise, we assume to be a Majorana fermion. Note that we have omitted couplings to left-handed leptons here which are more strongly constrained (though not ruled out) by collider searches and electroweak precision test Liu et al. (2013). We also do not consider the scalar potential for since these terms are irrelevant to our discussion. Finally, we disregard the vertex from the kinetic term of because it is higher order in the coupling constant and will thus be phenomenologically negligible.

The simplified model (1) has been studied previously for instance in Cao et al. (2009), and it has been shown in Bringmann et al. (2012) that the model could explain the 135 GeV feature in the Fermi-LAT data. The fit from Bringmann et al. (2012) results in a preferred DM mass of and an annihilation cross section . (Here, is the relative velocity of the two annihilating DM particles, and the average is taken over .)

The interactions in eq. (1) lead to annihilation of DM particles into pairs of SM leptons via -channel exchange of the charged scalar . This process can be decomposed into an -wave part and a -wave part, the latter of which can usually be neglected because it is suppressed by the square of the small velocity km/s of DM particles in the Milky Way. The -wave contribution is unsuppressed for Dirac DM, while for Majorana DM, it is helicity-suppressed by the small mass of the final state lepton Bringmann et al. (2012). This can be understood by noting that DM annihilation through the Yukawa interaction in Eq. (1) produces two leptons of the same chirality. For Majorana DM, however, Pauli blocking in the initial state requires the incoming DM particles to have opposite spin. Angular momentum conservation therefore requires a mass insertion on one of the final state lepton lines. Thus, for Majorana DM, higher order annihilation processes become important, in particular the process , with two charged leptons and a photon in the final state (see Fig. 1). Since the photon carries away one unit of angular momentum, it can lift the helicity suppression, see for instance Bergstrom (1989); Barger et al. (2009).

A helicity suppression of 2-body DM annihilation compared to the 3-body internal bremsstrahlung process exists also in models where the scalar mediator is replaced by a vector particle Barger et al. (2012b) and in models with scalar DM and fermionic mediators Barger et al. (2012b); Toma (2013); Giacchino et al. (2013). We will not consider these possibilities here, but will instead focus on the scenario from Eq. (1) as a representative for all internal bremsstrahlung models.

(a) | (b) |

If the mediator mass and the DM mass are nearly degenerate, the emission of an internal bremsstrahlung photon (first and fourth diagram in Fig. 1) is strongly peaked if the photon energy gets close to . The reason is that, in this case, one of the final state leptons is very soft, and one of the propagators gets close to the mass shell. In other words, internal bremsstrahlung with and can be viewed as DM annihilation into a lepton and a photon, with the emission of a soft lepton as a form of initial state radiation. While the spectral peak is thus due to internal bremsstrahlung only, it is important to take into account also the final state radiation diagrams to guarantee gauge invariance of the process. Note that, in contrast to gamma ray lines from DM annihilation to , or , the peaked signal from internal bremsstrahlung is not loop-suppressed, hence the cross section can be sizeable. The differential three-body cross section for in the case of Majorana DM has the following form Bringmann et al. (2012)

(2) |

with the electromagnetic fine structure constant , the number of final state lepton flavors , and with the definitions and . In Eq. (2), we have neglected the lepton mass and the DM velocity . is plotted in Fig. 2 for different values of . It is clear that, in order to have a distinct peak, a small degeneracy parameter is necessary. Integrating over , we immediately obtain also the full cross section Bringmann et al. (2012)

(3) |

Here, is the dilogarithm function.

The approximate expression for the relic density of Majorana DM in our toy model is Bringmann et al. (2012)

(4) |

for . For smaller , is smaller by an factor due to coannihilations Bringmann et al. (2012) (see Giacchino et al. (2013) Fig. VII, for a quantitative estimate of the effect of co-annihilations). We see that in the interesting parameter range , GeV, the model (1) naturally predicts a relic density comparable to the observed value . Here, the quoted upper limit on is taken from Planck Ade et al. (2013), whereas for the lower limit, we conservatively use the WMAP value Hinshaw et al. (2013). We thus account in a qualitative way for the uncertainty in from the yet unresolved tension between different measurements of the Hubble constant .

For Dirac DM, is not a
well-defined quantity in the limit due to infrared divergences
in the phase space region where the photon is soft or collinear with one of the
leptons. For , we can evaluate numerically, see
Fig. 2.^{2}^{2}2We have checked that the logarithms appearing
in the expression for
are sufficiently small for a perturbative treatment to be approximately valid.
We find that the spectrum is entirely dominated by
final state radiation and no internal bremsstrahlung peak is discernable at
. This means in particular that no sharp spectral
features are expected for Dirac DM. In the following, we will therefore use the
two-body annihilation cross section

(5) |

as a figure of merit for indirect detection of Dirac DM.

### ii.2 Extended models and connection to supersymmetry

A natural realization of scenario (1) is provided by the leptonic sector of supersymmetric extensions of the SM. There, the mediator is the lightest slepton and the DM candidate is the lightest neutralino, which is given in terms of its bino (), wino () and higgsino (, ) components as . Here, are elements of the neutralino mixing matrix. The next-to-lightest slepton, as well as the squarks, are assumed to be much heavier than . In the MSSM, the Yukawa coupling can be written in terms of the unit electric charge , the Weinberg angle , and the neutralino mixing matrix element as Haber and Kane (1985)

(6) |

If, instead of Eq. (1), we were considering couplings to left handed leptons and their corresponding sleptons, the Yukawa coupling in the MSSM would be given by . Since for conventional mechanisms of supersymmetry breaking, slepton masses of one chirality tend to be similar, we will also generalize (1) to include all three lepton flavors and slepton flavors of one chirality, where , , :

(7) |

Finally, we will also consider a more general model (which cannot be realized in the MSSM), in which couplings to both left-handed and right-handed fermions are included, and couplings are allowed to be flavor off-diagonal. The Lagrangian for this generalized toy model is

(8) |

Here, are the Yukawa matrices, and are the mass eigenstates of the scalar mediators, of which an arbitrary nmber could exist. The index runs over , , , while runs over all mass eigenstates.

Since our main motivation is the possibility of observing internal bremsstrahlung signals in future gamma ray observations, we will mostly focus on the case where the mass scale of the mediator(s), , is similar to the DM mass. The reason is that in this case the photon spectrum from internal bremsstrahlung is strongly peaked. Note that models with are notoriously difficult to probe at colliders because the charged leptons produced in slepton decays are very soft. In the supersymmetric context, a model with nearly degenerate neutralino and slepton masses has been studied with a different goal in Konishi et al. (2013). We will comment on this model also at the end of Sec. VI.

## Iii Electromagnetic form factors of dark matter and direct detection constraints

We now establish the connection between indirect gamma ray signatures of DM in our toy model and direct laboratory searches on Earth. Connecting the final state fermion lines in the internal bremsstrahlung and final state radiation diagrams from Fig. 1, we obtain an effective vertex coupling the DM particle to the photon through loops of the form shown in Fig. 3. At dimension 5 and 6, the most general form of this effective interaction for a neutral fermion is Fukushima and Kumar (2013)

(9) |

where is the magnetic dipole moment, is the electric dipole moment, and is the anapole moment. For Majorana DM, only the anapole term is nonzero Radescu (1985); Kayser and Goldhaber (1983), as can be seen by using the fact that a Majorana field is invariant under the charge conjugation operation , i.e. . Applying this identity to the fermion fields in Eq. (9), it is straightforward to show that the magnetic and electric dipole terms vanish.

Note that establishing a similar connection between DM annihilation and loop-induced electromagnetic form factors is also possible in internal bremsstrahlung models with scalar DM and fermionic mediators, or with Majorana DM and vector mediators Barger et al. (2012b); Frandsen et al. (2013). We have seen above that these scenarios are phenomenologically as interesting as our model with Majorana DM and a scalar mediator because internal bremsstrahlung dominates over DM annihilation to 2-body final states in all of them. The connection between gamma ray lines from DM annihilation and direct detection signals has been made also for models with loop-induced DM annihilation to photons in Frandsen et al. (2012).

### iii.1 One Loop Contribution to the Electromagnetic Moments

We will now compute the loop induced electromagnetic interactions for the DM particles in our toy model Eq. (1).

#### iii.1.1 Anapole moment for Majorana fermions

(a) | (b) |

We begin by evaluating the diagrams in Fig. 3 to obtain the anapole form factor in Eq. (9) for Majorana DM. For negligible 4-momentum transfer we find

() | (10) |

with , . Taking into account the behavior of the function when its argument approaches 1, it is easy to see that for , the anapole moment diverges logarithmically as . This behavior can be qualitatively understood by noting that, if and , all three propagators in the loops of Fig. 3 can be close to the mass shell simultaneously. In the limit , on the other hand, the leading term in is proportional to . Note that in this limit, the expression in Eq. (10) requires analytic continuation of the function into the complex plane. The dependence of on the degeneracy parameter is shown in Fig. 4 (a) for and GeV.

If , a case that is relevant for instance in DM–nucleus scattering through loops containing electrons, the approximation underlying Eq. (10) is not applicable. In this case, it is instead convenient to set and keep only to the leading term in , which leads to

(11) |

At very small or , one may wonder whether a calculation at fixed order in perturbation theory is still valid. However, in the case of interest to us, namely , the divergent logarithms in Eqs. (10) and (11) are at most of order 10 even for DM couplings to electrons.

#### iii.1.2 Dipole moment for Dirac fermions

If is a Dirac fermion rather than a Majorana particle, only the two diagrams on the left in Fig. 3 exist. They generate an anapole moment that is half as large as the one for Majorana DM, Eq. (10), and a magnetic dipole moment given by

(12) |

for . The dipole moment will turn out to be numerically much more important than the anapole moment in scattering processes involving Dirac DM. If is neglected compared to , i.e. , Eq. (12) simplifies to

(13) |

Note that, unlike the anapole moment , the dipole moment is not divergent for . For , the leading term in is proportional to . The behavior of as a function of is shown in Fig. 4 (b).

### iii.2 Direct detection signals

In this section we will discuss the experimental limits on dark matter scattering through anapole and magnetic dipole interactions. This has been done previously at the effective field theory level for instance in Refs. Sigurdson et al. (2004); Masso et al. (2009); Barger et al. (2011); Fitzpatrick and Zurek (2010); Banks et al. (2010); Del Nobile et al. (2012); Weiner and Yavin (2012b); Ho and Scherrer (2012); Del Nobile et al. (2013); Gresham and Zurek (2013); Del Nobile et al. (2014). Here, we carry out a similar analysis using the latest LUX Akerib et al. (2013) and XENON100 Aprile et al. (2012) data, and we then translate the resulting constraints into new limits on the expected indirect detection signals in our toy model. Since the differential DM–nucleus scattering cross section (where is the nuclear recoil energy) for anapole and dipole interactions differs from the conventional spin-independent or spin-dependent scenarios, we cannot directly use the published exclusion limits from LUX and XENON100, but instead have to fit the data at the event level. We do this by using a framework developed in Refs. Kopp et al. (2009, 2010, 2012), which we have extended by including LUX data and by implementing anapole and dipole interactions.

(a) | (b) |
---|

The differential cross section for DM–nucleus scattering through an anapole interaction is (cf. also Ho and Scherrer (2012); Del Nobile et al. (2014))

(14) |

while for dipole interactions we have Pospelov and ter Veldhuis (2000); Chang et al. (2010); Barger et al. (2011); Banks et al. (2010)

(15) |

In both equations, the first line corresponds to scattering on the nuclear
charge , while the second line describes scattering on the nuclear dipole
moment .^{3}^{3}3Note that the contributions from the nuclear charge and
from the nuclear dipole moment must be separated carefully. For instance, a
naive calculation involving the standard QED vertex for the nucleus would
correctly describe DM–charge scattering, but the contribution from DM–dipole
scattering would be correct only for a truly pointlike nucleus with magnetic
dipole moment . Here, instead, this spurious DM–dipole
scattering term must be subtracted out and replaced by the correct term for
scattering on the dipole moments of extended nuclei (second line of
Eqs. (14) and (15)).
The nuclear mass is denoted by , and is the velocity of
the incoming DM particle. We have also included the nuclear charge form factor
and the spin form factor . We parametrize
as Jungman et al. (1996) , where , fm, , fm (with the
nuclear mass number ). For , we use Banks et al. (2010)
for and , and otherwise. Here, . Note that
nuclear dipole moments are subdominant in many target materials, including
xenon, which we mostly focus on in this paper. The contribution from the nuclear
dipole moment may be comparable to the contribution from the nuclear charge for
instance in fluorine, sodium and iodine Chang et al. (2010). Note that
Eq. (14) can be integrated over to yield a total cross
section, while Eq. (15) has an infrared divergence, which
makes the total cross section for dipole interactions an ill-defined
quantity.

The differential DM–nucleus scattering rate per unit target mass is given by

(16) |

where GeV/cm is the local DM density, is the minimal DM velocity required to yield a recoil energy , is the reduced mass of the DM–nucleus system, and is the DM velocity distribution in the rest frame of the detector. We obtain by a Galilean transformation of the DM velocity distribution in the Milky Way rest frame, . For the latter, in turn, we assume the conventional Maxwell-Boltzmann form with a smooth cutoff, , with velocity dispersion km/s and escape velocity km/s. We expect the dependence of our results on this choice of velocity profile to be similar to what was found for DM scattering through contact interactions in the literature, see for instance McCabe (2010); Farina et al. (2011); Green (2012); Fairbairn et al. (2012).

In Fig. 5, we compare the differential reaction rates for anapole, dipole and spin-independent contact interactions, both with and without including nuclear form factor and detector effects. For easier comparison, all rates are normalized to a total rate of 1 event above 10 keV per kg per day before taking into account nuclear form factor and detector effects. We see that anapole and contact interactions lead to similar event spectra, while dipole interactions are strongly enhanced at low energies due to the dependence of the first term in Eq. (15). The nuclear form factor leads to a suppression of at higher energies. Note that at low energies, the scattering rate remains sizeable down to few keV even because such low energy events can occasionally produce a detectable number of photoelectrons due to Poisson statistics.

We conclude that with sufficient statistical power direct detection experiments could relatively easily distinguish dipole interactions from other interaction structures, while discriminating between anapole and contact interaction is challenging.

In the absence of a signal, we next derive limits on the anapole moment , the dipole moment and the total DM–nucleon scattering cross section for contact interactions, .

### iii.3 Constraints from direct detection data

(a) | (b) |
---|

In Fig. 6 we show the constraints on the anapole and magnetic dipole moments of dark matter from 85.3 days of LUX data Akerib et al. (2013) and from 225 days of XENON100 data Aprile et al. (2012). For the statistical analysis, we have used Yellin’s maximum gap method Yellin (2002). The code employed to derive limits has been developed in Kopp et al. (2009, 2010, 2012), and we have checked that it reproduces the XENON100 and LUX limits on standard spin-independent DM–nucleus scattering to very good accuracy. Note that the qualitative shape of the exclusion curves is similar to the well-known exclusion limit for scattering through contact interactions. At low DM mass, the loss in sensitivity is slightly less steep for dipole interactions due to the enhancement of the scattering rate at low energies (see Fig. 5).

We now derive our main results by translating the LUX constraint on the anapole moment from Fig. 6 (a) into a constraint on the annihilation cross section into two charged leptons plus an internal bremsstrahlung photon using Eqs. (10) and (3). Similarly, we convert the LUX limits on the dipole moment of Dirac DM from Fig. 6 (b) into bounds on the DM annihilation cross section into two charged leptons, using Eqs. (12) and (5). Note that the total cross section for the 3-body final state is ill-defined in the Dirac case due to infrared divergences. Moreover, annihilation into is subdominant for Dirac DM.

Our results are shown in Fig. 7 (a), (b) and (c) for Majorana DM, and in Fig. 7 (d) for Dirac DM. Figs. 7 (a) and (d) are for couplings to only one lepton species , while (b) and (c) are for flavor-universal couplings.

For Majorana DM, Fig. 7 (a) clearly reflects the increase in the anapole moment for small , which here translates into stronger limits on the model parameters and on for coupling to electrons than for coupling to or . We also clearly see the effect of degenerate and : for close to unity, the anapole moment is significantly larger than for well separated and (see Eq. (10) and Fig. 4). Comparing to the preferred parameter region from the gamma ray line search in Bringmann et al. (2012), we find that this region is still marginally compatible with direct detection constraints if . For , it is disfavored at the confidence level if DM has couplings to electrons and at the confidence level if DM couples only to muons. Comparing to the cross sections required for thermal relic DM (horizontal blue line in Fig. 7 (a)), we see that direct detection limits are just starting to probe this region. Note that our estimates for the thermal relic cross section are based on Eq. (4). They do not include the effect of co-annihilations Bringmann et al. (2012), which would move the thermal relic cross section to smaller values. Note also that our perturbative calculations become inaccurate close to the gray regions in Fig. 7, inside of which is larger than .

Comparing direct detection constraints to limits from gamma ray searches (Fig. 7 (b)), we find that for flavor-universal couplings and not too far from unity, direct searches are significantly more sensitive than continuum gamma ray searches in dwarf galaxies Bringmann et al. (2012) and competitive with the bounds from gamma ray line searches Garny et al. (2013). (Note that in Refs. Bringmann et al. (2012); Garny et al. (2013) these bounds are shown only for GeV, even though in principle, Fermi-LAT and H.E.S.S. are sensitive also to lower DM masses.) At GeV, direct detection limits are superseded by constraints from the anomalous magnetic moment of the electron and the muon (see Sec. V.1).

Looking into the future, Fig. 7 (c) illustrates that the sensitivity of direct detection experiments can be expected to improve by more than two orders of magnitude in the coming years thanks to the planned XENON1T Aprile (2012) and LUX-ZEPLIN (LZ) Malling et al. (2011) experiments. This will make direct DM searches highly sensitive to thermal relic DM. For XENON1T, we have assumed a total exposure of 2 200 kg yrs, while for LZ we use 10 000 kg yrs. In both cases, this corresponds to roughly 2 years of data taking. For comparison, we plot in Fig. 7 (c) also contours of constant (gray dot-dashed curves), where is the direct detection cross section averaged over the DM velocity distribution:

(17) |

Note that direct detection limits on are more than an order of magnitude weaker than direct detection limits on the cross section for DM–nucleon scattering through contact interactions. The reasons are the velocity dependence in as well as the fact that anapole interactions are proportional to the nuclear charge rather than the nuclear mass. As discussed in Sec. III.2, couplings to nuclear dipole moments are subdominant for the target material considered here.

For Dirac DM, Fig. 7 (d) shows that the qualitative picture is similar to Majorana DM, but the dependence on the lepton mass is less strong. Comparing the direct detection limits to constraints from the Fermi-LAT analysis of gamma ray signals from dwarf galaxies Ackermann et al. (2013), we find that for DM masses GeV, direct detection provides significantly stronger limits if and are not too different. In this case, also thermal production (horizontal blue band in Fig. 7 (d)) is excluded for GeV.

(a) | (b) |

(c) | (d) |

## Iv Collider searches for leptophilic dark matter

A set of constraints on leptophilic DM complementary to the limits from direct detection can be obtained from collider data. Since tree level production of DM at hadron colliders Cao et al. (2011); Beltran et al. (2010); Goodman et al. (2011); Bai et al. (2010); Goodman et al. (2010); Fox et al. (2012); Goodman and Shepherd (2011) is impossible in the leptophilic case, the strongest constraints are expected to come from mono-photon events at LEP Fox et al. (2011). In the future, mono-photon searches at a linear collider may improve on these bounds Dreiner et al. (2012).

Here, we apply the procedure described in Fox et al. (2011) to our toy model, Eq. (1). We simulate the process in CalcHEP 3.4 Belyaev et al. (2012) including the effect of initial state radiation and beamstrahlung (with default parameters) on the beam energy. We analyze the simulated events in a modified version of MadAnalysis 1.1.2 (from the MadGraph 4 package) Alwall et al. (2007) that implements the efficiencies and resolutions of the DELPHI detector at LEP Abdallah et al. (2005, 2009), see Fox et al. (2011) for details. We have checked that our simulation reproduces the predicted background from Abdallah et al. (2009) to very good accuracy. To set limits, we add our signal prediction to the background prediction from Abdallah et al. (2009), and compare to the DELPHI data from Fig. 1 of Abdallah et al. (2009), which corresponds to an integrated luminosity of 650 pb. Following Fox et al. (2011) we use a simple analysis to set limits on the Yukawa coupling as a function of and , and then convert these limits into constraints on , which are shown in Fig. 7 (b) and (d). Systematic uncertainties are subdominant compared to statistical uncertainties in DELPHI and are therefore neglected in our analysis.

We also estimate the sensitivity of a future linear collider with a center of mass energy GeV to leptophilic DM in our toy model. We simulate the signal and the dominant background in CalcHEP 3.4 Belyaev et al. (2012) while for the final state (with one photon escaping undetected) and for events (with an undetected pair) we follow Dreiner et al. (2012): we qualitatively include the background by simply increasing the background by 10%. For events, we reweight the spectrum by the energy dependent factor . Negative reweighting factors are excluded. The detector response of an ILC detector is modeled according to the information given in Abe et al. (2010); Bartels (2011); Dreiner et al. (2012). We assume an energy resolution of , where the notation means that the different terms correspond to separate, statistically independent Gaussian distributions. We restrict our analysis to the photon energy range (divided into 5 GeV bins) to remove events with on-shell production, and to the rapidity range . The detection efficiency is given by . We derive limits using a simple analysis, assuming an integrated luminosity of 50 fb and neglecting systematic uncertainties. Our projected ILC limits are included in Fig. 7 (c) and (d).

## V Constraints from Precision Experiments

### v.1 Lepton magnetic dipole moments

The extension of the SM by a DM particle and a charged mediator in our toy model Eq. (1) leads to a new contribution to the anomalous magnetic moment of leptons via the vertex correction loop shown in Fig. 8. This has been used previously in Bringmann et al. (2012); Fukushima and Kumar (2013) to constrain DM annihilation through charged mediators. In the case of complex Yukawa couplings, there can also be contributions to electric dipole moments, but we will not consider this possibility here. In the limit , , the anomalous magnetic moment of charged leptons is modified by Bringmann et al. (2012)

(18) |

For DM couplings to electrons, we compare Eq. (18) to the difference between the SM prediction for and the experimentally measured value, Aoyama et al. (2012) to derive the exclusion bound shown in Fig. 7 (b) for (lower edge of colored band) and for (upper edge of colored band).

For the of the muon, the difference between the measured best fit value and the theoretical prediction is Beringer et al. (2012). We add the experimental and theoretical uncertainties in quadrature. To account for the significant discrepancy between theory and experiment, we artificially inflate the error by linearly adding an ad-hoc uncertainty given by the central value of the discrepancy, . Note that the discrepancy has a sign opposite to the one predicted by Eq. (18). The resulting constraint on is shown in Fig. 7 (b).

We see that constraints are competitive with direct and indirect searches only at DM masses GeV.

### v.2 Positronium and muonium spectroscopy

Lepton–antilepton bound states such as positronium () and true muonium () are interesting laboratories for precision tests of QED because they can be studied accurately using spectroscopy, but are theoretically simpler than atoms. In particular, there are no nuclear effects that need to be taken into account. In our toy model for leptophilic DM, the box diagrams shown in Fig. 9 lead to an effective contact interaction of the form

(19) | |||

with | |||

(20) |

This contact interaction contributes to the electrostatic potential between the and , thus modifying the hyperfine splitting between the energy of the ortho-state (parallel spins, ) and the para-state (antiparallel spins, ). To obtain the new contribution to , we first calculate the new term in the Hamilton operator of the system by plugging explicit expressions for the and wave functions into (19), integrating over and adding a minus sign from the Legendre transform that converts the Lagrangian into the Hamiltonian as well as a factor 4 from the different ways in which the lepton fields can be contracted with the incoming and outgoing fermion states. The lepton wave functions are given by

(21) |

where is a non-relativistic particle or antiparticle Dirac spinor normalized to unity. We find that the energy of the ortho-state remains unchanged while the energy of the para-state is increased. The splitting between the two states is thus reduced, with

(22) |

For positronium, this implies

(23) |

which is an correction to Ishida et al. (2013), well below the experimental precision and the precision of the SM prediction. The reason for the low sensitivity is that positronium is a relatively large system, whereas the contact interaction is effective only at very short distance. The same is true for bound states.

More promising as a probe for contact interactions of the form of eq. (19), and of new physics in the lepton sector in general, seems to be “true muonium”, i.e. a bound state. Even though true muonium has never been directly produced and studied in the laboratory, precision experiments seem feasible Brodsky and Lebed (2009). For true muonium, we have

(24) |

which is only an correction to the leading term MHz Jentschura et al. (1997). Using Eq. (4) and comparing to Eq. (24), we obtain that to exclude thermal relic dark matter with GeV, , needs to be measured with an accuracy of 0.2 MHz.

### v.3 Lepton Flavor Violation

Even though in the simplest versions of our toy model motivated by supersymmetry, Eqs. (1) and (7), DM couplings to leptons are flavor diagonal, we now consider also the general Lagrangian Eq. (8) including flavor off-diagonal couplings. We derive constraints on these couplings from searches for the rare decays , and , which are mediated by the diagram shown in Fig. 8. Computing this diagram, we obtain for the decay rate

(25) |

where

(26) | ||||

(27) |

are Wilson coefficients in the effective Lagrangian

(28) |

and the loop functions , are given by

(29) | ||||

(30) |

We have used the definition , where are the masses of the charged mediators (see Eq. (8)).

We consider for illustrative purposes the special case where only three charged mediator , , exist, and where . This can be realized in supersymmetry if all left-handed sleptons are too heavy to be phenomenologically relevant. We obtain in this special case for the branching ratios (with the SM width )

(31) | ||||

(32) |

The expression for is identical to the one for
, with the replacements . With the current experimental limits Adam et al. (2013), Aubert et al. (2010) and Aubert et al. (2010), and using GeV, we then obtain the
following constraints on the elements of at :

Process | Coupling | Limit |
---|---|---|

7.0 | ||

6.1 |

We have seen in Eq. (4) that in our simplified model setup, at least one of the Yukawa couplings should be of order 0.1–1 to avoid DM overproduction. The above constraints show that flavor off-diagonal Yukawa couplings are therefore always subdominant. This justifies our neglecting them in the preceding sections.

We have also studied the decay , which constrains a different combination of Yukawa couplings because it also receives contributions from box diagrams similar to Fig. 9. If we assume that flavor-diagonal Yukawa couplings are , we obtain limits on the flavor off-diagonal couplings that are about a factor of 8 weaker than the limit from . To arrive at this estimate, we have used Ref. Kuno and Okada (2001) to express in terms of the Wilson coefficients of the effective operators in Eqs. (28) and (19). We have then compared the predicted branching ratio to the current experimental limit from Bellgardt et al. (1988); Beringer et al. (2012). Note that planned searches for will improve the limit on by up to four orders of magnitude Blondel et al. (2013).

## Vi Conclusions

In this paper, we have studied leptophilic dark matter models in which DM annihilation proceeds through a charged mediator and can therefore be accompanied by emission of a virtual internal bremsstrahlung photon. Such models are of great interest for indirect dark matter searches because internal bremsstrahlung can lead to spectral peaks in the gamma ray sky, a feature which is easily distinguishable from the large astrophysical gamma ray flux. Leptophilic DM models are also well motivated theoretically: they can be realized for instance in supersymmetric scenarios or radiative neutrino mass models, and in most cases, their parameter space is relatively unconstrained.

Here, we have established a connection between internal bremsstrahlung signals and loop-induced electromagnetic form factors of DM particles in leptophilic models. In particular, upon connecting the charged lepton lines in the internal bremsstrahlung diagrams in Fig. 1 to a loop, one immediately obtains the electromagnetic vertex corrections in Fig. 3. For Majorana DM, these lead to an anapole moment, while for Dirac DM, both anapole and magnetic dipole moments are generated, with the dipole moment being dominant in DM scattering processes. Interactions of the anapole and dipole moments with atomic nuclei then allow us to constrain the internal bremsstrahlung cross section using DM–nucleus scattering data from direct detection experiments. We have carried out this analysis for the most recent LUX and XENON100 data, and have found that direct detection constraints can be competitive with internal bremsstrahlung searches. This is true in particular if the mass splitting between the DM particle and the charged mediator is very small—the case which is also most interesting for internal bremsstrahlung searches due to the peaked gamma ray spectrum.

If DM is a Majorana fermion that couples universally to all charged leptons, direct detection limits are of the same order as limits from gamma ray line searches, and better than continuum gamma ray constraints from dwarf galaxies (see Fig. 7 (b)). Specifically, for small mass splitting , LUX constrains the internal bremsstrahlung cross section to be below cm/s at GeV. At DM masses of order 100 GeV, which have been invoked previously to explain a bump in Fermi-LAT gamma ray data Bringmann et al. (2012), LUX constraints imply that this interpretation is disfavored if DM couples to electrons or muons and if and differ by . If the last condition is significantly violated, however, the expected bump in the gamma ray spectrum becomes relatively broad, making line searches less sensitive. If , also direct searches for the charged mediator at colliders will impose important constraints, disfavoring GeV Bringmann et al. (2012); Liu et al. (2013). These constraints are ineffective if because the leptons from decay will be very soft in this case and thus hard to detect.

We note an interesting connection between our results and the scenario studied by Konishi et al. Konishi et al. (2013) to solve the cosmological lithium-7 problem in the Constrained Minimal Supersymmetric Standard Model (CMSSM) with sleptons that are nearly mass degenerate with the lightest neutralino. For the preferred mass range from Konishi et al. (2013), , this scenario would predict cm/s, well within the region testable by next generation direct detection experiments.

If DM is a Dirac fermion and the masses of and are of the same order of magnitude, but still differ by so that coannihilations are not yet relevant), direct detection constraints disfavor thermal relic production of DM for between 10–20 GeV and up to a few hundred GeV (see Fig. 7 (d)). For GeV, direct detection limits are also significantly stronger than astrophysical limits from gamma ray line searches and from continuum gamma rays searches in dwarf galaxies.

In the future, we expect the XENON1T and LUX-ZEPLIN experiments to improve these direct detection limits by about two orders of magnitude. These experiments will thus test the thermal relic hypothesis for DM masses of order , except for scenarios with a per cent level degeneracy between and , where coannihilations dominate in the early Universe. If a signal is detected, the spectrum of recoil events can be used to discriminate between anapole and dipole interaction and hence between Majorana and Dirac DM.

We have also studied constraints on our simplified model from low energy precision experiments. We confirm that bounds from the anomalous magnetic moment of the electron and the muon are weaker than the direct detection constraints at GeV. Searches for the lepton flavor violating decays , , and are very powerful in setting bounds on DM annihilation into flavor violating final states. Finally, we have studied the possibility of obtaining constraints from a future measurement of the hyperfine splitting in true muonium (a bound state). We have found such a measurement to be challenging for heavy DM ( GeV), where excluding thermal relic DM would require a measurement with a relative accuracy better than (see Eq. (24)). For lighter DM ( GeV), however, requirements are weaker and an interesting measurement may be possible.

In summary, our results show that direct dark matter searches are powerful tools to search for leptophilic DM even though DM–nucleus scattering occurs only at the loop level in this case. They are complementary to, and sometimes significantly superior to, indirect searches and precision experiments. Particularly in a scenario where a peak is observed in the cosmic gamma ray spectrum, but no other indirect hints for DM are found, virtual internal bremsstrahlung in a leptophilic DM model provides an attractive explanation. Our results show how this scenario can be confirmed in direct detection experiments by looking for the electromagnetic moment interactions of DM with nuclei. This illustrates once again that the search for Dark Matter is an interdisciplinary task, and that only a combination of different search strategies can yield optimal results.

## Acknowledgments

We would like to thank E. del Nobile, B. Kayser, T. Marrodán Undagoitia, H. Patel, T. Plehn, P. Schichtel, D. Schmeier, J. Tattersall, and C. Weniger for very helpful discussions. It is also a pleasure to thank Jonathan Schuster for his unusual but creative contributions to this work. JS acknowledges support from the IMPRS for Precision Tests of Fundamental Symmetries. JK would like to thank the Aspen Center for Physics (supported by NSF grant 1066293) for kind hospitality during part of this work.

## References

- Aharonian et al. (2012) F. Aharonian, D. Khangulyan, and D. Malyshev (2012), eprint 1207.0458.
- Bringmann et al. (2008) T. Bringmann, L. Bergstrom, and J. Edsjo, JHEP 0801, 049 (2008), eprint 0710.3169.
- Bell et al. (2011) N. F. Bell, J. B. Dent, T. D. Jacques, and T. J. Weiler, Phys.Rev. D83, 013001 (2011), eprint 1009.2584.
- Bringmann et al. (2012) T. Bringmann, X. Huang, A. Ibarra, S. Vogl, and C. Weniger (2012), eprint 1203.1312.
- Geringer-Sameth and Koushiappas (2011) A. Geringer-Sameth and S. M. Koushiappas, Phys.Rev.Lett. 107, 241303 (2011), eprint 1108.2914.
- Ackermann et al. (2013) M. Ackermann et al. (Fermi-LAT Collaboration) (2013), eprint 1310.0828.
- Hisano et al. (2011) J. Hisano, K. Ishiwata, and N. Nagata, Phys.Lett. B706, 208 (2011), eprint 1110.3719.
- Garny et al. (2012) M. Garny, A. Ibarra, M. Pato, and S. Vogl (2012), eprint 1207.1431.
- Garny et al. (2013) M. Garny, A. Ibarra, M. Pato, and S. Vogl (2013), eprint 1306.6342.
- Goodman et al. (2011) J. Goodman, M. Ibe, A. Rajaraman, W. Shepherd, T. M. Tait, et al., Phys.Lett. B695, 185 (2011), eprint 1005.1286.
- Bai et al. (2010) Y. Bai, P. J. Fox, and R. Harnik, JHEP 1012, 048 (2010), eprint 1005.3797.
- Fox et al. (2012) P. J. Fox, R. Harnik, J. Kopp, and Y. Tsai, Phys.Rev. D85, 056011 (2012), eprint 1109.4398.
- Lin et al. (2013) T. Lin, E. W. Kolb, and L.-T. Wang (2013), eprint 1303.6638.
- The CMS collaboration (2013) The CMS collaboration (2013), CMS-PAS-EXO-12-048.
- The ATLAS collaboration (2012) The ATLAS collaboration (2012), ATLAS-CONF-2012-147, ATLAS-COM-CONF-2012-190.
- Aad et al. (2013) G. Aad et al. (ATLAS Collaboration) (2013), eprint 1309.4017.
- Boyd (2013) J. Boyd (ATLAS), Overview of SUSY results from the ATLAS experiment (2013), talk given at the SUSY 2013 conference, slides available from http://susy2013.ictp.it.
- Richman (2013) J. D. Richman (CMS), Searches for Supersymmetry in the CMS Experiment (2013), talk given at the SUSY 2013 conference, slides available from http://susy2013.ictp.it.
- Adriani et al. (2009) O. Adriani et al. (PAMELA Collaboration), Nature 458, 607 (2009), eprint 0810.4995.
- Adriani et al. (2013) O. Adriani et al. (PAMELA Collaboration) (2013), eprint 1308.0133.
- Ackermann et al. (2012a) M. Ackermann et al. (Fermi LAT Collaboration), Phys.Rev.Lett. 108, 011103 (2012a), eprint 1109.0521.
- Aguilar et al. (2013) M. Aguilar et al. (AMS Collaboration), Phys.Rev.Lett. 110, 141102 (2013).
- Cirelli et al. (2009) M. Cirelli, M. Kadastik, M. Raidal, and A. Strumia, Nucl.Phys. B813, 1 (2009), eprint 0809.2409.
- Donato et al. (2009) F. Donato, D. Maurin, P. Brun, T. Delahaye, and P. Salati, Phys.Rev.Lett. 102, 071301 (2009), eprint 0810.5292.
- Nardi et al. (2009) E. Nardi, F. Sannino, and A. Strumia, JCAP 0901, 043 (2009), eprint 0811.4153.
- Bertone et al. (2009) G. Bertone, M. Cirelli, A. Strumia, and M. Taoso, JCAP 0903, 009 (2009), eprint 0811.3744.
- Fox and Poppitz (2009) P. J. Fox and E. Poppitz, Phys.Rev. D79, 083528 (2009), eprint 0811.0399.
- Evoli et al. (2011) C. Evoli, I. Cholis, D. Grasso, L. Maccione, and P. Ullio (2011), eprint 1108.0664.
- Kopp (2013) J. Kopp (2013), eprint 1304.1184.
- Cholis and Hooper (2013) I. Cholis and D. Hooper (2013), eprint 1304.1840.
- Bergstrom et al. (2013) L. Bergstrom, T. Bringmann, I. Cholis, D. Hooper, and C. Weniger, Phys.Rev.Lett. 111, 171101 (2013), eprint 1306.3983.
- Ibarra et al. (2013) A. Ibarra, A. S. Lamperstorfer, and J. Silk (2013), eprint 1309.2570.
- Adriani et al. (2010) O. Adriani et al. (PAMELA Collaboration), Phys.Rev.Lett. 105, 121101 (2010), eprint 1007.0821.
- Bartoli et al. (2012) B. Bartoli et al. (ARGO-YBJ Collaboration), Phys.Rev. D85, 022002 (2012), eprint 1201.3848.
- Hooper and Goodenough (2011) D. Hooper and L. Goodenough, Phys.Lett. B697, 412 (2011), eprint 1010.2752.
- Hooper and Linden (2011) D. Hooper and T. Linden, Phys.Rev. D84, 123005 (2011), eprint 1110.0006.
- Hooper (2012) D. Hooper (2012), eprint 1201.1303.
- Su et al. (2010) M. Su, T. R. Slatyer, and D. P. Finkbeiner, Astrophys.J. 724, 1044 (2010), eprint 1005.5480.
- Hooper and Slatyer (2013) D. Hooper and T. R. Slatyer (2013), eprint 1302.6589.
- Huang et al. (2013a) W.-C. Huang, A. Urbano, and W. Xue (2013a), eprint 1307.6862.
- Huang et al. (2013b) W.-C. Huang, A. Urbano, and W. Xue (2013b), eprint 1310.7609.
- Linden et al. (2011) T. Linden, D. Hooper, and F. Yusef-Zadeh, Astrophys.J. 741, 95 (2011), eprint 1106.5493.
- Kopp et al. (2009) J. Kopp, V. Niro, T. Schwetz, and J. Zupan, Phys. Rev. D80, 083502 (2009), eprint 0907.3159.
- Essig et al. (2012a) R. Essig, J. Mardon, and T. Volansky, Phys.Rev. D85, 076007 (2012a), eprint 1108.5383.
- Essig et al. (2012b) R. Essig, A. Manalaysay, J. Mardon, P. Sorensen, and T. Volansky (2012b), eprint 1206.2644.
- Schmidt et al. (2012) D. Schmidt, T. Schwetz, and T. Toma, Phys.Rev. D85, 073009 (2012), eprint 1201.0906.
- Ackermann et al. (2012b) M. Ackermann et al. (LAT Collaboration) (2012b), eprint 1205.2739.
- The Fermi-LAT collaboration (2013) The Fermi-LAT collaboration (2013), eprint 1305.5597.
- Boyarsky et al. (2012) A. Boyarsky, D. Malyshev, and O. Ruchayskiy (2012), eprint 1205.4700.
- Whiteson (2012) D. Whiteson (2012), eprint 1208.3677.
- Hektor et al. (2012) A. Hektor, M. Raidal, and E. Tempel (2012), eprint 1209.4548.
- Finkbeiner et al. (2012) D. P. Finkbeiner, M. Su, and C. Weniger (2012), eprint 1209.4562.
- Whiteson (2013) D. Whiteson (2013), eprint 1302.0427.
- Radescu (1985) E. Radescu, Phys.Rev. D32, 1266 (1985).
- Kayser and Goldhaber (1983) B. Kayser and A. S. Goldhaber, Phys.Rev. D28, 2341 (1983).
- Ho and Scherrer (2012) C. M. Ho and R. J. Scherrer (2012), eprint 1211.0503.
- Gresham and Zurek (2013) M. I. Gresham and K. M. Zurek (2013), eprint 1311.2082.
- Gao et al. (2013) Y. Gao, C. M. Ho, and R. J. Scherrer (2013), eprint 1311.5630.
- Del Nobile et al. (2014) E. Del Nobile, G. B. Gelmini, P. Gondolo, and J.-H. Huh (2014), eprint 1401.4508.
- Heo (2010) J. H. Heo, Phys.Lett. B693, 255 (2010), eprint 0901.3815.
- Masso et al. (2009) E. Masso, S. Mohanty, and S. Rao, Phys. Rev. D80, 036009 (2009), eprint 0906.1979.
- Del Nobile et al. (2012) E. Del Nobile, C. Kouvaris, P. Panci, F. Sannino, and J. Virkajarvi (2012), eprint 1203.6652.
- Weiner and Yavin (2012a) N. Weiner and I. Yavin (2012a), eprint 1209.1093.
- Barger et al. (2012a) V. Barger, W.-Y. Keung, D. Marfatia, and P.-Y. Tseng (2012a), eprint 1206.0640.
- Haisch and Kahlhoefer (2013) U. Haisch and F. Kahlhoefer (2013), eprint 1302.4454.
- Akerib et al. (2013) D. Akerib et al. (LUX Collaboration) (2013), eprint 1310.8214.
- Aprile et al. (2012) E. Aprile et al. (XENON100 Collaboration) (2012), eprint 1207.5988.
- Chun et al. (2010) E. J. Chun, J.-C. Park, and S. Scopel, JCAP 1002, 015 (2010), eprint arXiv:0911.5273.
- Buckley et al. (2013) M. R. Buckley, D. Hooper, and J. Kumar (2013), eprint 1307.3561.
- Liu et al. (2013) J. Liu, B. Shuve, N. Weiner, and I. Yavin (2013), eprint 1303.4404.
- Cao et al. (2009) Q.-H. Cao, E. Ma, and G. Shaughnessy, Phys.Lett. B673, 152 (2009), eprint 0901.1334.
- Bergstrom (1989) L. Bergstrom, Phys.Lett. B225, 372 (1989).
- Barger et al. (2009) V. Barger, Y. Gao, W. Y. Keung, and D. Marfatia, Phys.Rev. D80, 063537 (2009), eprint 0906.3009.
- Barger et al. (2012b) V. Barger, W.-Y. Keung, and D. Marfatia, Phys.Lett. B707, 385 (2012b), eprint 1111.4523.
- Toma (2013) T. Toma, Phys.Rev.Lett. 111, 091301 (2013), eprint 1307.6181.
- Giacchino et al. (2013) F. Giacchino, L. Lopez-Honorez, and M. H. Tytgat, JCAP 1310, 025 (2013), eprint 1307.6480.
- Ade et al. (2013) P. Ade et al. (Planck Collaboration) (2013), eprint 1303.5076.
- Hinshaw et al. (2013) G. Hinshaw et al. (WMAP), Astrophys.J.Suppl. 208, 19 (2013), eprint 1212.5226.
- Haber and Kane (1985) H. E. Haber and G. L. Kane, Phys.Rept. 117, 75 (1985).
- Konishi et al. (2013) Y. Konishi, S. Ohta, J. Sato, T. Shimomura, K. Sugai, et al. (2013), eprint 1309.2067.
- Fukushima and Kumar (2013) K. Fukushima and J. Kumar (2013), eprint 1307.7120.
- Frandsen et al. (2013) M. T. Frandsen, F. Sannino, I. M. Shoemaker, and O. Svendsen (2013), eprint 1312.3326.
- Frandsen et al. (2012) M. T. Frandsen, U. Haisch, F. Kahlhoefer, P. Mertsch, and K. Schmidt-Hoberg, JCAP 1210, 033 (2012), eprint 1207.3971.
- Sigurdson et al. (2004) K. Sigurdson, M. Doran, A. Kurylov, R. R. Caldwell, and M. Kamionkowski, Phys.Rev. D70, 083501 (2004), eprint astro-ph/0406355.
- Barger et al. (2011) V. Barger, W.-Y. Keung, and D. Marfatia, Phys.Lett. B696, 74 (2011), eprint 1007.4345.
- Fitzpatrick and Zurek (2010) A. Fitzpatrick and K. M. Zurek, Phys.Rev. D82, 075004 (2010), eprint 1007.5325.
- Banks et al. (2010) T. Banks, J.-F. Fortin, and S. Thomas (2010), eprint 1007.5515.
- Weiner and Yavin (2012b) N. Weiner and I. Yavin (2012b), eprint 1206.2910.
- Del Nobile et al. (2013) E. Del Nobile, G. Gelmini, P. Gondolo, and J.-H. Huh (2013), eprint 1306.5273.
- Kopp et al. (2010) J. Kopp, T. Schwetz, and J. Zupan, JCAP 1002, 014 (2010), eprint 0912.4264.
- Kopp et al. (2012) J. Kopp, T. Schwetz, and J. Zupan, JCAP 1203, 001 (2012), eprint 1110.2721.
- Aprile et al. (2011) E. Aprile et al. (XENON100 Collaboration), Phys.Rev.Lett. (2011), eprint 1104.2549.
- Pospelov and ter Veldhuis (2000) M. Pospelov and T. ter Veldhuis, Phys.Lett. B480, 181 (2000), eprint hep-ph/0003010.
- Chang et al. (2010) S. Chang, N. Weiner, and I. Yavin, Phys.Rev. D82, 125011 (2010), eprint 1007.4200.
- Jungman et al. (1996) G. Jungman, M. Kamionkowski, and K. Griest, Phys. Rept. 267, 195 (1996), eprint hep-ph/9506380.
- McCabe (2010) C. McCabe, Phys.Rev. D82, 023530 (2010), eprint 1005.0579.
- Farina et al. (2011) M. Farina, D. Pappadopulo, A. Strumia, and T. Volansky, JCAP 1111, 010 (2011), eprint 1107.0715.
- Green (2012) A. M. Green, Mod.Phys.Lett. A27, 1230004 (2012), eprint 1112.0524.
- Fairbairn et al. (2012) M. Fairbairn, T. Douce, and J. Swift (2012), eprint 1206.2693.
- Yellin (2002) S. Yellin, Phys. Rev. D66, 032005 (2002), eprint physics/0203002.
- Aprile (2012) E. Aprile (XENON1T collaboration) (2012), eprint 1206.6288.
- Malling et al. (2011) D. Malling, D. Akerib, H. Araujo, X. Bai, S. Bedikian, et al. (2011), see also http://www.hep.ucl.ac.uk/darkMatter/LZ.shtml, eprint 1110.0103.
- Cao et al. (2011) Q.-H. Cao, C.-R. Chen, C. S. Li, and H. Zhang, JHEP 1108, 018 (2011), eprint 0912.4511.
- Beltran et al. (2010) M. Beltran, D. Hooper, E. W. Kolb, Z. A. Krusberg, and T. M. Tait, JHEP 1009, 037 (2010), eprint 1002.4137.
- Goodman et al. (2010) J. Goodman, M. Ibe, A. Rajaraman, W. Shepherd, T. M. Tait, et al., Phys.Rev. D82, 116010 (2010), eprint 1008.1783.
- Goodman and Shepherd (2011) J. Goodman and W. Shepherd (2011), eprint 1111.2359.
- Fox et al. (2011) P. J. Fox, R. Harnik, J. Kopp, and Y. Tsai, Phys.Rev. D84, 014028 (2011), eprint 1103.0240.
- Dreiner et al. (2012) H. Dreiner, M. Huck, M. Kramer, D. Schmeier, and J. Tattersall (2012), eprint 1211.2254.
- Belyaev et al. (2012) A. Belyaev, N. D. Christensen, and A. Pukhov (2012), eprint 1207.6082.
- Alwall et al. (2007) J. Alwall, P. Demin, S. de Visscher, R. Frederix, M. Herquet, et al., JHEP 0709, 028 (2007), eprint 0706.2334.
- Abdallah et al. (2005) J. Abdallah et al. (DELPHI Collaboration), Eur.Phys.J. C38, 395 (2005), eprint hep-ex/0406019.
- Abdallah et al. (2009) J. Abdallah et al. (DELPHI Collaboration), Eur.Phys.J. C60, 17 (2009), eprint arXiv:0901.4486.
- Abe et al. (2010) T. Abe et al. (ILD Concept Group - Linear Collider Collaboration) (2010), eprint 1006.3396.
- Bartels (2011) C. Bartels (2011), URL http://www-library.desy.de/cgi-bin/showprep.pl?thesis11-034.
- Aoyama et al. (2012) T. Aoyama, M. Hayakawa, T. Kinoshita, and M. Nio, Phys.Rev.Lett. 109, 111807 (2012), eprint 1205.5368.
- Beringer et al. (2012) J. Beringer et al. (Particle Data Group), Phys.Rev. D86, 010001 (2012).
- Ishida et al. (2013) A. Ishida, T. Namba, S. Asai, T. Kobayashi, H. Saito, et al. (2013), eprint 1310.6923.
- Brodsky and Lebed (2009) S. J. Brodsky and R. F. Lebed, Phys.Rev.Lett. 102, 213401 (2009), eprint 0904.2225.
- Jentschura et al. (1997) U. Jentschura, G. Soff, V. Ivanov, and S. G. Karshenboim (1997), eprint hep-ph/9706401.
- Adam et al. (2013) J. Adam et al. (MEG Collaboration) (2013), eprint 1303.0754.
- Aubert et al. (2010) B. Aubert et al. (BaBar Collaboration), Phys.Rev.Lett. 104, 021802 (2010), eprint 0908.2381.
- Kuno and Okada (2001) Y. Kuno and Y. Okada, Rev.Mod.Phys. 73, 151 (2001), eprint hep-ph/9909265.
- Bellgardt et al. (1988) U. Bellgardt et al. (SINDRUM Collaboration), Nucl.Phys. B299, 1 (1988).
- Blondel et al. (2013) A. Blondel, A. Bravar, M. Pohl, S. Bachmann, N. Berger, et al. (2013), eprint 1301.6113.