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 19 records

Experimental and theoretical examples of the value and limitations of transition state theory

Value and limitations of transition-state theory (TST) are reviewed. TST analyses of the temperature dependence of the 'direct' reactions CH3 + CH3CHO yields CH4 + CH3CO(1) and O + CH4 yields OH + CH3(2) are presented in detail, and other examples of TST usefulness are recalled. Limitations are discussed for bimolecular processes in terms of 'complex' vs. 'direct' mechanisms. The reaction OH + CO yields CO2 + H is discussed in this context. Limitations for unimolecular processes seem to arise only for simple bond fission processes, and recent advances are noted.

Golden, D. M.↗

Semiclassical transition state theory through the lens of the restricted partition function

The wide adoption of transition state theory in resolving the rates of molecular processes relies on the simplification from reducing the formal and numerical expense of dynamics by a geometric constraint. Such a reduction is at odds with the uncertainty in localization that the uncertainty principle requires. While many forms of semiclassical transition state theory (SCTST) have been aimed at addressing this challenge, a popular approach has relied on resolving the underlying phase space structure of the exact rate formula to leverage Bohr–Sommerfeld quantization. Here, the Hernandez–Miller SCTST reframed the thermal rate formula into an integral of the so-called restricted partition function (RPF) over the action associated with the reaction. The density-of-state SCTST has reframed the rate formula in terms of the instanton’s density of states (DoS). Here, we show the relationship between the RPF-SCTST and the DoS-SCTST and derive the latter from the former. In this way, we help unify these branches of SCTST and provide a clearer formalism for future advances.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Microscopic Derivation of Transition-state Theory for Complex Quantum Systems

The decay of quantum complex systems through a potential barrier is often described with transition-state theory, also known as RRKM theory in chemistry. Here we derive the basic formula for transition-state theory based on a generic Hamiltonian as might be constructed in a configuration-interaction basis. Two reservoirs of random Hamiltonians from Gaussian orthogonal ensembles are coupled to intermediate states representing the transition states at a barrier. Under the condition that the decay of the reservoirs to open channels is large, an analytic formula for reaction rates is derived. Here, the transition states act as independent Breit–Wigner resonances which contribute additively to the total transition probability, as is well known for electronic conductance through resonant tunneling states. It is also found that the transition probability is independent of the decay properties of the states in the second reservoir over a wide range of decay widths.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Restricted Partition Function Semiclassical Transition State Theory (RPF-SCTST): Applications to Reactions Involving Hydrogen Cyanide, Hydrogen Peroxide, and Formaldehyde Oxide

We demonstrate the efficiency of the numerical calculation of thermal semiclassical transition state theory (SCTST) rates across several representative chemical reactions using the restricted partition function (RPF)─viz a function that depends only on the imaginary action associated with a reaction. Here, we treat the potential energy surfaces (PESs) through fourth order expansions around the saddle point and integrate the Hamiltonian up to second order in employing vibrational perturbation theory. We apply this formalism to uni- and bimolecular reactions with different─though relatively high─barrier heights to uncover the influence of the barrier properties, minimal energies, and separability of the rotational motion on rate constants and tunneling corrections. Although all modes are coupled within the RPF, we found in our numerical examples that the rotational component in the absence of significant rotational distortions can be separated from the vibrational contributions. Moreover, a classical treatment of the rotational contribution is adequate over the temperature range from approximately 100 to 1000 K. We also found that the choice of DFT basis set can lead to variations in rate constants of up to an order of magnitude. Lastly, different schemes for counting eligible energy levels result in rate constants that differ by no more than a factor of 2, with the remaining discrepancies attributed to the treatment of near-convergent levels in perturbation theory.

74 ATOMIC AND MOLECULAR PHYSICS↗

Use of scaled external correlation, a double many-body expansion, and variational transition state theory to calibrate a potential energy surface for FH2

