Importance of metal-oxide interfaces in heterogeneous catalysis: A combined DFT, microkinetic, and experimental study of water-gas shift on Au/MgO

Uncertainty-Aware First-Principles Exploration of Chemical Reaction Networks

Exploring large chemical reaction networks with automated exploration approaches and accurate quantum chemical methods can require prohibitively large computational resources. Here, we present an automated exploration approach that focuses on the kinetically relevant part of the reaction network by interweaving (i) large-scale exploration of chemical reactions, (ii) identification of kinetically relevant parts of the reaction network through microkinetic modeling, (iii) quantification and propagation of uncertainties, and (iv) reaction network refinement. Such an uncertainty-aware exploration of kinetically relevant parts of a reaction network with automated accuracy improvement has not been demonstrated before in a fully quantum mechanical approach. Uncertainties are identified by local or global sensitivity analysis. The network is refined in a rolling fashion during the exploration. Moreover, the uncertainties are considered during kinetically steering of a rolling reaction network exploration. We demonstrate our approach for Eschenmoser-Claisen rearrangement reactions. The sensitivity analysis identifies that only a small number of reactions and compounds are essential for describing the kinetics reliably, resulting in efficient explorations without sacrificing accuracy and without requiring prior knowledge about the chemistry unfolding.

First principles simulations of MgO(100) surface hydration at ambient conditions

Developing a better understanding of water ordering and hydroxylation at oxide mineral surfaces is important across a breath of application spaces. Recent vibrational sum frequency generation (vSFG) measurements on MgO(100) surfaces at ambient conditions showed that water dissociates and hydroxylates the surface yielding a non-hydrogen bonded hydroxyl species. Starting from previously determined water hydroxylation patterns on MgO(100), we performed ab initio thermodynamic calculations and vibrational analysis to compare with the vSFG observations. At ambient conditions (i.e., T = 298.15 K and pH2O = 32 mbar), the most thermodynamically favorable surface hydroxylation is found to be p(3 × 2) - 8H2O, involving a dissociation of 25% of the adsorbed water. Analysis of the vibrational density of states for this hydroxylation configuration yielded three different hydrogen bonding environments with the frequency of the peaks in very good agreement with the vSFG measurements. However, the non-H-bonded spectral feature on this surface is predicted to be similar to that expected for Mg(OH)2, a thermodynamically downhill alteration of the surface that must be independently ruled out before one can be fully confident in the apparent theory/vSFG agreement. Our study provides more insights into the ordering and structure of water monolayer at MgO(100) surface at ambient conditions and completes previous theoretical and experimental analysis performed at low temperature and ultra-high vacuum conditions.

Theoretical insight into the promotion effect of potassium additive on the water-gas shift reaction over low-coordinated Au catalysts

CONTEXT

How to elucidate the effect of alkali metal promoters on gold-catalyzed water-gas shift reaction intrinsically remains a challenging, because that the complex synergy effects such as strong metal-support interactions, interfacial effects, and charge transfer of supported metal catalysts makes people difficulty in the understanding the alkali promotion phenomenon in nature. Herein, we report a systematically study of whole water-gas shift reaction mechanism on pure and the K-modified defected-Au(211) (i.e., by removing one surface Au atom from perfect Au(211) and make one model with the Au-Au coordination number is six) by using the microkinetic modeling based on first principles. Our results indicate that the presence of K can increase the adsorption ability of oxygen-containing species via the attractive coulomb interaction, has no significant effect on the adsorption of H species, but inhibits the adsorption of CO due to the steric effect. K promoter stabilizes the water adsorption by ~0.3 eV, which results in one order increasing of whole reaction rate. Interestingly, the strong promotion effect of the K can be assigned to the significant direct space interaction between K and the adsorbate H2O* through the inducted electric field, which can be further confirmed by the posed negative electric field on the unpromoted D-Au(211). Microkinetic modeling results revealed that the carboxyl mechanism is the most likely to occur, redox mechanism is the next one, and the formate mechanism is the least likely to occur. For different kinds of alkali metal additives, the adsorption strength of water molecules gradually weakens from Li to Cs, but Na shows the best promoter behavior at the low temperature. By considering the effect of K contents on the reactivity of water-gas shift reaction, we found that the K with the medium coverage (~0.2~0.3 ML) has the strongest promoting effect. It is expected that the conclusion of this work can be extended to other WGSR catalytic systems like Cu(or Pt).

