DOE PAGES title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: Electrons to Reactors Multiscale Modeling: Catalytic CO Oxidation over RuO2

Abstract

First-principles kinetic Monte Carlo (1p-kMC) simulations for CO oxidation on two RuO2 facets, RuO2(110) and RuO2(111), were coupled to the computational fluid dynamics (CFD) simulations package MFIX, and reactor-scale simulations were then performed. 1p-kMC coupled with CFD has recently been shown as a feasible method for translating molecular scale mechanistic knowledge to the reactor scale, enabling comparisons to in situ and online experimental measurements. Only a few studies with such coupling have been published. This work incorporates multiple catalytic surface facets into the scale-coupled simulation, and three possibilities were investigated: the two possibilities of each facet individually being the dominant phase in the reactor, and also the possibility that both facets were present on the catalyst particles in the ratio predicted by an ab initio thermodynamics-based Wulff construction. When lateral interactions between adsorbates were included in the 1p-kMC simulations, the two surfaces, RuO2(110) and RuO2(111), were found to be of similar order-of-magnitude in activity for the pressure range of 1 × 10–4 bar to 1 bar, with the RuO2(110) surface-termination showing more simulated activity than the RuO2(111) surface-termination. Coupling between the 1p-kMC and CFD was achieved with a lookup table generated by the error-based modified Shepard interpolation scheme. Isothermal reactormore » scale simulations were performed and compared to two separate experimental studies, conducted with reactant partial pressures of ≤0.1 bar. Simulations without an isothermality restriction were also conducted and showed that the simulated temperature gradient across the catalytic reactor bed is <0.5 K, which validated the use of the isothermality restriction for investigating the reactor-scale phenomenological temperature dependences. The approach with the Wulff construction based reactor simulations reproduced a trend similar to one experimental data set relatively well, with the (110) surface being more active at higher temperaures; in contrast, for the other experimental data set, our reactor simulations achieve surprisingly and perhaps fortuitously good agreement with the activity and phenomenological pressure dependence when it is assumed that the (111) facet is the only active facet present. Lastly, the active phase of catalytic CO oxidation over RuO2 remains unsettled, but the present study presents proof of principle (and progress) toward more accurate multiscale modeling from electrons to reactors and new simulation results.« less

Authors:
ORCiD logo [1]; ORCiD logo [2]; ORCiD logo [1]; ORCiD logo [3]; ORCiD logo [4];  [5]; ORCiD logo [1]
  1. Oak Ridge National Lab. (ORNL), Oak Ridge, TN (United States)
  2. Technische Univ. Munchen, Garching (Germany)
  3. Oak Ridge National Lab. (ORNL), Oak Ridge, TN (United States); Fiat Chrysler Automobiles, Auburn Hills, MI (United States)
  4. Oak Ridge National Lab. (ORNL), Oak Ridge, TN (United States); SABIC, Houston, TX (United States)
  5. Free Univ., Berlin (Germany)
Publication Date:
Research Org.:
Oak Ridge National Lab. (ORNL), Oak Ridge, TN (United States)
Sponsoring Org.:
USDOE Office of Science (SC)
OSTI Identifier:
1436956
Grant/Contract Number:  
AC05-00OR22725
Resource Type:
Accepted Manuscript
Journal Name:
ACS Catalysis
Additional Journal Information:
Journal Volume: 8; Journal ID: ISSN 2155-5435
Publisher:
American Chemical Society (ACS)
Country of Publication:
United States
Language:
English
Subject:
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CHEMISTRY; computational fluid dynamics; interpolation; kinetic Monte Carlo; kmos; mfix; multiscale; ruthenium; ruthenium oxide

Citation Formats

