Microkinetic Model Research Articles

Rational Design of Cost-Effective Metal-Doped ZrO2 for Oxygen Evolution Reaction

HighlightsSurface energy and surface Pourbaix diagram reveal that ZrO2 (1¯11\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\overline{1 }11$$\\end{document}) is the most thermodynamically stable facet and is preferentially occupied by HO* at the equilibrium potential of oxygen evolution reaction (OER).Microkinetic modeling analyzed the OER activity of 40 single-metal doped ZrO2 and identified 16 metals exhibit improved catalytic activity, with Rh and Fe dopants showing the remarkable improvement.Thermodynamic free energy diagrams, density of states analysis, and ab initio molecular dynamics simulations further confirm that Fe–ZrO2 and Rh–ZrO2 are highly promising catalysts for OER, showcasing low ΔG for the rate-determining step, high conductivity, and exceptional stability.

Microkinetic Modeling to Decode Catalytic Reactions and Empower Catalytic Design

AbstractKinetic model development is integral for designing, redesigning, monitoring, and optimizing chemical processes. Of the various approaches used within this field, microkinetic modeling is a crucial tool that focuses on surface events to analyze overall and preferential reaction pathways. This work covers noticeable features of microkinetic modeling for three critical case studies: (i) ammonia to hydrogen, (ii) oxidative coupling of methane to chemicals, and (iii) carbon dioxide hydrogenation for methanol synthesis. We analyze how microkinetic modeling enables predicting and optimizing complex reaction networks, allowing the design of efficient and tailored catalysts with enhanced activity and selectivity.

The Importance of Reaction Energy in Predicting Chemical Reaction Barriers with Machine Learning Models.

Improving our fundamental understanding of complex heterocatalytic processes increasingly relies on electronic structure simulations and microkinetic models based on calculated energy differences. In particular, calculation of activation barriers, usually achieved through compute-intensive saddle point search routines, remains a serious bottleneck in understanding trends in catalytic activity for highly branched reaction networks. Although the well-known Brønsted-Evans-Polyani (BEP) scaling - a one-feature linear regression model - has been widely applied in such microkinetic models, they still rely on calculated reaction energies and may not generalize beyond a single facet on a single class of materials, e. g., a terrace sites on transition metals. For highly branched and energetically shallow reaction networks, such as electrochemical CO2 reduction or wastewater remediation, calculating even reaction energies on many surfaces can become computationally intractable due to the combinatorial explosion of states that must be considered. Here, we investigate the feasibility of activation barrier prediction without knowledge of the reaction energy using linear and nonlinear machine learning (ML) models trained on a new database of over 500 dehydrogenation activation barriers. We also find that inclusion of the reaction energy significantly improves both classes of ML models, but complex nonlinear models can achieve performance similar to the simplest BEP scaling when predicting activation barriers on new systems. Additionally, inclusion of the reaction energy significantly improves generalizability to new systems beyond the training set. Our results suggest that the reaction energy is a critical feature to consider when building models to predict activation barriers, indicating that efforts to reliably predict reaction energies through, e. g., the Open Catalyst Project and others, will be an important route to effective model development for more complex systems.

Facet-Dependent Oxygen Evolution Reaction Activity of IrO2 from Quantum Mechanics and Experiments.

The diversity of chemical environments present on unique crystallographic facets can drive dramatic differences in catalytic activity and the reaction mechanism. By coupling experimental investigations of five different IrO2 facets and theory, we characterize the detailed elemental steps of the surface redox processes and the rate-limiting processes for the oxygen evolution reaction (OER). The predicted complex evolution of surface adsorbates and the associated charge transfer as a function of applied potential matches well with the distinct redox features observed experimentally for the five facets. Our microkinetic model from grand canonical quantum mechanics (GC-QM) calculations demonstrates mechanistic differences between nucleophilic attack and O-O coupling across facets, providing the rates as a function of applied potential. These GC-QM calculations explain the higher OER activity observed on the (100), (001), and (110) facets and the lower activity observed for the (101) and (111) facets. This combined study with theory and experiment brings new insights into the structural features that either promote or hinder the OER activity of IrO2, which are expected to provide parallels in structural effects on other oxide surfaces.

