Cosmic Microwave Background Spectral Distortions

Cosmic microwave background spectral distortions are small, frequency-dependent departures of the sky-averaged or direction-dependent cosmic microwave background (CMB) photon energy spectrum from an equilibrium Planckian blackbody distribution B_ν(T₀) = (2hν³/c²) / [exp(hν/(k_BT₀)) − 1] at present-day temperature T₀ = 2.72548 ± 0.00057 K. Formulated theoretically by Rashid A. Sunyaev and Yakov B. Zel'dovich (1969; 1970), Luigi Danese and Gianfranco De Zotti (1977; 1982), Wayne Hu and Joseph Silk (1993), and Jens Chluba and Rashid A. Sunyaev (2012), CMB spectral distortions encode a thermal fossil record of energetic processes occurring both before and after cosmological recombination (z ≈ 1,090)—spanning redshifts up to z ≈ 2 × 10⁶ (roughly two months after the Big Bang, when the cosmic plasma temperature was k_BT ≈ 0.5 keV). Governed by the Kompaneets equation coupled to photon-creating radiative processes—double Compton scattering (γ + e⁻ ↔ e⁻ + γ + γ) and thermal bremsstrahlung (e⁻ + p ↔ e⁻ + p + γ)—any injection or extraction of energy Δρ_γ/ρ_γ after z_th ≈ 2 × 10⁶ cannot be fully thermalized into a new blackbody because photon creation has frozen out. Instead, energy release between 5 × 10⁴ ≲ z ≲ 2 × 10⁶ drives the photon occupation number into a Bose–Einstein distribution characterized by a dimensionless chemical-potential distortion (μ-distortion), while energy release at z ≲ 5 × 10⁴ produces a Comptonized Zel'dovich–Sunyaev y-distortion (with intermediate residual i-type distortions generated during the transition epoch 1.5 × 10⁴ ≲ z ≲ 2 × 10⁵).

Photon Occupation Number and the Blackbody Baseline

In a homogeneous, isotropic Friedmann–Lemaître–Robertson–Walker (FLRW) universe, the phase-space distribution of the cosmic microwave background radiation field at frequency ν and temperature T is described by the dimensionless photon occupation number n(ν, t), related to specific intensity I_ν by I_ν = (2hν³/c²) n(ν, t). Defining the dimensionless redshift-invariant frequency variable:

an ideal Planckian blackbody spectrum corresponds to the Bose–Einstein occupation number with zero chemical potential (μ = 0):

Because both photon frequency ν ∝ (1 + z) and reference temperature T_z ∝ (1 + z) scale identically under cosmological expansion, the ratio x = hν/(k_BT_z) and the Planckian occupation number n_Pl(x) remain invariant in the absence of energy exchange or photon injection.

The Kompaneets Equation Plus Photon Creation Terms

When energy is injected into the baryon–photon plasma—creating a temperature difference between free electrons (T_e) and photons (T_γ)—the time evolution of the photon occupation number n(x, t) is governed by the generalized Kompaneets equation (Kompaneets, 1956; Weymann, 1965; Danese and De Zotti, 1977) supplemented by double Compton and bremsstrahlung emission and absorption terms:

where n_e is the free electron number density, σ_T = 6.652 × 10⁻²⁵ cm² is the Thomson cross section, m_e is the electron rest mass, and x_e = hν/(k_BT_e).

The three physical processes in this radiative transport equation play distinct roles in establishing thermodynamic equilibrium:

1. Elastic Compton Scattering (γ + e⁻ ↔ γ + e⁻): Described by the first diffusion term on the right-hand side, Compton scattering redistributes photon energies across frequency space via Doppler shifts (Δν/ν ~ (k_BT_e/m_ec²)^1/2) and electron recoil, driving n(x) toward a Bose–Einstein distribution n_BE(x_e) = 1 / [exp(x_e + μ) − 1] at the electron temperature T_e. Elastic Compton scattering conserves total photon number (∫ x² n(x) dx = const). If pure thermal energy ΔQ is added to the plasma without adding a matching number of photons (ΔN_γ/N_γ = (3/4) Δρ_γ/ρ_γ), Compton scattering alone cannot restore a Planck blackbody (μ = 0); it leaves a positive chemical potential μ > 0.
2. Double Compton Scattering (γ + e⁻ ↔ e⁻ + γ + γ): A relativistic radiative correction to Compton scattering whose cross section scales as σ_DC ~ (4α_em/3π) (k_BT/m_ec²)² σ_T, where α_em ≈ 1/137.036 is the fine-structure constant. Because the radiation density exceeds the baryon density by the photon-to-baryon ratio η⁻¹ = n_γ/n_b ≈ 1.6 × 10⁹, double Compton scattering dominates over bremsstrahlung as the primary photon-creation mechanism at high redshifts z ≳ 10⁵ (Lightman, 1981; Thorne, 1981).
3. Thermal Bremsstrahlung (Free–Free Emission/Absorption, e⁻ + p ↔ e⁻ + p + γ): Scales with the square of the baryon density (n_en_p ∝ (1 + z)⁶ T_e^−7/2) and dominates photon creation and absorption at lower redshifts z ≲ 10⁵ and very low frequencies x ≪ 1.