Sutton, Jonathan E., Lorenzi, Juan M., Krogel, Jaron T., Xiong, Qingang, Pannala, Sreekanth, Matera, Sebastian, and Savara, Aditya. Electrons to Reactors Multiscale Modeling: Catalytic CO Oxidation over RuO2. United States: N. p., 2018. Web. doi:10.1021/acscatal.8b00713.
Sutton, Jonathan E., Lorenzi, Juan M., Krogel, Jaron T., Xiong, Qingang, Pannala, Sreekanth, Matera, Sebastian, & Savara, Aditya. Electrons to Reactors Multiscale Modeling: Catalytic CO Oxidation over RuO2. United States. https://doi.org/10.1021/acscatal.8b00713
Sutton, Jonathan E., Lorenzi, Juan M., Krogel, Jaron T., Xiong, Qingang, Pannala, Sreekanth, Matera, Sebastian, and Savara, Aditya. Fri . "Electrons to Reactors Multiscale Modeling: Catalytic CO Oxidation over RuO2". United States. https://doi.org/10.1021/acscatal.8b00713. https://www.osti.gov/servlets/purl/1436956.
@article{osti_1436956,
title = {Electrons to Reactors Multiscale Modeling: Catalytic CO Oxidation over RuO2},
author = {Sutton, Jonathan E. and Lorenzi, Juan M. and Krogel, Jaron T. and Xiong, Qingang and Pannala, Sreekanth and Matera, Sebastian and Savara, Aditya},
abstractNote = {First-principles kinetic Monte Carlo (1p-kMC) simulations for CO oxidation on two RuO2 facets, RuO2(110) and RuO2(111), were coupled to the computational fluid dynamics (CFD) simulations package MFIX, and reactor-scale simulations were then performed. 1p-kMC coupled with CFD has recently been shown as a feasible method for translating molecular scale mechanistic knowledge to the reactor scale, enabling comparisons to in situ and online experimental measurements. Only a few studies with such coupling have been published. This work incorporates multiple catalytic surface facets into the scale-coupled simulation, and three possibilities were investigated: the two possibilities of each facet individually being the dominant phase in the reactor, and also the possibility that both facets were present on the catalyst particles in the ratio predicted by an ab initio thermodynamics-based Wulff construction. When lateral interactions between adsorbates were included in the 1p-kMC simulations, the two surfaces, RuO2(110) and RuO2(111), were found to be of similar order-of-magnitude in activity for the pressure range of 1 × 10–4 bar to 1 bar, with the RuO2(110) surface-termination showing more simulated activity than the RuO2(111) surface-termination. Coupling between the 1p-kMC and CFD was achieved with a lookup table generated by the error-based modified Shepard interpolation scheme. Isothermal reactor scale simulations were performed and compared to two separate experimental studies, conducted with reactant partial pressures of ≤0.1 bar. Simulations without an isothermality restriction were also conducted and showed that the simulated temperature gradient across the catalytic reactor bed is <0.5 K, which validated the use of the isothermality restriction for investigating the reactor-scale phenomenological temperature dependences. The approach with the Wulff construction based reactor simulations reproduced a trend similar to one experimental data set relatively well, with the (110) surface being more active at higher temperaures; in contrast, for the other experimental data set, our reactor simulations achieve surprisingly and perhaps fortuitously good agreement with the activity and phenomenological pressure dependence when it is assumed that the (111) facet is the only active facet present. Lastly, the active phase of catalytic CO oxidation over RuO2 remains unsettled, but the present study presents proof of principle (and progress) toward more accurate multiscale modeling from electrons to reactors and new simulation results.},
doi = {10.1021/acscatal.8b00713},
journal = {ACS Catalysis},
number = ,
volume = 8,
place = {United States},
year = {Fri Apr 20 00:00:00 EDT 2018},
month = {Fri Apr 20 00:00:00 EDT 2018}
}

Journal Article:
Free Publicly Available Full Text
Publisher's Version of Record

Citation Metrics:
Cited by: 29 works
Citation information provided by
Web of Science

Save / Share:

Works referenced in this record:

On the CO-Oxidation over Oxygenated Ruthenium
journal, January 2009

  • Rosenthal, Dirk; Girgsdies, Frank; Timpe, Olaf
  • Zeitschrift für Physikalische Chemie, Vol. 223, Issue 1-2
  • DOI: 10.1524/zpch.2009.6032

Kinetic Evidence for a Non-Langmuir-Hinshelwood Surface Reaction: H/D Exchange over Pd Nanoparticles and Pd(111)
journal, April 2013

  • Savara, Aditya; Ludwig, Wiebke; Schauermann, Swetlana
  • ChemPhysChem, Vol. 14, Issue 8
  • DOI: 10.1002/cphc.201300179

