Issue 
A&A
Volume 540, April 2012



Article Number  A124  
Number of page(s)  10  
Section  Cosmology (including clusters of galaxies)  
DOI  https://doi.org/10.1051/00046361/201118194  
Published online  11 April 2012 
Does BoseEinstein condensation of CMB photons cancel μ distortions created by dissipation of sound waves in the early Universe?
^{1} Max Planck Institut für Astrophysik, KarlSchwarzschildStr. 1, 85741 Garching, Germany
email: khatri@mpagarching.mpg.de
^{2} Space Research Institute, Russian Academy of Sciences, Profsoyuznaya 84/32, 117997 Moscow, Russia
^{3} Institute for Advanced Study, Einstein Drive, Princeton, New Jersey 08540, USA
^{4} Canadian Institute for Theoretical Astrophysics, 60 St George Street, Toronto, ON M5S 3H8, Canada
Received: 3 October 2011
Accepted: 12 February 2012
The difference in the adiabatic indices of photons and nonrelativistic baryonic matter in the early Universe causes the electron temperature to be slightly lower than the radiation temperature. Comptonization of photons with colder electrons results in the transfer of energy from photons to electrons and ions, owing to the recoil effect (spontaneous and induced). Thermalization of photons with a colder plasma results in the accumulation of photons in the RayleighJeans tail, aided by stimulated recoil, while the higher frequency spectrum tries to approach Planck spectrum at the electron temperature T_{γ}^{final} = T_{e} < T_{γ}^{initial}; i.e., BoseEinstein condensation of photons occurs. We find new solutions of the Kompaneets equation describing this effect. No actual condensate is, in reality, possible since the process is very slow and photons drifting to low frequencies are efficiently absorbed by bremsstrahlung and double Compton processes. The spectral distortions created by BoseEinstein condensation of photons are within an order of magnitude (for the present range of allowed cosmological parameters), with exactly the same spectrum but opposite in sign, of those created by diffusion damping of the acoustic waves on small scales corresponding to comoving wavenumbers 45 < k < 10^{4} Mpc^{1}. The initial perturbations on these scales are completely unobservable today due to their being erased completely by Silk damping. There is partial cancellation of these two distortions, leading to suppression of μ distortions expected in the standard model of cosmology. The net distortion depends on the scalar power index n_{s} and its running dn_{s}/dlnk, and may vanish for special values of parameters, for example, for a running spectrum with, n_{s} = 1,dn_{s}/dlnk = − 0.038. We arrive at an intriguing conclusion: even a null result, nondetection of μtype distortion at a sensitivity of 10^{9}, gives a quantitative measure of the primordial smallscale power spectrum.
Key words: radiation mechanisms: thermal / scattering / cosmic background radiation / cosmological parameters / cosmology: theory / early Universe
© ESO, 2012
1. Introduction
In the early Universe we have two strongly interacting fluids: cosmic microwave background (CMB) radiation and plasma (=ions + electrons). Electrons and CMB are strongly coupled due to the Compton process with the rate of energy transfer from CMB to electrons given by R_{compton} ~ 1.4 × 10^{20}(1 + z)^{4} s^{1}. This is much faster than the expansion rate , H(z) ~ 2.1 × 10^{20}(1 + z)^{2} s^{1}^{1}, and keeps T_{γ} − T_{e} ≪ T_{γ} at redshifts z ≳ 200, where T_{γ} is the effective temperature of radiation (to be defined below) and T_{e} is the electron temperature. Ions and electrons exchange energy via Coulomb collisions at a rate R_{coulomb} ~ 3 × 10^{9}(1 + z)^{3/2} s^{1} (Lifshitz & Pitaevskii 1981) keeping T_{e} − T_{i} ≪ T_{e}, where T_{i} is the temperature of ions. Thus we have T_{γ} ≈ T_{e} ≈ T_{i} to high precision at z > 200. Nonrelativistic matter has an adiabatic index of 5/3, while the adiabatic index of radiation is 4/3. Adiabatic cooling due to the expansion of the Universe thus makes the matter cool faster than the radiation, while Comptonization tries to maintain matter and radiation at the same temperature by transferring energy from radiation to matter. As a result, the energy of radiation decreases, while the number of photons is conserved (neglecting absorption of photons by bremsstrahlung and double Compton). This effect was discussed by Zeldovich et al. (1968) and Peebles (1968) and leads to a decrease in the temperature of plasma compared to the radiation temperature at z ≲ 200, when Comptonization becomes inefficient because of depletion of electrons from recombination. It also drives the 21cm line spin temperature to below the radiation temperature, raising the possibility of 21cm absorption of CMB photons by neutral hydrogen. This transfer of energy from CMB to matter happens at all redshifts, as long as there is significant free electron number density, and results in deviations in CMB spectrum from the blackbody. Recently, Chluba & Sunyaev (2012) have numerically studied these deviations from the blackbody in the CMB spectrum. Below we demonstrate that the problem of spectral deviations of CMB due to loss of energy to plasma, in the case of small spectral deviations and small Compton parameter y, has a simple analytic solution, when photon production and destruction is neglected. The resulting flow of photons towards low frequencies as the spectrum tries to approach the Planck spectrum is in fact the BoseEinstein condensation of photons. The BoseEinstein condensation of photons unfortunately does not progress very far in reality, because lowfrequency photons are efficiently absorbed by bremsstrahlung and double Compton processes, and the Compton process freezes out as the electron density falls owing to the expansion of the Universe and recombination.
Early papers about energy release in the early Universe were concerned with exotic sources such as annihilation of matter and antimatter, primordial turbulence, decay of new unstable particles, unwinding of topological defects like domain walls, and cosmic strings. Experiments beginning with COBE FIRAS (Fixsen et al. 1996) have been unable to find any significant distortions in the CMB from blackbody. New proposals like Pixie (Kogut et al. 2011) are demonstrating that a tremendous increase in the sensitivity is possible in the future experiments. Pixie is proposed to be 1000 times more sensitive than COBE FIRAS. There is hope that this is still not the last word, and even higher sensitivity might become possible in the future. Under these circumstances we decided to check what minimum levels of deviations from blackbody spectrum are expected in the standard cosmological model (Komatsu et al. 2011). In this paper we only consider the spectral distortions arising before the end of the dark ages and beginning of reionization. We demonstrate that, in the absence of decay and annihilation of new unknown particles or any other new physics beyond the standard model, there are only three key reasons for significant global spectral distortions in the CMB.

1.
BoseEinstein condensation of CMB photons due to the difference in the adiabatic indices of nonrelativistic plasma and photons.

2.
The energy release due to dissipation of sound waves and initial perturbations due to Silk damping (Silk 1968) and the secondorder Doppler effect due to nonzero peculiar velocity of electrons and baryons in the CMB rest frame, both leading to the superposition of blackbodies (Zeldovich et al. 1972) in the electron rest frame.

