Issue 
A&A
Volume 535, November 2011



Article Number  A26  
Number of page(s)  10  
Section  Cosmology (including clusters of galaxies)  
DOI  https://doi.org/10.1051/00046361/201116914  
Published online  27 October 2011 
WMAP7 and future CMB constraints on annihilating dark matter: implications for GeVscale WIMPs
^{1} Tartu Observatory, Tõravere 61602, Estonia
email: gert@aai.ee
^{2}
MaxPlanckInstitut für Astrophysik, KarlSchwarzschildStr. 1, 85741 Garching, Germany
^{3}
Canadian Institute for Theoretical Astrophysics, 60 St. George Street, Toronto, ON M5S 3H8, Canada
email: jchluba@cita.utoronto.ca
^{4} National Institute of Chemical Physics and Biophysics, Tallinn 10143, Estonia
email: [andi.hektor;martti.raidal]@cern.ch
Received: 17 March 2011
Accepted: 30 August 2011
Aims. We calculate constraints from current and future cosmic microwave background (CMB) measurements on annihilating dark matter (DM) with masses below the electroweak scale: m_{DM} = 5 − 100 GeV. In particular, we assume the Swave annihilation mode to be dominant, and focus our attention on the lower end of this mass range, as DM particles with masses m_{DM} ~ 10 GeV have recently been claimed to be consistent with the CoGeNT and DAMA/LIBRA results, while also providing viable DM candidates to explain the measurements of Fermi and WMAP haze. We study the model (in)dependence of the CMB power spectra on particle physics DM models, largescale structure formation and cosmological uncertainties. We attempt to find a simple and practical recipe for estimating current and future CMB bounds on a broad class of DM annihilation models.
Methods. We use a modelindependent description for DM annihilation into a wide set of Standard Model particles simulated by PYTHIA Monte Carlo. Our Markov chain Monte Carlo calculations used for finding model constraints involve realistic CMB likelihoods and assume a standard 6parameter ΛCDM background cosmological model, which is extended by two additional DM annihilation parameters: m_{DM} and ⟨ σ_{A}υ ⟩ /m_{DM}.
Results. We show that in the studied DM mass range the CMB signal of DM annihilations is independent of the details of largescale structure formation, distribution, and profile of DM halos and other cosmological uncertainties. All particle physics models of DM annihilation can be described with only one parameter, the fraction of energy carried away by neutrinos in DM annihilation. As the main result we provide a simple and rather generic fitting formula for calculating CMB constraints on the annihilation cross section of light WIMPs. We show that thermal relic DM in the CoGeNT, DAMA/LIBRA favored mass range is in a serious conflict with present CMB data for the annihilation channels with few neutrinos, and will definitely be tested by the Planck mission for all possible DM annihilation channels. Also, our findings strongly disfavor the claim that thermal relic DM annihilations with m_{DM} ~ 10 GeV and ⟨ σ_{A}υ ⟩ ~ 9 × 10^{25} cm^{3} s^{1} could be a cause of Fermi and WMAP haze.
Key words: cosmic background radiation / cosmological parameters / dark matter / cosmology: theory / elementary particles / diffuse radiation
© ESO, 2011
1. Introduction
According to our current understanding, the energymatter content of the Universe is dominated by dark components: dark energy (DE) (~73%) and dark matter (DM) (~22%), with ordinary baryonic matter accounting for only ~5% of the total (~critical) density (e.g., see Komatsu et al. 2011). In contrast to the poorly understood DE – the substance causing the Universe to expand in an accelerated fashion – we have physically wellmotivated models for the DM. Among those, the most promising scenario states that the DM of the Universe consists of thermal relic density of stable weakly interacting massive particles (WIMPs). It is quite miraculous that having particles with masses and annihilation cross sections set by the electroweak scale automatically provide the right DM density after freezeout (Jungman et al. 1996; Bertone et al. 2005).
The WIMP hypothesis, along with its potentially observable phenomenology, has initiated strong effort in the particle and astrophysics communities to try to find other than purely gravitational manifestations of DM. So far, we have good knowledge of DM only through its gravitational effects, starting from the scale of galaxies and galaxy clusters, up to the cosmologically largest observable scales (Jungman et al. 1996; Bertone et al. 2005; Einasto 2009). However, as there is already an impressive list of ongoing and upcoming direct DM detection experiments along with various indirect means of detection (see Feng 2010, for overview), the hopes are very high that in the nearest future the mystery of DM might at last be solved. Indeed, the first signals from DM particles could potentially have already been detected: the (expected) annual modulation signal from DAMA/LIBRA (Bernabei et al. 2010), signals from the CDMS (Ahmed et al. 2011) and CoGeNT (Aalseth et al. 2011) nuclear recoil experiments, anomalies of the cosmic ray positrons as revealed by PAMELA satellite (Adriani et al. 2009), or positrons+electrons as obtained by the Fermi satellite (Abdo et al. 2009) and HESS atmospheric Cherenkov telescope (Aharonian et al. 2009). While the cosmic ray positron anomaly can possibly be explained by TeVscale DM (Bergström et al. 2008; Barger et al. 2009; Cirelli et al. 2009b; ArkaniHamed et al. 2009; Fox & Poppitz 2009), the signal from CoGeNT calls for light WIMPs within the mass range 7 − 12 GeV (Aalseth et al. 2011). Consistent analyses of combined data from CoGeNT and DAMA/LIBRA determine the light DM mass to be 6 − 8 GeV (Hooper et al. 2010). Although this mass range is probed by the CDMS (Ahmed et al. 2011), XENON10 (Angle et al. 2009), and XENON100 (Aprile et al. 2010; XENON100 Collaboration et al. 2011) experiments, interpretation of those results (Kopp et al. 2010; Schwetz 2010) requires an ability to reliably reconstruct nuclear recoils at very low energy (Savage et al. 2011), as well as precise knowledge of DM distribution and velocity in the local halo. Therefore the CoGeNT and DAMA/LIBRA hints of light DM cannot be ruled out unambiguously.
In addition, there is an independent positive claim of the existence of O(10) GeV mass DM. A recent study by Dobler et al. (2011) also suggests that annihilating DM with similarly low masses (m_{DM} = 1 − 20 GeV) may give a good match to the observed Fermi and WMAP haze (Dobler et al. 2010; Dobler & Finkbeiner 2008). However, all those claims depend strongly on the knowledge of the profile of the DM halo of our Galaxy and precise knowledge of local DM density and halo substructure.
Thus the several interesting claims of the existence of O(10) GeV mass DM call for modelindependent tests of the light DM scenario. Since lower DM particle masses imply higher number densities (n_{DM} ∝ Ω_{DM}h^{2}/m_{DM}), and as the energy input from annihilations scales as $\mathrm{\propto}{\mathit{n}}_{\mathrm{DM}}^{\mathrm{2}}$, one might expect strong constraints on annihilation cross section, which might possibly reach below the standard thermal production value of ⟨ σ_{A}υ ⟩ ≃ 3 × 10^{26} cm^{3} s^{1} (Jungman et al. 1996). The constraints from gammaray measurements (Bell & Jacques 2009; Bertone et al. 2009; Bergström et al. 2009; Cirelli & Panci 2009; Meade et al. 2010; Cirelli et al. 2010; Papucci & Strumia 2010; Hütsi et al. 2010; Baxter et al. 2010; Arina & Tytgat 2011; Vincent et al. 2010; Zavala et al. 2011) along with CMB bounds (Galli et al. 2009; Slatyer et al. 2009; Cirelli et al. 2009a; Hütsi et al. 2009; Kanzaki et al. 2010), indicate that this might indeed be the case.
In this paper we investigate how well DM annihilation cross sections for WIMPs with masses m_{DM} = 5 − 100 GeV can be constrained with CMB measurements, in particular focusing on the lower end of this mass range. The main advantage of CMB over other indirect probes, like gamma ray and cosmic ray measurements, is that it is practically insensitive to the complications caused by the nonlinear evolution of the cosmic density field. As our analysis shows, for any realistic structure formation scenario the CMB bounds on annihilating DM arise solely around the redshifts of z ~ 1000, while the contribution from lower redshift cosmic structures is completely negligible.
CMB constraints on annihilating DM have been obtained in several earlier studies: e.g., Padmanabhan & Finkbeiner (2005); Mapelli et al. (2006); Zhang et al. (2006); Galli et al. (2009); Slatyer et al. (2009); Cirelli et al. (2009a); Hütsi et al. (2009); Kanzaki et al. (2010). Even though most of these analyses have assumed a simple “on the spot” approximation for the energy deposition^{1}, more recent studies (Slatyer et al. 2009; Hütsi et al. 2009) followed the energy transport problem including various energyloss mechanisms in a more realistic way. Compared to the analysis of Slatyer et al. (2009), which partially relies on the previously derived “on the spot” results of Galli et al. (2009), in this paper we perform a more elaborate treatment for the energy deposition joined to the analysis of the CMB data via Markov chain Monte Carlo calculations that incorporate the most recent WMAP likelihood code. Also, we make an attempt to unify the results from various annihilation channels and provide a simple and rather generic fitting formula for calculating CMB constraints on annihilation cross section ⟨ σ_{A}υ ⟩ for a broad range of annihilating DM models.
Our paper is organized as follows. In Sect. 2 we give a brief description of the energy input from DM annihilation and provide a simple treatment for its propagation. The effect on CMB temperature and polarization fluctuations is investigated in Sect. 3. Section 4 presents our main results about current and future CMB constraints. Our summary is given in Sect. 5.
2. Input signals and their propagation
As in Cirelli et al. (2009b, 2011), we treat our input signals from DM annihilation in an as model independent a way as possible. In Cirelli et al. (2011) the twoparticle annihilation channels to all Standard Model (SM) particles were considered: leptons, quarks, photons, gluons, weakinteraction gauge bosons, Higgs boson, and neutrinos. In addition, annihilations to four leptons via an intermediate new boson V were considered. Such a treatment can be considered as model independent since realistic models can always be decomposed into these basic channels where the particular branching ratios between the channels are given by the underlying theoretical particle physics model.
Since in our work we focus on DM particle masses below 100 GeV, out of all of the above channels the following remain: DM DM → SM $\overline{)\mathrm{SM}}$, where SM = { e,μ,τ,q,c,b,γ,g}; plus 4lepton channels via V. Here q denotes the light quarks u, d, and s. Because the masses of interest in this work are mostly below the masses of the electroweak gauge bosons, W^{ ± } and Z, those channels are left out. For the same reason we also do not need to distinguish between left and righthanded particles. Also, the neutrino channels in this case provide only trivial output; i.e. 100% of the energy is carried away by neutrinos, which escape freely at the redshifts of interest, so are not treated any further. Even though the channels γ and g are included in our modelindependent approach, these are strongly suppressed for realistic models since DM should not carry color or interact electromagnetically.
For all channels, the spectra of the emerging stable particles, e, p, γ, ν, after treatment of several decays, parton showers, and hadronization were calculated using PYTHIA Monte Carlo^{2} (Sjöstrand et al. 2008). All input spectra are downloadable from http://www.marcocirelli.net/PPPC4DMID.html. For more details, and in particular for a discussion on the level of possible uncertainties, we refer the reader to Cirelli et al. (2011).
Among the stable output particles, neutrinos propagate freely at the redshifts of interest, while e^{ ± } immediately interact with the ubiquitous CMB photons and upscatter those to the gammaray energy range via the inverse Compton (IC) mechanism. The total output energy in hadrons (p and d) is typically quite negligible. Only in gluon and quark channels does it reach up to ~15%, and that only for the highest DM particle masses considered in this paper. We therefore did not model this component in detail. However, in Sect. 4 we give some arguments that at the redshifts of interest, i.e. z ~ 1000, probably the majority of energy released in p and d channels is directly converted into heat by the cosmic medium. Thus, we only need a detailed treatment for the photons: (A) the energetic ones originating directly from the annihilation event (prompt photons); (B) the softer IC photons created by e^{ ± } upscattering CMB photons.
Fig. 1 The redshift z′ where the optical depth for photons reaches unity (i.e., τ(z,z′) = 1) for several “observer’s redshifts”: z = 0,10,500,1000. Here the energy plotted is the photon energy at redshift z. The light gray region corresponds to the DM mass interval considered in this paper. 
The processes and the corresponding cross sections relevant to the propagation of photons through the cosmic medium were taken from Zdziarski & Svensson (1989). Starting from the lowest of energies these include (i) photoionization; (ii) Compton losses (on both bound and free electrons); (iii) pair production on matter; (iv) photonphoton scattering; and (v) pair production on ambient photon fields.
In Fig. 1 we show the redshifts z′ where the optical depth τ(z,z′) = 1 for various redshifts of the observer: z = 0,10,500,1000. As can be seen, for intermediate photon energies (depending on the redshift) there is a wellknown Xray/gammaray energy window where the photons can propagate freely over cosmologically large distances (e.g., see Chen & Kamionkowski 2004). For the cosmological radiation transfer it is crucial that this “transparency window” is properly modeled. Once the photon gets outside of this window (we take it to happen after the first interaction) we assume that the following cascade will be locally absorbed in a very short time. Moreover, the fractions (1 − f_{ion})/3 and (1 + 2f_{ion})/3 of the total absorbed energy, are assumed to be going for ionization and heating, respectively. Here f_{ion} is the fraction of ionized hydrogen atoms, and a similar expression for helium can be used. Excitations of atoms are neglected. This approximation, motivated by the work of Shull & van Steenberg (1985), has been widely used in several subsequent papers, e.g., Chen & Kamionkowski (2004); Padmanabhan & Finkbeiner (2005); Mapelli et al. (2006); Zhang et al. (2006); Natarajan & Schwarz (2009). However, it is clear that these simple expressions only provide a rough estimate for the correct energy deposition efficiencies. As mentioned by Chluba (2010), the precise redshift dependence of the efficiency factors and the fraction of energy that goes into excitations of hydrogen and helium atoms need more careful consideration of the radiative transfer processes, including secondary lowenergy photons and their feedback. These extra photons have the potential of further delaying recombination, hence affect the last scattering surface and CMB anisotropies (Peebles et al. 2000). We leave a more detailed investigation of these ambiguities to a future paper; however, later on in Sect. 4 we briefly comment on how much the omission of excitations using the rough prescription of Chen & Kamionkowski (2004) changes our final results.
In Fig. 1 we also see that at high redshifts, as the density of the environment becomes much higher, the Xray/gammaray transparency window starts to close. Thus at sufficiently high z, we would expect all the produced annihilation energy (excluding the energy stored in neutrinos, since these can freely leak out) to be absorbed locally. In the following we call the ratio of the locally produced to the locally absorbed energy the fparameter. At high redshifts (but well after the neutrino decoupling) we expect the fparameter to asymptote to the value given by (1 − f_{ν}), where f_{ν} is the fraction of energy carried away by neutrinos. An example for the fparameters in the case of μ annihilation channel are shown in the upper lefthand panel of Fig. 2. Since ~60% of the energy is carried away by neutrinos in the μchannel the expected asymptotic highredshift fparameter should be ~0.4, which is indeed the case. With this in mind, we see that robust modelindependent results from the CMB analyses can be obtained for the DM masses below ~100 GeV. This is the reason we concentrate on light WIMPS in this work. For heavier WIMPs, the computation depends on more complicated details for energy absorption.
3. Effect on CMB
Fig. 2 Lefthand column from top to bottom: (i) fparameters for the μchannel assuming m_{DM} = 10 GeV (dark solid lines) and m_{DM} = 100 GeV (light solid lines). In both cases, the lowest curve out of the triple of lines corresponds to only the smooth background, the middle one includes halos with a lower mass cutoff of 10^{6} M_{⊙}, while the top one has 10^{9} M_{⊙}. The dotted lines represent highredshift fits (valid for z ≳ 170) as given by Slatyer et al. (2009). The dashed line shows the CMB visibility function. Here zg(z)dlnz gives the probability that the CMB photon last scattered in the logarithmic redshift interval dlnz centered on redshift z. (ii) Fraction of free electrons as a function of redshift for m_{DM} = 10 GeV and 100 GeV, along with two values for ⟨ σ_{A}υ ⟩ /( ⟨ σ_{A}υ ⟩ _{std}m_{DM} [ GeV] ): 1 and 10. Here ⟨ σ_{A}υ ⟩ _{std} = 3 × 10^{26} cm^{3} s^{1} is the standard thermal cross section. For the meaning of each line, see the description given in the legend. The lowest dotted line represents the standard ΛCDM case with no additional energy input from the annihilating DM. (iii) Matter temperature as a function of redshift. The models shown are exactly the same as in the panel above. The long dashed line compares the temperature of CMB. (iv) The CMB visibility function g(z) times redshift z for the same models as shown in the above two panels. Righthand column from top to bottom: (i) angular power spectra of CMB temperature fluctuations for the same models as already given above. (ii) Temperature and Emode polarization crossspectra. (iii) Emode polarization spectra. In all of the righthand panels the points with errorbars show the 7year measurements by the WMAP space mission. 
Because we are able to calculate fparameters, we can go on to calculate the effect on cosmological recombination. To this end we modify the cosmological recombination code RECFAST (Seager et al. 1999) along the lines presented in Padmanabhan & Finkbeiner (2005). For more details see also Hütsi et al. (2009). Although more advanced cosmological recombination codes have recently been released^{3} (Chluba & Thomas 2010; AliHaïmoud & Hirata 2011), at this point we do not include any of the recently discuss corrections to the cosmological recombination process arising from detailed radiative transfer and atomic physics (see Fendt et al. 2009; RubiñoMartín et al. 2010, and references therein for overview).
The main effect of DM annihilations is a delay in recombination at z ~ 1100 and an increase in the lowredshift freezeout tail. This changes the position and width of the last scattering surface and thus affects the CMB temperature and polarization anisotropies (see, for example, Chen & Kamionkowski 2004; Padmanabhan & Finkbeiner 2005). Since the input annihilation power scales as $\mathrm{\propto}{\mathit{\rho}}_{\mathrm{DM}}^{\mathrm{2}}$, it is clear that structure formation leads to a significant boost over the average density squared $\mathit{\rho \u0305}\begin{array}{}\mathrm{2}\\ \mathrm{DM}\end{array}$, i.e. $\mathrm{\u27e8}{\mathit{\rho}}_{\mathrm{DM}}^{\mathrm{2}}\mathrm{\left(}\mathit{z}\mathrm{\right)}\mathrm{\u27e9}\mathrm{=}\mathit{B}\mathrm{\left(}\mathit{z}\mathrm{\right)}\mathit{\rho \u0305}\begin{array}{}\mathrm{2}\\ \mathrm{DM}\end{array}\mathrm{\left(}\mathit{z}\mathrm{\right)}$. The details of how we calculated the structure boost factors B(z) can be found in Hütsi et al. (2009). The onset of structure formation at z ~ 100 is clearly visible in the upper lefthand panel of Fig. 2. Here the three dark (light) solid curves, representing the results for the μ annihilation channel, correspond to m_{DM} = 10 GeV (100 GeV). From the bottom to top, the trio of lines in each case represent: (i) no structure formation, i.e. only a smooth background; (ii) including structures with the lower halo mass cutoff of 10^{6} M_{⊙}; (iii) structures with the cutoff of 10^{9} M_{⊙}. The halo mass function was assumed to have an analytic ShethTormen form (Sheth & Tormen 1999) and the halo massconcentration relation followed Macciò et al. (2008) description. For the full details of this calculations, again see Hütsi et al. (2009). The dotted lines represent the highredshift fits (valid for z ≳ 170) for the fparameters as given by Slatyer et al. (2009). After adjusting slightly the highredshift normalizations we see that the shape of their fparameters at high redshifts agrees remarkably well with our results, even though our treatment is somewhat more simplified. The sharply peaked dashed curve at z ~ 1000 in Fig. 2 shows the Thomson visibility function (Sunyaev & Zeldovich 1970). More specifically, the quantity gives the probability that the CMB photon last scattered in the logarithmic redshift interval dlnz centered on redshift z. For the CMB calculations the values of the fparameter matter only for the redshift range over which the visibility function is large, and in this range our results agree with the Slatyer et al. (2009) calculations to better than 5% accuracy. We also tested the results for the τchannel and found very good agreement. Thus, our somewhat more simplistic calculation of the fparameters compared to Slatyer et al. (2009) seems to be justified.
The other lefthand panels of Fig. 2 (from top to bottom) show the fraction of free electrons, matter temperature, and the Thomson visibility function as a function of redshift for two DM particle masses, 10 GeV and 100 GeV, along with annihilation cross sections { 10,100 } times, and { 100,1000 } times the standard thermal cross section ⟨ σ_{A}υ ⟩ _{std} = 3 × 10^{26} cm^{3} s^{1}. See the legend for the line definitions. The dotted lines show the standard behavior of ionization fraction, matter temperature and visibility function without any additional energy input from DM annihilation. In the panel for the matter temperature we also show the line for the temperature of the CMB.
Because the annihilation power scales in proportion to $\mathrm{\u27e8}{\mathit{\sigma}}_{\mathrm{A}}\mathit{\upsilon}\mathrm{\u27e9}{\mathit{n}}_{\mathrm{DM}}^{\mathrm{2}}{\mathit{m}}_{\mathrm{DM}}\mathrm{\propto}\mathrm{\u27e8}{\mathit{\sigma}}_{\mathrm{A}}\mathit{\upsilon}\mathrm{\u27e9}\mathit{/}{\mathit{m}}_{\mathrm{DM}}$, we would expect the case with m_{DM} = 10 GeV, and ⟨ σ_{A}υ ⟩ = 10 ⟨ σ_{A}υ ⟩ _{std} (100 ⟨ σ_{A}υ ⟩ _{std}) give comparable results to the case m_{DM} = 100 GeV and ⟨ σ_{A}υ ⟩ = 100 ⟨ σ_{A}υ ⟩ _{std} (1000 ⟨ σ_{A}υ ⟩ _{std}), as long as the fparameter does not depend strongly on m_{DM}. This is indeed approximately so, as can be seen from Fig. 2. Even though the lines for the ionization fraction, matter temperature, and CMB visibility function are easily separated at lower redshifts, this simple scaling with ⟨ σ_{A}υ ⟩ /m_{DM} seems to hold very well near the peak of the visibility function, and so one would not expect to be able to clearly distinguish the models with the same value of ⟨ σ_{A}υ ⟩ /m_{DM} via CMB measurements alone^{4}. Also, the additional annihilation boost from the structure formation is not expected to influence the CMB measurements. Only in the very extreme case where instead of Macciò et al. (2008) massconcentration relation we use a simple powerlaw extrapolation down to the very low halo masses one is able to cause an additional significantly high lowredshift peak in the visibility function. However, for any realistic massconcentration relation, along with annihilation cross sections ⟨ σ_{A}υ ⟩ that do not violate the CMB data, the contribution from the structure formation to the CMB signal is completely negligible. Also, the contribution to the lowredshift ionization fraction is very mild, and thus the annihilating DM models which are compatible with CMB measurements could only play a marginal role in reionizing the lowz Universe. These results support our similar findings previously reported in Hütsi et al. (2009).
In the righthand panels of Fig. 2, we show the temperaturetemperature (TT), temperatureEmode polarization (TE), and EmodeEmode (EE) angular power spectra. The points with error bars show the WMAP sevenyear measurements (Larson et al. 2011). The dotted lines in all of the panels show theoretical predictions for the concordance ΛCDM model (Komatsu et al. 2011). As we might already expect, remembering the visibility function behavior as shown in the lowest lefthand panel, the models with the same values for ⟨ σ_{A}υ ⟩ /m_{DM} are hardly distinguishable. This is because the redshift dependence of the fparameter over the width of the visibility function is small, and so only the average value around z ~ 1100 really matters, in agreement with the statements of Slatyer et al. (2009).
As a final note in this section we point out that, if one increases ⟨ σ_{A}υ ⟩ to high enough values and adds the structure formation boost, so that the lowredshift ionization fraction starts to get significantly approach one, the original RECFAST code, along with its DVERK solver, is not properly able to deal with the underlying set of stiff ODEs and simply breaks down. To be able to numerically treat these cases, one can use the significantly improved Recfast++ code^{5} (Chluba & Thomas 2010), along with its stiff ODE solver. However, in our calculations the original RECFAST works fine, since the typically allowed values of ⟨ σ_{A}υ ⟩ /m_{DM} compatible with the CMB data are low enough.
4. CMB constraints
Fig. 3 The WMAP7 parameter constraints for the μ annihilation channel. Along with two annihilation parameters log _{10}(m_{DM}) and ⟨ σ_{A}υ ⟩ /m_{DM}, we also show the constraints on the other 6parameter ΛCDM background model: Ω_{B}h^{2}, Ω_{DM}h^{2}, H_{0}, τ, n_{S}, log _{10}(A_{S}). From the inside out the colored 2D areas show the 1sigma and 2sigma regions after marginalization over the other parameters, and the dashed lines show the same regions for the basic 6parameter ΛCDM without annihilating DM. The topmost panels in each column plot the marginalized 1D probability distributions for all of the parameters. 
Now we have all the ingredients available to calculate current and future CMB constraints on annihilating DM. We notice that the CMB constraints for the decaying DM are not competitive (e.g., when compared to the bounds available from the lowredshift gammaray measurements) due to ∝ ρ_{DM} scaling when compared to the $\mathrm{\propto}{\mathit{\rho}}_{\mathrm{DM}}^{\mathrm{2}}$ scaling for the annihilating DM, and thus are not being discussed in this work.
For the currently existing CMB data, we used the latest power spectra, i.e. sevenyear spectra, from the WMAP space mission (Larson et al. 2011). In reality the current CMB bounds can be somewhat tightened if one includes results from other smaller scale CMB measurements, in particular the measurements of the temperature power spectrum from ACBAR (Reichardt et al. 2009) and temperature and polarization spectra from QUaD (Pryke et al. 2009). In addition, data from ACT (Das et al. 2011) and SPT (Lueker et al. 2010) could be added to the analysis, providing additional leverage on small scales^{6}. However, to keep our analysis clearer we decided to restrict ourselves to the WMAP7 data only. Regarding future CMB results we made predictions for the bounds available from the currently ongoing Planck mission, and also for the idealized noisefree experiment up to the multipole of ℓ_{max} = 2000. As in Zaldarriaga et al. (2008), for Planck we assume 80% sky coverage, beam size θ_{FWHM} = 7 arcmin, and noisebias ${\mathit{N}}_{\mathit{\ell}}^{\mathrm{TT}}\mathrm{\simeq}{\mathrm{10}}^{4}$ μK^{2} and ${\mathit{N}}_{\mathit{\ell}}^{\mathrm{EE}}\mathrm{\simeq}\mathrm{3.5}$ × 10^{4} μK^{2} for TT and EE, respectively. For the ideal noisefree experiment, we assume a full sky coverage and the underlying uncertainties of the CMB fluctuations are solely due to a finite number of available fluctuation modes on the sky, i.e. due to cosmic variance, and so we call this type of idealized experiment “cosmic variance limited” (CVL) in the following. First steps towards CVL CMB measurements in both temperature and polarization down to small scales will become available from a combinations of SPTpol^{7} (McMahon et al. 2009) and ACTPol^{8} (Niemack et al. 2010) along with Planck data, so that the CVL case discussed here provides a good guideline to what could become possible in the near future.
We generated the synthetic data for the Planck and CVL experiments as described in http://cosmocoffee.info/^{9}. Along with WMAP7 likelihood code^{10}, they will be used in our Markov chain Monte Carlo (MCMC) parameter estimation. For the MCMC engine we used the publicly available CosmoMC tools^{11}(Lewis & Bridle 2002), where the standard recombination modules were modified in order to allow additional energy input from annihilating DM. Our baseline concordance ΛCDM model (Komatsu et al. 2011) is parameterized by six free parameters: Ω_{B}h^{2}, Ω_{DM}h^{2}, H_{0}, τ, n_{S}, log _{10}(A_{S}), where Ω_{B} and Ω_{DM} are density parameters for baryons and DM, H_{0} = 100 h km s^{1} Mpc^{1} is the Hubble constant, τ is the reionization optical depth, and A_{S} and n_{S} are the amplitude and spectral index for adiabatic scalar perturbations, respectively. This minimal sixparameter set is extended by two additional parameters describing the annihilating DM: its mass m_{DM}, and thermally averaged annihilation cross section ⟨ σ_{A}υ ⟩ . Since MCMC will give better results if model parameters are as uncorrelated as possible, we used the parameter pair { m_{DM}, ⟨ σ_{A}υ ⟩ /m_{DM} } in place of { m_{DM}, ⟨ σ_{A}υ ⟩ } .
As an example, we show the WMAP7 constraints on all of the eight parameters, assuming μ annihilation channel in Fig. 3. With dashed lines we show parameter constraints for the sixparameter ΛCDM model without annihilating DM. The colored 2D areas in all of the panels display 1sigma and 2sigma regions marginalized over all of the other parameters. For all of the parameters the topmost panels in each column plot the marginalized 1D probability distributions. It is reassuring to see that after introducing two additional parameters, all the previous six parameters, albeit with small shifts, are still as precisely determined. The strongest shifts compared to the baseline sixparameter model are quite understandably seen for parameters A_{S} and n_{S}, which both lead to the rise of the CMB power spectra if increased, and thus counterbalance the damping of the spectrum caused by the extra scattering off the additional free electrons created by the energy input from the DM annihilation.
Here it is important to note that these small shifts towards higher values of A_{S} and n_{S} are directed in the same direction as recently discussed recombination corrections from previously neglected standard physics processes (RubiñoMartín et al. 2010; Shaw & Chluba 2011). Although biases in A_{S} and n_{S} arising from recombination physics are not significant for WMAP7 (Shaw & Chluba 2011), parameters derived from Planck data will be affected by several standard deviations if recombination corrections are neglected. This indicates that precise constraints on models with annihilating DM should probably account for these recombination corrections simultaneously.
Fig. 4 Upper panel: fparameters for all of the channels considered in this work and for several DM particle masses in the interval m_{DM} = 5 − 100 GeV. For all of the cases, the asymptotic highredshift fparameter values have been renormalized to be equal to one. Lower panel: the same as above, after dividing with a “typical shape”, (f_{max}(z) + f_{min}(z))/2, of the fparameter curve. 
In Fig. 3 we only present the constraints for one particular channel. However, we would like to provide easily usable recipes to calculate bounds for all of the channels described in Sect. 2. As seen before, even though the fparameters vary somewhat if one changes the mass from 10 GeV to 100 GeV, the CMB spectra as shown in Fig. 2 were fairly insensitive to this change as long as the ratio ⟨ σ_{A}υ ⟩ /m_{DM} was kept the same. In the upper panel of Fig. 4 we plot the fparameters for all of the channels and for several masses in the range m_{DM} = 5 − 100 GeV. Here we have renormalized the highz values to be equal to one. In the redshift range where the visibility function peaks (see the lowest lefthand panel of Fig. 2), the variation among different DM models is not too large. To see this better, in the lower panel of Fig. 4 we have divided the above lines by a “typical shape” of the fparameter given by (f_{max}(z) + f_{min}(z))/2, where f_{max}(z) and f_{min}(z) are the extremal values of f(z) for a particular z. We see that, at most relevant redshifts, the variation among all the different models is around 15%. This indeed suggests that it might be possible to provide a unified approximate scheme for calculating the bounds for all of the annihilation channels.
To bracket all the possibilities, we calculated the bounds assuming the fparameter to have a typical shape (f_{max}(z) + f_{min}(z))/2 along with two extremal values f_{min}(z) and f_{max}(z). Of course, instead of two annihilation parameters we then only have to deal with one additional parameter: ⟨ σ_{A}υ ⟩ /m_{DM}. We took the asymptotic values for the fparameter to have values of 0.25,0.50,0.75, and 1.00. This asymptotic value, as explained above, is simply (1 − f_{ν}), where f_{ν} is the fraction of the total energy carried away by neutrinos. It turned out that the resulting 1sigma and 2sigma constraints on (1 − f_{ν}) ⟨ σ_{A}υ ⟩ /m_{DM} for the above four values of (1 − f_{ν}) are the same with (1 − 2)% accuracy. Our final results (averaged over these small variations due to changing (1 − f_{ν})) can thus be given in the form $\mathrm{(}\mathrm{1}\mathrm{}{\mathit{f}}_{\mathit{\nu}}\mathrm{)}\frac{\mathrm{\u27e8}{\mathit{\sigma}}_{\mathrm{A}}\mathit{\upsilon}\mathrm{\u27e9}\hspace{0.17em}\mathrm{[}\mathrm{3}\mathrm{\times}{\mathrm{10}}^{26}{\mathrm{cm}}^{\mathrm{3}}{\mathrm{s}}^{1}\mathrm{]}}{{\mathit{m}}_{\mathrm{DM}}\hspace{0.17em}\left[\mathrm{GeV}\right]}\mathit{}\mathit{r,}$(1)where the values of r for 1sigma and 2sigma bounds are given in Table 1. One can see that the approximate bounds we provide are typically accurate at about the 20 − 30% level. We tested that, indeed, the ranges of r given in Table 1 fully cover the values of r for all of the channels considered in this work. Two examples for the case of 1sigma upper bounds are shown in Fig. 5. Here the upper panel corresponds to the μchannel (f_{ν} ≃ 0.61) and the lower one to the echannel (f_{ν} ≃ 0.02). We chose the above two channels because these are the two extreme cases among all of the channels treated in this paper. Consequently, they bracket the expected results for any realistic annihilating DM, given as a superposition over the basis channels. The solid lines show the bounds calculated directly through the MCMC analysis, while the shaded regions represent the ranges as obtained from Eq. (1) and Table 1. Indeed, the solid lines are fully covered by the shaded regions, as it should be. The vertical gray stripe indicates the range of WIMP masses (m_{DM} = 6 − 8 GeV) that provide a good fit to CoGeNT and DAMA/LIBRA data (Hooper et al. 2010). The vertical dotted line marks the lowest DM particle mass 5 GeV used in our PYTHIA simulations. This cutoff is not physical, but occurs because PYTHIA does not work below that energy. Therefore, the extrapolations of our results are shown below 5 GeV DM mass. Note that directly calculated lines for the μ and echannel have slightly steeper slopes than given by the shaded regions, which increase as ∝ m_{DM}. As the fparameters generally fall off more slowly for lower m_{DM}, this behavior is also typical of the other channels. Thus, for lower m_{DM} values, one should actually get slightly stronger bounds than calculated directly from Eq. (1), and so our extrapolations shown in Fig. 5 are somewhat conservative.
Fig. 5 WMAP, Planck, and CVL 1sigma constraints on the ⟨ σ_{A}υ ⟩ − m_{DM} plane for μ (upper panel) and e (lower panel) annihilation channels. The solid lines show the upper bounds on annihilation cross section as determined directly through full MCMC calculations. The shaded regions around solid lines show the results from the simple recipe of Eq. (1) with values of r taken from Table 1. The vertical gray stripe shows the range of WIMP masses (m_{DM} = 6 − 8 GeV) that provide a good fit to CoGeNT and DAMA/LIBRA data (Hooper et al. 2010). The vertical dotted line marks the lowest DM particle mass 5 GeV available for PYTHIA simulations. Extrapolations are shown below that value. 
We see that for WMAP7, depending on the annihilation model, the limiting DM particle mass below which the upper bound on the annihilation cross section drops below the standard thermal production value is in the range 4.5 − 10 GeV, while the corresponding numbers reachable for the Planck and CVL experiment are 19 − 43 GeV and (45 − 100) GeV, respectively. Since the μchannel represents our most conservative case, one can instead say that, according to currently available CMB data, the annihilation cross section should be below the standard value of 3 × 10^{26} cm^{3} s^{1} as long as m_{DM} ≲ 5 GeV^{12}. Thus, for the CoGeNT and DAMA/LIBRA bestfit mass region 6 − 8 GeV, the standard thermal production cross section is still compatible with the CMB measurements. This could quite possibly be changed soon as Planck results become available. Of course one should keep in mind that the cross section bounds given here directly apply to redshifts of z ~ 1000, and if the cross section depends on velocity (as in the case of Pwave annihilation or Sommerfeldenhanced scenario), one should be careful in converting these numbers to the values relevant at z = 0. Also, note that the standard annihilation cross section 3 × 10^{26} cm^{3} s^{1}, with respect to what we are comparing our CMB bounds to, provides the desired thermal relic density (i.e., Ω_{DM} ~ 0.3) only if WIMPs are annihilating through Swave processes.
Except for a typical assumption of Swave annihilation our results are largely modelindependent, so all the particlephysics scenarios motivated by DAMA/LIBRA and CoGeNT results (Fitzpatrick et al. 2010; Andreas et al. 2010; Chang et al. 2010; Foot 2010; Barger et al. 2010b; Hooper et al. 2010; Fitzpatrick & Zurek 2010; Essig et al. 2010; Barger et al. 2010a; Cline et al. 2011; Buckley et al. 2011) including theoretically wellmotivated particle physics models that predict light DM, such as the MSSM (Feldman et al. 2010; Kuflik et al. 2010; Fornengo et al. 2011; Bottino et al. 2010) and the NMSSM (Kang et al. 2011; Belikov et al. 2011; Gunion et al. 2010; Draper et al. 2011), are stringently tested by WMAP and Planck.
A recent study by Dobler et al. (2011) suggests that annihilating DM with masses in the range m_{DM} = 1 − 20 GeV (i.e. compatible with CoGeNT and DAMA bestfit region), along with annihilation cross section ⟨ σ_{A}υ ⟩ ~ 9 × 10^{25} cm^{3} s^{1} (i.e. 30 times the standard thermal relic cross section), may give a good match to the observed Fermi and WMAP haze. From Fig. 5 it is clear that, for the DM particle masses in this range, such high annihilation cross sections certainly conflict with available CMB data. If the required boost by a factor of ~ 30 is obtained via the Sommerfeld enhancement (Cirelli et al. 2009b; ArkaniHamed et al. 2009; Slatyer 2010), then at z ~ 1000 the expected cross section is at least as large as the value given above, so the model could be in even bigger trouble. The alternative claims (Goodenough & Hooper 2009; Hooper & Goodenough 2011) of light DM annihilations in our Galaxy will also be stringently tested.
Finally, it would be interesting to compare the bounds given in Slatyer et al. (2009), which are based on the results obtained by Galli et al. (2009), with our values. However, since Galli et al. (2009) used “on the spot” approximation for the energy deposition, the comparison can only be approximate. The WMAP5 2sigma bound as given by Eq. (6) of Slatyer et al. (2009) can be cast in the form ⟨ σ_{A}υ ⟩ ^{ [ }3 × 10^{26} cm^{3} s^{1}^{] }/m_{DM} [GeV] < 0.12/f. To compare this with our WMAP7 results, we have to set f_{ν} = 0 in Eq. (1) and, in the above relation, use a typical value for the fparameter at z ~ 1000, which from Fig. 4 is f ~ 0.85. Thus, the value for r that should be directly comparable to the 2sigma WMAP value given in Table 1 is 0.12/0.85 ~ 0.14. Considering the differences in the treatment for the energy deposition, this value agrees reasonably well with our result. However, there are significantly greater differences if one compares the forecasts for the Planck and CVL experiments. In our case Planck and CVL would tighten the 2sigma bound of WMAP by a factor of ~6 and ~14, respectively. The corresponding numbers (~13 and ~40, respectively) from Galli et al. (2009) are certainly more optimistic. For CVL, some of this discrepancy is surely due to the higher ℓ_{max} value assumed in Galli et al. (2009): ℓ_{max} = 2500 compared to our ℓ_{max} = 2000. Even though we are comparing here the results derived from WMAP5 and WMAP7, the difference is expected to be negligible because the accuracy of the cosmological parameters with an additional two years of WMAP data improves only very moderately (typically somewhere between 5 − 15%), and thus cannot be the cause of the above discrepancy.
NOTE: After the first version of this work was submitted, a new paper by Galli et al. (2011) appeared where the authors have performed a similar analysis (now also using WMAP7 data), this time finding results that agree with ours more closely: typical deviations are now within factors 1.2 − 2, with their bounds on annihilation cross section being somewhat stronger. These relatively small deviations are actually quite negligible keeping in mind several approximations used in both analyses.
To see how much the omission of the additional excitations of atoms could change our results, we also performed several calculations where the exciting Lyα photons are treated along the lines presented in Chen & Kamionkowski (2004). In agreement with the claim in Galli et al. (2011), we find that the bounds on annihilation cross section are getting tighter only up to ~10%; i.e., the excitations seem to have only a relatively weak effect with the current prescription.
Note on protons
Although for most of the channels and, in particular, for the lower end of the considered DM particle mass range, the number of protons produced is completely negligible, for quark and gluon channels the contribution can reach ~15% of the total energy input if one has masses at the higher end of the considered range. As this is still only a relatively mild contribution, we did not attempt any detailed modeling for the proton component and simply assumed that all of this energy is absorbed as heat by the cosmic medium. A simple justification is the following. The energy loss rate for protons with the energies of interest in this paper is dominated by protonproton scattering. Thus the loss rate $\mathrm{\Gamma}\mathrm{\equiv}\mathrm{}\frac{\mathrm{1}}{\mathit{E}}\frac{\mathrm{d}\mathit{E}}{\mathrm{d}\mathit{t}}\mathrm{\simeq}{\mathit{n}}_{\mathrm{p}}\mathit{c}{\mathit{\sigma}}_{\mathrm{pp}}{\mathit{K}}_{\mathrm{pp}}$, where n_{p} is the number density of target protons, σ_{pp} the scattering cross section, and K_{pp} the inelasticity parameter. The cross section σ_{pp} depends only weakly on proton energy with typical values of 30 − 40 mbarn and inelasticity parameter K_{pp} ≃ 0.5 for the energies of interest in this work (see e.g., Wommer et al. 2008). At redshifts z ~ 1000 we therefore get Γ ≃ 10^{13} s^{1} for the energy loss rate. Comparing this to the expansion rate H at the same redshift, H ≃ 0.4 × 10^{13} s^{1}, we see that Γ ≳ H, so one might expect a significant fraction of the proton energy to be absorbed.
5. Summary
In this paper we have calculated the existing and future CMB constraints on annihilating DM assuming the Swave annihilation mode to be dominant. Our results can be summarized as follows.

