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 577 records · Page 32
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Dynamics of HCN, HNC, and HNCO Formation in the 193 nm Photodissociation of Formamide

Formamide (NH 2 CHO) is the simplest molecule containing a peptide linkage [-NH-C($=$O)-], and it plays an essential role in the study of prebiotic chemistry. Exposure to UV irradiation allows formamide to decompose and act as a prebiotic feedstock in the formation of nucleobases and other necessary starting materials. The photodissociation mechanism of gaseous formamide at 193 nm is studied using (a) chirped-pulse Fourier transform millimeter-wave spectroscopy in the 260-290 GHz spectral region in a room-temperature flow-tube reactor at 1 mu bar pressure, (b) a combination of electronic structure theory, transition state theory, and quasiclassical trajectories, and (c) the Active Thermochemical Tables. The HCN and HNC photoproducts of hydrogenated (NH 2 CHO) and deuterated (NH 2 CDO and ND 2 CHO) formamide precursors are examined to gain insight into the photodissociation mechanism. Here, the theoretical investigation has characterized the main pathway leading to each of the HCN/HNC isomers from the precursor isotopologues. The theoretical branching ratio [HNC]/[HCN] = 2.1 for nascent photofragments agrees with the experiment. The effect of the postphotolysis HNC ↔ HCN isomerization on the [HNC]/[HCN] ratio is predicted. We report the experimental branching ratio [HNCO]: ([HNC] + [HCN]) = 12 ± 3 and propose that most of HNCO originates from dissociation on the S-1 electronic state of formamide.

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

Hydrolysis Reaction Pathways of Thorium Oxide Nanoclusters

Density functional theory benchmarked by correlated molecular orbital theory is used to develop a fundamental and predictive understanding of the interaction of thorium oxide nanoclusters with gas phase water to provide insight into nuclear-waste storage, production of thorium nuclear reactor fuels, and reprocessing of spent fuel. The structures of Th n O 2n (n = 3 – 6) clusters and their interactions with water have been studied at the B3LYP, MP2, and CCSD(T) levels. Hydrolysis is initiated by the formation of Lewis acid-base adducts, with relative H 2 O binding energies (physisorption) ranging from −15 kcal/mol to −22 kcal/mol. The initial H 2 O physisorption energy is ca. −21 kcal/mol regardless of the cluster size and is consistent with the experimentally obtained initial adsorption energy on a thorium dioxide surface. The physisorption enthalpies for additional water molecules can be affected by the presence of terminal groups OH groups generated by proton transfer to a Th-O near the site of adsorption. The hydrolysis products (chemisorption) form either bridging or terminal hydroxides. More exothermic hydrolysis steps were predicted for the formation of terminal hydroxides as compared to the formation of bridging hydroxides. Here, the calculated transition state barriers for transfer of protons from bound water complexes to form the chemisorption products are very low. Overall, water readily reacts with thorium oxide clusters preferring hydroxide products over hydrated complexes. First and second order fits were predicted for the combined physisorption and chemisorption energies for the hydrolysis of thorium oxide clusters. Finally, ionization energies and electron affinities were calculated as were HOMO-LUMO gaps to provide additional insights into the properties of the thorium oxide and hydroxide clusters.

Adsorption↗
DOE OSTI2025

PySIDT: Subgraph Isomorphic Decision Trees for Molecular Property Prediction