kmos: A lattice kinetic Monte Carlo framework
journal, July 2014

  • Hoffmann, Max J.; Matera, Sebastian; Reuter, Karsten
  • Computer Physics Communications, Vol. 185, Issue 7
  • DOI: 10.1016/j.cpc.2014.04.003

Error propagation in first-principles kinetic Monte Carlo simulation
journal, April 2017


Multiscale modeling for emergent behavior, complexity, and combinatorial explosion
journal, March 2012


CO oxidation over ruthenium: identification of the catalytically active phases at near-atmospheric pressures
journal, January 2012

  • Gao, Feng; Goodman, D. Wayne
  • Physical Chemistry Chemical Physics, Vol. 14, Issue 19
  • DOI: 10.1039/c2cp40121e

Scale-up and multiphase reaction engineering
journal, August 2015


First-Principles Analysis of Structure Sensitivity in NO Oxidation on Pt
journal, January 2015

  • Bray, J. M.; Schneider, W. F.
  • ACS Catalysis, Vol. 5, Issue 2
  • DOI: 10.1021/cs501783q

Accurate coverage-dependence incorporated into first-principles kinetic models: Catalytic NO oxidation on Pt (111)
journal, February 2012


Efficient interpolation of precomputed kinetic data employing reduced multivariate Hermite Splines
journal, March 2017


SQERTSS: Dynamic rank based throttling of transition probabilities in kinetic Monte Carlo simulations
journal, October 2017

  • Danielson, Thomas; Sutton, Jonathan E.; Hin, Céline
  • Computer Physics Communications, Vol. 219
  • DOI: 10.1016/j.cpc.2017.05.016

Building large microkinetic models with first-principles׳ accuracy at reduced computational cost
journal, January 2015


Understanding the Structural Deactivation of Ruthenium Catalysts on an Atomic Scale under both Oxidizing and Reducing Conditions
journal, December 2004

  • Aßmann, Jens; Crihan, Daniela; Knapp, Marcus
  • Angewandte Chemie International Edition, Vol. 44, Issue 6
  • DOI: 10.1002/anie.200461805

Transport limitations and bistability for in situ CO oxidation at RuO 2 ( 110 ) : First-principles based multiscale modeling
journal, August 2010


A bubbling fluidization model using kinetic theory of granular flow
journal, April 1990


Structure-Activity Correlations for the Oxidation of CO over Polycrystalline RuO 2 Powder Derived from Steady-State and Transient Kinetic Experiments
journal, July 2005

  • Narkhede, Vijay; Aßmann, Jens; Muhler, Martin
  • Zeitschrift für Physikalische Chemie, Vol. 219, Issue 7-2005
  • DOI: 10.1524/zpch.219.7.979.67092

Kinetic modelling of heterogeneous catalytic systems
journal, November 2014


Coupling CFD with detailed microkinetic modeling in heterogeneous catalysis
journal, June 2013


Coupling of kinetic Monte Carlo simulations of surface reactions to transport in a fluid for heterogeneous catalytic reactor modeling
journal, February 2013

  • Schaefer, C.; Jansen, A. P. J.
  • The Journal of Chemical Physics, Vol. 138, Issue 5
  • DOI: 10.1063/1.4789419

Hierarchical multi-scale model reduction in the simulation of catalytic converters
journal, April 2013


Advanced coal gasifier designs using large-scale simulations
journal, July 2009


Parallel kinetic Monte Carlo simulation framework incorporating accurate models of adsorbate lateral interactions
journal, December 2013

  • Nielsen, Jens; d’Avezac, Mayeul; Hetherington, James
  • The Journal of Chemical Physics, Vol. 139, Issue 22
  • DOI: 10.1063/1.4840395

Multivariate interpolation of large sets of scattered data
journal, June 1988

  • Renka, Robert J.
  • ACM Transactions on Mathematical Software, Vol. 14, Issue 2
  • DOI: 10.1145/45054.45055

