 Solution

Write $n=n_0+y\mathcal{L}{n}_0+O(y^2)$, with $n_0=(e^x-1{)}^{-1}$ and $\mathcal{L}=x^{-2}{\partial }_x(x^4{\partial }_x)$. Differentiating the Planck photon distribution givesDividing by $n_0$ proves the first-order thermal Sunyaev-Zeldovich effect:In the Rayleigh-Jeans law limit, $x{e}^x/(e^x-1)\to 1$ and $x\coth (x/2)\to 2$, so $\Delta n/n_0\to -2y$. In the Wien approximation, the prefactor approaches $x$ and the bracket approaches $x-4$, so $\Delta n/n_0\sim yx(x-4)\sim y{x}^2$. The last limit is an asymptotic description of the first-order coefficient: a perturbative prediction at large $x$ still requires $y{x}^2\ll 1$.

Hot electrons transfer energy to the cosmic microwave background, removing low-frequency photons and enhancing the high-frequency tail. The first-order distortion crosses zero where $x\coth (x/2)=4$, at $x\simeq 3.830$; for today's CMB this is near $217\mathrm{GHz}$. Photon conservation follows directly from the boundary termIn contrast, integration by parts gives $\int x^3\Delta ndx=4y\int x^3{n}_{0}dx$, so radiation gains a fraction $4y$ of its original energy. The thermal Sunyaev-Zeldovich effect is thus a spectral distortion, not simply a new blackbody temperature.

A galaxy cluster with bulk peculiar velocity along the line of sight additionally produces the kinetic Sunyaev-Zeldovich effect. Define $v_{\parallel }>0$ for recession and $\tau =\int n_e{\sigma }_{T}d\ell$. To first order in velocity and optical depth, $\Delta T/T=-\tau v_{\parallel }/c$: receding gas gives a decrement and approaching gas an increment. Its occupation number perturbation is $(\Delta T/T)x{e}^x/(e^x-1{)}^2$, the spectrum of a small blackbody temperature change. It remains nonzero at the thermal null and therefore complicates measurements near that frequency.