 Solution

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.