Gas‐Phase Thermodynamics as a Validation of Computational Catalysis on Surfaces: A Case Study of Fischer–Tropsch Synthesis

Density functional theory has become a valuable tool to study surface catalysis. However, due to the scarcity of clean and reliable experimental data on surfaces, the theoretical methods employed to explore heterogeneous catalytic mechanisms are usually less well validated than those for gas‐phase reactions. We argue herein that gas‐phase reactions and the corresponding surface reactions are related through the Born–Haber cycle and computational catalysis on surfaces will be less meaningful if gas‐phase behavior cannot first be suitably determined. In this contribution, we have constructed a set of gas‐phase reactions relevant to the Fischer–Tropsch synthesis as a case study. With this set, we have tested the validity of the widely used PBE and B3LYP functionals and found that neither of them are capable of describing all kinds of gas‐phase reactions properly, such that some surface reactions may be biased falsely against the others. Significantly, XYG3, which is a double‐hybrid functional that includes Hartree–Fock‐like exchange and many‐body perturbation correlation effects, presents a significant improvement for all of the gas‐phase reactions, holding promise for further development for surface catalysis.     

Theoretical models of X–H bonds breaking (X = C,O, and H) over metal surfaces: Used for simulation of catalytic methane steam reforming

A short review of the works where theoretical models for describing kinetics of catalytic X–H bonds breaking reactions (X = C, O, and H) over metal surfaces were developed on the basis of concepts of the Dogonadze–Kuznetsov–Levich quantum mechanical theory of chemical processes. Numerical values of the rate constants of these reactions over (111) surfaces of nickel, platinum and rhodium, which are considered as steps of a complex catalytic process of methane steam reforming (MSR) are calculated and compared with experimental data. These rate constants are used for simulations of microkinetic models of the MSR reactions on the catalysts. Effects of external parameters on the MSR rates and on isotope effects are described.     

Rate‐Reactivity Model: A New Theoretical Basis for Systematic Kinetic Characterization of Heterogeneous Catalysts

Non‐steady‐state kinetic measurements contain a wealth of information about catalytic reactions and other gas–solid chemical interactions, which is extracted from experimental data via kinetic models. The standard mathematical framework of microkinetic models, which are typically used in computational catalysis and for advanced modeling of steady‐state data, encounters multiple challenges when applied to non‐steady‐state data. Robust phenomenological models, such as the steady‐state Langmuir–Hinshelwood–Hougen–Watson equations, are presently unavailable for non‐steady‐state data. Herein, a novel modeling framework is proposed to fulfill this need. The rate‐reactivity model (RRM) is formulated in terms of experimentally observable quantities including the gaseous transformation rates, concentrations, and surface uptakes. The model is linear with respect to these quantities and their pairwise products, and it is also linear in terms of its parameters (reactivities). The RRM parameters have a clear physicochemical meaning and fully characterize the kinetic behavior of a specific catalyst state, but unlike microkinetic models that rely on hypothetical surface intermediates and specific reaction networks, the RRM does not require any assumptions regarding the underlying mechanism. The systematic RRM‐based procedure outlined in this paper enables an effective comparison of various catalysts and the construction of more detailed microkinetic models in a rational manner. The model was applied to temporal analysis of products pulse‐response data as an example, but it is more generally applicable to other non‐steady‐state techniques that provide time‐resolved rates and concentrations. Several numerical examples are given to illustrate the application of the model to simple model reactions.     

不同金属催化水煤气变换反应活性的Monte Carlo模拟研究

运用BOC-MP方法对Cu(110),Cu(111),Pd(111)和Au(111)等过渡金属催化的WGS反应的可能微观动力学步骤进行了详尽的能学数据计算,并结合MonteCarlo方法对WGS反应的表面氧化还原机理进行了计算机模拟。结果表明,Cu的催化活性优于Pd,Au的催化活性,并获得了相应金属上WGS反应的表观活化能及动力学指前因子(相对值);在此基础上,对该反应的结构敏感性进行了研究,发现该反应为一结构敏感反应,与实验结果相符。     

