SEARCH · Search NASA
Results for “Microkinetic Modeling”
Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
From microkinetic model to process: understanding the role of the boron nitride surface and gas phase chemistry in the oxidative dehydrogenation of propane
In the oxidative dehydrogenation of propane over boron nitride, propylene selectivity is lost due to unselective gas-phase chemistry.
Mapping the Degree of Rate Control Using Automated Construction of Microkinetic Models with Rmg-Cat.
Abstract not provided.
Automated Construction of Microkinetic Models with RMG-Cat for Mapping the Degree of Rate Control.
Abstract not provided.
Mechanisms of Ethylene Epoxidation over Silver from Machine Learning-Accelerated First-Principles Modeling and Microkinetic Simulations
In this work we employed machine learning-augmented density functional theory (DFT) thermodynamic calculations to assess the stability of different AgO x structures under catalytic ethylene epoxidation reaction conditions. We found that there are multiple AgO x surface motifs that could co-exist under the relevant conditions. These included Ag surface oxides (e.g., AgO_p(4 × 4) and Ag 1.83 O) and atomic oxygen-covered Ag(111) surfaces. Furthermore, we employed DFT calculations to evaluate the energetics of different reaction mechanisms by which ethylene and oxygen can react on these surfaces. These studies revealed several energetically viable reaction pathways for ethylene epoxidation. Microkinetic modeling analyses, based on the DFT-calculated reaction pathways, showed that ethylene epoxidation can proceed on all surfaces and that multiple pathways, including those involving Langmuir–Hinshelwood and Eley–Rideal mechanisms, could be involved in selective and unselective reactions. The diversity of mechanisms that we discovered in the context of the relatively simple ethylene epoxidation reaction on Ag suggests that the richness and complexity of surface chemistry are most likely a rule rather than an exception in heterogeneous catalytic chemical transformations on metal surfaces and that the concept of a single or even a dominant mechanism and reaction intermediates might need to be revisited for many reactions.
Pore-Scale Transport Effects in Electrochemical CO 2 Reduction on Gold via Coupled Microkinetic-Transport Modeling
A pore-resolved modeling framework is developed to quantify how pore-scale transport affects the intrinsic microkinetics of CO 2 -to-CO on Au. A DFT-informed microkinetic model is coupled self-consistently to a Generalized-Modified Poisson–Nernst–Planck (GMPNP) transport description in a single, electrolyte-filled cylindrical pore, allowing local concentrations and potential to feed back into site-specific reaction rates. FIB-SEM is used to determine pore sizes within realistic electrode materials. Across pore diameters, d p = 10–6000 nm, the surface-averaged CO 2 reduction rate is systematically reduced relative to the ideal microkinetic baseline where mass transport is not accounted for; the effectiveness factor 𝜂 𝑠,CO 2 , which quantifies this ratio, decreases rapidly at more negative potentials and is about 1% near −1.0 V vs SHE due to reactant depletion. Spatial maps reveal pore-bulk alkalization that emerges at higher cathodic bias, with a small, near-wall pH dip due to electrostatic repulsion of hydroxide at the cathode interface. For a fixed aspect ratio L p /d p , narrower pores exhibit larger 𝜂 𝑠,CO 2 by shortening diffusion paths, whereas variations in the aspect ratio L p /d p play a secondary role. A dimensionless analysis (surface/bulk Damköhler numbers) delineates operating regimes. In conclusion, this work offers a concept for incorporating microkinetic models into homogenized porous-electrode models through effectiveness factors and pore-size distribution.
A Theoretical Study on Reaction Mechanisms and Kinetics of Thiophene Hydrodesulfurization over MoS 2 Catalysts
In the present work, thiophene hydrodesulphurization (HDS) over the Mo-edge, the S-edge, and the Mo-S connection edge of MoS 2 catalyst with 50% sulfur coverage was studied using first-principles based microkinetic modeling. Two parallel HDS routes, i.e., direct desulfurization (DDS) and hydrogenation (HYD) were taken into account. It has been found that the major reaction route of thiophene HDS on the Mo- and the Mo-S edges is temperature dependent. In the low temperature range of 500~600 K, the HYD route is dominant, leading to the C 4 H 8 formation. As the temperature increases, the DDS route becomes competitive with the HYD route. At the temperature above 650 K, the DDS route will be the dominant HDS reaction route on the Mo- and the Mo-S edges. While the DDS route that leads to the formation of C4H6 is the major thiophene HDS reaction route on the S-edge in the entire temperature range of 500~750 K. The microkinetic modeling results show the overall HDS activity on the S-edge is lower than it on the Mo- and the Mo-S edges. The Mo-S edge also provides a preferential reaction pathway, which facilitates 2-hydrothiophene migration from the Mo-edge to the S-edge, followed by remaining elementary steps with lower activation barriers in the DDS route. The work was financially supported by the National Natural Science Foundation of China (No. 21476012, 21571012, 21407007 and 91534201). The computing time was granted by a scientific theme user proposal in the William R. Wiley Environmental Molecular Sciences Laboratory (EMSL), which is a U.S. Department of Energy national scientific user facility located at Pacific Northwest National Laboratory (PNNL) in Richland, Washington. Part of computing time was provided the National Energy Research Scientific Computing Center (NERSC). D. Mei is supported by the US Department of Energy, Office of Sciences, Office of Basic Energy Sciences, Division of Chemical Sciences, Geosciences & Biosciences.
Microkinetic insights into the role of catalyst and water activity on the nucleation, growth, and dissolution during COF-5 synthesis
The chemical pathway for synthesizing covalent organic frameworks (COFs) involves a complex medley of reaction sequences over a rippling energy landscape that cannot be adequately described using existing theories. Even with the development of state-of-the-art experimental and computational tools, identifying primary mechanisms of nucleation and growth of COFs remains elusive. Other than empirically, little is known about how the catalyst composition and water activity affect the kinetics of the reaction pathway. In this work, for the first time, we employ time-resolved in situ Fourier transform infrared spectroscopy (FT-IR) coupled with a six-parameter microkinetic model consisting of ~10 million reactions and over 20 000 species. The integrated approach elucidates previously unrecognized roles of catalyst pK a on COF yield and water on growth rate and size distribution. COF crystalline yield increases with decreasing pK a of the catalysts, whereas the effect of water is to reduce the growth rate of COF and broaden the size distribution. The microkinetic model reproduces the experimental data and quantitatively predicts the role of synthesis conditions such as temperature, catalyst, and precursor concentration on the nucleation and growth rates. Furthermore, the model also validates the second-order reaction mechanism of COF-5 and predicts the activation barriers for classical and non-classical growth of COF-5 crystals. The microkinetic model developed here is generalizable to different COFs and other multicomponent systems.
Unraveling the CO Oxidation Mechanism over Highly Dispersed Pt Single Atom on Anatase TiO 2 (101)
Catalysts with noble metals deposited as single atoms on metal oxide supports have recently been studied extensively due to their maximized metal utilization and potential for performing difficult chemical conversions owing to their unique electronic properties. Understanding of the reaction mechanisms on supported single-metal atoms is still limited but is highly important for designing more efficient catalysts. In this study, we report the complexity of the CO oxidation reaction mechanism on Pt single atoms supported on anatase TiO 2 (Pt SA /a-TiO 2 ) by coupling density functional theory (DFT) calculations and microkinetic analysis with kinetic measurements, in situ/operando infrared, and X-ray absorption spectroscopies. Starting from the adsorbed Pt SA occupying an O vacancy induced by reductive pretreatment, we show that CO oxidation follows a complex mechanism consisting of initiation steps to reorganize the active site and multibranch reactive cycles, with the Pt SA /a-TiO 2 catalyst not returning to its initial configuration. The initiation step consists of CO and O 2 adsorption healing the O vacancy, followed by CO oxidation using gas-phase CO to form Pt(CO). The reactive cycle alternates O 2 adsorption and dissociation to oxidize the catalyst to Pt(O)(O)(CO) and branching pathways of competing Langmuir–Hinshelwood (LH)- or Eley–Rideal (ER)-type CO oxidation steps to reduce it again to Pt(CO). In situ/operando infrared experiments, including cryogenic CO adsorption and isotopic CO exchange, confirm the combined involvement of strongly adsorbed CO and gas-phase CO in an Eley–Rideal step along the reaction cycle. Microkinetic modeling shows that Pt single atoms are present in a mixture of Pt(CO), Pt(CO)(O 2 ), Pt(O)(CO)(O 2 ), and Pt(CO)(CO 3 ) structures as the main intermediates during steady-state CO oxidation, all having the C–O vibrational stretch close to the experimentally observed value of 2115 cm –1 . Microkinetic modeling also shows that the fractional orders of CO and O 2 measured experimentally originate from multiple steps with a high degree of rate control and not from a simple competitive adsorption. The results demonstrate the complex reaction pathways that even CO oxidation on a simple single-atom system can follow, providing mechanistic insights for designing efficient Pt-based single-atom catalysts. Here, we further show that microkinetic modeling results are sensitive to changes in energies of intermediate and transition states within errors of density functional theory, which can ultimately lead to incorrect conclusions regarding the reaction pathways and most abundant reaction intermediates if not accounted for by experiments.
A self-adjusting platinum surface for acetone hydrogenation
We show that platinum displays a self-adjusting surface that is active for the hydrogenation of acetone over a wide range of reaction conditions. Reaction kinetics measurements under steady-state and transient conditions at temperatures near 350 K, electronic structure calculations employing density-functional theory, and microkinetic modeling were employed to study this behavior over supported platinum catalysts. The importance of surface coverage effects was highlighted by evaluating the transient response of isopropanol formation following either removal of the reactant ketone from the feed, or its substitution with a similarly structured species. The extent to which adsorbed intermediates that lead to the formation of isopropanol were removed from the catalytic surface was observed to be higher following ketone substitution in comparison to its removal, indicating that surface species leading to isopropanol become more strongly adsorbed on the surface as the coverage decreases during the desorption experiment. This phenomenon occurs as a result of adsorbate–adsorbate repulsive interactions on the catalyst surface which adjust with respect to the reaction conditions. Reaction kinetics parameters obtained experimentally were in agreement with those predicted by microkinetic modeling when the binding energies, activation energies, and entropies of adsorbed species and transition states were expressed as a function of surface coverage of the most abundant surface intermediate (MASI, C 3 H 6 OH*). It is important that these effects of surface coverage be incorporated dynamically in the microkinetic model (e.g., using the Bragg–Williams approximation) to describe the experimental data over a wide range of acetone partial pressures.
Enhanced Selectivity for C 2 H 4 Production from C 2 H 6 on Partially Chlorinated IrO 2 (110) Surfaces
Modifying metal oxide surfaces to limit their oxidizing activity can provide a means of improving catalytic selectivity toward the partial oxidation of light alkanes. Here, in this study, we investigated the oxidation of C 2 H 6 on Cl-modified IrO 2 (110) surfaces using temperature-programmed reaction spectroscopy (TPRS) and first-principles microkinetic modeling. We find that substituting Cl for O in the IrO 2 (110) surface enhances the selectivity for C 2 H 6 conversion to C 2 H 4 during TPRS by suppressing extensive oxidation to CO x products, while also either enhancing C 2 H 4 production or altering it to a lesser extent, depending on the initial C 2 H 6 coverage. The C 2 H 4 selectivity increased with increasing C 2 H 6 and Cl coverage, but reached a limiting value below 50%. The Cl coverage changed negligibly during C 2 H 6 oxidation, and the surface reactivity decreased only marginally for Cl coverages up to 0.5 ML (monolayer). TPRS simulations using a microkinetic model predict C 2 H 4 and CO x product yields as a function of the Cl coverage that agree closely with the experimental results. According to the simulations, C 2 H 6 conversion to C 2 H 4 occurs on Cl-IrO 2 (110) by the hydrogenation of C 2 H 3 * species adsorbed in blocked states, in which neighboring sites are occupied only by unreactive HO and Cl species. The microkinetic modeling shows that H-hopping away from surface HO groups provides a relatively efficient route for C 2 H 3 * to escape blocked configurations and dehydrogenate, and that this process can limit the C 2 H 4 selectivity on Cl-IrO 2 (110) under the conditions studied. Overall, our results demonstrate that Cl-substitution into IrO 2 (110) enhances the selectivity for C 2 H 4 production from C 2 H 6 and provides insights into the reaction mechanism that can guide strategies to further improve the C 2 H 4 selectivity.
Automated Generation of Microkinetics for Heterogeneously Catalyzed Reactions Considering Correlated Uncertainties**
Abstract 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.
Resolving the Coverage Dependence of Surface Reaction Kinetics with Machine Learning and Automated Quantum Chemistry Workflows
Microkinetic models for catalytic systems require estimation of many thermodynamic and kinetic parameters that can be calculated for isolated species and transition states using ab initio methods. However, the presence of nearby coadsorbates on the surface can dramatically alter these thermodynamic and kinetic parameters causing them to be dependent on species coverage fractions. As there are combinatorially many coadsorbed configurations on the surface, computing the coverage dependence of these parameters is far less straightforward. We present a framework for generating and applying machine learning models to predict coverage-dependent parameters for microkinetic models. Our toolkit enables automatic calculation and evaluation of coadsorbed configurations allowing us to sample 2,000 coadsorbed adsorbates and transition states (TSs) for a diverse set of 9 reactions on Cu(111), a challenging surface, with four possible coadsorbates. This dataset was then used to train subgraph isomorphic decision trees (SIDTs) to predict the stability and association energy of configurations. We were able to achieve mean absolute errors (MAEs) of 0.106 eV on adsorbates, 0.172 eV on TSs, and due to natural error cancellation in SIDTs for relative properties, 0.130 eV on reaction energies and 0.180 eV on activation barriers. In conclusion, we describe how to use these models to predict coverage-dependent corrections for adsorbates and TSs and demonstrate on H*, HO*, and O* comparing the generated SIDT model with an iteratively refined version.
Investigation of solvent effects in the hydrodeoxygenation of levulinic acid to γ-valerolactone over Ru catalysts
Liquid phase, reductive deoxygenation of biomass derived platform chemicals over transition metal surfaces constitutes an efficient scheme for upgrading lignocellulosic biomass. The solvation effects on the reaction kinetics of the hydrodeoxygenation (HDO) of levulinic acid (LA) towards the formation of γ-valerolactone (GVL) over Ru(0001) has been studied in three condensed phase media, i.e., liquid water, methanol, and 1,4-dioxane. Detailed microkinetic models have been developed incorporating various catalytic pathways including formation of 4-hydroxypentanoic acid (HPA) and α-angelicalactone (AGL) to simulate the catalytic activity of Ru(0001) under various reaction conditions of solvent, temperature, and partial pressures. Our microkinetic models suggest that direct catalytic conversion with alkoxy formation is the preferred reaction mechanism in all reaction environments. Here, we find that water facilitates the reaction kinetics significantly and that the solvent effect is strongest at lower temperatures (T<373 K). Here, rate increases due to liquid water solvation effects of 2-4 orders of magnitude are observed. All solvents increase the rate of reaction relative to the gas phase; however, solvation effects decrease with decrease in polarity. 1,4-dioxane increases the rate only minimally due to competitive adsorption of the solvent molecules despite facilitating the partially rate controlling step of the LA hydrogenation to an alkoxy intermediate.
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. Here, 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.
Combining Computational Modeling with Reaction Kinetics Experiments for Elucidating the In Situ Nature of the Active Site in Catalysis
Microkinetic modeling based on density functional theory (DFT) derived energetics is important for addressing fundamental questions in catalysis. The quantitative fidelity of microkinetic models (MKMs), however, is often insufficient to conclusively infer the mechanistic details of a specific catalytic system. Here, this can be attributed to a number of factors such as an incorrect model of the active site for which DFT calculations are performed, deficiencies in the hypothesized reaction mechanism, inadequate consideration of the surface environment under reaction conditions, and intrinsic errors in the DFT exchange-correlation functional. Despite these limitations, we aim at developing a rigorous understanding of the reaction mechanism and of the nature of the active site for heterogeneous catalytic chemistries under reaction conditions. By achieving parity between experimental and modeling outcomes through robust parameter estimation and by ensuring coverage-consistency between DFT calculations and MKM predictions, it is possible to systematically refine the mechanistic model and, thereby, our understanding of the catalytic active site in situ.
Pynta-An Automated Workflow for Calculation of Surface and Gas-Surface Kinetics
Many important industrial processes rely on heterogeneous catalytic systems. However, given all possible catalysts and conditions of interest, it is impractical to optimize most systems experimentally. Automatically generated microkinetic models can be used to efficiently consider many catalysts and conditions. However, these microkinetic models require accurate estimation of many thermochemical and kinetic parameters. Manually calculating these parameters is tedious and error prone, involving many interconnected computations. Here, we present Pynta, a workflow software for automating the calculation of surface and gas–surface reactions. Pynta takes the reactants, products, and atom maps for the reactions of interest, generates sets of initial guesses for all species and saddle points, runs all optimizations, frequency, and IRC calculations, and computes the associated thermochemistry and rate coefficients. It is able to consider all unique adsorption configurations for both adsorbates and saddle points, allowing it to handle high index surfaces and bidentate species. Pynta implements a new saddle point guess generation method called harmonically forced saddle point searching (HFSP). HFSP defines harmonic potentials based on the optimized adsorbate geometries and which bonds are breaking and forming that allow initial placements to be optimized using the GFN1-xTB semiempirical method to create reliable saddle point guesses. This method is reaction class agnostic and fast, allowing Pynta to consider all possible adsorbate site placements efficiently. We demonstrate Pynta on 11 diverse reactions involving monodenate, bidentate, and gas-phase species, many distinct reaction classes, and both a low and a high index facet of Cu. Our results suggest that it is very important to consider reactions between adsorbates adsorbed in all unique configurations for interadsorbate group transfers and reactions on high index surfaces.