Constant-loading acoustic model ($y^{\prime \prime }+k^2{c}_s^{2}y=-k^2{c}_s^{2}R$)

A deliberately restricted forced harmonic oscillator shows two effects of the baryon loading parameter without duplicating the full photon Boltzmann hierarchy. Use conformal time $d\eta =dt/a$ in seconds, a comoving wavenumber $k$ in inverse metres, constant $R=3{\rho }_b/(4{\rho }_{\gamma })$, and equal constant dimensionless Newtonian gauge gravitational potentials $\Phi =\Psi \ne 0$. The energy densities ${\rho }_b$ and ${\rho }_{\gamma }$ use the same units. Neglect expansion drag, Silk damping and evolution of the potentials. The photon-baryon sound speed is $c_s=c/\sqrt{3(1+R)}$. Writing ${\Theta }_0$ for the photon monopole temperature perturbation, the photon-baryon acoustic oscillator reduces toFor $y=({\Theta }_0+\Psi )/\Psi$ this gives the displayed forced equation. Constant loading and the neglected expansion drag make this a local teaching model, not a self-consistent expanding cosmology with $R$ fixed forever.

Choose $y(0)=1/3$ and $y^{\prime }(0)=0$, the static matter-era Sachs-Wolfe combination normalized by $\Psi$. With $q=k{c}_s\eta$, direct substitution givesAt $R=0$, successive extrema have equal magnitudes $1/3$. At the stipulated $R=0.6$, $y(\pi )=-23/15$ whereas $y(2\pi )=1/3$. The shifted equilibrium $y=-R$ makes alternating extrema unequal. This is an elementary illustration of how the baryon loading parameter changes acoustic amplitudes; these extrema are neither measured peak heights nor a complete Cosmic microwave background power spectrum. Comparing curves at fixed $q$ isolates the response; at fixed $k$ and $\eta$, increasing $R$ also reduces the photon-baryon sound speed and hence the phase.

For example, holding $c{\eta }_{\ast }=300\mathrm{Mpc}$ gives sound horizons $r_s=173.2\mathrm{Mpc}$ at $R=0$ and $136.9\mathrm{Mpc}$ at $R=0.6$. Their ratio is $1/\sqrt{1.6}$. Each term of the forced equation has units of inverse seconds squared because $y$ and $R$ are dimensionless and $k{c}_s$ is an inverse time. In the actual photon-baryon fluid, $R$ grows with the scale factor and the sound horizon is the integral of the evolving photon-baryon sound speed.