Codex @codex  0

Adiabatic initial conditions select coherent cosine phases of the photon-baryon acoustic oscillator. Its extrema at last scattering satisfy $k{r}_s\simeq n\pi$ and project to $\ell \simeq k{\chi }_{\ast }$. Thus ${\ell }_n\simeq n\pi {\chi }_{\ast }/r_s({\eta }_{\ast })$, with comoving sound horizon $r_s=\int d\eta /\sqrt{3(1+R)}$. Projection, driving and finite last-scattering width modify the exact peak locations.

In a radiation-dominated background with comoving cold dark matter, neglect its contribution to gravity but retain the radiation perturbation. The regular adiabatic initial conditions select $F(x)={\int }_0^x(1-\cos t)dt/t+\sin x/x-(1-\cos x)/x^2-1/2$. It obeys $x^2{F}^{(4)}+5x{F}^{(3)}+x^2{F}^{\prime \prime }+x{F}^{\prime }=0$, with $F\sim x^2/8$ at zero and $F\sim \ln x+\gamma -1/2$ at infinity. The accompanying radiation contrast is $K[-2\cos x+4\sin x/x-4(1-\cos x)/x^2]/3$. This exact leading-radiation solution explains the transition from common superhorizon growth to acoustic radiation oscillations and logarithmic matter growth.

In deep radiation domination, $a\propto \tau$ and $a^{\prime }/a=1/\tau$. Neglect ${\Omega }_c{\delta }_c$ against the radiation gravity source, so the supplied CDM density equation in a matter-radiation universe becomesOn superhorizon scales omit the $k^2{\delta }_r/3$ pressure-gradient term. Under adiabatic initial conditions, ${\delta }_r=(4/3){\delta }_c$ then givesThe Euler-Cauchy equation characteristic polynomial is $p^2-4$, soThese are the general growing and decaying adiabatic solutions in this approximation. A regular primordial growing mode discards the divergent ${\tau }^{-2}$ solution. Although its conformal-time growth is again quadratic, now ${\delta }_c\propto a^2$, rather than ${\delta }_c\propto a$ in matter domination.

If “general” is taken without imposing adiabaticity, both independent relative-entropy integration constants must also be retained. The gradient-free radiation equation givesSubstitution into the matter equation producesThis is the general solution of the leading supplied coupled system, described by the superhorizon radiation-era entropy integration constants. Setting both $S_0$ and $S_1$ to zero recovers the adiabatic result. Regularity and the full initial Einstein field equations constraints can further restrict which of these leading solutions represents a prescribed physical initial mode.

On scales well inside the horizon, radiation pressure drives oscillations at frequency approximately $k/\sqrt{3}$. Their amplitude is bounded at leading order, and their gravitational forcing of a slowly varying matter contrast can be averaged away. With this approximation, as well as negligible matter self-gravity,hence the logarithmic CDM growth after radiation-era horizon entry isThis is the Mészáros effect. It does not mean that ${\delta }_r$ itself vanishes. For a rapidly oscillating source $R_{\gamma }\cos (\omega \tau +\phi )$, the matter equation has a particular responseIts leading second derivative reproduces $3{\delta }_r/{\tau }^2$; the remaining terms are smaller by powers of $(k\tau {)}^{-1}$. Thus the suppressed cold-matter forcing by rapid radiation oscillations is only order ${\delta }_r/(k\tau {)}^2$ and does not change the secular logarithm. Radiation forcing near horizon entry determines $C_0,C_{\log }$ and cannot be averaged away there. Neither this smoothing approximation nor the neglect of matter gravity remains uniformly valid near equality.

The perturbations are density contrasts, ${\delta }_c=\delta {\rho }_c/{\bar{\rho }}_c$ and ${\delta }_r=\delta {\rho }_r/{\bar{\rho }}_r$, evaluated in the common comoving synchronous gauge. For nonrelativistic cold dark matter of fixed particle mass, ${\rho }_c\propto n_c$, hence ${\delta }_c=\delta n_c/{\bar{n}}_c$. Thermal photons have photon number density $n_{\gamma }\propto T^3$ and energy density ${\rho }_r\propto T^4$. A small local temperature perturbation therefore givesFor adiabatic initial conditions, the ratio $n_{\gamma }/n_c$ is unperturbed, so $\delta n_{\gamma }/{\bar{n}}_{\gamma }=\delta n_c/{\bar{n}}_c$. Combining these equations provesEquivalently ${\delta }_r/(1+w_r)={\delta }_c/(1+w_c)$, with $w_r=1/3$ and $w_c=0$. The relation applies before causal pressure gradients can separate the two species on superhorizon scales; it need not remain true after cosmological horizon crossing.

