Search NASA⌕ Search

SEARCH · Search NASA

Results for “atomic kinetics”

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 469 records · Page 26
All sourcesNASA NTRSDOE OSTI
NASA NTRS1976

Absolute rate parameters for the reaction of ground state atomic oxygen with dimethyl sulfide and episulfide

It is pointed out that the investigated reaction of oxygen with dimethyl sulfide may play an important role in photochemical smog formation and in the chemical evolution of dense interstellar clouds. Kinetic data were obtained with the aid of the flash photolysis-resonance fluorescence method. The photodecomposition of molecular oxygen provided the oxygen atoms for the experiments. The decay of atomic oxygen was studied on the basis of resonance fluorescence observations. Both reactions investigated were found to be fast processes. A negative temperature dependence of the rate constants for reactions with dimethyl sulfide was observed.

Lee, J. H.↗
DOE OSTI2021

Atomistic modeling of meso -timescale processes with $\mathrm{SEAKMC}$: A perspective and recent developments

On-the-fly kinetic Monte Carlo (kMC) methods have recently garnered significant attentions after successful applications to various atomic-scale problems using a timescale outside the reach of classical molecular dynamics. These methods play a critical role in modeling atomistic meso-timescale processes, and it is therefore essential to further improve their capabilities. In this report we review one of the on-the-fly kMC methods, Self-Evolving Atomistic kinetic Monte Carlo (SEAKMC) and propose two schemes that considerably enhance the efficiency of saddle point searches (SPSs) during the simulations. The performance of these schemes is tested using the diffusion of point defects in bcc Fe. In addition, we discuss approaches to significantly mitigate limitations of these schemes, which further improves their efficiencies. Importantly, these schemes improve the SPS efficiency not only for SEAKMC but also for other on-the-fly kMC methods, broadening the applications of on-the-fly kMC simulations to complex meso-timescale problems.

36 MATERIALS SCIENCE↗
NASA NTRS1986

Hot hydrogen atom reactions moderated by H2 and He

Photolysis experiments were performed on the H2-CD4-NH3 and He-CD4-NH3 systems. The photolysis (1849 A) involved only NH3. Mixtures of H2:CD4:NH3 included all combinations of the ratios (200,400,800):(10,20,40):4. Two He:CD4:NH3 mixtures were examined where the ratios equalled the combinations 100:(10,20):4. Abstraction of a D from CD4 by the photolytically produced hot hydrogen from ammonia was monitored by mass spectrometric determination of HD. Both experiment and semiempirical hot-atom theory show that H2 is a very poor thermalizer of hot hydrogens with excess kinetic energy of about 2 eV. Applications of the hard-sphere collision model to the H2-CD4-NH3 system resulted in predicted ratios of net HD production to NH3 decomposition that were two orders of magnitude smaller than the experimental ratios. On the other hand, helium is found to be a very efficient thermalizer; here, the classical model yields reasonable agreement with experiments. Application of a semiempirical hot-atom program gave quantitative agreement with experiment for either system.

Aronowitz, S.↗
DOE OSTI2024

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.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Atomic-scale structure of ZrO 2 : Formation of metastable polymorphs

Metastable phases can exist within local minima in the potential energy landscape when they are kinetically “trapped” by various processing routes, such as thermal treatment, grain size reduction, chemical doping, interfacial stress, or irradiation. Despite the importance of metastable materials for many technological applications, little is known about the underlying structural mechanisms of the stabilization process and atomic-scale nature of the resulting defective metastable phase. Investigating ion-irradiated and nanocrystalline zirconia with neutron total scattering experiments, we show that metastable tetragonal ZrO 2 consists of an underlying structure of ferroelastic, orthorhombic nanoscale domains stabilized by a network of domain walls. The apparent long-range tetragonal structure that can be recovered to ambient conditions is only the configurational ensemble average of the underlying orthorhombic domains. This structural heterogeneity with a distinct short-range order is more broadly applicable to other nonequilibrium materials and provides insight into the synthesis and recovery of functional metastable phases with unique physical and chemical properties.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2025

Quantum dynamics of the temporary capture of light atoms by superfluid helium nanodroplets at very low collision energies (≈1–13 meV): the case of the hydrogen atom and its isotopes

