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 433 records · Page 24
All sourcesNASA NTRSDOE OSTI
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↗
DOE OSTI2025

Kinetic Isotope Effect in the Radical Step-Growth Termination of On-Surface Synthesized Graphene Nanoribbons

Bottom-up on-surface synthesis has emerged as a versatile tool to access and finely tune the electronic structure of nanographenes. The controlled generation of reactive intermediates catalyzed and stabilized by a supporting substrate has enabled the design, assembly, and characterization of a wide range of exotic tailor-made quantum materials. Even under these tightly controlled conditions, the growth of extended structures remains limited by termination processes and undesired side reactions. Here, we identify an H atom transfer as one principal contributor to the radical step growth termination that limits the on-surface growth of N = 7 armchair graphene nanoribbons (7-AGNRs) on Au(111) surfaces. Analysis of 7-AGNR lengths grown from protiated and deuterated molecular precursors reveals a primary kinetic isotope effect of KIE ∼ 1.4. First-principles density functional theory calculations suggest that a concerted H atom transfer mechanism that involves the breaking of a C–H/D bond in the transition state is associated with radical chain termination.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2023

Synergic Effects of Surface Chemistry and Applied Potentials on the Kinetics of the Electrocatalytic Oxygen Evolution Reaction in IrO 2

Improving the efficacy of the oxygen evolution reaction (OER) through water oxidation is critical for advancing photoelectrochemical water splitting. Among the catalysts, IrO 2 exhibits a high OER catalytic activity and stability under acidic conditions. The OER mechanism in this system has been a topic of intense research; however, many mechanistic understandings are lacking. In this work, we carried out first-principles calculations with an implicit solvation model at a constant potential to study the OER behavior on the IrO 2 (110) surface. We find that the surface hydrogen coverage has a significant effect on OER kinetics and transition states. We then develop a microkinetics model that accounts for the continuous evolution of both OER activation energy and hydrogen coverage as a function of an applied potential. We show that this inclusion leads to significant improvement in the simulated Tafel plot compared to available experiments. In conclusion, our results point to a complex interplay between surface chemistry and the applied potential on OER kinetics.

08 HYDROGEN↗
DOE OSTI2021

Reactive Species and Reaction Pathways for the Oxidative Cleavage of 4-Octene and Oleic Acid with H 2 O 2 over Tungsten Oxide Catalysts

The oxidative cleavage of carbon-carbon double bonds (C=C) in alkenes and fatty acids produces aldehydes and acids valued as chemical intermediates. Solid tungsten oxide catalysts are low cost, non-toxic, and selective for the oxidative cleavage of C=C bonds with hydrogen peroxide (H 2 O 2 ) and are, therefore, a promising option for continuous processes. Despite the relevance of these materials, the elementary steps involved and their sensitivity to the form of W sites present on surfaces have not been described. In this work, we combine in situ spectroscopy and rate measurements to identify significant steps in the reaction, the reactive species present on catalysts, and examine differences between the kinetics of this reaction on isolated W atoms grafted to alumina and on those exposed on crystalline WO 3 nanoparticles. Raman spectroscopy shows that W-peroxo complexes (W-(η 2 -O 2 ) formed from H 2 O 2 ) react in a kinetically relevant step to produce epoxides, which undergo hydrolysis at protic surface sites. Subsequently, the CH 3 CN solvent deprotonates diols to form alpha hydroxy ketones that react to form aldehydes and water following nucleophilic attack of H 2 O 2 . Turnover rates for oxidative cleavage, determined by in situ site titrations, on WO x -Al 2 O 3 are 75% greater than those on WO 3 at standard conditions. These differences reflect activation enthalpies (ΔH ‡ ) for oxidative cleavage of 4-octene that are much lower for the isolated WO x sites (36±3 and 60±6 kJ∙mol -1 for WO x -Al 2 O 3 and WO 3 , respectively) and correlate strongly with the difference between enthalpies of adsorption for epoxyoctane (ΔH ads,epox ), which resembles the transition state for epoxidation. The WO x -Al 2 O 3 catalysts mediate oxidative cleavage of oleic acid with H 2 O 2 following a mechanism comparable to that for 4-octene. The WO 3 materials, however, form only the epoxide and do not cleave the C-C bond or produce aldehydes and acids. These differences reflect the distinct site requirements for these reaction pathways and indicate that acid sites required for diol formation are strongly inhibited by oleic acids and epoxides on WO 3 whereas the Al 2 O 3 support provides sites competent for this reaction and increase the yield of the oxidative cleavage products.

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

