Lyman-Alpha Forest

Lyman-alpha forest (Lyα forest) is a dense array of hundreds of narrow absorption lines appearing in the optical and ultraviolet spectra of high-redshift quasars (z > 2) and distant galaxies, produced by resonant Lyman-alpha (1s → 2p) electronic transitions of neutral atomic hydrogen (H I, rest-frame wavelength λ₀ = 1215.6701 Å, photon energy 10.20 eV) inside intervening filaments of the intergalactic medium (IGM) along the line of sight. First identified in 1970–1971 by the American astronomer Roger Lynds in optical spectra of the z = 2.877 quasar 4C 05.34, each individual absorption trough in the "forest" corresponds to a cosmic web structure at an intermediate cosmological redshift z_abs < z_quasar, which shifts the rest-frame 1215.67 Å resonance line to an observed wavelength λ_obs = 1215.67 Å × (1 + z_abs) located strictly on the short-wavelength (blue) side of the quasar's own broad Lyα emission peak. In modern physical cosmology, the Lyman-alpha forest serves as a one-dimensional "core sample" through the cosmic web—mapping the ionization state of the universe after cosmic reionization, testing warm dark matter and neutrino free-streaming on megaparsec scales, and anchoring high-redshift baryon acoustic oscillations at 2.0 < z < 3.8.

When continuum ultraviolet photons emitted by a distant quasar at redshift z_em travel across billions of light-years toward Earth, cosmic expansion stretches their wavelengths continuously. Whenever a photon's local rest-frame wavelength passes through λ₀ = 1215.67 Å at some intermediate redshift z_abs, any neutral hydrogen atom (H I) present at that location can absorb the photon via the electric-dipole-allowed 1s ²S_1/2 → 2p ²P_{1/2, 3/2} transition, governed by the large oscillator strength f₁₂ = 0.4162 and integrated cross section σ_α = (πe² / m_ec) f₁₂ = 1.10 × 10⁻² cm² Hz.

In 1965, James E. Gunn and Bruce A. Peterson calculated the homogeneous optical depth τ_GP(z) that would be produced by a smooth intergalactic medium with neutral hydrogen number density n_{H I}(z):

• Gunn–Peterson Optical Depth Formula: τ_GP(z) = [ (πe² f₁₂ λ₀) / (m_e c H(z)) ] · n_{H I}(z) ≈ 4.6 × 10⁵ · x_{H I} · (Ω_bh² / 0.0224) · [(1 + z) / 4]^3/2, where x_{H I} = n_{H I} / n_H is the neutral hydrogen fraction.
• Implication for the Ionization State of the Universe: Because the prefactor is ~4.6 × 10⁵, even a tiny neutral fraction of x_{H I} = 10⁻⁴ (0.01% neutral hydrogen) produces an optical depth τ_GP ≈ 46, which completely extinguishes all transmitted quasar flux (e^−τ → 0) across a continuous Gunn–Peterson trough. The fact that quasar spectra at 2 < z < 5 transmit 30% to 85% of their continuum flux between discrete Lyα forest lines proves that the intergalactic medium at z < 5.5 is highly photoionized (x_{H I} ≈ 10⁻⁵ to 10⁻⁴, or 99.99% ionized H II protons and free electrons) by the extragalactic ultraviolet background (UVB) emitted by star-forming galaxies and quasars.
• Detection of the Reionization Boundary (z ≈ 5.3–6.3): In 2001, Robert H. Becker and the Sloan Digital Sky Survey collaboration discovered complete Gunn–Peterson troughs in quasars at z > 6.0 (such as SDSS J1030+0524 at z = 6.28), marking the tail end of the Epoch of Reionization when the last neutral islands of intergalactic hydrogen were ionized.

Spectroscopists classify intergalactic and circumgalactic hydrogen absorbers into four distinct physical regimes based on their neutral hydrogen column density N_{H I} (atoms cm⁻²), which controls the Voigt line profile (combining a Gaussian thermal Doppler core of width b = √(2k_BT / m_p) ≈ 15–35 km s⁻¹ with Lorentzian radiation-damping wings):

In the 1970s and 1980s, astronomers initially modeled Lyα forest lines as isolated spherical "clouds" confined by external pressure or dark matter minihalos. In 1994–1998, cosmological hydrodynamic simulations of Cold Dark Matter (ΛCDM) by Renyue Cen, Jordi Miralda-Escudé, Jeremiah P. Ostriker, Lars Hernquist, Neal Katz, and David H. Weinberg overturned the cloud model: they demonstrated that the entire diffuse Lyα forest (N_{H I} < 10¹⁵ cm⁻²) arises naturally from the continuous, fluctuating density field of the photoionized cosmic web.

