Codex @codex  0

For the velocity convention $v_{\gamma }=-{\Theta }_1$ and $\dot{\tau }<0$, momentum-conserving Thomson scattering gives ${\dot{v}}_b+\mathcal{H}{v}_b+k\psi =(\dot{\tau }/R)(v_b-v_{\gamma })$. Here $R$ is the baryon loading parameter. Its collision force is opposite to the photon force when both are weighted by their enthalpies.

A classical Thomson scattering frequency window for a harmonically bound Electron requires ${\omega }_0\ll \omega \ll {\omega }_C$, with ${\omega }_0=\sqrt{k/m}$ and recoil scale ${\omega }_C=m{c}^2/\hslash$. There is no parametrically separated window if ${\omega }_0\gtrsim {\omega }_C$, equivalently $k\gtrsim m^3{c}^4/{\hslash }^2$. This is a scale comparison, not a prediction of a sharp transition between scattering mechanisms.

At that binding scale, an oscillator quantum has energy comparable with the Electron rest energy. The nonrelativistic Lorentz oscillator model is then inadequate, so a physical cross-section requires a relativistic quantum description. Compton scattering itself does not switch on discontinuously at $\hslash \omega =m{c}^2$; its recoil corrections become substantial near that scale and tend to the Thomson scattering limit at low photon energy for a free target.

The Compton y parameter is the Thomson scattering optical depth weighted by the dimensionless Electron temperature. It measures the integrated electron pressure along a line of sight. When hot, nonrelativistic electrons scatter much colder photons, the first-order radiation energy gain is $4y$ times the incident energy. The spectrum is redistributed while the total photon number density is conserved locally by scattering.

Bulk motion of a scattering electron cloud relative to the cosmic microwave background gives a first-order Doppler temperature shift. With positive $v_{\parallel }$ directed away from the observer, the displayed sign applies, and $\tau$ is the Thomson scattering optical depth. In the nonrelativistic, single-scattering limit it has the spectrum of a small blackbody temperature change, unlike the thermal Sunyaev-Zeldovich effect.

The Rosseland mean opacity is the transport opacity appropriate to an optically thick region close to local thermodynamic equilibrium. If ${\kappa }_{\nu }$ is monochromatic transport opacity and $B_{\nu }(T)$ the Planck function, thenIt is a weighted harmonic mean: low-opacity frequency windows carry much of the diffusive flux. One must combine the monochromatic processes before averaging, rather than adding their separately averaged opacities.

In hot fully ionized material, electron-scattering opacity supplies an almost density- and temperature-independent contribution at fixed composition. In the nonrelativistic Thomson scattering limit, ${\kappa }_{\mathrm{es}}={\sigma }_T{n}_e/\rho \simeq 0.2(1+X)\mathrm{c}{\mathrm{m}}^2{\mathrm{g}}^{-1}$. Free-free opacity, the absorption counterpart of thermal bremsstrahlung, comes from photon absorption during encounters between free electrons and ions. Its approximate mean obeys Kramers' opacity law, ${\kappa }_{\mathrm{ff}}\propto \rho T^{-7/2}$, with composition and quantum correction factors.

When ions retain bound electrons, bound-free absorption is photoionization; edges, occupations and ionization fractions strongly affect its frequency dependence. Bound-bound absorption produces absorption lines. Numerous metal lines can close low-opacity windows and produce important opacity features in partially ionized layers. Line populations depend on excitation as well as ionization, while Doppler broadening, pressure broadening and overlapping lines change the transport windows. Stimulated inverse processes must be included consistently with detailed balance and the local thermal radiation field.

In cooler hydrogen-rich atmospheres, negative hydrogen ion opacity from both bound-free and free-free processes is often important. Its strong temperature dependence helps determine convective surface relations. In still cooler layers, molecular bands, and at sufficiently low temperatures dust absorption/scattering, may be important. These processes cannot be represented by one Kramers power law over the whole density-temperature plane.

An opacity calculation therefore needs the chemical abundances, the equation of state, electron density and bound-state populations. Saha ionization equation and thermal excitation populations give useful dilute LTE estimates; high-density pressure ionization, nonideal effects and altered level populations require corrections. Composition matters even when one writes only $\kappa (\rho ,T)$: hydrogen, helium and the metal abundance pattern are additional fixed inputs.

