Past exam of the mathematics course of the University of Cambridge / 2005 / iii / Paper 74 / 2

The Lyman-alpha forest is the large collection of narrow Lyman-alpha absorption features in the spectrum of a distant quasar. Neutral hydrogen in foreground gas absorbs at its own redshifted $1s\to 2p$ resonance, so an absorber at $z_{\mathrm{abs}}$ places its feature atOrdinary intervening absorbers have $z_{\mathrm{abs}}<z_Q$, putting their Lyman-alpha lines blueward of the quasar's own broad Lyman-alpha emission. Many different absorbing patches along the same line of sight make the forest; they are not different transitions of one cloud. At sufficiently short wavelengths higher Lyman-series transitions can overlap the forest, so an actual analysis identifies these blends rather than treating every dip as isolated Ly-alpha.

The upper panel is an original schematic spectrum for $z_Q=3$, with absorption features at smaller cosmological redshifts and a broad emission feature near $4863$ Å. The lower panel anticipates the column-density discussion: both curves are teaching sketches, not measured spectra or fitted survey distributions. The ordinary Lyman-alpha forest mainly traces the fluctuating, photoionized intergalactic medium, rather than a collection of identical isolated, neutral clouds.

Hydrogen is the most abundant element, and its Lyman-alpha line is a strong resonance from the ground state. Very small neutral fractions in an otherwise ionized intergalactic medium therefore produce detectable Lyman-alpha absorption. A long quasar sightline samples a great many modest density enhancements in the cosmic gas, with a much larger covering area than the compact high-column regions that produce damped Lyman-alpha systems.

Heavy-element absorption additionally requires metal enrichment and enough ions in the relevant ionization state. Many weak Lyman-alpha forest absorbers contain too few such ions to produce detectable metal lines at the same sensitivity, while high-column Lyman limit systems and damped Lyman-alpha systems are intrinsically rarer. Thus abundant hydrogen, a sensitive resonance and the large interception area of diffuse gas together explain the forest's numerical dominance. The comparison is between absorber classes at specified observational sensitivity; it does not mean that a metal-rich system has only one spectral transition.

Typical weak Lyman-alpha forest absorbers have neutral hydrogen column densities roughly $10^{12}$–$10^{16}{\mathrm{cm}}^{-2}$, with stronger forest absorption extending towards $10^{17}{\mathrm{cm}}^{-2}$. Common Doppler broadening parameters are a few tens of ${\mathrm{km}\mathrm{s}}^{-1}$, consistent with gas at approximately $10^4$ K together with bulk velocity gradients. The thermal contribution isA line width is consequently a temperature upper bound if all other broadening is neglected, not a unique thermometer. The gas commonly lies in diffuse sheets and filaments at modest density contrast; only the strongest systems need substantial self-shielding or a galactic environment.

For optically thin gas illuminated by the cosmic ionizing background, photoionization equilibrium requireswhere the last approximation uses almost complete hydrogen ionization. For example, $n_H=10^{-5}{\mathrm{cm}}^{-3}$, $\alpha \simeq 4\times 10^{-13}{\mathrm{cm}}^3{\mathrm{s}}^{-1}$ and ${\Gamma }_{\mathrm{HI}}=10^{-12}{\mathrm{s}}^{-1}$ give $x_{\mathrm{HI}}\simeq 4\times 10^{-6}$. A $100$ kpc path then has $N_{\mathrm{HI}}\sim 10^{13}{\mathrm{cm}}^{-2}$ despite its tiny neutral fraction. This illustrates the Lyman-alpha forest photoionization scaling: low density makes the two-body radiative recombination rate weak relative to the one-body photoionization rate, but a long path and a strong resonance still produce measurable absorption.

Independent information is important because a column density alone is not a volume density. Similar absorption in close pairs of quasars indicates substantial transverse coherence, while an ionizing optical depth $N_{\mathrm{HI}}{\sigma }_{912}\ll 1$ shows that these weak systems do not shield themselves from the background. The continuum-shielding threshold is $1/{\sigma }_{912}\simeq 1.6\times 10^{17}{\mathrm{cm}}^{-2}$. Together with the observed line widths and ionization balance, these facts favor extended, low-density, highly ionized gas over small cold neutral clouds. The physical interpretation and characteristic scales are consistent with Rauch's forest analysis and the measured column and line-width distributions.

Specify the path convention before integrating the neutral-hydrogen column-density distribution. Write $f_z={\partial }^2\mathcal{N}/(\partial N\partial z)$ and $f_X={\partial }^2\mathcal{N}/(\partial N\partial X)$, with absorption distance satisfyingThe lower panel of the preceding sketch shows the approximate decrease of $f$ with $N=N_{\mathrm{HI}}$ near $z=3$. Ordinary forest columns are roughly described by $f\propto N^{-\beta }$, $\beta \sim 1.5$, but the slope changes through partially shielded Lyman limit systems and high-column damped Lyman-alpha systems. No single slope should be extrapolated through all these regimes. In logarithmic bins the number is $N\ln (10)f$, whereas their neutral mass contribution is proportional to $N^{2}f$; confusing these two weightings would give a wrong baryon inventory.

