# The effect of quark interactions on dark matter kinetic decoupling and the mass of the smallest dark halos

## Abstract

The kinetic decoupling of dark matter (DM) from the primordial plasma sets the size of the first and smallest dark matter halos. Studies of the DM kinetic decoupling have hitherto mostly neglected interactions between the DM and the quarks in the plasma. Here we illustrate their importance using two frameworks: a version of the Minimal Supersymmetric Standard Model (MSSM) and an effective field theory with effective DM-quark interaction operators. We connect particle physics and astrophysics obtaining bounds on the smallest dark matter halo size from collider data and from direct dark matter search experiments. In the MSSM framework, adding DM-quark interactions to DM-lepton interactions more than doubles the smallest dark matter halo mass in a wide range of the supersymmetric parameter space.

###### pacs:

95.35.+d## I Introduction

The nature of dark matter is still an open question, despite the growing evidence in support of its existence. Weakly interacting massive particles (WIMPs) are among the favorite candidates for dark matter. One of the unique features of WIMPs is that due to their weak interaction and heavy mass, they can lose thermal contact with the heat bath at a relatively early stage in the history of the Universe. The thermal, or as often called kinetic, decoupling of dark matter from the relativistic plasma sets the scale of the smallest dark matter halos, which are the first to form. Kinetic decoupling could provide a powerful cosmological probe on the properties of dark matter, in a way analogous to baryon decoupling, which has been unveiling the nature of our Universe through the baryon acoustic oscillations and the cosmic microwave background. Kinetic decoupling has duly received a fair amount of attention, being a compelling interdisciplinary avenue to explore from both particle physics and astrophysics perspectives (1); (3); (6); (4); (8); (13); (11); (12); (14); (16); (17); (15); (7); (2); (5).

Kinetic decoupling is a distinct process from chemical decoupling. Kinetic equilibrium between DM and the primordial plasma is maintained by rapid momentum exchange through scattering. Chemical equilibrium, on the other hand, holds when reactions that change the number of DM particles are active (for WIMPs, these reactions are typically annihilation processes). Kinetic equilibrium lasts much longer than chemical equilibrium due to the much higher number density of relativistic plasma particles available for scattering compared to the non-relativistic WIMPs necessary for annihilation. At chemical decoupling the annihilation/creation rates fall well below the Hubble expansion rate, and the dark matter particle number freezes out to a constant value per comoving volume. Still the dark matter remains in thermal equilibrium with the relativistic plasma through frequent elastic scattering. At kinetic decoupling, even such scattering processes become slower than the Hubble expansion rate, and the dark matter becomes free from the heat bath and begins to stream freely. The dark matter kinetic decoupling determines the low-mass cutoff scale for the size of the dark matter halos (protohalos), which is of importance for structure formation in the Universe.

In this paper, we pay particular attention to the role of DM-quark interactions in the process of kinetic decoupling. The role of these interactions has not been fully explored (to the extent of our knowledge, only Ref. (14) includes them in the numerical analysis), contrary to the numerous studies that include scattering of DM and leptons. In addition to the lack of due attention to DM-quark interactions, another strong motivation for looking into the quark interactions stems from the unprecedented wealth of data from the Large Hadron Collider (LHC) and direct dark matter search experiments, which provide us with a direct probe of the DM-quark interactions. These particle physics experiments would in principle put useful astrophysical bounds on the size of dark matter protohalos.

The main results presented in this paper are (i) a kinetic-decoupling Fokker-Planck equation that extends the expression in Bertschinger (11) to general WIMP models (it was first obtained in the PhD thesis work of J. Kasahara under the direction of one of the authors (PG) (12)) and (ii) a connection between the smallest mass of the dark matter protohalos, particle searches at hadron colliders, and dark matter direct detection experiments (we illustrate this connection in the context of effective DM-quark interaction operators and of the Minimal Supersymmetric Standard Model).

The layout of the paper is as follows. First, in Sec. II, we outline how we estimate the kinetic decoupling temperature; here we generalize the formalism presented in Ref. (11) and give a general expression for the momentum exchange momentum relaxation rate that can be applied to any non-relativistic WIMP. Then in Sec. III we apply this method to the DM-quark effective operators and show that the current collider and dark matter experiments set an upper bound on the smallest allowed protohalo mass. Finally in Sec. IV we quantify the importance of the DM-quark scattering relative to the DM-lepton scattering in the MSSM.