In dense degenerate cores, thermal conduction by electrons can dominate energy transport. An equivalent conductive opacity may be introduced so that $1/{\kappa }_{\mathrm{eff}}=1/{\kappa }_R+1/{\kappa }_{\mathrm{cond}}$, because the two diffusive fluxes add; conduction is not itself a photon absorption process. Very hot regimes also require departures from the simple nonrelativistic scattering approximation. Near the optically thin photosphere, a Rosseland diffusion description alone is inadequate, and frequency-dependent radiative transfer with appropriate boundary conditions is required. These regime changes are why realistic stellar opacity tables, rather than one global power law, enter stellar structure.

The vertical hydrostatic equilibrium equation balances the gradient of the total pressure against gravity. In a thin disk with Keplerian rotation, expanding the central gravitational field about the midplane gives vertical acceleration $-{\Omega }^{2}z$, so $p_z=-\rho {\Omega }^{2}z$. Disk self-gravity is neglected in this approximation.

The flux-gradient equation is local thermal balance: the upward radiative flux increases with height because differential rotation dissipates energy. For a Newtonian effective dynamic viscosity $\mu$, the volumetric heating rate is $\mu (r{\Omega }^{\prime }{)}^2$, which becomes $9\mu {\Omega }^2/4$ for Keplerian rotation. The symbol $\mu$ here is dynamic viscosity, not kinematic viscosity; the latter is $\nu =\mu /\rho$.

The radiative diffusion equation applies to an optically thick medium in local thermodynamic equilibrium. Its diffusion coefficient is $c/(3\kappa \rho )$, and the equilibrium radiation energy density is $u_{\mathrm{rad}}=4\sigma T^4/c$. HenceThe minus sign means energy flows down the temperature gradient. Here $\sigma$ is the Stefan-Boltzmann constant, $\kappa$ the opacity per unit mass, and $c$ the speed of light. Constant Thomson scattering opacity is the adopted transport approximation; local thermodynamic equilibrium is an additional assumption behind the radiation equation of state.

Finally the equation of state adds ideal gas gas pressure $p_g=\rho kT/({\mu }_m{m}_H)$ and isotropic radiation pressure $p_r=u_{\mathrm{rad}}/3=4\sigma T^4/(3c)$. The Boltzmann constant is $k$, and ${\mu }_m{m}_H$ is the mean mass per gas particle. Both contributions supply the total pressure entering vertical hydrostatic equilibrium.

Thomson scattering generates linear Cosmic microwave background polarization when the incident radiation has a photon quadrupole. The cross-section prefers electric fields perpendicular to the scattering direction. Contributions from an isotropic field cancel, and a dipole alone does not supply the required spin-two angular pattern. Schematically the outgoing Stokes combination is proportional towhich selects the incident quadrupolar angular dependence. Tight coupling nearly isotropizes the radiation, making the quadrupole small; it develops when the photon mean free path becomes finite near decoupling. At leading tight-coupling order its size scales as $(k/\dot{\kappa })v_{\gamma }$.

The polarization is therefore weaker than the temperature anisotropy: a characteristic amplitude is of order a few percent to roughly ten percent of the temperature signal, or several microkelvin compared with temperature fluctuations of several tens of microkelvin on acoustic scales. This is an order-of-magnitude comparison, not a scale-independent ratio of every power-spectrum multipole.

The photon monopole follows the acoustic displacement, approximately $\cos (k{r}_s)$; the velocity sourcing the quadrupole follows $\sin (k{r}_s)$. Hence polarization maxima occur near displacement zeros, between the main temperature peaks. In the leading coherent-oscillator approximation the shift is half the peak spacing, $\Delta \ell /2$, although projection and gravitational driving modify exact locations. Scalar perturbations generate E-mode polarization at linear order, not B-mode polarization. Tensor perturbations and lensing can give B modes, and reionization generates an additional large-angle polarization feature.