Optimizing the Solvent Reorganization Free Energy by Metal Substitution for Nanocage Catalysis

The supramolecular capsule Ga 4 L 6 12– has been found to catalyze a number of important chemical reactions, but the close proximity of a poorly organized external solvent environment can serve as an impediment to the reaction chemistry. Here, using ab initio molecular dynamic calculations and electric field analyses, we find that metal substitution of indium for gallium lowers the total and electrostatic activation barriers by ~3–4 kcal/mol in water solvent. Using an energy decomposition analysis of the interaction of water solvent with the metal vertices, we find that Pauli repulsion between the metal and water decreases, allowing the water coordination with the metal to increase, upon substitution of In for Ga. We therefore determine that the stabilization of the transition state is due to a better arrangement of water molecules around the metal vertices of In 4 L 6 12– and provides a proof of concept of how to redesign the nanocage scaffold in order to reduce the solvent reorganization energy.

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

Modeling the Effect of Surface Platinum–Tin Alloys on Propane Dehydrogenation on Platinum–Tin Catalysts

Uncertainty analysis, reported experimental literature data, and density functional theory were synthesized to model the effect of surface tin coverage on platinum-based catalysts for nonoxidative propane dehydrogenation to propylene. Here, this study tests four different platinum–tin skin surface models as potential catalytic sites, Pt 3 Sn/Pt(100), PtSn/Pt(100), Pt 3 Sn/Pt(111), and Pt 2 Sn/Pt(211), and compares them to the corresponding pure Pt surface sites using an uncertainty analysis methodology that uses BEEF-vdW with its ensembles (BMwE) to generate the uncertainty for the energies of the intermediates and transition states. One experimental data set with two experimental observations, selectivity to propylene and turnover frequency of propylene, was used as a calibration data set to evaluate the impact of the experimental data on informing the models. This study finds that the prior model for Pt 3 Sn/Pt(100) is the most active and Pt 2 Sn/Pt(211) is the most selective toward propylene. Active sites on the (100) facet have the highest probability of being responsible for C 1 and C 2 product formations (C–C bond cleavage). Increasing the Sn coverage on the (100) surface facet to a PtSn/Pt(100) active site leads to a significantly reduced rate and might explain the experimentally observed higher selectivity of Sn-doped catalysts relative to pure Pt catalysts. Next, this study finds that for all surfaces, except PtSn/Pt(100), the rate-controlling steps are the initial dehydrogenation steps alongside some partially rate-controlling second dehydrogenation steps. For PtSn/Pt(100), only the initial terminal dehydrogenation step to CH 3 CH 2 CH 2 * and second dehydrogenation steps are rate-controlling. Next, the calibrated models for all surfaces were found to be selective toward propylene production and model the reported turnover frequency successfully. Nevertheless, Pt 2 Sn/Pt(211) emerges as the active site with some (minor) evidence as the main active site based on Jeffreys’ scale interpretation of Bayes factors. This observation agrees with prior studies that also found step sites to be most likely the most relevant active sites for pure Pt catalysts. Overall, the results indicate that tin, in addition to affecting the binding strength of the adsorbed species, prevents deeper dehydrogenation (reducing coking) and cracking reactions through increasing activation barriers for unwanted side reactions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