Analyzing the Case for Bifunctional Catalysis
journal, March 2016

  • Andersen, Mie; Medford, Andrew J.; Nørskov, Jens K.
  • Angewandte Chemie International Edition, Vol. 55, Issue 17
  • DOI: 10.1002/anie.201601049

High throughput multiscale modeling for design of experiments, catalysts, and reactors: Application to hydrogen production from ammonia
journal, January 2010


Multiscale optimization using hybrid PDE/kMC process systems with application to thin film growth
journal, December 2005


Ab Initio Thermodynamics and First-Principles Microkinetics for Surface Catalysis
journal, January 2016


A multiregion operator-splitting CFD approach for coupling microkinetic modeling with internal porous transport in heterogeneous catalytic reactors
journal, January 2016

  • Maffei, Tiziano; Gentile, Giancarlo; Rebughini, Stefano
  • Chemical Engineering Journal, Vol. 283
  • DOI: 10.1016/j.cej.2015.08.080

Heterogeneous oxidation catalysis on ruthenium: bridging the pressure and materials gaps and beyond
journal, April 2008


Exploring Pretreatment-Morphology Relationships: Ab Initio Wulff Construction for RuO 2 Nanoparticles under Oxidising Conditions
journal, July 2013


Assessing the reliability of calculated catalytic ammonia synthesis rates
journal, July 2014


When atomic-scale resolution is not enough: Spatial effects on in situ model catalyst studies
journal, November 2012


Role of RuO 2 (100) in surface oxidation and CO oxidation catalysis on Ru(0001)
journal, January 2016

  • Flege, Jan Ingo; Lachnitt, Jan; Mazur, Daniel
  • Physical Chemistry Chemical Physics, Vol. 18, Issue 1
  • DOI: 10.1039/C5CP05807D

Effect of local metal microstructure on adsorption on bimetallic surfaces: Atomic nitrogen on Ni/Pt(111)
journal, May 2013

  • Guo, Wei; Vlachos, Dionisios G.
  • The Journal of Chemical Physics, Vol. 138, Issue 17
  • DOI: 10.1063/1.4803128

Addressing global uncertainty and sensitivity in first-principles based microkinetic models by an adaptive sparse grid approach
journal, January 2018

  • Döpking, Sandra; Plaisance, Craig P.; Strobusch, Daniel
  • The Journal of Chemical Physics, Vol. 148, Issue 3
  • DOI: 10.1063/1.5004770

Effects of correlated parameters and uncertainty in electronic-structure-based chemical kinetic modelling
journal, February 2016

  • Sutton, Jonathan E.; Guo, Wei; Katsoulakis, Markos A.
  • Nature Chemistry, Vol. 8, Issue 4
  • DOI: 10.1038/nchem.2454

Silicon Chemical Vapor Deposition on macro and submicron powders in a fluidized bed
journal, March 2009


Uncertainty Quantification Framework Applied to the Water–Gas Shift Reaction over Pt-Based Catalysts
journal, May 2016

  • Walker, Eric; Ammal, Salai Cheettu; Terejanu, Gabriel A.
  • The Journal of Physical Chemistry C, Vol. 120, Issue 19
  • DOI: 10.1021/acs.jpcc.6b01348

Wavelet-based surrogate time series for multiscale simulation of heterogeneous catalysis
journal, April 2016


Implications of coverage-dependent O adsorption for catalytic NO oxidation on the late transition metals
journal, January 2014

  • Frey, Kurt; Schmidt, David J.; Wolverton, C.
  • Catal. Sci. Technol., Vol. 4, Issue 12
  • DOI: 10.1039/C4CY00763H

Composition, structure, and stability of RuO 2 ( 110 ) as a function of oxygen pressure
journal, December 2001


Unraveling the Complexity of Catalytic Reactions via Kinetic Monte Carlo Simulation: Current Status and Frontiers
journal, November 2012

  • Stamatakis, Michail; Vlachos, Dionisios G.
  • ACS Catalysis, Vol. 2, Issue 12
  • DOI: 10.1021/cs3005709

Evidence for the Active Phase of Heterogeneous Catalysts through In Situ Reaction Product Imaging and Multiscale Modeling
journal, June 2015


