AVS 66th International Symposium & Exhibition
    Fundamental Discoveries in Heterogeneous Catalysis Focus Topic Monday Sessions
       Session HC+SS-MoM

Paper HC+SS-MoM6
First-Principles Kinetic Monte Carlo Simulation of CO Oxidation on PdO(101): Role of Oxygen Vacancies

Monday, October 21, 2019, 10:00 am, Room A212

Session: Utilization of Theoretical Models, Machine Learning, and Artificial Intelligence for Heterogeneously-Catalyzed Reactions
Presenter: Minkyu Kim, The Ohio State University
Authors: M. Kim, The Ohio State University
A. Asthagiri, The Ohio State University
Correspondent: Click to Email

CO oxidation on transition metal (TM) oxide surfaces has been widely studied both experimentally and theoretically; however, a healthy debate continues on the coupling between changes in oxide phase and surface reactivity. In this study, we investigated CO oxidation over PdO(101) surface, which has been proposed to be important in Pd oxidation catalysis. In contrast to earlier studies of CO oxidation on TM oxide surfaces, we incorporate neighbor effects of oxygen vacancies on all the elementary surface processes. We employ density functional theory (DFT) to map out the kinetics of 50 elementary surface processes. We find that barriers for elementary steps such as O vacancy, O2 dissociation, and CO2 formation can be decreased by 20-40% in the presence of O vacancies, while O2 desorption is increased by 70%. Using the DFT-derived rate constants, we have developed a lattice-based kinetic Monte Carlo (kMC) framework that can simulate CO oxidation under both ultra-high vacuum (UHV) and reaction conditions.

Initially, the kMC simulations were performed under UHV conditions (low CO partial pressure, no O2pressure) as a function of increasing surface temperature. As the CO oxidation proceeds and the surface lattice O atom coverage is depleted, the CO oxidation rate decreases rapidly at 400 K; however, the rate sharply increases at temperatures above 450 K. At 450 K, we find a new complex pathway to CO2 formation that is activated in the presence of O vacancies and is the source of the increased CO oxidation rate despite the depletion of surface oxygen atoms. These kMC results match UHV isothermal experiments under similar conditions. KMC simulations at steady state conditions of low CO and O2pressures (PCO: 5x10-9 Torr / PO2: 1.5x10-8Torr), show that the rate-limiting step is O2 dissociation and this elementary step requires the presence of adjacent oxygen vacancies to be activated at temperatures below 500 K. Without the addition of O vacancy neighbor effects to the kMC model, the PdO(101) surface would be inactive to CO oxidation because surface oxygens cannot be healed by gaseous O2. In addition, we will discuss kMC simulations under reaction conditions at varying CO/O2partial pressures.