Role of oxygen vacancies in CO2 methanation over zirconia: A mechanistic DFT and microkinetic study

The use of CO2 as a source of carbon for energy-storage through their transformation into methane represents an attractive approach in modern days. The development of stable and active catalysts as well as a deeper understanding of reaction mechanisms is essential in that field. Herein, the CO2 methanation reaction over zirconia was thoroughly investigated via DFT calculations and microkinetic modeling. The catalytic surface was represented by a monoclinic ZrO2 surface where CO2 adsorption and elementary hydrogenation steps have been systematically explored. It was demonstrated that the participation of an active oxygen vacancy (reduced surface) is crucial to activate the stable bonds of CO2 for further hydrogenation steps; the stoichiometric surface has a minor contribution to CH4 production. Two main reaction pathways were investigated: i) CO2 dissociation with later hydrogenations; (ii) H-assisted (early hydrogenation) mechanisms. It was demonstrated that the latter is kinetically favorable passing through a formate intermediate. Microkinetic simulations indicate that the initial adsorption configuration plays a central role to methane formation: more weakly adsorbed geometry happens to be more important. Thus, the present findings provide new insights and critical information on the role of oxygen vacancies to methanation mechanisms.

Coverage- and Facet-Dependent Multiscale Modeling of O* and H* Adsorption on Pt Catalytic Nanoparticles.

Within computational heterogeneous catalysis, two critical factors exist-coverage and multifaceted effects-which are challenging to incorporate and contribute to differences between the results obtained from computational and experimental studies. Such disparities exist when significant adsorbate-adsorbate interactions are present, particularly when coupled with computationally limited facet sampling. Here, we designed a study to demonstrate the significance of coverage and facet effects on the predicted coverages for O* and H* on Pt nanoparticles. This is accomplished by employing multiscale modeling techniques using a three-pronged approach consisting of density functional theory (DFT), ab initio phase diagrams, and mean-field microkinetic models. Overall, adsorbate-adsorbate interactions are repulsive and far stronger for O*/Pt than for H*/Pt. For O* on Pt(111), repulsive interactions are both two- and three-body, but on Pt(100) and Pt(110), they are predominantly three-body. Through benchmarks to existing experimental literature, we demonstrate that experimentally observed coverages and desorption temperatures can be accurately estimated by computational models when the adsorbate-adsorbate interactions are included. Finally, by combining microkinetic models at the equilibrium limit with kubic harmonic interpolation, we model the impacts of the treatment of adsorbate-adsorbate interactions on the predicted O* and H* coverages over multifaceted Pt nanoparticles. Omitting adsorbate-adsorbate interactions over Pt nanoparticles leads to an overestimation of equilibrium coverages of the adsorbates (maximum of 0.29 and 0.11 ML for O* and H*, respectively) across a wide range of temperatures and pressures relevant to heterogeneous catalysis. Altogether, our work demonstrates that the inclusion of coverage and facet effects increases the accuracy of computational models for heterogeneous catalysts.

A Microkinetic Model for Isomerizing Ethenolysis

A Microkinetic Model for Isomerizing Ethenolysis

Steady-State and transient kinetic investigations of the oxidation of NO over Pt/SiO2

The reaction mechanism of the oxidation of NO under conditions of high NO concentration over a 2 wt% Pt/SiO2 was studied using both kinetic and transient isotopic tracing experiments. The reaction mechanism was found to involve both NO adsorption and oxygen adsorption, contrary to the situation under low NO concentration where NO reacts from the gas phase. The kinetic rate data could be fitted to the Langmuir-Hinshelwood model permitting the estimation of rate and equilibrium constants of elementary steps. These steps were investigated using transient responses following independent reactant gas isotopic switches for tracing the nitrogen path and the oxygen path during NO oxidation. These transients could be sufficiently described by a microkinetic model based on two pools of NO intermediates, one pool each for the adsorbed molecular oxygen and atomically adsorbed oxygen. Surface atomic oxygen was found to be formed by two different routes: via a direct dissociation of molecularly adsorbed oxygen and also via an assisted pathway that leads to NO2 formation. The latter surface reaction between adsorbed NO and adsorbed molecular oxygen comprised the kinetically relevant step. Insights from the two experimental and model approaches were in good agreement and provided a coherent reaction mechanism for NO oxidation under high NO reactant concentration.

