Abstract:Context: Molecular hydrogen ($\rm{H_2}$) is crucial in galaxy formation and evolution, serving as the main fuel for star formation (SF). In metal-enriched environments, $\rm{H_2}$ primarily forms on interstellar dust grain surfaces. However, due to the complexities of modelling this process, SF in cosmological simulations often relies on empirical or theoretical frameworks validated only in the Local Universe to estimate the abundance of $\rm{H_2}$. Aims: This study aims to model the connection between star, dust, and $\rm{H_2}$ formation processes in cosmological simulations. Methods: We include $\rm{H_2}$ formation on dust grain surfaces and account for molecule destruction and radiation shielding into the SF and feedback model MUPPI. Results: The model reproduces key properties of observed galaxies for stellar, dust, and $\rm{H_2}$ components. The cosmic density of $\rm{H_2}$ ($\rho_{\rm{H2}}$) peaks around $z=1.5$, then decreases by half towards $z=0$, showing milder evolution than observed. The $\rm{H_2}$ mass function since $z=2$ also shows gentler evolution. Our model successfully recovers the integrated molecular Kennicutt-Schmidt (mKS) law between surface star formation rate ($\Sigma_{\rm SFR}$) and surface $\rm{H_2}$ density ($\Sigma_{\rm H2}$) at $z=0$, already evident at $z=2$ with a higher normalization. We find hints of a broken power law with a steeper slope at higher $\Sigma_{\rm H2}$, aligning with some observational findings. Additionally, the $\rm{H_2}$-to-dust mass ratio in galaxies shows a decreasing trend with gas metallicity and stellar mass. The $\rm{H_2}$-to-dust mass fraction for the global galaxy population is higher at higher redshifts. The analysis of the atomic-to-molecular transition on a particle-by-particle basis suggests that gas metallicity cannot reliably substitute the dust-to-gas ratio in models simulating dust-promoted $\rm{H_2}$.

Vibrational level population of H$_2$ and H$_2^+$ in the early Universe

Non-equilibrium H$_2$ formation in the early Universe: energy exchanges, rate coefficients and spectral distortions

Quantum mechanical study of the high-temperature $\mathrm{H}^+ + \mathrm{HD} \to \mathrm{D}^+ + \mathrm{H}_2$ reaction for the primordial universe chemistry

Atomic and molecular gas from the epoch of reionization down to redshift 2

Influence of a quasi-molecular mechanism of recombination on the formation of hydrogen in the early universe

Vibrational Distribution of Hydrogen Molecular Ions in High-Energy Ionization Processes

Black-body radiation induced photodissociation and population redistribution of weakly bound states in H$_2^+$

H2 Formation in Low Metallicity Galaxies

Uncertainties in H2 and HD Chemistry and Cooling and their Role in Early Structure Formation

Molecular hydrogen formation on dust grains in the high redshift universe

Excitation mechanisms in newly discovered H2-bearing Damped Lyman-alpha clouds: systems with low molecular fractions

Modeling Molecular Hydrogen in Low-metallicity Galaxies

Dynamic Equilibrium Sets of the Atomic Content of Galaxies Across Cosmic Time

Molecular hydrogen in damped Ly-alpha systems: clues to interstellar physics at high-redshift

Intertwined Formation of $\rm{H_2}$, Dust, and Stars in Cosmological Simulations

H2 distribution during formation of multiphase molecular clouds

Modeling Molecular Hydrogen and Star Formation in Cosmological Simulations

Using a modified version of the Tavis-Cummings-Hubbard model to simulate the formation of neutral hydrogen molecule

Radiative cooling functions for primordial molecules

Modeling molecular hydrogen in low metallicity galaxies

Isotopic and vibrational-level dependence of H$_2$ dissociation by electron impact