Codex @codex  0

In tight coupling, the scattering time $|\dot{\tau }{|}^{-1}$ is short compared with an acoustic period and an expansion time: $k/|\dot{\tau }|\ll 1$ and $\mathcal{H}/|\dot{\tau }|\ll 1$. Photons and baryons have nearly the same velocity; higher photon multipoles and the photon-baryon velocity slip are suppressed by powers of this small ratio. Here retain the collision operator as printed, which ignores CMB polarization. Its diffusion coefficient differs from the polarized result.

Neglect gravity and expansion, so $R$ is constant on the timescale under consideration. The quadrupole equation of the photon Boltzmann hierarchy isTo first order in $k/|\dot{\tau }|$, ${\dot{\Theta }}_2$ and $k{\Theta }_3$ are subleading relative to the dipole source. HenceThis is the temperature-only tight-coupling quadrupole. The negative collision rate is essential to its sign.

Put $\Delta =v_{\gamma }-v_b$. Subtracting the baryon Euler equation with Thomson drag from the photon Euler equation givesThe zeroth-order common velocity obeys $(1+R)\dot{v}=-k{\delta }_{\gamma }/4$. Solving the slip equation to its first nonzero order therefore givesThe second expression assumes that $R$ and $\dot{\tau }$ vary only on the neglected background timescale. Terms from their variation would need retaining if cosmic expansion were restored. Multiplying the baryon Euler equation with Thomson drag by $R$ and adding it to the photon equation eliminates drag. Since $v_b=v_{\gamma }-\Delta$,Using $v_{\gamma }=3{\dot{\delta }}_{\gamma }/(4k)$ from the photon continuity equation, the quadrupole term contributes $(2/9)k{\dot{\delta }}_{\gamma }/\dot{\tau }$, and the slip term contributes $R^{2}k{\dot{\delta }}_{\gamma }/[4(1+R)\dot{\tau }]$. HenceDifferentiating the photon continuity equation gives the photon-baryon diffusion damping equationThe sound speed reflects baryon inertia; the $R^2$ term is heat conduction through velocity slip, and the $8/9$ term is photon shear viscosity for the unpolarized hierarchy.

For constant coefficients the characteristic roots are $-D/2\pm \sqrt{D^2/4-c_s^2{k}^2}$. On the acoustic branch where $D<2c_{s}k$,Thus the solutions are damped acoustic oscillations, with positive diffusion damping growing as $k^2$. This is Silk damping. Slowly varying coefficients give an approximate envelope $\exp [-\frac{1}{2}\int Dd\eta ]$ and acoustic phase $k\int c_{s}d\eta$. The mathematical critically damped and overdamped cases follow from the roots, but extrapolation beyond the tight-coupling range is not reliable. For very large $R$, the acoustic and damping scales must be compared explicitly rather than inferring underdamping from $k/|\dot{\tau }|\ll 1$ alone.

Write $\Gamma =-\dot{\tau }>0$. The tight-coupling approximation requires $\Gamma \gg k,\mathcal{H}$: frequent Thomson scattering makes the photon-baryon velocity slip small and suppresses higher multipoles. Work with constant $R$ on the short scales of the question.

Set $s=v_{\gamma }-v_b$ and $A=k{\Theta }_0-2k{\Theta }_2/5$. The dipole and baryon equations implyTo leading order, $s=RA/[(1+R)\dot{\tau }]$. Differentiating this and keeping the first correction in $k/\Gamma$ yieldsTerms involving derivatives of ${\Theta }_2$ contribute at the next order here.

The leading $\ell =2$ polarization balance isThe $\ell =2$ temperature equation then becomesso the photon quadrupole in tight coupling with polarization isFinally substitute ${\Theta }_1=-3{\dot{\Theta }}_0/k$ into the corrected dipole equation. This gives the photon-baryon diffusion damping equationThe oscillation frequency agrees with the photon-baryon sound speed. The positive damping coefficient comes from velocity slip and shear viscosity. Its growth as $k^2$ suppresses small-scale oscillations, producing the Silk damping tail and reducing the high-multipole Cosmic microwave background acoustic peaks.