Direct HCN synthesis via plasma-assisted conversion of methane and nitrogen.

Hydrogen cyanide (HCN) is synthesized from ammonia (NH3) and methane (CH4) at ~1200°C over a Pt catalyst. Ammonia synthesis entails several complex, highly emitting processes. Plasma-assisted HCN synthesis directly from CH4 and nitrogen (N2) could be pivotal for on-demand HCN production. Here, we evaluate the potential of dielectric barrier discharge (DBD) N2/CH4 plasma for decentralized catalyst-free selective HCN production. We demonstrate a single-step conversion of methane and nitrogen to HCN with a 72% yield at <300°C. HCN is favored at low CH4 concentrations with ethane (C2H6) as the secondary product. We propose a first-principles microkinetic model with few electron impact reactions. The model accurately predicts primary product yields and elucidates that methyl radical (·CH3) is a common intermediate in HCN and C2H6 synthesis. Compared to current industrial processes, N2/CH4 DBD plasma can achieve minimal CO2 emissions.

Theoretical Investigation of the Adsorbate and Potential-Induced Stability of Cu Facets During Electrochemical CO2 and CO Reduction.

The activity and product selectivity of electrocatalysts for reactions like the carbon dioxide reduction reaction (CO2RR) are intimately dependent on the catalyst's structure and composition. While engineering catalytic surfaces can improve performance, discovering the key sets of rational design principles remains challenging due to limitations in modeling catalyst stability under operating conditions. Herein, we perform first-principles density functional calculations adopting implicit solvation methods with potential control to study the influence of adsorbates and applied potential on the stability of different facets of model Cu electrocatalysts. Using coverage dependencies extracted from microkinetic models, we describe an approach for calculating potential and adsorbate-dependent contributions to surface energies under reaction conditions, where Wulff constructions are used to understand the morphological evolution of Cu electrocatalysts under CO2RR conditions. We identify that CO*, a key reaction intermediate, exhibits higher kinetically and thermodynamically accessible coverages on (100) relative to (111) facets, which can translate into an increased relative stabilization of the (100) facet during CO2RR. Our results support the known tendency for increased (111) faceting of Cu nanoparticles under more reducing conditions and that the relative increase in (100) faceting observed under CO2RR conditions is likely attributed to differences in CO* coverage between these facets.

Regulating the Electronic Synergy of Asymmetric Atomic Fe Sites with Adjacent Defects for Boosting Activity and Durability toward Oxygen Reduction

AbstractThe oxygen reduction reaction (ORR) plays a fundamental role in sustainable energy technologies. However, the creation of non‐precious metal electrocatalysts with high ORR activity and durability under all pH conditions is of great significance but remains challenging. Herein, the aim is to overcome this challenge by creating a Fe single atom catalyst on a 2D defect‐containing nitrogen‐doped carbon support (Fe1/DNC) via a microenvironment engineering strategy. Microkinetic modeling reveals that FeN4(OH) moieties are the real active sites under reaction conditions. Due to the synergistic promotion effect of denser accessible FeN4(OH) moieties and defect‐induced electronic properties, Fe1/DNC catalyst achieves extraordinary ORR activity under alkaline, acidic, and neutral conditions, with half‐wave potentials of 0.95, 0.82, and 0.70 V, respectively. Moreover, a negligible performance decay is observed with this Fe catalyst in stability and methanol tolerance tests. Zn‐air battery employing Fe1/DNC delivers remarkable peak power density and long‐term operational durability. Theoretical analysis provides compelling evidence that the defects adjacent to FeN4(OH) moieties can endow an inductive effect to reshape electronic properties to balance the OOH* formation and OH* reduction. This work offers insight into the regulation of asymmetric coordination structure and electronic properties of metal sites for boosting electrocatalytic activity and stability.

Understanding the Reverse Water Gas Shift Reaction over Mo2C MXene Catalyst: A Holistic Computational Analysis

