Realistic multisite lattice-gas modeling and KMC simulation of catalytic surface reactions: Kinetics and multiscale spatial behavior for CO-oxidation on metal (100) surfaces

Abstract A realistic molecular-level description of catalytic reactions on single-crystal metal surfaces can be provided by stochastic multisite lattice-gas (msLG) models. This approach has general applicability, although in this report, we will focus on the example of CO-oxidation on the unreconstructed fcc metal (1 0 0) or M(1 0 0) surfaces of common catalyst metals M = Pd, Rh, Pt and Ir (i.e., avoiding regimes where Pt and Ir reconstruct). These models can capture the thermodynamics and kinetics of adsorbed layers for the individual reactants species, such as CO/M(1 0 0) and O/M(1 0 0), as well as the interaction and reaction between different reactant species in mixed adlayers, such as (CO + O)/M(1 0 0). The msLG models allow population of any of hollow, bridge, and top sites. This enables a more flexible and realistic description of adsorption and adlayer ordering, as well as of reaction configurations and configuration-dependent barriers. Adspecies adsorption and interaction energies, as well as barriers for various processes, constitute key model input. The choice of these energies is guided by experimental observations, as well as by extensive Density Functional Theory analysis. Model behavior is assessed via Kinetic Monte Carlo (KMC) simulation. We also address the simulation challenges and theoretical ramifications associated with very rapid diffusion and local equilibration of reactant adspecies such as CO. These msLG models are applied to describe adsorption, ordering, and temperature programmed desorption (TPD) for individual CO/M(1 0 0) and O/M(1 0 0) reactant adlayers. In addition, they are also applied to predict mixed (CO + O)/M(1 0 0) adlayer structure on the nanoscale, the complete bifurcation diagram for reactive steady-states under continuous flow conditions, temperature programmed reaction (TPR) spectra, and titration reactions for the CO-oxidation reaction. Extensive and reasonably successful comparison of model predictions is made with experimental data. Furthermore, we discuss the possible transition from traditional mean-field-type bistability and reaction kinetics for lower-pressure to multistability and enhanced fluctuation effects for moderate- or higher-pressure. Behavior in the latter regime reflects a stronger influence of adspecies interactions and also lower diffusivity in the higher-coverage mixed adlayer. We also analyze mesoscale spatiotemporal behavior including the propagation of reaction–diffusion fronts between bistable reactive and inactive states, and associated nucleation-mediated transitions between these states. This behavior is controlled by complex surface mass transport processes, specifically chemical diffusion in mixed reactant adlayers for which we provide a precise theoretical formulation. The msLG models together with an appropriate treatment of chemical diffusivity enable equation-free heterogeneous coupled lattice-gas (HCLG) simulations of spatiotemporal behavior. In addition, msLG + HCLG modeling can describe coverage variations across polycrystalline catalysts surfaces, pressure variations across catalyst surfaces in microreactors, and could be incorporated into a multiphysics framework to describe mass and heat transfer limitations for high-pressure catalysis.