A fast species redistribution approach to accelerate the kinetic Monte Carlo simulation for heterogeneous catalysis

Achieving Theory-Experiment Parity for Activity and Selectivity in Heterogeneous Catalysis Using Microkinetic Modeling

Microkinetic modeling based on density functional theory (DFT) energies plays an essential role in heterogeneous catalysis because it reveals the fundamental chemistry for catalytic reactions and bridges the microscopic understanding from theoretical calculations and experimental observations. Microkinetic modeling requires building a set of ordinary differential equations (ODEs) based on the calculation results of thermodynamic properties of adsorbates and kinetic parameters for the reaction elementary steps. Solving a microkinetic model can extract information on catalytic chemistry, including critical reaction intermediates, reaction pathways, the surface species distribution, activity, and selectivity, thus providing vital guidelines for altering catalysts.

However, the quantitative reliability of traditional microkinetic models is often insufficient to conclusively extrapolate the mechanistic details of complex reaction systems. This can be attributed to several factors, the most important of which is the limitation of obtaining an accurate estimation of the energy inputs via traditional calculation methods. These limitations include the difficulty of using static DFT methods to calculate reaction energies of adsorption/desorption processes, often rate-controlling or selectivity-determining steps, and the inadequate consideration of surface coverage effects. In addition, the robust microkinetic software is rare, which also complicates the resolution of complex catalytic systems.

In this Account, we review our recent works toward refining the predictions of microkinetic modeling in heterogeneous catalysis and achieving theory–experiment parity for activity and selectivity. First, we introduce CATKINAS, a microkinetic software developed in our group, and show how it disentangles the problem that traditional microkinetic software has and how it can now be applied to obtain kinetic results for more sophisticated reaction systems. Second, we describe a molecular dynamics method developed recently to obtain the free-energy changes for the adsorption/desorption process to fill in the missing energy inputs. Third, we show that a rigorous consideration of surface coverage effects is pivotal for building more realistic models and obtaining accurate kinetic results. Following a series of studies on acetylene hydrogenation reactions on Pd catalysts, we demonstrate how this new approach can provide an improved quantitative understanding of the mechanism, active site, and intrinsic structural sensitivity. Finally, we conclude with a brief outlook and the remaining challenges in this field.

Kinetic Monte Carlo simulations for heterogeneous catalysis: Fundamentals, current status, and challenges

Kinetic Monte Carlo (KMC) simulations in combination with first-principles (1p)-based calculations are rapidly becoming the gold-standard computational framework for bridging the gap between the wide range of length scales and time scales over which heterogeneous catalysis unfolds. 1p-KMC simulations provide accurate insights into reactions over surfaces, a vital step toward the rational design of novel catalysts. In this Perspective, we briefly outline basic principles, computational challenges, successful applications, as well as future directions and opportunities of this promising and ever more popular kinetic modeling approach.

Perspective on computational reaction prediction using machine learning methods in heterogeneous catalysis

Heterogeneous catalysis plays a significant role in the modern chemical industry. Towards the rational design of novel catalysts, understanding reactions over surfaces is the most essential aspect. Typical industrial catalytic processes such as syngas conversion and methane utilisation can generate a large reaction network comprising thousands of intermediates and reaction pairs. This complexity not only arises from the permutation of transformations between species but also from the extra reaction channels offered by distinct surface sites. Despite the success in investigating surface reactions at the atomic scale, the huge computational expense of ab initio methods hinders the exploration of such complicated reaction networks. With the proliferation of catalysis studies, machine learning as an emerging tool can take advantage of the accumulated reaction data to emulate the output of ab initio methods towards swift reaction prediction. Here, we briefly summarise the conventional workflow of reaction prediction, including reaction network generation, ab initio thermodynamics and microkinetic modelling. An overview of the frequently used regression models in machine learning is presented. As a promising alternative to full ab initio calculations, machine learning interatomic potentials are highlighted. Furthermore, we survey applications assisted by these methods for accelerating reaction prediction, exploring reaction networks, and computational catalyst design. Finally, we envisage future directions in computationally investigating reactions and implementing machine learning algorithms in heterogeneous catalysis.

Perspective on theoretical methods and modeling relating to electro-catalysis processes

Theoretical methods and models for the description of thermodynamics and kinetics in electro-catalysis, including solvent effects, externally applied potentials, and many-body interactions, are discussed.

An approach to calculate the free energy changes of surface reactions using free energy decomposition on ab initio brute-force molecular dynamics trajectories

To understand the mechanisms and kinetics of catalytic reactions in heterogeneous catalysis, ab initio molecular dynamics is one of the powerful methods used to explore the free energy surface (FES) of surface elementary steps.