AbstractPristine Mo2C MXene has been proposed as an heterogeneous catalysis of the reverse water gas shift (RWGS) reaction. The present computational study aims at understanding its catalytic performance and reaction mechanisms tackling its thermodynamics, kinetics, and surface dynamic effects, combining Gibbs free energy profiles gained by density functional theory (DFT), mean‐field kinetics by microkinetic modeling, and rare‐event steps by kinetic Monte Carlo (kMC). The RWGS endergonicity goes for the use of high temperatures and reactants partial pressures to make the reaction exergonic. Gibbs free energy profiles show a preference for redox mechanism, whereas microkinetic simulations favor a low‐temperature preference of formate mechanism. The kMC reveals simultaneous operating redox and formate pathways, where surface coverage disfavors redox favoring the formate pathway. A peak performance is found at 700 K, in line with reported experiments, where the formation of surface O2* is found to be key, acting as a reservoir for O* adatoms while freeing surface sites upon O2* formation. Even though high turnover frequencies are predicted, the system could benefit from swing operando conditions, alternating CO production steps with H2 reduction regeneration steps, and/or ways to reduce the surface O2* and so to have more active catalytic sites.

High-Pressure oxidative coupling of methane on alkali metal catalyst – Microkinetic analysis and operando thermal visualization

To introduce promotional H2O effects for both CH4 rate and C2 selectivity, the OH radical formation, catalyzed through H2O activation with O2 surface species, was critical for modeling selective Mn-K2WO4/SiO2 catalysts. Based on our reported experimental evidence, which demonstrates the formation of H2O2 through surface alkali peroxide intermediate, the elementary reactions that account for the OH-mediated pathway were added into the microkinetic model. The advanced model adeptly replicated the promotional H2O effects on both OCM rate and selectivity. The data from a low-pressure microkinetic study were treated isothermally, and extended for near-industrially relevant pressures up to 901 kPa. Thermal visualization using an infrared camera found substantial temperature increases at undiluted high-pressure conditions which caused C2 selectivity to drop significantly. When the furnace temperatures were decreased after ignition, side reactions after O2 depletion (e.g., hydrocarbon reforming) were suppressed, obtaining 13.7 (11.8) % yields at 19.9 % CH4 conversion with 68.6 (59.1) % selectivities for C2-4 (C2) at 901 kPa. The temperature was found to be the determining factor of C2 yield which was perturbed by varying space velocity or CH4/O2 ratios. The optimum temperature for high-pressure conditions was predicted as 885 °C at 901 kPa. The study provides mechanistic and industrially relevant understandings for further OCM catalyst design and system application.

Electrokinetic Analyses Uncover the Rate‐Determining Step of Biomass‐Derived Monosaccharide Electroreduction on Copper

AbstractElectrochemical biomass conversion holds promise to upcycle carbon sources and produce valuable products while reducing greenhouse gas emissions. To this end, deep insight into the interfacial mechanism is essential for the rational design of an efficient electrocatalytic route, which is still an area of active research and development. Herein, we report the reduction of dihydroxyacetone (DHA)—the simplest monosaccharide derived from glycerol feedstock—to acetol, the vital chemical intermediate in industries, with faradaic efficiency of 85±5 % on a polycrystalline Cu electrode. DHA reduction follows preceding dehydration by coordination with the carbonyl and hydroxyl groups and the subsequent hydrogenation. The electrokinetic profile indicates that the rate‐determining step (RDS) includes a proton‐coupled electron transfer (PCET) to the dehydrated intermediate, revealed by coverage‐dependent Tafel slope and isotopic labeling experiments. An approximate zero‐order dependence of H+ suggests that water acts as the proton donor for the interfacial PCET process. Leveraging these insights, we formulate microkinetic models to illustrate its origin that Eley–Rideal (E−R) dominates over Langmuir–Hinshelwood (L−H) in governing Cu‐mediated DHA reduction, offering rational guidance that increasing the concentration of the adsorbed reactant alone would be sufficient to promote the activity in designing practical catalysts.

Interpretation of H2-TPR from Cu-CHA Using First-Principles Calculations.