New ab initio results and a double many-body expansion formalism have been used to parameterize a new FH2 potential energy surface with improved properties near the saddle point and in the region of long-range attraction. The functional form of the new surface includes dispersion forces by a double many-body expansion. Stationary point properties for the new surface are calculated along with the product-valley barrier maxima of vibrationally adiabatic potential curves for F + H2 - HF(nu-prime = 3) + H, F + HD - HF(nu-prime = 3) + D, and F + D2 - DF(nu-prime = 4) + D. The new surface should prove useful for studying the effect on dynamics of a low, early barrier with a wide, flat bend potential.

Lynch, Gillian C.↗

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↗

Association Kinetics for Perfluorinated n -Alkyl Radicals

Radical-radical reaction channels are important in the pyrolysis and oxidation chemistry of perfluoroalkyl substances (PFAS). In particular, unimolecular dissociation reactions within unbranched n-perfluoroalkyl chains, and their corresponding reverse barrierless association reactions, are expected to be significant contributors to the gas-phase thermal decomposition of families of species such as perfluorinated carboxylic acids and perfluorinated sulfonic acids. Unfortunately, experimental data for these reactions are scarce and uncertain. Furthermore, obtaining reliable theoretical predictions for such reactions is a laborious and computationally intensive task. Here, in this work, the chemical kinetics of the various association/decomposition reactions producing/decomposing the C 2 -C 4 series of unbranched n-perfluoroalkanes (C 2 F 6 , C 3 F 8 , and C 4 F 10 ) are examined using state-of-the-art ab initio transition-state-theory-based master-equation calculations. The variable-reaction-coordinate transition-state theory (VRC-TST) formalism is employed in computing the microcanonical and canonical rates for the association reactions. Reaction thermochemistry is obtained via composite quantum chemistry calculations and the laddering of error-canceling reaction schemes via a connectivity-based hierarchy approach employing ANL1/ANL0-style reference energies. Lennard-Jones collision model parameters for the considered systems were estimated by a direct dynamics approach, and collisional energy transfer parameters were obtained from analogies to systems of similar size and heavy-atom connectivity. A one-dimensional master equation approach was used to convert the microcanonical rate coefficients from the VRC-TST analysis into temperature- and pressure-dependent rate constants for the association reactions and the reverse dissociation reactions. The data are reported in standardized formats for usage in comprehensive chemical kinetic models for PFAS thermal destruction.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

An Ab initio based OH initiated oxidation kinetics of glycerol carbonate: A promising biofuel component

The global energy demand is steadily increasing because of the population explosion and economic growth. Fossil fuels supply around 85 % of global primary energy demand. On one hand, fulfilling the increasing energy demands is a big challenge for the next few decades. On the other hand, the continued burning of fossil fuels leads to higher CO2 emissions, severely impacting global warming. Therefore, the policymakers vow to shift from conventional fuels to renewable resources for economic, environmental, and future energy security reasons. In this context, biofuels from lignocellulosic biomass and/or carbon-neutral fuels produced in the sustainable carbon cycle can close the carbon cycle and reach net zero-carbon emission. Recently, glycerol carbonate has been proposed as a promising fuel or fuel additive for future sustainability. Therefore, we investigated the hydrogen abstraction reactions of glycerol carbonate (GC) by OH radicals using high-level ab initio and variational transition state theory calculations. We mapped out the potential energy surface using the CCSD(T)/cc-pV(D, T)Z//MP2/cc-pVTZ level of theory. Here we used the ab initio parameters to obtain the site-specific rate coefficients by employing the variational transition state theory. We observed that every hydrogen atom in GC displays a unique reactivity with OH radicals. We derived branching ratio of each channel that are difficult to access experimentally. The overall rate coefficients exhibit a strong non-Arrhenius behaviour, which can be represented as: $k^{CVT/SCT}_{ov}$ (T) = 3.39 x 10 -20 x T 2.659 x e$\frac{-750.0 Jmol^{-1}}{RT}$ $\frac{cm^{3}}{molecule s}$ This is the first reported rate data for the glycerol carbonate and OH radicals reaction.