The capture dynamics of a H atom and isotopic variants [D, T and Q (hypothetical isotope of mass equal to four times the mass of H)] by a superfluid helium nanodroplet (HeND) has been investigated theoretically. The HeND (T = 0.37 K) is ( 4 He) N=400 and a mean field quantum hybrid approach [TDDFT (helium) + quantum wave packet (H, D, T or Q)] at zero angular momentum, is used to explore a rather wide range of very low initial kinetic energies (E k,0 ≈ 10–150 K). The analysis of the capture mechanism shows the existence of a dynamical barrier and a dynamical minimum that play key roles to understand the time evolution of the capture, especially the former property. In general, the H atom shows a different behavior from the other isotopes, with the behavior of T and Q being very similar to each other and the D atom behaving inbetween H and T. Besides, it is worth noting that, in principle, at the very low initial kinetic energies considered only “short” and “long” lived atom⋯HeND collision complexes are formed, i.e., in the atom-helium nanodroplet collision only the temporary capture of the atom takes place. The different behaviors observed have been interpreted considering the faster motion of the H atom when colliding with ( 4 He) N=400 and the more quantum character of the H behavior both due to its significantly lower mass. As far as we know, this is the first quantum dynamics study carried out on the collision of light atoms with HeNDs at very low energies.

Sternberg, Michael [Argonne National Laboratory (A↗
DOE OSTI2021

Catalytic Site Requirements for N2O Decomposition on Cu-, Co-, and Fe-SSZ-13 Zeolites

N2O decomposition is investigated on Cu, Co and Fe-exchanged SSZ-13 zeolite catalysts at relatively low metal loadings. The catalysts are synthesized by solution ion exchange, and subjected to X-ray diffraction (XRD), temperature-programed-reduction by H2 (H2-TPR), temperature-programed-reaction of N2O (N2O-TPR) coupled with in-situ transmission FTIR, and finally steady-state flow reaction tests. At low N2O pressures (< 0.05 kPa), all catalysts display pseudo first-order kinetics. From Arrhenius analysis, Cu and Fe-SSZ-13 display very different apparent activation energies but similar pre-exponential factors, suggesting their similar reaction mechanisms. N2O decomposition follows a dual-site mechanism, occurring on dimeric M-O-M sites in these catalysts, and O2 is formed by the combination of two O ad-atoms from two vicinal metal sites. Under low N2O pressure (0.05 kPa) and first-order kinetic regime, the reaction is limited by N-O cleavage on bare metal active sites. In comparison to Cu-SSZ-13, the much higher N2O decomposition rate over Fe-SSZ-13 is attributed to the much lower activation barriers for the N-O cleavage step. N2O decomposition occurs on isolated Co2+ ions in Co-SSZ-13. The rate-limiting step is N-O cleavage on an O-occupied Co site in the low-pressure first order kinetic regime. This single-site mechanism leads to much higher pre-exponential factors as compared to the dual-site mechanism. This beneficial factor for reaction rate enhancement, however, is compromised by the much higher activation barriers over this catalyst.

Lin, Fan↗
NASA NTRS1976

A temperature dependence kinetics study of the reactions of Cl/2-P-3/2/ with O3, CH4, and H2O2

The temperature dependence of two chlorine atom reactions of considerable fundamental importance to stratospheric chemistry was studied using the technique of flash photolysis-resonance fluorescence. The reactions of interest were: (1) Cl + O3 yields ClO + O2 studied at 220-350 K, and (2) Cl + CH4 yields CH3 + HCl studied at 218-401 K. In addition, the reaction Cl + H2O2 yields HCl + HO2 was studied at 300 K. The corresponding rate constants are provided for the three reactions. The new rate data implies the need to revise downward by a factor of 2.4-3 the magnitude of the ozone perturbation due to the presence of ClO/x/ species in the stratosphere, predicted by earlier model calculations.

Watson, R.↗
NASA NTRS1997

Temperature-Dependent Kinetics Studies of the Reactions Br((sup 2)P3/2) + H2S yields SH + HBr and Br((sup 2)P3/2) + CH3SH yields CH3S + HBr. Heats of Formation of SH and CH3S Radicals

Time resolved resonance fluorescence detection of Br(sup 2)P3/2) atom disappearance or appearance following 266-nm laser flash photolysis of CF2Br2/H2S/H2/N2, CF2Br2/CH3SH/H2/N2, Cl2CO/H2S/HBr/N2, and CH3SSCH3/HBr/H2/N2 mixtures has been employed to study the kinetics of the reactions Br((sup 2)P3/2) + H2S = SH + HBr (1,-1) and Br((sup2)P3/2) + CH3SH = CH3S + HBr (2, -2) as a function of temperature over the range 273-431K. Arrhenius expressions in units of 10(exp -12) cu cm/molecule/s which describe the results are k1 = (14.2 +/- 3.4) exp[(-2752 +/- 90)/T],(k-1) = (4.40 +/- 0.92) exp[(-971 +/- 73)/T],k(2) = (9.24 +/- 1.15) exp[(-386 +/- 41)/T], and k(-2) = (1.46 +/-0.21) exp[(-399 +/-41)/T; errors are 2 sigma and represent precision only. By examining Br((sup 2)P3/2) equilibrium kinetics following 355nm laser flash photolysis of Br2/CH3SH/H2/N2 mixtures, a 298 K rate coefficient of (1.7 +/- 0.5) x 10(exp -10) cu cm/molecule/s has been obtained for the reaction CH3S + Br2 yields CH3SBr + Br. To our knowledge, these are the first kinetic data reported for each of the reactions studied. Measured rate coefficients, along with known rate coefficients for similar radical + H2S, CH3SH, HBr,Br2 reactions are considered in terms of possible correlations of reactivity with reaction thermochemistry and with IP - EA, the difference between the ionization potential of the electron donor and the electron affinity of the electron acceptor. Both thermochemical and charge-transfer effects appear to be important in controlling observed reactivities. Second and third law analyses of the equilibrium data for reactions 1 and 2 have been employed to obtain the following enthalpies of reaction in units of kcal/mol: for reaction 1, Delta-H(298) = 3.64 +/- 0.43 and Delta-H(0) = 3.26 +/-0.45; for reaction 2, Delta-H(298) = -0.14 +/- 0.28 and Delta-H(0) = -0.65 +/- 0.36. Combining the above enthalpies of reaction with the well-known heats of formation of Br, HBr, H2S, and CH3SH gives the following heats of formation for the RS radicals in units of kcal/mol: Delta-H(sub f)(sub 0)(SH) = 34.07 +/- 0.72, Delta-H(sub f)(sub 298)(SH) = 34.18 +/- 0.68, Delta-H(sub f)(sub 0)(CH3S) = 31.44 +/- 0.54, Delta-H(sub f)(sub 298)(CH3S) = 29.78 +/- 0.44; errors are 2 sigma and represent estimates of absolute accuracy. The SH heat of formation determined from our data agrees well with literature values but has reduced error limits compared to other available values. The CH3S heat of formation determined from our date is near the low end of the range of previous estimates and is 3-4 kcal/mol lower than values derived from recent molecular beam photofragmentation studies.

Nicovich, J. M.↗
DOE OSTI2020

Experimental and theoretical study of the soot-forming tendencies of furans as potential biofuels

Sustainably-produced biofuels, including furans derived from lignocellulosic biomass, are being studied as potential next-generation fuels. Biofuels have unique structures that are absent in conventional fuels and they potentially emit fewer particulates than conventional fuels because of the presence of oxygen atoms. In this work, we focus on quantification and measurement of sooting tendency. A set of furans was selected and their sooting tendencies were characterized by Yield Sooting Index (YSI), which was obtained by measuring the line-of-sight spectral radiance (LSSR) of a coflow laminar diffusion flame doped with 1000 ppm of each test fuel. The measured YSIs are lower than traditional gasoline fuels and benzenoid aromatics (e.g., YSIbenzene = 100). We also performed numerical simulations with a detailed chemical kinetic mechanism to analyze the flame chemistry of 2-methylfuran (2MF). 4→5 H-atom migration in the furan ring causes 31.0% of the total decomposition and is the dominant fuel consumption pathway leading to soot. DFT calculations via the G4 composite method were carried out to obtain activation energies (Ea) and bond dissociation energies (BDE). The results showed that the 4→5 H-atom migration in the furan ring has the second lowest Ea and subsequently forms a 1,2-diene-substituted ketone derivative as a stable species. The ring-opening products of the 4→5 H-atom migration are propyne + a linear ketene for 2-n-alkylfurans, and 1-butyne + a methyl-substituted ketene for 2-t-butylfuran. These results can explain the trends of the measured YSIs of 2-alkylfurans.

09 BIOMASS FUELS↗
DOE OSTI2023

Low Catalyst Loading Enhances Charge Accumulation for Photoelectrochemical Water Splitting

Abstract Solar water oxidation is a critical step in artificial photosynthesis. Successful completion of the process requires four holes and releases four protons. It depends on the consecutive accumulation of charges at the active site. While recent research has shown an obvious dependence of the reaction kinetics on the hole concentrations on the surface of heterogeneous (photo)electrodes, little is known about how the catalyst density impacts the reaction rate. Using atomically dispersed Ir catalysts on hematite, we report a study on how the interplay between the catalyst density and the surface hole concentration influences the reaction kinetics. At low photon flux, where surface hole concentrations are low, faster charge transfer was observed on photoelectrodes with low catalyst density compared to high catalyst density; at high photon flux and high applied potentials, where surface hole concentrations are moderate or high, slower surface charge recombination was afforded by low‐density catalysts. The results support that charge transfer between the light absorber and the catalyst is reversible; they reveal the unexpected benefits of low‐density catalyst loading in facilitating forward charge transfer for desired chemical reactions. It is implied that for practical solar water splitting devices, a suitable catalyst loading is important for maximized performance.

Liu, Tianying↗
DOE OSTI2023

Low Catalyst Loading Enhances Charge Accumulation for Photoelectrochemical Water Splitting

Abstract Solar water oxidation is a critical step in artificial photosynthesis. Successful completion of the process requires four holes and releases four protons. It depends on the consecutive accumulation of charges at the active site. While recent research has shown an obvious dependence of the reaction kinetics on the hole concentrations on the surface of heterogeneous (photo)electrodes, little is known about how the catalyst density impacts the reaction rate. Using atomically dispersed Ir catalysts on hematite, we report a study on how the interplay between the catalyst density and the surface hole concentration influences the reaction kinetics. At low photon flux, where surface hole concentrations are low, faster charge transfer was observed on photoelectrodes with low catalyst density compared to high catalyst density; at high photon flux and high applied potentials, where surface hole concentrations are moderate or high, slower surface charge recombination was afforded by low‐density catalysts. The results support that charge transfer between the light absorber and the catalyst is reversible; they reveal the unexpected benefits of low‐density catalyst loading in facilitating forward charge transfer for desired chemical reactions. It is implied that for practical solar water splitting devices, a suitable catalyst loading is important for maximized performance.

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

Insertion products in the reaction of carbonyl oxide Criegee intermediates with acids: Chloro(hydroperoxy)methane formation from reaction of CH 2 OO with HCl and DCl

The reactions of carbonyl oxides (“Criegee intermediates”) with organic and inorganic acids have been postulated to proceed predominantly by an insertion mechanism, and the hydroperoxy-substituted ester products from reactions with organic acids have been spectroscopically observed. Here this work characterizes the products from one of the simplest reactions of carbonyl oxides with inorganic acids, that of formaldehyde oxide, CH 2 OO, with hydrogen chloride, which occurs via a 1,2-insertion in the H–Cl bond. Reactions of both HCl and DCl isotopologues yield product signal at the mass of the insertion product chloro(hydroperoxy)methane and a dissociative ionisation peak at the mass of the protonated (or deuteronated) Criegee intermediate. The isotopic composition of the insertion product has been measured for reaction mixtures where both HCl isotopologues are present, and the H/D ratio of the product is consistently higher (by a factor of 1.6 ± 0.3) than that of the reactants. This isotope selectivity in the products has smaller uncertainty than the ratio of measured rate coefficients and suggests a normal (k H > k D ) kinetic isotope effect in the reaction. Theoretical kinetics calculations predict a small normal kinetic isotope effect for the overall reaction (k H / k D = 1.35 at 20 Torr N 2 and k H / k D = 1.2 for 1 atm N 2 ) but predict a substantial inverse kinetic isotope effect (k D > k H ) for the stabilisation fraction, in disagreement with the experimental observation.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2026

Phase-field model of freeze casting

Directional solidification of water-based solutions has emerged as a versatile technique for templating hierarchical porous materials. However, the underlying mechanisms of pattern formation remain incompletely understood. In this work, we present a detailed derivation and analysis of a quantitative phase-field model for simulating this nonequilibrium process. The phase-field model extends the thin-interface formulation of dilute binary alloy solidification with antitrapping to incorporate the highly anisotropic energetic and kinetic properties of the partially faceted ice-water interface. This interface is faceted in the basal plane normal to the ⟨0001⟩ directions and atomically rough in other directions within the basal plane. On the basal plane, the model reproduces a linear or nonlinear relationship between the interface growth rate and the kinetic undercooling that can be linked to experimental measurements. In both cases, spontaneous parity breaking of the solidification front is observed when the preferred growth direction is aligned with the temperature gradient. This phenomenon leads to the formation of partially faceted ice lamellae that drift laterally in one of the ⟨0001⟩ directions. Here, we demonstrate that the drifting velocity of the ice lamellae is controlled by the kinetics on the basal plane and converges as the thickness of the diffuse solid-liquid interface decreases. Furthermore, we examine the effect of the form of the kinetic anisotropy, which is chosen here such that the inverse of the kinetic coefficient varies linearly from a finite value in the ⟨0001⟩ directions to zero in all other directions within the basal plane, consistent with the assumption that the interface grows in local thermodynamic equilibrium in this plane. Our results indicate that the drifting velocity of ice lamellae is not affected by the slope of this linear relation, and the radius and undercooling at the tip of an ice lamella converge at relatively small slope values. Consequently, the phase-field simulations remain quantitative with computationally tractable choices of both the interface thickness and the slope assumed in the form of the kinetic anisotropy.

Materials science↗
DOE OSTI2022

Hydrogen Storage Performance of Preferentially Oriented Mg/rGO Hybrids

Chemical interactions on the surface of a functional nanoparticle are closely related to its crystal facets, which can regulate the corresponding energy storage properties like hydrogen absorption. In this study, we reported a one-step growth of magnesium (Mg) particles with both close- and nonclose-packed facets, that is, {0001} and {$2\overline{116}$} planes, on atomically thin reduced graphene oxide (rGO). The detailed microstructures of Mg/rGO hybrids were revealed by X-ray diffraction, selected-area electron diffraction, high-resolution transmission electron microscopy, and fast Fourier transform analysis. Hydrogen storage performance of Mg/rGO hybrids with different orientations varies: Mg with preferential high-index {$2\overline{116}$} crystal surface shows remarkably increased hydrogen absorption up to 6.2 wt % compared with the system exposing no preferentially oriented crystal surfaces showing inferior performance of 5.1 wt % within the first 2 h. First-principles calculations revealed improved hydrogen sorption properties on the {$2\overline{116}$} surface with a lower hydrogen dissociation energy barrier and higher stability of hydrogen atoms than those on the {0001} basal plane, supporting the hydrogen uptake experiment. In addition, the hydrogen penetration energy barrier is found to be much lower than that of {0001} because of low surface atom packing density, which might be the most critical process to the hydrogenation kinetics. Here, the experimental and calculation results present a new handle for regulating the hydrogen storage of metal hydrides by controlled Mg facets.

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

Molecular Dynamics Simulation of the Structures, Dynamics, and Aggregation of Dissolved Organic Matter

Dissolved organic matter (DOM) plays a significant role in the transport and transformation of pollutants in the aquatic environment. However, the experimental characterization of DOM has been limited mainly to bulk properties, and the molecular-level interactions among various components of DOM remain to be fully characterized. In this study, we use molecular dynamics (MD) simulations to probe the structural properties of model DOM systems at atomic detail. The 200 ns simulations, validated by available experimental data, reveal processes and mechanisms by which chemical species (cations, peptides, lipids, lignin, carbohydrates, and some low-molecular-weight aliphatic and aromatic compounds) aggregate to form complex DOM. The DOM aggregates are dynamic, consisting of a hydrophobic core and amphiphilic exterior. The lipid tails and other hydrophobic fragments form the core, with hydrophilic and amphiphilic groups exposed to water, making DOM accessible to both polar and nonpolar species. Thus, the lipid component acts as a nucleator, whereas cations (especially Ca 2+ ) connect the molecular fragments on the surface by coordinating with the O-containing functional groups of DOM. The structural details revealed here provide new insights including surface accessible atoms, overall assemblage, and interactions among the molecules of DOM for understanding the kinetics and mechanisms through which DOM interacts with metal and other contaminants.

54 ENVIRONMENTAL SCIENCES↗
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.↗