METHODS

All calculations were performed by using the plane-wave based periodic method implemented in Vienna ab initio simulation package (VASP, version 5.4.4), where the ionic cores are described by the projector augmented wave (PAW) method. The exchange and correlation energies were computed using the Perdew, Burke and Ernzerhof functional with the vdw correction (PBE-D3). The transition states (TSs) were searched using the climbing image nudged elastic band (CI-NEB) method. Some electronic structure properties like work function was predicated by the DS-PAW software. Microkinetic simulation was carried out using MKMCXX software.

Characterization and Simulation of Nanoscale Catastrophic Failure of Metal/Ceramic Interfaces

The catastrophic failure of metal/ceramic interfaces is a complex process involving the energy transfer between accumulated elastic strain energy and many types of energy dissipation. To quantify the contribution of bulk and interface cohesive energy to the interface cleavage fracture without global plastic deformation, we characterized the quasi-static fracture process of both coherent and semi-coherent fcc-metal/MgO(001) interface systems using a spring series model and molecular static simulations. Our results show that the theoretical catastrophe point and spring-back length by the spring series model are basically consistent with the simulation results of the coherent interface systems. For defect interfaces with misfit dislocations, atomistic simulations revealed an obvious interface weakening effect in terms of reduced tensile strength and work of adhesion. As the model thickness increases, the tensile failure behaviors show significant scale effects-thick models tend to catastrophic failure with abrupt stress drop and obvious spring-back phenomenon. This work provides insight into the origin of catastrophic failure at metal/ceramic interfaces, which highlights a pathway by combining the material and structure design to improve the reliability of layered metal-ceramic composites.

Metal-Support Interactions in Heterogeneous Catalysis: DFT Calculations on the Interaction of Copper Nanoparticles with Magnesium Oxide

Oxide supports play an important role in enhancing the catalytic properties of transition metal nanoparticles in heterogeneous catalysis. How extensively interactions between the oxide support and the nanoparticles impact the electronic structure as well as the surface properties of the nanoparticles is hence of high interest. In this study, the influence of a magnesium oxide support on the properties of copper nanoparticles with different size, shape, and adsorption sites is investigated using density functional theory (DFT) calculations. By proposing simple models to reduce the cost of the calculations while maintaining the accuracy of the results, we show using the nonreducible oxide support MgO as an example that there is no significant influence of the MgO support on the electronic structure of the copper nanoparticles, with the exception of adsorption directly at the Cu-MgO interface. We also propose a simplified methodology that allows us to reduce the cost of the calculations, while the accuracy of the results is maintained. We demonstrate in addition that the Cu nanowire model corresponds well to the nanoparticle model, which reduces the computational cost even further.

Theoretical Calculations on Metal Catalysts Toward Water-Gas Shift Reaction: a Review

Water-gas shift (WGS) reaction offers a dominating path to hydrogen generation from fossil fuel, in which heterogeneous metal catalysts play a crucial part in this course. This review highlights and summarizes recent developments on theoretical calculations of metal catalysts developed to date, including surface structure (e. g., monometallic and polymetallic systems) and interface structure (e. g., supported catalysts and metal oxide composites), with special emphasis on the characteristics of crystal-face effect, alloying strategy, and metal-support interaction. A systematic summarization on reaction mechanism was performed, including redox mechanism, associative mechanism as well as hybrid mechanism; the development on chemical kinetics (e. g., molecular dynamics, kinetic Monte Carlo and microkinetic simulation) was then introduced. At the end, challenges associated with theoretical calculations on metal catalysts toward WGS reaction are discussed and some perspectives on the future advance of this field are provided.

Enhancing Tungsten Oxide Gasochromism with Noble Metal Nanoparticles: The Importance of the Interface