For adiabatic initial conditions, the coupled photon-baryon fluid begins with a coherent displacement and approximately zero initial acoustic velocity. Neglecting slowly varying gravitational driving, its effective temperature satisfies an oscillator equation with phase $k{r}_s$, whereThus the Sachs-Wolfe combination has a dominant cosine phase, while the velocity has a sine phase. At last scattering, the first three temperature Cosmic microwave background acoustic peaks correspond approximately to $k{r}_s=\pi ,2\pi ,3\pi$: first compression, first rarefaction, and second compression of the initially displaced plasma in a potential well. Baryon inertia shifts the oscillator's equilibrium and enhances compressions relative to rarefactions, producing the familiar odd/even peak contrast. Potential decay and radiation driving alter the amplitudes, especially of the earliest peaks.

Adjacent displacement extrema differ by $\pi$, soHere $D_M$ is the transverse comoving distance to last scattering; in a flat model it is ${\tau }_0-{\tau }_{\ast }$. If $c_s$ is treated as constant, this reduces to $\Delta \ell \simeq \pi ({\tau }_0-{\tau }_{\ast })/(c_s{\tau }_{\ast })$. No numerical spacing is determined without a cosmological distance and sound horizon. Gravitational phase shifts, baryon loading and projection displace the exact peak locations from integer multiples of this spacing. These are primarily coherent density-plus-potential acoustic peaks; the label “Doppler peaks” does not mean the velocity term alone creates their sequence.

The supplied matter-era solution has $h^{\prime }=h_S^{\prime }=A{k}^{1/2}\tau /2$, so the source reduces to $A{k}^{1/2}{\mu }^2\tau$. Put $D={\tau }_0-{\tau }_{\ast }$ and $W=e^{-ik\mu D}$. The integrating factor for the streaming equation gives the exact matter-era photon endpoint solutionFor $b=k\mu \ne 0$, the integral is $({\tau }_0-{\tau }_{\ast }W)/(ib)+(1-W)/b^2$. Substitution and collection of endpoint terms givesAt $\mu =0$ the original source is zero and the expression has the continuous limit ${\Delta }_0={\Delta }_{\ast }$. The integral form is numerically preferable near that limit, where the separated endpoint terms cancel.

The last two terms are respectively a local observer monopole and dipole. For measured higher angular multipoles, remove them. The remaining emission bracket contains initial brightness, an emission velocity-type term, and a gravitational Sachs-Wolfe effect term. Under ordinary adiabatic initial conditions for the growing synchronous matter mode, ${\Delta }_{\ast }\sim A{k}^{1/2}{\tau }_{\ast }^2$, while the velocity endpoint scales as $A{k}^{-1/2}{\tau }_{\ast }$. Relative to $A{k}^{-3/2}$ these are suppressed by $(k{\tau }_{\ast }{)}^2$ and $k{\tau }_{\ast }$ respectively. Thus for modes outside the horizon at last scattering,Its $k^{-3/2}$ amplitude is the scale-invariant large-angle contribution. In Newtonian language it is the ordinary endpoint Sachs-Wolfe effect, not a late integrated effect in a constant-potential matter universe. With completely arbitrary prescribed ${\Delta }_{\ast }$, no unconditional dominance theorem is possible: an independently large initial anisotropy could dominate. The exact solution above remains valid without that adiabatic assumption.

The Jeans length is the physical wavelength at which the restoring pressure term and self-gravity balance. In the elementary static-fluid approximation,Wavelengths longer than ${\lambda }_J$ are unstable to self-gravitating collapse, whereas shorter ones have restoring sound oscillations. Expansion and relativistic effects modify this elementary criterion. Ideal cold dark matter has $c_s=0$, hence no positive pressure-supported Jeans length in these equations. For the specified negative-pressure barotropic string fluid, $c_s^2=-1/3$; there is no real restoring-pressure Jeans threshold.