Temperature-programmed reduction and oxidation are used to obtain information on the presence and abundance of different species in complex catalytic materials. The interpretation of the temperature-programmed reaction profiles is, however, often challenging. One example is H2 temperature-programmed reduction (H2-TPR) of Cu-chabazite (Cu-CHA), which is a material used for ammonia assisted selective catalytic reduction of NOx (NH3-SCR). The TPR profiles of Cu-CHA consist generally of three main peaks. A peak at 220 °C is commonly assigned to ZCuOH, whereas peaks at 360 and 500 °C generally are assigned to Z2Cu, where Z represents an Al site. Here, we analyze H2-TPR over Cu-CHA by density functional theory calculations, microkinetic modeling, and TPR measurements of samples pretreated to have a dominant Cu species. We find that H2 can react with Cu ions in oxidation state +2, whereas adsorption on Cu ions in +1 is endothermic. Kinetic modeling of the TPR profiles suggests that the 220 °C peak can be assigned to Z2CuOCu and ZCuOH, whereas the peaks at higher temperatures can be assigned to paired Z2Cu and Z2CuHOOHCu species (360 °C) or paired Z2Cu and Z2CuOOCu (500 °C). The results are in good agreement with the experiments and facilitate the interpretation of future TPR experiments.

Electrokinetic Analyses Uncover the Rate-Determining Step of Biomass-Derived Monosaccharide Electroreduction on Copper.

Electrochemical biomass conversion holds promise to upcycle carbon sources and produce valuable products while reducing greenhouse gas emissions. To this end, deep insight into the interfacial mechanism is essential for the rational design of an efficient electrocatalytic route, which is still an area of active research and development. Herein, we report the reduction of dihydroxyacetone (DHA)-the simplest monosaccharide derived from glycerol feedstock-to acetol, the vital chemical intermediate in industries, with faradaic efficiency of 85±5 % on a polycrystalline Cu electrode. DHA reduction follows preceding dehydration by coordination with the carbonyl and hydroxyl groups and the subsequent hydrogenation. The electrokinetic profile indicates that the rate-determining step (RDS) includes a proton-coupled electron transfer (PCET) to the dehydrated intermediate, revealed by coverage-dependent Tafel slope and isotopic labeling experiments. An approximate zero-order dependence of H+ suggests that water acts as the proton donor for the interfacial PCET process. Leveraging these insights, we formulate microkinetic models to illustrate its origin that Eley-Rideal (E-R) dominates over Langmuir-Hinshelwood (L-H) in governing Cu-mediated DHA reduction, offering rational guidance that increasing the concentration of the adsorbed reactant alone would be sufficient to promote the activity in designing practical catalysts.

A Brief Review of Recent Theoretical Advances in Fe-Based Catalysts for CO2 Hydrogenation.

Catalytic hydrogenation presents a promising approach for converting CO2 into valuable chemicals and fuels, crucial for climate change mitigation. Iron-based catalysts have emerged as key contributors, particularly in driving the reverse water-gas shift and Fischer-Tropsch synthesis reactions. Recent research has focused on enhancing the efficiency and selectivity of these catalysts by incorporating alkali metal promoters or transition metal dopants, enabling precise adjustments to their composition and properties. This review synthesizes recent theoretical advancements in CO2 hydrogenation with iron-based catalysts, employing density functional theory and microkinetic modeling. By elucidating the underlying mechanisms involving metallic iron, iron oxides, and iron carbides, we address current challenges and provide insights for future sustainable CO2 hydrogenation developments.

Ni-based bimetallic catalysts for hydrogen production via (sorption-enhanced) steam methane reforming