Crystalline tungsten trioxide (WO₃ ) thin films covered by noble metal (gold and platinum) nanoparticles are synthesized via wet chemistry and used as optical sensors for gaseous hydrogen. Sensing performances are strongly influenced by the catalyst used, with platinum (Pt) resulting as best. Surprisingly, it is found that gold (Au) can provide remarkable sensing activity that tuned out to be strongly dependent on the nanoparticle size: devices sensitized with smaller nanoparticles display better H₂ sensing performance. Computational insight based on density functional theory calculations suggested that this can be related to processes occurring specifically at the Au nanoparticle-WO₃ interface (whose extent is in fact dependent on the nanoparticle size), where the hydrogen dissociative adsorption turns out to be possible. While both experiments and calculations single out Pt as better than Au for sensing, the present work reveals how an exquisitely nanoscopic effect can yield unexpected sensing performance for Au on WO₃ , and how these performances can be tuned by controlling the nanoscale features of the system.

Differences in solvation thermodynamics of oxygenates at Pt/Al2O3 perimeter versus Pt(111) terrace sites

A prominent role of water in aqueous-phase heterogeneous catalysis is to modify free energies; however, intuition about how is based largely on pure metal surfaces or even homogeneous solutions. Using multiscale modeling with explicit liquid water molecules, we show that the influence of water on the free energies of adsorbates at metal/support interfaces is different than that on pure metal surfaces. We specifically compute free energies of solvation for methanol and its constituents on a Pt/Al₂O₃ catalyst and compare the results to analogous values calculated on a pure Pt catalyst. We find that the more hydrophilic Pt/Al₂O₃ interface leads to smaller (more positive) free energies of solvation due to an increased entropy penalty resulting from the additional work necessary to disrupt the interfacial water structure and accommodate the interfacial species. The results will be of interest in other fields, including adsorption and proteins.

Predicting Frequency from the External Chemical Environment: OH Vibrations on Hydrated and Hydroxylated Surfaces

Robust correlation curves are essential to decipher structural information from IR-vibrational spectra. However, for surface-adsorbed water and hydroxides, few such correlations have been presented in the literature. In this paper, OH vibrational frequencies are correlated against 12 structural descriptors representing the quantum mechanical or geometrical environment, focusing on those external to the vibrating molecule. A nonbiased fitting procedure based on Gaussian process regression (GPR) was used alongside simple analytical functional forms. The training data consist of 217 structurally unique OH groups from 38 water/metal oxide interface systems for MgO, CaO and CeO₂, all optimized at the DFT level, and the fully anharmonic and uncoupled OH vibrational signatures were calculated. Among our results, we find the following: (i) The intermolecular R(H···O) hydrogen bond distance is particularly strong, indicating the primary cause of the frequency shift. (ii) Similarly, the electric field along the H-bond vector is also a good descriptor. (iii) Highly detailed machine learning descriptors (ACSF, SOAP) are less intuitive but were found to be more capable descriptors. (iv) Combinations of geometric and QM descriptors give the best predictions, supplying complementary information.

Bifunctional hydroformylation on heterogeneous Rh-WOx pair site catalysts

Metal-catalysed reactions are often hypothesized to proceed on bifunctional active sites, whereby colocalized reactive species facilitate distinct elementary steps in a catalytic cycle^1-8. Bifunctional active sites have been established on homogeneous binuclear organometallic catalysts^9-11. Empirical evidence exists for bifunctional active sites on supported metal catalysts, for example, at metal-oxide support interfaces^2,6,7,12. However, elucidating bifunctional reaction mechanisms on supported metal catalysts is challenging due to the distribution of potential active-site structures, their dynamic reconstruction and required non-mean-field kinetic descriptions^7,12,13. We overcome these limitations by synthesizing supported, atomically dispersed rhodium-tungsten oxide (Rh-WO_x) pair site catalysts. The relative simplicity of the pair site structure and sufficient description by mean-field modelling enable correlation of the experimental kinetics with first principles-based microkinetic simulations. The Rh-WO_x pair sites catalyse ethylene hydroformylation through a bifunctional mechanism involving Rh-assisted WO_x reduction, transfer of ethylene from WO_x to Rh and H₂ dissociation at the Rh-WO_x interface. The pair sites exhibited >95% selectivity at a product formation rate of 0.1 g_propanal cm^-3 h^-1 in gas-phase ethylene hydroformylation. Our results demonstrate that oxide-supported pair sites can enable bifunctional reaction mechanisms with high activity and selectivity for reactions that are performed in industry using homogeneous catalysts.

