Ethylene Epoxidation on Ag(100), Ag(110), and Ag(111): A Joint Ab Initio and Kinetic Monte Carlo Study and Comparison with Experiments

Speeding up the Detection of Adsorbate Lateral Interactions in Graph-Theoretical Kinetic Monte Carlo Simulations

Kinetic Monte Carlo (KMC) has become an indispensable tool in heterogeneous catalyst discovery, but realistic simulations remain computationally demanding on account of the need to capture complex and long-range lateral interactions between adsorbates. The Zacros software package (https://zacros.org) adopts a graph-theoretical cluster expansion (CE) framework that allows such interactions to be computed with a high degree of generality and fidelity. This involves solving a series of subgraph isomorphism problems in order to identify relevant interaction patterns in the lattice. In an effort to reduce the computational burden, we have adapted two well-known subgraph isomorphism algorithms, namely, VF2 and RI, for use in KMC simulations and implemented them in Zacros. To benchmark their performance, we simulate a previously established model of catalytic NO oxidation, treating the O* lateral interactions with a series of progressively larger CEs. For CEs with long-range interactions, VF2 and RI are found to provide impressive speedups relative to simpler algorithms. RI performs best, giving speedups reaching more than 150× when combined with OpenMP parallelization. We also simulate a recently developed methane cracking model, showing that RI offers significant improvements in performance at high surface coverages.

Enhanced Selectivity of the Propylene Epoxidation Reaction on a Cu Monolayer Surface via Eley-Rideal Mechanism

The aerobic oxidation of propylene to selectively achieve propylene oxide (PO) is a challenging reaction in catalysis. Therefore, an active catalyst which shows enhanced PO selectivity is extremely desired. In the present investigation, an attempt has been made to explore the catalytic activity of a mono-atomically thin two-dimensional (2D) hexagonal (HX) Cu layer for selective propylene epoxidation using molecular O₂ with the aid of density functional theory calculations. The results reveal that the conversion of propylene to PO via Eley-Rideal mechanism is an exoergic and barrierless reaction on the O₂ pre-adsorbed Cu monolayer. The Pauli energy component plays a decisive role for barrierless reaction whereas the electrostatic and orbital contribution governs the energetic stability of PO. Car-Parrinello molecular dynamics (CPMD) simulation reinforces the outcomes of climbing image nudged elastic band (CI-NEB) calculation. Further, the formation of oxametallacycle OMC-2 (0.47 eV) is kinetically favourable over OMC-1 (0.87 eV) and AHS (0.50 eV) on O pre-adsorbed 2D HX Cu. Interestingly, the energy barrier for the conversion of OMC-2 to PO (0.70 eV) is considerably low in comparison with the acetone formation (0.90 eV). Therefore, it is worth to mention that the 2D HX Cu surface provides a promising platform for selective propylene epoxidation.

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.

Directing transition metal-based oxygen-functionalization catalysis

This review presents the recent progress of oxygen functionalization reactions based on non-electrochemical (conventional organic synthesis) and electrochemical methods. Although both methods have their advantages and limitations, the former approach has been used to synthesize a broader range of organic substances as the latter is limited by several factors, such as poor selectivity and high energy cost. However, because electrochemical methods can replace harmful terminal oxidizers with external voltage, organic electrosynthesis has emerged as greener and more eco-friendly compared to conventional organic synthesis. The progress of electrochemical methods toward oxygen functionalization is presented by an in-depth discussion of different types of electrically driven-chemical organic synthesis, with particular attention to recently developed electrochemical systems and catalyst designs. We hope to direct the attention of readers to the latest breakthroughs of traditional oxygen functionalization reactions and to the potential of electrochemistry for the transformation of organic substrates to useful end products.

Ag (111) surface for ambient electrolysis of nitrogen to ammonia

In this paper, the reaction process of N₂ convert to NH₃ catalyzed by Ag (111) surface was obtained through the construction of Ag (111) surface and computational simulation. The charge transfer in the reaction process and the change of N≡N bond length are described. Since the N₂ reduction reaction (NRR) usually occurs under alkaline solution conditions, we calculated and described the coexistence of OH* and N₂. At the same time, the co-adsorption structure of OH* and N₂ at different adsorption sites was studied.

Revitalizing silver nanocrystals as a redox catalyst by modifying their surface with an isocyanide-based compound

Silver is an excellent catalyst for oxidation reactions such as ethylene epoxidation, but it shows limited activity toward reduction reactions. Here we report a strategy to revitalize Ag nanocrystals as a redox catalyst for the production of an aromatic azo compound by modifying their surface with an isocyanide-based compound. We also leverage in situ fingerprint spectroscopy to acquire molecular insights into the reaction mechanism by probing the vibrational modes of all chemical species at the catalytic surface with surface-enhanced Raman spectroscopy. We establish that binding of isocyanide to Ag nanocrystals makes it possible for Ag to extract the oxygen atoms from the nitro-groups of nitroaromatics and then use these atoms to oxidize isocyanide to isocyanate. Concurrently, the coupling between two adjacent deoxygenated nitroaromatic molecules leads to the formation of an aromatic azo compound.

Identification of an AgS2 Complex on Ag(110)

Adsorbed sulfur has been investigated on the Ag(110) surface at two different coverages, 0.02 and 0.25 monolayers. At the lower coverage, only sulfur adatoms are present. At the higher coverage, there are additional bright features which we identify as linear, independent AgS₂ complexes. This identification is based upon density functional theory (DFT) and its comparison with experimental observations including bias dependence and separation between complexes. DFT also predicts the absence of AgS₂ complexes at low coverage, and the development of AgS₂ complexes around a coverage of 0.25 monolayers of sulfur, as is experimentally observed. To our knowledge, this is the first example of an isolated linear sulfur-metal-sulfur complex.

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.