The Population of Damped Lyman-alpha and Lyman Limit Systems in the Cold Dark Matter Model

Lyman limit and damped Lyman-alpha absorption systems probe the distribution of collapsed, cold gas at high redshift. Numerical simulations that incorporate gravity and gas dynamics can predict the abundance of such absorbers in cosmological models. We develop a semi-analytical method to correct the numerical predictions for the contribution of unresolved low mass halos, and we apply this method to the Katz et al. (1996) simulation of the standard cold dark matter model ($\Omega=1$, $h=0.5$, $\Omega_b=0.05$, $\sigma_8=0.7$). Using this simulation and higher resolution simulations of individual low mass systems, we determine the relation between a halo's circular velocity $v_c$ and its cross section for producing Lyman limit or damped absorption. We combine this relation with the Press-Schechter formula for the abundance of halos to compute the number of absorbers per unit redshift. The resolution correction increases the predicted abundances by about a factor of two at z=2, 3, and 4, bringing the predicted number of damped absorbers into quite good agreement with observations. Roughly half of the systems reside in halos with circular velocities $v_c>100\kms$ and half in halos with $35\kms=2$ only for column densities exceeding $N_{HI}=10^{22} cm^{-2}$.