Rational design of highly efficient MXene-based catalysts for the water-gas-shift reaction

Water molecules linked by hydrogen bonds are responsible for the high efficiency of bi-functional catalysts for the water-gas-shift (WGS) reaction because water can act as a proton transfer medium. Herein, we propose an associative pathway for the WGS reaction assisted by water to realize hydrogen production. Based on this pathway, we show by first-principles calculations that a large family of oxygen-terminated two-dimensional transition metal carbides and nitrides (MXenes) deposited on Au clusters are promising catalysts for the WGS reaction. Remarkably, the rate-determining barriers for *CO → *COOH on Au/M_n+1X_nO₂ are in the range from 0.15 eV to 0.39 eV, indicating that WGS can occur at much lower temperatures. Furthermore, a comprehensive microkinetic model is constructed to describe the turnover frequencies (TOF) for the product under the steady-state conditions. More importantly, there is a perfect linear scaling relationship between the rate-determining barriers of the WGS and the free energy of the adsorbed hydrogen. Besides, the potential energy diagrams for CO reforming reveal that the F terminations introduced in experiments have only a slight influence on the catalytic performance of the oxygen-terminated MXenes. Our work not only opens a new avenue towards the WGS reaction but also provides many ideal catalysts for hydrogen production.

Influence of a Cu-zirconia interface structure on CO2 adsorption and activation

CO₂ adsorption and activation on a catalyst are key elementary steps for CO₂ conversion to various valuable products. In the present computational study, we screened different Cu-ZrO₂ interface structures and analyzed the influence of the interface structure on CO₂ binding strength using density functional theory calculations. Our results demonstrate that a Cu nanorod favors one position on both tetragonal and monoclinic ZrO₂ surfaces, where the bottom Cu atoms are placed close to the lattice oxygens. In agreement with previous calculations, we find that CO₂ prefers a bent bidentate configuration at the Cu-ZrO₂ interface and the molecule is clearly activated being negatively charged. Straining of the Cu nanorod influences CO₂ adsorption energy but does not change the preferred nanorod position on zirconia. Altogether, our results highlight that CO₂ adsorption and activation depend sensitively on the chemical composition and atomic structure of the interface used in the calculations. This structure sensitivity may potentially impact further catalytic steps and the overall computed reactivity profile.

Vibrational energy relaxation of interfacial OH on a water-covered α-Al2O3(0001) surface: a non-equilibrium ab initio molecular dynamics study

Vibrational relaxation of adsorbates is a sensitive tool to probe energy transfer at gas/solid and liquid/solid interfaces. The most direct way to study relaxation dynamics uses time-resolved spectroscopy. Here we report on a non-equilibrium ab initio molecular dynamics (NE-AIMD) methodology to model vibrational relaxation of OH vibrations on a hydroxylated, water-covered α-Al₂O₃(0001) surface. In our NE-AIMD approach, after exciting selected O-H bonds their coupling to surface phonons and to the water adlayer is analyzed in detail, by following both the energy flow in time, as well as the time-evolution of Vibrational Density of States (VDOS) curves. The latter are obtained from Time-dependent Correlation Functions (TCFs) and serve as prototypical, generic representatives of time-resolved vibrational spectra. As most important results, (i) we find a few-picosecond lifetime of the excited modes and (ii) identify both hydrogen-bonded aluminols and water molecules in the adsorbed water layer as main dissipative channels, while the direct coupling to Al₂O₃ surface phonons is of minor importance on the timescales of interest. Our NE-AIMD/TCF methodology is powerful for complex adsorbate systems, in principle even reacting ones, and opens a way towards time-resolved vibrational spectroscopy.

Microkinetic Modeling: A Tool for Rational Catalyst Design