09 BIOMASS FUELS↗

From Coherence to Function: Exploring the Connection in Chemical Systems

The role of quantum mechanical coherences or coherent superposition states in excited state processes has received considerable attention in the last two decades largely due to advancements in ultrafast laser spectroscopy. These coherence effects hold promise for enhancing the efficiency and robustness of functionally relevant processes, even when confronted with energy disorder and environmental fluctuations. Understanding coherence deeply drives us to unravel mechanisms and dynamics controlled by order and synchronization at a quantum mechanical level, envisioning optical control of coherence to enhance functions or create new ones in molecular and material systems. In this frontier, the interplay between electronic and vibrational dynamics, specifically the influence of vibrations in directing electronic dynamics, has emerged as the leading principle. Here, two energetically disparate quantum degrees of freedom work in-sync to dictate the trajectory of an excited state reaction. Moreover, with the vibrational degree being directly related to the structural composition of molecular or material systems, new molecular designs could be inspired by tailoring certain structural elements. In the realm of chemical kinetics, our understanding of the dynamics of chemical transformations is underpinned by fundamental theories, such as transition state theory, activated rate theory, and Marcus theory. These theories elucidate reaction rates by considering the energy barriers that must be overcome for reactants to transform into products. Those barriers are surmounted by the stochastic nature of energy gap fluctuations within reacting systems, emphasizing that the reaction coordinate, the pathway from reactants to products, is not rigidly defined by a specific vibrational motion but encompasses a diverse array of molecular motions. While less is known about the involvement of specific intramolecular vibrational modes, their significance in certain cases cannot be overlooked. In this Account, we summarize key experimental findings that offer deeper insights into the complex electronic–vibrational trajectories encompassing excited states afforded from state-of-the-art ultrafast laser spectroscopy in three exemplary processes: photoinduced electron transfer, singlet–triplet intersystem crossing, and intramolecular vibrational energy flow in molecular systems. We delve into the rapid decoherence, or loss of phase and amplitude correlations, of vibrational coherences along promoter vibrations during subpicosecond intersystem crossing dynamics in a series of binuclear platinum complexes. This rapid decoherence illustrates the vibration-driven reactive pathways from the Franck–Condon state to the curve crossing region. We also explore the generation of new vibrational coherences induced by impulsive reaction dynamics rather than by the laser pulse in these systems, which sheds light on specific energy dissipation pathways and thereby on the progression of the reaction trajectory in the vicinity of the curve crossing on the product side. Another property of vibrational coherences, amplitude, reveals how energy can flow from one vibration to another in the electronic excited state of a terpyridine–molybdenum complex hosting a nonreactive dinitrogen substrate. In conclusion, a slight change in vibrational energy triggers a quasi-resonant interaction, leading to constructive wavepacket interference and ultimately intramolecular vibrational redistribution from a Franck–Condon active terpyridine vibration to a dinitrogen stretching vibration, energizing the dinitrogen bond.

Electrical energy↗

Nuclear–Electronic Orbital General Rate Theory: Predicting Hydrogen Kinetic Isotope Effects in the Deep Tunneling Regime

