Abstract:Simulations of chemical reaction probabilities in gas surface dynamics require the calculation of ensemble averages over many tens of thousands of reaction events to predict dynamical observables that can be compared to experiments. At the same time, the energy landscapes need to be accurately mapped, as small errors in barriers can lead to large deviations in reaction probabilities. This brings a particularly interesting challenge for machine learning interatomic potentials, which are becoming well-established tools to accelerate molecular dynamics simulations. We compare state-of-the-art machine learning interatomic potentials with a particular focus on their inference performance on CPUs and suitability for high throughput simulation of reactive chemistry at surfaces. The considered models include polarizable atom interaction neural networks (PaiNN), recursively embedded atom neural networks (REANN), the MACE equivariant graph neural network, and atomic cluster expansion potentials (ACE). The models are applied to a dataset on reactive molecular hydrogen scattering on low-index surface facets of copper. All models are assessed for their accuracy, time-to-solution, and ability to simulate reactive sticking probabilities as a function of the rovibrational initial state and kinetic incidence energy of the molecule. REANN and MACE models provide the best balance between accuracy and time-to-solution and can be considered the current state-of-the-art in gas-surface dynamics. PaiNN models require many features for the best accuracy, which causes significant losses in computational efficiency. ACE models provide the fastest time-to-solution, however, models trained on the existing dataset were not able to achieve sufficiently accurate predictions in all cases.

Ridge-based Bias Potentials to Accelerate Molecular Dynamics.

Hyper Molecular Dynamics with A Local Bias Potential

Local bias potential in hyper molecular dynamics method

Automated construction of effective potential via algorithmic implicit bias

Constrained Broyden Dimer Method with Bias Potential for Exploring Potential Energy Surface of Multistep Reaction Process.

Enhanced Sampling of Multidimensional Free-Energy Landscapes Using Adaptive Biasing Forces

Application of the Accelerated Molecular Dynamics Simulations to the Folding of A Small Protein

Accelerated Sampling of Rare Events using a Neural Network Bias Potential

Transition rates, survival probabilities, and quality of bias from time-dependent biased simulations

Free energy calculations: An efficient adaptive biasing potential method

Good rates from bad coordinates: the exponential average time-dependent rate approach

Refining Potential Energy Surface through Dynamical Properties via Differentiable Molecular Simulation

Adaptively accelerating reactive molecular dynamics using boxed molecular dynamics in energy space

Advantages of multi-dimensional biasing in accelerated dynamics: application to the calculation of the acid $pK_a$ for acetic acid

Reinforced Dynamics for Enhanced Sampling in Large Atomic and Molecular Systems

Benchmarking of machine learning interatomic potentials for reactive hydrogen dynamics at metal surfaces

Neural Network and Nearest Neighbor Algorithms for Enhancing Sampling of Molecular Dynamics

Universal Aspects of Barrier Crossing Under Bias

Uncertainty-biased molecular dynamics for learning uniformly accurate interatomic potentials

Reweighted Jarzynski sampling: Acceleration of rare events and free energy calculation with a bias potential learned from nonequilibrium work