Importance of equivariant features in machine-learning interatomic potentials for reactive chemistry at metal surfaces

Abstract

Reactive chemistry of molecular hydrogen at surfaces, notably dissociative sticking and hydrogen evolution, play a crucial role in energy storage, fuel cells, and chemical synthesis. Copper is a particularly interesting metal for studying these processes due to its widespread use as both a catalyst in industry and a model catalyst in fundamental research. Theoretical studies can help to decipher underlying mechanisms and reaction design, but studying these systems computationally is challenging due to the complex electronic structure of metal surfaces and the high sensitivity towards reaction barriers. In addition, ab initio molecular dynamics, based on density functional theory, is too computationally demanding to explicitly simulate reactive sticking or desorption probabilities. A promising solution to such problems can be provided through high-dimensional machine learning-based interatomic potentials (MLIPs). Despite the remarkable accuracy and fidelity of MLIPs, particularly in molecular and bulk inorganic materials simulations, their application to different facets of hybrid systems and the selection of appropriate representations remain largely unexplored. This paper addresses these issues and investigates how feature equivariance in MLIPs impacts adaptive sampling workflows and data efficiency. Specifically, we develop high-dimensional MLIPs to investigate reactive hydrogen scattering on copper surfaces and compare the performance of various MLIPs that use equivariant features for atomic representation (PaiNN) with those that use invariant representations (SchNet). Our findings demonstrate that using equivariant features can greatly enhance the accuracy and reliability of MLIPs for gas surface dynamics and that this approach should become the standard in this field