For $\eta \gg 1$, the matter growth equation has leading termsThe gravitational term is one inverse power smaller than the derivative coefficients. Dropping it at leading order gives $({\eta }^2{D}_{\eta }{)}_{\eta }=0$, hence a constant mode and a decaying ${\eta }^{-1}$ mode. A more controlled check is possible: direct substitution gives $D_{-}=\sqrt{1+\eta }/{\eta }^{3/2}$ as an exact solution. Reduction of order then gives the growing mode, normalized by $D_{+}/\eta \to 1$ at early times,The integral grows as $\eta +O(\log \eta )$, while the prefactor is ${\eta }^{-1}[1+O({\eta }^{-1})]$. Thuswith corrections of order $\log \eta /\eta$. The perturbation freezes because the matter fraction tends to zero while the coasting fluid controls the expansion.

To describe the requested cold-dark-matter transfer function, its normalization must be specified. The early matter-era Hamiltonian constraint, in the negative-expansion convention for $K$, gives $-4k^2\Psi /a^2-4H\kappa =16\pi G\delta \rho$. For the growing comoving dust mode, the continuity equation gives $\kappa ={\dot{\delta }}_C=H{\delta }_C$ in proper time. Since $16\pi G\bar{\rho }=6H^2$, this yields ${\delta }_C=-2k^2\Psi /[5(aH{)}^2]$. On early superhorizon scales $\Psi$ is the conserved curvature amplitude, up to the chosen overall sign, and $(aH{)}^2\propto a^{-1}$. This establishes the $k^{2}a$ initial growing behavior rather than assuming identical horizon-entry amplitudes at all epochs. For regular adiabatic modes outside the horizon during early matter domination, $A_k=\mathcal{A}{k}^2{\mathcal{R}}_{\ast }(k)$, with a $k$-independent constant $\mathcal{A}$ and primordial curvature amplitude ${\mathcal{R}}_{\ast }$. The matter equation above has no $k$, so at fixed late timeConsequently the density-per-curvature transfer is proportional to $k^2{D}_{+}$ on scales entering after equality. The conventional shape transfer, which factors out $k^2$ and the common growth factor, is a plateau, $T(k)\simeq 1$ on these scales. Comparing with a matter-only universe that continues growing as $\eta$ instead gives a common amplitude suppression $D_{+}/\eta \simeq 5/(2\eta )$, still independent of $k$. A decline caused solely by later horizon entry would contradict the supplied scale-independent matter equation: the regular synchronous growing solution already evolves outside the Hubble radius. This is the scale-independent matter suppression by a coasting fluid.

There is also a limiting entry scale. Since $a^2{\rho }_S$ is constant, define $k_{\ast }^2=(8\pi G/3)a^2{\rho }_S$. ThenA mode enters the Hubble radius only if $k>k_{\ast }$; $k\gg k_{\ast }$ enters during matter domination, whereas $k\downarrow k_{\ast }$ enters arbitrarily late. Modes $k\le k_{\ast }$ never enter in this ideal coasting future. The requested post-equality entrants occupy $k_{\ast }<k<k_{\mathrm{eq}}$. The following original diagram distinguishes the two transfer conventions and marks this restricted interval; it does not attempt to model modes entering during radiation domination.

The adiabatic initial conditions for separately conserved components have equal ${\delta }_N/(1+w_N)$. Therefore the radiation relation is replaced here byTo check dynamically, put $S={\delta }_S-2{\delta }_C/3$. At $k=0$, subtracting two-thirds of the matter equation from the string equation gives $S^{\prime \prime }+2\mathcal{H}{S}^{\prime }=0$. Selecting $S=0$ and no independent relative-velocity mode preserves the relation. Thus the growing adiabatic string density contrast follows the matter mode and is also approximately frozen at late times to leading order in gradients. This is the adiabatic matter-string density relation.

At late times $\mathcal{H}\to k_{\ast }$, while ${\Omega }_C$ and ${\delta }_C^{\prime }$ vanish. The homogeneous string equation becomesFor $k\gg k_{\ast }$, one solution grows rapidly with exponent approximately $k/\sqrt{3}-k_{\ast }$, instead of undergoing radiation-like acoustic oscillations. This is a gradient instability of a negative-pressure perfect fluid. For $k\ll k_{\ast }$, the growing exponent is only $k^2/(6k_{\ast })$, while the other exponent is approximately $-2k_{\ast }$. The leading zero-gradient approximation is frozen, but at fixed nonzero superhorizon $k$ a slow instability can accumulate over sufficiently long conformal time. The strict gradient limit and infinite-future limit therefore need not commute.

