Generalized Temporal Acceleration Scheme for Kinetic Monte Carlo Simulations of Surface Catalytic Processes by Scaling the Rates of Fast Reactions

Comparative simulations of methanol steam reforming on PdZn alloy using kinetic Monte Carlo and mean-field microkinetic model

Methanol steam reforming (MSR) is an attractive route for producing clean energy hydrogen. PdZn alloys are extensively studied as potential MSR catalysts for their stability and high CO2 selectivity. Here, we investigated the reaction mechanism using density functional calculations, mean-field microkinetic modeling (MF-MKM), and kinetic Monte Carlo (kMC) simulations. To overcome the over-underestimation of CO2 selectivity by log-kMC, an ads-kMC algorithm is proposed in which the adsorption/desorption rate constants were reduced under certain requirements and the diffusion process was treated by redistributing surface species each time an event occured. The simulations show that the dominant pathway to CO2 at low temperatures is CH3OH → CH3O → CH2O → H2COOH → H2COO → HCOO → CO2. The ads-kMC predicted OH coverage is 2-3 times that of MF-MKM, while they produce similar coverage for other species. Analyses indicate that surface OH promotes the dehydrogenation of CH3OH, CH3O, and H2COOH significantly and plays a key role in the MSR process. The dissociation of water/methanol is the most important rate-limiting/rate-inhibiting step. The CO2 selectivity obtained by the two methods is close to each other and consistent with the experimental trend with temperature. Generally, the ads-kMC results agree with the MF-MKM ones, supporting the previous finding that kMC and MF-MKM predict similar results if the diffusion is very fast and adsorbate interactions are neglected. The present study sheds light on the MSR process on PdZn alloys, and the proposed scheme to overcome the stiff problems in kMC simulations is worthy of being extended to other systems.

On-the-fly kinetic Monte Carlo simulations with neural network potentials for surface diffusion and reaction

We develop an adaptive scheme in the kinetic Monte Carlo simulations, where the adsorption and activation energies of all elementary steps, including the effects of other adsorbates, are evaluated "on-the-fly" by employing the neural network potentials. The configurations and energies evaluated during the simulations are stored for reuse when the same configurations are sampled in a later step. The present scheme is applied to hydrogen adsorption and diffusion on the Pd(111) and Pt(111) surfaces and the CO oxidation reaction on the Pt(111) surface. The effects of interactions between adsorbates, i.e., adsorbate-adsorbate lateral interactions, are examined in detail by comparing the simulations without considering lateral interactions. This study demonstrates the importance of lateral interactions in surface diffusion and reactions and the potential of our scheme for applications in a wide variety of heterogeneous catalytic reactions.

Multiscale Investigation of the Mechanism and Selectivity of CO2 Hydrogenation over Rh(111)

CO2 hydrogenation over Rh catalysts comprises multiple reaction pathways, presenting a wide range of possible intermediates and end products, with selectivity toward either CO or methane being of particular interest. We investigate in detail the reaction mechanism of CO2 hydrogenation to the single-carbon (C1) products on the Rh(111) facet by performing periodic density functional theory (DFT) calculations and kinetic Monte Carlo (kMC) simulations, which account for the adsorbate interactions through a cluster expansion approach. We observe that Rh readily facilitates the dissociation of hydrogen, thus contributing to the subsequent hydrogenation processes. The reverse water-gas shift (RWGS) reaction occurs via three different reaction pathways, with CO hydrogenation to the COH intermediate being a key step for CO2 methanation. The effects of temperature, pressure, and the composition ratio of the gas reactant feed are considered. Temperature plays a pivotal role in determining the surface coverage and adsorbate composition, with competitive adsorption between CO and H species influencing the product distribution. The observed adlayer configurations indicate that the adsorbed CO species are separated by adsorbed H atoms, with a high ratio of H to CO coverage on the Rh(111) surface being essential to promote CO2 methanation.

Kinetic Trapping of Charge-Transfer Molecules at Metal Interfaces