3.
The cosmological recombination radiation from hydrogen and helium (Dubrovich 1975; Peebles 1968; Zeldovich et al. 1968).
We omit the distortions caused by the cosmological recombination process (e.g., see Chluba & Sunyaev 2006; Sunyaev & Chluba 2009, for more details and references therein) in the discussion below. They have narrow features and different continuum spectrum and can be distinguished from y and μ distortions from the first two mechanisms. Cosmological recombination radiation can, however, become comparable to or larger than those discussed here (Chluba & Sunyaev 2012), if the distortions from the first and second processes partially cancel. We also note that for any experiment with finite beam size there will be a mixing of blackbodies in the beam owing to the temperature fluctuations on the last scattering surface on scales smaller than the beam size, leading to inevitable y distortions of magnitude ~ (ΔT/T)^{2} ~ 10^{9} − 10^{10} (Chluba & Sunyaev 2004).
In this paper we find new solutions for the Kompaneets equation describing the CMB spectral distortions arising from process 1 described above. A second mechanism results in a ydistortion and heating of the electrons at low redshifts. At high redshifts Comptonization of CMB on the hotter electrons converts the ytype distortion to a μtype distortion (Sunyaev & Zeldovich 1970b; Illarionov & Sunyaev 1975b). We demonstrate below that thesetwo sources of distortions work against each other, individually resulting in spectral distortions with opposite signs but the same spectral shape in any part of the CMB spectrum.
The cooling of electrons and the corresponding spectral distortions are easy to calculate and only depend on the standard cosmological parameters, such as baryon to photon number density n_{B}/n_{γ} and helium fraction, which decide the amount of energy losses by the CMB, the Hubble constant (H_{0}), and densities of constituents of the Universe, which in turn determine the expansion rate and thus the efficiency of Comptonization. On the other hand, the energy released by dissipation of sound waves crucially depends on the power in smallscale fluctuations, and thus the spectral index (n_{s}) in the standard cosmological model, in addition to the other wellmeasured parameters of the standard cosmological model (Komatsu et al. 2011). It is interesting to note that for a spectrum with constant n_{s} the energy release due to dissipation exceeds the energy losses to adiabatic cooling of baryons, and there is net heating of electrons. For a primordial spectrum with a running spectral index, the role of Comptonization with colder electrons and BoseEinstein condensation can become dominant, with the spectral distortions changing sign, and there can be net cooling of electrons and a corresponding decrease of CMB entropy per baryon (specific entropy). Additional increase in entropy during the recombination epoch due to superposition of blackbodies, free streaming as well as Silk damping and the secondorder Doppler effect only produce ytype distortions and can be distinguished from the μ distortions created in the earlier epoch.
2. Thermodynamic equilibrium in the early Universe
We can obtain the energy losses due to adiabatic cooling of baryons under the assumption that Compton scattering, bremsstrahlung, and double Compton scattering can maintain full thermodynamic equilibrium between the electrons and photons. This is true to high accuracy at z > 10^{6} while at lower redshifts, although bremsstrahlung and double Compton scattering cannot help in restoring a Planck spectrum at all frequencies, Comptonization keeps the electron temperature equal to the effective radiation temperature (cf. Eq. (16)) to high precision until z ~ 200. In the thermodynamic equilibrium, entropy is conserved, and we can use this to calculate the energy losses due to the expansion of the Universe without referring to Comptonization, since it is not really important which physical process is responsible for the equilibrium.
We start with the thermodynamic relation, (1)where S is the entropy, U total thermal energy, P pressure, V volume, M chemical potential, and N the number of particles. T_{γ} is the common temperature of photons, ions, and electrons. We can ignore the last term for ions and electrons since their number is fixed after big bang nucleosynthesis ends^{2}. For photons as well we can ignore the last term if we assume that their chemical potential is 0, which is true in full thermodynamic equilibrium. With this assumption we can write the equation for total entropy per baryon, σ = s/n_{B}, where s is the entropy density and n_{B} is the baryon number density: Here a_{R} is the radiation constant and N the number of nonrelativistic particles per baryon, so that N n_{B} = n_{eb}. We now integrate Eq. (2)to obtain (5)where C is an arbitrary constant of integration, which is not important for our calculation. The first term above is the contribution from photons, and it is clear that in the absence of second term, T_{γ} ∝ (1 + z). The second term is the contribution of nonrelativistic particles and causes the temperature to drop slightly faster in order to conserve the entropy per baryon. We can write T_{γ} = T(1 + t), where T = T_{i}(1 + z)/(1 + z_{i}) is the background temperature proportional to 1 + z, and t ≪ 1 is the fractional deviation from this law. We can take the initial deviation t(z_{i}) = 0 at initial redshift z_{i} and then calculate the subsequent energy losses to adiabatic cooling of baryons at later redshifts. Thus we can write σ, correct to leading order in t, as (6)Equating the above to the initial entropy σ_{i} = σ(z_{i}) with t ≡ 0,T ≡ T_{i}, we get upon solving for t (and ignoring the last term in Eq. (6)compared to the second term) (7)with ΔE/E = 4t, since E_{γ} ∝ T^{4}, and . We note that t is negative since z_{i} > z and there is net energy loss. ln(z_{max}/z_{min}) ~ 10 for the redshift range of interest, 200 < z < 2 × 10^{6}, where a mixture of μ and ytype distortions are produced. This result for the total energy extraction is consistent with the estimate obtained by Chluba (2005) and Chluba & Sunyaev (2012).
To estimate the magnitude of total distortions, μ and ytype, we take z_{max} = z_{i} = 2 × 10^{6}, since at greater redshifts spectral distortions are rapidly destroyed by bremsstrahlung and double Compton processes and Comptonization. The latter redistributes photons over the whole spectrum. For z_{min} = 200 we get ΔE/E = 5.4 × 10^{9}. This transfer of energy from radiation to baryons due to a difference in their adiabatic indices results in an inevitable distortion and BoseEinstein condensation of photons. We show below that the distortion has a magnitude of Y_{BEC} ~ (1/4) ΔE/E ~ 10^{9}. We define parameter Y_{BEC} and prove this result below. The corresponding deviation of the electron temperature from the equilibrium temperature of electrons in the radiation field is described by Eq. (16)and is (T_{γ} − T_{e})/T_{γ} ≈ 10^{12} at z = 10^{6} growing to ~0.01 at z = 500.
In the Universe there are also inevitable processes like the dissipation of acoustic waves and SZ effect from reionization, which would lead to a normal SZ effect of similar or greater magnitude. Due to the additivity of small spectral distortions, these effects will cancel the Y_{BEC} distortions, and this diminishes the hope that this BoseEinstein condensation will ever be observed. Nevertheless it is important to stress that in the standard cosmology there was an epoch 10^{6} > z > 1000 when BoseEinstein condensation of CMB photons due to the difference in the adiabatic indices of matter and radiation was able to create a peculiar deviation of CMB spectrum from the blackbody.
3. Alternative direct calculation of energy losses in CMB
We would like to remind the reader that adiabatic cooling of baryons leads to the cooling of radiation, but Comptonization itself conserves the number of photons. We can write the equation for the evolution of average thermal energy density E_{B} in baryons before recombination as (8)The derivative in the first term on the righthand side, dn_{eb}/dz = 3n_{eb}/(1 + z), is just the decrease in number density of particles (ions and electrons) with redshift z due to the Hubble expansion, with n_{eb} = ρ_{b}/μ_{mol}, where ρ_{b} is the baryon mass density and μ_{mol} the mean molecular weight. The change in temperature has contributions from adiabatic cooling and also from energy gained from CMB by Comptonization: (9)where S_{Compton} = 4k_{B}E_{γ}n_{e}σ_{T}(T_{γ} − T_{e})/m_{e}cH(z)(1 + z) is the energy transfer rate per unit volume from radiation to baryons by Compton scattering, n_{e} is the number density of free electrons, σ_{T} the Thomson scattering cross section, k_{B} the Boltzmann constant, E_{γ} the energy density of radiation, c the speed of light, m_{e} the mass of electron, and H(z) the Hubble parameter. That redshift decreases with increasing time, so the terms with “ + ” sign are cooling terms and terms with “ − ” sign are heating terms. The change in photon energy density can therefore be written as (10)where the first term is just the adiabatic cooling due to the expansion of the Universe. We can estimate the energy transfer, S_{Compton}, by noting that the baryon temperature T_{e} = T_{γ} ∝ 1 + z to a high accuracy until z ~ 200, where T_{γ} = 2.725(1 + z) K is the CMB temperature. Using this to evaluate the total derivative of electron temperature dT_{e}/dz on the lefthand side of Eq. (9)and the first term on the righthand side, we get (11)Thus we have fractional rate at which energy is lost by radiation to baryons (12)We can now calculate the total fractional energy losses of radiation that contribute to the spectral distortions between redshifts z_{min} and z_{max}, (13)We have integrated dE/E, since immediate distortions are proportional to ΔE/E, and the distortions can be added linearly if they are small. This is the same result as Eq. (7).
4. Kompaneets equation
The interaction of radiation with electrons through Compton scattering or Comptonization is described by the Kompaneets equation (Kompaneets 1956) in the FokkerPlanck approximation, when the energy transfer in each scattering is small compared to temperature, and the incoming photon distribution is wide compared to the width of the scattering kernel: (14)where we have defined the Compton parameter (15)which is convenient to use instead of time or redshift. We start our calculation at the reference redshift z_{max}. One can furthermore introduce the Compton equilibrium electron temperature in a radiation field (Levich & Sunyaev 1971; Zel’Dovich & Levich 1970) (16)Here T_{BB} is a reference temperature which is equal to the radiation blackbody temperature if the initial radiation field is a Planck spectrum, and x = hν/k_{B}T_{BB} is the dimensionless frequency with ν the frequency of photons. The equations written in this way factor out the expansion of the Universe and are applicable to cosmology, as well as nonexpanding astrophysical systems. This equation and the Kompaneets equation form a coupled system to be solved simultaneously.
The three terms in the inner brackets of Eq. (14)describe recoil (n), induced recoil (n^{2}), and the Doppler effect of the thermal motion of electrons (T_{e}/T_{BB}∂n/∂x). For y ≪ 1 and an initial blackbody spectrum with temperature T_{BB}, n_{BB} = 1/(e^{x} − 1), a simple firstorder correction to the blackbody spectrum can be found (Zeldovich & Sunyaev 1969). By approximating^{3}n + n^{2} ≈ − ∂n/∂x in Eq. (14)we get (17)where Δ ≡ T_{e}/T_{BB} − 1. This particular form of equation emphasizes that the Doppler term dominates over recoil and gives rise to the spectral distortion by boosting lowenergy photons to higher frequency. The above equation can be reduced to the diffusion or heat equation (Zeldovich & Sunyaev 1969) by changing variables, for which the Green’s function is well known.
For small distortions, the solution takes a particularly simple form and can be found by substituting n_{BB} on the righthand side of Eqs. (14)or (17): (18)where we have defined . The SZ parameter is related to fractional energy release by the formula Y_{SZ} = (1/4)ΔE_{γ}/E_{γ}. This can be obtained by integrating the above spectrum over all frequencies.
5. New solution of Kompaneets equation for Comptonization of photons with colder plasma
The Kompaneets equation describes, to lowest order, the competition between the Doppler boosting of lowenergy photons to high energies and the down scatter of highenergy photons to low energies by recoil and stimulated recoil. The initial radiation spectrum is blackbody with temperature T_{BB}. If the electrons are also at the same temperature, T_{e} = T_{BB}, then these two effects cancel each other out exactly. If T_{e} > T_{BB}, the Doppler boosting becomes stronger than recoil and we have the normal SZ effect in the limit of small y parameter. If T_{e} < T_{BB} then Doppler boosting is weaker and we have a net movement of photons from high to low energies, which is the complete opposite of the SZ effect. But since at T_{e} = T_{BB} the two effects exactly balance each other out, at linear order the spectral distortions for T_{e} < T_{BB} would be exactly the same as for T_{e} > T_{BB}, just with the opposite sign.
For T_{e} < T_{BB}, we can write the approximate Kompaneets equation by approximating ∂n/∂x ≈ − (n + n^{2}) in Eq. (14), (19)Analytic solutions can be obtained by the standard method of characteristics for the recoil (n) and induced recoil (n^{2}) terms individually, and they consist of photons moving down the frequency axes at a speed proportional to x^{2} for the recoil term (Arons 1971; Illarionov & Syunyaev 1972) and at a speed proportional to the x^{2}n for the induced recoil term (Zel’Dovich & Levich 1969; Syunyaev 1971). To get the linear solution, we can just substitute n(x) = n_{BB}(x) on the righthand side of Eq. (19)and define the parameter for the amplitude of distortion in the case T_{e} < T_{BB} as , the result is Eq. (18)with Y_{SZ} = − Y_{BEC}. We should emphasize that the linear solution, Eq. (18), can be arrived at in a straightforward way by just substituting n(x) = n_{BB}(x) on the righthand side of the full Kompaneets equation, Eq. (14), without any restriction on T_{e}. However the dominant physical effects when T_{e} > T_{BB} (Doppler) and when T_{e} < T_{BB} (recoil and induced recoil) are completely different and that in equilibrium they balance each other out exactly gives us the same mathematical form of the linear solution in both cases.
Fig. 1 Difference in intensity from the blackbody radiation for normal SZ effect with Y_{SZ} = 10^{9} and negative SZ effect with Y_{BEC} = 10^{9}. 
Fig. 2 Difference in intensity from the blackbody radiation for μ distortion defined by the equation with n(ν) = 1/(e^{hν/kBTμ + μ} − 1). T_{μ} is the temperature to which the spectrum relaxes for an initial energy addition/loss of ΔE/E with ΔT_{μ}/T = 0.64ΔE/E = ± 2.5 × 10^{9}. ΔE/E = ± 4 × 10^{9} is the energy addition/loss that gives rise to SZ distortion of 10^{9}. The dimensionless chemical potential μ is given by μ = 2.2ΔT/T = ± 5.6 × 10^{9}. 
In the RayleighJeans part of the spectrum we now have an increase in the brightness temperature of ΔT/T = +2Y_{BEC}, which is independent of frequency and thus maintains the RayleighJeans shape of the spectrum. In the Wien part we have a decrease in intensity ΔI/I = Δn/n ≈ − x^{2}Y_{BEC}. Figure 1 shows the difference in intensity from the blackbody radiation for normal SZ effect with Y_{SZ} = 10^{9} and the negative SZ effect with Y_{BEC} = 10^{9}. Figure 2 shows the spectrum that would be achieved at a high value of the Compton parameter y. It is a BoseEinstein spectrum with the dimensionless chemical potential^{4}μ defined by the equation with n(ν) = 1/(e^{hν/kBTμ + μ} − 1), where T_{μ} is the temperature to which the spectrum relaxes for an initial energy addition/loss of ΔE/E with ΔT_{μ}/T = 0.64ΔE/E = ±2.5 × 10^{9} (Illarionov & Sunyaev 1975a), and ΔE/E = ±4 × 10^{9} is the energy addition/loss that gives rise to SZ distortion of 10^{9}. The chemical potential μ is given by^{5}μ = 2.2ΔT/T = ± 5.6 × 10^{9}. This is the spectrum that an initial spectrum with Y_{BEC},Y_{SZ} = 10^{9} will approach at high y at x ≫ μ ~ 10^{9}. The frequency at which the distortion crosses zero is at x = 2.19 compared to x = 3.83 for the SZ distortion in Fig. 1.
The ytype and μtype distortions expected in the early Universe (calculated in the later sections) are compared with the cosmological recombination spectrum (RubiñoMartín et al. 2006) in Fig. 3. Clearly μ type distortions have a different spectral shape than the recombination radiation (both from hydrogen (Chluba & Sunyaev 2006) and helium (RubiñoMartín et al. 2008)) and ytype distortions and can be distinguished from the last two. This is very important because the information in the μ type distortions about the early Universe physics can be extracted. On the other hand, the ytype distortions from the early Universe get swamped by the much larger ytype distortions from the low redshifts and the two contributions are difficult to separate. The μ distortions expected from the early Universe also have higher magnitude than the recombination spectrum in the RayleighJeans part of the spectrum but has no quasiperiodic structure like the cosmological recombination radiation.
Fig. 3 Comparison of positive and negative μtype distortions expected in the early Universe (to be calculated in later sections) with the ytype distortions before reionization and the cosmological recombination spectrum (taken from Chluba & Sunyaev 2006; RubiñoMartín et al. 2008) for illustration. 
Fig. 4 Evolution of initial spectrum at y = 0 given by Eq. (18)with − Y_{SZ} = Y_{BEC} = 10^{5}. ΔT/T ≡ T_{eff}/T − 1 is plotted, where T_{eff} is the temperature of a blackbody spectrum with the corresponding intensity at a particular frequency. BoseEinstein spectrum defined by n(ν) = 1/(e^{hν/kBTe + μ} − 1) with negative μ is labeled “ve BE”, red dashed, doubledotted line in the figure. The distortions are positive in the RayleighJeans region (labeled “+ve”) on the left side and negative in the Wien region (labeled “ve”) on the right side of the figure. 
This flow of photons toward lower frequencies as the spectrum tries to approach the Planck spectrum due to recoil, and induced recoil is BoseEinstein condensation of photons (Illarionov & Sunyaev 1975a; Landau & Lifshitz 1980). We show the evolution of the spectrum (solution of the coupled system of Eqs. (14)and (16)) starting with the initial distortion given by Eq. (18)with − Y_{SZ} = Y_{BEC} = 10^{5} in Fig. 4. In the RayleighJeans region intensity I_{ν} ∝ T_{eff} and therefore ΔT/T = ΔI/I. The initial evolution is similar to that of a spectrum with positive Y_{SZ} with the photon distribution approaching a BoseEinstein spectrum defined by n(ν) = 1/(e^{hν/kBTe + μ} − 1) with negative μ (marked “ve BE” in the figure). We should emphasize that there is no singularity in the actual solutions of the Kompaneets equation plotted above. The singularity is just in the mathematical formula, which correctly describes the spectrum at high frequencies, x ≫ μ. The actual spectrum deviates from the BoseEinstein spectrum near the singularity (positive everywhere in the RayleighJeans region) and can be described by a chemical potential decreasing in magnitude with decreasing frequency. The evolution at y > 1 is therefore very different from the positive Y_{SZ} case. For x ≫ μ the spectrum is the BoseEinstein spectrum at the electron temperature. For x ≪ μ there is an excess of photons compared to the BoseEinstein spectrum which grows with time, a feature of BoseEinstein condensation. The fractional change in temperature ΔT/T takes a very simple form in the RayleighJeans region (x ≪ 1) for a BoseEinstein spectrum and can be understood as follows. It is given, in terms of frequency referred to the electron temperature, x_{e} ≡ hν/kT_{e}, by (20)For μ ≪ x_{e} ≪ 1 we have ΔT/T = − μ/x_{e}. Thus the fractional temperature deviation is positive for negative μ and negative for the positive μ. For x_{e} greater than and close to μ,μ < 0 the spectral distortions can become very large and exceed unity. Since the highfrequency spectrum is in equilibrium at the electron temperature by y = 10, the subsequent evolution is very slow, and μ decreases while the electron temperature increases slowly as the photon distribution approaches a Planck spectrum with the extra photons accumulating at low frequencies; i.e., BoseEinstein condensation happens. The lowfrequency spectrum approaches the stationary solution of Kompaneets equation with only the induced recoil term, n(x) ∝ 1/x^{2}, with a continuous flow of photons towards x = 0 (Syunyaev 1971). Then, T_{e}/T − 1 just increases from −2.558 × 10^{5} to −2.555 × 10^{5} in going from y = 10 to y = 100 for the chosen energy losses of ΔE/E = 4 × 10^{5}. The effect discussed below in the real Universe is of magnitude ~10^{9}.
Fig. 5 Regions of the blackbody photosphere probed by different experiments, COBE FIRAS with sensitivity μ = 10^{5}, Pixie with sensitivity μ = 10^{8} and a hypothetical experiment with sensitivity of μ = 10^{10}. Shaded regions below the curves are allowed and inaccessible to these experiments. We also plot the ratio of final μ distortion today to the initial μ distortion as a function of energy injection redshift. We define the blackbody surface as the energy injection redshift such that μ^{final}/μ^{initial} = 1/e, z_{bb} = 1.98 × 10^{6}. 
6. BoseEinstein condensation of CMB photons in the early Universe
In the early Universe, before z ~ 2 × 10^{6} double Compton and Comptonization destroy any spectral distortions and maintain the Planck spectrum of CMB. For small distortions due to single, quasiinstantaneous episode of energy release, we can write the ratio of finaltoinitial μ as an exponential function of redshift defined by the square root of the product of Comptonization and absorption rates (Sunyaev & Zeldovich 1970b). For a double Compton process, this formula gives (Danese & de Zotti 1982) (21)where z_{dc} ≈ 1.98 × 10^{6} defines the “surface of the blackbody photosphere”. We call this function G(z) the blackbody visibility function for spectral distortions. We call the region z > z_{dc} where the initial μ distortion can be reduced by a factor of more than e the “blackbody photosphere”. Thus inside the blackbody photosphere, a Planckian spectrum can be established efficiently. Figure 5 shows how the μ distortion at high redshifts decreases due to the double Compton and Compton scattering. The regions of the blackbody photosphere allowed by different experiments are shaded.
At redshifts 10^{5} < z < 2 × 10^{6}, we have the Compton parameter y > 1, and the BoseEinstein spectrum with a negative chemical potential is established at x ≳ 0.01^{6}. At lower frequencies, bremsstrahlung and double Compton create a Planck spectrum corresponding to the electron temperature. At z < 10^{4} the Compton parameter y < 0.01 and the distortions created can be described by the linear solution of Eq. (18). Thus the net distortions created in the early Universe owing to adiabatic cooling of baryons are a linear superposition of distortions corresponding to different values of y, and a few of them are shown in Fig. 6. We show the evolution of actual spectrum (Eqs. (14) and (16)) including the effects of bremsstrahlung and double Compton processes (Burigana et al. 1991; Chluba & Sunyaev 2012; Hu & Silk 1993; Illarionov & Sunyaev 1975a,b) in the ΛCDM cosmology with WMAP parameters (Komatsu et al. 2011). The reference spectrum is the blackbody CMB spectrum at z = 10^{7}. At x < 0.01 the spectrum is dominated by bremsstrahlung and double Compton processes that destroy the condensing lowfrequency photons. The effect of bremsstrahlung compared to Comptonization becomes stronger with decreasing redshift and can be felt at successively higher frequencies.
At z ≳ 5 × 10^{4} the spectral distortions in the x ≳ 0.01 region can relax to a BoseEinstein spectrum due to high value of parameter y. Contributions from lower redshifts are described by the analytic solution of Eq. (18). The distortions created at z < 10^{4} are just the negative of the usual SZ effect and would be completely overwhelmed by the SZ effect from reionization at z ~ 10 which is expected to be Y_{SZ} = (k_{B}T_{e}/m_{e}c^{2})τ_{ri} ~ 10^{6} × 0.1 = 10^{7}, where τ_{ri} ~ 0.1 is the optical depth to the last scattering surface due to reionization and T_{e} ~ 10^{4} K is the average electron temperature during reionization. The distortions created by energy loss at z > 10^{5} can, however, be destroyed only by energy injection at z > 10^{5} since at lower redshift Comptonization cannot create μtype distortions.
Fig. 6 Evolution of spectral distortions in CMB in ΛCDM cosmology. The values of total Compton parameter y = y(100,z) for different curves are (z = 1.5 × 10^{6}, y = 109), (z = 4.6 × 10^{5}, y = 10), (z = 1.5 × 10^{5}, y = 1), (z = 5.9 × 10^{4}, y = 0.16), (z = 1.2 × 10^{4}, y = 0.005), (z = 1760,y = 3 × 10^{5}), (z = 840, y = 1.8 × 10^{8}) and (z = 220, y = 1.2 × 10^{10}). Reference spectrum is the blackbody CMB spectrum at z = 10^{7}. At low frequencies the spectrum is dominated by double Compton and bremsstrahlung which maintain Planck spectrum at the electron temperature. The zero point at low frequencies moves to the right at low redshifts because bremsstrahlung dominates over the Compton scattering at higher and higher frequencies and absorbs photons, although both are slower than the expansion rate. 
7. Energy release from dissipation of acoustic waves during radiation domination
In the standard cosmology the main source of energy injection at high redshifts is the dissipation of acoustic waves in the baryonphoton plasma (Sunyaev & Zeldovich 1970b) due to Silk damping (Silk 1968). Sunyaev & Zeldovich (1970a) first proposed using the resulting spectral distortions to measure the spectral index and power in the smallscale fluctuations, which themselves do not survive. Later during recombination the secondorder Doppler effect due to nonzero electron/baryon velocity in the CMB rest frame also becomes important, but at that time only ytype distortions similar to the SZ effect arise, and μ distortions cannot be created. We can estimate the amount of μ distortions created by dissipation of acoustic waves following Hu et al. (1994) by calculating how the total power in the density fluctuations changes with time because of photon diffusion. The comoving energy density in acoustic waves in the photonbaryon plasma (neglecting the baryon energy density) is given by , where ρ_{γ} is the comoving photon energy density, is sound speed squared, δ_{γ} is the photon density perturbation at position x, and angular brackets denote the ensemble average: (22)where , and is the initial power spectrum. It was shown by Mukhanov & Chibisov (1981) that primordial fluctuations from inflation can have a spectrum deviating from the scaleinvariant HarrisonZeldovich spectrum, with spectral index n_{s} < 1. In the cyclic ekpyrotic models n_{s} > 1 is also possible (Lehners et al. 2007). WMAP gives the constraint on curvature perturbation in comoving gauge, ζ. This is related to the gravitational perturbation ψ in the radiation era (assuming neutrinos are free streaming) by the relation ψ = ζ/(2/5R_{ν} + 1.5), where R_{ν} = ρ_{ν}/(ρ_{γ} + ρ_{ν}) ≈ 0.4, ρ_{ν} is the neutrino energy density and (Ma & Bertschinger 1995). Thus and , k_{0} = 0.002 Mpc^{1}, A_{ζ} = 2.4 × 10^{9} (Keisler et al. 2011; Komatsu et al. 2011), r = dn_{s}/dlnk is the running of the index. The transfer function for modes well inside the horizon before recombination is given by (23)and the diffusion scale by (Kaiser 1983; Weinberg 2008) where R ≡ 3ρ_{b}/4ρ_{γ}, ρ_{b} is the baryon energy density.
Replacing cos^{2}(kr_{s}) with its average value over an oscillation of 1/2, we get the energy release per unit redshift (26)For a running index, the above integration must be done numerically. We are interested in radiationdominated epoch where μ distortions are generated, and we have and n_{e}(z) = (n_{H0} + 2n_{He0})(1 + z)^{3} ≡ n_{e0}(1 + z)^{3}. In addition we have and d(1/k_{D}^{2})/dz = − 3A_{D}(1 + z)^{4}, and the above equation simplifies to where (27)The redshift and n_{s} dependence obtained above matches that of Hu et al. (1994). The diffusion scale k_{D} = 46 Mpc^{1} at z = 5 × 10^{4} and k_{D} = 10^{4} Mpc^{1} at z = 2 × 10^{6}. Thus we are probing the primordial fluctuations on very small scales that are not accessible in any other way.
Thus for a HarrisonZeldovich spectrum with n_{s} = 1 (28)We note the surprising fact that this redshift dependence is exactly the same as for the energy losses due to adiabatic cooling Eq. (13), as already pointed out by Chluba & Sunyaev (2012). For the currently preferred value of n_{s} = 0.96 we have (29)Figure 7 compares the rate of energy release multiplied by the blackbody visibility function, G(z)(1 + z)dQ/dz/E_{γ} for different initial power spectra with the energy losses due to the adiabatic cooling of baryons, (1 + z)S_{Compton}/E_{γ} obtained by solving Eq. (9)for the standard recombination history (Seager et al. 2000) calculated using the effective multilevel approach and taking recent corrections into account (AliHaïmoud & Hirata 2011; Chluba & Thomas 2011). The spike due to increase in Silk damping during recombination is clearly noticeable^{7}. The smaller spikes due to HeII and HeI are also noticeable. Recombination of each species leads to a decrease in the number of particles, hence to adiabatic losses. But the photon diffusion length, and the associated energy release, increases owing to a decrease in the electron fraction. The transition from radiation to matter domination results in unnoticeable change in the slope of the curves. Even after recombination photons continue to mix on horizon scales because of free streaming, while the energy losses from Comptonization drop sharply as a result of depletion of electrons. We switch to the freestreaming solution described in Appendix A when diffusion length 2π/k_{D} becomes equal to the comoving particle horizon. In the blackbody photosphere region, the blackbody visibility function G(z) makes the curves drop sharply. Global energy release and ytype distortions during the free streaming epoch after recombination are due to the superposition of blackbody spectra on horizon scales. The energy injected into ytype distortion during recombination in the peak (800 < z < 1500) is ΔE/E ~ 10^{8} for n_{s} = 0.96 (see also footnote 7). The y distortion due to mixing of blackbodies during free streaming after recombination and up to reionization, 20 < z < 800, calculated using the result in the Appendix A, is ΔE/E ~ 10^{10}. These are, however, difficult to separate from much the larger thermal and nonlinear Doppler y distortion from reionization of Y_{SZ} ~ 10^{7}.
Fig. 7 Sketch of fractional rate of energy release due to Silk damping and free streaming for different initial power spectra. Also shown for comparison is the rate of energy loss due to adiabatic cooling of baryonic matter. 
Energy injection in μ distortions during 5 × 10^{4} < z < 2 × 10^{6} for different initial power spectra without running compared with energy losses due to BoseEinstein condensation.
Energy injection in μ distortions during 5 × 10^{4} < z < 2 × 10^{6} for different initial power spectra with running compared with energy losses due to BoseEinstein condensation.
Energy releases for different power spectra, including the ones with running spectral index, are summarized in Tables 1 and 2. The energy losses from adiabatic cooling of baryons during the same time period is ~2.2 × 10^{9}. This equals the energy release, for example, for a spectral index of n_{s} = 1.0, dn_{s}/dlnk = − 0.038 for a running spectrum. The net μ distortion as a function of spectral index n_{s} without running is shown in Fig. 8 and for n_{s} = 1.0,1.02,1.05 as a function of the running of the index in Fig. 9. For low values of n_{s} and high negative values of running, BoseEinstein condensation dominates, and net distortion approaches a constant value of −3 × 10^{9}. A detection or nondetection at a sensitivity of 10^{9} of μtype distortion is thus a very sensitive probe of primordial power spectrum on small scales. Given the importance of the effective heating rate for the amplitude of the net distortion, it is important to carry out a more careful and refined calculation of the problem, which we present elsewhere (Chluba et al. 2012).
Fig. 8 μ = 1.4ΔE/E as a function of spectral index n_{s} without running. 
Fig. 9 μ = 1.4ΔE/E as a function of running of the spectral index r for three values of spectral index n_{s}. 
Over the years a number of ground and balloonbased experiments that aimed at measuring the absolute CMB brightness have been conducted before COBE (Johnson & Wilkinson 1987; Levin et al. 1992) and after it (Bersanelli et al. 1994; Fixsen et al. 2011, 2004; Raghunathan & Subrahmnayan 2000; Staggs et al. 1996; Zannoni et al. 2008). The best limits on CMB μtype spectral distortions from these experiments are at the level of 10^{4}−10^{5} (Fixsen et al. 1996; Gervasi et al. 2008; Seiffert et al. 2011). Achieving a sensitivity of 10^{9} would require subtracting foregrounds due to synchrotron emission, free free emission, dust emission, and spinning dust emission at the same precision level. The proposed experiment Pixie (Kogut et al. 2011) aims to achieve this goal by using 400 effective channels between frequencies of 30 GHz to 6 THz. Their simulations indicate that an accuracy of 1nK in foreground subtraction is achievable. There is, however, some uncertainty in our understanding of the foregrounds and possible systematics as indicated by an observed excess signal at 3 GHz by ARCADE (Fixsen et al. 2011) experiment, which is not completely explained by the current galactic and extragalactic emission models; nevertheless achievement of a goal of 10^{9} in sensitivity, although challenging, seems possible in the near future.
8. Conclusions
Comptonization of CMB photons with slightly cooler electrons results in spectral distortions that then evolve (in the absence of any other physical process) towards the BoseEinstein condensate solution, which is the equilibrium solution and consists of a blackbody spectrum with zero chemical potential, with the extra photons accumulating at zero frequencies (Illarionov & Sunyaev 1975a; Landau & Lifshitz 1980). This is because the radiation has more photons than can be accommodated in a BoseEinstein spectrum with nonnegative chemical potential at the electron temperature. However, in reality, bremsstrahlung and double Compton scattering destroy extra photons at low frequencies. Thus in practice we have the highfrequency spectrum evolving towards a blackbody with the extra photons slowly moving down in frequency to be eaten up by bremsstrahlung and double Compton scattering. This is what is seen in the numerical solution in Fig. 6.
The difficulty of conserving photon number is the reason BoseEinstein condensation (or accumulation of photons at low frequencies) is not observed in astrophysical systems. BoseEinstein condensation of photons has only been achieved in the laboratory very recently by Klaers et al. (2010). In the real Universe we also have inevitable sources of energy injection into the CMB, such as the SZ effect after reionization, which creates distortions of similar magnitude but with opposite signs canceling and overwhelming this effect. The situation is much more interesting for μ distortions, which can be created only at high redshifts and thus have no foregrounds at lower redshifts. In standard cosmology, positive μ distortions are created by dissipation of acoustic waves on small scales due to photon diffusion. For the currently allowed values of cosmological parameters, there is a possibility that the positive μ distortions from photon diffusion can almost exactly cancel the negative μ distortions from BoseEinstein condensation of CMB, leading to a net distortion that is much smaller. BoseEinstein condensation can even dominate, leading to a net distortion μ ~ −3 × 10^{9} to which Silk damping contributes, but only as a small perturbation. Nevertheless, we must emphasize that the photon number conserving Comptonization, along with differences in adiabatic indices of radiation and matter, creates a unique system in the early Universe in which photons can begin to evolve towards a BoseEinstein condensate. It is a remarkable coincidence that a completely unrelated physical process, diffusion damping of primordial perturbations, can produce distortions of almost exactly the same or greater magnitude, leading to the suppression of this unique effect in astrophysics. We thus arrive at an intriguing conclusion: even a null result, nondetection of μtype distortion at a sensitivity of 10^{9}, rather than just placing a upper limit, actually gives a quantitative measure of the primordial smallscale power spectrum. The importance of improving the experimental sensitivity to reach this critical value cannot be overemphasized.
This definition has a sign difference with respect to the usual definition of the chemical potential in thermodynamics (Landau & Lifshitz 1980) used in Eq. (1).
Strictly speaking, a BoseEinstein spectrum is established at y ~ few as can be seen from Fig. 4 However at y ~ 0.25−1 the spectrum is already very close to BoseEinstein at high frequencies. We use this to divide the energy release into μtype and ytype estimates. Exact results are presented elsewhere (Chluba et al. 2012).
Acknowledgments
We would like to thank Matias Zaldarriaga for important remarks on the manuscript. We would also like to thank Yacine AliHaïmoud for careful reading and comments on the manuscript.
References
 AliHaïmoud, Y., & Hirata, C. M. 2011, Phys. Rev. D, 83, 043513 [NASA ADS] [CrossRef] [Google Scholar]
 Arons, J. 1971, ApJ, 164, 437 [NASA ADS] [CrossRef] [Google Scholar]
 Bersanelli, M., Bensadoun, M., de Amici, G., et al. 1994, ApJ, 424, 517 [NASA ADS] [CrossRef] [Google Scholar]
 Burigana, C., Danese, L., & de Zotti, G. 1991, A&A, 246, 49 [NASA ADS] [Google Scholar]
 Chluba, J. 2005, Ph.D. Thesis, LMU München [Google Scholar]
 Chluba, J., & Sunyaev, R. A. 2004, A&A, 424, 389 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Chluba, J., & Sunyaev, R. A. 2006, A&A, 458, L29 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Chluba, J., & Thomas, R. M. 2011, MNRAS, 412, 748 [NASA ADS] [Google Scholar]
 Chluba, J., & Sunyaev, R. A. 2012, MNRAS, 419, 1294 [NASA ADS] [CrossRef] [Google Scholar]
 Chluba, J., Khatri, R., & Sunyaev, R. A. 2012, MNRAS, submitted [arXiv:1202.0057] [Google Scholar]
 Danese, L., & de Zotti, G. 1982, A&A, 107, 39 [NASA ADS] [Google Scholar]
 Dodelson, S. 2003, Modern cosmology Scott Dodelson (Amsterdam: Netherlands Academic Press) [Google Scholar]
 Dubrovich, V. K. 1975, Sov. Astron. Lett., 1, 196 [Google Scholar]
 Fixsen, D. J., Cheng, E. S., Gales, J. M., et al. 1996, ApJ, 473, 576 [NASA ADS] [CrossRef] [Google Scholar]
 Fixsen, D. J., Kogut, A., Levin, S., et al. 2004, ApJ, 612, 86 [NASA ADS] [CrossRef] [Google Scholar]
 Fixsen, D. J., Kogut, A., Levin, S., et al. 2011, ApJ, 734, 5 [NASA ADS] [CrossRef] [Google Scholar]
 Gervasi, M., Zannoni, M., Tartari, A., Boella, G., & Sironi, G. 2008, ApJ, 688, 24 [NASA ADS] [CrossRef] [Google Scholar]
 Hu, W., & Silk, J. 1993, Phys. Rev. D, 48, 485 [NASA ADS] [CrossRef] [Google Scholar]
 Hu, W., Scott, D., & Silk, J. 1994, ApJ, 430, L5 [NASA ADS] [CrossRef] [Google Scholar]
 Illarionov, A. F., & Syunyaev, R. A. 1972, Sov. Ast., 16, 45 [Google Scholar]
 Illarionov, A. F., & Sunyaev, R. A. 1975a, Sov. Ast., 18, 413 [Google Scholar]
 Illarionov, A. F., & Sunyaev, R. A. 1975b, Sov. Ast., 18, 691 [Google Scholar]
 Johnson, D. G., & Wilkinson, D. T. 1987, ApJ, 313, L1 [NASA ADS] [CrossRef] [Google Scholar]
 Kaiser, N. 1983, MNRAS, 202, 1169 [NASA ADS] [Google Scholar]
 Keisler, R., Reichardt, C. L., Aird, K. A., et al. 2011, ApJ, 743, 28 [NASA ADS] [CrossRef] [Google Scholar]
 Klaers, J., Schmitt, J., Vewinger, F., & Weitz, M. 2010, Nature, 468, 545 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
 Kogut, A., Fixsen, D. J., Chuss, D. T., et al. 2011, JCAP, 07, 025 [NASA ADS] [CrossRef] [Google Scholar]
 Komatsu, E., Smith, K. M., Dunkley, J., et al. 2011, ApJS, 192, 18 [NASA ADS] [CrossRef] [Google Scholar]
 Kompaneets, A. S. 1956, Zh. Eksp. Teor. Fiz., 31, 876 [Google Scholar]
 Landau, L. D., & Lifshitz, E. M. 1980, Statistical physics, Part. 1 (Oxford: ButterworthHeinemann) [Google Scholar]
 Lehners, J.L., McFadden, P., Turok, N., & Steinhardt, P. J. 2007, Phys. Rev. D, 76, 103501 [NASA ADS] [CrossRef] [Google Scholar]
 Levich, E. V., & Sunyaev, R. A. 1971, Sov. Ast., 15, 363 [Google Scholar]
 Levin, S., Bensadoun, M., Bersanelli, M., et al. 1992, ApJ, 396, 3 [NASA ADS] [CrossRef] [Google Scholar]
 Lifshitz, E. M., & Pitaevskii, L. P. 1981, Physical Kinetics (Oxford: ButterworthHeinemann) [Google Scholar]
 Ma, C.P., & Bertschinger, E. 1995, ApJ, 455, 7 [NASA ADS] [CrossRef] [Google Scholar]
 Mukhanov, V. F., & Chibisov, G. V. 1981, Sov. J. Exp. Theor. Phys. Lett., 33, 532 [Google Scholar]
 Peebles, P. J. E. 1968, ApJ, 153, 1 [NASA ADS] [CrossRef] [Google Scholar]
 Raghunathan, A., & Subrahmnayan, R. 2000, JApA, 20, 1 [NASA ADS] [Google Scholar]
 RubiñoMartín, J. A., Chluba, J., & Sunyaev, R. A. 2006, MNRAS, 371, 1939 [NASA ADS] [CrossRef] [Google Scholar]
 RubiñoMartín, J. A., Chluba, J., & Sunyaev, R. A. 2008, A&A, 485, 377 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Seager, S., Sasselov, D. D., & Scott, D. 2000, ApJS, 128, 407 [NASA ADS] [CrossRef] [Google Scholar]
 Seiffert, M., Fixsen, D. J., Kogut, A., et al. 2011, ApJ, 734, 6 [NASA ADS] [CrossRef] [Google Scholar]
 Silk, J. 1968, ApJ, 151, 459 [NASA ADS] [CrossRef] [Google Scholar]
 Staggs, S. T., Jarosik, N. C., Wilkinson, D. T., & Wollack, E. J. 1996, ApJ, 458, 407 [NASA ADS] [CrossRef] [Google Scholar]
 Syunyaev, R. A. 1971, Sov. Ast., 15, 190 [Google Scholar]
 Sunyaev, R. A., & Chluba, J. 2009, Astron. Nachr., 330, 657 [NASA ADS] [CrossRef] [Google Scholar]
 Sunyaev, R. A., & Zeldovich, Y. B. 1970a, Ap&SS, 9, 368 [NASA ADS] [CrossRef] [Google Scholar]
 Sunyaev, R. A., & Zeldovich, Y. B. 1970b, Ap&SS, 7, 20 [Google Scholar]
 Weinberg, S. 2008, Cosmology (Oxford: Oxford University Press) [Google Scholar]
 Zannoni, M., Tartari, A., Gervasi, M., et al. 2008, ApJ, 688, 12 [NASA ADS] [CrossRef] [Google Scholar]
 Zel’Dovich, Y. B., & Levich, E. V. 1969, Sov. J. Exp. Theor. Phys., 28, 1287 [NASA ADS] [Google Scholar]
 Zel’Dovich, Y. B., & Levich, E. V. 1970, Sov. J. Exp. Theor. Phys. Lett., 11, 35 [NASA ADS] [Google Scholar]
 Zeldovich, Y. B., & Sunyaev, R. A. 1969, Ap&SS, 4, 301 [NASA ADS] [CrossRef] [Google Scholar]
 Zeldovich, Y. B., Kurt, V. G., & Syunyaev, R. A. 1968, Zh. Eksp. Teor. Fiz., 55, 278 [Google Scholar]
 Zeldovich, Y. B., Illarionov, A. F., & Syunyaev, R. A. 1972, Sov. J. Exp. Theor. Phys., 35, 643 [NASA ADS] [Google Scholar]
