Past exam of the mathematics course of the University of Cambridge / 2004 / iii / Paper 67 / 4

Define neutral-hydrogen column density along a ray by $N=\int n_{\mathrm{HI}}d\ell$, in atoms per unit area. The neutral-hydrogen column-density distribution per redshift isso $f_{z}dNdz$ is the expected number of absorbers in those intervals along one sightline. A logarithmic column bin instead has distribution $N\ln (10)f_z$; the bin convention must be specified.

Let $n_c$ be a constant comoving number density and ${\sigma }_p$ a constant physical interception area. The corresponding physical number density is $n_p=n_c(1+z{)}^3$. Proper light-path length satisfies $d\ell =c|dt|=cdz/[(1+z)H(z)]$. The expected intersection count, or first-order probability in infinitesimal $dz$, is thereforeThe Friedmann equation, the constant-equation-of-state density scaling for pressureless matter, and the present cosmological density parameters giveCombining them yields the absorber incidence and comoving number density relationHere the constant cross-section is physical, as required for this result; a fixed comoving cross-section would instead decrease in physical area as $(1+z{)}^{-2}$. It is convenient to define absorption distance by $dX/dz=(H_0/H(z))(1+z{)}^2$ and $f_X=f_z/(dX/dz)$, so a population with these non-evolution assumptions has constant incidence per $X$.

A useful description of the low-column Lyman-alpha forest at $z\simeq 3$ is a power law $f_X(N,z)=A(z)N^{-\beta }$, with $\beta$ around $1.4$–$1.5$ over the ordinary unsaturated forest range. The same column dependence holds for $f_z$ at a fixed redshift. It is not a universal unbroken power law: saturation, increased neutral fraction and changing absorber populations produce breaks. Lyman limit systems become optically thick at $N\simeq 10^{17.2}{\mathrm{cm}}^{-2}$ because $N{\sigma }_{912}\simeq 1$ with ${\sigma }_{912}\simeq 6.3\times 10^{-18}{\mathrm{cm}}^2$. Damped Lyman-alpha systems have $N\gtrsim 2\times 10^{20}{\mathrm{cm}}^{-2}$ and strong radiative damping wings. They are rare in number but contain substantial neutral mass. The highest-column distribution steepens or cuts off rather than extending indefinitely.

This sketch shows representative shapes and transition scales, not an observational fit; its ordinate is distribution per unit column, not counts per logarithmic bin. In an Einstein-de Sitter universe, $H=H_0(1+z{)}^{3/2}$, so a fixed comoving population with fixed physical cross-section above a fixed column threshold has $d\mathcal{N}/dz\propto (1+z{)}^{1/2}$. ThusThe observed redshift evolution of the Lyman-alpha forest is much stronger for commonly selected forest lines. For example, Kim, Cristiani and D'Odorico's absorption survey finds an index about $2.19$ in $d\mathcal{N}/dz\propto (1+z{)}^{\gamma }$ for $10^{13.64}<N/{\mathrm{cm}}^{-2}<10^{16}$ over $1.5<z<4$, corresponding to an increase by approximately $2^{2.19}\simeq 4.56$. This number is selection-dependent. Expansion changes the gas density, and the evolving cosmic ionizing background changes its neutral fraction through photoionization equilibrium, moving systems across a fixed neutral-column threshold. Structure also evolves. A non-evolving population with fixed physical cross-section is therefore not an adequate physical model of the observed forest.

Use absorption distance $X$ and its distribution $f_X$. Summing hydrogen columns per unit $X$ gives $\int N{f}_{X}dN$. The physical light-path length per $X$ isThe mean neutral-hydrogen number density is therefore $(H_0/c)(1+z{)}^3\int N{f}_{X}dN$. Divide by $(1+z{)}^3$ to obtain its comoving density, multiply by hydrogen atomic mass, and divide by today's mass critical density ${\rho }_{\mathrm{crit},0}=3H_0^2/(8\pi G)$. This derives the neutral gas density from an absorption distribution:If “neutral gas” includes helium, replace $m_H$ by ${\mu }_H{m}_H$, with a stated abundance correction approximately ${\mu }_H\simeq 1.3$. Equivalently, for the distribution per redshift the prefactor inside the integral is $m_{H}H(z)/[c{\rho }_{\mathrm{crit},0}(1+z{)}^2]$ multiplying $N{f}_z$. In the PDF's energy-density convention, both the mass density and critical density are multiplied by $c^2$, leaving their ratio unchanged. The usual quoted abundance uses comoving mass density relative to today's critical density; a ratio of physical density to the epoch's critical density instead multiplies this answer by $(1+z{)}^3{H}_0^2/H(z{)}^2$.