Despite the common expectation that conjugated organic molecules on metals adsorb in a flat-lying layer, several recent studies have found coverage-dependent transitions to upright-standing phases, which exhibit notably different physical properties. In this work, we argue that from an energetic perspective, thermodynamically stable upright-standing phases may be more common than hitherto thought. However, for kinetic reasons, this phase may often not be observed experimentally. Using first-principles kinetic Monte Carlo simulations, we find that the structure with lower molecular density is (almost) always formed first, reminiscent of Ostwald's rule of stages. The phase transitions to the upright-standing phase are likely to be kinetically hindered under the conditions typically used in surface science. The simulation results are experimentally confirmed for the adsorption of tetracyanoethylene on Cu(111) using infrared and X-ray photoemission spectroscopy. Investigating both the role of the growth conditions and the energetics of the interface, we find that the time for the phase transition is determined mostly by the deposition rate and, thus, is mostly independent of the nature of the molecule.

Multiscale modelling of biomolecular corona formation on metallic surfaces

In the realm of food industry, the choice of non-consumable materials used plays a crucial role in ensuring consumer safety and product quality. Aluminum is widely used in food packaging and food processing applications, including dairy products. However, the interaction between aluminum and milk content requires further investigation to understand its implications. In this work, we present the results of multiscale modelling of the interaction between various surfaces, that is (100), (110), and (111), of fcc aluminum with the most abundant milk proteins and lactose. Our approach combines atomistic molecular dynamics, a coarse-grained model of protein adsorption, and kinetic Monte Carlo simulations to predict the protein corona composition in the deposited milk layer on aluminum surfaces. We consider a simplified model of milk, which is composed of the six most abundant milk proteins found in natural cow milk and lactose, which is the most abundant sugar found in dairy. Through our study, we ranked selected proteins and lactose adsorption affinities based on their corresponding interaction strength with aluminum surfaces and predicted the content of the naturally forming biomolecular corona. Our comprehensive investigation sheds light on the implications of aluminum in food processing and packaging, particularly concerning its interaction with the most abundant milk proteins and lactose. By employing a multiscale modelling approach, we simulated the interaction between metallic aluminum surfaces and the proteins and lactose, considering different crystallographic orientations. The results of our study provide valuable insights into the mechanisms of lactose and protein deposition on aluminum surfaces, which can aid in the general understanding of protein corona formation.

Lithium-ion battery multi-scale modeling coupled with simplified electrochemical model and kinetic Monte Carlo model

The multi-scale modeling of lithium-ion battery (LIB) is difficult and necessary due to its complexity. However, it is difficult to capture the aging behavior of batteries, and the coupling mechanism between multiple scales is still incomplete. In this paper, a simplified electrochemical model (SEM) and a kinetic Monte Carlo (KMC)-based solid electrolyte interphase (SEI) film growth model are used to study the multi-scale characteristics of LIBs. The single-particle SEM (SP-SEM) is described for macro scale, and a simple and self-consistent multi-particle SEM (MP-SEM) is developed. Then, the KMC-based SEI model is established for micro-scale molecular evolution. And, the two models are coupled to construct the full-cycle multi-scale model. After modeling, validation is performed by using a commercial 18650-type LIB. Finally, the effect of parameters on the SEI model is studied, including qualitative trend analysis and quantitative sensitivity analysis. The growth of SEI film with different particle sizes is studied by MP-SEM coupling simulation.

Exact distributed kinetic Monte Carlo simulations for on-lattice chemical kinetics: lessons learnt from medium- and large-scale benchmarks

Kinetic Monte Carlo (KMC) simulations have been instrumental in multiscale catalysis studies, enabling the elucidation of the complex dynamics of heterogeneous catalysts and the prediction of macroscopic performance metrics, such as activity and selectivity. However, the accessible length- and time-scales have been a limiting factor in such simulations. For instance, handling lattices containing millions of sites with 'traditional' sequential KMC implementations is prohibitive owing to large memory requirements and long simulation times. We have recently established an approach for exact, distributed, lattice-based simulations of catalytic kinetics which couples the Time-Warp algorithm with the Graph-Theoretical KMC framework, enabling the handling of complex adsorbate lateral interactions and reaction events within large lattices. In this work, we develop a lattice-based variant of the Brusselator system, a prototype chemical oscillator pioneered by Prigogine and Lefever in the late 60s, to benchmark and demonstrate our approach. This system can form spiral wave patterns, which would be computationally intractable with sequential KMC, while our distributed KMC approach can simulate such patterns 15 and 36 times faster with 625 and 1600 processors, respectively. The medium- and large-scale benchmarks thus conducted, demonstrate the robustness of the approach, and reveal computational bottlenecks that could be targeted in further development efforts. This article is part of a discussion meeting issue 'Supercomputing simulations of advanced materials'.