To derive the neutral gas density from an absorption distribution, a redshift interval $dz$ has proper light-path length $d\ell =cdz/[(1+z)H]$. The sum of neutral hydrogen columns per path length estimates the mean proper number density:The proper mass density is $m_H{\overline{n}}_{\mathrm{HI}}$; division by $(1+z{)}^3$ converts it to a comoving mass density. Thus, with today's critical density ${\rho }_{\mathrm{crit},0}=3H_0^2/(8\pi G)$,A representative finite-sample estimator replaces the last integral by ${\sum }_j{N}_j/\Delta X$, with completeness and selection corrections. This defines a comoving inventory divided by today's critical density; a fraction of the instantaneous critical density has a different normalization. Neutral gas including accompanying helium is obtained by multiplying by a specified mass correction, approximately $1.3$; hydrogen-only ${\Omega }_{\mathrm{HI}}$ must not silently include it.

At $z\sim 3$ the observed neutral inventory is of order $10^{-3}$, far below the cosmological baryon density of order $0.04$–$0.05$ for $h\sim 0.7$. High-column damped Lyman-alpha systems contribute most of this neutral mass although weak forest lines dominate the counts; for $\beta =1.5$ the mass per logarithmic column interval increases as $N^{1/2}$ in the forest regime. The small neutral inventory therefore does not imply a small total gas inventory: the intergalactic medium contains a large ionized baryon reservoir invisible to a neutral-only count. Recovering that reservoir requires an ionization correction and a gas-density model. The survey normalization and high-column contribution are described in the neutral-gas inventory analysis.

For a fixed column density or rest-equivalent-width threshold, the number of ordinary Lyman-alpha forest features per unit redshift generally increases towards higher redshift. Over much of $z\sim 1$–$4$, a rough description is $d\mathcal{N}/dz\propto (1+z{)}^{\gamma }$ with an index of order $2$–$3$, dependent on selection; the evolution is much slower at low redshift. Near reionization, blending and widespread opacity make a decomposition into individually detected lines increasingly inappropriate.

First separate physical evolution from geometrical path length. For comoving absorber density $n_c$ and proper interception area ${\sigma }_p$, absorber incidence and comoving number density givesEven a non-evolving population is not constant per unit redshift. Using absorption distance removes this particular geometrical effect.

The remaining redshift evolution of the Lyman-alpha forest reflects several competing processes. Cosmic expansion dilutes the gas; at fixed overdensity, photoionization equilibrium gives $n_{\mathrm{HI}}\propto \alpha (T)(1+z{)}^6/\Gamma$, so dilution can move many features below a fixed threshold. The cosmic ionizing background evolves as quasars, stellar sources and ionizing-photon absorption evolve; a decline of its intensity at late times partly offsets dilution. Photoheating changes temperature and radiative recombination, while gravitational instability redistributes gas into denser structures and hot shocked phases and changes its velocity field. Thus evolving line counts are not a direct count of conserved clouds; the density, thermal state, radiation field, geometry and observational selection must all be accounted for. The high-redshift increase and its physical interpretation are documented in the forest observations and models.

Near the background quasar, the quasar proximity effect typically reduces the number and strength of Lyman-alpha forest features relative to regions farther along the sightline. It is an ionization effect, not evidence that the quasar lies in an empty region. At fixed gas density and temperature, photoionization equilibrium givesIf the unperturbed neutral-hydrogen column-density distribution is $A{N}^{-\beta }$, $\beta >1$, detecting a column above $N_{\min }$ now requires its unperturbed column to exceed $(1+\omega )N_{\min }$. Integrating that distribution shows that the number above threshold is reduced by $(1+\omega {)}^{1-\beta }$. This supplies an explicit model for the decline of line counts.

The known quasar luminosity estimates its photoionization rate at proper separation $R$:These expressions assume unattenuated, steady, isotropic quasar emission and a specified ionizing spectrum. Fitting the measured suppression as a function of separation determines the ratio ${\Gamma }_Q/{\Gamma }_{\mathrm{bg}}$ and hence constrains the cosmic ionizing background. A weaker background gives a larger region in which the quasar dominates. For fixed background spectral shape, this constrains its specific intensity near the Lyman limit system ionization edge, not an arbitrary spectral function. The proximity-effect estimator of the ionizing background is biased if the host gas is overdense, the quasar's past luminosity differs from its observed luminosity, its radiation is anisotropic, or systemic-redshift errors distort the inferred separation; absorption between source and gas also changes ${\Gamma }_Q$. These effects explain why the method needs environmental and source modeling, as discussed in the proximity-effect measurements.