Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “transition state”

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.

At least 397 records · Page 22
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

An automated cluster surface scanning method for exploring reaction paths on metal-cluster surfaces

Metal-cluster surfaces present a wide variety of unique coordination environments. This complexity makes it difficult to manually probe the surface reactivity of such clusters. Here, we present a simple and automated method to systematically discover reaction pathways on cluster surfaces, based on the automated cluster surface scanning (ACSS) technique for mapping out potential energy surfaces. We showcase our method on 55-atom icosahedral Cu and Ag clusters, where we determine the activation energies of four elementary steps common in heterogeneous catalysis – hydrogen recombination (H* + H* → H 2 * + *), oxygen recombination (O* + O* → O 2 * + *), water formation (OH* + H* → H 2 O* + *), and CO oxidation (CO* + O* → CO 2 * + *) – with density functional theory calculations (DFT-PBE + D3). We show that the ACSS method requires significantly less human effort than the established manually performed climbing-image nudged elastic band (MP + CI-NEB) technique and locates transition states with comparable accuracy (root-mean-squared error of 0.10 eV) and similar computational cost. Rigorous sampling of the potential energy surface with the ACSS method allows one to locate all lowest-energy reaction pathways obtained via the MP+CI-NEB approach, as well as alternative pathways that one may have missed with the MP+CI-NEB approach due to the many possible pathways available on these clusters. The accuracy and efficiency afforded by the ACSS method could enable high-throughput exploration of the diverse reactivity of metal clusters.

36 MATERIALS SCIENCE↗
DOE OSTI2023

QMMM 2023: A program for combined quantum mechanical and molecular mechanical modeling and simulations

Combined quantum mechanical and molecular mechanical (QM/MM) methods play an important role in multiscale modeling and simulations. QMMM 2023 is a general-purpose program for single-point calculations, geometry optimizations, transition-state optimizations, and molecular dynamics (MD) at the QM/MM level. It calls a QM package and an MM package to perform the required single-level calculations and combines them into a QM/MM energy by a variety of schemes. QMMM 2023 supports GAMESS-US, Gaussian, and ORCA as QM packages and Tinker as the MM package. Four types of treatments are available for embedding the QM subsystem in the MM environment: mechanical embedding with gas-phase calculations of the QM region, electronic embedding that allows polarization of the QM region by the MM environment, polarizable embedding for mutual polarization of the QM and MM regions, and flexible embedding for both mutual polarization and partial charge transfer between the QM and MM regions. Boundaries between QM and MM regions that pass through covalent bonds can be treated by several methods, including the redistributed charge (RC) scheme, redistributed charge and dipole (RCD) scheme, balanced-RC scheme, balanced-RCD scheme, screened charge scheme that takes account of charge penetration effects, and smeared charge scheme that delocalizes the MM charges near the QM–MM boundary. Geometry optimization can be done using the optimizer implemented in QMMM 2023 or the Berny optimizer in Gaussian through external calls to Gaussian. Molecular dynamics simulations can be performed at the pure-MM level, pure-QM level, fixed-partitioning QM/MM level, and adaptive-partitioning QM/MM level. As a result, the adaptive-partitioning treatments permit on-the-fly relocation of the QM–MM boundary by dynamically reclassifying atoms or groups into the QM or MM subsystems.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Preparation of a uranium monocarbide anode and electrochemical characterization in molten LiCl-KCl-UCl 3

Porous uranium carbide (UC) pellets possessing moderate electrical conductivity were synthesized by reaction of UO 2 with graphite at temperatures up to 1550°C under rough vacuum. Conversions as high as 98% were achieved at soak times of 2-4 hours. The electrochemistry of the UC pellets in molten LiCl-KCl-6.5 wt% UCl 3 was explored using a variety of techniques including DC polarization methods, cyclic voltammetry, chronopotentiometry and bulk electrolysis. Here, the electrode reaction for anodic dissolution was found to be kinetically controlled by dissociation of UC to a transition state complex that was hypothesized to consist of a uranium atom partially complexed by chloride ions. Precise measurements of current efficiencies using chronopotentiometry indicated upper limits of 90.9 ± 3.4% and 98.3 +1.7/-3.7% for anode and cathode, respectively, when operating at anodic overpotentials near +300 mV. Bulk electrolysis of a UC pellet performed by passing 98% of the theoretical charge resulted in nearly complete recovery of its uranium content as highly pure metal at the cathode.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Fission yeast Duf89 and Duf8901 are cobalt/nickel-dependent phosphatase–pyrophosphatases that act via a covalent aspartyl–phosphate intermediate

Domain of Unknown Function 89 (DUF89) proteins are metal-dependent phosphohydrolases. Exemplary DUF89 enzymes differ in their metal and phosphosubstrate preferences. Here, we interrogated the activities and structures of two DUF89 paralogs from fission yeast—Duf89 and Duf8901. We find that Duf89 and Duf8901 are cobalt/nickel-dependent phosphohydrolases adept at hydrolyzing p-nitrophenylphosphate and PP i . Crystal structures of metal-free Duf89 and Co 2+ -bound Duf8901 disclosed two enzyme conformations that differed with respect to the position of a three-helix module, which is either oriented away from the active site in Duf89 or forms a lid over the active site in Duf8901. Lid closure results in a 16 Å movement of Duf8901 Asp195, vis-à-vis Asp199 in Duf89, that brings Asp195 into contact with an octahedrally coordinated cobalt. Reaction of Duf8901 with BeCl 2 and NaF in the presence of divalent cations Co 2+ , Ni 2+ , or Zn 2+ generated covalent Duf8901-(Asp248)–beryllium trifluoride (BeF 3 )•Co 2+ , Duf8901-(Asp248)–BeF 3 •Ni 2+ , or Duf8901-(Asp248)–BeF 3 •Zn 2+ adducts, the structures of which suggest a two-step catalytic mechanism via formation and hydrolysis of an enzyme-(aspartyl)–phosphate intermediate. Alanine mutations of Duf8901 Asp248, Asn249, Lys401, Asp286, and Asp195 that interact with BeF 3 •Co 2+ squelched p-nitrophenylphosphatase activity. A 1.8 Å structure of a Duf8901-(Asp248)–AlF 4 –OH 2 •Co 2+ transition-state mimetic suggests an associative mechanism in which Asp195 and Asp363 orient and activate the water nucleophile. Whereas deletion of the duf89 gene elicited a phenotype in which expression of phosphate homeostasis gene pho1 was derepressed, deleting duf8901 did not, thereby hinting that the DUF89 paralogs have distinct functional repertoires in vivo.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Acetone to isobutene conversion on Zn x Ti y O z : Effects of TiO 2 facet

In this study, liquid-phase chemical grafting method was used to graft Zn onto TiO 2 with preferentially exposed (1 0 1) or (0 0 1) facet. The obtained Zn x Ti y O z materials were characterized using various techniques (e.g., XRD, Raman, DRIFTS etc.) and evaluated for the acetone-to-isobutene reaction. It was found that over TiO 2 (0 0 1), both terminal and bridging hydroxyls were readily titrated by Zn deposition, whereas a substantial amount of bridging hydroxyls on TiO 2 (1 0 1) remained. Although dominant Zn-O-terminated surface was obtained on two Zn x Ti y O z samples, bridging hydroxyls with high H-D exchange reactivity were observed on Zn x Ti y O z (1 0 1) compared with Zn x Ti y O z (0 0 1). The bridging hydroxyls showing rapid proton transfer efficiently stabilizes a transition state of diacetone alcohol intramolecular rearrangement for isobutene production as opposed to the diacetone alcohol dehydration.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Understanding the influence of solvents on the Pt-catalyzed hydrodeoxygenation of guaiacol

Bio-oil derived from fast pyrolysis of lignocellulosic biomass needs to be deoxygenated to become a substitute for petroleum fuels. Here, we study the hydrodeoxygenation mechanism of guaiacol, a bio-oil model compound derived from the lignin fraction of biomass, on Pt(1 1 1) terrace sites in the presence of water, diethyl ether, 1-butanol, and n-hexane as solvent. Using first-principles periodic density functional theory (DFT) calculations and mean-field microkinetic reactor modeling, a detailed reaction mechanism is investigated targeting various products such as catechol, phenol, anisole, benzene, cyclohexanone, and cyclohexanol. Solvent phase DFT outcomes are mostly similar to that of the vapor phase; however, microkinetic modeling results suggest that rate controlling species and transition states differ somewhat in the various reaction environments. Catechol was found to be the major aromatic product across all reaction environments. Over Pt(1 1 1), unsaturated monooxygenate production from catechol is unlikely.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

A constitutive framework for rocks undergoing solid dissolution

Here, we formulate a time-dependent damage theory for rocks subjected to mechanical deformation and solid dissolution. The constitutive description is inspired by the transition state theory, which states that the rate of dissolution is a function of the reactive surface area measured through the crack density in the volume. We use a gradient-enhanced damage framework in which damage depends on the deformation of the material as well as on the amount of solid mass dissolved over time. The gradient-enhanced formulation is characterized by a three-field variational formulation with the solid displacement, nonlocal equivalent strain, and nonlocal rate of solid dissolution as the basic state variables. Traditionally, time-independent damage theories have only allowed damage to increase with increasing external load. In the proposed framework, the degree of damage may increase due to solid dissolution even when the external load is held fixed. In this way, solid dissolution is viewed as a process that is responsible for bringing about rate-dependent effects such as creep and stress-relaxation, which are two common features of geomaterial behavior.

42 ENGINEERING↗
DOE OSTI2025

Flexibility of oxygen sublattice and hydrogen bond length predict proton mobility in ternary metal oxides

Discovery of fast proton conductors is important for advancing clean energy technologies. This requires a better understanding of proton migration mechanisms. While structural and chemical traits of ternary metal oxides have been related to proton migration barriers, lattice dynamical effects have not been resolved quantitatively. Here, in this work, we introduce a phonon-based dynamic descriptor, termed ‘‘thermal O…O fluctuation,’’ quantifying the flexibility of donor-acceptor oxide-ion pairs. This enables direct comparison of O-sublattice flexibility across diverse metal oxides. Using regression models, we ranked physical descriptors as predictors of proton mobility, finding that H-bond length and thermal O…O fluctuation were the strongest descriptors. Further analysis revealed a critical O…O spacing of 2.4 A˚ at the transition state, which is easier to reach by more flexible donor-acceptor pairs, enabling facile proton transfer. Our results demonstrate oxygen sublattice flexibility as a dynamic descriptor and provide guiding principles for enhancing proton mobility in ternary metal oxides.

diffusivity↗
DOE OSTI2026

Synthesis and application of thermally responsive nanofiber coatings for overtemperature monitoring

This study presents a one-pot synthesis route to organometallic nanofibers based on copper thiolate, exhibiting distinctive chemical and physical characteristics. Electron microscopy analysis of morphology and composition revealed 2-10 μm-long, 50-90 nm-diameter hollow and non-hollow fibers composed of copper, sulfur, oxygen, hydrocarbon, and chlorine. Thermogravimetric analysis showed a pronounced mass loss within 120°C-135°C. To elucidate the thermal responsive pathways, the nanofibers were characterized before and after heating. X-ray photoelectron spectroscopy indicates that an initially mixed Cu(I)/Cu(II) oxidation states transition to predominantly Cu(I) upon heating. A layer of nanofiber was coated on battery pouch foil and evaluated as a candidate thermally sensitive coating. At elevated temperature (100-130°C), nanofiber coating released volatile organic compounds, sulfide and sulfur dioxide as detected using multiple gas sensors. This thermally responsive gas release/sensing approach provides a potential large-area temperature monitoring strategy, which is particularly relevant where direct temperature measurements of individual batteries is impractical. The results established proof of concept for nanofiber-coated battery pouch foil as overtemperature warning platform that can provide alerts when surface temperatures exceed a critical threshold. More broadly, the ability to form interconnected fiber networks positions copper thiolate nanofiber coatings as promising materials for advanced applications.

