ADAPTIVE MESH REFINEMENT SIMULATIONS OF THE IONIZATION STRUCTURE AND KINEMATICS OF DAMPED Lyα SYSTEMS WITH SELF-CONSISTENT RADIATIVE TRANSFER

摘要

We use high-resolution Eulerian hydrodynamics simulations to study kinematic properties of the low-ionization species in damped Lyα systems at redshift z = 3. Our adaptive mesh refinement simulations include most key ingredients relevant for modeling neutral gas in high column density absorbers: hydrodynamics, gravitational collapse, continuum radiative transfer above the hydrogen Lyman limit, and gas chemistry, but no star formation. We model high-resolution Keck spectra with unsaturated low ion transitions in two Si II lines (1526 and 1808 A) and compare simulated line profiles to the data from the SDSS DLA survey. We find that with increasing grid resolution the models show a trend in convergence toward the observed distribution of H I column densities. While in our highest resolution model we recover the cumulative number of DLAs per unit absorption distance, none of our models predicts DLA velocity widths as high as indicated by the data, suggesting that feedback from star formation might be important. Atz = 3 a nonnegligible fraction of DLAs with column densities below 10~(21) cm~(-2) is caused by filamentary structures in more massive halo environments. Lower column density absorbers with N_(H I) < 10~(21.4) cm~(-2) are sensitive to changes in the UV background, resulting in a 10% reduction of the cumulative number of DLAs for twice the quasar background relative to the fiducial value and nearly a 40% reduction for 4 times the quasar background. We find that the mass cutoff below which a large fraction of dwarf galaxies cannot retain gas after reionization is ～7 x 10~7 solar mass, lower than the previous estimates. Finally, we show that models with self-shielding commonly used in the literature produce significantly lower DLA velocity widths than the full radiative transfer runs, which essentially renders these self-shielded models obsolete.