Utilizing Data-Driven Optimization to Automate the Parametrization of Kinetic Monte Carlo Models

Kinetic Monte Carlo (kMC) simulations are a popular tool to investigate the dynamic behavior of stochastic systems. However, one major limitation is their relatively high computational costs. In the last three decades, significant effort has been put into developing methodologies to make kMC more efficient, resulting in an enhanced runtime efficiency. Nevertheless, kMC models remain computationally expensive. This is in particular an issue in complex systems with several unknown input parameters where often most of the simulation time is required for finding a suitable parametrization. A potential route for automating the parametrization of kinetic Monte Carlo models arises from coupling kMC with a data-driven approach. In this work, we equip kinetic Monte Carlo simulations with a feedback loop consisting of Gaussian Processes (GPs) and Bayesian optimization (BO) to enable a systematic and data-efficient input parametrization. We utilize the results from fast-converging kMC simulations to construct a database for training a cheap-to-evaluate surrogate model based on Gaussian processes. Combining the surrogate model with a system-specific acquisition function enables us to apply Bayesian optimization for the guided prediction of suitable input parameters. Thus, the amount of trial simulation runs can be considerably reduced facilitating an efficient utilization of arbitrary kMC models. We showcase the effectiveness of our methodology for a physical process of growing industrial relevance: the space-charge layer formation in solid-state electrolytes as it occurs in all-solid-state batteries. Our data-driven approach requires only 1-2 iterations to reconstruct the input parameters from different baseline simulations within the training data set. Moreover, we show that the methodology is even capable of accurately extrapolating into regions outside the training data set which are computationally expensive for direct kMC simulation. Concluding, we demonstrate the high accuracy of the underlying surrogate model via a full parameter space investigation eventually making the original kMC simulation obsolete.

Concepts Relevant for the Kinetic Analysis of Reversible Reaction Systems

The net rate of a reversible chemical reaction is the difference between unidirectional rates of traversal along forward and reverse reaction paths. In a multistep reaction sequence, the forward and reverse trajectories, in general, are not the microscopic reverse of one another; rather, each unidirectional route is comprised of distinct rate-controlling steps, intermediates, and transition states. Consequently, traditional descriptors of rate (e.g., reaction orders) do not reflect intrinsic kinetic information but instead conflate unidirectional contributions determined by (i) the microscopic occurrence of forward/reverse reactions (i.e., unidirectional kinetics) and (ii) the reversibility of reaction (i.e., nonequilibrium thermodynamics). This review aims to provide a comprehensive resource of analytical and conceptual tools which deconvolute the contributions of reaction kinetics and thermodynamics to disambiguate unidirectional reaction trajectories and precisely identify rate- and reversibility-controlling molecular species and steps in reversible reaction systems. The extrication of mechanistic and kinetic information from bidirectional reactions is accomplished through equation-based formalisms (e.g., De Donder relations) grounded in principles of thermodynamics and interpreted in the context of theories of chemical kinetics developed in the past 25 years. The aggregate of mathematical formalisms detailed herein is general to thermochemical and electrochemical reactions and encapsulates a diverse body of scientific literature encompassing chemical physics, thermodynamics, chemical kinetics, catalysis, and kinetic modeling.

Large-scale benchmarks of the time-warp/graph-theoretical kinetic Monte Carlo approach for distributed on-lattice simulations of catalytic kinetics

