The Population of Damped Lyα and Lyman Limit Systems in the Cold Dark Matter Model

Open Access
Publisher Website
Google Scholar
Abstract
Lyman limit and damped Lyα absorption systems probe the distribution of collapsed, cold gas at high redshift. Numerical simulations that incorporate gravity and gasdynamics can predict the abundance of such absorbers in cosmological models. We develop a semianalytical method to correct the numerical predictions for the contribution of unresolved low-mass halos, and we apply this method to the Katz et al. simulation of the standard cold dark matter model (Ω = 1, h = 0.5, Ωb = 0.05, and σ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 vc and its cross section for producing Lyman limit or damped Lyα absorption. We combine this relation with the Press-Schechter formula for the abundance of halos—itself calibrated against the simulated halo population—in order to compute the number of absorbers per unit redshift. The resolution correction increases the predicted abundances by about a factor of 2 at z = 2, 3, and 4, bringing the predicted number of damped Lyα absorbers into quite good agreement with observations. Roughly half of these systems reside in halos with circular velocities vc ≥ 100 km s-1, and half in halos with 35 km s-1 ≤ vc ≤ 100 km s-1. Halos with vc > 150 km s-1 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 3 below observational estimates, signifying either a failure of the standard cold dark matter model or a failure of these simulations to resolve most of the systems responsible for Lyman limit absorption. By comparing simulations with and without star formation, we find that the depletion of the gas supply by star formation affects absorption-line statistics at z ≥ 2 for column densities exceeding NH I = 1022 cm-2 only, even though half of the cold, collapsed gas has been converted to stars by z = 2.
All Related Versions