The catalytic performance of a monometallic Ni/Al2O3 and three bimetallic catalysts (Ni3M1/Al2O3, with M = Cu, Fe, and Ge) for the (sorption-enhanced) steam methane reforming reaction was evaluated. Ni3Cu1/Al2O3 was found to be the optimal catalyst in terms of methane conversion, hydrogen yield, and purity. Ge also has a promoting effect on the monometallic Ni catalyst, whereas the addition of Fe negatively influenced its performance. Physico-chemical characterization of the materials indicated the formation of alloys upon activation of the materials with hydrogen. The addition of Cu increased the surface area and metal dispersion, and improved the overall morphology of the catalyst. The experimental observations were also supported by a numerical study combining Density Functional Theory-based calculations and Microkinetic modelling of the SMR process. Ni3Cu1 and Ni3Ge1 were calculated to have a similar level of catalytic activity as Ni, whereas Ni3Fe1 was unsuitable for the reaction. The SMR reaction was further improved by adding calcium oxide as the CO2 sorbent, which increased methane conversion, CO selectivity, hydrogen yield, and hydrogen purity. The highest methane conversion of 97 % was achieved by Ni/Al2O3 and Ni3Cu1/Al2O3 at 700 °C.

Reaction of (N4Py)Fe with H2O2 and the relevance of its Fe(IV)=O species during and after H2O2 disproportionation

AbstractThe catalytic disproportionation of by non‐heme Fe(II) complexes of H2O2 the ligand N4Py (1,1‐bis(pyridin‐2‐yl)‐N,N‐bis(pyridin‐2‐ylmethyl)methanamine) and the formation and reactivity of Fe(III)‐OOH and Fe(IV)=O species is studied by UV/Vis absorption, NIR luminescence, (resonance) Raman and headspace Raman spectroscopy, 1O2 trapping and DFT methods. Earlier DFT studies indicated that disproportionation of H2O2 catalysed by Fe(II)‐N4Py complexes produce only 3O2, however, only the low‐spin state pathway was considered. In the present study, DFT calculations predict two pathways for the reaction between Fe(III)‐OOH and H2O2, both of which yield 3O2/H2O2 and involve either the S=1/2 or the S=3/2 spin state, with the latter being spin forbidden. The driving force for both pathways are similar, however, a minimal energy crossing point (MECP) provides a route for the formally spin forbidden reaction. The energy gap between the reaction intermediate and the MECP is lower than the barrier across the non‐adiabatic channel. The formation of 3O2 only is confirmed experimentally in the present study through 1O2 trapping and NIR luminescence spectroscopy. However, attempts to use the 1O2 probe ( ‐terpinene) resulted in initiation of auto‐oxidation rather than formation of the expected endoperoxide, which indicated formation of OH radicals from Fe(III)‐OOH, e. g., through O−O bond homolysis together with saturation of methanol with 3O2. Microkinetic modelling of spectroscopic data using rate constants determined earlier, reveal that there is another pathway for Fe(III)‐OOH decomposition in addition to competition between the reaction of Fe(III)‐OOH with H2O2 and homolysis to form Fe(IV)=O and hydroxyl radical. Notably, after all H2O2 is consumed the decay of the Fe(III)‐OOH species is predominantly through a second order self reaction (with Fe(III)‐OOH). The conclusion reached is that the rate of O−O bond homolysis in the Fe(III)‐OOH species to form Fe(IV)=O and an hydroxyl radical is too low to be responsible for the observed oxidation of organic substrates.

Unravelling the reaction mechanism for H2 production via formic acid decomposition over Pd: Coverage-Dependent microkinetic modeling

As a key liquid organic hydrogen carrier, investigating the decomposition of formic acid (HCOOH) on the Pd (111) transition metal surface is imperative for harnessing hydrogen energy. Despite a multitude of studies, the major mechanisms and key intermediates involved in the dehydrogenation process of formic acid remain a great topic of debate due to ambiguous adsorbate interactions. In this research, we develop an advanced microkinetic model based on first-principles calculations, accounting for adsorbate–adsorbate interactions. Our study unveils a comprehensive mechanism for the Pd (111) surface, highlighting the significance of coverage effects in formic acid dehydrogenation. Our findings unequivocally demonstrate that H coverage on the Pd (111) surface renders formic acid more susceptible to decompose into H2 and CO2 through COOH intermediates. Consistent with experimental results, the selectivity of H2 in the decomposition of formic acid on the Pd (111) surface approaches 100 %. Considering the influence of H coverage, our kinetic analysis aligns perfectly with experimental values at a temperature of 373 K.