Hydrogen formation is a key process for the physics and the chemistry of interstellar clouds. Molecular hydrogen is believed to form on the carbonaceous surface of dust grains, and several mechanisms have been invoked to explain its abundance in different regions of space, from cold interstellar clouds to warm photondominated regions. Here, we investigate direct (Eley-Rideal) recombination including lattice dynamics, surface corrugation, and competing H-dimers formation by means of ab initio molecular dynamics. We find that Eley-Rideal reaction dominates at energies relevant for the interstellar medium and alone may explain observations if the possibility of facile sticking at special sites (edges, point defects, etc.) on the surface of the dust grains is taken into account.graphite-graphene | interstellar chemistry | hydrogen recombination | density functional theory T he complicated interaction between hydrogen atoms and graphite has been widely studied in the last decade because of its relevance in various fields, from hydrogen storage (1) to nuclear fusion (2), from graphene technology (3-6) to interstellar chemistry (7).In the chemistry of the interstellar medium (ISM), H 2 is a precursor for more complex molecules through its protonated form H þ 3 , shields the clouds from the radiation field, and acts as primary cooling agent during the gravitational collapse, which ultimately leads to star formation (7). Its appearance in the early universe steered the development of the first stars and triggered galaxy formation (8).In the harsh ISM conditions, molecular hydrogen is continuously dissociated by stellar UV radiation and cosmic rays; hence efficient synthetic pathways are needed to explain the observed abundances (7). These are known to involve the carbonaceous surface of interstellar dust grains (9), but the detailed mechanisms behind H 2 formation on graphitic surfaces have yet to be fully uncovered.Physisorbed H atoms can only be invoked in cold molecular clouds, where the surface temperature is low enough to prevent desorption (T s ≤ 40 K) from the shallow physisorption well (E ∼ 4 kJ·mol −1 ) (10, 11). In this case, hydrogen molecules may follow from either an Eley-Rideal (ER)/hot-atom (HA) or a LangmuirHinshelwood (LH) pathway, because H atoms are very mobile even at extremely low temperatures (11). Stable (chemisorbed) species on the basal graphitic plane, however, require energetic projectile atoms to overcome an activation barrier to sticking (E ∼ 20 kJ·mol −1 ) (12, 13), which is due to the carbon sp 2 -sp 3 rehybridization needed by the bond formation process, although tunneling might help in this context (14). Chemically bound H atoms are immobile on the surface [they desorb rather than diffusing (15)]; hence LH pathway is ruled out in sufficiently warm environments where only chemisorbed species can be found. No matter how a first chemisorbed species is formed, secondary H atoms undergo facile sticking and preferentially form ortho and para (16, 17) dimers because of the aromatic nature of the...