Effect of errors in linear scaling relations and Brønsted–Evans–Polanyi relations on activity and selectivity maps
journal, June 2016


Kinetics of carbon monoxide oxidation on single-crystal palladium, platinum, and iridium
journal, September 1988

  • Berlowitz, Paul J.; Peden, Charles H. F.; Goodman, D. Wayne
  • The Journal of Physical Chemistry, Vol. 92, Issue 18
  • DOI: 10.1021/j100329a030

Structure sensitivity in oxide catalysis: First-principles kinetic Monte Carlo simulations for CO oxidation at RuO 2 (111)
journal, November 2015

  • Wang, Tongyu; Reuter, Karsten
  • The Journal of Chemical Physics, Vol. 143, Issue 20
  • DOI: 10.1063/1.4936354

Effect of multiscale model uncertainty on identification of optimal catalyst properties
journal, July 2011


In situ adaptive tabulation for the CFD simulation of heterogeneous reactors based on operator-splitting algorithm
journal, August 2016

  • Bracconi, Mauro; Maestri, Matteo; Cuoci, Alberto
  • AIChE Journal, Vol. 63, Issue 1
  • DOI: 10.1002/aic.15441

CO oxidation on ruthenium: The nature of the active catalytic surface
journal, October 2007


Reduced order modeling and dynamic optimization of multiscale PDE/kMC process systems
journal, September 2008


Kinetics of carbon monoxide oxidation over ruthenium(0001)
journal, March 1986

  • Peden, Charles H. F.; Goodman, D. Wayne
  • The Journal of Physical Chemistry, Vol. 90, Issue 7
  • DOI: 10.1021/j100398a031

Oscillatory Behavior in the CO-oxidation over Bulk Ruthenium Dioxide — the Effect of the CO/O 2 Ratio
journal, January 2011

  • Rosenthal, Dirk; Girgsdies, Frank; Timpe, Olaf
  • Zeitschrift für Physikalische Chemie, Vol. 225, Issue 1
  • DOI: 10.1524/zpch.2011.5515

Beyond mean-field approximations for accurate and computationally efficient models of on-lattice chemical kinetics
journal, July 2017

  • Pineda, M.; Stamatakis, M.
  • The Journal of Chemical Physics, Vol. 147, Issue 2
  • DOI: 10.1063/1.4991690

Does phenomenological kinetics provide an adequate description of heterogeneous catalytic reactions?
journal, May 2007

  • Temel, Burcin; Meskine, Hakim; Reuter, Karsten
  • The Journal of Chemical Physics, Vol. 126, Issue 20
  • DOI: 10.1063/1.2741556

Universal transition state scaling relations for (de)hydrogenation over transition metals
journal, January 2011

  • Wang, S.; Petzold, V.; Tripkovic, V.
  • Physical Chemistry Chemical Physics, Vol. 13, Issue 46
  • DOI: 10.1039/c1cp20547a

Comment on “CO oxidation on ruthenium: The nature of the active catalytic surface” by D.W. Goodman, C.H.F. Peden, M.S. Chen
journal, December 2007


Reply to comment on “CO oxidation on ruthenium: The nature of the active catalytic surface” by H. Over, M. Muhler, A.P. Seitsonen
journal, December 2007


Experiment-Based Kinetic Monte Carlo Simulations: CO Oxidation over RuO 2 (110)
journal, December 2011

  • Farkas, A.; Hess, F.; Over, H.
  • The Journal of Physical Chemistry C, Vol. 116, Issue 1
  • DOI: 10.1021/jp204703p

Solving the master equation without kinetic Monte Carlo: Tensor train approximations for a CO oxidation model
journal, June 2016

  • Gelß, Patrick; Matera, Sebastian; Schütte, Christof
  • Journal of Computational Physics, Vol. 314
  • DOI: 10.1016/j.jcp.2016.03.025

Molecular origins of surface poisoning during CO oxidation over RuO2(110)
journal, August 2009


A review of multiscale modeling of metal-catalyzed reactions: Mechanism development for complexity and emergent behavior
journal, October 2011

  • Salciccioli, M.; Stamatakis, M.; Caratzoulas, S.
  • Chemical Engineering Science, Vol. 66, Issue 19
  • DOI: 10.1016/j.ces.2011.05.050

