Skip to main content
Sandbox Physics
Cosmology from expansion to inference

L2 · Chapter 3

Why primordial sound leaves a series of peaks

Before recombination, photon pressure pushed baryons out of compressed regions and gravity pulled them back into potential wells. Different wavelengths reached different phases by last scattering. Modes near compression or rarefaction extrema left peaks in the angular power spectrum.

Photons push and baryons add inertia

Frequent scattering couples photons to electrons, and electromagnetic interactions couple electrons to nuclei. In tight coupling the photon-baryon fluid moves almost together. Photon pressure supplies the restoring force; baryon mass resists acceleration. Define the loading R=3ρb/(4ργ)R=3\rho_b/(4\rho_\gamma). Both densities are energy per physical volume.

The sound speed is cs=c/3(1+R)c_s=c/\sqrt{3(1+R)}. More baryons lower it. Baryons and photons dilute at different rates, so RR grows as the universe expands. Sound therefore travels a different comoving distance in each interval.

rs(z∗)=∫0η∗cs(η) dηr_s(z_*)=\int_0^{\eta_*}c_s(\eta)\,d\eta
The sound horizon is the comoving distance sound travels before last scattering. Here conformal time retains time units and sound speed retains length per time.

Reading: Hu & Sugiyama

Solve one oscillator inside a fixed well

Let Θ0=δT/T\Theta_0=\delta T/T be the photon temperature monopole and use the Newtonian-gauge time potential Ψ\Psi. Set light speed to one. With constant potential and baryon loading, and omitting fluid-drag terms, the compression mode obeys the forced-oscillator equation below. Primes denote conformal-time derivatives.

Gravitational redshift contributes to the observed temperature, so define X=Θ0+ΨX=\Theta_0+\Psi. Moving the constant force into a shifted equilibrium gives an ordinary cosine when the initial velocity vanishes.

Θ0′′+k2cs2Θ0=−k23Ψ,X′′+k2cs2(X+RΨ)=0\Theta_0''+k^2c_s^2\Theta_0=-\frac{k^2}{3}\Psi,\qquad X''+k^2c_s^2(X+R\Psi)=0
In this equation sound speed is measured in units of light speed. The fixed-potential, fixed-loading approximation isolates the equilibrium shift.
X(η)+RΨ=[X(0)+RΨ]cos⁡(krs)X(\eta)+R\Psi=[X(0)+R\Psi]\cos(k r_s)
This is the zero-initial-velocity solution. The effective temperature oscillates about minus the baryon loading times the time potential.

Reading: Hu & Sugiyama

Calculate the two turning points

Choose a fixed well with Ψ=−10−5\Psi=-10^{-5}, R=1/2R=1/2, initial effective temperature X(0)=Ψ/3X(0)=\Psi/3 and zero initial velocity. The potential in the well is negative, so the baryon-induced equilibrium effective temperature is positive. Substitution into the cosine solution gives the dimensionless expression below.

After half an oscillation, krs=πkr_s=\pi and X=(4/3)×10−5X=(4/3)\times10^{-5}. After a full oscillation, X=−(1/3)×10−5X=-(1/3)\times10^{-5}. Their absolute amplitudes differ by a factor of four. This shows the effect of the shifted equilibrium for the specified initial conditions. A temperature spectrum also combines velocity, driving and projection, as the next section explains.

X10−5=12−56cos⁡(krs)\frac{X}{10^{-5}}=\frac12-\frac56\cos(kr_s)
This restricted example calculates the equilibrium shift and both turning-point amplitudes for the chosen initial conditions.

Baryons make compression and rarefaction unequal

Without baryon loading, XX oscillates around zero. Adding baryons shifts the equilibrium and deepens compression inside potential wells. For the usual adiabatic initial conditions, compression peaks rise relative to rarefaction peaks, producing the familiar odd-even temperature-peak pattern.

The complete temperature spectrum also contains Doppler velocities, changing potentials and photon diffusion. Increasing baryon density changes recombination, sound speed and the sound horizon together. A realistic comparison therefore needs the coupled calculation. The frozen CAMB grid retains these simultaneous effects.

Reading: CAMB source and documentation · Hu & Sugiyama

Estimate the separation between peaks

Modes with krs≃nπk r_s\simeq n\pi reach oscillator extrema at last scattering. Projection places them roughly at ℓ≃kDM\ell\simeq kD_M. Combining these relations gives the main acoustic spacing.

For an illustrative rs=144  Mpcr_s=144\;\mathrm{Mpc} and DM=14000  MpcD_M=14000\;\mathrm{Mpc}, the angle is about 0.0103  rad0.0103\;\mathrm{rad} and the acoustic multipole scale is about ℓA=305\ell_A=305. This estimates a spacing scale. Driving, phase shifts and projection move the first temperature peak away from that simple value.

θ∗=rsDM,ℓA=πθ∗,Δℓ≈ℓA\theta_* = \frac{r_s}{D_M},\qquad \ell_A=\frac{\pi}{\theta_*},\qquad \Delta\ell\approx\ell_A
Both lengths are comoving. The next chapter replaces the rough projection with the transfer-function integral.
ObserverEqual comoving sound horizons

Increasing comoving distance to the ruler →

θ∗=rs(z∗)/DM(z∗)\theta_* = r_s(z_*) / D_M(z_*)
Schematic geometry, with both lengths comoving. Temperature peaks also carry phase shifts and projection effects.

Try the prediction

Connect compression to a Boltzmann spectrum

  1. Locate the first three peaks in the reference preset and predict which respond most to more baryons.
  2. Change baryon density and compare with the fixed reference.
  3. Change late-time expansion separately to distinguish altered peak heights from a shift of the peak sequence.

Exercises and explanations

  1. At fixed sound horizon, increase transverse comoving distance to last scattering by ten percent. How do the acoustic angle and multipole scale change?

    Read the explanation

    The angle becomes ten-elevenths as large and the multipole scale rises by ten percent. The same physical pattern looks smaller and projects to higher multipoles.

  2. Why do peak positions alone fail to determine baryon density and late expansion uniquely?

    Read the explanation

    They mainly constrain a ratio of sound horizon to transverse distance. Different parameter combinations can preserve that ratio. Relative heights, damping, polarization, lensing and external distances provide additional information.

Sources and calculation scope

  1. Hu & Sugiyama, Small-Scale Cosmological Perturbations: An Analytic Approach, arXiv:astro-ph/9510117v2 (1996)

    Tight-coupling acoustic physics and baryon loading. The constant-potential, constant-baryon-loading oscillator in chapter three is a restricted analytic example; the interactive spectra come from the frozen CAMB calculation.

  2. CAMB source and documentation

    The acoustic-grid generator pins CAMB 1.6.6, saves its full cosmological settings and checks interpolation at off-grid points. The browser interpolates a fixed grid; it does not run a Boltzmann hierarchy.

FeedbackSomething broken, confusing, or missing on this page?Tell us
What kind of feedback?

Sent with this page’s address, its Lab ID, the site language and your device type. Please leave out personal details. Feedback is kept for up to 12 months.Privacy policy