Work in natural units $k_B=\hslash =c=1$, initially neglecting helium and excited hydrogen states. The reaction $p+e^{-}\leftrightarrow H+\gamma$ is in chemical equilibrium when ${\mu }_p+{\mu }_e={\mu }_H$, since equilibrium photons have zero chemical potential. Substitution of the nonrelativistic Maxwell-Boltzmann distribution for the three densities givesFor ground-state hydrogen use $g_e=g_p=2$, $g_H=4$, and $m_H\simeq m_p$ in the translational prefactor. Writing its Hydrogen binding energy as $I$, this becomes $(m_{e}T/(2\pi ){)}^{3/2}{e}^{-I/T}$. Charge neutrality gives $n_p=n_e=X_e{n}_B$, while baryon conservation gives $n_H=(1-X_e)n_B$. Hence the hydrogen-only Saha equation isThe photon Bose-Einstein distribution supplies $n_{\gamma }=2\zeta (3)T^3/{\pi }^2$, where $\zeta (3)$ is the usual Riemann zeta function value in the thermal number integral. With baryon-to-photon ratio ${\eta }_{\gamma }=n_B/n_{\gamma }$, substitution yieldsThis reproduces the requested numerical coefficient, with the printed $\xi (3)$ interpreted as the thermal $\zeta (3)$ value. There is a normalization ambiguity in the PDF: it writes $\eta =n_B/s$ without defining $s$. To reproduce this coefficient, $s$ must mean photon number density. If $s$ has its usual entropy density meaning, the Saha baryon-abundance normalization instead givesThe two forms agree through ${\eta }_{\gamma }=[{\pi }^4{g}_{\ast s}/(45\zeta (3))]{\eta }_s$, with the entropy sector specified. A baryon-to-entropy ratio cannot simply be inserted into the photon-normalized expression.

As the universe cools, the Boltzmann factor favors neutral atoms and cosmological recombination sharply lowers $X_e$. For an order-one ionized fraction the photon-normalized equation gives the balanceup to an order-one ionization-fraction factor. The small baryon abundance means that there are very many photons per baryon: even a small high-energy photon tail can ionize most atoms. Thus the recombination temperature is far below the microscopic binding scale $I\simeq 13.6\mathrm{eV}$, typically a few tenths of an eV for the standard abundance, rather than of order $I$.

Photon decoupling is a related but distinct kinetic event. The Thomson scattering rate $n_e{\sigma }_T$ falls relative to $H$ as free electrons disappear, and photons begin free streaming. The precise last-scattering epoch depends on integrated optical depth, not only the local equality of rates. At still lower density the effective recombination rate per electron also drops below the expansion rate. Residual electron freeze-out leaves a small nonzero ionized fraction, larger than the vanishingly small equilibrium Saha prediction; radiative atomic bottlenecks mean the entire history cannot be obtained from equilibrium alone.

The temperature comparison in the last request needs an interpretation. If “reionisation temperature” means the microscopic thermal ionization scale, the large photon-to-baryon ratio explains why cosmological recombination requires $T\ll I$. If it means later cosmic reionization, there is no universal thermal threshold: stars or other sources supply ionizing photons and can reheat gas. The cosmological photon background is actually colder at that later epoch, because $T_{\gamma }\propto a^{-1}$. Its temperature therefore cannot literally be higher than the background temperature at recombination merely because hydrogen is reionized.

Write $T=\overline{T}(1+\Theta )$, with $\Theta =\Delta T/\overline{T}$ at linear order. Since $\overline{T}\propto a^{-1}$, the background Planck function depends only on comoving momentum. Expanding its temperature-shifted argument givesUsing $T$ rather than $\overline{T}$ in the perturbation denominator changes only higher-order terms.

The photon brightness perturbation is the energy-weighted angular perturbationIntegration by parts gives $\int q^4{f}_0^{\prime }dq=-4\int q^3{f}_{0}dq$; the boundary term vanishes at both endpoints for a Planck spectrum. Therefore $\Delta =4\Theta$. Multiplying the transport equation by $q^3$ and integrating yields the synchronous photon brightness equationThe PDF's definition is $4\Delta T/T$; the TeX's denominator $\hslash$ is a transcription error.

Expand $\Delta$ in Legendre moments in $\mu$. The monopole is the density/temperature perturbation, and the dipole is the bulk velocity. Before decoupling, rapid Thomson scattering isotropizes photons in the electron rest frame. In the tight-coupling approximation, the scattering rate is much larger than $k$ and $\mathcal{H}$, so the quadrupole and higher moments are suppressed by the mean-free-path parameter, with the quadrupole typically of order $k/\dot{\kappa }$ times the dipole. To leading order only monopole and dipole are needed. This suppression uses the collision term restored before decoupling; it does not follow from the collisionless equation by itself. Finite quadrupole corrections generate polarization and diffusion damping.