CO oxidation over Ru(0001) at near-atmospheric pressures: From chemisorbed oxygen to RuO2
journal, April 2009


Patched bimetallic surfaces are active catalysts for ammonia decomposition
journal, October 2015

  • Guo, Wei; Vlachos, Dionisios G.
  • Nature Communications, Vol. 6, Issue 1
  • DOI: 10.1038/ncomms9619

The gap-tooth method in particle simulations
journal, September 2003


A multiscale scheme for modeling catalytic flow reactors
journal, January 2006

  • Majumder, Debarshi; Broadbelt, Linda J.
  • AIChE Journal, Vol. 52, Issue 12
  • DOI: 10.1002/aic.11030

Efficient implementation of detailed surface chemistry into reactor models using mapped rate data
journal, April 2009


Adlayer inhomogeneity without lateral interactions: Rationalizing correlation effects in CO oxidation at RuO 2 (110) with first-principles kinetic Monte Carlo
journal, February 2011

  • Matera, Sebastian; Meskine, Hakim; Reuter, Karsten
  • The Journal of Chemical Physics, Vol. 134, Issue 6
  • DOI: 10.1063/1.3553258

To address surface reaction network complexity using scaling relations machine learning and DFT calculations
journal, March 2017

  • Ulissi, Zachary W.; Medford, Andrew J.; Bligaard, Thomas
  • Nature Communications, Vol. 8, Issue 1
  • DOI: 10.1038/ncomms14621

A Theoretical and Computational Analysis of Linear Free Energy Relations for the Estimation of Activation Energies
journal, July 2012

  • Sutton, Jonathan E.; Vlachos, Dionisios G.
  • ACS Catalysis, Vol. 2, Issue 8
  • DOI: 10.1021/cs3003269

Application of Wavelet-Based Methods for Accelerating Multi-Time-Scale Simulation of Bistable Heterogeneous Catalysis
journal, February 2017

  • Gur, Sourav; Frantziskonis, George N.; Pannala, Sreekanth
  • Industrial & Engineering Chemistry Research, Vol. 56, Issue 9
  • DOI: 10.1021/acs.iecr.6b04407

Multiscale modelling couples patches of non-linear wave-like simulations
journal, October 2015

  • Cao, Meng; Roberts, A. J.
  • IMA Journal of Applied Mathematics, Vol. 81, Issue 2
  • DOI: 10.1093/imamat/hxv034

Rationalizing the Relation between Adlayer Structure and Observed Kinetics in Catalysis
journal, February 2016


Local-metrics error-based Shepard interpolation as surrogate for highly non-linear material models in high dimensions
journal, October 2017

  • Lorenzi, Juan M.; Stecher, Thomas; Reuter, Karsten
  • The Journal of Chemical Physics, Vol. 147, Issue 16
  • DOI: 10.1063/1.4997286

SBH10: A Benchmark Database of Barrier Heights on Transition Metal Surfaces
journal, September 2017

  • Mallikarjun Sharada, Shaama; Bligaard, Thomas; Luntz, Alan C.
  • The Journal of Physical Chemistry C, Vol. 121, Issue 36
  • DOI: 10.1021/acs.jpcc.7b05677

Works referencing / citing this record:

A Critical Look at Direct Catalytic Hydrogenation of Carbon Dioxide to Olefins
journal, August 2019

  • Ronda‐Lloret, Maria; Rothenberg, Gadi; Shiju, N. Raveendran
  • ChemSusChem, Vol. 12, Issue 17
  • DOI: 10.1002/cssc.201900915

Electronic Structure, Morphological Aspects, Optical and Electrochemical Properties of RuO2 Nanocrystals
journal, July 2019


SQERT-T: alleviating kinetic Monte Carlo (KMC)-stiffness in transient KMC simulations
journal, June 2018

  • Savara, Aditya; Sutton, Jonathan E.
  • Journal of Physics: Condensed Matter, Vol. 30, Issue 29
  • DOI: 10.1088/1361-648x/aacb6d