Neutral-hydrogen column-density distribution ($f(N,z)={\partial }^2\mathcal{N}/(\partial N\partial z)$)

The neutral-hydrogen column-density distribution counts Lyman-alpha absorption systems per unit column density and redshift. A count in logarithmic column bins is $N\ln (10)f(N,z)$. At low forest columns an approximate power law $f\propto N^{-\beta }$ with $\beta$ near $1.5$ is useful, but the slope and normalization depend on column range and epoch. Lyman limit systems are optically thick at the hydrogen photoionization edge; damped Lyman-alpha systems have strong radiative damping wings. A distribution per absorption distance instead removes the geometrical path factor.

A neutral-hydrogen column density $N$ traces total hydrogen column $N/x_{\mathrm{HI}}$ only after specifying its neutral fraction $x_{\mathrm{HI}}$. Integrating the latter against the distribution per absorption distance gives the comoving gas mass density in the displayed expression. The factor ${\mu }_H$ includes helium if desired. In photoionization equilibrium for optically thin, almost fully ionized hydrogen, $x_{\mathrm{HI}}\simeq \alpha (T)n_H/\Gamma$, so gas density, temperature, photoionization rate and absorber size are required. An observed distribution alone does not fix the total baryon inventory.

Let $f_X={\partial }^2\mathcal{N}/(\partial N\partial X)$ count neutral hydrogen absorbers per column density and absorption distance. The summed column per unit proper light-path length is the mean neutral number density. Since $d\ell =cdz/[(1+z)H(z)]$ and $dX/dz=H_0(1+z{)}^2/H(z)$, that density is $(H_0/c)(1+z{)}^3\int N{f}_{X}dN$. Dividing the proper mass density by $(1+z{)}^3$ gives the comoving density in the displayed formula. The denominator is today's critical density; including neutral helium requires a stated mass correction. For a distribution per redshift, replace $H_0{f}_X$ by $H(z)f_z/(1+z{)}^2$.