# Overcoming the challenges of simulating surface reactions with complex interaction models

** Published:**

The complexity of first-principles-based Kinetic Monte Carlo (KMC) models to treat chemical reactions at surfaces has risen dramatically over the past few years. Lateral interactions have been disregarded in kinetic models of surface reactions for a long time, in part due to the inability to obtain reliable parameters from experiments, in part due to the belief that these interactions play only a minor role in the kinetics of surface reactions, and also, in part, due to the extensive computational effort associated with computing them on the level of electronic structure theory and simulating reactions with lateral interactions in KMC.

Lateral interactions in KMC simulations are often modeled in a Cluster Expansion (CE) Hamiltonian, where the total energy of the lattice and adsorption energies are represented as a sum of interaction clusters, where each cluster represents a pattern in which adsorbates are arranged on the surface, for instance, pairwise, three-body, four body, and larger clusters. The Cluster Expansion is, in principle, exact, but is always truncated in practical application due to constraints in the number of parameters that can be determined (and meaningfully fitted) on a first principles level, and in the evaluation of the Hamiltonian.

While lateral interaction models for surface coadsorbate systems have steadily gained popularity and grown in terms of complexity, their use in chemical kinetics has been impeded by the low performance of current KMC algorithms. The origins of this additional computational cost are traced back to the more elaborate Cluster Expansion Hamiltonian, the more extensive rate updating, and to the impracticality of rate-catalog-based algorithms for interacting adsorbate systems. However, these challenges can be overcome by dedicated optimization strategies and we have made significant progress in developing new KMC algorithms using Cluster Expansion Hamiltonians. The code and benchmark have recently been released on Github and in a J. Comp. Chem. paper.

### Related publications

Probing the Activity of Different Oxygen Species in the CO Oxidation over RuO

_{2}(110) by Combining Transient Reflection–Absorption Infrared Spectroscopy with Kinetic Monte Carlo Simulations F. Hess, C. Sack, D. Langsdorf, H. Over.*ACS Catal.*7 (**2017**) 8420-8428.10.1021/acscatal.7b02838 Rate-Determining Step or Rate-Determining Configuration? The Deacon Reaction over RuO

_{2}(110) Studied by DFT-Based KMC Simulations F. Hess, H. Over.*ACS Catal.*7 (**2017**) 128-138.10.1021/acscatal.6b02575 Combined experiment and theory approach in surface chemistry: Stairway to heaven? K.S. Exner, F. Hess, H. Over, A.P. Seitsonen.

*Surf. Sci.*640 (**2015**) 165-180.10.1016/j.susc.2015.01.006 Perspective article: Kinetic Monte Carlo simulations of heterogeneously catalyzed oxidation reactions F. Hess, H. Over.

*Catal. Sci. Technol.*4 (**2014**) 583-598.10.1039/C3CY00833A One-dimensional confinement in heterogeneous catalysis: Trapped oxygen on RuO

_{2}(110) model catalysts F. Hess, P.P.T. Krause, S.F. Rohrlack, J.P. Hofmann, H. Over.*Surf. Sci.*606 (**2012**) L69-L73.10.1016/j.susc.2012.04.019 “First-Principles” Kinetic Monte Carlo simulations revisited: CO oxidation over RuO

_{2}(110) F. Hess, A. Farkas, A.P. Seitsonen, H. Over.*J. Comp. Chem.*33 (**2012**) 757-766.10.1002/jcc.22902 Experiment-Based Kinetic Monte Carlo Simulations: CO Oxidation over RuO

_{2}(110) A. Farkas, F. Hess, H. Over.*J. Phys. Chem. C.*116 (**2012**) 581-591.10.1021/jp204703p