Skip to content

Power Spectra, Horizon Crossing, and Freeze-Out

A power spectrum is a late-time two-point function, whereas “horizon crossing” names an intermediate epoch. Replacing evolution by evaluation at k=aHk=aH is accurate only when background parameters vary slowly and the relevant long-wavelength mode rapidly approaches a constant attractor solution.

Required background. Mukhanov–Sasaki scalar modes supplies scalar normalization; tensor modes supplies both helicities; and vacuum choice and initial-state effects supplies the state condition.

Helpful background. Mode matching across cosmological eras supplies regulated junction and transfer methods.

Define the scalar dimensionless spectrum by

ζkζk=(2π)3δ3(k+k)2π2k3Pζ(k).\langle\zeta_{\mathbf k}\zeta_{\mathbf k'}\rangle =(2\pi)^3\delta^3(\mathbf k+\mathbf k') \frac{2\pi^2}{k^3}\mathcal P_\zeta(k).

For a slowly varying single-clock attractor with an adiabatic initial state,

Pζ(k)=Hs28π2ϵscs,sMPl2[1+O(ϵ,ηH,s)],csk=aHat ts,\mathcal P_\zeta(k)= \frac{H_s^2}{8\pi^2\epsilon_s c_{s,s}M_{\rm Pl}^2} \left[1+O(\epsilon,\eta_H,s)\right], \qquad c_sk=aH\quad\text{at }t_s,

where ηH=ϵ˙/(Hϵ)\eta_H=\dot\epsilon/(H\epsilon) and s=c˙s/(Hcs)s=\dot c_s/(Hc_s). The tensor spectrum summed over the two helicities is

Pt(k)=2Ht2π2MPl2[1+O(ϵ)],k=aHat tt.\mathcal P_t(k)=\frac{2H_t^2}{\pi^2M_{\rm Pl}^2} \left[1+O(\epsilon)\right], \qquad k=aH\quad\text{at }t_t.

The scalar and tensor crossing times differ when cs1c_s\ne1. At leading order in the minimal theory, transporting each mode to its constant regime gives r=Pt/Pζ=16ϵcsr=\mathcal P_t/\mathcal P_\zeta=16\epsilon c_s, with the parameters evaluated consistently. The systematic slow-roll expansion and its evaluation prescription are derived in Stewart and Lyth 1993, §§II–III, Eqs. (20)–(42).

Spectral tilts are derivatives of the late-time spectra,

ns1=dlnPζdlnk,nt=dlnPtdlnk.n_s-1=\frac{d\ln\mathcal P_\zeta}{d\ln k}, \qquad n_t=\frac{d\ln\mathcal P_t}{d\ln k}.

At lowest slow-roll order, ns1=2ϵηHsn_s-1=-2\epsilon-\eta_H-s and nt=2ϵn_t=-2\epsilon. These are model- and hierarchy-dependent truncations, not exact symmetry identities.

The canonical frequency ωk2=cs2k2z/z\omega_k^2=c_s^2k^2-z''/z changes character near the sound horizon. The turning region sets the amplitude, but the superhorizon solution

ζ=C1+C2dηz2\zeta=C_1+C_2\int\frac{d\eta}{z^2}

determines whether it freezes. A reliable first application compares three objects for constant slow roll: the exact Hankel solution at late time, its value at csk=aHc_sk=aH, and the slow-roll crossing formula. Their difference is the controlled correction that the star notation compresses.

A sharp transition between eras should be treated as a limit of smooth backgrounds. Integrating the action through the transition preserves ζ\zeta and its canonical momentum z2ζz^2\zeta' unless a surface operator supplies a jump. Matching an arbitrarily chosen field and its ordinary derivative can violate the constraint system.

For a smooth feature, evolution can be summarized by a transfer matrix acting on the constant and second long-wavelength solutions. Its determinant is fixed by the conserved symplectic product, while its off-diagonal entries quantify mode conversion. This is more informative than assigning a single “crossing time”: a feature can occur after csk=aHc_sk=aH and still convert the second solution into the final constant mode. Numerical transport should therefore continue until both dlnP/dNd\ln\mathcal P/dN and the second-mode contribution fall below declared tolerances.

The same distinction matters for observational pivot scales. Parameters quoted at a pivot kk_* are background quantities evaluated through a model-dependent mapping between kk_* and the number of e-folds before the end of inflation. Reheating history changes that mapping without changing the primordial mode equation. A power-spectrum prediction should keep mode evolution, background-to-pivot mapping, and later transfer functions as separate uncertainty sources.

The structure map places spectrum evaluation after state normalization and before interaction corrections.

A normalized subhorizon mode crosses its effective turning region and evolves through constant and second superhorizon solutions before defining the late-time spectrum

Crossing estimates the amplitude near the turning region; freeze-out is the later dynamical suppression of the second long-wavelength solution. Schematic; not to scale.

In ultra-slow-roll evolution, ϵa6\epsilon\propto a^{-6} and the nominal second solution grows as a3a^3. Evaluating Pζ\mathcal P_\zeta at k=aHk=aH then misses the subsequent enhancement. The adversarial check evolves the exact mode through the end of the nonattractor phase and verifies continuity of the canonical data. If the late amplitude depends on the duration or exit profile, it must be reported as a transfer calculation rather than a universal horizon-crossing formula.

An excited initial state gives a different failure signature: it changes oscillatory phase and amplitude already before crossing but need not spoil later conservation. Varying the state and the attractor history independently distinguishes an initial-condition correction from superhorizon transport.

See the chapter’s domain and failure conditions. The validity map distinguishes slow-roll truncation error, matching error, state error, and genuine superhorizon evolution.

Rapid background change, an excited state, incorrect junction data, or a growing superhorizon mode separates horizon crossing from the final power spectrum

The crossing formula fails when slow variation, adiabatic state preparation, canonical matching, or attractor freeze-out fails; exact mode transport is then required. Schematic; not to scale.

  • Stewart, E. D., and D. H. Lyth, “A More Accurate Analytic Calculation of the Spectrum of Cosmological Perturbations Produced during Inflation,” Physics Letters B 302, 171–175 (1993), doi:10.1016/0370-2693(93)90222-S.