Anisotropic stress 2026-10-06
In a local orthonormal tetrad, anisotropic stress is the trace-free spatial part of the stress-energy tensor, conventionally denoted . Some linear cosmological perturbation theory conventions instead use and write . The sign must be specified before using a photon Euler equation. The photon quadrupole determines the linear photon anisotropic stress.
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 to
which 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 .
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 ; the velocity sourcing the quadrupole follows . 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, , 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.
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 . A first-order boost by the fluid velocity adds a temperature photon dipole ; it does not generate a photon quadrupole until second order in that boost. Thus the photon brightness perturbation at instantaneous decoupling is
Higher 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 and . Multiplication of the synchronous photon brightness equation by its integrating factor gives
Integration from to and division by four yield
For the real-space observation point , the corresponding background ray is . Fourier inversion gives the synchronous Sachs-Wolfe line-of-sight formula
Every source is evaluated at its appropriate point on the ray, rather than at the observation position. Here 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 , while the projection distance is . Modes project roughly to photon temperature multipole , and the acoustic scale is , with successive peak spacing of order . 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, the Thomson-scattering rate is much faster than both spatial streaming and cosmological evolution: and . Frequent scattering locks the photon and baryon velocities and suppresses higher angular multipoles. The photon dipole is leading order, while the photon quadrupole is first order in the mean-free-time expansion.
Use the quadrupole equation. Its time derivative is smaller than its collision term by an additional expansion parameter, and starts one order beyond . The leading balance is therefore
Thus for the temperature-only collision operator supplied here,
The temperature-only tight-coupling quadrupole is not numerically interchangeable with a hierarchy that includes polarization feedback. The latter changes the leading collision coefficient and gives the existing photon quadrupole in tight coupling with polarization instead. Here the stated temperature-only approximation is the one being expanded.
For , scalar axisymmetry leaves only in the photon quadrupole: with the normalization used above. The two spin-weighted spherical harmonics satisfy . Translation from last scattering supplies . Thus the two helicity combinations have the same scalar source:
where is a common normalization constant independent of . In this meridian polarization basis the scalar contribution has ; this does not mean that every rotated polarization basis has zero .
Apply the spin-raising and lowering operators to an axisymmetric scalar potential . The first action is . In the second action, the intermediate spin is respectively or ; both give
It follows that
Subtracting the two equations gives . The possible affine consists solely of the unobservable kernel of the twice-applied spin operators. Defining polarization potentials to have no such components sets the physical scalar B mode to zero, . This is the low-multipole kernel of polarization potentials, not a physical dipole polarization.
For the remaining equation integrates to . Remove the same monopole/dipole kernel and use the plane-wave expansion. Thus the scalar E-mode radial projection is
The overall minus sign is absorbed in the stated proportionality constant. Although the unprojected antiderivative has , its divergent monopole and dipole pieces are discarded; the physical limit is finite because .
The temperature-only tight-coupling quadrupole is suppressed by . For regular superhorizon scalar modes, the photon dipole is itself proportional to , so the recombination quadrupole and the resulting polarization are especially small at large angular scales. Last-scattering polarization is generated by the small departure from perfect tight coupling, not by an isotropic monopole. Reionization can generate a separate large-angle signal after recombination and is absent from the instantaneous-last-scattering approximation here. The no-B-mode conclusion applies to linear scalar sources; it does not include tensor sources or conversion by lensing.
The tight-coupling approximation requires the positive scattering rate to exceed both and the conformal Hubble parameter. The photon-baryon velocity slip and photon quadrupole are then small, and the leading fluid motion has . One must combine the equations before setting the slip to zero: a small slip times a large collision rate can exert a finite force.
Add the photon Euler equation to times the baryon Euler equation with Thomson drag. The collisions cancel and, neglecting at this order, the result is
Since , the common velocity obeys
The inertia of the baryons reduces the photon-baryon sound speed to .
Differentiate the photon continuity equation and substitute this velocity equation:
Thus the photon-baryon acoustic oscillator is
The terms on the right drive the oscillation gravitationally; the first-derivative term comes from changing baryon inertia. Diffusion damping enters only beyond this leading tight-coupling approximation.
Let and let differentiate along the unperturbed photon trajectory. Use to decompose the Thomson scattering gain term:
The Orthogonality of Legendre polynomials identifies as a pure photon quadrupole. Dropping is an approximation, justified when tight coupling makes the scattering quadrupole small; the supplied unpolarized equation already neglects polarization corrections. The Cosmic microwave background line-of-sight solution requested here omits that quadrupole source.
Set . Adding to the Photon Boltzmann equation with Thomson scattering yields
Along , the integrating factor is because . Consequently
An early optically thick boundary makes vanish for finite initial perturbations, and . Therefore
Here and are partial time derivatives, not total derivatives along the ray. We work at first order in scalar cosmological perturbations, evaluate sources along the background radial null geodesic in FLRW spacetime, and ignore observer peculiar velocity. Perturbed paths and lensing affect this first-order source solution only at higher perturbative order. A moving observer adds its local Doppler dipole.
The original PDF, unlike the transcribed TeX, also asks for the instantaneous-last-scattering limit and its interpretation. For with no later scattering, and a matter-era growing mode with constant and , the Integrated Sachs-Wolfe effect vanishes. Hence
The four terms are the intrinsic emission temperature, the gravitational potential at emission, the observer's potential, and the Doppler CMB anisotropy. The two potentials describe the gravitational redshift; the observer term is independent of direction and is absorbed into the measured temperature monopole. The intrinsic temperature and emission potential form the Sachs-Wolfe combination. On superhorizon scales for adiabatic initial conditions in matter domination, , giving the ordinary Sachs-Wolfe effect . The Doppler sign is positive with the paper's propagation direction ; the viewing direction is .
Write and use the Legendre polynomial recurrence relation. The definition in the question is inverted by
For Newtonian-gauge potentials the scanned geodesic equation is the standard identity
Its monopole and dipole are and . Projecting the photon Boltzmann hierarchy onto and therefore gives
The first is the photon continuity equation; the second is the photon Euler equation, with the photon quadrupole providing the anisotropic-stress term. Indeed, with photon density contrast and velocity divergence ,
Project the photon Boltzmann equation onto Legendre polynomials. The angular average uses and , so the collision term has no monopole and
For the dipole, the recurrence relation makes free streaming couple to and . The gravitational term contributes , while projection of gives the baryon-velocity source. Thus
Neglecting the photon quadrupole gives
When gravity and expansion can be neglected and is constant, the photon monopole satisfies
The first damping contribution comes from photon-baryon velocity slip; the second comes from the photon quadrupole. This is a local form of Silk damping.
For the expansion without a factor, the leading polarization relation is . The photon Boltzmann hierarchy then gives . This small photon quadrupole supplies the shear-viscosity part of Silk damping.
For a scalar photon quadrupole emitted at distance , the meridian-basis source is . Integrating its scalar-potential equation and removing the low-multipole kernel of polarization potentials yields , with no scalar . The apparent divergence lies in discarded low multipoles; makes the physical small- limit regular.