## Ii The dark matter kinetic decoupling and the smallest protohalos

In this section we outline the formalism to estimate the kinetic decoupling temperature and the corresponding smallest dark matter protohalo mass.

Let us first start with a heuristic order of magnitude argument for the kinetic decoupling process before jumping to the Boltzmann and Fokker-Planck equations. The momentum transfer per collision between the plasma (heat bath) at temperature and the heavy dark matter particles of mass is of order , much smaller than the average momentum of the dark matter particles, which is of order as follows from the fact that the average kinetic energy of the dark matter particles, , is of order . Many collisions, , are required for the dark matter to transfer a large part of its momentum to the plasma or to acquire it from the plasma. The momentum relaxation rate is thus , where is the elastic collision rate. The dark matter is in thermal equilibrium with the plasma when the momentum relaxation rate is larger than the Hubble expansion rate. Thermal decoupling occurs when the relaxation rate becomes of the order of the Hubble expansion rate, and this defines the kinetic decoupling temperature.

We can estimate the relaxation rate and consequently the decoupling temperature more accurately through the Fokker-Planck equation without approximating the dark matter as a perfect fluid or fully collisionless gas. Ref. (12) has re-derived the Fokker-Planck equation first discussed in Ref. (11) without assuming a specific cross section for WIMP-lepton scattering and allowing for massive particles in the plasma. The latter generalization is important at temperatures of the order of the electron or muon or quark masses. The former generalization goes beyond the zero-momentum transfer approximation of the previous literature (13); (14); (4); (5); (6) (Mandelstam variable ), and allows the study of general particle models. Indeed the only assumptions in Ref. (12) are that the dark matter particle is heavy ( or any other mass/energy scale) and that to momentum transfer is small. Under these assumptions, the Boltzmann equation for the time-dependence of the dark matter particle occupation number in a Friedmann-Robertson-Walker universe reduces to the Fokker-Planck equation

(1) |

with momentum relaxation rate

(2) |

Here the sum extends over the species in the relativistic plasma, with mass , occupation number , and statistical degrees of freedom (e.g. for the number of spin, color, and particle-antiparticle states of each flavor of Standard Model (SM) quark and antiquark). The + sign in and corresponds to bosons (stimulated emission), the - sign to fermions (Pauli blocking). And is the differential scattering cross section for the elastic scattering of and , written as a function of the Mandelstam variable and of the center-of-mass momentum , which in the heavy mass limit () equals the incoming momentum of particle in the plasma rest frame. Also, in the same limit,

(3) |

where is the invariant scattering amplitude and an overline indicates the usual sum over final polarizations and average over initial polarizations.

The Fokker-Planck equation (1) automatically conserves the number of particles per comoving volume, since the time derivative of is a total divergence in momentum space. The relaxation rate is an average over the thermal distribution of the mean square momentum transfer times the collision rate.

The expression for the momentum relaxation rate in Eq. 2 was first presented in Ref. (12) and it reproduces the formula for in Ref. (11) which only considered the bino scattering off massless leptons with a simplified invariant amplitude (the formula in Ref. (13); (14), which uses the forward scattering cross section in place of in Eq. (2), gives a value which is 20% larger). The formula given in Eq. 2 can be applied to a generic scattering amplitude and is thus of wide use. It has been implemented (12) in an extension of the DarkSUSY computer code for particle dark matter (18).

Multiplying the Fokker-Planck equation by the kinetic energy , and integrating in neglecting the stimulated emission or Pauli blocking factors (), leads to an equation for the kinetic temperature , defined as 2/3 of the average kinetic energy,

(4) |

where is the plasma temperature. Refs. (13); (11) present analytic solutions for the case proportional to a power of . At temperatures greater than the kinetic decoupling temperature , the particles are coupled to the plasma and . At temperatures smaller than , , as appropriate for non-relativistic particles of momenta that expand freely decoupled from the rest of the universe. In general a numerical solution of Eq. (4) is necessary. Here we content ourselves with estimating the kinetic decoupling temperature as the solution of (11); (12)