Appendix A: Erasure of perturbations due to free streaming
We can derive the solution for erasure if perturbations on horizon scales due to free streaming as follows. We start with the Boltzmann hierarchy, ignoring quadrupole and higher order moments. We define multipole (ℓ) moments of the temperature anisotropy , where is the Legendre polynomial, and the angle between photon momentum and comoving wavenumber. In the limit of zero Thomson optical depth we get where η = dt/a is the conformal time. Using the fact that in the matterdominated era the potential is suppresses by a factor of 9/10 compared to primordial value on large scales, we have Θ_{0} = − 2ψ/3 (see for example Dodelson 2003) and looking for a solution with time dependence of type Θ_{ℓ} ∝ e^{i∫dηω} we get ω = ik/ Thus free streaming damps the perturbations by a factor of . Finally we want to mention that this solution is very approximate, and in principle we should also take higher ℓ modes into account. For example, on going up to ℓ = 3 we get for the damping factor . Also for free streaming we have super position of blackbodies, and the y distortion is given by (Chluba & Sunyaev 2004) , and so the equivalent energy release is .
All Tables
Energy injection in μ distortions during 5 × 10^{4} < z < 2 × 10^{6} for different initial power spectra without running compared with energy losses due to BoseEinstein condensation.
Energy injection in μ distortions during 5 × 10^{4} < z < 2 × 10^{6} for different initial power spectra with running compared with energy losses due to BoseEinstein condensation.
All Figures
Fig. 1 Difference in intensity from the blackbody radiation for normal SZ effect with Y_{SZ} = 10^{9} and negative SZ effect with Y_{BEC} = 10^{9}. 