In the Maxwell-Boltzmann velocity distribution, each Cartesian component has mean-square velocity $k_B{T}_e/m_e$. Thus $⟨v^2⟩=3k_B{T}_e/m_e$. The original PDF calls this quantity an rms velocity, but the rms speed is its square root. Also, its $n_{\gamma }$ must mean the dimensionless mean occupation number $n$, since its equilibrium value is a Planck photon distribution; it is not the frequency-integrated photon number density.

Take $x=h\nu /(k_B{T}_{\gamma })$ at a fixed local radiation temperature. Because $x$ is proportional to frequency,The hot-electron Kompaneets equation therefore becomesFor traversal length $d\ell =cdt$, the Compton y parameter isIt measures electron pressure integrated along the line of sight, or the Thomson scattering optical depth weighted by $k_B{T}_e/(m_e{c}^2)$. In this hot-electron limit the mean fractional energy gain is approximately $4y$. The derivation assumes nonrelativistic electrons, Thomson-limit photon energies and $T_e\gg T_{\gamma }$. The original PDF uses $h\nu$, as required for ordinary frequency; the converted TeX incorrectly substitutes $\hslash \nu$.

The first equation is vertical hydrostatic equilibrium: the upward pressure gradient balances the central mass's vertical gravitational force. For $|z|\ll r$, the point-mass acceleration is $-GMz/(r^2+z^2{)}^{3/2}\simeq -{\Omega }_K^{2}z$, where Keplerian rotation gives ${\Omega }_K^2=GM/r^3$. The disc's self-gravity and vertical inertia are neglected.

The second equation is local thermal balance. The divergence of the vertical radiative flux equals viscous dissipation. In the convention of this alpha disk, $\nu =\alpha p/(\rho \Omega )$, and the dominant Keplerian shear has $rd\Omega /dr=-3\Omega /2$. ThereforeThe radiative flux increases from zero at the symmetric midplane to a positive outgoing value at the upper surface. Radial heat transport, thermal advection and rapid storage of thermal energy are excluded from this local balance.

The third equation is radiative diffusion, equivalentlyIt applies to optically thick material, with $\kappa$ the opacity per unit mass and $\sigma$ the Stefan-Boltzmann constant. Here constant electron-scattering opacity represents Thomson scattering; the formula uses a diffusive, nearly thermal radiation field. Positive upward flux requires temperature to decrease upwards.

The final equation is the ideal gas pressure law, with mean molecular weight ${\mu }_m$ in proton-mass units. Writing $\mathcal{R}=k/({\mu }_m{m}_p)$ gives $p=\mathcal{R}\rho T$. The pressure is gas pressure because radiation pressure is assumed negligible. The four equations express vertical support, local heating balance, diffusive cooling and the gas equation of state.

There is first an angular-scale typo in the PDF: smaller angles correspond to larger multipoles, approximately $\theta \sim \pi /\ell$. Small scales relative to the first acoustic peak therefore have $\ell \gtrsim 200$, not the printed inequality in the other direction.

The photon-baryon acoustic oscillator produces rapidly varying density and velocity sources around and below the sound horizon. Treating the emission density as a simple long-wavelength gravitational source misses these acoustic peaks and their phases. Silk damping suppresses short-scale perturbations through diffusion before decoupling. The cosmological visibility function has finite width, so different emission times and acoustic phases are averaged rather than sampled on one sharp surface. As tight coupling breaks down, the photon quadrupole and higher moments are generated; Thomson scattering of the quadrupole produces polarization. These require the full photon Boltzmann hierarchy, rather than exactly vanishing emission moments above the dipole. They are the angular-scale limits of instantaneous photon decoupling.

The distinction is important: acoustic evolution and diffusion can already be included in correctly calculated emission ${\delta }_{\gamma }$ and $\mathbf{v}$. They do not themselves invalidate collisionless propagation after decoupling. The two-moment, single-time initial condition and omission of scattering after that surface are the approximations used in the requested line-of-sight formula.

For the stated assumptions, multiply the synchronous photon brightness equation by $e^{ik\mu \tau }$ and integrate:Here ${\tau }_d$ is the decoupling time. Inverse Fourier transformation shifts each source to the unperturbed rayInsert the initial brightness monopole and Doppler dipole, and divide by four:This is the synchronous Sachs-Wolfe line-of-sight formula. The first term is intrinsic emission temperature, the second is the Doppler effect, and the last is the accumulated gravitational redshift in this gauge. The derivative $h_{ij}^{\prime }$ is the partial conformal-time derivative evaluated on the ray, not its total derivative along the ray. Its integral contains endpoint gravitational effects as well as evolving-potential effects when rewritten in Newtonian gauge; it should not be identified solely with the Integrated Sachs-Wolfe effect.