In full agreement with our earlier findings, first presented in Hütsiet al. (2009), we confirm that in caseof any realistic halo massconcentration relation, along withannihilation cross sections ⟨ σ_{A}υ ⟩ that do not violate the CMB data, the contribution from structure formation to the CMB signal is completely negligible. Also, the contribution to the lowredshift ionization fraction is very mild, and thus the annihilating DM models compatible with CMB measurements could only play a minor role in helping to reionize the lowz Universe. Therefore for lowmass DM the CMB constraints can be considered almost free of cosmological uncertainties.

At the same time, the particle physics uncertainties on the DM annihilation channels can all be described with one single parameter f_{ν}, the energy fraction carried away by neutrinos. Therefore the DM annihilations to two muons represents the least stringently tested DM scenario.

As our main results, in the form of Eq. (1) and Table 1, we provided a simple recipe for estimating the upper bounds on annihilation cross section, valid for a broad range of DM models. Two examples for 1sigma upper bounds for the μ and echannels (which are the two extremal cases among all of the channels considered in this paper; i.e., those upper bounds bracket the expected results for more general annihilating DM models, given as a superposition over the basis channels) are shown in Fig. 5.

For the DM particle masses m_{DM} = 6 − 8 GeV, which give best fits to CoGeNT and DAMA/LIBRA data, current CMB data is still compatible with a standard thermal relic annihilation cross section 3 × 10^{26} cm^{3} s^{1} only if the annihilations are dominantly into the μ or τ (and corresponding 4 lepton) channels. All other annihilation channels already now conflict with the CMB data.