Both double Compton scattering and bremsstrahlung emit and absorb low-frequency (x ≪ 1) photons at a rate proportional to x⁻³ [1 − n(x)(exp(x_e) − 1)].

Depending on the redshift z at which energy is injected into the cosmic plasma, the competition between the Hubble expansion rate H(z), the Compton energy-exchange rate t_C⁻¹ = (4k_BT_e/m_ec²) n_eσ_Tc, and the double-Compton photon-production rate divides cosmic history into three canonical eras.

1. The Blackbody / Complete Thermalization Era (z ≳ 2 × 10⁶)

At redshifts z > z_th ≈ 1.98 × 10⁶ (cosmic time t ≲ 7 × 10⁶ s, T ≳ 5.4 × 10⁶ K), double Compton scattering at low frequencies (x ≪ 1) produces new soft photons while elastic Compton scattering upscatters them across the thermal peak (x ≈ 2.82). The blackbody visibility function 𝒢_th(z) formulated by Sunyaev and Zel'dovich (1970) and Danese and De Zotti (1982) equals:

For any energy injection at z ≫ z_th, the suppression factor exp[−(z/z_th)^5/2] erases any spectral distortion, relaxing the radiation field to a Planckian blackbody B_ν(T′) at shifted temperature T′ = T[1 + (1/4)(Δρ_γ/ρ_γ)].

2. The Chemical-Potential / μ-Distortion Era (5 × 10⁴ ≲ z ≲ 2 × 10⁶)

Between z_μ,y ≈ 5 × 10⁴ and z_th ≈ 2 × 10⁶, double Compton and bremsstrahlung processes are too slow compared to Hubble expansion to restore the photon-to-energy ratio across the thermal peak, except at low frequencies x ≤ x_c ≈ 5 × 10⁻³. However, the elastic Comptonization parameter:

remains much greater than unity (y_γ ≫ 1). Elastic Compton scattering establishes kinetic equilibrium between photons and electrons, driving the photon occupation number into a Bose–Einstein distribution n_BE(x) = 1 / [exp(x + μ(x)) − 1], where μ(x) = μ₀ exp(−x_c/x).