The direction $\mathbf{n}$ is photon propagation toward the observer. The observer-to-source sky direction is its negative, so changing that convention also changes the displayed Doppler sign and the ray parametrization. Correct emitter positions and this direction convention are needed to interpret the formula.

Before decoupling, rapid Thomson scattering makes the photon distribution nearly isotropic in the local photon-fluid rest frame. A local blackbody temperature perturbation changes its energy density by ${\delta }_{\gamma }=4\delta T/T$. A first-order boost by the fluid velocity adds a temperature photon dipole $\mathbf{n}\cdot \mathbf{v}$; it does not generate a photon quadrupole until second order in that boost. Thus the photon brightness perturbation at instantaneous decoupling isHigher angular moments are suppressed by the photon mean-free-path to perturbation-length ratio in the tight-coupling approximation. The idealized equilibrium assumption neglects those finite-mean-free-path corrections, polarization effects and finite thickness of the last-scattering surface. It is not an assertion that collisionless photons remain isotropic after decoupling. The distinction between a fluid treatment and the subsequent hierarchy is discussed in the primary synchronous and Newtonian perturbation treatment.

Let $\mu =\hat{\mathbf{k}}\cdot \mathbf{n}$ and $W(\tau )=e^{ik\mu (\tau -{\tau }_0)}$. Multiplication of the synchronous photon brightness equation by its integrating factor givesIntegration from ${\tau }_{\ast }={\tau }_{\mathrm{dec}}$ to ${\tau }_0$ and division by four yieldFor the real-space observation point ${\mathbf{x}}_0$, the corresponding background ray is $\mathbf{x}(\tau )={\mathbf{x}}_0-\mathbf{n}({\tau }_0-\tau )$. Fourier inversion gives the synchronous Sachs-Wolfe line-of-sight formulaEvery source is evaluated at its appropriate point on the ray, rather than at the observation position. Here $\mathbf{n}$ is photon propagation direction, opposite to the outward observer-to-source sky direction; changing that convention changes the written sign of the Doppler projection consistently.

The first term is the intrinsic photon temperature at last scattering. On wavelengths larger than the sound horizon, it retains primordial adiabatic information; inside that horizon the photon-baryon fluid has acoustic energy density oscillations, producing the acoustic temperature peaks. The second term is the Doppler effect from the fluid's bulk motion at emission. A regular long-wavelength velocity is gradient-suppressed; the Doppler contribution is important near and below the acoustic horizon, with its acoustic phase displaced from that of the energy density oscillation.

The integral is the gravitational energy shift produced by the changing spatial metric tensor along the photon path. In synchronous gauge it includes endpoint gravitational redshifts as well as evolving-potential effects. It must not be identified solely with the Integrated Sachs-Wolfe effect: in Newtonian gauge the same total separates into the ordinary last-scattering Sachs-Wolfe effect, local observer terms, and an integral of time-varying potentials. The ordinary contribution is important on scales outside the sound horizon at decoupling. A late integrated contribution is strongest on large angular scales when dark energy or curvature changes the potentials; evolution around matter-radiation equality also contributes near the first acoustic scales.

The relevant comoving sound scale is $r_s({\tau }_{\ast })={\int }^{{\tau }_{\ast }}{c}_s(\tau )d\tau$, while the projection distance is $D_{\ast }\simeq {\tau }_0-{\tau }_{\ast }$. Modes project roughly to photon temperature multipole $\ell \sim k{D}_{\ast }$, and the acoustic scale is $\ell \sim D_{\ast }/r_s$, with successive peak spacing of order $\pi D_{\ast }/r_s$. At much shorter wavelengths Silk damping and the finite last-scattering thickness suppress energy density and velocity anisotropies; the instantaneous-decoupling approximation alone does not model that cutoff. These scale distinctions explain why the three terms cannot be assigned one common range of important wavelengths.

