The effects of the smallscale behaviour of dark matter power spectrum on CMB spectral distortion
Abstract
After numerous astronomical and experimental searches, the precise particle nature of dark matter is still unknown. The standard Weakly Interacting Massive Particle(WIMP) dark matter, despite successfully explaining the largescale features of the universe, has longstanding smallscale issues. The spectral distortion in the Cosmic Microwave Background(CMB) caused by Silk damping in the prerecombination era allows one to access information on a range of small scales , whose dynamics can be precisely described using linear theory. In this paper, we investigate the possibility of using the Silk damping induced CMB spectral distortion as a probe of the smallscale power. We consider four suggested alternative dark matter candidates—Warm Dark Matter (WDM), Late Forming Dark Matter (LFDM), Ultra Light Axion (ULA) dark matter and Charged Decaying Dark Matter (CHDM); the matter power in all these models deviate significantly from the CDM model at small scales. We compute the spectral distortion of CMB for these alternative models and compare our results with the CDM model. We show that the main impact of alternative models is to alter the subhorizon evolution of the Newtonian potential which affects the latetime behaviour of spectral distortion of CMB. The parameter diminishes by a few percent as compared to the CDM model for a range of parameters of these models: LFDM for formation redshift (7%); WDM for mass (2%); CHDM for decay redshift (5%); ULA for mass (3%). This effect from the prerecombination era can be masked by orders of magnitude higher distortions generated by latetime sources, e.g. the Epoch of Reionization and tSZ from the cluster of galaxies. We also briefly discuss the detectability of this deviation in light of the upcoming CMB experiment PIXIE, which might have the sensitivity to detect this signal from the prerecombination phase.
1,2]Abir Sarkar, 1]Shiv.K.Sethi, 3]Subinoy Das
Prepared for submission to JCAP
The effects of the smallscale behaviour of dark matter power spectrum on CMB spectral distortion