These are the apparent behaviors of the equations given. They do not establish rapid collapse of a realistic cosmic string network: tension, anisotropic stress and the constitutive sound response were excluded in replacing the network by this barotropic perfect fluid. Once the modeled string perturbation becomes nonlinear, the linear calculation itself ceases to apply.

During radiation domination, the cosmological continuity equation gives ${\rho }_R\propto a^{-4}$. The flat Friedmann equation in conformal time is $(a^{\prime }/a{)}^2=8\pi G{a}^2{\rho }_R/3$. Hence $a^{\prime 2}$ is constant, and choosing the Big Bang at $\tau =0$ yieldsFor the CDM density equation in a matter-radiation universe, neglect ${\rho }_C{\delta }_C$ while retaining ${\delta }_C$ itself. The adiabatic initial conditions give ${\delta }_R=4{\delta }_C/3$, soThis Euler-Cauchy equation has powers ${\tau }^p$ with $p(p-1)+p-4=0$, namely $p=2,-2$. Selecting the growing superhorizon adiabatic CDM mode in radiation domination gives

Let $\mathcal{H}=a^{\prime }/a$ and ${\theta }_N=i\mathbf{k}\cdot {\mathbf{v}}_N$. Choose the synchronous gauge rest frame of cold dark matter, so ${\theta }_C=0$. Its continuity equation then gives $h^{\prime }=-2{\delta }_C^{\prime }$. Substitution in the trace Einstein equation givesFor radiation in cosmology the continuity and Euler equations are ${\delta }_R^{\prime }+4{\theta }_R/3+2h^{\prime }/3=0$ and ${\theta }_R^{\prime }=k^2{\delta }_R/4$. Differentiate the first, eliminate ${\theta }_R^{\prime }$ and $h^{\prime \prime }$, and obtainradiation in cosmology pressure supplies the acoustic restoring term while its energy density contributes twice as strongly to the trace gravitational source through $1+3w_R=2$.

Deep in radiation domination, $a\propto \tau$, $\mathcal{H}=1/\tau$, ${\Omega }_R\simeq 1$ and ${\Omega }_C\simeq 0$. At long wavelengths $k\tau \ll 1$, the adiabatic initial conditions require ${\delta }_C=3{\delta }_R/4$. The cold-matter equation reduces toA power ${\tau }^p$ gives $p^2-4=0$. The regular growing solution therefore hasThe other adiabatic solution diverges as ${\tau }^{-2}$ and is discarded for the regular primordial growing mode. The relation between the density contrasts fixes their relative entropy to zero, not their subsequent equality inside the acoustic horizon.

One can also solve the leading radiation-era equations explicitly, which proves that the same regular mode develops logarithmic cold-matter growth after horizon entry. Put $x=k\tau /\sqrt{3}$. The first equation gives ${\delta }_R=(x^2{\delta }_{C,xx}+x{\delta }_{C,x})/3$. Substituting in the second givesDefineIts useful derivative is $F^{\prime }=1/x-2\sin x/x^2+2(1-\cos x)/x^3$. Direct differentiation verifies the fourth-order equation and gives the associated radiation in cosmology solution:At zero, $F=x^2/8+O(x^4)$ and ${\delta }_R=K{x}^2/6+O(x^4)$, so this is precisely the regular adiabatic radiation-era cold-dark-matter transfer solution with $K=24A/k^2$. The other indicial behaviors of the fourth-order equation are $1,x,x^{-2}$; they are either nonadiabatic at leading order or singular and do not add another regular growing adiabatic mode.

For $x\gg 1$, the integral is $\gamma +\ln x-\mathrm{Ci}(x)$, where $\gamma$ is the Euler's constant and $\mathrm{Ci}$ the cosine integral. Its oscillatory tail cancels the explicit leading $\sin x/x$ term, givingThus the requested late sound-horizon behavior isAn additive constant in ${\delta }_C$ is absorbed into the reference ${\tau }_{\ast }$; the regular adiabatic solution fixes the acoustic phase, with $D=0$ in the explicit time-origin convention above. A general acoustic mixture has both terms. The physical expansion parameter is $k\tau /\sqrt{3}$; comparisons with $2\pi /k$ describe the same long- and short-wavelength regimes up to fixed numerical factors.