(5) |

The kinetic decoupling temperature is important to establish the mass of the smallest dark matter protohalos. As pointed out in Refs. (4); (7); (11), acoustic oscillations due to the coupling between the dark matter and the plasma damp the amplitude of fluctuations at scales smaller than the horizon size at decoupling. After decoupling, dark matter particles can stream freely without interacting with the plasma, and this process erases fluctuations up to the distance to which they can stream from the time of kinetic decoupling. The mass of the smallest protohalo is determined by the larger of the DM mass inside the horizon at kinetic decoupling and the DM mass within the free streaming length ,

(6) |

The DM mass contained within the free streaming length is

(7) |

where and is the dark matter density at the present time. The comoving free streaming scale

(8) |

where after kinetic decoupling, grows logarithmically during the radiation era () and saturates during the matter domination era (). We find (see Appendix)

(9) |

where

(10) |

is the total energy density at the time of matter-radiation equality. Finally, is a correction term that takes into account the change in the effective number of degrees of freedom between kinetic decoupling and the time of matter-radiation equality,

(11) |

where and are the energy and entropy degrees of freedom, respectively, and is a temperature much smaller than the electron-positron annihilation temperature and much larger than the temperature at equality, (we take keV for concreteness (see Appendix for the details)). For the dark matter velocity at decoupling , we use , with the coefficient obtained numerically (11).

Acoustic damping is characterized by the dark matter mass inside the horizon at decoupling

(12) | |||||

(13) |

where is the present temperature, is the present DM density, and the Hubble parameter at decoupling can be obtained from the Friedmann equation with total energy density at kinetic decoupling.

For the effective energy and entropy degrees of freedom and we adopted the model of Ref. (19) (equation of state B) as implemented in DarkSUSY (18). The factor , which was not present in Ref. (11), takes into account the change of comoving volume due to the entropy increase in the radiation. This ratio can be bigger than 10 for temperatures of the order of the QCD scale ( MeV in the model of our choice (19)).

The acoustic damping and free streaming scales are plotted in Fig. 1. For a wide range of parameters, the acoustic damping scale is larger than the free streaming scale and thus determines the cutoff scale of the smallest halo size. The free streaming length becomes more important than the acoustic damping scale when becomes large as seen in the figure. This behavior is expected because the acoustic damping length scales as , while the free streaming length roughly scales as . Notice that the scale of the smallest dark protohalo decreases with increasing kinetic decoupling temperatures.

## Iii Quark interactions in the effective field theory

It is illustrative to consider the effective field theory approach to see the essential features of the DM-quark interactions and their relations to the kinetic decoupling temperature. A great feature of studying such effective DM-quark operators is that one can study the range of the kinetic decoupling temperature allowed by the recent data from the LHC and dark matter direct detection experiments that directly probe DM-quark interactions.

The effective field theory approach has been studied extensively and we refer the readers to the existing literature for a complete survey of the effective operators (30); (25); (22); (28); (21); (27); (23); (24); (20); (26); (29) ^{1}

(14) | |||||

(15) |

These operators lead respectively to spin-independent and spin-dependent interaction whose interaction strength is set by the effective cutoff scale . The other DM-quark interaction operators besides the scalar and axial-vector operators vanish in the non-relativistic limit for Majorana fermion dark matter.

We make the simplifying assumptions that the DM couples universally to the Standard Model (SM) up and down type quarks via , and that it couples in through the quark mass suppression factor implied by chirality breaking.

The spin-independent cross section per proton reads

(16) |

where is the proton-DM reduced mass and is the fraction of the proton mass coupled to the scalar operator , which is equal to the mass fraction carried by quarks plus 6/27 of the mass fraction carried by gluons (for the default values of the nucleon parameters in DarkSUSY, which we use, ). The spin-dependent cross section per nucleon is

(17) |

where we use the value for the spin fraction carried by quarks in the nucleon (31); (32) and is the nuclear spin (1/2 for a free nucleon). For the dark matter direct search experiment constraints, we use the data from SIMPLE (33) for the spin-dependent cross section and from XENON100 (34) for the spin-independent interaction, which are currently the most sensitive direct detection experiments at the high end of the DM mass.