In the text 
Fig. 2 Difference in intensity from the blackbody radiation for μ distortion defined by the equation with n(ν) = 1/(e^{hν/kBTμ + μ} − 1). T_{μ} is the temperature to which the spectrum relaxes for an initial energy addition/loss of ΔE/E with ΔT_{μ}/T = 0.64ΔE/E = ± 2.5 × 10^{9}. ΔE/E = ± 4 × 10^{9} is the energy addition/loss that gives rise to SZ distortion of 10^{9}. The dimensionless chemical potential μ is given by μ = 2.2ΔT/T = ± 5.6 × 10^{9}. 

In the text 
Fig. 3 Comparison of positive and negative μtype distortions expected in the early Universe (to be calculated in later sections) with the ytype distortions before reionization and the cosmological recombination spectrum (taken from Chluba & Sunyaev 2006; RubiñoMartín et al. 2008) for illustration. 

In the text 
Fig. 4 Evolution of initial spectrum at y = 0 given by Eq. (18)with − Y_{SZ} = Y_{BEC} = 10^{5}. ΔT/T ≡ T_{eff}/T − 1 is plotted, where T_{eff} is the temperature of a blackbody spectrum with the corresponding intensity at a particular frequency. BoseEinstein spectrum defined by n(ν) = 1/(e^{hν/kBTe + μ} − 1) with negative μ is labeled “ve BE”, red dashed, doubledotted line in the figure. The distortions are positive in the RayleighJeans region (labeled “+ve”) on the left side and negative in the Wien region (labeled “ve”) on the right side of the figure. 

