Connecting cosmological accretion to strong Ly α absorbers

We present an analytical model for the cosmological accretion of gas on to dark matter haloes, based on a similarity solution applicable to spherical systems. Performing simplified radiative transfer, we compute how the accreting gas turns increasingly neutral as it self-shields from the ionizing background, and obtain the column density, $N_{\rm H\,{\small I}}$ , as a function of impact parameter. The resulting column-density distribution function (CDDF) is in excellent agreement with observations. The analytical expression elucidates (1) why haloes over a large range in mass contribute about equally to the CDDF as well as (2) why the CDDF evolves so little with redshift in the range z = 2-5. We show that the model also predicts reasonable damped Lyman-αabsorber(DLA) line widths (v₉₀), bias, and molecular fractions. Integrating over the CDDF yields the mass density in neutral gas, $\Omega _{\rm H\,{\small I}}$ , which agrees well with observations. $\Omega _{\rm H\,{\small I}}(z)$ is nearly constant even though the accretion rate on to haloes evolves. We show that this occurs because the fraction of time that the inflowing gas is neutral depends on the dynamical time of the halo, which is inversely proportional to the accretion rate. Encapsulating results from cosmological simulations, the simple model shows that most Lyman-limit systems and DLAs are associated with the cosmological accretion of gas on to haloes.