Motivated by the need to perform large-scale kinetic Monte Carlo (KMC) simulations, in the context of unravelling complex phenomena such as catalyst reconstruction and pattern formation, we extend the work of Ravipati et al. [S. Ravipati, G. D. Savva, I.-A. Christidi, R. Guichard, J. Nielsen, R. Réocreux and M. Stamatakis, Comput. Phys. Commun., 2022, 270, 108148] in benchmarking the performance of a distributed-computing, on-lattice KMC approach. The latter, implemented in our software package Zacros, combines the graph-theoretical KMC framework with the Time-Warp algorithm for parallel discrete event simulations, and entails dividing the lattice into subdomains, each assigned to a processor. The cornerstone of the Time-Warp algorithm is the state queue, to which snapshots of the simulation state are saved regularly, enabling historical KMC information to be corrected when conflicts occur at subdomain boundaries. Focusing on three model systems, we highlight the key Time-Warp parameters that can be tuned to optimise performance. The frequency of state saving, controlled by the state saving interval, δ_snap, is shown to have the largest effect on performance, which favours balancing the overhead of re-simulating KMC history with that of writing state snapshots to memory. Also important is the global virtual time (GVT) computation interval, Δτ_GVT, which has little direct effect on the progress of the simulation but controls how often the state queue memory can be freed up. We also find that pre-allocating memory for the state queue data structure favours performance. These findings will guide users in maximising the efficiency of Zacros or other distributed KMC software, which is a vital step towards realising accurate, meso-scale simulations of heterogeneous catalysis.

Ensemble Effects in Adsorbate-Adsorbate Interactions in Microkinetic Modeling

Adsorbates on a surface experience lateral interactions that result in a distribution of adsorption energies. The adsorbate-adsorbate interactions are known to affect the kinetics of surface reactions, which motivates efforts to develop models that accurately account for the interactions. Here, we use density functional theory (DFT) calculations combined with Monte Carlo simulations to investigate how the distribution of adsorbates affects adsorption and desorption of CO from Pt(111). We find that the mean of the average adsorption energy determines the adsorption process, whereas the desorption process can be described by the low energy part of the adsorbate stability distribution. The simulated results are in very good agreement with calorimetry and temperature-programmed desorption experiments and provide a guideline of how to include adsorbate-adsorbate interactions in DFT-based mean-field kinetic models.

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.

Effects of surface diffusion in electrocatalytic CO2 reduction on Cu revealed by kinetic Monte Carlo simulations

Kinetic Monte Carlo (KMC) methods are frequently used for mechanistic studies of thermally driven heterogeneous catalysis systems but are underused for electrocatalysis. Here, we develop a lattice KMC approach for electrocatalytic CO₂ reduction. The work is motivated by a prior experimental report that performed electroreduction of a mixed feed of ¹²CO₂ and ¹³CO on Cu; differences in the ¹³C content of C2 products ethylene and ethanol (Δ¹³C) were interpreted as evidence of site selectivity. The lattice KMC model considers the effect of surface diffusion on this system. In the limit of infinitely fast diffusion (mean-field approximation), the key intermediates ¹²CO* and ¹³CO* would be well mixed on the surface and no evidence of site selectivity could have been observed. Using a simple two-site model and adapting a previously reported microkinetic model, we assess the effects of diffusion on the relative isotope fractions in the products using the estimated surface diffusion rate of CO* from literature reports. We find that the size of the active sites and the total surface adsorbate coverage can have a large influence on the values of Δ¹³C that can be observed. Δ¹³C is less sensitive to the CO* diffusion rate as long as it is within the estimated range. We further offer possible methods to estimate surface distribution of intermediates and to predict intrinsic selectivity of active sites based on experimental observations. This work illustrates the importance of considering surface diffusion in the study of electrochemical CO₂ reduction to multi-carbon products. Our approach is entirely based on a freely available open-source code, so will be readily adaptable to other electrocatalytic systems.

Real-time dynamics and structures of supported subnanometer catalysts via multiscale simulations

Understanding the performance of subnanometer catalysts and how catalyst treatment and exposure to spectroscopic probe molecules change the structure requires accurate structure determination under working conditions. Experiments lack simultaneous temporal and spatial resolution and could alter the structure, and similar challenges hinder first-principles calculations from answering these questions. Here, we introduce a multiscale modeling framework to follow the evolution of subnanometer clusters at experimentally relevant time scales. We demonstrate its feasibility on Pd adsorbed on CeO2(111) at various catalyst loadings, temperatures, and exposures to CO. We show that sintering occurs in seconds even at room temperature and is mainly driven by free energy reduction. It leads to a kinetically (far from equilibrium) frozen ensemble of quasi-two-dimensional structures that CO chemisorption and infrared experiments probe. CO adsorption makes structures flatter and smaller. High temperatures drive very rapid sintering toward larger, stable/metastable equilibrium structures, where CO induces secondary structure changes only.

KIMERA: A Kinetic Montecarlo Code for Mineral Dissolution