The sensitivity of Planck space mission allows one to test all the DM annihilation channels definitively. If the Planck mission will not find signals of annihilating DM in the CMB spectrum, all the light DM scenarios motivated by DAMA/LIBRA, CoGeNT, Fermi haze, and WMAP haze will be stringently tested if DM is a thermal relic.

Our findings strongly disfavor a claim that annihilating DM with m_{DM} ~ 1−20 GeV and ⟨ σ_{A}υ ⟩ ~ 9 × 10^{25} cm^{3} s^{1} (i.e. ~30 times above the standard value) could be a cause for Fermi and WMAP haze.

The last conclusion applies to all DM scenarios with large boost factors from the Sommerfeld enhancement. Because CMB is sensitive to DM annihilations at z ~ 1000, at that time the Sommerfeld enhancement of annihilation cross section should have been at least as large as today. Therefore, for light DM scenarios the boost of DM annihilation in our Galaxy can only come form the halo substructure.

Our MCMC calculations assumed standard 6parameter ΛCDM cosmology, which was extended by two additional parameters describing the annihilating DM: m_{DM} and ⟨ σ_{A}υ ⟩ /m_{DM}. It turns out that, after introducing these 2 additional degrees of freedom, all the previous 6 parameters, albeit with small shifts (the most noticeable of which being the ones for A_{S} and n_{S}), were still as precisely determined.