The design of heterogeneous catalysts relies on understanding the fundamental surface kinetics that controls catalyst performance, and microkinetic modeling is a tool that can help the researcher in streamlining the process of catalyst design. Microkinetic modeling is used to identify critical reaction intermediates and rate-determining elementary reactions, thereby providing vital information for designing an improved catalyst. In this review, we summarize general procedures for developing microkinetic models using reaction kinetics parameters obtained from experimental data, theoretical correlations, and quantum chemical calculations. We examine the methods required to ensure the thermodynamic consistency of the microkinetic model. We describe procedures required for parameter adjustments to account for the heterogeneity of the catalyst and the inherent errors in parameter estimation. We discuss the analysis of microkinetic models to determine the rate-determining reactions using the degree of rate control and reversibility of each elementary reaction. We introduce incorporation of Brønsted-Evans-Polanyi relations and scaling relations in microkinetic models and the effects of these relations on catalytic performance and formation of volcano curves are discussed. We review the analysis of reaction schemes in terms of the maximum rate of elementary reactions, and we outline a procedure to identify kinetically significant transition states and adsorbed intermediates. We explore the application of generalized rate expressions for the prediction of optimal binding energies of important surface intermediates and to estimate the extent of potential rate improvement. We also explore the application of microkinetic modeling in homogeneous catalysis, electro-catalysis, and transient reaction kinetics. We conclude by highlighting the challenges and opportunities in the application of microkinetic modeling for catalyst design.

Computational Methods in Heterogeneous Catalysis

The unprecedented ability of computations to probe atomic-level details of catalytic systems holds immense promise for the fundamentals-based bottom-up design of novel heterogeneous catalysts, which are at the heart of the chemical and energy sectors of industry. Here, we critically analyze recent advances in computational heterogeneous catalysis. First, we will survey the progress in electronic structure methods and atomistic catalyst models employed, which have enabled the catalysis community to build increasingly intricate, realistic, and accurate models of the active sites of supported transition-metal catalysts. We then review developments in microkinetic modeling, specifically mean-field microkinetic models and kinetic Monte Carlo simulations, which bridge the gap between nanoscale computational insights and macroscale experimental kinetics data with increasing fidelity. We finally review the advancements in theoretical methods for accelerating catalyst design and discovery. Throughout the review, we provide ample examples of applications, discuss remaining challenges, and provide our outlook for the near future.

Selective Halogenation of Pyridines Using Designed Phosphine Reagents

Halopyridines are key building blocks for synthesizing pharmaceuticals, agrochemicals, and ligands for metal complexes, but strategies to selectively halogenate pyridine C-H precursors are lacking. We designed a set of heterocyclic phosphines that are installed at the 4-position of pyridines as phosphonium salts and then displaced with halide nucleophiles. A broad range of unactivated pyridines can be halogenated, and the method is viable for late-stage halogenation of complex pharmaceuticals. Computational studies indicate that C-halogen bond formation occurs via an S_NAr pathway, and phosphine elimination is the rate-determining step. Steric interactions during C-P bond cleavage account for differences in reactivity between 2- and 3-substituted pyridines.

Water Poisons H2 Activation at the Au-TiO2 Interface by Slowing Proton and Electron Transfer between Au and Titania

Understanding the dynamic changes at the active site during catalysis is a fundamental challenge that promises to improve catalytic properties. While performing Arrhenius studies during H₂ oxidation over Au/TiO₂ catalysts, we found different apparent activation energies (E_app) depending on the feedwater pressure. This is partially attributed to changing numbers of metal-support interface (MSI) sites as water coverage changes with temperature. Constant water coverage studies showed two kinetic regimes: fast heterolytic H₂ activation directly at the MSI (E_app ∼ 25 kJ/mol) and significantly slower heterolytic H₂ activation mediated by water (E_app ∼ 45 kJ/mol). The two regimes had significantly different kinetics, suggesting a complicated mechanism of water poisoning. Density functional theory (DFT) showed water has minor effects on the reaction thermodynamics, primarily attributable to intrinsic differences in surface reactivity of different Au sites in the DFT model. The DFT model suggested significant surface restructuring of the TiO₂ support during heterolytic H₂ adsorption; evidence for this phenomenon was observed during in situ infrared spectroscopy experiments. A monolayer of water on the hydroxylated TiO₂ surface increased the H₂ dissociation activation barrier by ∼0.2 eV, in good agreement the difference in experimentally measured values. DFT calculations suggested H₂ activation goes through a proton-coupled electron-transfer-like mechanism. During proton transfer to a basic support hydroxyl group, electron density is distributed through the gold nanorod and partially localized on the protonated support hydroxyl group. Water slows H₂ activation by slowing this H⁺ transfer, forcing negative charge buildup on the Au and increasing the transition state energy.