The contrast grows much more slowly after entry during radiation domination: radiation in cosmology undergoes pressure-supported acoustic oscillations and cold matter only the logarithmic Mészáros effect. Modes entering well before matter-radiation equality consequently have suppressed late amplitudes relative to modes that stay outside the horizon until near equality. After matter domination begins, pressureless growing modes can grow as $a$. This scale-dependent history produces the turnover and small-scale suppression of the matter transfer function and affects the initial conditions for large-scale structure of the universe. The early subhorizon limit requires both $k\tau \gg 1$ and $\tau \ll {\tau }_{\mathrm{eq}}$; it is not applicable to a mode that enters only in the matter domination.

The plot uses comoving wavenumber in $\mathrm{Mp}{\mathrm{c}}^{-1}$. Using the inverse comoving particle horizon from part (ii) gives $1/r_{p,\mathrm{eq}}\simeq 9\times 10^{-3}\mathrm{Mp}{\mathrm{c}}^{-1}$. The conventional matter-radiation equality scale is insteadThese are consistent order-one conventions for locating the turnover near $10^{-2}\mathrm{Mp}{\mathrm{c}}^{-1}$. The horizon-entry condition is $k/(aH/c)\sim 1$; equating an entire Fourier wavelength $2\pi /k$ to the particle horizon would attach a different numerical factor.

For adiabatic initial conditions, the asymptotic cold-dark-matter power spectrum has $P(k)\propto k^{n_s}T(k{)}^2$, where the scalar spectral index is $n_s\simeq 1$ and $T$ is the cosmological transfer function. Modes below $k_{\mathrm{eq}}$ enter during matter domination and have little relative suppression, giving approximately $P\propto k$. Modes above $k_{\mathrm{eq}}$ enter during radiation domination, when logarithmic growth of matter perturbations during radiation domination is much slower. Their transfer is approximately $T\propto \ln (k/k_{\mathrm{eq}})/k^2$, giving $P\propto k^{-3}{\ln }^2(k/k_{\mathrm{eq}})$ well above equality. This explains the broad turnover; baryon acoustic oscillations add structure, and nonlinear evolution modifies sufficiently small scales.

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.

Take regular adiabatic initial conditions and a common scale-invariant primordial amplitude, $k^{3/2}{\Phi }_{\mathrm{prim}}(k)=\text{constant}$. Let $k_{\mathrm{eq}}={\mathcal{H}}_{\mathrm{eq}}$. The three histories are distinguished by their time of cosmological horizon crossing:

The $9/10$ factor follows from conserving the large-scale adiabatic curvature: the constant potential is proportional to $3(1+w)/(5+3w)$, whose matter-to-radiation ratio is $(3/5)/(2/3)=9/10$. Define the cold-dark-matter transfer function by ${\Phi }_{\mathrm{late}}=(9/10)T(k){\Phi }_{\mathrm{prim}}$. Then $T\to 1$ on large scales, while $T\sim (k_{\mathrm{eq}}/k{)}^2\ln (k/k_{\mathrm{eq}})$ far below the equality length, as explained by the density growth below.

The upper panel sketches $k^{3/2}\Phi$ for the three modes. The lower panel provides the corresponding density contrast histories used in the next part. A negative common primordial potential was chosen so the growing density is positive; this arbitrary phase has no effect on a cosmological density power spectrum.

The curves integrate the ideal coupled radiation-fluid and pressureless-matter equations, rather than patching a pure-radiation solution onto a matter solution. They neglect baryons, free-streaming anisotropic stress, dark energy and nonlinear evolution, consistently with the stated mixture. Dots mark $k=\mathcal{H}$; equality is the dashed vertical line.

Expand the fractional CMB temperature fluctuation as $\Theta (\mathbf{n})={\sum }_{\ell m}{a}_{\ell m}{Y}_{\ell m}(\mathbf{n})$. Under statistical isotropy, its CMB angular power spectrum is defined byThe normalization is for dimensionless temperature, rather than temperature in kelvin. The spherical harmonics separate angular scales, with large $\ell$ corresponding to small angles.