Hydrogen transfer is a critical component of many chemical and biological processes. The ratio of rate constants for hydrogen and deuterium transfer defines the H/D kinetic isotope effect (KIE), which is a powerful tool for elucidating hydrogen transfer mechanisms. Interpretation of experimental H/D KIEs relies on accurate and affordable computational methods. However, due to their light mass, hydrogen and deuterium can undergo tunneling, which is challenging to describe in multidimensional molecular systems. Herein, we introduce the nuclear–electronic orbital general rate theory (NEO-GRT), which enables the efficient prediction of H/D KIEs based on full-dimensional molecular quantum chemistry calculations. The NEO-GRT approach describes the hydrogen transfer rate constant with a general expression that spans the vibrationally adiabatic and nonadiabatic hydrogen tunneling regimes. The input quantities are computed using NEO density functional theory, which treats the transferring hydrogen or deuterium nucleus quantum mechanically on the same level as the electrons. We investigate two intramolecular proton transfer reactions in organic molecules at temperatures down to 50 K to evaluate the performance of NEO-GRT by comparison to transition state theory and ring-polymer instanton theory. The KIEs computed with NEO-GRT agree with those calculated using ring-polymer instanton theory for the full-dimensional molecular systems at the same level of electronic structure theory. This agreement indicates that NEO-GRT captures the deep hydrogen tunneling effects, in contrast to transition state theory, which neglects such effects. Given its relatively low computational cost, NEO-GRT is a promising approach for predicting H/D KIEs in large organic and organometallic systems.

Hydrogen↗

Kinetic conversion of CO to CH4 in the Solar System

Some of the most interesting chemistry in the Solar System involves changes in the oxidation state of the simple carbon species. The chemical pathways for the conversion of CH4 to CO and CO2 are for the most part known. The reverse process, the reduction of CO to CH4, is, however, poorly understood. This is surprising in view of the importance of the reduction process in the chemistry of the Solar System. Recently we investigated the chemical kinetics of a hitherto unsuspected reaction. It is argued that the formation of the methoxy radical (CH3O) from H+H2CO may play an essential role in the reduction of CO to CH4. The rate coefficient for this reaction has been estimated using the approximate theory of J. Troe and transition state theory. We will discuss the implications of this reaction for the chemistry of CO on Jupiter, in the solar nebula, for interpreting the laboratory experiments of A. Bar-Nun and A. Shaviv and A. Bar-Nun and S. Chang, and for organic synthesis in the prebiotic terrestrial atmosphere. The possible relation of CO reduction in the solar nebula and polyoxymethylene observed in comet Halley will be discussed.

Yung, Y. L.↗

A theoretical kinetic study of ĊH 3 + ṄH 2 : From electronic structure to NH 3 /CH 4 combustion modelling implications

Carbon–nitrogen interaction reactions play an important role in governing the reactivity of ammonia blended fuels. However, there remains uncertainties regarding their detailed reaction pathways and rate constants, hampering the development of high-fidelity chemical kinetic models. In this study, the kinetics of ĊH 3 + ṄH 2 , a key C–N interaction reaction in ammonia/methane blend combustion have been investigated. The potential energy surface has been explored using the high-level ANL0F method, yielding highly accurate stationary point energies that agree with ATcT values within 0.1 kcal mol –1 . Variable reaction coordinate transition state theory is used to treat the barrierless association and decomposition reaction channels, based on directly sampled radical-radical interaction energies at the CASPT2-F12(2e,2o)/cc-pVTZ-F12 level of theory. The minimum transitional mode numbers of states obtained are then coupled with the RRKM/master equation to calculate temperature- and pressure-dependent rate constants. Our a priori calculations capture available experimental measurements from the literature very well. The calculated rate constants have been incorporated into an NH 3 /CH 4 chemical kinetic model currently under development at the University of Galway. The effect of the updated kinetic data for ĊH 3 + ṄH 2 on model predicted NH 3 /CH 4 fuel reactivity is elucidated.

ab initio↗

Analysis of differential scanning calorimetry data for aged plutonium

Differential scanning calorimetry data for samples of a 52 year old plutonium alloy with 3.3 at. % Ga that were heated beyond the melting point is analyzed using transition state theory to find activation energies for the δ to ε and ε to liquid phase transitions. A Bayesian statistical method involving a Gaussian process model is used to find mean values and confidence intervals for the activation energies. The activation energy for the δ to ε phase transition increases by 3.3 ± 3.8% per decade, relative to the case when all age related plutonium lattice point defects have been removed through annealing. The corresponding increase in activation energy for the ε to liquid transition is shown to be 7.1 ± 1.8% per decade. It is postulated that the change in activation energy with age for both phase transitions is caused, in part, by the accumulation of the same type of lattice point defects associated with the observed increase in elastic bulk modulus over time.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Estimation of the reaction rate for the formation of CH3O from H + H2CO - Implications for chemistry in the solar system