A systematic theoretical study of the water gas shift reaction on the Pt/ZrO2 interface and Pt(111) face: key role of a potassium additive

K can enhance the activity of the WGSR on the Pt₄₀/ZrO₂ model by reducing both the H₂O and COOH dissociation barriers.

A low-temperature water–gas shift reaction catalyzed by hybrid NiO@NiCr-layered double hydroxides: catalytic property, kinetics and mechanism investigation

The realization of a high efficiency water gas shift reaction (WGSR) at low temperatures has always been a research hotspot and is difficult to achieve.

An electronic perturbation in TiC supported platinum monolayer catalyst for enhancing water-gas shift performance: DFT study

The water-gas shift (WGS) reaction behaviors over the TiC(0 0 1) supported Pt monolayer catalyst (Pt_ML/TiC(0 0 1)) are investigated by using the spin-unrestricted density functional theory calculations. Importantly, we find that the Pt_ML/TiC(0 0 1) system exhibits a much lower density of Pt-5d states nearby the Fermi level compared with that for Pt(1 1 1), and the monolayer Pt atoms undergo an electronic perturbation when in contact with TiC(0 0 1) support that would strongly improve the WGS activity of supported Pt atoms. Our calculations clearly indicate that the dominant reaction path follows a carboxyl mechanism involving a key COOH intermediate, rather than the common redox mechanism. Furthermore, through the detailed comparisons, the results demonstrate that the strong interactions between the monolayer Pt atoms and TiC(0 0 1) support make Pt_ML/TiC(0 0 1) a highly active catalyst for the low-temperature WGS reaction. Following the route presented by Bruix et al (2012 J. Am. Chem. Soc. 134 8968-74), the positive effect derived from strong metal-support interaction in the metal/carbide system is revealed.

Vibrational spectra of dissociatively adsorbed D2O on Al-terminated α-Al2O3(0001) surfaces from ab initio molecular dynamics

Water can adsorb molecularly or dissociatively onto different sites of metal oxide surfaces. These adsorption sites can be disentangled using surface-sensitive vibrational spectroscopy. Here, we model Vibrational Sum Frequency (VSF) spectra for various forms of dissociated, deuterated water on a reconstructed, Al-terminated α-Al₂O₃(0001) surface at submonolayer coverages (the so-called 1-2, 1-4, and 1-4' modes). Using an efficient scheme based on velocity-velocity autocorrelation functions, we go beyond previous normal mode analyses by including anharmonicity, mode coupling, and thermal surface motion in the framework of ab initio molecular dynamics. In this way, we calculate vibrational density of states curves, infrared, and VSF spectra. Comparing computed VSF spectra with measured ones, we find that relative frequencies of resonances are in quite good agreement and linewidths are reasonably well represented, while VSF intensities coincide not well. We argue that intensities are sensitively affected by local interactions and thermal fluctuations, even at such low coverage, while absolute peak positions strongly depend on the choice of the electronic structure method and on the appropriate inclusion of anharmonicity.

Abundant Ce3+ Ions in Au-CeOx Nanosheets to Enhance CO2 Electroreduction Performance

The electroreduction of CO₂ to CO provides a potential way to solve the environmental problems caused by excess fossil fuel utilization. Loading transition metals on metal oxides is an efficient strategy for CO₂ electroreduction as well as for reducing metal usage. However, it needs a great potential to overcome the energy barrier to increase CO selectivity. This paper describes how 8.7 wt% gold nanoparticles (NPs) loaded on CeO_x nanosheets (NSs) with high Ce³⁺ concentration effectively decrease the overpotential for CO₂ electroreduction. The 3.6 nm gold NPs on CeO_x NSs containing 47.3% Ce³⁺ achieve CO faradaic efficiency of 90.1% at -0.5 V in 0.1 m KHCO₃ solution. Furthermore, the CO₂ electroreduction activity shows a strong relationship with the fractions of Ce³⁺ on Au-CeO_x NSs, which has never been reported. In situ surface-enhanced infrared absorption spectroscopy shows that Au-CeO_x NSs with high Ce³⁺ concentration promote CO₂ activation and *COOH formation. Theoretical calculations also indicate that the improved performance is attributed to the enhanced *COOH formation on Au-CeO_x NSs with high Ce³⁺ fraction.