Automated Generation of Microkinetics for Heterogeneously Catalyzed Reactions Considering Correlated Uncertainties**

The study presents an ab-initio based framework for the automated construction of microkinetic mechanisms considering correlated uncertainties in all energetic parameters and estimation routines. 2000 unique microkinetic models were generated within the uncertainty space of the BEEF-vdW functional for the oxidation reactions of representative exhaust gas emissions from stoichiometric combustion engines over Pt(111) and compared to experiments through multiscale modeling. The ensemble of simulations stresses the importance of considering uncertainties. Within this set of first-principles-based models, it is possible to identify a microkinetic mechanism that agrees with experimental data. This mechanism can be traced back to a single exchange-correlation functional, and it suggests that Pt(111) could be the active site for the oxidation of light hydrocarbons. The study provides a universal framework for the automated construction of reaction mechanisms with correlated uncertainty quantification, enabling a DFT-constrained microkinetic model optimization for other heterogeneously catalyzed systems.     

First Principles Calculations for Hydrogenation of Acrolein on Pd and Pt: Chemoselectivity Depends on Steric Effects on the Surface

The chemoselective hydrogenation of acrolein on Pt(111) and Pd(111) surfaces is investigated employing density functional theory calculations. The computed potential energy surfaces together with the analysis of reaction mechanisms demonstrate that steric effects are an important factor that governs chemoselectivity. The reactions at the C=O functionality require more space than the reactions at the C=C functionality. Therefore the formation of allyl alcohol is more favorable at low coverage, while the reduction of the C=C bond and the formation of propanal becomes kinetically more favorable at higher coverage. The elementary reaction steps are found to follow different reaction mechanisms, which are identified according to terminology typically used in organometallic catalysis. The transition state scaling (TSS) relationship is demonstrated and the origin of multiple TSS lines is linked to variation of an internal electronic structure of a carbon skeleton.     

Sonochemical Reaction of Bifunctional Molecules on Silicon (111) Hydride Surface

While the sonochemical grafting of molecules on silicon hydride surface to form stable Si–C bond via hydrosilylation has been previously described, the susceptibility towards nucleophilic functional groups during the sonochemical reaction process remains unclear. In this work, a competitive study between a well-established thermal reaction and sonochemical reaction of nucleophilic molecules (cyclopropylamine and 3-Butyn-1-ol) was performed on p-type silicon hydride (111) surfaces. The nature of surface grafting from these reactions was examined through contact angle measurements, X-ray photoelectron spectroscopy (XPS) and atomic force microscopy (AFM). Cyclopropylamine, being a sensitive radical clock, did not experience any ring-opening events. This suggested that either the Si–H may not have undergone homolysis as reported previously under sonochemical reaction or that the interaction to the surface hydride via a lone-pair electron coordination bond was reversible during the process. On the other hand, silicon back-bond breakage and subsequent surface roughening were observed for 3-Butyn-1-ol at high-temperature grafting (≈150 °C). Interestingly, the sonochemical reaction did not produce appreciable topographical changes to surfaces at the nano scale and the further XPS analysis may suggest Si–C formation. This indicated that while a sonochemical reaction may be indifferent towards nucleophilic groups, the surface was more reactive towards unsaturated carbons. To the best of the author’s knowledge, this is the first attempt at elucidating the underlying reactivity mechanisms of nucleophilic groups and unsaturated carbon bonds during sonochemical reaction of silicon hydride surfaces.     

Dynamic interplay between diffusion and reaction: nitrogen recombination on Rh{211} in car exhaust catalysis

The adsorption and diffusion of atomic nitrogen on Rh{211} as well as formation and desorption of molecular nitrogen from this surface have been investigated by means of density functional theory (DFT) calculations. The elementary step reaction mechanism derived from this comprehensive DFT study forms the foundation of a detailed microkinetic model including diffusion, recombination, and desorption of nitrogen species. It will be shown that nitrogen formation on a stepped rhodium surface is a dynamic interplay of atomic nitrogen diffusion and reaction. Moreover, evidence will be presented that not one but several on-step recombination reactions are responsible for dinitrogen formation and desorption.     