A complementary set of constraints comes from collider experiments. Collider constraints do not suffer from the astrophysical uncertainties, such as the local galactic dark matter density and velocity distribution, that afflict direct search experiments. The recent LHC data release from the CMS Collaboration (35) presented 1142 observed events in a mono jet analysis with leading jet transverse momentum GeV, pseudorapidity 2.4 and missing transverse energy GeV, to be compared with a Standard Model prediction of for the data sample of 4.7/fb total integrated luminosity at a center-of-mass energy of 7 TeV. We set a 2 collider lower bound on the effective coupling scale . We implement the effective operators by treating the scalar and axial-vector operators separately (one at a time) in Madgraph/Madevent. We use Pythia for the hadronization and the initial/final state radiation, treating the jets with the pycell subroutine.

The lower bounds on from the direct dark matter search and collider experiments are shown in Fig. 2 as a function of the dark matter mass , along with contours of the kinetic decoupling temperature . The LHC sensitivity decreases toward higher dark matter masses due to kinematic reasons. The collider bounds show the limits for the mono-jet events without additional nearby jets, in accord with our simplifying assumption of perfect efficiency. The direct detection experiments on the other hand suffer from the threshold to detect the DM recoil, and consequently, for the scalar operator, the collider constraints become stronger for light DM ( 10 GeV). The spin-dependent DM-nucleus scattering is not enhanced by the square of the atomic mass of the target nucleus, in contrast to the spin-independent interaction rate, and hence direct search experiments lead to bounds on the scale of the axial-vector operator that are weaker than the current LHC constraints .

The combination of direct detection and collider bounds in Fig. 2 forces the kinetic decoupling temperature to be larger than MeV, a regime in which DM-quark scattering becomes important. The power-law behavior of the kinetic decoupling temperature contours in the figure emerges from the relation and power-law dependence of the low-momentum transfer relaxation rate on , and .

Our effective-operator analysis provides us with upper bounds on the smallest dark protohalo mass directly from the current LHC and direct dark matter search limits, without having to scan the parameters space of a specific particle model. For instance, we find that for GeV the smallest allowed kinetic decoupling temperature is 350 MeV for the scalar operator and 150 MeV for the axial-vector operator, corresponding to upper limits on the smallest protohalo mass of and respectively. Note that here the protohalo cutoff mass scales as . Independently of the DM mass and spin-dependent or -independent interaction, there is an absolute lower bounds on the kinetic decoupling temperature of MeV and a corresponding absolute upper bound of (of order of the Earth’s mass) on the mass of the smallest dark protohalos.

In the model-independent analysis of this section, each quark was decoupled at its mass scale but the QCD phase transition was not taken account of which is currently heavily model dependent. Indeed our simplified analysis adding quarks to the plasma as a free gas is not really accurate during the QCD phase transition, when hadrons are also present, whose treatment would require more input from lattice simulations. Nonetheless, our findings show that the current LHC and direct dark matter search data impose lower bounds on that are well in the QCD phase transition regime, supporting our statement that DM-quark interactions, hitherto neglected in studies of kinetic decoupling, are important and must be included. The inclusion of DM-quark interactions will become even more important as forthcoming LHC/dark matter experiments probe DM-quark interactions further, potentially pushing the lower bound on up and above the QCD transition regime.

## Iv Quark interactions in a concrete model

An advantage of specifying a concrete model rather than effective DM-quark operators is that in a concrete model all the DM interaction terms are present, including the interference terms, which are cumbersome to specify in a model independent approach treating the different effective operators separately. In this section, we use a specific model to examine the relative significance of DM-quark interactions in comparison to DM interactions with other particles.

As a concrete model, we choose the Minimal Supersymmetric Standard Model (MSSM), in which at the varying of the model parameters the values of are known to span a wide range, from tens of MeV to more than a GeV (1); (3); (6); (4); (8); (13); (11); (14); (16); (17); (15); (7); (2); (5). Because of this, the MSSM provides a good theoretical proving ground to check the relative importance of the DM-quark interactions compared with the DM scattering off leptons that has been discussed extensively in the previous literature.