To infer total gas, one needs the neutral fraction $x_{\mathrm{HI}}=N_{\mathrm{HI}}/N_H$ for each absorber. The ionized gas density from an absorption distribution isThis is an ionization-corrected integral, not an inference available from the observed columns alone. For optically thin, almost fully ionized hydrogen, photoionization equilibrium gives $\Gamma n_{\mathrm{HI}}=\alpha (T)n_e{n}_{\mathrm{HII}}\simeq \alpha (T)n_H^2$. If the absorbing depth $L$ is specified, $N=\alpha (T)n_H^{2}L/\Gamma$ andThus photoionization rate, temperature and depth are all needed.

One useful physical closure sets $L$ to a local Jeans length, $L\propto (f_{g}T/n_H{)}^{1/2}$, with local gas fraction $f_g$. Solving $N\propto [\alpha (T)/\Gamma ](f_{g}T{)}^{1/2}{n}_H^{3/2}$ then gives the Jeans-length ionization correction for the Lyman-alpha forest:This illustrates why weak, highly ionized forest lines can dominate the total gas inventory even though damped Lyman-alpha systems dominate the neutral-column-weighted inventory. It is a model for locally confined optically thin gas; optically thick Lyman limit systems and damped Lyman-alpha systems require a separate radiative transfer treatment. The physical basis and observational application of this closure are developed in Schaye's absorption-based matter-density calculation.

No measured normalization of the neutral-hydrogen column-density distribution, ionizing field or gas temperature is supplied, so the paper's formulas do not determine numerical densities on their own. With observational data and suitable ionization corrections, the high-redshift Lyman-alpha forest gives a total gas density of order a few $10^{-2}$ relative to today's critical density, comparable to the cosmological baryon density ${\Omega }_b\sim 0.04$–$0.05$ for $H_0\sim 70\mathrm{km}{\mathrm{s}}^{-1}\mathrm{Mp}{\mathrm{c}}^{-1}$. The conclusion is that a large fraction of baryons at $z\simeq 3$ reside in diffuse ionized gas, rather than that every baryon has been individually accounted for. Schaye's analysis and Kim and collaborators' survey support that interpretation under their stated gas and ionizing-background assumptions.

The neutral-gas reservoir traced predominantly by damped Lyman-alpha systems is much smaller, of order $10^{-3}$ in the same convention. It is broadly comparable in order of magnitude to today's luminous stellar matter, whose cosmic stellar mass density is a few $10^{-3}$, but exact comparisons depend on cosmology, helium correction and which stellar components are counted. Storrie-Lombardi, McMahon and Irwin's neutral-gas census finds especially close comparison with present stellar disks. Most high-redshift gas is ionized; the smaller neutral reservoir has a mass of roughly the same order as today's visible stellar population. This comparison does not imply that the measured neutral reservoir alone simply becomes all present-day stars: gas accretion, recycling, ionization changes and outflows continue during galaxy evolution.

The neutral gas density from an absorption distribution is comparatively robust for high columns where the hydrogen is mostly neutral, but saturated lines, dust-obscured quasars, incomplete surveys and small numbers of rare high-column systems affect the integral. The ionized gas density from an absorption distribution has additional systematic uncertainty in the cosmic ionizing background, temperature, absorber depth, helium treatment and small-scale density inhomogeneity. At low columns, line blending and incompleteness matter; at high columns, attenuation makes an optically thin ionization correction invalid. A single power law extrapolated outside the measured column range can also produce a misleading mass integral. Hot gas, molecular gas and gas missed by the selected absorption population are not included automatically. These are observational and modelling qualifications to a baryon census, not failures of the column-mass accounting itself.