Damped Lyman-$\alpha$ absorbers and atomic hydrogen in galaxies: the view of the GAEA model

TL;DR

This study uses the GAEA semi-analytic model to analyze the relationship between Damped Lyman-alpha systems and atomic hydrogen in galaxies, comparing predictions with observations across different redshifts.

Contribution

It introduces a detailed simulation framework for DLAs using the GAEA model, improving agreement with observations by adjusting galaxy HI masses and sizes.

Findings

The model reproduces the shape of the column density distribution function at 2<z<3.

Increasing HI masses and disk radii improves the match with observed DLA statistics.

The typical DLA host halo mass is around 10^11 solar masses.

Abstract

Using the GAEA semi-analytic model, we analyse the connection between Damped Lyman-$\alpha$ systems (DLAs) and HI in galaxies. Our state-of-the-art semi-analytic model is tuned to reproduce the local galaxy HI mass function, and that also reproduces other important galaxy properties, including the galaxy mass - gas metallicity relation. To produce catalogs of simulated DLAs we throw $10^5$ random lines of sight in a composite simulated volume: dark matter haloes with log$(\frac{M_{200}}{ M_{\odot}}) \geq 11.5$ are extracted from the Millennium Simulation, while for $9.2 \leq \log(\frac{M_{200}}{ M_{\odot}})<11.5$ we use the Millennium II, and for $8 \leq \log(\frac{M_{200}}{M_{\odot}}) < 9.2$ a halo occupation distribution model. At $2 < z < 3$, where observational data are more accurate, our fiducial model predicts the correct shape of the column density distribution function, but its normalization falls short of the observations, with the discrepancy increasing at higher redshift. The agreement with observations is significantly improved increasing both the HI masses and the disk radii of model galaxies by a factor 2, as implemented 'a posteriori' in our $2M-2R$ model. In the redshift range of interest, haloes with $M_{200} \geq {10}^{11} M_{\odot}$ give the major contribution to $\Omega_{\rm DLA}$, and the typical DLA host halo mass is $\sim {10}^{11} M _{\odot}$. The simulated DLA metallicity distribution is in relatively good agreement with observations, but our model predicts an excess of DLAs at low metallicities. Our results suggest possible improvements for the adopted modelling of the filtering mass and metal ejection in low-mass haloes.