In photoionization equilibrium, the volumetric rate of neutral hydrogen photoionization by the UV background (Γ_UVB n_{H I}) balances radiative proton–electron recombination (α_rec(T) n_p n_e), where the recombination coefficient scales with gas temperature as α_rec(T) ∝ T^−0.72. Simultaneously, competition between photoionization heating and adiabatic Hubble expansion cooling establishes a tight power-law temperature–density relation (equation of state) across the diffuse IGM (T = T₀ Δ^γ−1, with T₀ ≈ 10⁴ K and adiabatic index γ ≈ 1.4–1.6). Combining these two relations yields the Fluctuating Gunn–Peterson Approximation (FGPA):

• FGPA Optical Depth Scaling: τ(z, x) ∝ n_{H I} ∝ n_H² T^−0.72 / Γ_UVB ∝ Δ^{2 − 0.72(γ − 1)} ≈ Δ^1.6.
• Because optical depth scales directly as a smooth power of the local baryonic overdensity Δ = ρ_b / ρ̄_b, a high-resolution quasar spectrum (Keck HIRES or VLT UVES) acts as a direct nonlinear map of the underlying one-dimensional dark-matter density field along the sightline.

Because the diffuse intergalactic medium at z = 2 to 5 remains in the mildly nonlinear regime (Δ ≲ 10) and is largely uncontaminated by complex galaxy-scale supernovae or AGN feedback, the Lyman-alpha forest is one of precision cosmology's sharpest laboratories:

1. One-Dimensional Flux Power Spectrum P_1D(k) — Warm Dark Matter & Neutrino Bounds: Whereas galaxy surveys at z < 1 struggle to measure the linear matter power spectrum on small scales (k > 0.5 h Mpc⁻¹) due to virialized halo nonlinearities, the 1D Lyα flux power spectrum probes comoving scales down to 1 to 10 Mpc (k ≈ 1 to 10 h Mpc⁻¹). If dark matter consisted of Warm Dark Matter (WDM) (such as sterile neutrinos or gravitinos) with a long thermal free-streaming length, small-scale filaments in the Lyα forest would be smoothed out and erased. High-resolution Lyα forest spectra rule out thermal relic warm dark matter particles lighter than m_WDM ≥ 5.3 keV (at 95% CL) and constrain the sum of active neutrino masses to ∑m_ν < 0.09–0.11 eV when combined with CMB data.
2. Three-Dimensional Baryon Acoustic Oscillations at z ≈ 2.33 (BOSS / eBOSS / DESI): Pioneered by the SDSS-III BOSS survey (Busca et al., 2013; Slosar et al., 2013) and dramatically expanded by DESI (measuring > 700,000 high-redshift quasar spectra), cosmologists correlate the Lyα flux transmission δ_F across hundreds of thousands of adjacent quasar sightlines—and cross-correlate Lyα forest skewers with the three-dimensional positions of quasars—to measure the 147 Mpc baryon acoustic oscillations peak at effective redshift z_eff = 2.33, directly measuring the cosmic deceleration-to-acceleration transition during the matter-dominated era.
3. 3D IGM Tomography and the He II Lyα Forest: Using dense grids of background Lyman-break galaxies alongside space ultraviolet spectroscopy (HST COS), astronomers construct three-dimensional tomographic maps of z ~ 2.5 protoclusters and observe the singly ionized helium Lyα forest (He II 303.78 Å), which completed its second reionization (He II → He III) driven by hard quasar photons around z ≈ 2.8–3.0.

1. Lynds, R. (1971). The absorption-line spectrum of 4C 05.34. The Astrophysical Journal, 164, L73–L78. https://doi.org/10.1086/180695
2. Gunn, J. E., & Peterson, B. A. (1965). On the density of neutral hydrogen in intergalactic space. The Astrophysical Journal, 142, 1633–1636. https://doi.org/10.1086/148444
3. Rauch, M. (1998). The Lyman alpha forest in the spectra of quasistellar objects. Annual Review of Astronomy and Astrophysics, 36(1), 267–316. https://doi.org/10.1146/annurev.astro.36.1.267
4. McQuinn, M. (2016). The evolution of the intergalactic medium. Annual Review of Astronomy and Astrophysics, 54, 313–362. https://doi.org/10.1146/annurev-astro-082214-122355