Raman Research Institute, CV Raman Ave Sadashivnagar,
Bengaluru, Karnataka 560080, India 
Indian Institute of Science,CV Raman Ave, Devasandra Layout,
Bengaluru, Karnataka 560012, India 
Indian Institute of Astrophysics,100 Feet Rd, Madiwala, 2nd Block,
Koramangala, Bengaluru, Karnataka 560034, India
Contents
1 Introduction
After intensive searches throughout the world over decades, the nature of dark matter is yet to be confirmed. The existence of this component, whose presence is revealed only through its gravitational interaction, is well established by many observations covering a broad range of length scales and epochs of the universe. This list contains cosmic microwave background (CMB) anisotropy experiments [1, 2, 3], large scale structure surveys [4, 5, 6], the study of the galaxy rotation curves [7], cosmological weak gravitational lensing observations [8, 9], etc.
One of the leading candidates for dark matter, the Weakly Interacting Massive Particle (WIMP) or the traditional cold dark matter (CDM), is inspired by the wellknown WIMP miracle[10]. The supersymmetric extension of the standard model of particle physics gives rise to a particle with selfannihilation crosssection and mass in the range of 100 GeV, which coincidentally produces the correct present abundance of dark matter. Inspired by this discovery, a lot of direct [11, 12, 13], indirect [14, 15, 16, 17] and collider [18, 19] searches have been performed worldwide but none of these experiments have yet succeeded in providing consistent information about the particle nature of dark matter. In fact, the results from many of these experiments are found to be in conflict [20] with each other.
Besides, there also exist some longstanding astrophysical problems with the WIMP. One of them is the cuspcore problem [21], indicated by the discrepancy between increasing dark matter halo profile (cusp) towards the centre of the galaxy from CDM Nbody simulations [22], while observationally relatively flat density profiles are found [23]. Another issue with WIMP is the missing satellite problem[24, 25]; Nbody simulations of structure formation with CDM produce much more satellite halos of a MilkyWay type galaxy than observed. Another issue is the “too big to fail” [26, 27] problem, which underlines the fact, based on Nbody simulations, that a majority of the most massive subhalos of the Milky Way are too dense to host any of its bright satellites. Some recent works claim that even when the effects of small scale baryonic physics are included, these issues may still persist [28, 29, 30]. All these problems have inspired a drive to go beyond the standard picture of dark matter and consider alternative candidates, which differ from CDM on galactic scales but must reproduce its success on cosmological scales.
We consider four alternative dark matter models in this work: Warm Dark Matter (WDM), Late Forming Dark Matter (LFDM), Ultra Light Axion Dark Matter (ULADM), dark matter produced by the decay of a heavy charged particle in the early universe (CHDM). The best studied alternative dark matter model is the WDM model (e.g. [31, 32]); we consider only thermally produced WDM in this work. This model has been proposed to alleviate some of the smallscale problems of the CDM model [33]. WDM has been extensively tested against observations at small scales: number of satellites based on Nbody simulations [34, 35, 36]; The cuspcore problem [37, 38]; Lyman forest flux power spectrum [32]; structure formation constraints based on hydrodynamic simulations [39]; the “too big to fail” issue [36]. These detailed studies prefer WDM masses (or lower limit on the mass) to lie in the range: . These constraints do not reveal a unique picture with respect to thermal WDM being a possible dark matter candidate and underline the need to investigate alternative models. The dark matter produced by the decay of a heavy charged particle in the early universe(CHDM) has also been proposed to address issues with smallscale power of the CDM model [40]. LFDM has its origin in extended neutrino physics [41] while the ULA dark matter originates from the string axiverse scenario [42]. In these cases, the dark matter starts behaving as CDM after the epoch of Big Bang Nucleosynthesis(BBN) and before the epoch of MatterRadiation Equality (MRE) and have similar features (suppression in small scale power followed by damped oscillations) in matter power spectra. For LFDM and CHDM the suppression happens at scales that enter the horizon before the dark matter is formed. For WDM and ULADM models, the scales that suffer suppression are determined by either freestreaming of the particle (WDM) or variable sound speed of the scalar field (ULADM). Both of these scales depend on the mass of the particle/scalar field. Recent Nbody simulations based on the LFDM and ULA models suggest that these models offer acceptable solutions to the cuspcore issue while being consistent with largescale clustering data ([43] [44, 45, 46]). Cosmological constraints favour the following lower limits on the mass of ULA: [47, 48].
All the models we consider deviate from the CDM model at small scales. One possible probe of this small scale deviation is the spectral distortion of the CMB caused by Silk damping in the prerecombination era. This damping pumps entropy into the thermal plasma and can distort the CMB spectra if the energy injection occurs after [49, 50, 51, 52]. The acoustic waves on scales get dissipated before the recombination which means that the CMB spectral distortion can be used to constrain the matter power at these scales (for details e.g. [53, 54, 55, 56, 57, 58, 59, 60, 61, 62, 63, 64, 65, 66, 67, 68]). The COBEFIRAS results gave the current upper bounds on the spectral distortion parameters: and [69]. Silk damping generally results in spectral distortion amplitudes many orders of magnitude smaller than these limits. The upcoming experiment PIXIE [70] will be able to measure and , assuming zero foreground contamination. The foregrounds depend on several physical parameters and marginalisation over all of them will degrade the sensitivity[71].
In this work, we study CMB distortion parameters for alternative dark matter models that differ from the CDM model at small scales as possible probes of these models. For these models, we compute the distortion parameters using the tightcoupling approximation in the prerecombination era.
The paper is organised as follows. In Section 2, we provide a brief description of all the models considered in this work. In section 3, we review the physical picture of the creation and evolution of spectral distortion in the CMB and describe the relation between the dissipation of acoustic wave and CMB spectral distortion. In section 4, we discuss in detail the physical basis of how a change in the dark matter model impacts Silk damping in the tightcoupling approximation. Section 5 contains the main results of this work. Section 6 is reserved for the conclusion and prospects of this paper.
Throughout this paper, unless specified otherwise, we use the bestfit cosmological parameters given by Planck [1]: spatially flat universe with , (both at the present epoch), and .
2 Models of dark matter
In this section, we describe the origins of the models considered in this work by motivating their particle physics origins and briefly discuss their cosmological signatures.
2.1 Wdm
Warm dark matter (WDM) particles with a common mass of , inspired by particle physics models of sterile neutrinos, have been advocated as a solution to the smallscale anomalies of CDM. A WDM particle is mainly produced nonthermally before the epoch of BBN by active to sterile neutrino oscillation [72] in the early universe. However, there are models of thermal keV mass Sterile neutrinos which were in thermal equilibrium with some hidden sector particle and decouple with correct dark matter relic density [73].
Sterile neutrinos in this mass range cannot be detected in standard WIMP searches at least with current experimental capabilities, but in galactic Xray data, its imprint can be captured if a DM sterile neutrino decays to a photon and relativistic active neutrino. This photon would have keV energy which in principal can be detected as a galactic Xray excess. Recent Xray anomalies from XMMNewton and Chandra data can be explained by the decay of a 7 keV WDM[74].
In this work, we have considered only thermally produced WDM. In the early universe, WDM is highly relativistic, and its freestreaming scale is the horizon at that time. As time proceeds, it cools down and passing through its semirelativistic phase it finally becomes nonrelativistic CDM. The freestreaming scale during its semirelativistic phase is given by [33]
(2.1) 
Here is the relic density of WDM and is the mass. In the matter power spectrum produced by WDM, the power is cut at the WDM freestreaming scale as compared to the standard CDM model and the decrement is smooth in nature owing to the fact that the thermal velocity of WDM decays as and is not small enough to be clustered at small scales for a long time (Figure 1). Clearly, lighter WDM will cut power at larger scales (e.g. [33]). This freestreaming also suppresses the formation of lowmass halos or subhalos and its finite phasespace density prevents the development of density cusps[32, 75]. In this work we consider WDM particles in the mass range: .
2.2 Lfdm
In this scenario, the dark matter is formed as a consequence of phase transition of a scalar field. The scalar field, which was initially trapped in a metastable minimum by thermal effects or due to its interaction with other scalar fields resulting in a hybrid potential in the early universe, makes a phase transition (which is generally between the redshift of BBN and MRE in this model) to the true minima. After this, the scalar field starts oscillating coherently and its equation of state changes from to , making it clump exactly like the CDM. The dynamics of the phase transition is very similar to the case of the hybrid inflation [76], where the scalar field plays the role of the waterfall field and behaves like dark matter after the phase transition. The hybrid potential for two scalar fields in the context of LFDM was derived in [41, 77] and is given by
(2.2) 
For a high value of in the early universe, the scalar field is trapped in a metastable minimum behaving like the cosmological constant. Later, at a critical redshift , becomes Tachyonic and rolls down to the true minima. After reaching the true minima, exhibits coherent oscillations and behaves exactly like the CDM. For details of this model, we refer the reader to [41].
LFDM can also originate from a Fermionic field which has been proposed in [78]. Here, a relativistic Fermionic fluid like the massless neutrino stops freestreaming at a certain redshift due to the onset of a strong fifth force mediated by a scalar field. As the fermions are massive, eventually the fifth force binds all the fermions within a Compton volume of the mediating scalar, forming nuggets. These nuggets are much heavier than the original fermions, and since their formation, they behave exactly like the CDM.
For both of the scenarios mentioned above, LFDM gets its initial conditions for evolution from the massless neutrino. If LFDM is formed at , power is suppressed at a scale that entered the horizon at that redshift. Unlike WDM, the suppression is sharp as the dark matter is assumed to form via an almost instantaneous phase transition (Figure 1). At , the matter power spectrum carries damped oscillations, a typical feature exhibited by the massless neutrino. This feature is in contrast to WDM where the power falls monotonically without any oscillation. Earlier the dark matter forms, smaller is the scale where the decrement of power occurs. The cosmology of this model is studied in [79] and it is found that dark matter should form deep inside the radiation dominated era and before . In the present work, we consider LFDM having formation redshift in the range .
2.3 Ula
Another very well studied dark matter is the dark matter from ultra light axion (ULA) fields in the context of string axiverse [42]. An axionlike particle as a dark matter candidate can be described by [46, 80, 81, 82] a twoparameter model, whose action is given by:
(2.3) 
where is a dimensionless and periodic scalar field, represented as . and are the two parameters
of the model. For sufficient small value of (which is the case for dark matter), it can be shown that mass of the scalar is given by . For a cosmologically and astrophysically acceptable dark matter candidate, a reasonable value for the mass is eV. All models of particle physics derived from string theory have several periodic scalar fields such as and it has been argued that such a low mass is quite reasonable from particle physics perspective.
ULA obtains its initial conditions after spontaneous symmetry breaking in the early universe and behaves like a coherent scalar field. Early when is high, the friction term dominates and the field is stuck at some random initial value and behaves like the cosmological constant. Later, when at a certain redshift, the field rolls off and start oscillating coherently around the nearest minima of the periodic potential and starts behaving like CDM. The adiabatic perturbations in the scalar field have a momentumdependent and thus massdependent effective sound speed. At scales below the effective sound horizon, perturbations are washed out due to freestreaming, and the matter power spectrum features very similar to LFDM (Figure 1) are found in the matter power spectrum. The freestreaming scale is given by [46]
(2.4) 
This means that lighter axions will push the scale of suppression to a higher value. The cosmologically relevant mass of ULA ranges from and can be as tiny as eV. Some recent works have put a lower limit on the mass of ULA: [83, 48]. In this paper, we consider the mass range: .
2.4 Charged decaying dark matter
This model, its variants, and their cosmological implications have been investigated in detail [40, 84, 85, 86, 87, 88]. We consider a model in which a heavy negatively charged particle of mass decays into a heavy neutral particle of mass and a relativistic electron (supersymmetric models in which a selectron decays into an electron and a gravitino might achieve this scenario [40]). These two masses and the decay time parameterize the model. The decay time and the mass difference between the two heavy particles are tightly constrained because the relativistic electron thermalizes with electronphoton coupled system, thereby causing spectral distortion of CMB if the decay time corresponds to redshifts . The threemomentum of the relativistic electron and in the limit all this relativistic energy is transferred to the photon gas, we get the fractional energy increase:
(2.5) 
Using the current bounds on and parameters, we get for decay times in the redshift range .
This constrains the energy density of the relativistic electrons to be dynamically unimportant and allows us to assume that the masses of the two heavy particles are the same. This means that the main difference between models such as the LFDM model and the decay charged particle model is that whereas the initial conditions (density and velocity perturbations) in the former case arise from massless neutrino, they are inherited from the baryonphoton fluid in the latter case. ^{1}^{1}1This assumption needs further explanation. For a massive charged particle to share the bulk velocity of the baryonphoton fluid, it must be tightly coupled to this fluid which means the time scale of energy exchange between this particle and the photonbaryons fluid should be far smaller than the expansion rate. This time scale which is much shorter than the expansion rate at redshifts of interest. Another complication, in this case, arises from the possibility that the charged dark matter particle can form an atom with a proton. This hydrogenic atom would have binding energy on the order of 24 keV (). These atoms, in turn, could be converted into an ion containing the charged dark matter particle and doubly ionized helium through charge exchange reactions(for details e.g. [89]). Only a tiny fraction of these atoms will form until the epoch there are not enough number of energetic CMB photons on the tail of black body spectrum to ionize the atom; this epoch corresponds roughly to a temperature . This could be an additional source of spectral distortion which we neglect in this paper. This results in a qualitative difference between the two cases as seen in Figure 1. While density and velocity perturbations of massless neutrinos decay after horizon entry, these perturbations oscillate with nearly constant amplitude for the photonbaryon fluid after horizon entry for .
It is seen in Figure 1 that the matter power spectrum in this case oscillates for scales that are subhorizon during the predecay phase but its value can exceed the matter power spectrum for the CDM model for the same cosmological parameters.
This can be understood as follows. In the CDM model, the density perturbations of the CDM component, in Synchronous gauge, or they are completely determined by metric perturbations and are independent of velocity perturbations which are zero at all times, . In this case, the CDM density perturbations at subhorizon scales grow logarithmically during the radiationdominated era. However, in the decaying charged particle model, the initial velocity perturbations of the CDM component (the postdecay neutral particle) are derived from the photonbaryon fluid and constitute an additional source of density perturbations. For wavenumbers at which the velocity perturbations of the initial conditions combine in phase with density perturbations of the CDM component, the density perturbations can overshoot perturbations in the CDM model. On the other hand, this effect also serves to increase the decrement as compared to the CDM model for wavenumbers at which velocity perturbations act to suppress the growth of density perturbations.
All the models we consider here leave the matterradiation equality unchanged and their matter energy densities at the present epoch are normalized to . In Figure 1, the power spectra of all the four models are displayed and compared to that of CDM. In the Figure, we have plotted models which induce cut in power at nearly the same scale . The power spectra shown in the figure have been computed using the publiclyavailable codes CMBFAST (WDM, CHDM, and LFDM by modifying this code) and axionCAMB (ULA).
2.5 Evolution of the transfer functions of different dark matter candidates
In this subsection, we briefly describe how the dark matter transfer function evolves for different dark matter candidates considered in this work. The transfer functions are given in Figure 2 at . The model specifications are identical to that used in Figure 1. This allows us to motivate the
discussion in section 4.
WDM: In this work, we have considered keVmass thermally produced WDM candidates, which become nonrelativistic at –. So in the very early universe, when they are still relativistic, the transfer function at subhorizon scales will be similar to that of a relativistic particle like the massless neutrino. As time progresses they become nonrelativistic but still cause suppression as compared to the CDM model at small scales owing to their freestreaming velocity that decays at a slower pace, as . We have shown the evolution of transfer functions for WDM with of mass 0.33 keV in Figure 2, which is relativistic until , so the subhorizon features at and
are essentially the same. WDM with higher masses become nonrelativistic earlier, before the onset of the distortion era (discussed in the next section).
ULADM: The ULADM forms due to spontaneous symmetry breaking in the early universe and behaves like a coherent scalar field. Once the mass of the field drops below the Hubble constant, the field starts oscillating coherently around the true minima of the potential and behaves like cold dark matter. We have plotted the evolution of transfer function for . ULADM with this mass decouples from the Hubble drag at and the suppression occurs at scales smaller than . So both at the onset of distortion (Figure 3 in section 3) and at earlier times the subhorizon transfer functions are the same, as seen from the Figure 2.
LFDM: As mentioned in subsection 2.2, LFDM is formed due to phase transition in the neutrino sector and gets its initial conditions from the massless neutrino. In Figure 2, we show
the evolution of the transfer function for . We notice that the transfer functions of LFDM and the WDM with mass 0.33 keV are identical at and because at high redshifts both behave like massless neutrinos.
At , the WDM has already become nonrelativistic but LFDM continues to be relativistic which explains greater suppression in the LFDM power spectrum. After , the LFDM behaves as CDM, but WDM continues to cause suppression
of power at small scales owing to its slowly decaying but still significant freestreaming velocities. At more recent time, i.e., , the LFDM transfer function is identical to that of CDM at scales , the scale that enters the horizon at the formation redshift, but smaller scales still carry the massless neutrinolike features.
CHDM: This model is qualitatively different from the other models as the dark matter candidate, in this case, is charged and coupled to the photonbaryon plasma at early times. The time evolution of CHDM transfer function is identical to that of baryons at high redshifts with as explained in subsection 2.4. At lower redshifts and at scales larger than the scale that entered the horizon at , its transfer function is identical to CDM. At smaller scales, the transfer function can exceed the matter power of the CDM model, as seen in Figure 2 and explained in subsection 2.4.
3 Physics of Spectral Distortion
In this section, the physics of spectral distortion in the CMB will be discussed in more detail. In Figure 3 different epochs related to the evolution of CMB distortion are shown. There are four phases in this regard, and we provide a brief physical picture of each of these.
In the early universe, there are two classes of physical processes that can cause energy exchange between electrons and photons and therefore act to equilibrate photon distribution function to a black body. One is photon conserving process, dominated by inverse Compton scattering and photon producing processes which are mediated by double Compton scattering and freefree emission with double Compton scattering dominating the rate of production of photons in the early universe. If both these processes act on time scales shorter than the expansion time scale of the universe, any injection of energy in the universe is rapidly shared between the photons and charged particles, and the CMB spectrum relaxes to a black body ( see e.g. [51, 90, 91, 92]).
This condition is obtained in the first era . Any injection of energy during this era, e.g. annihilation and BBN, leaves no trace on the CMB spectrum.
For the rate of photon production is not large enough to equilibriate the photon distribution function to a black body except at low photon frequencies. This redshift thus approximately marks the time after which the signature of energy injection can no longer be erased. As Compton scattering remains an efficient process for energy exchange between photons and electrons, the resultant distribution function relaxes to a BoseEinstein distribution with but with a nonzero chemical potential, . This arises because Compton scattering is a photon conserving process that cannot cause a transition between two Black bodies at two different temperatures. During this phase, the evolution of the photon distribution function, with Compton scattering being the only energyexchange process, is given by the wellknown Kompaneets equation:
(3.1) 
Here is the electron temperature, is the number density of electrons, is the mass of the electron and , being the frequency of the photon, is the dimensionless frequency. The general equilibrium solution of this equation is given by [93, 49]
(3.2) 
The chemical potential can be related to the fractional energy absorbed by the photon gas: (e.g. [51, 49]).
For , the Comptonization process is not efficient enough to bring the photon distribution function into equilibrium, and the equilibrium solution of Kompaneets equation (Eq. 3.2) is not valid anymore. This redshift marks the end of distortion era. This era is followed by “distortion” and “distortion” eras. We first discuss the physics of distortion.
The distribution function in Eq. 3.2 satisfies the following identity : where . Assuming small departure from equilibrium, Eq. 3.1 can be written as follows:
(3.3) 
The time variable of Eq. 3.3 can be modified as a new parameter which denotes the temperature difference between electrons and photons:
\[ https://www.sharelatex.com/project/5666c560182bc84705fe62a3dy=\frac{n_{e}% \sigma_{T}k_{B}(T_{e}T)}{m_{e}c}dt \]  (3.4) 
Eq. 3.3 can be solved to give [93, 94]:
(3.5) 
In this case the fractional change in photon energy can be related to the parameter as: [94].
Between  and distortion era there exists another era called the “intermediate (type) distortion era”. In this time, the Comptonization time scale is not short enough to relax the spectrum to equilibrium. Instead, the system settles into a state where the distortion is given by the sum of , and some residual distortion [65, 95]. The thermalization in this era is approximated as a weighted combination of pure  and distortion[95]. The residual distortion (type) is between 1030 of the total distortion. This type distortion depends sensitively on the time of energy injection, which is not the case for pure  or distortion. A different approach was taken in [60] to quantitatively describe the distortion. The distortion can be characterised by a modified parameter defined as:
(3.6) 
The Kompaneets equation is written in terms of this new time parameter and then expanded about . The solution shows that the distorted spectrum and thus the distortion is dependent on . According to Eq. 3.6, is sensitive to the time when the energy was injected into the system through . This makes the distortion able to estimate not only the amount of energy injected but also the time of injection, something that can’t be estimated by observing the or type distortion.
After , becomes very small and distortion epoch commences and lasts up to present time. Late time phenomena of the universe like reionization [96, 60], heating of photons by warmhot intergalactic medium [97, 98, 99] and SZ effects from groups and clusters of galaxies ([64, 55] and references therein) also contribute to the distortion. As long as the amount of distortion is small, the three different distortions (++) can be linearly added to give the final spectrum.
Various processes like decay of massive particles[65, 100, 101], annihilation of particles [58, 102, 103], dissipation of acoustic wave [58, 50, 54, 56, 57, 104], adiabatic cooling of electrons [58, 59] and cosmological recombination radiation [105, 62, 63, 106, 107] can contribute energy into or extract energy from the photonbaryon plasma leading to spectral distortion in the early universe. In this work, we study the spectral distortion caused by dissipation of acoustic wave (Silk Damping) for different dark matter models. A detailed description of that process is given in next section.
3.1 Energy released due to Silk Damping
At , the photon and baryons are tightly coupled to each other via Compton scattering and behave like a single fluid. Adiabatic perturbations in this fluid behave like standing waves inside the sound horizon , where is the speed of sound in the plasma and . At scales much smaller than the sound horizon, photon diffusion causes damping of density perturbation and bulk motion of this fluid [108]. This process (Silk damping) can be modelled by expanding the evolution of perturbations in the coupled photonbaryon fluid to second order in the mean free path of photons [109, 110] or as the dissipation of energy of sounds waves owing to radiative viscosity and also on thermal conduction at late times [111]. Silk damping causes the injection of entropy into the thermal plasma.
As this is a continuous energy injection process, it can lead to type, type or type distortion depending on the era when the energy was injected. In the redshift range, , the structures corresponding to scales are completely wiped out due to this damping. Thus, conversely, one possible way to study and constrain the initial power spectrum at these scales is by observing the spectral distortion that is imprinted in the CMB spectrum due to this process. It should be noted that this is the only known probe of linear structures for such a wide range of small scales. In comparison, the observed CMB temperature anisotropies from Planck probe scales [1] and the smallest scales probed by galaxy clustering data correspond to nearly linear scales at the present: (e.g. [4]).
The damping of adiabatic perturbations at small scales has been well studied in the literature. Recently a precise calculation of  and type distortion due to Silk damping has been performed using the second order perturbation theory[57]. It has been shown that tightcoupling approximation provides a good approximation for modelling the Silk damping in the prerecombination era. In this approximation, the source function for heating due to Silk damping is given by [57]^{2}^{2}2the source function given in Eq. (3.7) can be understood as arising from radiative viscosity. In this case the energy pumped into the thermal plasma is proportion to the square of the product of the photon mean free path with the velocity shear field (e.g. [53]) which is proportional to . Alternatively, the dissipation of the energy can be modelled in terms of photon monopole [56]. Both these approaches give results within a factor of 3/4 of each other [57]. While numerically computing CMB spectral distortion, we use both the methods and find reasonable agreement:
(3.7) 
Where is the CMB dipole anisotropy, is the derivative of Compton scattering optical depth with respect to the conformal time and is the frequency dependent function representing the distortion defined in Eq. 3.5. The first term in the square bracket comes from heat conduction and the second one is due to radiative viscosity. At high redshift, the radiative viscosity dominates as . As any time, the average source function can be obtained by integrating over all wavenumbers:
(3.8) 
To obtain effective heating rate needs to be multiplied by and integrated over all frequencies. The integration over frequencies, , yields . This result provides the heating rate as a function of conformal time which is converted to a function of redshift by dividing it by . Thus, the final expression of heating rate is given by
(3.9) 
The and parameters are related to heating rates as follows [57, 112]
(3.10) 
and
(3.11) 
and are the redshift of decoupling and the redshift denoting the end of distortion era, respectively.
4 Spectral distortion and dark matter models
In this section, we discuss how altering the dark matter model impacts the generation of entropy owing to Silk damping in the prerecombination era. As noted in the previous section, the dynamics of this energy pumping can be captured by the time evolution of the photon dipole, .
We discuss here how the photon dipole is altered when the dark matter model is changed. We use Newtonian conformal gauge for motivating our discussion as the underlying physics is more transparent in this gauge. In the tightcoupling approximation, the photon dipole can be expressed as (e.g. [113]):
(4.1) 
Here is the photon monopole and and are the two gravitational potentials in the Newtonian gauge. If the small contribution from neutrinos is neglected, . ^{3}^{3}3 with ; for three massless neutrino degrees of freedom. Two of the models we consider here can alter . This means that the first term in Eq. (4.1), which corresponds to initial conditions, can also be used to distinguish between different dark matter models [57]. In WDM models, the dark matter particle is initially relativistic and therefore would increase . However, even during this phase the contribution of this particle to neutrino relativistic degrees of freedom is negligible as compared to the standard model neutrino; for , the particle contributes less than 1% of the energy density of a standard model neutrino. In LFDM models, a tiny fraction of neutrino energy density is used to create the dark matter particle during a phase transition (e.g. [79]). For both these models, the effect of initial conditions is generally negligible. We take into account these effects in our detailed modelling using CMBFAST but drop it in this subsection to uncover the main determinants of the impact of dark matter models on the photon dipole. and . In Eq (4.1) we have neglected and its evolution in the prefactors of the equation. ^{4}^{4}4For the charged decay model, at early times as the dark matter particle is charged and tightly coupled to the baryonphoton fluid. For this model, this is an additional effect that determines the evolution of . In the results presented in this section, we take into account this effect using CMBFAST. corresponds to the scale that undergoes Silk damping at any time:
(4.2) 
Here . Under these approximations, the evolution of can be completely determined by the Newtonian potential and the Silk damping scale . As discussed in the next section, the dissipation of energy at a given scale occurs predominantly at a time such that . As we argue below, the main effect we seek depends on the time evolution of in the time range, , i.e. from the horizon entry of a scale to the time at which the mechanical energy at this scale dissipates.
All the models we consider in this paper are based on altering the nature of CDM with an aim to suppress power at small scales. While these models impact the radiation content of the universe in the early universe (e.g. WDM model is based on a particle of mass which is relativistic in the early universe) they all leave unchanged the matterradiation equality epoch. This implies that their impact is proportional to the ratio of dark matter to the radiation energy density, , which is much smaller than unity at early times. This means that the impact of changing the dark matter change is essentially captured by the second term of Eq. (4.1), which depends on the evolution of the Newtonian potential .
In Figure 4 we show the evolution of for two of the models we consider here for a range of wavenumbers . The Newtonian potential is computed from CMBFAST code, which is written in Synchronous gauge, using the following transformation:
(4.3) 
Here and are the potentials in Synchronous Gauge (Eq. 18 of [114]). Some notable features of the comparison between the models we consider and the CDM models are: (a) the difference between the potentials is negligible of for large scales (small ). This is expected as the models we consider agree with the CDM model on large scales, (b) As the scale gets smaller, the difference between alternative DM models and the CDM model becomes more significant at larger . We discuss the nature of this deviation below.
LFDM model: In this model, the CDM forms at a redshift from a tiny fraction of relativistic neutrinos due to a phase transition. It inherits the density and velocity perturbations of the neutrino component, resulting in a sharp reduction in matter power for scales , where is the scale that enters the horizon at .
To understand Figure 4 we consider scales smaller and larger than : (a) : these scales are outside the horizon when the dark matter forms. These scales evolve outside the horizon in a purely radiation dominated universe. The perturbations at these scales are not affected by LFDM physics except in determining the initial conditions which are not necessary for the following reasons: the photon perturbations outside the horizon are constant and are set by the potential. For scales outside the horizon, the potential changes by a factor 9/10 in making a transition from radiation to the matterdominated era (e.g. [113] and references therein). If the era of LFDM lies deep inside RD era, the change is negligible. This explains the large scale behaviour of potential in Figure 4.
(b) : these scales enter the horizon before the formation of the cold dark matter. As these scales evolve subhorizon in the radiation dominated era. In radiation dominated era, the potential decays as for if perturbations in radiation determine the evolution of the potential (e.g. [113]). At subhorizon scales, neutrino perturbation decay exponentially while photon density perturbations equal baryon perturbations which oscillate with nearly constant amplitude for . The CDM perturbations, on the other hand, grow logarithmically during this phase. This means matter perturbations can determine the evolution of the potential even in deep radiation dominated era. This behaviour is seen the evolution of the potential for in Figure 4 as flattening of the potential for larger . For LFDM models, the potential decay is sharper because of the CDM forms late. This behaviour along with other features for seen in Figure 4 explains how the evolution of in LFDM models differs from the CDM model. It is important to note that this difference is most pronounced for scales that enter the horizon at . For even smaller scales, the effect of matter perturbations is less important as these scales enter the horizon when is lower and therefore matter perturbations have a smaller impact on their evolution. Equivalently, it is seen from Eq. (4.1) that the photon dipole depends on the integral of the potential and as the potential decays inside the horizon, most of the contribution comes from epochs close to the horizon entry. It is also seen in the evolution of the potential for . The potential deviates insignificantly from the CDM model for but shows sharp deviation at later times. The evolution of perturbations at these scales (from to ) correspond roughly to the case of no dark matter at all time. The potential for this case falls exponentially for large which occurs because the only source of potential in the LFDM model is photon and neutrino perturbations which both decay exponentially for . It is of interest to note that the potential at is dominated by matter perturbations for .
For the LFDM model shown in Figure 4, the scales of interest for measuring the deviation from the CDM model lie in the range . These scales enter the horizon before, but not significantly before, the era of matter formation and the perturbations at these scales suffer Silk damping after .
For models such as WDM and ULA, the preceding discussion is directly applicable. We briefly discuss the charged decaying particle case below.
Charged Decaying particle: For charged decay model, the variation in the dynamics of closely follows the previous discussion for LFDM model: nearly follows the potential for CDM model at large scales. At very small scales the effect is suppressed because the matter density is small when they enter the horizon (Figure 4). The main impact is captured by intermediate scales that enter the horizon around but before the time of decay.
Some of the salient differences between the two cases can be seen in Figure 4. For the Charged decaying particle model, the potential can exceed the CDM values for a range of wavenumbers, e.g. in the figure.
It follows from our discussion that the impact of alternative DM models on declines if the transition (decay) redshift moves deeper into the radiation dominated era. We note that the tightcoupling approximation while capturing the essential physics, tends to break down close to the recombination era; this is relevant for the computation of distortion [115, 57]. However, this approximation allows us to compare our results in two different gauges and identify the main determinants of altering CMB distortion for alternative dark matter models.
In Figure 5 we show the dynamics of for a range of alternative DM models. We compute directly from CMBFAST code but also verify that Eq. (4.1), which shows that depends on the history of the variation of , explains the main difference expected for alternative DM models. We establish it by running a smaller set of equations for density and velocity perturbations in photonbaryon fluid and dark matter along with the evolution of in conformal Newtonian gauge [114].
In Figure 6, we show the values of at of different models at the same redshift. All the models shown correspond roughly to cases where the decrement in the power spectrum occurs at .
The preceding discussion allows us to establish that the main impact of alternative dark matter models is to alter distortion and not  or distortion. The distortion era occurs in the redshift range . In this era the scales that dissipate their energy and therefore cause the CMB distortion lie in the range . In the range of redshifts of interest, the matter power for alternative models at these scales is highly suppressed as compared to the CDM model, as seen in Figure 2. However, as argued above, this does not have a significant impact on because it depends on the history of the evolution of the potential (Eq. (4.1)). As the potential decays after horizon entry (Figure 4), the photon dipole is governed by the initial condition of the potential and its evolution closer to the time of horizon entry of the scale () rather than its value at . The scale (which decays at the end of the distortion era) enters the horizon at . As the ratio of radiation to matter energy density at this redshift , the potential is insensitive to even a substantial change in the matter power at this scale at early times. The initial conditions for different models are not the same as the CDM model. However, as discussed above, requiring these models to agree with large scale observations forces the impact of the change in initial conditions on the potential to be negligible.
This also implies that the dynamics of scales that determine distortion of CMB are minimally affected by a change in the dark matter model. As we showed above, the main effect in the photon dipole is caused by scales in the range , this also means the impact of alternative dark matter on the distortion era is also negligible. We test this in our study by computing the photon dipole for scales in the range for alternative models numerically and using an analytic approximation for even smaller scales to discern the impact of initial conditions (the first term of Eq. (4.1)). For all the models we consider, the  and distortion is less than .
5 Results
Using Eqs. (4.1) and (3.7) one can analytically show that the dissipation at any given time is dominated by scales such that . Such an estimate is based on approximating by the first term of Eq. (4.1) (e.g. [57]). This gives:
(5.1) 
Here is the power spectrum of the potential at the initial time. Using (e.g. [57, 113]) where is the scalar spectral index and integrating over gives:
(5.2) 
This integral in the equation is dominated by which means the energy dissipation at any given time is governed by this condition. This result holds even when both the terms in Eq. (4.1) are retained.^{5}^{5}5While the first term of Eq. (4.1) gives a reasonable analytic estimate, the second term generally dominates the first term even for the CDM model. The factor proportional to in the second term dominates over the first term and cancels it but the remainder has the form of the first term and has a similar order of magnitude. The factor proportional to in the second term is generally subdominant. These two factors contain information of the evolution of potential which is germane to our work.
This means that the dissipation is dominated by photon dipole at at any time , . In Figure 7 we show the redshift evolution of , the difference between for different dark matter models and that of CDM. The set of model parameters we have chosen is provided in Table 1. The difference decreases as the redshift of formation for LFDM, and the decay redshift for CHDM and masses of WDM and Axion are increased. One further point to note is that the difference for alternative models increases as the redshift decreases. These results are in agreement with our discussion in the previous section which suggested that the alternative dark matter models yield substantial difference only during the distortion era.
Model  Parameter  difference of from CDM  

CDM  [1]  
LFDM  

keV  
keV  
WDM  keV  
keV  
keV  
CHDM  
eV  
eV  
ULA DM  eV  
eV 
Using Eqs. (3.7), (3.8) and (3.11), we calculate parameter for different models; is normalized using the condition on mass dispersion at scale at the present epoch: [1]. The list of parameters for different models is provided in Table 1. We note that the parameter can vary by up to for a range of models. This behaviour is in line with the variation of . The parameter approaches its value for CDM model as for LFDM and for CHDM models are increased, and also when the mass of WDM and axion is increased as expected from our discussion above. Furthermore, we can use Eq. (3.5) to determine the shape of distorted CMB spectrum.
Can the distortion be used to distinguish between different dark matter models? We address this question by comparing models where the power is cut at nearly the same scale. In Figure 8, we show the evolution of and the spectral shape of the distorted spectrum for models for which the power has been cut at . Figure 8 also shows the distorted spectrum from distortion for these cases; the spectra have the same shape (Eq. (3.5)) but differ owing to slightly different values of the parameter. A more detailed analysis based on understanding degeneracies between parameters of different models would be needed to quantify the results shown in Figure 8. We also note that this effect can be masked by several latetime phenomena like the Epoch of Reionization, thermal SunyaevZel’dovich(tSZ) effects in the galaxy groups and clusters that give rise to the spectral shape that arises from type distortion with orders of magnitude higher parameters.
Many of the models we consider are already tightly constrained by cosmological observations. From galaxy clustering and CMB anisotropy data and observed neutral hydrogen (HI) abundance at high redshifts, the LFDM models require [79, 48]. Lyman data puts even stronger constraints on such models [79]. Recent studies of ULA models suggest that eV is consistent with current observational data [47]. Other studies based on the abundance of HI at high redshift constrain the axion mass [48]. These constraints suggest that some of the models shown in Table 1 are ruled out and in particular the ones that show largest deviations from the CDM model. However, there is significant variation in the predictions of different models, and we might expect to see up to 5% deviation (e.g. ) within the framework of constraints that arise from galaxy clustering [43] and CMB anisotropies [116].
The charged decaying particle model offers a more complicated scenario of spectral distortion of the CMB. As noted in section 2.4, the relativistic electron released as the decay product rapidly thermalizes its energy with the thermal plasma. Eq. 2.5 gives the fractional increase in the photon energy density due to this process. This constrains the mass difference between the two heavy particles to be tiny. This scenario therefore presents two different ways of distinguishing this model from the CDM model. For instance, for , the decay will cause type spectral distortion in the CMB (whose amplitude will depend on ) while Silk damping causes a difference in distortion. In principle, given their distinct spectral signatures, these two phases of distortions are distinguishable from each other[100, 101].
The upcoming experiment Primordial Inflation Explorer (PIXIE) [70] is likely to improve the FIRAS bounds on CMB spectral distortion by many orders of magnitude. Its projected sensitivity corresponds to: and . However, unlike the  and distortion, the distortion created during the prerecombination era can be masked by CMB distortion in the postrecombination era. For instance, the epoch of reionization is likely to produce global CMB distortion corresponding to , though the most dominating source of postrecombination distortion is the tSZ effect in galaxy groups and clusters corresponding to [64]. These sources would constitute a strong foreground to the prerecombination distortion [70]. In principle, it might be possible to distinguish the prerecombination distortion from reionization signal using spatial information (e.g. [104, 57]) but it would be a challenge.
6 Conclusion
The largescale behaviour of dark matter has been well established by a range of cosmological observables such as CMB anisotropies and clustering of galaxies. However, there remain issues with this model at small scales, and the nature of the dominant fraction of dark matter in the universe remains a mystery.
The evolution of spectral distortion in the prerecombination era allows us to study scales in the range as density perturbations at these scale decay and leave observable signatures on CMB spectrum. These perturbations are in the linear regime of their growth in the prerecombination and therefore can be theoretically modelled very accurately.
We consider four alternative dark matter models which give significant deviations of smallscale power as compared to the CDM model and study the spectral distortion of CMB owing to Silk damping for these models. These models are motivated by different aspects of physics in the early universe: phase transition (LFDM), freestreaming of massive particles (WDM), the decay of massive charged particles(CHDM), and dynamics of a scalar field with nonzero effective mass (ULA).
We show that main impact of the models we consider is to alter the late time spectral distortion history by lowering the parameter by a few percent for an acceptable range of parameters for these models (Table 1).
In this work, we only consider models that leave invariant the matterradiation equality. This excludes models such as decaying dark matter particles that create relativistic decay products thereby delaying the matterradiation equality epoch (e.g. [117, 118]). We also excluded models that are based on extra relativistic degrees of freedom (for details see [1] and references therein) or are based on a change in initial conditions (e.g. [67, 68]).
Our analysis suggests that all class of models that give suppression of power as compared to the CDM model should result in late time spectral distortion. In section 4 we present general arguments which show that the change in spectral distortion is dominated by the evolution of potential after the scale enters the horizon. This effect scales as the ratio of energy densities of dark matter and radiation and is negligible at early times. This might also allow one to distinguish the impact of modifying the initial matter power spectrum on CMB spectral distortion from a change in the dark matter model that affects the matter power at small scales, as the former would cause all the three forms of CMB distortions while the latter would not.
Future experiments such as PIXIE and CMBS4 have the potential to establish the nature of dark matter (e.g. [116]) by unprecedented improvement in mapping the CMB spectral and spatial structures. Our work is one step in that direction and points out the challenges involved in such an endeavour.
Acknowledgements
We thank the anonymous referee for sending useful comments and suggestions that can make the paper better. The authors would like to acknowledge Rishi Khatri, Jens Chluba, David J E Marsh, Ayuki Kamada, Sunny Vagnozzi and Farinaldo Queiroz for their useful comments and suggestions.
References
 [1] Planck Collaboration, P. A. R. Ade et al., Planck 2015 results. XIII. Cosmological parameters, Astron. Astrophys. 594 (2016) A13, [arXiv:1502.01589].
 [2] WMAP Collaboration, G. Hinshaw et al., NineYear Wilkinson Microwave Anisotropy Probe (WMAP) Observations: Cosmological Parameter Results, Astrophys.J.Suppl. 208 (2013) 19, [arXiv:1212.5226].
 [3] Atacama Cosmology Telescope Collaboration, J. L. Sievers et al., The Atacama Cosmology Telescope: Cosmological parameters from three seasons of data, JCAP 1310 (2013) 060, [arXiv:1301.0824].
 [4] BOSS Collaboration, F. Beutler et al., The clustering of galaxies in the completed SDSSIII Baryon Oscillation Spectroscopic Survey: Baryon Acoustic Oscillations in Fourierspace, Mon. Not. Roy. Astron. Soc. (2016) [arXiv:1607.03149].
 [5] SDSS Collaboration, M. Tegmark et al., Cosmological Constraints from the SDSS Luminous Red Galaxies, Phys.Rev. D74 (2006) 123507, [astroph/0608632].
 [6] SDSS Collaboration, M. Tegmark et al., Cosmological parameters from SDSS and WMAP, Phys.Rev. D69 (2004) 103501, [astroph/0310723].
 [7] K. Begeman, A. Broeils, and R. Sanders, Extended rotation curves of spiral galaxies: Dark haloes and modified dynamics, Mon.Not.Roy.Astron.Soc. 249 (1991) 523.
 [8] M. Bartelmann and P. Schneider, Weak gravitational lensing, Phys.Rept. 340 (2001) 291–472, [astroph/9912508].
 [9] D. Clowe, M. Bradac, A. H. Gonzalez, M. Markevitch, S. W. Randall, C. Jones, and D. Zaritsky, A direct empirical proof of the existence of dark matter, Astrophys. J. 648 (2006) L109–L113, [astroph/0608407].
 [10] N. Craig and A. Katz, The Fraternal WIMP Miracle, JCAP 1510 (2015), no. 10 054, [arXiv:1505.07113].
 [11] G. Angloher, M. Bauer, I. Bavykina, A. Bento, C. Bucci, et al., Results from 730 kg days of the CRESSTII Dark Matter Search, Eur.Phys.J. C72 (2012) 1971, [arXiv:1109.0702].
 [12] XENON100 Collaboration Collaboration, E. Aprile et al., First Dark Matter Results from the XENON100 Experiment, Phys.Rev.Lett. 105 (2010) 131302, [arXiv:1005.0380].
 [13] CDMSII Collaboration Collaboration, Z. Ahmed et al., Results from a LowEnergy Analysis of the CDMS II Germanium Data, Phys.Rev.Lett. 106 (2011) 131302, [arXiv:1011.2482].
 [14] PAMELA Collaboration Collaboration, O. Adriani et al., PAMELA results on the cosmicray antiproton flux from 60 MeV to 180 GeV in kinetic energy, Phys.Rev.Lett. 105 (2010) 121101, [arXiv:1007.0821].
 [15] Fermi LAT Collaboration Collaboration, M. Ackermann et al., Measurement of separate cosmicray electron and positron spectra with the Fermi Large Area Telescope, Phys.Rev.Lett. 108 (2012) 011103, [arXiv:1109.0521].
 [16] HEAT Collaboration Collaboration, S. Barwick et al., Measurements of the cosmic ray positron fraction from 1GeV to 50GeV, Astrophys.J. 482 (1997) L191–L194, [astroph/9703192].
 [17] AMS01 Collaboration Collaboration, M. Aguilar et al., Cosmicray positron fraction measurement from 1 to 30GeV with AMS01, Phys.Lett. B646 (2007) 145–154, [astroph/0703154].
 [18] J. Goodman, M. Ibe, A. Rajaraman, W. Shepherd, T. M. Tait, et al., Constraints on Light Majorana dark Matter from Colliders, Phys.Lett. B695 (2011) 185–188, [arXiv:1005.1286].
 [19] P. J. Fox, R. Harnik, J. Kopp, and Y. Tsai, Missing Energy Signatures of Dark Matter at the LHC, Phys.Rev. D85 (2012) 056011, [arXiv:1109.4398].
 [20] D. Hooper, Revisiting XENON100’s Constraints (and Signals?) For LowMass Dark Matter, JCAP 1309 (2013) 035, [arXiv:1306.1790].
 [21] W. J. G. de Blok, The CoreCusp Problem, Adv. Astron. 2010 (2010) 789293, [arXiv:0910.3538].
 [22] J. F. Navarro, C. S. Frenk, and S. D. M. White, The Structure of cold dark matter halos, Astrophys. J. 462 (1996) 563–575, [astroph/9508025].
 [23] M. G. Walker and J. Penarrubia, A Method for Measuring (Slopes of) the Mass Profiles of Dwarf Spheroidal Galaxies, Astrophys. J. 742 (2011) 20, [arXiv:1108.2404].
 [24] A. A. Klypin, A. V. Kravtsov, O. Valenzuela, and F. Prada, Where are the missing Galactic satellites?, Astrophys.J. 522 (1999) 82–92, [astroph/9901240].
 [25] B. Moore, S. Ghigna, F. Governato, G. Lake, T. R. Quinn, et al., Dark matter substructure within galactic halos, Astrophys.J. 524 (1999) L19–L22, [astroph/9907411].
 [26] S. GarrisonKimmel, M. BoylanKolchin, J. S. Bullock, and E. N. Kirby, Too Big to Fail in the Local Group, arXiv:1404.5313.
 [27] M. BoylanKolchin, J. S. Bullock, and M. Kaplinghat, Too big to fail? The puzzling darkness of massive Milky Way subhaloes, Mon.Not.Roy.Astron.Soc. 415 (2011) L40, [arXiv:1103.0007].
 [28] M. S. Pawlowski, B. Famaey, D. Merritt, and P. Kroupa, On the persistence of two smallscale problems in CDM, arXiv:1510.08060.
 [29] J. Oñorbe, M. BoylanKolchin, J. S. Bullock, P. F. Hopkins, D. Kereš, C.A. FaucherGiguère, E. Quataert, and N. Murray, Forged in FIRE: cusps, cores and baryons in lowmass dwarf galaxies, Mon.Not.Roy.Astron.Soc. 454 (Dec., 2015) 2092–2106, [arXiv:1502.02036].
 [30] T. Sawala et al., Local Group galaxies emerge from the dark, arXiv:1412.2748.
 [31] A. Dolgov and S. Hansen, Massive sterile neutrinos as warm dark matter, Astropart.Phys. 16 (2002) 339–344, [hepph/0009083].
 [32] M. Viel, J. Lesgourgues, M. G. Haehnelt, S. Matarrese, and A. Riotto, Constraining warm dark matter candidates including sterile neutrinos and light gravitinos with WMAP and the Lymanalpha forest, Phys.Rev. D71 (2005) 063534, [astroph/0501562].
 [33] P. Bode, J. P. Ostriker, and N. Turok, Halo formation in warm dark matter models, Astrophys. J. 556 (2001) 93–107, [astroph/0010389].
 [34] E. Polisensky and M. Ricotti, Constraints on the Dark Matter Particle Mass from the Number of Milky Way Satellites, Phys. Rev. D83 (2011) 043506, [arXiv:1004.1459].
 [35] D. Anderhalden, J. Diemand, G. Bertone, A. V. Maccio, and A. Schneider, The Galactic Halo in Mixed Dark Matter Cosmologies, JCAP 1210 (2012) 047, [arXiv:1206.3788].
 [36] M. R. Lovell, V. Eke, C. S. Frenk, L. Gao, A. Jenkins, T. Theuns, J. Wang, D. M. White, A. Boyarsky, and O. Ruchayskiy, The Haloes of Bright Satellite Galaxies in a Warm Dark Matter Universe, Mon. Not. Roy. Astron. Soc. 420 (2012) 2318–2324, [arXiv:1104.2929].
 [37] A. V. Maccio, S. Paduroiu, D. Anderhalden, A. Schneider, and B. Moore, Cores in warm dark matter haloes: a Catch 22 problem, Mon. Not. Roy. Astron. Soc. 424 (2012) 1105–1112, [arXiv:1202.1282].
 [38] A. Schneider, R. E. Smith, A. V. Maccio, and B. Moore, Nonlinear Evolution of Cosmological Structures in Warm Dark Matter Models, Mon. Not. Roy. Astron. Soc. 424 (2012) 684, [arXiv:1112.0330].
 [39] J. Baur, N. PalanqueDelabrouille, C. Yèche, C. Magneville, and M. Viel, Lymanalpha Forests cool Warm Dark Matter, JCAP 1608 (2016), no. 08 012, [arXiv:1512.01981].
 [40] K. Sigurdson and M. Kamionkowski, Charged  particle decay and suppression of small  scale power, Phys. Rev. Lett. 92 (2004) 171302, [astroph/0311486].
 [41] S. Das and N. Weiner, Late Forming Dark Matter in Theories of Neutrino Dark Energy, Phys. Rev. D84 (2011) 123511, [astroph/0611353].
 [42] A. Arvanitaki, S. Dimopoulos, S. Dubovsky, N. Kaloper, and J. MarchRussell, String Axiverse, Phys. Rev. D81 (2010) 123530, [arXiv:0905.4720].
 [43] S. Agarwal, P. S. Corasaniti, S. Das, and Y. Rasera, Small scale clustering of late forming dark matter, arXiv:1412.1103.
 [44] H.Y. Schive, T. Chiueh, and T. Broadhurst, Cosmic Structure as the Quantum Interference of a Coherent Dark Wave, Nature Phys. 10 (2014) 496–499, [arXiv:1406.6586].
 [45] D. J. E. Marsh and A.R. Pop, Axion dark matter, solitons and the cuspcore problem, Mon. Not. Roy. Astron. Soc. 451 (2015), no. 3 2479–2492, [arXiv:1502.03456].
 [46] D. J. E. Marsh, Axion Cosmology, Phys. Rept. 643 (2016) 1–79, [arXiv:1510.07633].
 [47] L. Hui, J. P. Ostriker, S. Tremaine, and E. Witten, On the hypothesis that cosmological dark matter is composed of ultralight bosons, arXiv:1610.08297.
 [48] A. Sarkar, R. Mondal, S. Das, S. Sethi, S. Bharadwaj, and D. J. E. Marsh, The effects of the smallscale DM power on the cosmological neutral hydrogen (HI) distribution at high redshifts, JCAP 1604 (2016), no. 04 012, [arXiv:1512.03325].
 [49] R. A. Sunyaev and Y. B. Zeldovich, The interaction of matter and radiation in the hot model of the Universe, II, Astrophys. Space Sci. 7 (Apr., 1970) 20–30.
 [50] R. A. Sunyaev and Y. B. Zeldovich, Small scale entropy and adiabatic density perturbations  Antimatter in the Universe, Astrophys. Space Sci. 9 (Dec., 1970) 368–382.
 [51] A. F. Illarionov and R. A. Siuniaev, Comptonization, the backgroundradiation spectrum, and the thermal history of the universe, Exp. Theor. Phys. 18 (June, 1975) 691–699.
 [52] C. Burigana, G. de Zotti, and L. Danese, Analytical description of spectral distortions of the cosmic microwave background., Astron.Astrophys. 303 (Nov., 1995) 323.
 [53] L. Danese and G. de Zotti, The relic radiation spectrum and the thermal history of the Universe, Nuovo Cimento Rivista Serie 7 (Sept., 1977) 277–362.
 [54] R. A. Daly, Spectral distortions of the microwave background radiation resulting from the damping of pressure waves, Astrophys. J. 371 (Apr., 1991) 14–28.
 [55] C. Burigana, L. Danese, and G. de Zotti, Formation and evolution of early distortions of the microwave background spectrum  A numerical study, Astron. Astrophys. 246 (June, 1991) 49–58.
 [56] W. Hu, D. Scott, and J. Silk, Power spectrum constraints from spectral distortions in the cosmic microwave background, Astrophys. J. 430 (1994) L5–L8, [astroph/9402045].
 [57] J. Chluba, R. Khatri, and R. A. Sunyaev, CMB at 2x2 order: The dissipation of primordial acoustic waves and the observable part of the associated energy release, Mon. Not. Roy. Astron. Soc. 425 (2012) 1129–1169, [arXiv:1202.0057].
 [58] J. Chluba, E. Dimastrogiovanni, M. A. Amin, and M. Kamionkowski, Evolution of CMB spectral distortion anisotropies and tests of primordial nonGaussianity, arXiv:1610.08711.
 [59] R. Khatri, R. A. Sunyaev, and J. Chluba, Does BoseEinstein condensation of CMB photons cancel distortions created by dissipation of sound waves in the early Universe?, Astron. Astrophys. 540 (2012) A124, [arXiv:1110.0475].
 [60] R. Khatri and R. A. Sunyaev, Beyond y and : the shape of the CMB spectral distortions in the intermediate epoch, , JCAP 1209 (2012) 016, [arXiv:1207.6654].
 [61] R. Khatri and R. A. Sunyaev, Creation of the CMB spectrum: precise analytic solutions for the blackbody photosphere, JCAP 6 (June, 2012) 038, [arXiv:1203.2601].
 [62] J. A. RubinoMartin, J. Chluba, and R. A. Sunyaev, Lines in the Cosmic Microwave Background Spectrum from the Epoch of Cosmological Hydrogen Recombination, Mon. Not. Roy. Astron. Soc. 371 (2006) 1939–1952, [astroph/0607373].
 [63] J. Chluba and R. A. Sunyaev, Prerecombinational energy release and narrow features in the CMB spectrum, Astron. Astrophys. 501 (2009) 29–47, [arXiv:0803.3584].
 [64] J. C. Hill, N. Battaglia, J. Chluba, S. Ferraro, E. Schaan, and D. N. Spergel, Taking the Universe’s Temperature with Spectral Distortions of the Cosmic Microwave Background, Phys. Rev. Lett. 115 (2015), no. 26 261301, [arXiv:1507.01583].
 [65] J. Chluba and R. A. Sunyaev, The evolution of CMB spectral distortions in the early Universe, Mon. Not. Roy. Astron. Soc. 419 (2012) 1294–1314, [arXiv:1109.6552].
 [66] J. B. Dent, D. A. Easson, and H. Tashiro, Cosmological constraints from CMB distortion, Phys. Rev. D86 (2012) 023514, [arXiv:1202.6066].
 [67] J. Chluba, A. L. Erickcek, and I. BenDayan, Probing the inflaton: Smallscale power spectrum constraints from measurements of the CMB energy spectrum, Astrophys. J. 758 (2012) 76, [arXiv:1203.2681].
 [68] S. Clesse, B. Garbrecht, and Y. Zhu, Testing Inflation and Curvaton Scenarios with CMB Distortions, JCAP 1410 (2014), no. 10 046, [arXiv:1402.2257].
 [69] D. J. Fixsen, E. S. Cheng, J. M. Gales, J. C. Mather, R. A. Shafer, and E. L. Wright, The Cosmic Microwave Background spectrum from the full COBE FIRAS data set, Astrophys. J. 473 (1996) 576, [astroph/9605054].
 [70] A. Kogut et al., The Primordial Inflation Explorer (PIXIE): A Nulling Polarimeter for Cosmic Microwave Background Observations, JCAP 1107 (2011) 025, [arXiv:1105.2044].
 [71] M. H. Abitbol, J. Chluba, J. C. Hill, and B. R. Johnson, Prospects for Measuring Cosmic Microwave Background Spectral Distortions in the Presence of Foregrounds, arXiv:1705.01534.
 [72] X.D. Shi and G. M. Fuller, A New dark matter candidate: Nonthermal sterile neutrinos, Phys. Rev. Lett. 82 (1999) 2832–2835, [astroph/9810076].
 [73] S. Das and K. Sigurdson, Cosmological Limits on Hidden Sector Dark Matter, Phys. Rev. D85 (2012) 063510, [arXiv:1012.4458].
 [74] A. Boyarsky, O. Ruchayskiy, D. Iakubovskyi, and J. Franse, Unidentified Line in XRay Spectra of the Andromeda Galaxy and Perseus Galaxy Cluster, Phys. Rev. Lett. 113 (2014) 251301, [arXiv:1402.4119].
 [75] P. Colin, V. AvilaReese, and O. Valenzuela, Substructure and halo density profiles in a warm dark matter cosmology, Astrophys. J. 542 (2000) 622–630, [astroph/0004115].
 [76] A. D. Linde, Hybrid inflation, Phys. Rev. D49 (1994) 748–754, [astroph/9307002].
 [77] R. Fardon, A. E. Nelson, and N. Weiner, Supersymmetric theories of neutrino dark energy, JHEP 03 (2006) 042, [hepph/0507235].
 [78] S. Das, Sterile neutrino, hidden dark matter and their cosmological signatures, J. Phys. Conf. Ser. 405 (2012) 012011.
 [79] A. Sarkar, S. Das, and S. K. Sethi, How Late can the Dark Matter form in our universe?, JCAP 1503 (2015), no. 03 004, [arXiv:1410.7129].
 [80] D. J. E. Marsh and P. G. Ferreira, UltraLight Scalar Fields and the Growth of Structure in the Universe, Phys. Rev.D. D82 (2010) 103528, [arXiv:1009.3501].
 [81] W. Hu, R. Barkana, and A. Gruzinov, Cold and fuzzy dark matter, Phys. Rev. Lett. 85 (2000) 1158–1161, [astroph/0003365].
 [82] L. Amendola and R. Barbieri, Dark matter from an ultralight pseudoGoldsoneboson, Phys. Lett. B642 (2006) 192–196, [hepph/0509257].
 [83] R. Hlozek, D. Grin, D. J. E. Marsh, and P. G. Ferreira, A search for ultralight axions using precision cosmological data, Phys. Rev. D91 (2015), no. 10 103512, [arXiv:1410.2896].
 [84] M. Kaplinghat, Dark matter from early decays, Phys. Rev. D72 (2005) 063510, [astroph/0507300].
 [85] S. Profumo, K. Sigurdson, P. Ullio, and M. Kamionkowski, A Running spectral index in supersymmetric darkmatter models with quasistable charged particles, Phys. Rev. D71 (2005) 023518, [astroph/0410714].
 [86] K. Kohri and T. Takahashi, Cosmology with LongLived Charged Massive Particles, Phys. Lett. B682 (2010) 337–341, [arXiv:0909.4610].
 [87] A. Kamada, K. Kohri, T. Takahashi, and N. Yoshida, Effects of electrically charged dark matter on cosmic microwave background anisotropies, arXiv:1604.07926.
 [88] R. Foot and S. Vagnozzi, Dissipative hidden sector dark matter, Phys. Rev. D 91 (Jan., 2015) 023512, [arXiv:1409.7174].
 [89] B. T. Draine, Physics of the Interstellar and Intergalactic Medium. Princeton University Press, 2011.
 [90] G. Rybicki and A. Lightman, Radiative processes in astrophysics. WILEYVCH Verlag GmbH Co. KGaA, 1979.
 [91] A. P. Lightman, Double Compton emission in radiation dominated thermal plasmas, Astrophys. J. 244 (Mar., 1981) 392–405.
 [92] W. Hu and J. Silk, Thermalization and spectral distortions of the cosmic background radiation, Phys. Rev. D48 (1993) 485–502.
 [93] J. Bernstein and S. Dodelson, Aspects of the Zel’dovichSunyaev mechanism, Phys. Rev.D. 41 (Jan., 1990) 354–373.
 [94] Y. Zeldovich and R. Sunyaev, The Interaction of Matter and Radiation in a HotModel Universe, Astrophys. Space Sci. 4 (July, 1969) 301–316.
 [95] J. Chluba, Green’s function of the cosmological thermalization problem, Mon. Not. Roy. Astron. Soc. 434 (2013) 352, [arXiv:1304.6120].
 [96] W. Hu, D. Scott, and J. Silk, Reionization and cosmic microwave background distortions: A Complete treatment of second order Compton scattering, Phys. Rev. D49 (1994) 648–670, [astroph/9305038].
 [97] B. B. Nath and J. Silk, Heating of the intergalactic medium due to structure formation, Mon. Not. Roy. Astron. Soc. 327 (2001) L5, [astroph/0107394].
 [98] E. Komatsu and U. Seljak, The SunyaevZel’dovich angular power spectrum as a probe of cosmological parameters, Mon. Not. Roy. Astron. Soc. 336 (2002) 1256, [astroph/0205468].
 [99] N. Sehgal, P. Bode, S. Das, C. HernandezMonteagudo, K. Huffenberger, Y.T. Lin, J. P. Ostriker, and H. Trac, Simulations of the Microwave Sky, Astrophys. J. 709 (2010) 920–936, [arXiv:0908.0540].
 [100] J. Chluba, Distinguishing different scenarios of early energy release with spectral distortions of the cosmic microwave background, Mon. Not. Roy. Astron. Soc. 436 (2013) 2232–2243, [arXiv:1304.6121].
 [101] J. Chluba and D. Jeong, Teasing bits of information out of the CMB energy spectrum, Mon. Not. Roy. Astron. Soc. 438 (2014), no. 3 2065–2082, [arXiv:1306.5751].
 [102] N. Padmanabhan and D. P. Finkbeiner, Detecting dark matter annihilation with CMB polarization: Signatures and experimental prospects, Phys. Rev.D. 72 (July, 2005) 023508, [astroph/0503486].
 [103] J. Chluba, Could the cosmological recombination spectrum help us understand annihilating dark matter?, Mon. Not. Roy. Astron. Soc. 402 (Feb., 2010) 1195–1207, [arXiv:0910.3663].
 [104] E. Pajer and M. Zaldarriaga, A hydrodynamical approach to CMB distortion from primordial perturbations, JCAP 1302 (2013) 036, [arXiv:1206.4479].
 [105] P. J. E. Peebles, Recombination of the Primeval Plasma, Astrophys. J. 153 (July, 1968) 1.
 [106] J. Chluba and R. A. Sunyaev, Cosmological hydrogen recombination: influence of resonance and electron scattering, Astron. Astrophys. 503 (Aug., 2009) 345–355, [arXiv:0904.2220].
 [107] R. Sunyaev and J. Chluba, Signals From the Epoch of Cosmological Recombination, in Reviews in Modern Astronomy (S. Röser, ed.), vol. 21 of Reviews in Modern Astronomy, pp. 1–38, 2009.
 [108] J. Silk, Cosmic BlackBody Radiation and Galaxy Formation, Astrophys. J. 151 (Feb., 1968) 459.
 [109] P. J. E. Peebles and J. T. Yu, Primeval Adiabatic Perturbation in an Expanding Universe, Astron. Astrophys. 162 (Dec., 1970) 815.
 [110] N. Kaiser, Smallangle anisotropy of the microwave background radiation in the adiabatic theory, Mon. Not. Roy. Astron. Soc. 202 (Mar., 1983) 1169–1180.
 [111] S. Weinberg, Entropy Generation and the Survival of Protogalaxies in an Expanding Universe, Astron. Astrophys. 168 (Sept., 1971) 175.
 [112] W. Hu and J. Silk, Thermalization and spectral distortions of the cosmic background radiation, Phys. Rev.D. 48 (July, 1993) 485–502.
 [113] S. Dodelson, Modern cosmology. Academic Press, 2003.
 [114] C.P. Ma and E. Bertschinger, Cosmological perturbation theory in the synchronous and conformal Newtonian gauges, Astrophys. J. 455 (1995) 7–25, [astroph/9506072].
 [115] R. Khatri, R. A. Sunyaev, and J. Chluba, Mixing of blackbodies: entropy production and dissipation of sound waves in the early Universe, “ Astron. Astrophys.” 543 (July, 2012) A136, [arXiv:1205.2871].
 [116] R. Hložek, D. J. E. Marsh, D. Grin, R. Allison, J. Dunkley, and E. Calabrese, Future CMB tests of dark matter: ultralight axions and massive neutrinos, arXiv:1607.08208.
 [117] S. Bharadwaj and S. K. Sethi, Decaying Neutrinos and LargeScale Structure Formation, “Astrophys.J” 114 (Jan., 1998) 37–57, [astroph/9707143].
 [118] K. Ichiki, M. Oguri, and K. Takahashi, Constraints from the Wilkinson Microwave Anisotropy Probe on Decaying Cold Dark Matter, Physical Review Letters 93 (Aug., 2004) 071302, [astroph/0403164].