Accurately modeling gas-surface collision dynamics presents a great challenge for theory, especially in the low energy (or temperature) regime where quantum effects are important. Here, a path integral based non-equilibrium ring polymer molecular dynamics (NE-RPMD) approach is adapted to calculate dissociative initial sticking probabilities (S0) of H2 on Cu(111) and D2O on Ni(111), revealing distinct quantum nature in the two benchmark surface reactions. NE-RPMD successfully captures quantum tunneling in H2 dissociation at very low energies, where the quasi-classical trajectory (QCT) method suddenly fails. Additionally, QCT substantially overestimates S0 of D2O due to severe zero point energy (ZPE) leakage, even at collision energies higher than the ZPE-corrected barrier. Instead, NE-RPMD predicts S0 values of D2O in much improved agreement with reference results obtained by the quantum wavepacket method with reasonable corrections of the thermal contribution. Our results suggest NE-RPMD as a promising approach to model quantum effects in gas-surface reactions.
TOC graphic3 Dissociative adsorption of small molecules at surfaces is the initial and often rate-limiting step in many interfacial processes such as heterogeneous catalysis and corrosion. Initial sticking probability (S0) is an important observable to reveal the adsorption mechanism and dynamics at solid surfaces, which can be now precisely measured as a function of incidence energy by molecular beam experiments 1-2 . S0 can be 10 -5 or even lower at low energies in many highly activated systems, e.g. methane and water dissociative chemisorption on metal surfaces 1, 3 , representing indispensable steps in methane steaming reforming. Given the large number of degrees of freedom (DOFs), however, it is very challenging to accurately predict S0 from first-principles calculations.Ideally, given an accurate global potential energy surface (PES), exact quantum initial sticking probabilities can be extracted by fully-coupled quantum dynamical (QD) methods, including both time-dependent 4-6 and time-independent 7 implementations. While such molecule-surface PESs can be now routinely developed 8-9 , high-dimensional QD calculations remain extremely difficult and are limited to involve at most nine molecular DOFs so far 10 , due to their poor scaling with dimensionality. Alternatively, the quasiclassical trajectory (QCT) method is much more efficient to model surface reaction dynamics and visualize the associated mechanism 4 . In QCT calculations, the zero point energy (ZPE) of the reactant is approximately included, while the atomic motion is evolved classically. Such QCT applications in H2 activated dissociation on various metal surfaces do reproduce the QD calculated S0 values above the ZPE-corrected barrier quite well 4, 11 .When explicit PESs are unavailable, ab initial molecular dynamics (AIMD) simulations, in which the energies and forces are calculated on-the-fly, have also been performed to study surface reactions 12 . Although AIMD is computationally expen...