The elastic scattering of neutralinos off the fermions in the MSSM occurs through the exchange of gauge bosons, Higgs bosons, and sfermions. We extended the modified DarkSUSY first developed in Ref. (12) so as to include all possible dark matter interactions and interference terms in the MSSM for dark matter scattering off fermions in the primordial plasma. For the effective number of degrees of freedom, we adopt the model of Ref. (19) with QCD phase transition temperature of 154 MeV (equation of state B (19), the default model in DarkSUSY (18)). For the scattering of DM off quarks, the Fokker-Planck equation is numerically fully solved down to 154 MeV including all the relevant interactions and we simply turn off the DM-quark scatterings below 154 MeV. Even though the temperature at which the asymptotically free quark description becomes valid is expected to be higher than the QCD phase transition temperature, this simplification suffices for our purpose of showing the potential significance of DM-quark scattering in estimating the dark matter protohalo mass. Choosing a higher temperature for the threshold of DM-quark scattering to ensure the validity of the free quark description does not change the conclusions of this section (adding DM-quark scattering besides DM-lepton scattering can increase the protohalo mass estimation by a factor of 2 or more).

We scanned the MSSM-7 parameter space, characterized (36) by seven parameters specified at the electroweak scale: 2 trilinear A-terms, a soft sfermion mass parameter, a gaugino mass parameter, and three Higgs-sector parameters. The trilinear A-terms and the soft sfermion mass matrices were assumed to be diagonal to avoid flavor-changing-neutral-current issues and they were parameterized as , . All the soft sfermion mass matrices are times the identity matrix. The gaugino mass GUT relation was assumed and the SU(2) gaugino mass was chosen as the free parameter. The Higgs sector was parameterized by the CP-odd Higgs boson mass , the Higgsino mass parameter and the ratio of Higgs vacuum expectation values. These seven parameters were randomly scanned over the ranges , TeV and for the remaining mass parameters , , , and . We applied the phenomenological bounds available in the DarkSUSY subroutines (18) for the theoretical consistency of the model and the experimental constraints. We excluded parameter sets that have charged or colored vacua, that have potentials unbounded from below, that do not have the neutralino as the lightest supersymmetric particle, and that violate experimental bounds on supersymmetric masses or rare processes such as .

Because our purpose here is to show quantitatively the significance of the DM-quark scattering we do not require that the neutralino thermal relic abundance should match the observed cold dark matter value (actually, through non-thermal production, entropy production, or non-standard expansion history, the neutralino thermal relic abundance can be smaller or larger than the cosmological value , even if the neutralinos comprise all of the cold dark matter (37)).

In Fig. 3 we compare the neutralino kinetic decoupling temperatures with and without including DM-quark interactions for a random scan in the MSSM-7 parameter space. The open circles represent the values computed including the DM-quark interactions, and the crosses are the values computed without including the DM-quark interactions.
We first observe that kinetic decoupling can occur above the QCD phase transition scale, as first pointed out in Ref. (8), which performed an analogous numerical scan without including the DM-quark interactions.
For , the quarks are bounded inside the hadrons and the pions were much less abundant than the light leptons in the thermal bath. We numerically checked that the inclusion of pion scattering affects the kinetic decoupling temperature by less than one per cent in our scenarios ^{2}

This conclusion on the importance of including DM-quark interactions would not change even if we had switched off DM-quark scattering up to temperatures four or five times higher than the QCD phase transition temperature, as to guarantee the presence of free quarks, because our MSSM-7 parameter scan already covers such regime of well above the QCD phase transition temperature ^{3}

## V Discussion/Conclusion

Studies of particle dark matter can provide us with a unique link between particle physics and astrophysics. To illustrate one such connection, we examined the importance of DM-quark interactions in the kinetic decoupling of particle dark matter and the consequent mass of the smallest dark protohalos. We focus on a model independent analysis through effective DM-quark interaction operators and a model specific analysis within a seven-parameter MSSM.