Accurate molecular property prediction is important across all fields of chemistry. Deep neural networks (DNNs) have become increasingly popular due to their ability to train automatically, avoiding the incredibly tedious process of constructing and extending traditional property estimation schemes. However, DNNs require large amounts of training data, are challenging to interpret, require large amounts of memory to load even during inference, and have severe difficulties incorporating qualitative chemical knowledge, which are often desired for molecular property prediction tasks. Here, in this study, we present PySIDT (https://github.com/zadorlab/PySIDT), a software for training and running inference on Subgraph Isomorphic Decision Trees (SIDTs). SIDTs are graph-based decision trees made of nodes associated with molecular substructures. Inference is done by descending target molecular structures down the decision tree to nodes with matching subgraph isomorphic substructures and making predictions based on the final (most specific) nodes matched. SIDTs scale down well to dataset sizes much smaller than is feasible for DNNs. As trees of molecular substructures, SIDTs are inherently readable and easy to visualize, making them easy to analyze. They are also straightforward to extend and retrain, facilitate uncertainty estimation, and enable easy integration of expert knowledge. We demonstrate the SIDT approach discussing its application to a diverse range of molecular prediction tasks: rate coefficient estimation, diffusion coefficient estimation, thermochemistry estimation, transition state bond stretch prediction, p K a prediction, stability of molecular structures, stability of surface structures, and prediction of surface lateral interaction energetics. Additionally, we demonstrate the power of the SIDT algorithms in two direct learning curve vanilla comparisons with the popular DNN-based software Chemprop and the popular gradient boosted trees-based software XGBoost on enthalpy of formation and rate coefficient prediction tasks. In particular, in the enthalpy of formation case, vanilla PySIDT is able to outperform vanilla Chemprop and XGBoost across the full range of training/validation set sizes out to 11,560 data points.

Johnson, Matthew Sean [Sandia National Laboratorie↗
DOE OSTI2026

Unimolecular Dynamics of Partially Deuterated Hydroperoxyalkyl Intermediates (•QOOD) in Cyclohexane and Cyclopentane Oxidation

The hydroperoxyalkyl radicals (•QOOH) formed in cyclohexane and cyclopentane oxidation with carbon radical center at the β site relative to the hydroperoxy group have recently been observed through their infrared fingerprint and time- and energy-resolved unimolecular dissociation dynamics to hydroxyl (OH) radical and cyclic ether products. Partial deuteration shifts the IR transitions associated with OH/OD motion to lower frequencies, providing a new window on the IR spectroscopy and unimolecular dynamics of these transient intermediates. Specifically, the overtone OD stretch (2ν OD ) has been identified at 5215.0 cm –1 and 5216.0 cm –1 for the β-QOOD intermediates in cyclohexane and cyclopentane oxidation, respectively. Unimolecular decay rates of the •QOOD intermediates at these energies have been obtained from the time-resolved appearance of OD products. The experimental rates are compared with statistical microcanonical rates evaluated using RRKM theory, including heavy-atom tunneling associated with simultaneous O–O bond elongation and C–O–O angle contraction along the reaction pathway. The experimental rates provide a further test of the previously determined transition state (TS) barriers, which were computed utilizing a benchmarking approach that builds on higher-level reference calculations for the smaller ethane oxidation system. Here, the experimental rate measurements for both •QOOD intermediates agree well with the computed rates after accounting for small changes in zero-point energy and a minor empirical adjustment of the TS barriers in both systems.

Ethers↗
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 OSTI2020

Resolving Heterogeneous Dynamics of Excess Protons in Aqueous Solution with Rate Theory

Rate theories have found great utility across the chemical sciences by providing a physically transparent way to analyze dynamical processes. Here we demonstrate the benefits of using transition state theory and Marcus theory to study the rate of proton transfer in HCl solutions. By using long ab initio molecular dynamics simulations we show that good agreement is obtained between these two different formulations of rate theory and how they can be used to study the pathways and life-time of proton transfer in aqueous solution. Since both rate theory formulations utilize identical sets of molecular data, this provides a self-consistent theoretical picture of the rates of aqueous phase proton transfer. Specifically, we isolate and quantify the rates of proton transfer, ion-pair dissociation, and solvent exchange in concentrated HCl solutions. Our analysis predicts a concentration dependence to both proton transfer and ion-pairing. Moreover, our estimate of the life-time for the Zundel species is 0.8 ps and 1.3 ps for 2M and 8M HCl, respectively. In this paper, we demonstrate that concentration effects, can indeed be quantified through the combination of state-of-the-art simulation and theory and provides a picture of the important correlations between the cation (hydronium) and the counter-ion in acid solutions.

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

Simulations of the Biodegradation of Citrate Based Polymers for Artificial Scaffolds Using Accelerated Reactive Molecular Dynamics

Here, we explore the reactivity and mechanical properties of poly(1,6-hexanediol-co-Citric acid) via ReaxFF molecular dynamics simulations. We implement an accelerated scheme within the ReaxFF framework to study the hydrolysis reaction of the polymer which is provided with sufficient amount of energy known as the restrain energy after a suitable pre-transition state configuration is obtained to overcome the activation energy barrier and the desired product is obtained. The validity of the ReaxFF force field is established by comparing the ReaxFF energy barriers of ester and ether hydrolysis with benchmark DFT values in literature. We perform chemical and mechanical degradation of polymer chain bundles at 300 K. We find that ester hydrolyzes faster than ether due to lower activation energy barrier of the reaction. The selectivity of the bond-boost scheme has been demonstrated by lowering the boost parameters of the accelerated simulation which almost stops the ether hydrolysis. Mechanical degradation of pre-hydrolyzed and intermittent hydrolyzed polymer bundles is performed along the longitudinal direction at two different strain rates. We find that the tensile modulus of the polymers are increases with increase in strain rates which shows that polymers show a strain dependent behavior. The tensile modulus of the Polyester-ether is higher than Polyester yet reaches yield-stress faster than Polyester. This makes Polyester more ductile than Polyester-ether.

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

Molecular Properties and Chemical Transformations Near Interfaces

The properties of bulk water and aqueous solutions are known to change in the vicinity of an interface and/or in a confined environment, including the thermodynamics of ion selectivity at interfaces, transition states and pathways of chemical reactions, and nucleation events and phase growth. Here we describe joint progress in identifying unifying concepts about how air, liquid, and solid interfaces can alter molecular properties and chemical reactivity compared to bulk water and multicomponent solutions. Furthermore, we also discuss progress made in interfacial chemistry through advancements in new theory, molecular simulation, and experiments.

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

Theoretical Investigation of Solvent Effects on the Hydrodeoxygenation of Propionic Acid over a Ni(111) Catalyst Model

The effect of two solvents, liquid water and 1,4-dioxane, has been studied from first principles on the hydrodeoxygenation of propionic acid over a Ni (111) catalyst surface model. A mean-field microkinetic model was developed to investigate these effects at a temperature of 473 K. Under all reaction conditions, a decarbonylation mechanism is favored significantly over a decarboxylation pathway. Although no significant solvent effects were observed on the decarbonylation rate, a substantial solvent stabilization of two key surface intermediates in the decarboxylation mechanism, CH 3 CCOO and CH 3 CHCOO, lead to a notable increase of the decarboxylation rate by two orders of magnitude in liquid water and by one order of magnitude in liquid 1,4-dioxane. Furthermore, a significant solvent stabilization of the transition state of C-H bond cleavage of the α-carbon of CH 3 CHCO, relative to the stabilization of the C-C bond cleavage of the α-carbon of CH 3 CHCO, leads to a change in dominant pathway in the liquid phase environments. Finally, a sensitivity analysis shows that the C-OH bond cleavage of propionic acid and C-C bond cleavage of the α-carbon of CH 3 CHCO are the most rate controlling states in the gas phase. In contrast, in solvents the dehydrogenation of CH 3 CHCO becomes the most influential step. This shift in rate controlling state is attributed to the solvent effect on the dehydrogenation of CH 3 CHCO, which is facilitated in aqueous phase. Altogether, it is likely that the investigated (111) facet of Ni is not active for the hydrodeoxygenation of propionic acid in neither the gas nor liquid phase and other Ni facets or phases must be responsible for the experimentally observed kinetics.

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

Comparative Study on the Machine Learning-Based Prediction of Adsorption Energies for Ring and Chain Species on Metal Catalyst Surfaces

Computation of adsorption and transition state energies for a large number of surface intermediates for numerous active site models pose significant computational overhead in computational screening of catalysts. Machine learning (ML) techniques can be used to predict part of these energies. To predict the energies, ML models need to be fed appropriate metal and species descriptors. For complex surface chemistries, the structures of the intermediate species can vary greatly. In this paper, working with the hydrodeoxygenation of succinic acid on six different metal surfaces, we have studied the effect of linear and non-linear ML models used along with pen-and-paper based species descriptors and two categories of metal descriptors on two different categories of intermediate species: chain and ring. More specifically, our computations include the prediction of chain species when trained on only chain species and also when trained on both chain and ring species. Similar computations were performed for predictions of ring species. In each case, results of linear ML models were compared with kernel based non-linear models. Our results indicate that ring species data does not improve the prediction of chain species. Similarly, chain species data does not improve the prediction of ring species. The use of non-linear ML models, however, did help to minimize the prediction errors compared to the linear models. Furthermore, the study also shows that electronic or adsorption energy based metal descriptors along with bond count based species fingerprints can achieve a mean absolute error (MAE) of less than 0.2 eV for complex chain molecules when used with an appropriate machine learning model.

Adsorption↗
DOE OSTI2023

Ethylene Hydrogenation Molecular Mechanism on MoC y Nanoparticles

Ethylene hydrogenation catalyzed by MoC y nanoparticles has been studied by means of density functional theory methods and several models. These include MetCar (Mo 8 C 12 ), Nanocube (Mo 14 C 13 ), and Mo 12 C 12 nanoparticles as representatives of experimental MoC y nanostructures. The effect of hydrogen coverage has been studied in detail by considering low-, intermediate-, and high-hydrogen regimes. The calculated enthalpy and energy barriers show that ethylene hydrogenation is feasible on the MetCar, Mo 12 C 12 , and Nanocube but at low, medium, and high hydrogen coverages, respectively. An additional step, related to the H* migration from a Mo to a C site in the nanoparticle, has been found to be the key to establishing the best hydrogenation system. In most cases, the reactions are exothermic, featuring low hydrogenation energy barriers, especially for the Nanocube at high hydrogen coverage. In addition, the calculated adsorption Gibbs free energy shows that, for this system, the C 2 H 4 adsorption is feasible in the 300–400 K temperature range and pressures from 10 –10 to 2 atm. For the hydrogenation steps, calculated transition state theory rates show that the overall process is limited by the first hydrogenation step (C 2 H 4 → C 2 H 5 ) at temperatures of 330–400 K. However, at the lower temperatures of 300–320 K, the reaction rates are comparable for the two steps. The present results indicate that the Mo 14 C 13 Nanocube models of MoC y nanoparticles exhibit appropriate thermodynamic and kinetic features to catalyze ethylene hydrogenation at a high-hydrogen-coverage regime. The present findings provide a basis for understanding the chemistry of active MoC y catalysts, suggest appropriate working conditions for the reaction to proceed, and provide a basis for future experimental studies.

03 NATURAL GAS↗
DOE OSTI2023

First-Principles Study of n -Butane Monomolecular Cracking and Dehydrogenation on Two-Dimensional-Zeolite Model Systems: Reaction Mechanisms and Effects of Spatial Confinement

Two-dimensional (2D) ultrathin (~0.5 nm) aluminosilicate bilayer films, consisting of hexagonal prisms (a.k.a. double 6-membered rings D6R) with acidic bridging hydroxyl groups exposed on the surface, have been previously synthesized on a Ru(0001) surface as a zeolite model system. These structures are helpful for mimicking zeolite catalysts with D6R building blocks, such as chabazite. We performed density functional theory calculations to investigate the monomolecular cracking and dehydrogenation of n-butane molecules over the acidic hydroxyl groups of the 2D model system and compared the reaction energetics with that in bulk chabazite. The intrinsic activation energy barrier is the highest for dehydrogenation and lowest for central C–C bond cracking in bulk chabazite. The trend of intrinsic energy barriers for dehydrogenation and terminal and central C–C bond cracking is reproduced on the 2D aluminosilicate film. Overall, the activation barriers are higher on the 2D film than in bulk chabazite due to the lack of confinement in the former. We further explored the effects of the zeolite channel size on the n-butane adsorption and monomolecular cracking using different bulk nanoporous zeolite frameworks (TON, MEL, MEI, and VFI). We found that as the confinement of channels decreases, n-butane adsorption becomes weaker, and the intrinsic energy barrier of terminal C–C cracking increases. The activation energy barriers (dehydrogenation and terminal and central C–C cracking) on the 2D bilayer film surface, which may be considered as zeolite cages at the infinite cage size limit, are close to that in VFI with a relatively large channel size. Comparing the reaction pathway of n-butane terminal C–C cracking in 3D nanocages and on the surface of the 2D aluminosilicate film revealed that stabilizing the transition states in the 3D nanocages is responsible for the decrease in the intrinsic energy barriers for bulk zeolites.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Toward Understanding and Controlling Organic Reactions on Metal Oxide Catalysts

Metal oxides have structurally complex surfaces on which a variety of adsorption site types can occur, including cation sites, anion sites, oxygen vacancy sites, and Brønsted acid sites. These sites can catalyze the catalytic transformation of organic molecules via diverse routes, thus enabling H abstraction, O abstraction, C–C bond formation, and other reactions. This Perspective provides an update on recent advances and future directions for various organic reactions on metal oxide catalyst surfaces, particularly for C–H activation of alkanes and for C–C bond formation with organic oxygenate reactants. Here, we put emphasis on the molecular scale details, on the active site structures required to enable the formation of kinetically relevant transition states, energetic descriptors, as well as contemporary ideas to enable low activation energies. This progress has been enabled by specialized experiments and the increased capabilities of modern electronic structure calculations.

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

Single-Step Selective Oxidation of Methane by Iron-Oxo Species in the Metal–Organic Framework MFU-4l

The direct and selective conversion of methane to methanol can be considered a holy grail for catalysis research. In this work, we study a metal-organic framework known as MFU-4l, modified by design to include highly reactive iron-oxo species for the catalytic C-H bond activation of methane. We investigate the oxidation of methane and the further potential oxidation of the product methanol using N 2 O as an oxygen source and map the potential energy landscape of these reactions using density functional theory calculations. We show that the highest energy barrier encountered during the methane oxidation process is not the C-H bond breaking, but the activation of the iron center by N 2 O. Furthermore, the potential energy landscape for the C-H bond activation exhibits a large, high-energy plateau region instead of a sharp transition state, thus differing from the traditional radical rebound mechanism. Furthermore, this insight offers interesting potential routes to enhance the catalytic activity of the catalyst, to hinder unwanted deactivation pathways, and to reduce the activity towards the over-oxidation of the product.

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

Active Site Engineering via Optimizing the Heterogeneous Support Structure for Single-Atom Catalysis

Supported single-atom catalysts show a large range of activities and selectivities that depend on the local environment of the catalytic sites. Here a theory-based optimization strategy is presented that is based on a density functional theory determination of the transition states and intermediates for a low-dimensional coordinate representation of the heterogeneity of the active sites. The approach is applied to a vanadium catalyst on an amorphous SiO 2 support that involves a large kinetic network described using a full chemistry model. Without assuming a priori scaling relations or mechanism reduction, the optimal state of heterogeneity is found to lie at atomic configurations where the activation energies for two distinct key chemical processes are equal. It is found a posteriori that the behavior of the system is consistent with linear free energy scaling relations in the randomness parameters. The energetic span theory proves quite useful in reducing the full chemistry model to a small number of key reactions. The use of a nonlinear optimization algorithm in combination with energetic span theory provides significant simplification in treating disordered systems.

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

Direct Deoxygenation of Phenol over Fe-Based Bimetallic Surfaces Using On-the-Fly Surrogate Models

We present an accelerated nudged elastic band (NEB) study of phenol direct deoxygenation (DDO) on Fe-based bimetallic surfaces using a recently developed Gaussian process regression (GPR) calculator. Our test calculations demonstrate that the GPR calculator achieves up to 3 times speedup compared to conventional density functional theory calculations while maintaining high accuracy, with energy barrier errors below 0.015 eV. Using GPR-NEB, we systematically examine the DDO mechanism on pure Fe(110) and surfaces modified with Co and Ni in both top and subsurface layers. Our results show that subsurface Co and Ni substitutions preserve favorable thermodynamics and kinetics for both C–O bond cleavage and C–H bond formation, comparable to those on the pure Fe(110) surface. In contrast, top-layer substitutions generally increase the C–O bond cleavage barrier, render the step endothermic, and result in significantly higher reverse reaction rates, making DDO unfavorable on these surfaces. This work demonstrates the effectiveness of GRR-accelerated transition state searches for complex surface reactions and provides insights into rational design of bimetallic catalysts for selective deoxygenation.

Aromatic compounds↗
DOE OSTI2026

Diffusion and Adsorption of 2-Methylpentane and 3-Methylpentane in MFI-Type Zeolite Crystals

Due to the intimate contact of guest molecules with the inner surface of microporous adsorbents, slight changes in the structure of the guest molecules may lead to remarkable changes in their microdynamics. As an impressive example, in this work, diffusion measurements via IR microimaging (IRM) show that the transport diffusivity of 2-methylpentane (2MP) exceeds the value attained for 3-methylpentane (3MP) by a factor of 2−4.5 in silicalite- 1, an MFI-type zeolite crystal and a technically relevant porous host system. This experimental finding is rationalized by comparison with the outcome of dynamically corrected transition state theory (dcTST) simulations.

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

Kinetics and Mechanism of the Singlet Oxygen Atom Reaction with Dimethyl Ether

Here, we combine in situ laser spectroscopy, quantum chemistry, and kinetic calculations to study the reaction of a singlet oxygen atom with dimethyl ether. Infrared laser absorption spectroscopy and Faraday rotation spectroscopy are used for the detection and quantification of the reaction products OH, H 2 O, HO 2 , and CH 2 O on submillisecond time scales. Fitting temporal profiles of products with simulations using an in-house reaction mechanism allows product branching to be quantified at 30, 60, and 150 Torr. The experimentally determined product branching agrees well with master equation calculations based on electronic structure data and transition state theory. The calculations demonstrate that the dimethyl peroxide (CH 3 OOCH 3 ) generated via O-insertion into the C–O bond undergoes subsequent dissociation to CH 3 O + CH 3 O through energetically favored reactions without an intrinsic barrier. This O-insertion mechanism can be important for understanding the fate of biofuels leaking into the atmosphere and for plasma-based biofuel processing technologies.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