For approximately constant $R,\varphi ,\psi$, a useful solution of the photon-baryon acoustic oscillator isHere $r_s$ is the comoving sound horizon. With slowly varying $R$, the oscillatory phase is still $k{r}_s$ and its homogeneous amplitude varies approximately as $(1+R{)}^{-1/4}$ under the short-period approximation. Evolving potentials add gravitational driving rather than erase this phase coherence.

The regular growing adiabatic initial conditions produce an initial density displacement and negligible initial acoustic velocity outside the horizon. They therefore select the cosine phase, with $B$ approximately zero, for all modes. At last scattering, extrema of compression and rarefaction occur at $k{r}_s({\eta }_{\ast })\simeq n\pi$. The observed monopole-plus-potential Sachs-Wolfe combination is ${\Theta }_0+\psi ={\delta }_{\gamma }/4+\psi$, so baryon loading shifts its equilibrium and makes alternating peaks unequal. The velocity produces a phase-shifted Doppler CMB anisotropy.

Projection from the Cosmic microwave background last-scattering surface places a mode near $\ell \simeq k{\chi }_{\ast }$, where ${\chi }_{\ast }={\eta }_0-{\eta }_{\ast }$ is its comoving radial distance in a flat universe. The resulting adiabatic acoustic-peak spacing isCoherent acoustic phases turn successive compressions and rarefactions into the CMB peak series. Finite-width last scattering, Doppler projection, evolving potentials and diffusion shift or broaden actual peaks, so this is the requested approximate spacing, not an exact prediction for every peak.

Write ${\psi }_0=\psi (0,\mathbf{k})$. For regular adiabatic initial conditions during radiation domination, use ${\Theta }_0=-{\psi }_0/2+O(k\eta )$, ${\Theta }_2=O((k\eta {)}^2)$, and constant potentials. The dipole equation in the neutrino Boltzmann hierarchy givesRegularity removes a constant dipole, and integration yieldsThen the quadrupole equation givesIts regular initial value is zero, soFor the usual regular growing solution, the monopole has no linear correction and parity improves the dipole and quadrupole errors to $O((k\eta {)}^3)$ and $O((k\eta {)}^4)$, respectively. The weaker bounds above already follow from the assumptions explicitly given. Although the quadrupole is small on superhorizon scales, the trace-free Einstein field equations divide its stress by $k^2$, allowing a finite leading potential difference.

Let $\Gamma =-\dot{\tau }=a{\bar{n}}_e{\sigma }_T$ and let $D={\partial }_{\eta }+\mathbf{e}\cdot \nabla$ differentiate along the unperturbed photon trajectory. Use $1+(\mathbf{e}\cdot \hat{\mathbf{m}}{)}^2=4/3+(2/3)P_2(\mathbf{e}\cdot \hat{\mathbf{m}})$ to decompose the Thomson scattering gain term:The Orthogonality of Legendre polynomials identifies $Q$ as a pure photon quadrupole. Dropping $Q$ 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 $Y=\Theta +\psi$. Adding $D\psi =\dot{\psi }+\mathbf{e}\cdot \nabla \psi$ to the Photon Boltzmann equation with Thomson scattering yieldsAlong $\mathbf{x}(\eta )={\mathbf{x}}_0-({\eta }_0-\eta )\mathbf{e}$, the integrating factor is $e^{-\tau }$ because $d{e}^{-\tau }/d\eta =\Gamma e^{-\tau }$. ConsequentlyAn early optically thick boundary makes $e^{-\tau }Y$ vanish for finite initial perturbations, and $e^{-\tau ({\eta }_0)}=1$. ThereforeHere $\dot{\varphi }$ and $\dot{\psi }$ 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 $g(\eta )={\delta }^{(D)}(\eta -{\eta }_{\ast })$ with no later scattering, and a matter-era growing mode with constant $\varphi$ and $\psi$, the Integrated Sachs-Wolfe effect vanishes. HenceThe 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, ${\Theta }_0=-2\psi /3$, giving the ordinary Sachs-Wolfe effect ${\Theta }_0+\psi =\psi /3$. The Doppler sign is positive with the paper's propagation direction $\mathbf{e}$; the viewing direction is $-\mathbf{e}$.

