Search NASA⌕ Search

SEARCH · Search NASA

Results for “transition state theory”

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

Model associated with: "Thermodynamic control on the decomposition of organic matter across different electron acceptors"

This model data package is associated with the publication “Thermodynamic control on the decomposition of organic matter across different electron acceptors” submitted to Soil Biology and Biochemistry (Zheng et al., 2023; https://doi.org/10.1016/j.soilbio.2024.109364).In this research, a thermodynamic modeling framework is built to flexibly incorporate both organic matter (OM) molecules and electron acceptors for estimating potential free energy release from various redox reactions and to further predict reaction rates based on Microbial Transition State Theory. The model package includes scripts for thermodynamic modeling and postprocessing. Input Fourier-transform ion cyclotron resonance (FTICR) data are from a previous experimental study (Boye et al., 2018), and model outputs are free energy predictions and stoichiometric coefficients associated with all possible redox reactions.This data package is associated with the project GitHub repository found at MM_bioenergetic_modeling.This data package contains four folders (Input_FTICR, Model, Output, and Output_processing), a file-level metadata (FLMD) csv, and a data dictionary (dd) csv. Please see Zheng_bioenergetic_modeling_flmd.csv for a list of all files contained in this data package and descriptions for each. The Zheng_bioenergetic_modeling_dd.csv file describes the csv column headers. The “Model” folder contains scripts to run energy balance calculations for each electron acceptor. The “Output” folder contains csv files with stoichiometric information from model simulations. And the "Output_processing" folder contains scripts for reaction rate calculations and to generate plots.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1973

Estimation of rate constants of elementary processes - A review of the state of the art.

'Thermochemical kinetics,' the codification and extrapolation of empirical observations, as applied to certain elementary reactions of importance to combustion studies, is described. This approach allows the critical scrutiny of experimental data in areas where sufficient previous data exist, while, at the same time, illuminating those key areas where more experimentation is crucial. It is shown that combination of transition-state theory with an understanding of the molecular basis of entropy puts fairly rigid constraints on the values of the Arrhenius A-factor for most reactions. This, in turn, means that the activation energy is often the key datum that is missing, and that such data can be obtained with some degree of confidence, even from measurements of rate constants at only one temperature. In complex mechanisms, it is often possible to distinguish among alternate pathways and pinpoint key processes.

Golden, D. M.↗
NASA NTRS1980

A molecular beam investigation of the oxidation of CO on Pt/9/111/x/100//

The CO oxidation on Pt/9(111)x(100)/ was studied by molecular beam relaxation spectroscopy (MBRS). The reaction proceeded via the reaction of adsorbed CO and adsorbed oxygen. No evidence for direct reactive collisions between gaseous CO and adsorbed atomic oxygen was seen. The second order rate constant was measured by linearizing the reaction system to be 10 to the -7th exp/-(9700 kcal/mole)RT/ per sq cm s. The relatively low pre-exponential factor was explained by transition state theory on the basis of a high partition function for adsorbed carbon monoxide obtained previously from studies of CO desorption on this surface.

Fair, J. A.↗
NASA NTRS1986

Kinetic study of the reaction CH (X 2Pi) + H2 yields CH2 (X 3B1) + H in the temperature range 372 to 675 K

The kinetics of the reversible reaction CH (X 2Pi) + H2 yields CH2 (X 3B1) + H at 372-675 K and total pressure 100 torr (mainly Ar) is investigated experimentally. The ground-state CH radicals are produced by photolysis of CHBr3 using 10-mJ 266-nm laser pulses (repetition rate 10 Hz) and monitored by measuring the fluorescence induced by a 429.8-nm dye laser, in the apparatus described by Berman et al. (1982) and Berman and Lin (1984). The results are presented in tables and graphs and characterized. The absolute rate constants for the forward and reverse reactions are determined, and their temperature dependence is given by Arrhenius expressions and formulas obtained in transition-state-theory calculations. The heat of formation of CH2 at 0 K is estimated (assuming that the recombination reaction CH2 + H has zero activation energy) as 92.6 + or - 0.5 kcal/mol.

Zabarnick, S.↗
NASA NTRS1988

Computed potential energy surfaces for chemical reactions

The minimum energy path for the addition of a hydrogen atom to N2 is characterized in CASSCF/CCI calculations using the (4s3p2d1f/3s2p1d) basis set, with additional single point calculations at the stationary points of the potential energy surface using the (5s4p3d2f/4s3p2d) basis set. These calculations represent the most extensive set of ab initio calculations completed to date, yielding a zero point corrected barrier for HN2 dissociation of approx. 8.5 kcal mol/1. The lifetime of the HN2 species is estimated from the calculated geometries and energetics using both conventional Transition State Theory and a method which utilizes an Eckart barrier to compute one dimensional quantum mechanical tunneling effects. It is concluded that the lifetime of the HN2 species is very short, greatly limiting its role in both termolecular recombination reactions and combustion processes.

Walch, Stephen P.↗
NASA NTRS1990

Theoretical characterization of the 5Pi and 3Pi potential energy surfaces for NH + O yields N + OH

The reactant, product, and saddle point regions of the 5Pi and 3Pi potential energy surfaces for the reaction NH + O yields N + OH have been characterized using complete active space self consistent field/externally contracted configuration interaction calculations with large atomic natural orbital basis sets. The computed barrier heights are 5.6 and 11.7 kcal/mol on the 5Pi and 3Pi surfaces, respectively. Transition state theory with an Eckart tunneling correction is used to estimate the rate constant on the 5Pi surface.

Walch, Stephen P.↗
NASA NTRS1993

Isomerization reaction dynamics and equilibrium at the liquid-vapor interface of water. A molecular-dynamics study

The gauche-trans isomerization reaction of 1,2-dichloroethane at the liquid-vapor interface of water is studied using molecular-dynamics computer simulations. The solvent bulk and surface effects on the torsional potential of mean force and on barrier recrossing dynamics are computed. The isomerization reaction involves a large change in the electric dipole moment, and as a result the trans/gauche ratio is considerably affected by the transition from the bulk solvent to the surface. Reactive flux correlation function calculations of the reaction rate reveal that deviation from the transition-state theory due to barrier recrossing is greater at the surface than in the bulk water. This suggests that the system exhibits non-Rice-Ramsperger-Kassel-Marcus behavior due to the weak solvent-solute coupling at the water liquid-vapor interface.

Benjamin, Ilan↗
NASA NTRS2002

Tensile Strength of Carbon Nanotubes Under Realistic Temperature and Strain Rate

Strain rate and temperature dependence of the tensile strength of single-wall carbon nanotubes has been investigated with molecular dynamics simulations. The tensile failure or yield strain is found to be strongly dependent on the temperature and strain rate. A transition state theory based predictive model is developed for the tensile failure of nanotubes. Based on the parameters fitted from high-strain rate and temperature dependent molecular dynamics simulations, the model predicts that a defect free micrometer long single-wall nanotube at 300 K, stretched with a strain rate of 1%/hour, fails at about 9 plus or minus 1% tensile strain. This is in good agreement with recent experimental findings.

Wei, Chen-Yu↗
NASA NTRS2001

Reactions of SiCl2 and SiHCl with H and Cl Atoms

Calculations have been carried out for the reaction of SiCl2 and SiHCl with H and Cl atoms. In each case, the stationary point geometries and harmonic frequencies were characterized using CASSCF/derivative methods and the cc-pVDZ basis set. Accurate energetics were obtained by combining the CCSD(T) results using the a-cc-pVTZ basis set with an extrapolation to the basis set limit using the a-cc-pVDZ, a-cc-pVTZ, and a-cc-pVQZ basis sets at the MP2 level. The geometries, energetics, and harmonic frequencies were used to obtain rate constants using conventional transition state theory or a Gorin-like model. In each case we find direct abstraction pathways compete with an addition elimination pathway. In the case of SiClH + H the two direct pathways are H abstraction which is barrierless and Cl abstraction with a barrier of 13.5 kcal/mol, while the addition elimination process has a barrier of 26.9 kcal/mol. In the case of SiCl2 + H the direct pathway is Cl abstraction with a barrier of 16.4 kcal/mol, while the addition elimination pathway has a barrier of 29.6 kcal/mol. In the case of SiClH + Cl the direct pathway is H abstraction which is barrierless and the addition elimination pathway has a barrier of 2.0 kcal/mol.

Walch, Stephen P.↗
NASA NTRS2002

Tensile Yielding of Multi-Wall Carbon Nanotube

The tensile yielding of multiwall carbon nanotubes (MWCNTs) has been studied using Molecular Dynamics simulations and a Transition State Theory based model. We find a strong dependence of the yielding on the strain rate. A critical strain rate has been predicted above/below which yielding strain of a MWCNT is larger/smaller than that of the corresponding single-wall carbon nanotubes. At experimentally feasible strain rate of 1% /hour and T = 300K, the yield strain of a MWCNT is estimated to be about 3-4 % higher than that of an equivalent SWCNT (Single Wall Carbon Nanotube), in good agreement with recent experimental observations.

Wei, Chenyu↗
NASA NTRS2008

Synchrotron Photoionization Mass Spectrometry Measurements of Kinetics and Product Formation in the Allyl Radical (H2CCHCH2)Self Reaction

Product channels for the self-reaction of the resonance-stabilized allyl radical, C3H5 + C3H5, have been studied with isomeric specificity at temperatures from 300-600 K and pressures from 1-6 Torr using time-resolved multiplexed photoionization mass spectrometry. Under these conditions 1,5-hexadiene was the only C6H10 product isomer detected. The lack of isomerization of the C6H10 product is in marked contrast to the C6H6 product in the related C3H3 + C3H3 reaction, and is due to the more saturated electronic structure of the C6H10 system. The disproportionation product channel, yielding allene + propene, was also detected, with an upper limit on the branching fraction relative to recombination of 0.03. Analysis of the allyl radical decay at 298 K yielded a total rate coefficient of (2.7 +/- 0.8) x 10(exp -11) cu cm/molecule/s, in good agreement with pre.vious experimental measurements using ultraviolet kinetic absorption spectroscopy and a recent theoretical determination using variable reaction coordinate transition state theory. This result provides independent indirect support for the literature value of the allyl radical ultraviolet absorption cross-section near 223 nm.

Selby, Talitha M.↗
NASA NTRS2012

Isotopic Ratios in Titan's Methane: Measurements and Modeling

The existence of methane in Titan's atmosphere (approx. 6% level at the surface) presents a unique enigma, as photochemical models predict that the current inventory will be entirely depleted by photochemistry in a timescale of approx 20 Myr. In this paper, we examine the clues available from isotopic ratios (C-12/C-13 and D/H) in Titan's methane as to the past atmosphere history of this species. We first analyze recent infrared spectra of CH4 collected by the Cassini Composite Infrared Spectrometer, measuring simultaneously for the first time the abundances of all three detected minor isotopologues: (13)CH4, (12)CH3D, and (13)CH3D. From these we compute estimates of C-12/C-13 = 86.5 +/- 8.2 and D/H = (1.59 +/- 0.33) x 10(exp -4) , in agreement with recent results from the Huygens GCMS and Cassini INMS instruments. We also use the transition state theory to estimate the fractionation that occurs in carbon and hydrogen during a critical reaction that plays a key role in the chemical depletion of Titan's methane: CH4 + C2H yields CH3 + C2H2. Using these new measurements and predictions we proceed to model the time evolution of C-12/C-13 and D/H in Titan's methane under several prototypical replenishment scenarios. In our Model 1 (no resupply of CH4), we find that the present-day C-12/C-13 implies that the CH4 entered the atmosphere 60-1600 Myr ago if methane is depleted by chemistry and photolysis alone, but much more recently-most likely less than 10 Myr ago-if hydrodynamic escape is also occurring. On the other hand, if methane has been continuously supplied at the replenishment rate then the isotopic ratios provide no constraints, and likewise for the case where atmospheric methane is increasing, We conclude by discussing how these findings may be combined with other evidence to constrain the overall history of the atmospheric methane.

Nixon, C. A.↗
DOE OSTI2021

Scaling of Transition State Vibrational Frequencies and Application of d -Band Theory to the Brønsted–Evans–Polanyi Relationship on Surfaces

Semiempirical energy relations provide a means of estimating thermodynamic properties. Specifically, linear scaling relationships (LSRs) and Brønsted–Evans–Polanyi (BEP) relationships correlate adsorption energies between adsorbates across surfaces and reaction energies with activation barriers, respectively. Although vibrational scaling relations (VSRs) exist between adsorbates at identical sites, scaling between vibrational frequencies of adsorbed local minima and transition states is lacking. Here, we present density functional theory calculations for AH X (A = C, N, O) diffusions on transition metal surfaces and reveal linear scaling between frequencies of local minima and the transition state between those minima. Using d-band theory and linear muffin tin orbital theory (LMTO), we derive the slopes of these transition state vibrational scaling relations (TSVSRs) and, in so doing, provide a rigorous theory extending the original BEP relations developed for solution chemistry to surface chemistry. Furthermore, with a single reference DFT calculation, we predict the slopes and quantify uncertainty in the predictions.

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

Nucleation Rate Theory for Coordination Number: Elucidating Water-Mediated Formation of a Zigzag Na 2 SO 4 Morphology

Predicting and controlling nanostructure formation during nucleation can pave the way to synthesizing novel energy materials via crystallization. However, such control over nucleation and crystallization remains challenging due to an inadequate understanding of critical factors that govern evolving atomistic structures and dynamics. Here, in this paper, we utilize coordination number as a reaction coordinate and rate theory to investigate how sodium sulfate, commonly known as a phase-change energy material, nucleates in a supersaturated aqueous solution. In conjunction with ab initio and force field-based molecular dynamics simulation, the rate theoretical analysis reveals that sodium sulfate from an initially dissolved metastable state transits to a heterogeneous mixture of prenucleated clusters and finally to a large cylindrical zigzag morphology. Measurements of Raman spectra and their ab initio modeling confirm that this nucleated morphology contains a few waters for every sulfate. Rate processes such as solvent exchange and desolvation exhibit high sensitivity to the evolving prenucleation/nucleation structures, providing a means to distinguish between critical nucleation precursors. Desolvation and forming the first-shell interionic coordination structure via monomer-by-monomer addition around sulfates are found to explain the formation of large nuclei. Thus, a detailed understanding of the step-by-step structure formation across scales has been achieved. This can be leveraged to predict nucleation-related structures and dynamics and potentially control the synthesis of novel phase-change materials for energy applications.

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

Supervised learning and the finite-temperature string method for computing committor functions and reaction rates

A central object in the computational studies of rare events is the committor function. Though costly to compute, the committor function encodes complete mechanistic information of the processes involving rare events, including reaction rates and transition-state ensembles. Under the framework of transition path theory, Rotskoff et al. [Proceedings of the 2nd Mathematical and Scientific Machine Learning Conference, Proceedings of Machine Learning Research (PLMR, 2022), Vol. 145, pp. 757–780] proposes an algorithm where a feedback loop couples a neural network that models the committor function with importance sampling, mainly umbrella sampling, which collects data needed for adaptive training. Here, in this work, we show additional modifications are needed to improve the accuracy of the algorithm. The first modification adds elements of supervised learning, which allows the neural network to improve its prediction by fitting to sample-mean estimates of committor values obtained from short molecular dynamics trajectories. The second modification replaces the committor-based umbrella sampling with the finite-temperature string (FTS) method, which enables homogeneous sampling in regions where transition pathways are located. We test our modifications on low-dimensional systems with non-convex potential energy where reference solutions can be found via analytical or finite element methods, and show how combining supervised learning and the FTS method yields accurate computation of committor functions and reaction rates. We also provide an error analysis for algorithms that use the FTS method, using which reaction rates can be accurately estimated during training with a small number of samples. The methods are then applied to a molecular system in which no reference solution is known, where accurate computations of committor functions and reaction rates can still be obtained.

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

High-Pressure Rate Rules for Ether Alkylperoxy Radical Isomerization

The first isomerization reaction of an alkylperoxy (RO 2 ) radical holds significant importance in low-temperature oxidation, as it governs the branching ratios of the hydroperoxyalkyl (QOOH) radicals, which influence the competition between the chain-propagation and chain-branching reactions. In this study, we systematically calculated high-pressure rate rules for the RO 2 isomerization reaction of monoethers, exploring 5-, 6-, 7-, and 8-membered ring transition states. Primary, secondary, and tertiary carbon sites, where both the abstracting peroxy group and the abstracted hydrogen are located, were considered, with particular emphasis on distinguishing between secondary carbons adjacent (alpha) and nonadjacent to the ether functional group. Using the G4//B3LYP/6-311++G(2df,2pd) level of theory and the transition state theory, we estimated the rate constants and the Arrhenius coefficient for over 120 possible isomerization reactions. We examined the effect of ring size and ring atoms, revealing that 6- and 7-membered ring isomerizations were generally the fastest. The impact of the ether functional group on transition states was investigated by comparing reactions with identical ring size, peroxy, and radical positions, but with the ether functional group positioned either outside (i.e., out) or inside (i.e., in) the transition state ring, leading to differences in the rate constants. When comparing to analogous alkane rate constants, differences of up to an order of magnitude were observed, underscoring the need for caution when assigning rate rules by analogy. We applied our rate constants in the di-iso-butyl ether kinetic model and evaluated their influence on low-temperature chemistry finding that they altered the branching ratios by up to a factor of 9, highlighting the significance of site-specific rate constants for more accurate low-temperature modeling.

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

Density functional theory-based surrogate kinetic models for heterogeneous reactions of hydrocarbon intermediates on silicon carbide

The increasing demand for high-performance materials in advanced technologies highlights the importance of achieving a fundamental understanding and potential control of silicon carbide (SiC) deposition processes. However, existing models often lack sufficient theoretical detail, relying heavily on empirical data and offering limited predictive capability. In particular, the complex surface chemistry governing SiC growth remains poorly understood. This study addresses these challenges by employing density functional theory (DFT) to investigate key heterogeneous reactions involving hydrocarbon intermediates on SiC surfaces, including dehydrogenation, hydrogenation, and carbon deposition. Transition state searches were conducted to identify reaction pathways and energy barriers. While first-principles calculations offer high accuracy, they are computationally intensive. To extend the utility of these first-principles results, vibrational analyses were performed using phonon-based statistical thermochemistry to compute temperature-dependent reaction rates which were used to develop Arrhenius-type surrogate kinetic models. Furthermore, the resulting framework provides a more rigorous, physically grounded basis for integrating atomistic insights into continuum-scale modeling, ultimately enabling improved prediction and optimization of SiC film growth in high-performance material systems.

Density Functional Theory↗