Troe's (1977) approximate theory is presently used in conjunction with transition state theory to estimate the rate coefficient of the reaction by which CO is reduced to CH4; attention is given to the role that may be played in the reduction process by the formation of the CH3O radical from H + H2CO. Attention is given to the implications of such a reaction (1) for the CO chemistry on Jupiter and within the solar nebula, (2) for the interpretation of such experimental results as those of Bar-Nun and Shaviv (1975) and Bar-Nun and Chang (1983), and (3) for organic synthesis in the prebiotic terrestrial atmosphere.

Yung, Yuk L.↗

The reaction NH 2 + CH 2 O: Kinetic measurements and computational studies

The reaction of amidogen with formaldehyde is relevant to astrochemistry and couples nitrogen and carbon chemistry in flames. The first measurements of the temperature dependence of the rate constant k were made, over 409–643 K, and are summarized as k = 7.1 × 10 −12 exp(−18.2 kJ mol −1 /RT) cm 3 molecule −1 s −1 with a confidence interval of ±20%. These results address a controversy over the presence of a significant barrier in the addition/elimination pathway leading to the production of formamide between theoretical models that include the zero-point energy of all modes at the transition state and a model that excludes transitional modes. The latter would lead to a negligible barrier, which is inconsistent with the experiments. Contrary to earlier claims, CCSDT(Q)-corrected energies combined with transition state theory yield quantitative accord with the measurements.

Astrochemistry↗

A Robust Methodology to Elucidate Kinetics of Room Temperature Electrochemical Propane Adsorption on Platinum

Electrocatalytic activation of alkanes can further decarbonize chemical manufacturing by leveraging affordable renewable electricity and readily available shale gas reserves in the United States. Earlier works have identified the unique role of Pt in adsorbing and activating alkanes, like propane, at room temperature in acidic, aqueous electrolytes, revealing spontaneous formation of deeply dehydrogenated propane-derived surface species with an intact C 3 - backbone. Although an adsorption mechanism was hypothesized, it has not been explicitly investigated to date, preventing the quantification of kinetic rate parameters. A robust methodology to investigate and benchmark propane adsorption kinetics on Pt is critical for the rational design of electrocatalysts that exhibit higher selectivity toward desired partially oxidized products. Herein, we analyze an oxidative current transience that appears during the adsorption of propane on Pt in aqueous electrochemical conditions and develop a methodology that elucidates the adsorption mechanism and enables quantification of rate parameters such as order dependences and apparent activation barriers. This method yields an expected first-order dependence with respect to propane concentration at low coverage and reveals a second-order dependence with respect to the concentration of surface active sites. Additionally, the apparent activation barrier for propane adsorption was calculated using an Arrhenius analysis of the current transience under temperature control. The experimentally measured activation barrier of 35 kJ mol –1 is in excellent agreement with the theoretical barrier calculated by density functional theory (DFT). The kinetic analysis was extended, via the use of transition state theory, to extract entropy and enthalpy of activation, yielding consistent results with the proposed two-step adsorption mechanism and DFT calculations. These results demonstrate reliable quantification of kinetic parameters for electrocatalytic activation of C–H bonds in alkanes that can be employed for rational catalyst development for a versatile range of electrocatalytic conditions.

alkane activation↗

Calculation of ion–ion mutual neutralization rate constants using Landau–Zener theory coupled with trajectory simulations for Ar + –Cl − , Br − , I −