KIMERA is a scientific tool for the study of mineral dissolution. It implements a reversible Kinetic Monte Carlo (KMC) method to study the time evolution of a dissolving system, obtaining the dissolution rate and information about the atomic scale dissolution mechanisms. KIMERA allows to define the dissolution process in multiple ways, using a wide diversity of event types to mimic the dissolution reactions, and define the mineral structure in great detail, including topographic defects, dislocations, and point defects. Therefore, KIMERA ensures to perform numerous studies with great versatility. In addition, it offers a good performance thanks to its parallelization and efficient algorithms within the KMC method. In this manuscript, we present the code features and show some examples of its capabilities. KIMERA is controllable via user commands, it is written in object-oriented C++, and it is distributed as open-source software.

The enhanced extended phenomenological kinetics method to deal with timescale disparity problem among different reaction pathways

Kinetic Monte Carlo method can provide valuable mechanistic insights for catalytic systems. Nonetheless, it suffers from the notorious problem of timescale disparity due to the existence of the complex catalytic network that consists of fast events and slow events. Previously, we have proposed the extended phenomenological kinetics (XPK) method that effectively deals with the timescale disparity problem between diffusion and reaction. However, it remains a great challenge to simulate systems with timescale disparity among different reaction pathways, which is important when selectivity is the major concern. In this study, we implement the enhanced XPK method to address this problem. The new algorithm works by identifying states connected through fast transitions and compressing them into a "superstate" when the chosen states satisfy a local steadystate condition. This state compression algorithm simplifies the reaction network by concealing the fast transitions. The accuracy and efficiency of the algorithm are demonstrated by two model systems: selective catalytic hydrogenation and selective catalytic decomposition. The enhanced XPK method is expected to be beneficial to the kinetic simulations of catalytic systems, especially those with complex reaction networks.

A Caching Scheme To Accelerate Kinetic Monte Carlo Simulations of Catalytic Reactions

Kinetic Monte Carlo (KMC) simulations have been instrumental in advancing our fundamental understanding of heterogeneously catalyzed reactions, with particular emphasis on structure sensitivity, ensemble effects, and the interplay between adlayer structure and adsorbate-adsorbate lateral interactions in shaping the observed kinetics. Yet, the computational cost of KMC remains high, thereby motivating the development of acceleration schemes that would improve the simulation efficiency. We present an exact such scheme, which implements a caching algorithm along with shared-memory parallelization to improve the computational performance of simulations incorporating long-range adsorbate-adsorbate lateral interactions. This scheme is based on caching information about the energetic interaction patterns associated with the products of each possible lattice process (adsorption, desorption, reaction etc.). Thus, every time a reaction occurs ("ongoing reaction"), it enables fast updates of the rate constants of "affected reactions", i.e., possible reactions in the region of influence of the "ongoing reaction". Benchmarks on KMC simulations of NO_x oxidation/reduction, yielded acceleration factors of up to 20, when comparing single-thread runs without caching to runs on 16 threads with caching, for simulations with a cluster expansion Hamiltonian that incorporates up to 8th-nearest-neighbor interactions.

Acceleration scheme for particle transport in kinetic Monte Carlo methods

Kinetic Monte Carlo (kMC) simulations are frequently used to study (electro-)chemical processes within science and engineering. kMC methods provide insight into the interplay of stochastic processes and can link atomistic material properties with macroscopic characteristics. Significant problems concerning the computational demand arise if processes with large time disparities are competing. Acceleration algorithms are required to make slow processes accessible. Especially, the accelerated superbasin kMC (AS-kMC) scheme has been frequently applied within chemical reaction networks. For larger systems, the computational overhead of the AS-kMC is significant as the computation of the superbasins is done during runtime and comes with the need for large databases. Here, we propose a novel acceleration scheme for diffusion and transport processes within kMC simulations. Critical superbasins are detected during the system initialization. Scaling factors for the critical rates within the superbasins, as well as a lower bound for the number of sightings, are derived. Our algorithm exceeds the AS-kMC in the required simulation time, which we demonstrate with a 1D-chain example. In addition, we apply the acceleration scheme to study the time-of-flight (TOF) of charge carriers within organic semiconductors. In this material class, time disparities arise due to a significant spread of transition rates. The acceleration scheme allows a significant acceleration up to a factor of 65 while keeping the error of the TOF values negligible. The computational overhead is negligible, as all superbasins only need to be computed once.

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

