Anisotropies of the cosmic microwave background (CMB) have proven to be a very powerful tool to constrain dark matter annihilation at the epoch of recombination. However, CMB constraints are currently derived using a number of reasonable but yet un-tested assumptions that could potentially lead to a misestimation of the true bounds. In this paper we examine the potential impact of these systematic effects. In particular, we separately study the propagation of the secondary particles produced by annihilation in two energy regimes; first following the shower from the initial particle energy to the keV scale, and then tracking the resulting secondary particles from this scale to the absorption of their energy as heat, ionization, or excitation of the medium. We improve both the high and low energy parts of the calculation, in particular finding that our more accurate treatment of losses to sub-10.2 eV photons produced by scattering of high-energy electrons weakens the constraints on particular DM annihilation models by up to a factor of two. On the other hand, we find that the uncertainties we examine for the low energy propagation do not significantly affect the results for current and upcoming CMB data. We include the evaluation of the precise amount of excitation energy, in the form of Lyman-alpha photons, produced by the propagation of the shower, and examine the effects of varying the Helium fraction and Helium ionization fraction. In the recent literature, simple approximations for the fraction of energy absorbed in different channels have often been used to derive CMB constraints: we assess the impact of using accurate versus approximate energy fractions. Finally we check that the choice of recombination code (between RECFAST v1.5 and COSMOREC), to calculate the evolution of the free electron fraction in the presence of dark matter annihilation, introduces negligible differences.

We study the development of high energy (E_in < 1TeV) cascades produced by a primary electron of energy E_in injected into the intergalactic medium (IGM). To this aim we have developed the new code MEDEA (Monte Carlo Energy DEposition Analysis) which includes Bremsstrahlung and Inverse Compton (IC) processes, along with H/He collisional ionizations and excitations, and electron-electron collisions. The cascade energy partition into heating, excitations and ionizations depends primarily on the IGM ionized fraction, x_e, but also on redshift, z, due to IC on CMB photons. While Bremsstrahlung is unimportant under most conditions, IC becomes largely dominant at energies E_in > 1MeV. The main effect of IC at injection energies E_in < 100MeV is a significant boost of the fraction of energy converted into low energy photons (h\nu < 10.2eV) which do not further interact with the IGM. For energies E_in > 1GeV CMB photons are preferentially upscattered within the X-ray spectrum ($h\nu > 10^4eV) and can free stream to the observer. Complete tables of the fractional energy depositions as a function of redshift, E_in and ionized fraction are given. Our results can be used in many astrophysical contexts, with an obvious application related to the study of decaying/annihilating Dark Matter (DM) candidates in the high-z Universe.

Beyond reionization epoch cosmic hydrogen is neutral and can be directly observed through its 21 cm line signal. If dark matter (DM) decays or annihilates the corresponding energy input affects the hydrogen kinetic temperature and ionized fraction, and contributes to the Ly_alpha background. The changes induced by these processes on the 21 cm signal can then be used to constrain the proposed DM candidates, among which we select the three most popular ones: (i) 25-keV decaying sterile neutrinos, (ii) 10-MeV decaying light dark matter (LDM) and (iii) 10-MeV annihilating LDM. Although we find that the DM effects are considerably smaller than found by previous studies (due to a more physical description of the energy transfer from DM to the gas), we conclude that combined observations of the 21 cm background and of its gradient should be able to put constrains at least on LDM candidates. In fact, LDM decays (annihilations) induce differential brightness temperature variations with respect to the non decaying/annihilating DM case up to Delta_delta T_b=8 (22) mK at about 50 (15) MHz. In principle this signal could be detected both by current single dish radio telescopes and future facilities as LOFAR; however, this assumes that ionospheric, interference and foreground issues can be properly taken care of.

We use numerical simulations of cosmic reionization and radiative processes related to the HI 21 cm emission line to produce synthetic radio maps as seen by next generation radio telescopes that will operate at low radio frequencies (e.g. LOFAR). Two different scenarios, in which the end of reionization occurs early (z approx 13) or late (z approx 8) depending on the Initial Mass Function (IMF) of the first stars and ionizing photon escape fraction, have been explored. For each of these models we produce synthetic HI 21 cm emission maps by convolving the simulation outputs with the provisional LOFAR sampling function in the frequency range 76-140 MHz. If reionization occurs late, LOFAR will be able to detect individual HI structures on arcmin scales, emitting at a brightness temperature of approx 35 mK as a 3-sigma signal in about 1000 hours of observing time. In the case of early reionization, the detection would be unlikely, due to decreased sensitivity and increased sky temperatures. These results assume that ionospheric, interference and foreground issues are fully under control.