Quantum dynamical investigation of the isotope effect in H2 formation on graphite at cold collision energies

The Eley–Rideal abstraction of hydrogen atoms on graphitic surfaces at cold collision energies was investigated using a time-dependent wave packet method within the rigid-flat surface approximation, with a focus on hydrogen–deuterium isotopic substitutions. It is found that the marked isotope effect of collinear collisions disappears when the full dimensionality of the problem is taken into account, thereby suggesting that abstraction is less direct than commonly believed and proceeds through glancing rather than head-on collisions. In contrast, a clear isotope effect is observed for “hot-atom” formation, which appears to be strongly favored for heavy projectiles because of their higher density of physisorbed states. Overall, the dynamics is essentially classical and reasonably well described by quasi-classical trajectory methods at all but the lowest energies (≲10 meV). A comparison of the results obtained in the (substrate) adiabatic and diabatic limits suggests that the reaction is only marginally affected by the lattice dynamics, but highlights the importance of including energy dissipation processes in order to accurately describe the internal excitation of the product molecules.