In the text 
Fig. 5 Regions of the blackbody photosphere probed by different experiments, COBE FIRAS with sensitivity μ = 10^{5}, Pixie with sensitivity μ = 10^{8} and a hypothetical experiment with sensitivity of μ = 10^{10}. Shaded regions below the curves are allowed and inaccessible to these experiments. We also plot the ratio of final μ distortion today to the initial μ distortion as a function of energy injection redshift. We define the blackbody surface as the energy injection redshift such that μ^{final}/μ^{initial} = 1/e, z_{bb} = 1.98 × 10^{6}. 

In the text 
Fig. 6 Evolution of spectral distortions in CMB in ΛCDM cosmology. The values of total Compton parameter y = y(100,z) for different curves are (z = 1.5 × 10^{6}, y = 109), (z = 4.6 × 10^{5}, y = 10), (z = 1.5 × 10^{5}, y = 1), (z = 5.9 × 10^{4}, y = 0.16), (z = 1.2 × 10^{4}, y = 0.005), (z = 1760,y = 3 × 10^{5}), (z = 840, y = 1.8 × 10^{8}) and (z = 220, y = 1.2 × 10^{10}). Reference spectrum is the blackbody CMB spectrum at z = 10^{7}. At low frequencies the spectrum is dominated by double Compton and bremsstrahlung which maintain Planck spectrum at the electron temperature. The zero point at low frequencies moves to the right at low redshifts because bremsstrahlung dominates over the Compton scattering at higher and higher frequencies and absorbs photons, although both are slower than the expansion rate. 

In the text 
Fig. 7 Sketch of fractional rate of energy release due to Silk damping and free streaming for different initial power spectra. Also shown for comparison is the rate of energy loss due to adiabatic cooling of baryonic matter. 

In the text 
Fig. 8 μ = 1.4ΔE/E as a function of spectral index n_{s} without running. 

In the text 
Fig. 9 μ = 1.4ΔE/E as a function of running of the spectral index r for three values of spectral index n_{s}. 

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.