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<v_c<100\kms$. Halos with $v_c>150\kms$ typically harbor two or more systems capable of producing damped absorption. Even with the resolution correction, the predicted abundance of Lyman limit systems is a factor of three below observational estimates, signifying either a failure of standard CDM or a failure of these simulations to resolve the systems responsible for most Lyman limit absorption. By comparing simulations with and without star formation, we find that depletion of the gas supply by star formation affects absorption line statistics at $z>=2$ only for column densities exceeding $N_{HI}=10^{22} cm^{-2}$. 
Jeffrey P. Gardner, Neal Katz, Lars Hernquist, and David H. Weinberg