The first-principles kinetic Monte Carlo (kMC) simulation has been demonstrated as a reliable multiscale modeling approach in silico to disclose the interplay among all the elementary steps in a complex reaction network for heterogeneous catalysis. Heterogeneous catalytic systems frequently contain fast surface diffusion processes of some adsorbates while the elementary steps in it would be much slower than those in fast diffusion. Consequently, the kMC simulation for these systems is easily trapped in the sub-basins of a super basin on a potential energy surface due to the continuous and repeated sampling of these fast processes, which would significantly increase the total accessible simulation time and even make it impossible to get the reasonable simulation results using the kMC simulation. In this work, we present an improved fast species redistribution (FSR) method for the kMC simulation to overcome the stiffness problem resulting from the low-barrier surface diffusion to accelerate the heterogeneous catalytic kMC simulation. Taking CO oxidations on Pt(111) and Pt(100) as examples, we demonstrate that the FSR approach can properly reproduce the results of an equivalent first-principles microkinetic model simulation with more reasonable reaction rates. The improved kMC simulation based on the FSR method can accurately incorporate the effect of the fast diffusion of species on the surface and provide several orders of magnitude of acceleration compared to the standard kMC simulation.

Kinetic Monte Carlo Modeling for the NO-CO Reaction Mechanism on Rh(100) and Rh(111)

The NO-CO reaction on Rh(100) and Rh(111) is a prototypical catalytic system with various practical applications, including the treatment of automotive gas exhausts. With parameters derived from first-principles calculations, the Brønsted-Evans-Polanyi (BEP) relation for the reaction steps of NO-CO on Rh(100) and Rh(111) surfaces is fitted, which is more accurate and practical for the calculation of the effect of interaction between adsorbates on activation energy compared to the basic BEP relation. Further, a kinetic Monte Carlo (kMC) model for the NO-CO reaction systems on Rh(100) and Rh(111) is constructed for the exploration of the system's reaction mechanism. Besides the temperature and pressure, the coverage and activation sites are essential factors for reaction kinetic of the NO-CO reaction system. Our results are beneficial for designing more efficient, economical, and environmentally friendly next-generation catalysts.

Dynamic vs static behaviour of a supported nanoparticle with reaction-induced catalytic sites in a lattice model

Modern literature shows a rapidly growing interest to the supported nanocatalysts with dynamic behaviour under reaction conditions. This new frontier of heterogeneous catalysis is recognized as one of the most challenging and worthy of consideration from all possible angles. In this context, a previously suggested lattice model is used to get an insight, by means of kinetic Monte Carlo, into the influence of the mobility of reaction-induced catalytic sites of a two-dimensional supported nanoparticle on the system behaviour. The results speak in favour of feasibility of dynamic nanocatalysts with self-organized structures capable of robust functioning. This approach, from the macroscopic end, is believed to be a useful complement to ever developing experimental and first principle approaches.

A DFT and KMC based study on the mechanism of the water gas shift reaction on the Pd(100) surface

We present a combined density functional theory (DFT) and Kinetic Monte Carlo (KMC) study of the water gas shift (WGS) reaction on the Pd(100) surface. We propose a mechanism comprising both the redox and the associative pathways for the WGS within a single framework, which consists of seven core elementary steps, which in turn involve splitting of a water molecule followed by the production of an H-atom and an OH-species on the Pd(100) surface. In the following steps, these intermediates then recombine with each other and with CO leading to the evolution of CO₂, and H₂. Seven other elementary steps, involving the diffusion and adsorption of the surface intermediate species are also considered for a complete description of the mechanism. The geometrical and electronic properties of each of the reactants, products, and the transition states of the core elementary steps are presented. We also discuss the analysis of Bader charges and spin densities for the reactants, transition states and the products of these elementary steps. Our study indicates that the WGS reaction progresses simultaneously via the direct oxidation and the carboxyl paths on the Pd(100) surface.

A Practical Guide to Surface Kinetic Monte Carlo Simulations