The weakest point of the analysis presented in this paper is the frequently used simple approximation about how the input energy of the cascading particles gets partitioned between ionizations and heating of the environment. For better treatment of this issue one possibly has to rely on detailed Monte Carlo calculations. This, however, we leave for a possible future study.
CosmoRec (www.Chluba.de/CosmoRec) already provides a simple module that accounts for the effect of DM annihilation or more general energy injection, but this module is still being validated.
To achieve better distinction between various models, one might try to use the differences in the behavior of the lowz matter temperature and its observable consequence on the 21 cm transition measurements of the neutral hydrogen. The effect of annihilating DM on the 21 cm signal has been discussed e.g., in Mapelli et al. (2006); Furlanetto et al. (2006); Valdés et al. (2007); Natarajan & Schwarz (2009).
Recfast++ is part of more advanced recombination code CosmoRec and can be downloaded from www.Chluba.de/CosmoRec
See also http://lpsc.in2p3.fr/perotto/ in case CMB lensing is needed.
Acknowledgments
We thank Marco Cirelli, Mario Kadastik, and Alessandro Strumia for discussions and our referee for useful comments and suggestions. G.H. acknowledges the support provided through a visiting scientist fellowship at the MPA. This work was supported by the ESF ERMOS Postdoctoral Research Grant 35, SF0060067s08, SF0690030s09, ESF8005, ESF8943, ESF8090, MTT8, MJD52, and by EU FP7INFRA20071.2.3 contract No. 223807. J.C. is also very grateful for additional financial support from the Beatrice D. Tremaine fellowship 2010.
References
 Aalseth, C. E., Barbeau, P. S., Bowden, N. S., et al. 2011, Phys. Rev. Lett., 106, 131301 [Google Scholar]
 Abdo, A. A., Ackermann, M., Ajello, M., et al. 2009, Phys. Rev. Lett., 102, 181101 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
 Adriani, O., Barbarino, G. C., Bazilevskaya, G. A., et al. 2009, Nature, 458, 607 [CrossRef] [PubMed] [Google Scholar]
 Aharonian, F., Akhperjanian, A. G., Anton, G., et al. 2009, A&A, 508, 561 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Ahmed, Z., Akerib, D. S., Arrenberg, S., et al. 2011, Phys. Rev. Lett., 106, 131302 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
 AliHaïmoud, Y., & Hirata, C. M. 2011, Phys. Rev. D, 83, 043513 [NASA ADS] [CrossRef] [Google Scholar]
 Andreas, S., Arina, C., Hambye, T., Ling, F., & Tytgat, M. H. G. 2010, Phys. Rev. D, 82, 043522 [NASA ADS] [CrossRef] [Google Scholar]
 Angle, J., Aprile, E., Arneodo, F., et al. 2009, Phys. Rev. D, 80, 115005 [NASA ADS] [CrossRef] [Google Scholar]
 Aprile, E., Arisaka, K., Arneodo, F., et al. 2010, Phys. Rev. Lett., 105, 131302 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
 Arina, C., & Tytgat, M. H. G. 2011, J. Cosmology Astropart. Phys., 1, 11 [NASA ADS] [CrossRef] [Google Scholar]
 ArkaniHamed, N., Finkbeiner, D. P., Slatyer, T. R., & Weiner, N. 2009, Phys. Rev. D, 79, 015014 [NASA ADS] [CrossRef] [Google Scholar]
 Barger, V., Keung, W., Marfatia, D., & Shaughnessy, G. 2009, Phys. Lett. B, 672, 141 [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
 Barger, V., Gao, Y., McCaskey, M., & Shaughnessy, G. 2010a, Phys. Rev. D, 82, 095011 [NASA ADS] [CrossRef] [Google Scholar]
 Barger, V., McCaskey, M., & Shaughnessy, G. 2010b, Phys. Rev. D, 82, 035019 [NASA ADS] [CrossRef] [Google Scholar]
 Baxter, E. J., Dodelson, S., Koushiappas, S. M., & Strigari, L. E. 2010, Phys. Rev. D, 82, 123511 [NASA ADS] [CrossRef] [Google Scholar]
 Belikov, A. V., Gunion, J. F., Hooper, D., & Tait, T. M. P. 2011, Phys. Lett. B, 705, 82 [NASA ADS] [CrossRef] [Google Scholar]
 Bell, N. F., & Jacques, T. D. 2009, Phys. Rev. D, 79, 043507 [NASA ADS] [CrossRef] [Google Scholar]
 Bergström, L., Bringmann, T., & Edsjö, J. 2008, Phys. Rev. D, 78, 103520 [CrossRef] [Google Scholar]
 Bergström, L., Bertone, G., Bringmann, T., Edsjö, J., & Taoso, M. 2009, Phys. Rev. D, 79, 081303 [NASA ADS] [CrossRef] [Google Scholar]
 Bernabei, R., Belli, P., Cappella, F., et al. 2010, Eur. Phys. J. C, 67, 39 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Bertone, G., Hooper, D., & Silk, J. 2005, Phys. Rep., 405, 279 [NASA ADS] [CrossRef] [Google Scholar]
 Bertone, G., Cirelli, M., Strumia, A., & Taoso, M. 2009, J. Cosmology Astropart. Phys., 3, 9 [Google Scholar]
 Bottino, A., Donato, F., Fornengo, N., & Scopel, S. 2010, Phys. Rev. D, 81, 107302 [NASA ADS] [CrossRef] [Google Scholar]
 Buckley, M. R., Hooper, D., & Tait, T. M. P. 2011, Phys. Lett. B, 702, 216 [NASA ADS] [CrossRef] [Google Scholar]
 Chang, S., Liu, J., Pierce, A., Weiner, N., & Yavin, I. 2010, J. Cosmology Astropart. Phys., 8, 18 [Google Scholar]
 Chen, X., & Kamionkowski, M. 2004, Phys. Rev. D, 70, 043502 [NASA ADS] [CrossRef] [Google Scholar]
 Chluba, J. 2010, MNRAS, 402, 1195 [NASA ADS] [CrossRef] [Google Scholar]
 Chluba, J., & Thomas, R. M. 2010, MNRAS, 1876 [Google Scholar]
 Cirelli, M., & Panci, P. 2009, Nucl. Phys. B, 821, 399 [Google Scholar]
 Cirelli, M., Iocco, F., & Panci, P. 2009a, J. Cosmology Astropart. Phys., 10, 9 [Google Scholar]
 Cirelli, M., Kadastik, M., Raidal, M., & Strumia, A. 2009b, Nucl. Phys. B, 813, 1 [NASA ADS] [CrossRef] [Google Scholar]
 Cirelli, M., Panci, P., & Serpico, P. D. 2010, Nucl. Phys. B, 840, 284 [Google Scholar]
 Cirelli, M., Corcella, G., Hektor, A., et al. 2011, J. Cosmology Astropart. Phys., 3, 51 [Google Scholar]
 Cline, J. M., Frey, A. R., & Chen, F. 2011, Phys. Rev. D, 83, 083511 [NASA ADS] [CrossRef] [Google Scholar]
 Das, S., Marriage, T. A., Ade, P. A. R., et al. 2011, ApJ, 729, 62 [NASA ADS] [CrossRef] [Google Scholar]
 Dobler, G., & Finkbeiner, D. P. 2008, ApJ, 680, 1222 [NASA ADS] [CrossRef] [Google Scholar]
 Dobler, G., Finkbeiner, D. P., Cholis, I., Slatyer, T., & Weiner, N. 2010, ApJ, 717, 825 [NASA ADS] [CrossRef] [Google Scholar]
 Dobler, G., Cholis, I., & Weiner, N. 2011, ApJ, 741, 25 [NASA ADS] [CrossRef] [Google Scholar]
 Draper, P., Liu, T., Wagner, C. E. M., Wang, L.T., & Zhang, H. 2011, Phys. Rev. Lett., 106, 121805 [NASA ADS] [CrossRef] [MathSciNet] [PubMed] [Google Scholar]
 Einasto, J. 2009 [arXiv:0901.0632] [Google Scholar]
 Essig, R., Kaplan, J., Schuster, P., & Toro, N. 2010 [arXiv:1004.0691] [Google Scholar]
 Feldman, D., Liu, Z., & Nath, P. 2010, Phys. Rev. D, 81, 117701 [NASA ADS] [CrossRef] [Google Scholar]
 Fendt, W. A., Chluba, J., RubiñoMartín, J. A., & Wandelt, B. D. 2009, ApJS, 181, 627 [NASA ADS] [CrossRef] [Google Scholar]
 Feng, J. L. 2010, ARA&A, 48, 495 [Google Scholar]
 Fitzpatrick, A. L., & Zurek, K. M. 2010, Phys. Rev. D, 82, 075004 [NASA ADS] [CrossRef] [Google Scholar]
 Fitzpatrick, A. L., Hooper, D., & Zurek, K. M. 2010, Phys. Rev. D, 81, 115005 [NASA ADS] [CrossRef] [Google Scholar]
 Foot, R. 2010, Phys. Lett. B, 692, 65 [NASA ADS] [CrossRef] [Google Scholar]
 Fornengo, N., Scopel, S., & Bottino, A. 2011, Phys. Rev. D, 83, 015001 [NASA ADS] [CrossRef] [Google Scholar]
 Fox, P. J., & Poppitz, E. 2009, Phys. Rev. D, 79, 083528 [NASA ADS] [CrossRef] [Google Scholar]
 Furlanetto, S. R., Oh, S. P., & Pierpaoli, E. 2006, Phys. Rev. D, 74, 103502 [NASA ADS] [CrossRef] [Google Scholar]
 Galli, S., Iocco, F., Bertone, G., & Melchiorri, A. 2009, Phys. Rev. D, 80, 023505 [NASA ADS] [CrossRef] [Google Scholar]
 Galli, S., Iocco, F., Bertone, G., & Melchiorri, A. 2011, Phys. Rev. D, 84, 027302 [NASA ADS] [CrossRef] [Google Scholar]
 Goodenough, L., & Hooper, D. 2009 [arXiv:0910.2998] [Google Scholar]
 Gunion, J. F., Belikov, A. V., & Hooper, D. 2010 [arXiv:1009.2555] [Google Scholar]
 Hooper, D., & Goodenough, L. 2011, Phys. Lett. B, 697, 412 [NASA ADS] [CrossRef] [Google Scholar]
 Hooper, D., Collar, J. I., Hall, J., McKinsey, D., & Kelso, C. M. 2010, Phys. Rev. D, 82, 123509 [NASA ADS] [CrossRef] [Google Scholar]
 Hütsi, G., Hektor, A., & Raidal, M. 2009, A&A, 505, 999 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Hütsi, G., Hektor, A., & Raidal, M. 2010, J. Cosmology Astropart. Phys., 7, 8 [Google Scholar]
 Jungman, G., Kamionkowski, M., & Griest, K. 1996, Phys. Rep., 267, 195 [NASA ADS] [CrossRef] [Google Scholar]
 Kang, Z., Li, T., Liu, T., Tong, C., & Yang, J. M. 2011, J. Cosmology Astropart. Phys., 1, 28 [NASA ADS] [CrossRef] [Google Scholar]
 Kanzaki, T., Kawasaki, M., & Nakayama, K. 2010, Prog. Theor. Phys., 123, 853 [NASA ADS] [CrossRef] [Google Scholar]
 Komatsu, E., Smith, K. M., Dunkley, J., et al. 2011, ApJS, 192, 18 [NASA ADS] [CrossRef] [Google Scholar]
 Kopp, J., Schwetz, T., & Zupan, J. 2010, J. Cosmology Astropart. Phys., 2, 14 [Google Scholar]
 Kuflik, E., Pierce, A., & Zurek, K. M. 2010, Phys. Rev. D, 81, 111701 [NASA ADS] [CrossRef] [Google Scholar]
 Larson, D., Dunkley, J., Hinshaw, G., et al. 2011, ApJS, 192, 16 [NASA ADS] [CrossRef] [Google Scholar]
 Lewis, A., & Bridle, S. 2002, Phys. Rev. D, 66, 103511 [NASA ADS] [CrossRef] [Google Scholar]
 Lueker, M., Reichardt, C. L., Schaffer, K. K., et al. 2010, ApJ, 719, 1045 [NASA ADS] [CrossRef] [Google Scholar]
 Macciò, A. V., Dutton, A. A., & van den Bosch, F. C. 2008, MNRAS, 391, 1940 [NASA ADS] [CrossRef] [Google Scholar]
 Mapelli, M., Ferrara, A., & Pierpaoli, E. 2006, MNRAS, 369, 1719 [NASA ADS] [CrossRef] [Google Scholar]
 McMahon, J. J., Aird, K. A., Benson, B. A., et al. 2009, in AIP Conf. Ser. 1185, ed. B. Young, B. Cabrera, & A. Miller, 511 [Google Scholar]
 Meade, P., Papucci, M., Strumia, A., & Volansky, T. 2010, Nucl. Phys. B, 831, 178 [Google Scholar]
 Natarajan, A., & Schwarz, D. J. 2009, Phys. Rev. D, 80, 043529 [NASA ADS] [CrossRef] [Google Scholar]
 Niemack, M. D., Ade, P. A. R., Aguirre, J., et al. 2010, in SPIE Conf. Ser., 7741 [Google Scholar]
 Padmanabhan, N., & Finkbeiner, D. P. 2005, Phys. Rev. D, 72, 023508 [NASA ADS] [CrossRef] [Google Scholar]
 Papucci, M., & Strumia, A. 2010, J. Cosmology Astropart. Phys., 3, 14 [Google Scholar]
 Peebles, P. J. E., Seager, S., & Hu, W. 2000, ApJ, 539, L1 [NASA ADS] [CrossRef] [Google Scholar]
 Pryke, C., Ade, P., Bock, J., et al. 2009, ApJ, 692, 1247 [NASA ADS] [CrossRef] [Google Scholar]
 Reichardt, C. L., Ade, P. A. R., Bock, J. J., et al. 2009, ApJ, 694, 1200 [NASA ADS] [CrossRef] [Google Scholar]
 RubiñoMartín, J. A., Chluba, J., Fendt, W. A., & Wandelt, B. D. 2010, MNRAS, 403, 439 [NASA ADS] [CrossRef] [Google Scholar]
 Savage, C., Gelmini, G., Gondolo, P., & Freese, K. 2011, Phys. Rev. D, 83, 055002 [NASA ADS] [CrossRef] [Google Scholar]
 Schwetz, T. 2010 [arXiv:1011.5432] [Google Scholar]
 Seager, S., Sasselov, D. D., & Scott, D. 1999, ApJ, 523, L1 [NASA ADS] [CrossRef] [Google Scholar]
 Shaw, J. R., & Chluba, J. 2011, MNRAS, 415, 1343 [NASA ADS] [CrossRef] [Google Scholar]
 Sheth, R. K., & Tormen, G. 1999, MNRAS, 308, 119 [NASA ADS] [CrossRef] [Google Scholar]
 Shull, J. M., & van Steenberg, M. E. 1985, ApJ, 298, 268 [NASA ADS] [CrossRef] [Google Scholar]
 Sjöstrand, T., Mrenna, S., & Skands, P. 2008, Comp. Phys. Commun., 178, 852 [Google Scholar]
 Slatyer, T. R. 2010, J. Cosmology Astropart. Phys., 2, 28 [NASA ADS] [CrossRef] [Google Scholar]
 Slatyer, T. R., Padmanabhan, N., & Finkbeiner, D. P. 2009, Phys. Rev. D, 80, 043526 [NASA ADS] [CrossRef] [Google Scholar]
 Sunyaev, R. A., & Zeldovich, Y. B. 1970, Ap&SS, 7, 3 [NASA ADS] [Google Scholar]
 Valdés, M., Ferrara, A., Mapelli, M., & Ripamonti, E. 2007, MNRAS, 377, 245 [NASA ADS] [CrossRef] [Google Scholar]
 Vincent, A. C., Xue, W., & Cline, J. M. 2010, Phys. Rev. D, 82, 123519 [NASA ADS] [CrossRef] [Google Scholar]
 Wommer, E., Melia, F., & Fatuzzo, M. 2008, MNRAS, 387, 987 [NASA ADS] [CrossRef] [Google Scholar]
 XENON100 Collaboration, Aprile, E., Arisaka, K., Arneodo, F., et al. 2011, Phys. Rev. D, 84, 052003 [NASA ADS] [CrossRef] [Google Scholar]
 Zaldarriaga, M., Colombo, L., Komatsu, E., et al. 2008 [arXiv:0811.3918] [Google Scholar]
 Zavala, J., Vogelsberger, M., Slatyer, T. R., Loeb, A., & Springel, V. 2011, Phys. Rev. D, 83, 123513 [NASA ADS] [CrossRef] [Google Scholar]
 Zdziarski, A. A., & Svensson, R. 1989, ApJ, 344, 551 [NASA ADS] [CrossRef] [Google Scholar]
 Zhang, L., Chen, X., Lei, Y.A., & Si, Z.G. 2006, Phys. Rev. D, 74, 103519 [NASA ADS] [CrossRef] [Google Scholar]
All Tables
All Figures
Fig. 1 The redshift z′ where the optical depth for photons reaches unity (i.e., τ(z,z′) = 1) for several “observer’s redshifts”: z = 0,10,500,1000. Here the energy plotted is the photon energy at redshift z. The light gray region corresponds to the DM mass interval considered in this paper. 

In the text 
Fig. 2 Lefthand column from top to bottom: (i) fparameters for the μchannel assuming m_{DM} = 10 GeV (dark solid lines) and m_{DM} = 100 GeV (light solid lines). In both cases, the lowest curve out of the triple of lines corresponds to only the smooth background, the middle one includes halos with a lower mass cutoff of 10^{6} M_{⊙}, while the top one has 10^{9} M_{⊙}. The dotted lines represent highredshift fits (valid for z ≳ 170) as given by Slatyer et al. (2009). The dashed line shows the CMB visibility function. Here zg(z)dlnz gives the probability that the CMB photon last scattered in the logarithmic redshift interval dlnz centered on redshift z. (ii) Fraction of free electrons as a function of redshift for m_{DM} = 10 GeV and 100 GeV, along with two values for ⟨ σ_{A}υ ⟩ /( ⟨ σ_{A}υ ⟩ _{std}m_{DM} [ GeV] ): 1 and 10. Here ⟨ σ_{A}υ ⟩ _{std} = 3 × 10^{26} cm^{3} s^{1} is the standard thermal cross section. For the meaning of each line, see the description given in the legend. The lowest dotted line represents the standard ΛCDM case with no additional energy input from the annihilating DM. (iii) Matter temperature as a function of redshift. The models shown are exactly the same as in the panel above. The long dashed line compares the temperature of CMB. (iv) The CMB visibility function g(z) times redshift z for the same models as shown in the above two panels. Righthand column from top to bottom: (i) angular power spectra of CMB temperature fluctuations for the same models as already given above. (ii) Temperature and Emode polarization crossspectra. (iii) Emode polarization spectra. In all of the righthand panels the points with errorbars show the 7year measurements by the WMAP space mission. 

In the text 
Fig. 3 The WMAP7 parameter constraints for the μ annihilation channel. Along with two annihilation parameters log _{10}(m_{DM}) and ⟨ σ_{A}υ ⟩ /m_{DM}, we also show the constraints on the other 6parameter ΛCDM background model: Ω_{B}h^{2}, Ω_{DM}h^{2}, H_{0}, τ, n_{S}, log _{10}(A_{S}). From the inside out the colored 2D areas show the 1sigma and 2sigma regions after marginalization over the other parameters, and the dashed lines show the same regions for the basic 6parameter ΛCDM without annihilating DM. The topmost panels in each column plot the marginalized 1D probability distributions for all of the parameters. 

In the text 
Fig. 4 Upper panel: fparameters for all of the channels considered in this work and for several DM particle masses in the interval m_{DM} = 5 − 100 GeV. For all of the cases, the asymptotic highredshift fparameter values have been renormalized to be equal to one. Lower panel: the same as above, after dividing with a “typical shape”, (f_{max}(z) + f_{min}(z))/2, of the fparameter curve. 

In the text 
Fig. 5 WMAP, Planck, and CVL 1sigma constraints on the ⟨ σ_{A}υ ⟩ − m_{DM} plane for μ (upper panel) and e (lower panel) annihilation channels. The solid lines show the upper bounds on annihilation cross section as determined directly through full MCMC calculations. The shaded regions around solid lines show the results from the simple recipe of Eq. (1) with values of r taken from Table 1. The vertical gray stripe shows the range of WIMP masses (m_{DM} = 6 − 8 GeV) that provide a good fit to CoGeNT and DAMA/LIBRA data (Hooper et al. 2010). The vertical dotted line marks the lowest DM particle mass 5 GeV available for PYTHIA simulations. Extrapolations are shown below that value. 

In the text 
Current usage metrics show cumulative count of Article Views (fulltext article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 4896 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.