In the tight-coupling approximation, repeated Thomson scattering drives the photon distribution toward isotropy in the local baryon rest frame. At linear order a temperature perturbation gives a monopole and a small boost gives only a dipole. A photon propagating along $\mathbf{n}$ has rest-frame energy $p(1-\mathbf{n}\cdot \mathbf{v})$, soSince the rest-frame density contrast is ${\delta }_{\gamma }=4{\Theta }_{\mathrm{rest}}$, the initial photon brightness perturbation isThe quadrupole and higher photon temperature multipoles are small when the conformal collision rate $\kappa =a{n}_e{\sigma }_T$ greatly exceeds both $k$ and $\mathcal{H}$. For example the quadrupole is suppressed by a factor of order $k/\kappa$ relative to the dipole. Instantaneous decoupling idealizes this strongly coupled state as the initial data for collisionless propagation; finite mean free path and polarization corrections are omitted from that initial truncation.

Multiply the synchronous photon brightness equation by $e^{ik\mu \tau }$ and integrate from ${\tau }_{\ast }$ to ${\tau }_0$. The line-of-sight solution for free-streaming photons isThe phase is essential: it specifies the emission point and the points traversed by the ray. In real space define $\mathbf{x}(\tau )={\mathbf{x}}_0-\mathbf{n}({\tau }_0-\tau )$ and ${\mathbf{x}}_{\ast }=\mathbf{x}({\tau }_{\ast })$. Dividing by four gives the synchronous Sachs-Wolfe line-of-sight formula,The spatial arguments suppressed in a shorthand version must be understood in this retarded sense; evaluating the emission terms at ${\mathbf{x}}_0$ instead would discard free streaming. Here $\mathbf{n}$ is propagation direction. The observer-to-source direction on the sky is $-\mathbf{n}$, explaining the opposite Doppler sign often used with that convention. The formula describes the synchronous-frame observer; an additional observer peculiar velocity contributes a dipole.

The first term is the intrinsic temperature perturbation at emission, since ${\delta }_{\gamma }=4\Theta$. The second is the Doppler CMB anisotropy from the plasma bulk velocity. The third is gravitational frequency shifting by the time-dependent spatial metric along the ray, and can include scalar and tensor metric contributions. In scalar Newtonian gauge, rearranging that metric contribution supplies the emission gravitational potential as well as the Integrated Sachs-Wolfe effect, up to the observer terms. It should not be identified solely with the integrated Sachs-Wolfe term: a constant Newtonian potential need not give a vanishing synchronous metric integral.

On large angles, modes are outside or comparable to the sound horizon at emission. For the growing adiabatic mode during matter domination, the intrinsic and emission-potential contributions give the ordinary Sachs-Wolfe effect, ${\Theta }_0+\Psi =\Psi /3$ in Newtonian variables. Velocity is a gradient effect, of order $k{\tau }_{\ast }\Psi$ on superhorizon scales, so its Doppler contribution is suppressed there. A scale-invariant primordial curvature spectrum produces approximately the familiar large-angle plateau in $\ell (\ell +1)C_{\ell }$. Evolving potentials during radiation-to-matter transition or late acceleration can add integrated gravitational anisotropy; constant matter-era scalar potentials do not generate a genuine scalar integrated Sachs-Wolfe source.

On smaller angles, the photon-baryon fluid has undergone acoustic oscillations. Its density and velocity are approximately in cosine and sine phases, producing acoustic structure in the intrinsic and Doppler sources and their projection onto angular multipoles. Silk damping suppresses sufficiently short wavelengths before decoupling, while the finite width of last scattering also smooths the real small-angle signal. These corrections refine the instantaneous-decoupling model rather than follow from a collisionless equation alone. The separate intrinsic, velocity and metric pieces depend on gauge convention; the complete observed anisotropy, including the appropriate observer terms and removal of the unobservable mean, is the physical quantity.

The collision rate $\kappa =a{n}_e{\sigma }_T$ is the scattering probability per unit conformal time. Thomson scattering approximately conserves photon number and energy in the electron rest frame, but redistributes propagation directions. Consequently the collision term drives the photon angular distribution toward an isotropic monopole plus the dipole associated with the baryon/electron bulk velocity. It damps the photon-baryon relative velocity and the quadrupole and higher directional anisotropies. Photon and baryon momentum-transfer terms have equal and opposite forces when weighted by their background enthalpies.

The metric terms describe the change in photon energy along a geodesic in a time-dependent spatial metric. The trace perturbation supplies an isotropic redshift, while the scalar trace-free perturbation supplies a direction-dependent redshift proportional to $3{\mu }^2-1$. Brightness is proportional to the fourth power of temperature, accounting for the factor four relative to the temperature shift. These are gravitational sources, not scattering terms.