This review article is intended as a practical guide for newcomers to the field of kinetic Monte Carlo (KMC) simulations, and specifically to lattice KMC simulations as prevalently used for surface and interface applications. We will provide worked out examples using the kmos code, where we highlight the central approximations made in implementing a KMC model as well as possible pitfalls. This includes the mapping of the problem onto a lattice and the derivation of rate constant expressions for various elementary processes. Example KMC models will be presented within the application areas surface diffusion, crystal growth and heterogeneous catalysis, covering both transient and steady-state kinetics as well as the preparation of various initial states of the system. We highlight the sensitivity of KMC models to the elementary processes included, as well as to possible errors in the rate constants. For catalysis models in particular, a recurrent challenge is the occurrence of processes at very different timescales, e.g., fast diffusion processes and slow chemical reactions. We demonstrate how to overcome this timescale disparity problem using recently developed acceleration algorithms. Finally, we will discuss how to account for lateral interactions between the species adsorbed to the lattice, which can play an important role in all application areas covered here.

SQERT-T: alleviating kinetic Monte Carlo (KMC)-stiffness in transient KMC simulations

Lattice based kinetic Monte Carlo (KMC) is often used for simulating the dynamics of systems at a supramolecular scale, based on molecular scale transitions. A common challenge in KMC simulations is rapid 'back-and-forth' reactions, which dominate the events executed during simulations and inhibit the ability for simulations to reach longer time scales. Such processes are fast frivolous processes (FFPs) and are one manifestation of a phenomenon referred to as KMC-stiffness. Here, an algorithm for staggered quasi-equilibrium rank-based throttling geared towards transient kinetics (SQERT-T) is presented. Within the SQERT-T methodology, a pace-restrictor reaction and an FFP floor are utilized along with throttling of the process transition rate constants to accelerate the KMC simulations while still retaining sufficient time resolution for sampling of the data. KMC simulations were performed for CO oxidation over RuO₂(1 1 0) and over RuO₂(1 1 1), and the results were compared to experimental data obtained using RuO₂ powders. The experiments and simulations were for transient conditions: the system was subjected to a temperature program which included temperatures in the range of 363 to 453 K. The timescales that were achieved during the KMC simulations in this study would not have been accessible without KMC acceleration, and were enabled by the use of SQERT-T.

Assessment of mean-field microkinetic models for CO methanation on stepped metal surfaces using accelerated kinetic Monte Carlo

First-principles screening studies aimed at predicting the catalytic activity of transition metal (TM) catalysts have traditionally been based on mean-field (MF) microkinetic models, which neglect the effect of spatial correlations in the adsorbate layer. Here we critically assess the accuracy of such models for the specific case of CO methanation over stepped metals by comparing to spatially resolved kinetic Monte Carlo (kMC) simulations. We find that the typical low diffusion barriers offered by metal surfaces can be significantly increased at step sites, which results in persisting correlations in the adsorbate layer. As a consequence, MF models may overestimate the catalytic activity of TM catalysts by several orders of magnitude. The potential higher accuracy of kMC models comes at a higher computational cost, which can be especially challenging for surface reactions on metals due to a large disparity in the time scales of different processes. In order to overcome this issue, we implement and test a recently developed algorithm for achieving temporal acceleration of kMC simulations. While the algorithm overall performs quite well, we identify some challenging cases which may lead to a breakdown of acceleration algorithms and discuss possible directions for future algorithm development.

A Lattice Kinetic Monte Carlo Solver for First-Principles Microkinetic Trend Studies

Mean-field microkinetic models in combination with Brønsted-Evans-Polanyi like scaling relations have proven highly successful in identifying catalyst materials with good or promising reactivity and selectivity. Analysis of the microkinetic model by means of lattice kinetic Monte Carlo promises a faithful description of a range of atomistic features involving short-range ordering of species in the vicinity of an active site. In this paper, we use the "fruit fly" example reaction of CO oxidation on fcc(111) transition and coinage metals to motivate and develop a lattice kinetic Monte Carlo solver suitable for the numerically challenging case of vastly disparate rate constants. As a result, we show that for the case of infinitely fast diffusion and absence of adsorbate-adsorbate interaction it is, in fact, possible to match the prediction of the mean-field-theory method and the lattice kinetic Monte Carlo method. As a corollary, we conclude that lattice kinetic Monte Carlo simulations of surface chemical reactions are most likely to provide additional insight over mean-field simulations if diffusion limitations or adsorbate-adsorbate interactions have a significant influence on the mixing of the adsorbates.