Interfacial energy exchange and reaction dynamics in collisions of gases on model organic surfaces

Molecular beam scattering experiments and molecular dynamics simulations have been combined to develop an atomic-level understanding of energy transfer, accommodation, and reactions during collisions between gases and model organic surfaces. The work highlighted in this progress report has been motivated by the scientific importance of understanding fundamental interfacial chemical reactions and the relevance of reactions on organic surfaces to many areas of environmental chemistry. The experimental investigations have been accomplished by molecular beam scattering from ω-functionalized self-assembled monolayers (SAMs) on gold. Molecular beams provide a source of reactant molecules with precisely characterized collision energy and flux; SAMs afford control over the order, structure, and chemical nature of the surface. The details of molecular motion that affect energy exchange and scattering have been elucidated through classical-trajectory simulations of the experimental data using potential energy surfaces derived from ab initio calculations. Our investigations began by employing rare-gas scattering to explore how alkanethiol chain length and packing density, terminal group relative mass, orientation, and chemical functionality influence energy transfer and accommodation at organic surfaces. Subsequent studies of small molecule scattering dynamics provided insight into the influence of internal energy, molecular orientation, and gas–surface attractive forces in interfacial energy exchange. Building on the understanding of scattering dynamics in non-reactive systems, our work has recently explored the reaction probabilities and mechanisms for O₃ and atomic fluorine in collisions with a variety of functionalized SAM surfaces. Together, this body of work has helped construct a more comprehensive understanding of reaction dynamics at organic surfaces.     

Hydrogenation of CO on a silica surface: an embedded cluster approach

The sequential addition of H atoms to CO adsorbed on a siliceous edingtonite surface is studied with an embedded cluster approach, using density functional theory for the quantum mechanical (QM) cluster and a molecular force field for the molecular mechanical (MM) cluster. With this setup, calculated QM/MM adsorption energies are in agreement with previous calculations employing periodic boundary conditions. The catalytic effect of the siliceous edingtonite (100) surface on CO hydrogenation is assessed because of its relevance to astrochemistry. While adsorption of CO on a silanol group on the hydroxylated surface did not reduce the activation energy for the reaction with a H atom, a negatively charged defect on the surface is found to reduce the gas phase barriers for the hydrogenation of both CO and H2C=O. The embedded cluster approach is shown to be a useful and flexible tool for studying reactions on (semi-)ionic surfaces and specific defects thereon. The methodology presented here could easily be applied to study reactions on silica surfaces that are of relevance to other scientific areas, such as biotoxicity of silica dust and geochemistry.     

STM studies of photochemistry and plasmon chemistry on metal surfaces

We review our recent studies of photochemistry and plasmon chemistry of dimethyl disulfide, (CH₃S)₂, molecules adsorbed on metal surfaces using a scanning tunneling microscope (STM). The STM has been used not only for the observation of surface structures at atomic spatial resolution but also for local spectroscopies. The STM combined with optical excitation by light can be employed to investigate chemical reactions of single molecules induced by photons and localized surface plasmons. This technique allows us to gain insights into reaction mechanisms at a single molecule level. The experimental procedures to examine the chemical reactions using the STM are briefly described. The mechanism for the photodissociation reaction of (CH₃S)₂ molecules adsorbed on metal surfaces is discussed based on both the experimental results obtained with the STM and the electronic structures calculated by density functional theory. The dissociation reaction of the (CH₃S)₂ molecule induced by the optically excited plasmon in the STM junction between a Ag tip and metal substrate is also described. The reaction mechanism and pathway of this plasmon-induced chemical reaction are discussed by comparison with those proposed in plasmon chemistry.     

Gas‐phase reaction mechanism of Pd+ with CH3CHO: A density functional theoretical study

