We present a comprehensive analysis of atomic hydrogen (HI) properties using a semi-analytical model of galaxy formation and N-body simulations covering a large cosmological volume at high resolution. We examine the HI mass function and the HI density, characterizing both their redshift evolution and their dependence on hosting halo mass. We analyze the HI content of dark matter haloes in the local Universe and up to redshift $z=5$, discussing the contribution of different galaxy properties. We find that different assembly history plays a crucial role in the scatter of this relation. We propose new fitting functions useful for constructing mock HI maps with HOD techniques. We investigate the HI clustering properties relevant for future $21$~cm Intensity Mapping (IM) experiments, including the HI bias and the shot noise level. The HI bias increases with redshift and it is roughly flat on the largest scales probed. The scale dependency is found at progressively larger scales with increasing redshift, apart from a dip feature at $z=0$. The shot-noise values are consistent with the ones inferred by independent studies, confirming that shot-noise will not be a limiting factor for IM experiments. We detail the contribution from various galaxy properties on the HI power spectrum and their relation to the halo bias. We find that HI poor satellite galaxies play an important role at the scales of the 1-halo term. Finally, we present the $21$~cm signal in redshift space, a fundamental prediction to be tested against data from future radio telescopes such as SKA.

We extend our previous work focused at $z\sim0$, studying the redshift evolution of galaxy dynamical properties using the state-of-the-art semi-analytic model GAEA: we show that the predicted size-mass relation for disky/star forming and quiescent galaxies is in good agreement with observational estimates, up to $z\sim2$. Bulge dominated galaxies have sizes that are offset low with respect to observational estimates, mainly due to our implementation of disk instability at high redshift. At large masses, both quiescent and bulge dominated galaxies have sizes smaller than observed. We interpret this as a consequence of our most massive galaxies having larger gas masses than observed, and therefore being more affected by dissipation. We argue that a proper treatment of quasar driven winds is needed to alleviate this problem. Our model compact galaxies have number densities in agreement with observational estimates and they form most of their stars in small and low angular momentum high-$z$ halos. GAEA predicts that a significant fraction of compact galaxies forming at high-$z$ is bound to merge with larger structures at lower redshifts: therefore they are not the progenitors of normal-size passive galaxies at $z=0$. Our model also predicts a stellar-halo size relation that is in good agreement with observational estimates. The ratio between stellar size and halo size is proportional to the halo spin and does not depend on stellar mass but for the most massive galaxies, where AGN feedback leads to a significant decrease of the retention factor (from about 80 per cent to 20 per cent).

We use a state-of-the-art semi-analytic model to study the size and the specific angular momentum of galaxies. Our model includes a specific treatment for the angular momentum exchange between different galactic components. Disk scale radii are estimated from the angular momentum of the gaseous/stellar disk, while bulge sizes are estimated assuming energy conservation. The predicted size--mass and angular momentum--mass relations are in fair agreement with observational measurements in the local Universe, provided a treatment for gas dissipation during major mergers is included. Our treatment for disk instability leads to unrealistically small radii of bulges formed through this channel, and predicts an offset between the size--mass relations of central and satellite early-type galaxies, that is not observed. The model reproduces the observed dependence of the size--mass relation on morphology, and predicts a strong correlation between specific angular momentum and cold gas content. This correlation is a natural consequence of galaxy evolution: gas-rich galaxies reside in smaller halos, and form stars gradually until present day, while gas-poor ones reside in massive halos, that formed most of their stars at early epochs, when the angular momentum of their parent halos is low. The dynamical and structural properties of galaxies can be strongly affected by a different treatment for stellar feedback, as this would modify their star formation history. A higher angular momentum for gas accreted through rapid mode does not affect significantly the properties of massive galaxies today, but has a more important effect on low-mass galaxies at higher redshift.

We update our recently published model for GAlaxy Evolution and Assembly (GAEA), to include a self-consistent treatment of the partition of cold gas in atomic and molecular hydrogen. Our model provides significant improvements with respect to previous ones used for similar studies. In particular, GAEA (i) includes a sophisticated chemical enrichment scheme accounting for non-instantaneous recycling of gas, metals, and energy; (ii) reproduces the measured evolution of the galaxy stellar mass function; (iii) reasonably reproduces the observed correlation between galaxy stellar mass and gas metallicity at different redshifts. These are important prerequisites for models considering a metallicity dependent efficiency of molecular gas formation. We also update our model for disk sizes and show that model predictions are in nice agreement with observational estimates for the gas, stellar and star forming disks at different cosmic epochs. We analyse the influence of different star formation laws including empirical relations based on the hydro-static pressure of the disk, analytic models, and prescriptions derived from detailed hydro-dynamical simulations. We find that modifying the star formation law does not affect significantly the global properties of model galaxies, neither their distributions. The only quantity showing significant deviations in different models is the cosmic molecular-to-atomic hydrogen ratio, particularly at high redshift. Unfortunately, however, this quantity also depends strongly on the modelling adopted for additional physical processes. Useful constraints on the physical processes regulating star formation can be obtained focusing on low mass galaxies and/or at higher redshift. In this case, self-regulation has not yet washed out differences imprinted at early time.

In this work, we study the basic statistical properties of HI-selected galaxies extracted from six different semi-analytic models, all run on the same cosmological N-body simulation. One model includes an explicit treatment for the partition of cold gas into atomic and molecular hydrogen. All models considered agree nicely with the measured HI mass function in the local Universe, with the measured scaling relations between HI and galaxy stellar mass, and with the predicted 2-point correlation function for HI rich galaxies. One exception is given by one model that predicts very little HI associated with galaxies in haloes above 10^12 Msun: we argue this is due to a too efficient radio-mode feedback for central galaxies, and to a combination of efficient stellar feedback and instantaneous stripping of hot gas for satellites. We demonstrate that treatment of satellite galaxies introduces large uncertainties at low HI masses. While models assuming non instantaneous stripping of hot gas tend to form satellite galaxies with HI masses slightly smaller than those of centrals with the same stellar mass, instantaneous gas stripping does not translate necessarily in lower HI masses. In fact, the adopted stellar feedback and star formation affect the satellites too. We analyze the relation between HI content and spin of simulated haloes: low spin haloes tend to host HI poor galaxies while high spin haloes are populated by galaxies in a wide range of HI mass. In our simulations, this is due to a correlation between the initial gas disk size and the halo spin.