When $\kappa \gg k,\mathcal{H}$, angular anisotropy is scattered away faster than free streaming or background expansion can regenerate it. In the photon Boltzmann hierarchy, a moment of order $\ell \ge 2$ has damping of order $\kappa {\Delta }_{\ell }$ and is sourced by neighboring moments with coefficients of order $k$. Successive moments are therefore suppressed by powers of $k/\kappa =k{\tau }_c$. This is the tight-coupling approximation. At leading perfect-fluid order the quadrupole can be neglected; its first nonzero correction must be retained to calculate the diffusion damping in part (iii).

The printed equations do not all share a consistent velocity convention. To fix the intended one, take the Fourier velocity to be $v^i=i{k}^i\theta$. Ignoring the metric source, the synchronous perfect-fluid density equation then givesThe PDF's baryon continuity sign is opposite to this convention. Keeping it literally would make the adiabatic ratio and common velocity in part (ii) incompatible for any nonstationary perturbation. Likewise, the momentum-transfer terms for these potentials are $-\kappa ({\theta }_{\gamma }-{\theta }_b)$ and $-R\kappa ({\theta }_b-{\theta }_{\gamma })$, without the printed $k^{-2}$.

For isotropic scattering in this convention the brightness source is $\kappa [{\delta }_{\gamma }+4ik\mu {\theta }_b-\Delta ]$. Its dipole relaxes toward the baryon velocity. The printed source involving ${\theta }_{\gamma }-{\theta }_b$ is not this source, and cannot yield the stated Euler drag by taking its dipole. The following solutions use the consistent continuity and momentum equations. These repairs do not change the requested tight-coupling reduction because the opposite drag forces cancel from its total-momentum equation. The velocity-divergence convention in the standard photon and baryon evolution equations is converted to this potential convention by $\vartheta =-k^2\theta$.

Use the corrected baryon density equation from part (ii). Subtracting $4{\delta }_b^{\prime }/3$ from the photon equation gives $({\delta }_{\gamma }-4{\delta }_b/3{)}^{\prime }=(4/3)k^2({\theta }_{\gamma }-{\theta }_b)$. Tight coupling therefore preserves adiabatic initial conditions to leading order: ${\theta }_b={\theta }_{\gamma }=\theta$ and ${\delta }_b=3{\delta }_{\gamma }/4$.

Multiply the photon Euler equation by the paper's inertia ratio $R$ and add the baryon Euler equation before setting the velocities equal. The Thomson scattering terms cancel exactly. With $\mathcal{H}=a^{\prime }/a$ and $c_b^2$ denoting the baryon sound speed, this yieldsNo derivative of $R$ occurs in this instantaneous weighted addition; writing a derivative of total momentum instead would also differentiate the background inertia. Define the photon-baryon sound speedThen the coupled equations areThe sound speed is radiation pressure divided by the total photon-baryon inertia. Baryons add inertia with little pressure, reducing it below $1/\sqrt{3}$. The combined formula is meaningful despite the printed collision-factor error because the same erroneous factor cancels; maintenance of the assumed adiabatic relation does require correcting the baryon continuity sign.

In the decoupled perfect fluid in general relativity limit, photons have $w=1/3$ and no anisotropic stress, while cold baryons have $w\simeq 0$ with a small pressure correction $c_b^2{\delta }_b$. The density equations express particle and energy conservation; Thomson scattering transfers momentum rather than leading-order rest-frame energy. The photon Euler equation consequently gains a drag toward the baryon velocity potential, and the baryon Euler equation gains the equal-and-opposite drag weighted by their inertias.

Let ${\Gamma }_T=a{n}_e{\sigma }_T>0$. With the velocity potential used in part (i), these forces are $-{\Gamma }_T({\theta }_{\gamma }-{\theta }_b)$ and $-R{\Gamma }_T({\theta }_b-{\theta }_{\gamma })$, with $R=4{\bar{\rho }}_{\gamma }/(3{\bar{\rho }}_b)$. Multiplication by the photon and baryon enthalpies, $4{\bar{\rho }}_{\gamma }/3$ and ${\bar{\rho }}_b$, makes their sum zero. Note that this paper's $R$ is the reciprocal of the usual baryon loading parameter. Thomson scattering also relaxes the photon quadrupole, giving the $-{\Gamma }_T{\sigma }_{\gamma }$ term in the stipulated simplified hierarchy. This quadrupole is an anisotropic stress, omitted by a perfect-fluid model; finite relaxation creates Silk damping.

