Multi-lattice Kinetic Monte Carlo Simulations from First Principles: Reduction of the Pd(100) Surface Oxide by CO

Aggregation of Size-Selected Oxide Clusters Deposited onto Au(111)

Kinetic Monte Carlo (kMC) simulations along with density functional theory (DFT) calculations were used to investigate the aggregation of size-selected Nb3Oy (y = 5, 6, 7) clusters deposited onto the Au(111) surface. Recent STM experiments showed that the cluster binding sites and sizes of the cluster assemblies on the Nb3Oy/Au(111) surfaces strongly depend on the stoichiometry of the clusters, i.e., the oxygen-to-niobium ratio. To better understand the origins of these differences, kMC simulations of the nucleation and growth of cluster assemblies were performed using energy barriers for diffusion and intercluster interactions estimated from DFT calculations of cluster binding and dimerization energies, respectively. Comparisons of the kMC simulations with STM images of the as-deposited Nb3Oy/Au(111) surfaces at RT and after high temperature annealing were used to further optimize the energetics and gauge the importance of nearest neighbor interactions. The kMC simulations demonstrate that the assembly of Nb3Oy clusters on Au(111) are largely controlled by the magnitude of the barriers for diffusion and interparticle-bond formation, while changes at higher temperatures are sensitive to the binding energies between nearest neighbors. Simulations for the Nb3O5 and Nb3O6 clusters, which exhibit smaller cluster assembly sizes in STM, required larger diffusion barriers as well as different barriers for interparticle binding, which reflected differences in DFT calculated dimerization energies. The results demonstrate the effectiveness of combined DFT and kMC calculations for understanding how the stoichiometry affects the aggregation of small oxide clusters on a metal surface.

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.

Computational Methods in Heterogeneous Catalysis

The unprecedented ability of computations to probe atomic-level details of catalytic systems holds immense promise for the fundamentals-based bottom-up design of novel heterogeneous catalysts, which are at the heart of the chemical and energy sectors of industry. Here, we critically analyze recent advances in computational heterogeneous catalysis. First, we will survey the progress in electronic structure methods and atomistic catalyst models employed, which have enabled the catalysis community to build increasingly intricate, realistic, and accurate models of the active sites of supported transition-metal catalysts. We then review developments in microkinetic modeling, specifically mean-field microkinetic models and kinetic Monte Carlo simulations, which bridge the gap between nanoscale computational insights and macroscale experimental kinetics data with increasing fidelity. We finally review the advancements in theoretical methods for accelerating catalyst design and discovery. Throughout the review, we provide ample examples of applications, discuss remaining challenges, and provide our outlook for the near future.

Vacancy assisted diffusion on single-atom surface alloys

Bimetallic surfaces can exhibit an improved catalytic activity through tailoring the concentration and/or the arrangement of the two metallic components. However, in order to be catalytically active, the active bimetallic surface structure has to be stable under operating conditions. Typically, structural changes in metals occur via vacancy diffusion. Based on the first-principles determination of formation energies and diffusion barriers we have performed kinetic Monte-Carlo (kMC) simulations to analyse the (meta-)stability of PtRu/Ru(0001), AgPd/Pd(111), PtAu/Au(111) and InCu/Cu(100) surface alloys. In a first step, here we consider single-atom alloys together with one vacancy per simulation cell. We will present results of the time evolution of these structures and analyse them in terms of the interaction between the constituents of the bimetallic surface.

Modeling nonadiabatic dynamics in condensed matter materials: some recent advances and applications

This review focuses on recent developments in the field of nonadiabatic molecular dynamics (NA-MD), with particular attention given to condensed-matter systems. NA-MD simulations for small molecular systems can be performed using high-level electronic structure (ES) calculations, methods accounting for the quantization of nuclear motion, and using fewer approximations in the dynamical methodology itself. Modeling condensed-matter systems imposes many limitations on various aspects of NA-MD computations, requiring approximations at various levels of theory-from the ES, to the ways in which the coupling of electrons and nuclei are accounted for. Nonetheless, the approximate treatment of NA-MD in condensed-phase materials has gained a spin lately in many applied studies. A number of advancements of the methodology and computational tools have been undertaken, including general-purpose methods, as well as those tailored to nanoscale and condensed matter systems. This review summarizes such methodological and software developments, puts them into the broader context of existing approaches, and highlights some of the challenges that remain to be solved.

Operando chemistry of catalyst surfaces during catalysis

The chemistry of a catalyst surface during catalysis is crucial for a fundamental understanding of the mechanisms of a catalytic reaction performed on the catalyst in the gas or liquid phase.