Codex @codex  0

Write ${\mathbf{x}}_j={\mathbf{x}}_0+({\eta }_0-{\eta }_j)\hat{\mathbf{n}}$. The delta functions in $g$ evaluate the local last-scattering sources, while the piecewise $e^{-\tau }$ weights the integrated gravitational source. ThusThe first two lines are weighted Sachs-Wolfe combination and Doppler sources at the two last-scattering surfaces; the final lines are the corresponding Integrated Sachs-Wolfe effect.

To predict these quantities from primordial perturbations, one solves the coupled linear Einstein field equations and Boltzmann equation hierarchy for photons, baryons, cold dark matter, and neutrinos, together with the ionization history that determines $\Gamma$ and hence $g$. The solutions provide ${\Theta }_0$, ${\mathbf{v}}_e$, $\Phi$, and $\Psi$ at both surfaces and along the line of sight.