The column densities of molecular gas across cosmic time: bridging observations and simulations

Observations of the cosmic evolution of different gas phases across time indicate a marked increase in the molecular gas mass density towards z similar to 2-3. Such a transformation implies an accompanied change in the global distribution of molecular hydrogen column densities (N-H(2)). Using observations by PHANGS-ALMA/SDSS and simulations by GRIFFIN/IllustrisTNG we explore the evolution of this H2 column density distribution function [f(N-H(2))]. The H-2 (and H I) column density maps for TNG50 and TNG100 are derived in post-processing and are made available through the IllustrisING online API. The shape and normalization of f(N-H2) of individual main-sequence star-forming galaxies are correlated with the star formation rate (SFR), stellar mass (M-*), and H-2 mass (M-H2) in both observations and simulations. TNG100, combined with 112 post-processing models, broadly reproduces observations, albeit with differences in slope and normalization. Also, an analytically modelled f(N), based on exponential gas discs, matches well with the simulations. The GRIFFIN simulation gives first indications that the slope of f(N-H(2)) might not majorly differ when including non-equilibrium chemistry in simulations. The f(N-H(2)) by TNG100 implies that higher molecular gas column densities are reached at z = 3 than at z = 0. Further, denser regions contribute more to the molecular mass density at z = 3. Finally, H-2 starts dominating compared to H I only at column densities above log(N-H(2)/cm(-2)) similar to 21.8-22 at both redshifts. These results imply that neutral atomic gas is an important contributor to the overall cold gas mass found in the ISM of galaxies including at densities typical for molecular clouds at z = 0 and 3.