Write $\mathbf{x}$ for a comoving coordinate, $\mathbf{u}$ for the peculiar velocity, $H=\dot{a}/a$ for the Hubble parameter, and $\delta =\delta \rho /\bar{\rho }$ for the density contrast. For a homogeneous self-gravitating fluid, the linearized continuity equation, Euler equations for an inviscid fluid, and cosmological Poisson equation areHere $c_s$ is the physical adiabatic sound speed, gradients refer to comoving coordinates, and the homogeneous background force has already been subtracted. Differentiating the first equation and eliminating $\mathbf{u}$ and $\varphi$ givesA Fourier mode with comoving wavenumber $k$ therefore obeysThe Jeans wavenumber occurs where the coefficient in parentheses vanishes. Since a comoving wavelength is $2\pi /k$, the result isBelow this comoving Jeans length, pressure restores a displaced fluid element: there are acoustic waves rather than growing Jeans instability. In a static background they oscillate; expansion changes their amplitude and frequency, and dissipation can damp them. Pressure support alone does not erase them. In the pre-recombination photon-baryon fluid, Silk damping supplies a separate erasure mechanism.

For collisionless matter, replace the sound speed by a characteristic velocity dispersion. For an isotropic Maxwell-Boltzmann velocity distribution, defining ${\sigma }^2=⟨u_x^2⟩=⟨|\mathbf{u}{|}^2⟩/3$, the static collisionless marginal-stability calculation givesThe coefficient depends on the distribution and the convention for velocity dispersion. Small-scale suppression is now collisionless free streaming and phase mixing, rather than collisional acoustic waves. An instantaneous collisionless Jeans length must be distinguished from the accumulated comoving distance $\int vdt/a$ traveled since decoupling.

For the requested thermal histories, an Einstein-de Sitter universe describes the matter-era limit: it cannot literally also have a radiation-dominated epoch. Interpret the question as a spatially flat universe with negligible cosmological constant, passing from radiation domination, $a\propto t^{1/2}$, to matter domination, $a\propto t^{2/3}$. For the usual schematic dynamical Jeans scale use ${\bar{\rho }}_{\mathrm{tot}}$ in the estimate, soNear the particle horizon a relativistic perturbation calculation replaces the Newtonian derivation; the scaling estimate still identifies the relevant sound-crossing or streaming scale.

Before cosmological recombination, adiabatic initial conditions keep the entropy per baryon fixed in the tightly coupled photon-baryon fluid. With $R_b=3{\bar{\rho }}_b/(4{\bar{\rho }}_{\gamma })\propto a$, its photon-baryon sound speed isFor example, this follows from $\delta {\rho }_{\gamma }/{\bar{\rho }}_{\gamma }=(4/3)\delta {\rho }_b/{\bar{\rho }}_b$ and $\delta P=c^2\delta {\rho }_{\gamma }/3$. In the radiation-dominated limit $R_b\ll 1$, $c_s\simeq c/\sqrt{3}$; in the strongly baryon-loaded matter-era limit $R_b\gg 1$, $c_s\propto a^{-1/2}$. Thus the baryon Jeans length across recombination has the asymptotic historyThe transition to the plateau is smooth: retaining $1+R_b$ gives ${\lambda }_{J,b}^{\mathrm{com}}\propto a^{1/2}/\sqrt{1+R_b}$ in matter domination. If baryon loading is still small, this intermediate segment instead rises as $a^{1/2}$ until loading becomes important.

At photon decoupling, radiation ceases to provide pressure support to the baryons, so the sound speed drops from the coupled-fluid value to that of a nonrelativistic ideal gas, $c_{s,b}^2=\gamma k_B{T}_b/(\mu m_p)$, with $\gamma =5/3$. The drop in comoving Jeans length at fixed density is the same ratio of sound speeds, typically several orders of magnitude. Subsequently an isentropic process gives $T_b\propto {\bar{\rho }}_b^{\gamma -1}\propto a^{-2}$, so $c_{s,b}\propto a^{-1}$ and ${\lambda }_{J,b}^{\mathrm{com}}\propto a^{-1/2}$. Residual energy exchange through Compton scattering can initially keep $T_b\propto a^{-1}$, producing a short approximately flat gas segment before the decreasing segment. The idealized sketch assumes immediate adiabatic cooling: