SEARCH · Search NASA

Results for “Reaction rate constants”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 307 records · Page 17

All sources NASA NTRS DOE OSTI

DOE OSTI2021

Identifying Support Effects in Au-Catalyzed CO Oxidation

Some catalytic oxide supports are more equal than others, with numerous variable properties ranging from crystal symmetry to surface chemistry and electronic structure. As a consequence, it is often very difficult to determine which of these act as the driver of performance changes observed in catalysis. In this work, we hold many of these variable properties constant with structurally similar LnScO 3 (Ln = La, Sm, and Nd) nanoparticle supports with cuboidal shapes and a common Sc-rich surface termination. Using CO oxidation over supported Au nanoparticles as a probe reaction, we observe higher activation energy and a slower rate using NdScO 3 as the support material. This change is found to correlate to the strength of CO 2 binding to the support surface, identified by temperature-programmed desorption measurements. Furthermore, the change is due to differences in the 4f electrons of the lanthanide cations, the cations’ Lewis acidity, and the inductive effect they impose.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2021

Electrochemical Modeling of Iodide Oxidation in Metal-Halide Molten Salts

The oxidation of iodide in NaI-AlBr 3 , NaI-AlCl 3 , and NaI-GaCl 3 molten salts was analyzed using simulation software to extract relevant kinetic parameters. The experimental oxidation potentials were ordered AlCl 3 < AlBr 3 < GaCl 3 , with higher oxidation potentials correlating with softer Lewis acidity of the metal halide. An iodide oxidation and metal halide speciation model was developed and simulated to fit the electrochemical response, enabling determination of electrochemical charge transfer parameters and chemical equilibrium constants. NaI-AlBr 3 displayed the fastest electron transfer rates yet showed the lowest current densities. All salts revealed smaller than expected current densities, explained by equilibrium between various species, where some are not electrochemically active at the studied potentials. These equilibrium reactions are due to the various metal halide species, controlling the reactant concentration of iodide and the resultant current. We hypothesize the electrochemically active iodide species, present as a metal halide monomer (MX 3 I − ), is decreased dramatically from the expected concentration, sequestered as a more stable metal halide dimer species (M 2 X 6 I − ) with a higher oxidation potential. Traditional Tafel analysis of the experimental data supports the validity of the simulations. These results increase understanding of iodide oxidation in low-temperature Lewis acidic molten salts and inform task-specific molten salt design.

Lee, Rose Y. (ORCID:0000000226587220)↗

DOE OSTI2024

Density functional theory (DFT) study of UF 6 hydrolysis: reaction pathways, spectroscopy, and chemical kinetics

Depleted uranium hexafluoride (UF 6 ), a stockpiled byproduct of the nuclear fuel cycle, reacts readily with atmospheric humidity, but the gas-phase reaction mechanism and associated chemical kinetics are poorly understood. During the performance period we undertook development of a state-of-the-art ab initio gas-phase chemical kinetics simulation workflow to model the hydrolysis of uranium hexafluroride (UF 6 ). In doing so, we addressed several outstanding issues in the theoretical treatment of uranium-containing systems. At the outset it was unclear how to generate accurate estimates of kinetic and thermodynamic data for U-containing chemical reactions. Generation of such data has been made routine. Prior to our work, the literature associated with UF 6 hydrolysis were disparate and inaccurate. This body of work provides a modern and comprehensive theoretical assessment of the reaction mechanism, molecular clustering towards deposition, and chemical kinetics. New methodological implementations and software integrations resulting from this work are also highlighted. As much as possible, our predictions were validated against experimental data including particle morphologies, vibrational spectroscopy, atomization enthalpies, and kinetic rate constants. Nevertheless, we were unable to reconcile kinetic measurements with high-accuracy simulations.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2026

The Reaction of Atomic Hydrogen with an Imine: Direct Kinetic Measurements and Computational Studies

Imines are of interest in the combustion modeling of nitrogen-based fuels such as ammonia, as intermediates in atmospheric oxidation of amines, and may offer a new class of environmentally benign substitutes for perhalogenated reagents. Elementary gas-phase rate constants for the consumption of H atoms by 1,1,1,3,3,3-hexafluoro-2-propanimine were measured over 294–736 K using the laser flash photolysis/atomic resonance fluorescence technique. The results are summarized as 3.0 × 10 –11 exp(−14.5 kJ mol –1 /RT) cm 3 molecule –1 s –1 with a confidence interval of ±13%. These are the first temperature-dependent kinetic experiments on an imine of the form R 2 C═NH and thereby provide the only validation of transition-state modeling for such systems. Several pathways are feasible; comparison of theory and experiment suggests the dominant channel is addition of H atoms to the nitrogen atom. On this basis, the observations are rationalized quantitatively. Computations also yield the thermochemistry of (CF 3 ) 2 CNH, with Δ f H 298 = −1232 kJ mol –1 .

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Iron (IV) Formation and the pH Dependent Kinetics of the Fenton Reaction

Despite its widespread importance for biological and environmental chemistry and decades of study, the mechanism underlying the Fenton reaction is still a matter of some controversy. Here, to elucidate the pH dependence of this complex reaction, a new kinetic model is developed to explain the increase in rate and mechanistic shift that occurs from acidic to neutral conditions. This mechanism originated from a re-analysis of a previously proposed model, which neglected explicit iron speciation, leading to unrealistic rate constants. Accounting for speciation suggests a much faster formation rate of Fe(IV), which is estimated to be on the order of 10 6 M −1 s −1 . Expanding on prior kinetic studies that include speciation under acidic conditions, we propose a unified kinetic model that captures the pH-dependent rate acceleration in Fe(II) oxidation by H 2 O 2 , which is a significant step toward resolving the long-standing mechanistic ambiguity of Fenton chemistry.

Cohen, Liron [Lawrence Berkeley National Laborator↗

DOE OSTI2021

Directing Energy Transfer in Halide Perovskite–Chromophore Hybrid Assemblies

Directing the flow of energy and the nature of the excited states that are produced in nanocrystal–chromophore hybrid assemblies is crucial for realizing their photocatalytic and optoelectronic applications. Using a combination of steady-state and time-resolved absorption and photoluminescence (PL) experiments, we have probed the excited-state interactions in the CsPbBr 3 –Rhodamine B (RhB) hybrid assembly. PL studies reveal quenching of the CsPbBr 3 emission with a concomitant enhancement of the fluorescence of RhB, indicating a singlet-energy-transfer mechanism. Transient absorption spectroscopy shows that this energy transfer occurs on the ~200 ps time scale. To understand whether the energy transfer occurs through a Förster or Dexter mechanism, we leveraged facile halide-exchange reactions to tune the optical properties of the donor CsPbBr 3 by alloying with chloride. This allowed us to tune the spectral overlap between the donor CsPb(Br 1–x Cl x ) 3 emission and acceptor RhB absorption. For CsPbBr 3 -RhB, the rate constant for energy transfer (k ET ) agrees well with Förster theory, whereas alloying with chloride to produce chloride-rich CsPb(Br 1–x Cl x ) 3 favors a Dexter mechanism. These results highlight the importance of optimizing both the donor and acceptor properties to design light-harvesting assemblies that employ energy transfer. Furthermore, the ease of tuning optical properties through halide exchange of the nanocrystal donor provides a unique platform for studying and tailoring excited-state interactions in perovskite–chromophore assemblies.

14 SOLAR ENERGY↗

DOE OSTI2025

In Situ Measurement of Stress Corrosion Crack Tip Properties During Environmentally Induced Cracking in Stainless Steel 304L

The influence of the external cathodic area, stress intensity, and electrochemical polarization level on the resultant crack growth rate and crack tip pH and potential was studied for stainless steel 304L in a sodium chloride environment. Before testing, the validity of in situ crack tip property measurements using microprobes and drilled holes was supported through modeling efforts and confirmed with experimental observations. Here, it was found that increasing the R-ratio during corrosion fatigue increased the crack growth rate and decreased the measured crack tip pH of the 304L sample. Decreasing the electrochemical potential (i.e., cathodic polarization) of the bulk sample increased the crack tip pH and stopped measured crack extension under constant K conditions. Conversely, increasing the potential of the bulk sample decreased the crack tip pH and induced cracking under constant K conditions. Decreasing the cathodic area external to the crack increased the measured pH at the crack tip and decreased the crack growth rate under fatigue conditions. For the given configuration, environment, and alloy studied, fatigue crack growth is not sustained by cathodic reactions solely occurring in the crack wake and notch. These results have implications for the ability to predict atmospheric crack growth rates from laboratory, full immersion experiments.

corrosion fatigue↗

DOE OSTI2024

Force-Free Identification of Minimum-Energy Pathways and Transition States for Stochastic Electronic Structure Theories

Here, the accurate mapping of potential energy surfaces (PESs) is crucial to our understanding of the numerous physical and chemical processes mediated by atomic rearrangements, such as conformational changes and chemical reactions, and the thermodynamic and kinetic feasibility of these processes. Stochastic electronic structure theories, e.g., Quantum Monte Carlo (QMC) methods, enable highly accurate total energy calculations that in principle can be used to construct the PES. However, their stochastic nature poses a challenge to the computation and use of forces and Hessians, which are typically required in algorithms for minimum-energy pathway (MEP) and transition state (TS) identification, such as the nudged elastic band (NEB) algorithm and its climbing image formulation. Here, we present strategies that utilize the surrogate Hessian line-search method, previously developed for QMC structural optimization, to efficiently identify MEP and TS structures without requiring force calculations at the level of the stochastic electronic structure theory. By modifying the surrogate Hessian algorithm to operate in path-orthogonal subspaces and at saddle points, we show that it is possible to identify MEPs and TSs by using a force-free QMC approach. We demonstrate these strategies via two examples, the inversion of the ammonia (NH 3 ) molecule and the nucleophilic substitution (S N 2) reaction F – + CH 3 F → FCH 3 + F – . We validate our results using Density Functional Theory (DFT)- and Coupled Cluster (CCSD, CCSD(T))-based NEB calculations. We then introduce a hybrid DFT-QMC approach to compute thermodynamic and kinetic quantities, free energy differences, rate constants, and equilibrium constants that incorporates stochastically optimized structures and their energies, and show that this scheme improves upon DFT accuracy. Our methods generalize straightforwardly to other systems and other high-accuracy theories that similarly face challenges computing energy gradients, paving the way for highly accurate PES mapping, transition state determination, and thermodynamic and kinetic calculations at significantly reduced computational expense.

Iyer, Gopal R.↗

DOE OSTI2020

Modified Isotherm Modeling to Predict Cs Exchange with Crystalline Silicotitanate in Tank Waste Simulants

The U.S. Department of Energy is working to expedite processing of Hanford tank waste supernate at the Hanford Waste Treatment and Immobilization Plant. To support this goal, Washington River Protection Solutions is designing a Tank Side Cesium Removal (TSCR) system for suspended solids and cesium (Cs/ 137 Cs) removal from Hanford tank waste supernate. The ion exchange media selected for Cs removal at TSCR is crystalline silicotitanate (CST) that is manufactured in a nearly spherical form by Honeywell UOP LLC (UOP; Des Plaines, IL) as product IONSIV® R9140-B (Na form). The Zheng Anthony Miller (ZAM) isotherm model (Zheng et al. 1997) is a multicomponent ion exchange model used to predict the exchange of Group I metals onto CST. The ZAM isotherm has historically been used to predict Cs distribution values from Hanford and Savannah River Site (SRS) tank waste simulants. Figure S.1 summarizes model predictions from the ZAM isotherm that indicate poor prediction of Cs distribution values for simple and complex simulants with the engineered form of CST, IONSIV® R9140-B, and IONSIV® R9120-B where the solid line indicates a perfect fit by the model. The dotted lines indicate ±20% error. Batch contact testing with Hanford tank waste complex and simple simulants was used in conjunction with SRS simulants to experimentally determine Cs distribution values using a modification to the original isotherm model. The experimentally determined maximum Cs capacity for IONSIV® R9140-B CST in both the simple and complex matrices was found to be 0.53±0.3 mmoles Cs/g of CST. This value is not drastically different from the maximum Cs capacity of 0.58 mmoles Cs/g TAM-5 reported by Zheng et al. (1997). However, it is important to note that TAM-5 (commercially IONSIV® IE-910) is a powder. Hamm et al. (2002) determined that a dilution factor was needed to account for the Zr(OH)2 binder in the engineered form of CST. Hamm et al. determined that a dilution factor of 0.68 was appropriate to account for binder contribution and correct overprediction of Cs exchange on the engineered form of CST in ZAM calculations. This reduced the total capacity from 0.58 mmol/g with TAM-5 to 0.39 mmol/g for the engineered form of CST (Hamm et al. 2002). Despite substituting the experimentally determined maximum Cs capacity of 0.55 mmoles Cs/g for the literature-reported capacity of 0.39 mmoles Cs/g, it was determined that additional modifications to the model’s equilibrium rate constants were necessary in refining the isotherm model. The modified model overpredicted K+ uptake by the CST when compared to digested CST results described by Campbell et al. (2019). The modified model was further revised to omit three of the five K+ exchange equilibrium reactions described by ZAM to reduce the additional K+ loading seen by the model. Figure S.2 summarizes the revised model isotherm predictions plotted against measured Kd values.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DOE OSTI2024

Formation of the oxyl’s potential energy surface by the spectral kinetics of a vibrational mode

One of the most reactive intermediates for oxidative reactions is the oxyl radical, an electron-deficient oxygen atom. The discovery of a new vibration upon photoexcitation of the oxygen evolution catalysis detected the oxyl radical at the SrTiO3 surface. The vibration was assigned to a motion of the sub-surface oxygen underneath the titanium oxyl (Ti–O●−) created upon hole transfer to (or electron extraction from) a hydroxylated surface site. Evidence for such an interfacial mode is derived from its spectral shape, which exhibited a Fano resonance—a coupling of a sharp normal mode to continuum excitations. Here, this Fano resonance is utilized to derive precise formation kinetics of the oxyl radical and its associated potential energy surface (PES). From the Fano lineshape, the formation kinetics are obtained from the anti-resonance (the kinetics of the coupling factor), the resonance (the kinetics of the coupled continuum excitations), and the frequency integrated spectrum (the kinetics of the normal mode’s cross-section). All three perspectives yield logistic function growth with a half-rise of 2.3 ± 0.3 ps and a time constant of 0.48 ± 0.09 ps. A non-equilibrium transient associated with photoexcitation is separated from the rise of the equilibrated PES. The logistic function characterizes the oxyl coverage at the very initial stages (t ∼ 0) to have an exponential growth rate that quickly decreases toward zero as a limiting coverage is reached. Such time-dependent reaction kinetics identify a dynamic activation barrier associated with the formation of a PES and quantify it for oxyl radical coverage.

Chemistry↗

DOE OSTI2022

The Reaction of HO2 and CH3O2: CH3OOH Formed from the Singlet Electronic State Surface

High-level coupled-cluster calculations in combination with two-dimensional master equation simulations were used to study the HO2 + CH3O2 reaction, which plays an important role in the oxidation of methane and hydrocarbons in the Earth’s atmosphere and low-temperature combustion. The main reaction pathways taking place on the lowest-lying triplet and singlet potential energy surfaces (PES) were characterized. Interestingly, methyl hydroperoxide (CH3OOH), the sole product, could be produced from both the triplet and singlet PESs, with a ratio of roughly 9:1. Formaldehyde is not made as a primary product, but can be formed via secondary chemistry. The formation of methyl tetraoxide (MTO) from the singlet PES is unimportant. The calculated reaction rate coefficients were found to be practically pressure-independent for p ≤ 760 Torr and can be given by k(T)=2.75×10−13×e+1.75 kcal mol−1/RT (in cm3/s), an expression useful for kinetics modeling over the range T = 200–800 K. The rate constant has a slight negative Arrhenius energy dependence of about −1.75 kcal mol–1, falling about a factor of 30 from 200 K to 800 K.

2DME↗

DOE OSTI2023

Mass transfer in catalytic depolymerization: External effectiveness factors and serendipitous processivity in stagnant and stirred melts

Several heterogeneous catalysts are being developed to recycle plastics. Most operate in viscous polymer melts, where external mass transfer effects could limit the supply of co-reactants to active sites. External mass transfer can also impede the diffusion of long chain products away from the catalyst after each cut. Product egress limitations could potentially confer unintentional processivity to catalyst operation, i.e. a tendency for the catalyst to repeatedly cut the same chain after an initial encounter. We formulate reaction–diffusion equations to quantify mass transfer effects on the co-reactant transport to the catalyst and the degree of serendipitous processivity. Results are developed for catalysts in stagnant or stirred melts, with simple expressions involving Damkohler, Peclet, and Sherwood numbers, i.e. dimensionless combinations of rate constants, catalyst particle size, polymer diffusivities, and shear rates (where applicable). In conclusion, we estimate the impact of these effects for a spherical core–shell catalyst.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Tracking Spatiotemporal Electric Potential in Batteries Using High-Resolution Operando X‑ray Transmission Imaging

The formation of significant concentration gradients across electrolytes in batteries affects the rate at which electrochemical reactions occur. In this work, we use high-resolution operando X-ray transmission imaging to capture spatiotemporal salt concentration profiles c(x,t) in a symmetric cell comprising a polymer electrolyte sandwiched between two lithium–indium alloy electrodes during a constant-current experiment followed by open-circuit relaxation. The decay of open-circuit potential is related to the concentration dependence of the potential across concentration cells, U. We show how operando c(x,t) data can be used to calculate the spatiotemporal electric potential “inside” the polarized electrolyte. We track the spatial- and time-dependent cell potential during the constant-current step and distinguish its two contributions: a concentration overpotential governed by U. and an ohmic contribution governed by ionic conductivity. Over most of the time window, the concentration overpotential dominates. At steady state, it is a factor of 7 larger than the ohmic contribution. Such findings indicate that efforts to design new polymer electrolytes should focus on minimizing concentration gradients.

Electrical conductivity↗

DOE OSTI2020

One-electron redox kinetics of aqueous transition metal couples Zn 2+/+ , Co 2+/+ , and Ni 2+/+ using pulse radiolysis

The one-electron redox potentials for aqueous metal couples Co 2+/+ and Ni 2+/+ have been investigated by using pulse radiolysis using their reactions with a series of reference compounds to establish the most positive upper limits of E 0 . Experiments with Zn + were also carried out to confirm the characteristic shape of the expected reduction kinetics. Both formate ions and t-BuOH were employed to scavenge ˙OH radicals and ˙H atoms. Kinetics and fitted first and second order reaction rates have been reported for reactions with methyl viologen, fluorescein, Ru(NH 3 ) 6 3+ , Co(en) 3 2+ , Co(sepulcrate) 3+ , Ru(bpy) 3 2+ , Cr(bpy)33+, and Ni(Me 6 [14]4,11-dieneN4) 2+ . Previous work demonstrated that both Co 2+ and Ni 2+ can be reduced by CO 2 ˙- radicals, giving a negative E 0 limit of -1.9 V vs. SHE. A definite reaction of Ni + with fluorescein di-anions provides a new upper limit of the Ni 2+/+ couple as -0.906 V vs. SHE. The reaction of Co + with Ru(bpy) 3 2+ has been confirmed, giving E 0 = -1.3 V vs. SHE as a rigorous upper limit of the Co 2+/+ couple. In the case of Co 2+/+ , kinetics were complicated by a self-catalyzed metal clustering phenomenon. Initiation rate constants of this process have also been reported.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2023

General Kinetic Model for pH Dependence of Proton-Coupled Electron Transfer: Application to an Electrochemical Water Oxidation System

The pH dependence of proton-coupled electron transfer (PCET) reactions, which are critical to many chemical and biological processes, is a powerful probe for elucidating their fundamental mechanisms. Herein, a general, multichannel kinetic model is introduced to describe the pH dependence of both homogeneous and electrochemical PCET reactions. According to this model, a weak pH dependence can arise from the competition among multiple sequential and concerted PCET channels involving different forms of the redox species, such as protonated and deprotonated forms, as well as different proton donors and acceptors. The contribution of each channel is influenced by the relative populations of the reactant species, which often depend strongly on pH, leading to complex pH dependence of PCET apparent rate constants. This model is used to explain the origins of the experimentally observed weak pH dependence of the electrochemical PCET apparent rate constant for a ruthenium-based water oxidation catalyst attached to a tin-doped In2O3 (ITO) surface. The weak pH dependence is found to arise from the intrinsic differences in the rate constants of participating channels and the dependence of their relative contributions on pH. This model predicts that the apparent maximum rate constant will become pH-independent at higher pH, which is confirmed by experimental measurements. Our analysis also suggests that the dominant channels are electron transfer at lower pH and sequential PCET via electron transfer followed by fast proton transfer at higher pH. Furthermore, this work highlights the importance of considering multiple competing channels simultaneously for PCET processes.

Catalysts↗

DOE OSTI2024

A generalization of the shock invariant relationship

Shock invariant relationship, which was conceived for inert shock waves to derive the 4th power relationship between shock pressure and maximum strain rate, is generalized for reactive shock waves such as Chapman–Jouget detonation and shock-induced vaporization. The generalization, based on the first-order reaction models, is a power function relationship between overall dissipated energy (Δe dis ) and reaction time Δτ such that Δe dis Δτ 1/α = constant, where the power coefficient α is found to be in the range of 2/3–4. Experimental data, though scarce, are consistent with the generalization. Implication of the generalization for inert shocks is also considered and suggests a broad range of the 4th power coefficient including an inequality equation that constrains the shock and particle velocity relationship.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2024

Deuterium Permeation Through a Self-Supported Palladium-Silver Membrane in Helium Gas Mixtures

Palladium membranes and membrane reactors can separate and purify tritium from impurities in the plasma exhaust processing system for the fusion energy fuel cycle. Membranes can also act as a continuous separation method to remove tritium from helium streams in the breeder blanket tritium extraction system, such as from the purge gas of solid breeders. To investigate the potential of these membranes for these applications, we performed a deuterium permeation campaign with a self-supported palladium-silver (Pd-25Ag wt%) tube of 0.15 m length, 76 µ m wall thickness, and 3.0 mm inner diameter. A gas mixture of 3.95% D 2 and a balance of He was delivered to the inside of the tube and D permeated radially outwards through the membrane into a vacuum chamber. Permeation experiments were conducted at 300 ° C, 350 ° C, 400° C, and 450°C, from 100 to 1000 sccm total flow rate, and with total pressures of 90, 150, 190, and 250 kPa. Further, deuterium permeation was consistently lower than predicted from diffusion-limited permeation models, thus we developed a transport model that included gas-phase mass transfer and surface reactions to model experimental results. The dissociation constant was optimized to fit the developed model to experimental data.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2026

Nuclear–Electronic Orbital General Rate Theory: Predicting Hydrogen Kinetic Isotope Effects in the Deep Tunneling Regime

Hydrogen transfer is a critical component of many chemical and biological processes. The ratio of rate constants for hydrogen and deuterium transfer defines the H/D kinetic isotope effect (KIE), which is a powerful tool for elucidating hydrogen transfer mechanisms. Interpretation of experimental H/D KIEs relies on accurate and affordable computational methods. However, due to their light mass, hydrogen and deuterium can undergo tunneling, which is challenging to describe in multidimensional molecular systems. Herein, we introduce the nuclear–electronic orbital general rate theory (NEO-GRT), which enables the efficient prediction of H/D KIEs based on full-dimensional molecular quantum chemistry calculations. The NEO-GRT approach describes the hydrogen transfer rate constant with a general expression that spans the vibrationally adiabatic and nonadiabatic hydrogen tunneling regimes. The input quantities are computed using NEO density functional theory, which treats the transferring hydrogen or deuterium nucleus quantum mechanically on the same level as the electrons. We investigate two intramolecular proton transfer reactions in organic molecules at temperatures down to 50 K to evaluate the performance of NEO-GRT by comparison to transition state theory and ring-polymer instanton theory. The KIEs computed with NEO-GRT agree with those calculated using ring-polymer instanton theory for the full-dimensional molecular systems at the same level of electronic structure theory. This agreement indicates that NEO-GRT captures the deep hydrogen tunneling effects, in contrast to transition state theory, which neglects such effects. Given its relatively low computational cost, NEO-GRT is a promising approach for predicting H/D KIEs in large organic and organometallic systems.

Hydrogen↗