For a fractional energy release Δρ_γ/ρ_γ ≪ 1 during the μ-era, conservation of photon number and energy across the Bose–Einstein spectrum fixes the chemical potential amplitude (Sunyaev and Zel'dovich, 1970; Chluba and Sunyaev, 2012):

Expanding n_BE(x) − n_Pl(x) to first order in μ₀ at constant photon number (ΔT/T = μ₀/2.1923) yields the μ-distortion intensity shape:

and I₀ = (2h/c²)(k_BT₀/h)³ ≈ 270.33 MJy sr⁻¹. A positive μ-distortion (μ₀ > 0) produces a photon deficit below ν_0,μ ≈ 124.5 GHz (x < 2.1923) and an excess above 124.5 GHz.

3. The Compton / y-Distortion Era (z ≲ 5 × 10⁴) and Intermediate i-Distortions

At redshifts z ≲ z_μ,y ≈ 5 × 10⁴, the Comptonization parameter drops below unity (y_γ ≲ 1). Heated electrons partially upscatter low-frequency CMB photons toward higher frequencies in the small-y diffusion regime of the Kompaneets equation (Zel'dovich and Sunyaev, 1969), yielding the Zel'dovich–Sunyaev y-type spectral distortion:

with dimensionless Compton parameter y = ∫ (k_B(T_e − T_γ)/m_ec²) n_eσ_Tc dt = (1/4)(Δρ_γ/ρ_γ)_y. A positive y-distortion shifts photons from the Rayleigh–Jeans region (ΔI_ν < 0 for ν < 217.6 GHz) into the Wien tail (ΔI_ν > 0 for ν > 217.6 GHz), with a null at ν_0,y = 217.6 GHz (x = 3.830). For energy injected during the transition epoch 1.5 × 10⁴ ≲ z ≲ 2 × 10⁵ (0.01 ≲ y_γ ≲ 10), partial Comptonization produces intermediate (i-type) residual distortions (Burigana, Danese, and De Zotti, 1991; Khatri and Sunyaev, 2012; Chluba and Sunyaev, 2012) that cannot be expressed as a sum of M(x) and Y_SZ(x), encoding the injection redshift z_inj.

Even in the standard ΛCDM cosmological model with zero exotic particle physics, four physical mechanisms generate CMB spectral distortions across amplitudes spanning y ~ 10⁻⁶ down to |μ| ~ 10⁻⁹:

1. Dissipation of Small-Scale Acoustic Waves (Silk Damping, μ ≈ 2 × 10⁻⁸): Prior to recombination, primordial adiabatic density perturbations on comoving scales 1 Mpc⁻¹ ≲ k ≲ 10⁴ Mpc⁻¹ enter the horizon and oscillate as baryon acoustic oscillations. Shear viscosity and thermal conduction in the photon–baryon fluid damp these sound waves (Silk, 1968), converting acoustic energy into isotropic heat (Sunyaev and Zel'dovich, 1970; Hu, Scott, and Silk, 1994; Chluba, Khatri, and Sunyaev, 2012). For a Planck-normalized scalar power spectrum (A_s = 2.1 × 10⁻⁹, n_s = 0.965), acoustic dissipation during 5 × 10⁴ ≲ z ≲ 2 × 10⁶ generates a chemical potential μ_ac ≈ 2.0 × 10⁻⁸, probing the primordial power spectrum at k ≈ 50–10⁴ Mpc⁻¹.
2. Adiabatic Cooling of Baryonic Matter (μ ≈ −2.7 × 10⁻⁹, y ≈ −5 × 10⁻¹⁰): Non-relativistic baryons and electrons cool with expansion as T_b ∝ (1 + z)², whereas relativistic CMB photons cool as T_γ ∝ (1 + z). Keeping electrons and ions thermally coupled to the radiation field (T_e ≈ T_γ) continuously extracts heat from the photon bath into the baryons (Chluba and Sunyaev, 2012), producing a negative distortion: μ_cool ≈ −2.7 × 10⁻⁹ and y_cool ≈ −5 × 10⁻¹⁰.
3. Cosmological Recombination Radiation (CRR, ΔI_ν/I₀ ~ 10⁻⁹): During helium II→I (z ≈ 6,000), helium I (z ≈ 2,000), and hydrogen recombination (z ≈ 1,090), atomic transitions (Lyman-α, Balmer Hα, Paschen, Brackett, and two-photon 2s→1s decays) inject ~5 non-thermal photons per baryon into the CMB across 1–1,000 GHz (Dubrovich, 1975; Rubiño-Martín, Chluba, and Sunyaev, 2006).
4. Reionization and Cosmic Structure Formation (y_reion+WHIM ≈ 1.77 × 10⁻⁶): At late times (z ≲ 15), ultraviolet radiation from stars and quasars photoionizes the intergalactic medium (T_e ~ 10⁴ K, y_reion ≈ 2 × 10⁻⁷), while shock heating inside galaxy groups, clusters, and the warm–hot intergalactic medium (WHIM, T_e ~ 10⁶–10⁸ K) generates a sky-averaged thermal Sunyaev–Zel'dovich effect of ⟨y⟩ ≈ (1.58–1.77) × 10⁻⁶ (Refregier et al., 2000; Hill et al., 2015).

The table below summarizes the physical epochs, spectral signatures, and observational constraints of the principal CMB spectral distortion modes.

The COBE/FIRAS Benchmark (1990–1996)

The absolute measurement of the CMB monopole spectrum was performed by the Far Infrared Absolute Spectrophotometer (FIRAS) aboard NASA's COBE satellite (Mather et al., 1990; Fixsen et al., 1996). Using a cryogenic Martin–Puplett polarizing Fourier-transform Michelson interferometer comparing the sky against an external reference blackbody calibrator (XCAL) across 68–2,910 GHz, Fixsen et al. (1996) established that the CMB spectrum matches a blackbody within 50 parts per million of peak brightness, yielding the 95% confidence upper limits:

Joint analyses combining COBE/FIRAS with ARCADE-2 and TRIS radiometer data tightened the chemical-potential bound to |μ| < 4.7 × 10⁻⁵ (Bianchini and Fabbian, 2022), constraining decaying relic particles with lifetimes 10⁷ s ≲ τ ≲ 10¹² s and evaporating primordial black holes in the mass window 10¹¹ g ≲ M_PBH ≲ 10¹³ g.

Next-Generation Fourier-Transform Spectrometers (PIXIE, Super-PIXIE, FOSSIL)

Detecting the ΛCDM sky-averaged Compton distortion (⟨y⟩ ≈ 1.77 × 10⁻⁶) and the Silk-damping chemical potential (μ ≈ 2.0 × 10⁻⁸) requires photon-noise-limited multimoded cryogenic Fourier-transform spectrometers operating in space:

• PIXIE (Primordial Inflation Explorer): Proposed by Kogut et al. (2011; 2024), PIXIE pairs two symmetric input ports with multimoded feedhorns cooled to 2.7 K, achieving σ(y) ≈ 2 × 10⁻⁹ and σ(μ) ≈ 10⁻⁸ across 15 GHz to 6 THz.
• Super-PIXIE and Voyage 2050 Concepts (FOSSIL / PRISM): ESA Voyage 2050 mission designs (Chluba et al., 2021) combine low-frequency (10–600 GHz) and high-frequency (300–3,000 GHz) interferometers to separate galactic synchrotron, free–free, thermal dust, and cosmic infrared background foregrounds, reaching σ(μ) ≈ (2–5) × 10⁻⁹.

Why can elastic Compton scattering never erase a μ-distortion on its own?

Elastic Compton scattering (γ + e⁻ ↔ γ + e⁻) transfers thermal energy between electrons and photons while conserving total photon number N_γ. In an equilibrium Planck blackbody B_ν(T), both energy density ρ_γ = a_RT⁴ and photon number density n_γ = b_RT³ are fixed by the single temperature T. If pure energy Δρ_γ > 0 is injected without adding photons, the radiation bath has too few photons for its new energy density. Elastic Compton scattering drives the photons to the maximum-entropy state at fixed (N_γ, ρ_γ), which is a Bose–Einstein distribution with positive chemical potential μ ≈ 1.401(Δρ_γ/ρ_γ) > 0. Erasing μ back to zero requires photon-creating processes—double Compton scattering and thermal bremsstrahlung—which freeze out below z_th ≈ 2 × 10⁶.

What is the difference between the sky-averaged y-distortion and the cluster Sunyaev–Zel'dovich effect?

Both phenomena obey the Kompaneets equation spectral function Y_SZ(x) = [x⁴e^x/(e^x − 1)²][x coth(x/2) − 4] with a null at 217.6 GHz. The cluster Sunyaev–Zel'dovich effect is an anisotropic (direction-dependent) y-distortion (y ~ 10⁻⁵–10⁻⁴) measured along lines of sight through individual galaxy clusters by differential CMB imagers (Planck, ACT, SPT) that subtract adjacent sky patches. The monopole y-distortion is the all-sky average ⟨y⟩ ≈ 1.77 × 10⁻⁶ integrated across all cosmic epochs (z ≲ 5 × 10⁴) and all diffuse ionized gas, measured by an absolute spectrometer comparing the sky to an internal blackbody reference load.

How do μ-distortions constrain primordial black holes and small-scale inflation?

Primary CMB temperature anisotropies measured by Planck probe comoving scales 10⁻⁴ Mpc⁻¹ ≲ k ≲ 0.3 Mpc⁻¹ because smaller scales are erased by Silk diffusion damping during recombination. However, the acoustic energy dissipated by Silk damping on small scales 1 Mpc⁻¹ ≲ k ≲ 10⁴ Mpc⁻¹ between z ≈ 2 × 10⁶ and z ≈ 5 × 10⁴ is converted into a chemical-potential distortion μ_ac ≈ 1.4 ∫ P_ℛ(k) W_μ(k) d ln k. Upper limits on μ therefore constrain any blue-tilted primordial power spectrum or density spikes that would collapse into stellar-mass or intermediate-mass primordial black holes.

1. Sunyaev, R. A., & Zel'dovich, Y. B. (1970). The interaction of matter and radiation in the hot model of the Universe, II. Astrophysics and Space Science, 7(1), 20–30. https://doi.org/10.1007/BF00653472
2. Danese, L., & De Zotti, G. (1982). Double Compton process and the spectrum of the microwave background. Astronomy and Astrophysics, 107(1), 39–42. https://ui.adsabs.harvard.edu/abs/1982A&A...107...39D
3. Hu, W., & Silk, J. (1993). Thermalization and spectral distortions of the cosmic background radiation. Physical Review D, 48(2), 485–502. https://doi.org/10.1103/PhysRevD.48.485
4. Fixsen, D. J., Cheng, E. S., Gales, J. M., Mather, J. C., Shafer, R. A., & Wright, E. L. (1996). The cosmic microwave background spectrum from the full COBE FIRAS data set. The Astrophysical Journal, 473(2), 576–587. https://doi.org/10.1086/178173
5. Chluba, J., & Sunyaev, R. A. (2012). The evolution of cosmic microwave background spectral distortions in the early Universe. Monthly Notices of the Royal Astronomical Society, 419(2), 1294–1314. https://doi.org/10.1111/j.1365-2966.2011.19786.x