The printed baryon density sign and collision normalization are inconsistent with part (i). The density equation must be ${\delta }_b^{\prime }-k^2{\theta }_b=-h^{\prime }/2$. As a direct check, with $h^{\prime }=0$, nonzero common velocity, and initially ${\delta }_{\gamma }=4{\delta }_b/3$, the printed plus sign gives $({\delta }_{\gamma }-4{\delta }_b/3{)}^{\prime }=8k^2\theta /3$, violating the tight-coupling adiabatic condition immediately. Likewise division of the collision forces by $k^2$ belongs to neither this potential convention nor the usual common velocity-divergence convention; the physical relaxation rate is ${\Gamma }_T$, independent of $k$. Both printed drag terms still cancel in the weighted sum used below, so that particular algebraic combination is unaffected by this normalization error.

For a nonrelativistic distribution, $T^{ij}=\int d^{3}p{p}^i{p}^{j}f/E$ is down by $v^2$ relative to the rest-mass energy density, and the corresponding anisotropic stress is also suppressed. The streaming terms coupling successive angular moments are proportional to $kv$, and stresses from the velocity distribution disappear in the cold limit. Thus, on scales where velocity dispersion is negligible, density and bulk velocity suffice; the photon Boltzmann hierarchy has no analogous small-speed suppression. This argument needs a cold distribution or a controlled small-dispersion approximation. Collisionless nonrelativistic matter with appreciable dispersion or multistreaming need not be an exact perfect fluid, and near photon decoupling the higher photon moments cannot all remain negligible.

For $R\gg 1$, $3{\tilde{c}}_s^2\simeq 1$, and the baryon-inertia corrections to the common Euler equation vanish. Differentiating the photon continuity equation givesSet aside the gravitational driving terms and balance the quadrupole streaming source against Thomson scattering. With ${\tau }_c={\Gamma }_T^{-1}$ and ${\sigma }_{\gamma }^{\prime }\simeq 0$,The positive friction is Silk damping; this particular truncation gives shear-only photon diffusion damping. Its coefficient follows the supplied simplified quadrupole equation; adding polarization and the full velocity-slip correction changes that coefficient, as in the general photon-baryon diffusion damping equation.

For completeness, put ${\omega }_0=k/\sqrt{3}$, $b=4k^2{\tau }_c/15$ and $D=\frac{1}{2}{\int }_{{\tau }_0}^{\tau }b(s)ds$. Substitution ${\delta }_{\gamma }=e^{-D}y$ gives the exact transformed equationFor constant ${\tau }_c$ the two independent solutions are $e^{-b(\tau -{\tau }_0)/2}\cos [\sqrt{{\omega }_0^2-b^2/4}(\tau -{\tau }_0)]$ and the corresponding sine. Their frequency correction begins at second order in ${\tau }_c$. For variable ${\tau }_c$, retaining first order alone does not automatically remove $b^{\prime }/2$: a general first-order representation isThe intended cosmological approximation additionally uses slow variation on an acoustic period, $\mathcal{H}/k\ll 1$, with tight coupling $k{\tau }_c\ll 1$. Then the $b^{\prime }$ correction to the frequency is subleading and the two independent photon-baryon acoustic oscillator solutions take the resummed damping formA shift of the time origin is absorbed in $A,B$. If the ionization fraction is constant, $n_e\propto a^{-3}$ gives ${\tau }_c\propto a^2\propto {\tau }^4$ during matter domination, and an initial origin at zero gives $k_D^{-2}=2\tau {\tau }_c/75$. Across recombination the ionization fraction varies, so the integral is the general answer. The PDF prints $O({\tau }_c^{-2})$, while the TeX aid prints $O({\tau }_c^2)$; the small-mean-free-time expansion requires the latter interpretation.

The intrinsic photon temperature fluctuation is $\delta T/T={\delta }_{\gamma }/4$. The cosine and sine produce the Cosmic microwave background acoustic peaks; the envelope suppresses their small-scale amplitudes. The corresponding contribution to the Cosmic microwave background power spectrum contains $e^{-2k^2/k_D^2}$, before angular projection and the other temperature sources are included. Thus Silk damping erases the high-multipole acoustic structure progressively rather than shifting every peak to a new frequency.