In the effective operator approach, in which we assumed DM-lepton interactions are negligible, we found that current direct dark matter search and collider constraints force the DM kinetic decoupling temperature to exceed MeV, whether the interaction is spin-dependent or spin-independent. This sets an upper limit of (of order of the Earth’s mass) to the mass of the smallest dark protohalo.

In the MSSM-7 study, we found that inclusion of neutralino-quark interactions can increase the smallest dark protohalo mass by more than a factor of 2 whenever the neutralino decoupling temperature exceeds the QCD transition temperature and acoustic damping dominates over free streaming.

The constraints on the kinetic decoupling temperature from the current LHC and direct dark matter search experiments turned out to be in the regime of the QCD phase transition. The forthcoming LHC and direct dark matter search data will most likely push this bound up, possibly beyond the quark-hadron transition, into a regime in which our approximation of a gas of free quarks is definitely applicable.

If DM-quark interactions are discovered at the LHC or in direct detection experiments, measurements of their strength can provide us with a lower bound of the mass of the smallest dark matter protohalos, once proper account is taken of possible additional interactions with leptons.

The survival of the smallest dark matter protohalos to the present time, and their observability, depend strongly on complicated tidal forces and other astrophysics such as stellar interactions, and have been the subject of vigorous debate. More input from computer simulations of structure formation in the universe would help clarify these issues, despite the vast dynamic range of masses and non-linear effects involved in simulating the relevant processes (42); (47); (45); (46); (43); (44); (48).

Although we discussed only the hadron collider (LHC) data that constrain the DM-quark interaction operators, an analogous exercise can be performed for the DM-lepton interactions considering a lepton collider (38); (39) which is crucial for a concrete model such as the MSSM where DM-lepton interactions cannot be neglected for the estimation of the kinetic decoupling. Inclusion of non-collider experiments and indirect search experiments also deserve further studies. For instance the constraints from the cosmic antiproton flux can put relatively tight constraints on the axial-vector operator (40); (41), and more systematic studies of the dark matter kinetic decoupling including additional experiments and additional effective operators beyond those considered here are worthy of future work.

Our results on the impact of collider and dark matter search experiments onto the formation of structure in the early universe calls for further exploration of the creation and evolution of dark matter protohalos, in view of their potential role in probing the nature of the dark matter.

## Acknowledgments

For the numerical estimation of kinetic decoupling temperature, we have modified the DarkSUSY where the momentum relaxation rate and the full complement of supersymmetric scattering cross sections were first implemented by Junya Kasahara for his PhD thesis project (12). This work was supported in part by a Grant-in-Aid for the Global COE program and for Scientific Research from the MEXT of Japan (JH and KK), and by NSF award PHY-1068111 (P.G.). We thank the Kavli Institute for Theoretical Physics China where part of this work was conducted. P.G. also thanks the Korean Institute for Advanced Studies, Nagoya University and Seoul National University for support and hospitality during the completion of this work.

## Appendix A Derivation of the free streaming length formula

Here we sketch the derivation of Eq. (9) for the free streaming length, which includes the change in the effective degrees of freedom before matter-radiation equality. Since the expressions for the energy and entropy degrees of freedom and usually do not extend down to temperatures , we divide the integration in

(19) |

into two parts, which we join at a temperature defined so that for . The temperature is thus much smaller than the electron-positron annihilation temperature but much larger than the temperature at equality.

For , we use an FRW model with matter and radiation,

(20) |

where is the Hubble parameter at the matter-radiation equality (), with equal to the total density at that time (the contribution to from the low-redshift cosmological constant term is negligible). We find

(21) | ||||

(22) | ||||

(23) | ||||

(24) |

where, on using for ,

(25) |

with

(26) |

For , we include the full temperature dependence of the degrees of freedom. Using conservation of total entropy , from which

(27) |

we find

(28) | ||||

(29) | ||||

(30) |

In the last step, we wrote in terms of , and then , using the matter and radiation model at where ,

(31) |

After adding Eqs. (24) and (30), we formally extend the integration in Eq. (24) from to using the definition

(32) |

and the mathematical identity

(33) |

Notice that since for , the variable in Eq. (32) is not the conformal time at plasma temperatures , as instead it is at , but is only a convenient mathematical variable. This gives

(34) |