Details on the reactions of: (1) Pd⁺ + CH₃CHO → PdCO⁺ + CH₄ and (2) Pd⁺ + CH₃CHO → PdH + CH₃CO⁺ in the gas phase were investigated using density functional theory (B3LYP), in conjunction with the LANL2DZ+6‐311+G(d) basis set. Three encounter complexes were located on the potential energy surfaces and the calculations indicated that both the C? C and aldehyde C? H bond activation of acetaldehyde could lead to the dominant demethanation reaction. The charge transfer process for PdH abstraction was caused by an intramolecular PdH rearrangement of the newly found η¹‐aldehyde attached complex. © 2010 Wiley Periodicals, Inc. Int J Quantum Chem, 2011     

硅苯与HX (X=F,OH, NH2)的1,2-及1,4-加成反应

采用密度泛函理论方法在B3LYP/6-311++G(d,p)水平上, 研究了硅苯与HX (X=F, OH, NH₂)的1,2-及1,4-加成反应的微观机理和势能剖面, 考察了Si 原子上的取代基及四氢呋喃溶剂对反应势能剖面的影响. 研究结果表明, 标题反应有两种可能的机理: (1) 硅苯与一个HX (X=F, OH, NH₂)分子先形成中间复合物, 然后经过四元环过渡态(机理1)生成最终产物; (2) 硅苯与两个HX分子先形成中间复合物, 然后经过六元环过渡态(机理2)生成另一中间复合物, 该中间复合物脱去一个HX分子形成最终产物. 机理2 在动力学上远较机理1 有利. 1,2-及1,4-加成产物哪种优先形成由动力学控制且与X基团的种类有关. HX在气相中参与加成反应从易到难的次序为: HF>H₂O>NH₃. Si 原子上具有较强供电子和吸电子性质的取代基, 在热力学和动力学上均有利于反应的进行, 但具有较大体积的2,4,6-三甲基苯基取代基对反应反而不利. 四氢呋喃溶剂在热力学上不利于硅苯与HX的1,2-及1,4-加成反应, 在动力学上对HF或H₂O作为加成试剂的反应也不利, 但对NH3作为加成试剂的反应反而有利.     

Interpolating DFT Data for 15D Modeling of Methane Dissociation on an fcc Metal

Detailed simulation of reactions occurring on and with the surfaces of crystalline materials usually require a continuous representation of the potential energy surface that describes the adsorbate–surface interaction. Only a few techniques are available to describe interactions with polyatomic adsorbates that respect all of the symmetries of the interactions. The modified Shepard interpolation has recently been reformulated to ensure symmetries are rigorously imposed. In this work, the modified Shepard interpolation is used to construct a 15D potential energy surface for the reaction of methane with the {100} surface of a face‐centered cubic metal, in the Born‐‐Oppenheimer static surface (BOSS) approximation. The energy of the system is calculated using density functional theory (DFT), and the geometries around which the potential is expanded are selected by quasi‐classical trajectory calculations. The energy of the resulting continuous potential energy surface exactly matched the DFT energy at these points; there is no fitting error. It is demonstrated that the classical reaction probability converges with a reasonable number of interpolation points for this 15D system.     

Classical, phenomenological analysis of the kinetics of reactions at the gas-exposed surface of mixed ionic electronic conductors

The kinetics of reactions occurring at the gas-exposed surfaces of charged mixed ionic electronic conductors (MIECs) are examined from theoretical first principles. Analysis based on the classical electrochemical potential-transition state theory model reveals that the nature of the reactions is electrochemical in general. However, the influence of the surface potential on the reaction rate is opposite for adsorption and incorporation reactions. Two-dimensional finite volume models of an MIEC as working electrode in a half-cell configuration are presented. The results for a simple, two-step reduction process show that the effect of the surface potential on the rate of reactions is minimal for incorporation-limited reactions but more influential for adsorption-limited reactions. An erratum to this article is available at .     

Analysis of Nascent Rotational Energy Distributions and Reaction Mechanisms of the Gas‐Phase Radical–Radical Reaction O(3P)+(CH3)2CH→C3H6+OH