Ihala Gamaralalage, Chanaka [ORNL] (ORCID:00000002↗
DOE OSTI2022

Experimental and modeling study of water time histories during H 2 S-N 2 O combustion in a shock tube

Here, with the interest in directly burning sour gas in gas turbines, and the fact that even small amounts of H 2 S or its combustion products can alter combustion characteristics, many research studies have been performed to better understand the combustion chemistry of H 2 S. In the present study, the water formation was followed by laser absorption with N 2 O as an oxidant, instead of O 2 . Nitrous Oxide being essentially consumed via N 2 O (+M) ⇌ N 2 + O (+M), the water formation via the H 2 S + O route can then be probed to further validate the models. Three H 2 S/N 2 O mixtures diluted in 98% Ar were studied to cover the following range of equivalence ratios: 0.5, 1.0, and 2.0, over a wide range of temperatures (1580–1940 K) around atmospheric pressure. A chemical kinetic model was then developed first by validating the base N 2 O kinetics mechanism, then by investigating the experimental results presented herein. The N 2 O kinetics mechanism was updated with recent work on the low- and high-pressure limits for N 2 O decomposition (N 2 O (+M) ⇌ N 2 + O (+M)) as well as a review of rate coefficients for N 2 O + H ⇌ N 2 + OH from the literature. Good agreement is shown for H 2 /N 2 O mixtures. Updates were then made to the H 2 S kinetics mechanism, specifically, an update from the literature on SO 2 + H (+M) ⇌ SO + OH (+M) and an adjustment to SO + SH ⇌ S 2 + OH. Additionally, reactions between SH and N 2 O were determined using W1BD data and transition state theory which required the addition of an NNS sub-mechanism. With these updates, the mechanism provides good agreement with the H 2 S/N 2 O experiments.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2023

From electronic structure to model application of key reactions for gasoline/alcohol combustion: Hydrogen-atom abstraction by $CH_3O\dot{O}$ radicals

Hydrogen atom abstraction by methyl peroxy ($CH_3O\dot{O}$) radicals can play an important role in gasoline/ethanol interacting chemistry for fuels that produce high concentrations of methyl radicals. Detailed kinetic reactions for hydrogen atom abstraction by $CH_3O\dot{O}$ radicals from the components of FGF-LLNL (a gasoline surrogate) including cyclopentane, toluene, 1-hexene, n -heptane, and isooctane have been systematically studied in this work. Here, the M06-2X/6-311 ++ G(d,p) level of theory was used to obtain the optimized structure and vibrational frequency for all stationary points and the low-frequency torsional modes. The 1-D hindered rotor treatment for low-frequency torsional modes was treated at M06-2X/6-31G level of theory. The UCCSD(T)-F12a/cc-pVDZ-F12 and QCISD(T)/CBS level of theory were used to calculate single point energies for all species. High pressure limiting rate constants for all hydrogen atom abstraction channels were performed using conventional transition state theory with unsymmetric tunneling corrections. Individual rate constants are reported in the temperature range from 298.15 to 2000 K. Our computed results show that the abstraction of allylic hydrogen atoms from 1-hexene is the fastest at low temperatures. When the temperature increases, the hydrogen atom abstraction reaction channel at the primary alkyl site gradually becomes dominant. Thermodynamics properties for all stable species and high-pressure limiting rate constants for each reaction pathway obtained in this work were incorporated into the latest gasoline surrogate/ethanol model to investigate the influence of the rate constants calculated here on model predicted ignition delay times.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2026

Pressure-dependent kinetics of cyclopentene consumption

Recent experiments have uncovered deficiencies in the ability of current kinetic models for cyclopentane combustion to predict species profiles as a function of O2 concentration. One likely source of these discrepancies flows from the limited availability of relevant kinetic data for key functionalized intermediates, requiring the use of distant analogies to estimate their rates. This work fills this gap for cyclopentene, the most prominent C5 intermediate in low-temperature cyclopentane combustion. Ab initio transition-state-theory based master equation methods are used to calculate pressure-dependent rate constants across four potential energy surfaces (C5H8 + OH, C5H8 + HO2, C5H7, and C5H7 + O2). Potential energy surfaces are characterized at the CCSD(T)-F12/cc-pVDZ-F12//revDSD-PBEP86-D3BJ/def2-TZVPP level of theory. Pressure-dependent rate constants and branching fractions are determined via 1D master equation calculations. The results are used to analyze the competition for radicals between cyclopentene and cyclopentane in the intermediate temperature region, and to identify the dominant pathways flowing from these initiation channels. On the OH initiation surface, we find good agreement with recently published experimental rate constants at high temperature, and we elucidate the complex temperature dependence below 800 K, where 𝜋𝜋 addition becomes prominent. On the HO2 initiation surface, we find an interesting alternative route to bicyclic ether formation via well-skipping over the HO2 adduct well. On the cyclopentenyl radical surface, the major product is cyclopentadiene plus H. On the cyclopentenyl + O2 surface, we find that resonance stabilization largely prevents the allylic radical from undergoing oxidation, whereas the alkylic radical has a somewhat elevated propensity for oxidation pathways relative to cyclopentyl.

ab initio calculations↗
DOE OSTI2021

On the structure sensitivity of and CO coverage effects on formic acid decomposition on Pd surfaces

Using density functional theory calculations, the Pd-catalyzed vapor-phase formic acid decomposition was studied, with a focus on the structure sensitivity and CO coverage effects. A comprehensive reaction network was developed on both the (111) and (100) facets of Pd, at CO coverages of 0 and 5/9 monolayer (ML). Pd(100) was determined to be more reactive than Pd(111) at both CO coverages. The introduction of 5/9 ML CO decreased the activity of both facets significantly, due to destabilization of the surface intermediates and transition states on the CO-decorated surfaces. Three reaction pathways were explored on the clean surfaces: the formate (HCOO) pathway, the carboxyl (COOH) pathway leading to the formation of CO 2 , and the COOH pathway leading to the formation of CO (COOH→CO). Based on the DFT-derived energetics alone, it appears that all three pathways contribute to the reaction on clean Pd, whereas the presence of 5/9 ML of CO inhibits the HCOO pathway on both facets and favors the COOH→CO pathway on the (111) facet, but the COOH→CO 2 one on the (100) facet. Moreover, at high CO coverages, alternative spectator CO-assisted adsorbate decomposition pathways were discovered, which could potentially play a role in formic acid decomposition on Pd catalysts under realistic reaction conditions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Effects of co-adsorbed water on different bond cleavages involved in acetic acid decomposition on Pt (111)

Acetic acid decomposition on Pt (111) in the presence of co-adsorbed water is a good model system for oxygenate decomposition on Pt (111) in aqueous phase, with application in hydrogen production and biomass conversion. Here, in this study, we present a density functional theory (DFT) theory calculation of how co-adsorbed water affects different bond cleavages of acetic acid decomposition on Pt (111). The presence of co-adsorbed water generally enhances O$-$H bond cleavage while inhibiting OC$-$O and OC$-$OH bond cleavage. The influence of co-adsorbed water on C$-$H bond cleavage varies the most and depends on the nature of the transition state and how co-adsorbed water stabilizes the initial and final state. Although these trends are useful as general guidance, they are not sufficient to predict the effect on a complex reaction network such as acetic acid decomposition on Pt (111) which has several parallel reaction paths with similar energies. In the absence of co-adsorbed water, the two lowest energy pathways are decarboxylation (DCX) and decarbonylation (DCN) pathways through a common CH 2 COO intermediate, in which the DCX pathway (C$-$H bond cleavage of CH 2 COO) is more favorable than DCN pathway (OC$-$O bond cleavage of CH 2 COO). In the presence of co-adsorbed water, the energy difference between Csingle bondH bond cleavage and OC$-$O bond cleavage of CH 2 COO increases, suggesting an increase in the favorability of acetic acid decarboxylation (formation of carbon dioxide) over acetic acid decarbonylation (formation of carbon monoxide) on Pt (111).

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

The P( 4 S) + NH( 3 Σ – ) and N( 4 S) + PH( 3 Σ – )reactions as sources of interstellar phosphorus nitride

Phosphorus nitride (PN) is believed to be one of the major reservoirs of phosphorus in the interstellar medium (ISM). For this reason, understanding which reactions produce PN in space and predicting their rate coefficients is important for modelling the relative abundances of P-bearing species and clarifying the role of phosphorus in astrochemistry. Here, in this work, we explore the potential energy surfaces of the $\textrm{P}(^4\textrm{S}) + \textrm{NH}(^3\Sigma^-)$ and $\textrm{N}(^4\textrm{S}) + \textrm{PH}(^3\Sigma^-)$ reactions and the formation of $\textrm{H}(^2\textrm{S}) + \textrm{PN}(^1\Sigma^+)$ through high accuracy ab initio calculations and the variable reaction coordinate transition state theory (VRC-TST). We found that both reactions proceed without an activation barrier and with similar rate coefficients that can be described by a modified Arrhenius equation ( $k(T)=\alpha\!\left( T/300 \right)^{\beta} \exp\!{(\!-\!\gamma/T)})$ with $\alpha=0.93\times 10^{-10}\rm cm^3\,s^{-1}$ , $\beta=-0.18$ and $\gamma=0.24\, \rm K$ for the $\textrm{P} + \textrm{NH} \longrightarrow \textrm{H} + \textrm{PN}$ reaction and $\alpha=0.88\times 10^{-10}\rm cm^3\,s^{-1}$ , $\beta=-0.18$ and $\gamma=1.01\, \rm K$ for the $\textrm{N} + \textrm{PH} \longrightarrow \textrm{H} + \textrm{PN}$ one. Both reactions are expected to be relevant for modelling PN abundances even in the cold environments of the ISM. Given the abundance of hydrogen in space, we have also predicted rate coefficients for the destruction of PN via H + PN collisions.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Ab Initio Direct Dynamics

The reactivity and dynamics of molecular systems can be explored computationally by classical trajectory calculations. The traditional approach involves fitting a functional form of a potential energy surface (PES) to the energies from a large number of electronic structure calculations and then integrating numerous trajectories on this fitted PES to model the molecular dynamics. The ever-decreasing cost of computing and continuing advances in computational chemistry software have made it possible to use electronic structure calculations directly in molecular dynamics simulations without first having to construct a fitted PES. In this “on-the-fly” approach, every time the energy and its derivatives are needed for the integration of the equations of motion, they are obtained directly from quantum chemical calculations. This approach started to become practical in the mid-1990s as a result of increased availability of inexpensive computer resources and improved computational chemistry software. The application of direct dynamics calculations has grown rapidly over the last 25 years and would require a lengthy review article. The present Account is limited to some of our contributions to methods development and various applications. To improve the efficiency of direct dynamics calculations, we developed a Hessian-based predictor-corrector algorithm for integrating classical trajectories. Hessian updating made this even more efficient. Furthermore, this approach was also used to improve algorithms for following the steepest descent reaction paths. For larger molecular systems, we developed an extended Lagrangian approach in which the electronic structure is propagated along with the molecular structure. Strong field chemistry is a rapidly growing area, and to improve the accuracy of molecular dynamics in intense laser fields, we included the time-varying electric field in a novel predictor-corrector trajectory integration algorithm. Since intense laser fields can excite and ionize molecules, we extended our studies to include electron dynamics. Specifically, we developed code for time-dependent configuration interaction electron dynamics to simulate strong field ionization by intense laser pulses. Our initial application of ab initio direct dynamics in 1994 was to CH 2 O → H 2 + CO; the calculated vibrational distributions in the products were in very good agreement with experiment. In the intervening years, we have used direct dynamics to explore energy partitioning in various dissociation reactions, unimolecular dissociations yielding three fragments, reactions with branching after the transition state, nonstatistical dynamics of chemically activated molecules, dynamics of molecular fragmentation by intense infrared laser pulses, selective activation of specific dissociation channels by aligned intense infrared laser fields, angular dependence of strong field ionization, and simulation of sequential double ionization.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

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. 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. Here, 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.

42 ENGINEERING↗