Estimating Free Energy Barriers for Heterogeneous Catalytic Reactions with Machine Learning Potentials and Umbrella Integration

Abstract

Predicting the rate constants of elementary reaction steps is key for the computational modeling of catalytic processes. Within transition state theory (TST), this requires an accurate estimation of the corresponding free energy barriers. While sophisticated methods for estimating free energy differences exist, these typically require extensive (biased) molecular dynamics simulations that are computationally prohibitive with the first-principles electronic structure methods that are typically used in catalysis research. In this contribution, we show that machine-learning (ML) interatomic potentials can be trained in an automated iterative workflow to perform such free energy calculations at a much reduced computational cost as compared to a direct density functional theory (DFT) based evaluation. For the decomposition of CHO on Rh(111), we find that thermal effects are substantial and lead to a decrease in the free energy barrier, which can be vanishingly small, depending on the DFT functional used. This is in stark contrast to previously reported estimates based on a harmonic TST approximation, which predicted an increase in the barrier at elevated temperatures. Since CHO is the reactant of the putative rate limiting reaction step in syngas conversion on Rh(111) and essential for the selectivity toward oxygenates containing multiple carbon atoms (C2+ oxygenates), our results call into question the reported mechanism established by microkinetic models.

Toward an ab Initio Description of Adsorbate Surface Dynamics

“…Although the calculation of reaction barriers remains cumbersome, recent work by us and others has demonstrated how machine learning potentials (MLPs) can be used to overcome this bottleneck. − For instance, using Cu-exchanged zeolites as a prototypical example, we have explicitly calculated the transition state geometries and reaction barriers of methane activation for thousands of [CuOCu] 2+ sites across 52 zeolites. While most sites show linear trends between the C–H activation barrier and the H binding energy, our analysis identifies several important factors (e.g., confinement and accessibility) that cause deviation from the expected universal scaling behavior .…”

Transferable Machine Learning Interatomic Potential for Bond Dissociation Energy Prediction of Drug-like Molecules

“…In recent years, machine learning interatomic potentials (MLIPs) have emerged as successful tools for accurately approximating ab initio potential energy surfaces at a significantly reduced computational cost: roughly milliseconds instead of minutes per single-point evaluation of a typical structure found in this work. − MLIPs have found applications in various computational chemistry problems, ranging from highly accurate spectra prediction long MD simulations of electrolytes, , path-integral MD of supramolecular complexes, to crystal structure prediction, heterogeneous catalysis, , and radical reaction networks . MLIPs offer an efficiency advantage by avoiding the need to solve the all-electron Schrödinger equation; instead, they predict energy and forces from the atom positions directly.…”