This paper reports on the gas‐phase radical–radical dynamics of the reaction of ground‐state atomic oxygen [O(³P), from the photodissociation of NO₂] with secondary isopropyl radicals [(CH₃)₂CH, from the supersonic flash pyrolysis of isopropyl bromide]. The major reaction channel, O(³P)+(CH₃)₂CH→C₃H₆ (propene)+OH, is examined by high‐resolution laser‐induced fluorescence spectroscopy in crossed‐beam configuration. Population analysis shows bimodal nascent rotational distributions of OH (X²Π) products with low‐ and high‐N′′ components in a ratio of 1.25:1. No significant spin–orbit or Λ‐doublet propensities are exhibited in the ground vibrational state. Ab initio computations at the CBS‐QB3 theory level and comparison with prior theory show that the statistical method is not suitable for describing the main reaction channel at the molecular level. Two competing mechanisms are predicted to exist on the lowest doublet potential‐energy surface: direct abstraction, giving the dominant low‐N′′ components, and formation of short‐lived addition complexes that result in hot rotational distributions, giving the high‐N′′ components. The observed competing mechanisms contrast with previous bulk kinetic experiments conducted in a fast‐flow system with photoionization mass spectrometry, which suggested a single abstraction pathway. In addition, comparison of the reactions of O(³P) with primary and tertiary hydrocarbon radicals allows molecular‐level discussion of the reactivity and mechanism of the title reaction.     

An unexpected pathway for the catalytic oxidation of methylidyne on Rh{111} as a route to syngas

This study investigates the adsorption properties of methylidyne (CH) on Rh{111}, its partial and full oxidation as well as its surface mobility, by means of plane-wave density functional theory (DFT) calculations. Besides investigating known oxidation pathways on rhodium, such as decomposition of CH and subsequent oxidation of the decomposition products, new pathways such as direct reaction of methylidyne and oxygen toward a surface aldyhyde-type species and the decomposition of this species are considered. The unexpected and novel pathway determined here by DFT is utilized for a microkinetic model of the formation of CO and CO2 from methylidyne. A comparison of this microkinetic study with experimental data shows that our novel mechanism can indeed describe the observations. This comparison strongly suggests that this new alternative route is the main reaction pathway for the conversion of methylidyne.     

Microkinetic model for NO–CO reaction: Model reduction

The objective of this work is to elucidate controlling mechanisms in NO_x reduction, develop reduced‐order reaction models, and analyze the reactor performance using the reduced‐order reaction model for the NO–CO reaction. We start with the microkinetic model on platinum, which describes the mechanism of catalytic reduction of NO by CO. The formation of the main product N₂O and the competitive formation of the side product N₂ are accounted for in the microkinetic model. Sensitivity and reaction path analysis have been carried out to determine the rate‐limiting steps as well as the most abundant reactive intermediates in the system. Owing to the differences between system performance at high and low temperatures, the model has been analyzed in detail in these temperature regimes. Two closed‐form expressions, corresponding to the two global reactions involved, have been derived. The characteristic features of the microkinetic model such as the sharp increase in NO conversion and the selectivity to N₂O are captured well by the reduced model. The reduced‐order model has been extended to the rhodium catalyst as well. © 2012 Wiley Periodicals, Inc. Int J Chem Kinet 44: 577–585, 2012     

On the mechanism of low-temperature water gas shift reaction on copper

Periodic, self-consistent density functional theory (DFT-GGA) calculations are used to investigate the water gas shift reaction (WGSR) mechanism on Cu(111). The thermochemistry and activation energy barriers for all the elementary steps of the commonly accepted redox mechanism, involving complete water activation to atomic oxygen, are presented. Through our calculations, we identify carboxyl, a new reactive intermediate, which plays a central role in WGSR on Cu(111). The thermochemistry and activation energy barriers of the elementary steps of a new reaction path, involving carboxyl, are studied. A detailed DFT-based microkinetic model of experimental reaction rates, accounting for both the previous and the new WGSR mechanism show that, under relevant experimental conditions, (1) the carboxyl-mediated route is the dominant path, and (2) the initial hydrogen abstraction from water is the rate-limiting step. Formate is a stable "spectator" species, formed predominantly through CO2 hydrogenation. In addition, the microkinetic model allows for predictions of (i) surface coverage of intermediates, (ii) WGSR apparent activation energy, and (iii) reaction orders with respect to CO, H2O, CO2, and H2.