with

(35) |

The approximation then gives the expression of in Eq. (11).

### Footnotes

- We do not consider the UV completion of the effective operators of our interest. Strictly speaking the effective theory breaks down if the 4-momentum transfer is comparable to or larger than the mass of a particle mediating the interaction (the typical momentum transfer is of order for DM kinetic decoupling). See for instance Refs. (23); (30) for further discussions on the validity of the effective operator approach in the case of light mediators.
- The leading order DM-pion coupling term can be estimated with a coupling constant (9) (18) where , with up or down quark in this case, is the coefficient of the effective scalar current operator in the MSSM (10). This form of was derived using the soft-pion technique in the framework of chiral perturbation theory, which is valid for an energy scale ( is the pion decay constant), and hence is applicable to the soft-pion case.
- That such high are common can be seen for instance in the MSSM scans in Ref. (14), which used a different QCD phase transition model and a different method of estimating the kinetic decoupling temperature.

### References

- C. Schmid, D. J. Schwarz and P. Widerin, Phys. Rev. D 59, 043517 (1999) [astro-ph/9807257].
- C. Boehm, P. Fayet and R. Schaeffer, Phys. Lett. B 518, 8 (2001) [astro-ph/0012504].
- X. -L. Chen, M. Kamionkowski and X. M. Zhang, Phys. Rev. D 64, 021302 (2001) [astro-ph/0103452].
- S. Hofmann, D. J. Schwarz and H. Stoecker, Phys. Rev. D 64, 083507 (2001) [astro-ph/0104173].
- V. Berezinsky, V. Dokuchaev and Y. Eroshenko, Phys. Rev. D 68, 103003 (2003) [astro-ph/0301551].
- A. M. Green, S. Hofmann and D. J. Schwarz, JCAP 0508, 003 (2005) [astro-ph/0503387].
- A. Loeb and M. Zaldarriaga, Phys. Rev. D 71, 103520 (2005) [astro-ph/0504112].
- S. Profumo, K. Sigurdson and M. Kamionkowski, Phys. Rev. Lett. 97, 031301 (2006) [astro-ph/0603373].
- G. Prezeau, A. Kurylov, M. Kamionkowski and P. Vogel, Phys. Rev. Lett. 91, 231301 (2003) [astro-ph/0309115].
- M. Drees and M. Nojiri, Phys. Rev. D 48, 3483 (1993) [hep-ph/9307208].
- E. Bertschinger, Phys. Rev. D 74, 063509 (2006) [astro-ph/0607319].
- J. Kasahara, Ph.D. dissertation, University of Utah (2009; ISBN 9781109295320; UMI microform 3368246); J. Kasahara and P. Gondolo, unpublished.
- T. Bringmann and S. Hofmann, JCAP 0407, 016 (2007) [hep-ph/0612238].
- T. Bringmann, New J. Phys. 11, 105027 (2009) [arXiv:0903.0189 [astro-ph.CO]].
- X. -J. Bi, P. -F. Yin and Q. Yuan, Phys. Rev. D 85, 043526 (2012) [arXiv:1106.6027 [hep-ph]].
- L. G. van den Aarssen, T. Bringmann and Y. CGoedecke, arXiv:1202.5456 [hep-ph].
- J. M. Cornell and S. Profumo, arXiv:1203.1100 [hep-ph].
- P. Gondolo, J. Edsjo, P. Ullio, L. Bergstrom, M. Schelke and E. A. Baltz, JCAP 0407, 008 (2004) [astro-ph/0406204].
- M. Hindmarsh and O. Philipsen, Phys. Rev. D 71, 087302 (2005) [hep-ph/0501232].
- M. Beltran, D. Hooper, E. W. Kolb, Z. A. C. Krusberg and T. M. P. Tait, JHEP 1009, 037 (2010) [arXiv:1002.4137 [hep-ph]].
- P. Agrawal, Z. Chacko, C. Kilic and R. K. Mishra, arXiv:1003.1912 [hep-ph].
- J. Goodman, M. Ibe, A. Rajaraman, W. Shepherd, T. M. P. Tait and H. -B. Yu, Phys. Lett. B 695, 185 (2011) [arXiv:1005.1286 [hep-ph]].
- Y. Bai, P. J. Fox and R. Harnik, JHEP 1012, 048 (2010) [arXiv:1005.3797 [hep-ph]].
- J. Goodman, M. Ibe, A. Rajaraman, W. Shepherd, T. M. P. Tait and H. -B. Yu, Phys. Rev. D 82, 116010 (2010) [arXiv:1008.1783 [hep-ph]].
- J. Goodman, M. Ibe, A. Rajaraman, W. Shepherd, T. M. P. Tait and H. -B. Yu, Nucl. Phys. B 844, 55 (2011) [arXiv:1009.0008 [hep-ph]].
- A. Rajaraman, W. Shepherd, T. M. P. Tait and A. M. Wijangco, Phys. Rev. D 84, 095013 (2011) [arXiv:1108.1196 [hep-ph]].
- P. J. Fox, R. Harnik, J. Kopp and Y. Tsai, arXiv:1109.4398 [hep-ph].
- K. Cheung, P. -Y. Tseng, Y. -L. S. Tsai and T. -C. Yuan, arXiv:1201.3402 [hep-ph].
- A. L. Fitzpatrick, W. Haxton, E. Katz, N. Lubbers and Y. Xu, arXiv:1203.3542 [hep-ph].
- J. March-Russell, J. Unwin and S. M. West, arXiv:1203.4854 [hep-ph].
- J. R. Ellis, K. A. Olive and C. Savage, Phys. Rev. D 77, 065026 (2008) [arXiv:0801.3656 [hep-ph]].
- M. Alekseev et al. [COMPASS Collaboration], Phys. Lett. B 660, 458 (2008) [arXiv:0707.4077 [hep-ex]].
- M. Felizardo, T. Girard, T. Morlat, A. C. Fernandes, F. Giuliani, A. R. Ramos, J. G. Marques and M. Auguste et al., arXiv:1106.3014 [astro-ph.CO].
- E. Aprile et al. [XENON100 Collaboration], Phys. Rev. Lett. 107, 131302 (2011) [arXiv:1104.2549 [astro-ph.CO]].
- The CMS collaboration, EXO11059, 4.7/fb (2012)
- L. Bergstrom and P. Gondolo, Astropart. Phys. 5, 263 (1996) [hep-ph/9510252].
- G. B. Gelmini and P. Gondolo, Phys. Rev. D 74, 023510 (2006) [hep-ph/0602230].
- P. J. Fox, R. Harnik, J. Kopp and Y. Tsai, Phys. Rev. D 84, 014028 (2011) [arXiv:1103.0240 [hep-ph]].
- Y. Mambrini and B. Zaldivar, JCAP 1110, 023 (2011) [arXiv:1106.4819 [hep-ph]].
- O. Adriani et al. [PAMELA Collaboration], Phys. Rev. Lett. 105, 121101 (2010) [arXiv:1007.0821 [astro-ph.HE]].
- K. Cheung, P. -Y. Tseng and T. -C. Yuan, JCAP 1101, 004 (2011) [arXiv:1011.2310 [hep-ph]].
- J. Diemand, B. Moore and J. Stadel, Nature 433, 389 (2005) [astro-ph/0501589].
- V. Berezinsky, V. Dokuchaev and Y. Eroshenko, Phys. Rev. D 73, 063504 (2006) [astro-ph/0511494].
- A. M. Green and S. P. Goodwin, Mon. Not. Roy. Astron. Soc. 375, 1111 (2007) [astro-ph/0604142].
- T. Goerdt, O. Y. Gnedin, B. Moore, J. Diemand and J. Stadel, Mon. Not. Roy. Astron. Soc. 375, 191 (2007) [astro-ph/0608495].
- S. M. Koushiappas, New J. Phys. 11, 105012 (2009) [arXiv:0905.1998 [astro-ph.CO]].
- A. Schneider, L. Krauss and B. Moore, Phys. Rev. D 82 (2010) 063525 [arXiv:1004.5432 [astro-ph.GA]].
- J. P. Ostriker and P. J. Steinhardt, Science 300, 1909 (2003) [astro-ph/0306402].