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 289 records · Page 16
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

H 2 O-assisted O 2 reduction by H 2 on Pt and PtAu bimetallic nanoparticles: Influences of composition and reactant coverages on kinetic regimes, rates, and selectivities

Hydrogen peroxide (H 2 O 2 ) can replace hazardous oxidants in industrial processes but is currently too expensive for many such applications. While direct synthesis of H 2 O 2 (H 2 + O 2 → H 2 O 2 ) may reduce costs in comparison to incumbent technology, current catalysts lack the requisite stability and selectivity. Here, we examine the direct synthesis of H 2 O 2 on bimetallic Pt 1 Au x (0 ≤ x ≤ 230) and Pt catalysts at steady-state in pure water and relate kinetic parameters for H 2 O 2 and H 2 O formation to possible active site structures informed by complementary characterization methods. X-ray photoelectron spectra show significant Pt surface enrichment compared to the bulk composition. Analysis of infrared spectra of mixed monolayers of 12 CO* and 13 CO* indicate that Pt and Au form substitutional surface alloys. The Pt 1 Au x nanoparticles with the greatest mole fractions of Au predominantly expose Pt monomers (i.e., isolated Pt atoms), yet Pt atoms exposed upon all these nanoparticles possess electronic structures distinct from bulk Pt. Despite these differences, rate measurements are consistent with product formation through proton-electron transfer pathways for all Pt 1 Au x catalysts. In situ XAS indicate that Pt remains metallic during H 2 O 2 synthesis. Under the most oxidizing conditions, selectivities toward H 2 O 2 increase strongly with the Au to Pt ratio from 2% for monometallic Pt to 85% for Pt 1 Au 170 . However, selectivities are similar among all catalysts within reducing conditions. Comparisons of apparent activation enthalpies for the formation of H 2 O 2 and H 2 O across these catalysts and the range of conditions suggest that Pt monomers within Au provide the greatest selectivities for H 2 O 2 formation, because these active sites present high barriers for O-O bond rupture. Further, selectivities decrease with increasing ratios of H 2 to O 2 pressures, because Pt atoms aggregate and form oligomers that readily dissociate dioxygen intermediates. The combined use of spectroscopy, kinetics, and concepts employed in reaching these conclusions take inspiration from the legacy of Prof. Michel Boudart, and specifically his elegant methods for interrogating bimetallic catalysts.

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

Laser-Irradiated Holey Graphene-Supported Single-Atom Catalyst towards Hydrogen Evolution and Oxygen Reduction

Single-atom catalysts (SAC) can boost the intrinsic catalytic activity of hydrogen evolution reaction (HER) and oxygen reduction reaction (ORR). However, the challenge remains due to the complex synthesis process and insufficient stability. A sustainable approach is applied to synthesizing SACs through laser irradiation and gaining mesoporous graphene oxide (MGO). The surface dangling bonds of nitrogen-doped MGO (NMGO) extract metal atoms species from Co or Fe metal foams and convert them to SAC via an appropriate synthesis approach. Notably, the Co-NMGO electrocatalyst requires low potentials of 146 mV to convey a current density of 10 mA cm -2 towards HER. Similarly, the Fe-NMGO electrocatalyst offers an onset of 0.79 V towards ORR in acidic solution. The individual metal atoms are confirmed via aberration-corrected scanning transmission electron microscopy, and X-ray absorption near-edge structure and extended X-ray absorption fine structure. Density functional theory calculations by applying the grand canonical potential kinetics model revealed that Co-NMGO shows the optimum free reaction energy of -0.17 eV at -0.1 V for HER, and Fe-NMGO has less limiting potential than that of Co-NMGO for ORR case. Furthermore, this work opens a new approach towards the synthesis of SAC and its mechanistic understandings.

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

Probing the antiknock effect of anisole through an ignition, speciation and modeling study of its blends with isooctane

In order to unravel the reaction pathways relevant to anisole co-oxidation within a fuel blend, a detailed study of isooctane/anisole blends was performed with the ULille RCM. Ignition delays as well as mole fraction profiles were measured during a two-stage ignition delay using sampling and GC techniques. These results are used to validate a kinetic model developed from ab initio calculations for the most relevant rate constants which included H-atom abstraction reactions from anisole, and reactions on the potential energy surfaces of methoxyphenyl + O 2 and anisyl + O 2 . Pressure dependent rate constants were computed for the methoxyphenyl + O 2 and anisyl + O 2 reactive systems using master equation code analysis. The new kinetic model shows good agreement with the experimental data. Dual brute-force sensitivity analysis was performed, on both first- and second-stages of ignition, allowing the identification of the most important reactions in the prediction of both ignition delays. It was observed that while pure anisole does not show NTC behavior, a 60/40 isooctane/anisole blend displays such behavior, as well as two-stage ignition. This suggests anisole addition may not be as beneficial to knock resistance as expected from its high octane number. The kinetic modeling results demonstrate the importance of H-abstraction reactions both from the methoxy group and from the aryl ring in ortho-position and the addition of the resultant radicals to O2, mostly leading to the formation of polar or non-aromatic products.

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

First-principles elucidation of defect-mediated Li transport in hexagonal boron nitride

Hexagonal boron nitride (hBN) is a promising candidate as a protective membrane or separator in Li-ion and Li–S batteries, given its excellent chemical stability, mechanical robustness, and high thermal conductivity. In addition, hBN can be functionalized by introducing defects and dopants, or be directly integrated into other active components of batteries, which further augments its appeal to the field. Here, we use first-principles simulations to evaluate the role of atomic defects in hBN in regulating the Li-ion diffusion mechanism and associated kinetics. Specifically, the following four distinct types of vacancy defects are considered: isolated single B and N vacancies, a B–N vacancy pair, and a B 3 N vacancy cluster. It is found that these defect sites generally favor Li intercalation and out-of-plane diffusion but slow down in-plane Li-ion diffusion due to a strong Li trapping effect at the defect sites. Such a trapping effect is, however, highly local such that it does not necessarily affect the overall Li-ion conductivity in defected hBN layers. The present systematic evaluation of the impact of atomic defects on Li ion migration and accompanied charge analysis of hBN lattice in response to Li-ion diffusion provide a mechanistic understanding of Li-ion transport behavior in defected hBN and highlight the potential of defect engineering to achieve optimal material performance.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
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↗
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↗
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.↗