First Principles Reactive Flux Theory for Surface Reactions: Multiple Channels and Recrossing Dynamics

Heterogenous reactions typically consist of multiple elementary steps and their rate coefficients are of fundamental importance in elucidating the mechanisms and micro-kinetics of these processes. Transition-state theory (TST) for calculating surface reaction rate coefficients often relies solely on the harmonic approximation of adsorbent vibrations and neglects recrossing dynamics. Here, we combine, for the first time, an efficient metadynamics enhanced sampling method with a more general reactive flux approach to calculate rate coefficients of surface reactions of any order and/or with multiple reaction coordinates, overcoming these limitations of TST. We apply this approach to a textbook surface reaction, CO oxidation on Pt(111), for which rate constants have been precisely measured, using a full-dimensional neural network potential energy surface constructed from first-principles data. An accurate multi-dimensional free-energy surface is obtained by incorporating three collective variables, yielding rate coefficients for both CO oxidation and the competing CO desorption that are in good agreement with experimental data. Interestingly, our results reveal significant dynamic recrossing in both channels, which however arises from distinct physical mechanisms. This approach represents an accurate and general framework for calculating rate coefficients of elementary surface processes from first-principles, which is vital for developing predictive kinetic models for heterogenous catalysis.