In this computational study, we self-consistently calculate the rate constants of mutual neutralization reactions by incorporating the electron transfer probability, using Landau–Zener state transition theory with inputs derived from ab initio quantum chemistry calculations, into classical trajectory simulations. Electronic structure calculations are done using correlation consistent basis sets with multi-reference configuration interaction to map all the molecular electronic states below the ion-dissociation limit as a function of the distance between the reacting species. Our electronic structure calculations have been significantly improved from our previous work through improved selection of molecular electronic configurations maintaining a fine grid of 1a 0 over a wide range of bond lengths and accurate treatment of spin–orbit couplings. Non-adiabatic coupling matrix elements are calculated with the three-point central difference method near each avoided crossing to estimate the exact crossing point R x and coupling parameter H if , which are inputs to the multi-channel Landau–Zener theory to calculate the electron transition probability. Our approach is applied to estimate the mutual neutralization rate constants for the following ion pairs: Ar + –Cl − , Ar + –Br − , Ar + –I − at ∼133 Pa. Furthermore, our predictions are compared against the experimental data reported. It is seen that the improvement in the electronic structure calculation results in excellent agreement between the simulation results and the available experimental data to within a factor of ∼2 or ∼±50%.

Complete-active space self-consistent field↗

ELECTRONIC STRUCTURE METHODS AND PROTOCOLS WITH APPLICATION TO DYNAMICS, KINETICS AND THERMOCHEMISTRY

Hydrocarbon combustion involves the reaction dynamics of a tremendous number of species beginning with many-component fuel mixtures and proceeding via a complex system of intermediates to form primary and secondary products. Combustion conditions corresponding to new advanced engines and/or alternative fuels rely increasingly on autoignition and low-temperature-combustion chemistry. In these regimes various transient radical species such as HO2, ROO·, ·QOOH, HCO, NO2, HOCO, and Criegee intermediates play important roles in determining the detailed as well as more general dynamics. A clear understanding and accurate representation of these processes is needed for effective modeling. Given the difficulties associated with making reliable experimental measurements of these systems, computation can play an important role in developing these energy technologies. Accurate calculations have their own challenges since even within the simplest dynamical approximations such as transition state theory, the rates depend exponentially on critical barrier heights and these may be sensitive to the level of quantum chemistry. Moreover, it is well-known that in many cases it is necessary to go beyond statistical theories and consider the dynamics. Quantum tunneling, resonances, radiative transitions, and non-adiabatic effects governed by spin-orbit or derivative coupling can be determining factors in those dynamics. Building upon progress made during a period of prior support through the DOE Early Career Program, this project combines developments in the areas of potential energy surface (PES) fitting and multistate multireference quantum chemistry to allow spectroscopically and dynamically/kinetically accurate investigations of key molecular systems (such as those mentioned above), many of which are radicals with strong multireference character and have the possibility of multiple electronic states contributing to the observed dynamics. An ongoing area of investigation is to develop general strategies for robustly convergent electronic structure theory for global multichannel reactive surfaces including diabatization of energy and other relevant surfaces such as dipole transition. Combining advances in ab initio methods with automated interpolative PES fitting allows the construction of high-quality PESs (incorporating thousands of high-level data) to be done rapidly through parallel processing on high-performance computing (HPC) clusters. In addition, new methods and approaches to electronic structure theory will be developed and tested through applications. This project will explore limitations in traditional multireference calculations (e.g., MRCI) such as those imposed by internal contraction, lack of high-order correlation treatment and poor scaling. Methods such as DMRG-based extended active-space CASSCF and various Quantum Monte Carlo (QMC) methods will be applied (including VMC/DMC and FCIQMC). Insight into the relative significance of different orbital spaces and the robustness of application of these approaches on leadership class computing architectures will be gained. Synergy with other components of this research program such as automated PES fitting and multireference quantum chemistry will be used to address challenges encountered by the standard approaches to computational thermochemistry (those being single-reference quantum chemistry and perturbative treatments of the anharmonic vibrational energy, which break down for some cases of electronic structure or floppy strongly coupled vibrational modes).

74 ATOMIC AND MOLECULAR PHYSICS↗