Search NASA⌕ Search

SEARCH · Search NASA

Results for “Reaction rate constants”

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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

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↗
DOE OSTI2025

Coupling to rotational manifolds to improve gas-phase pump–probe spectroscopic models

The physical picture of gas-phase optical transitions is normally presented as an isolated two-level system balanced by upward and downward processes. Isolated models assume a phenomenological treatment of collisional dephasing but do not strictly account for collisional population exchange with the rotational baths. While this assumption is valid under low-intensity conditions, where excitation is rate-limiting, isolated models can deviate from Beer’s Law at sufficient pressures and monochromatic intensities when both collisional broadening and power broadening are comparable to (or greater than) lifetime broadening, which are not uncommon conditions for cavity enhanced spectroscopies in the mid-IR spectral range. Although this problem has been addressed by rate-equation models for linear absorption measurements, a general treatment for multi-level quantum mechanical models suitable for non-linear absorption measurements (two-photon/two-color/pump–probe) is lacking. Isolated models require physical parameter inputs that disagree with expected values by at least an order of magnitude. These non-physical models undermine the ability to predict non-linear signal strengths under untested conditions and thereby limit the potential to optimize the sensitivity of non-linear spectroscopies and to expand their analytical applications (e.g., new analytes and/or buffer gases, changes in cavity free-spectral-range, changes in intracavity powers or wavelengths, and accurate investigation of physical phenomena). In this study, we derive bath-coupled models for gaseous pump–probe spectroscopy by application of the quantum Lindblad equation and detailed balance. Bath-coupled models are shown to fit data consistently across variations in intensity and agree with all physically expected values.

Cavity ring-down spectroscopy↗
DOE OSTI2026

Dynamics of iodine geminate recombination in supercritical xenon solvent: Caging effect

Understanding the dynamics of chemical reactions in solutions is vital, as their rates and kinetics are significantly affected by the solvent environment. Supercritical solvents offer extensive applications in chemical reactions by enabling the manipulation of the solution environment. Here, in this study, we investigate the geminate recombination of iodine in a supercritical xenon solvent by using ReaxFF-based molecular dynamics simulations. Our findings reveal that the highest iodine recombination rate occurs near supercritical conditions, while lower-pressure conditions lead to reduced collision rates and unstable recombination, and higher-pressure conditions hinder iodine diffusion, resulting in a lower recombination rate. Our analysis shows that the xenon local density at the time of recombination is at least 2.5 times higher than the global density, confirming the presence of xenon clusters surrounding the Iodine atoms. This observation is further supported by coordination number analysis, which confirms an elevated xenon local density during recombination. In addition, the correlation between the total energy of xenon atoms within a cluster and recombined iodine atoms underscores the kinetic energy transfer process, validating the occurrence of geminate recombination. The excess kinetic energy from the recombining iodine atoms is transferred to the surrounding xenon atoms. Our examination of geminate recombination demonstrates that iodine atoms confined within xenon clusters—whether through manual insertion of atoms or the fast dissociation of an iodine molecule within xenon clusters—are more likely to recombine as primary geminate recombination. However, extending the iodine molecule dissociation time allows iodine atoms to diffuse out of the cluster, and the recombination to shift toward secondary geminate recombination.

Cage effect↗
DOE OSTI2026

General framework for quantifying dissipation pathways in open quantum systems. III. Off-diagonal subsystem–bath couplings

This paper extends the previously reported theory of dissipation pathways [C. W. Kim and I. Franco, J. Chem. Phys. 160, 214111 (2024)] to incorporate off-diagonal subsystem–bath coupling, which is often required to model molecular systems where the environment directly influences transitions and couplings between subsystem states. We systematically derive master equations for both population transfer and dissipation into individual bath components, for which we also rigorously prove energy conservation and detailed balance. The approach is based on second-order perturbation theory with respect to the subsystem–bath couplings, whose form is not limited to any specific model. The accuracy of the developed method is tested by applying it to diverse model Hamiltonians involving linearly coupled harmonic oscillator baths and comparing the outcomes against the hierarchical equations of motion (HEOM) method. Overall, our method accurately quantifies the contributions of specific bath components to the overall dissipation while significantly reducing the computational cost compared to numerically exact methods such as HEOM, thus offering a path to examine how vibronic interactions steer non-adiabatic processes in realistic chemical systems.

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

Analytic rate theory of polariton relaxation that explains long polariton lifetime

Hybridization of a molecular exciton with a quantized photon creates a polariton. Despite extensive experimental investigations, the apparent lifetime of the exciton–polariton is not well-understood. For this work, we examined the steady-state population dynamics for a Holstein–Tavis–Cumming Hamiltonian to illuminate the long-term polaritonic dynamics and lifetime of the exciton–polariton in an optical cavity. For a realistic description of polariton relaxation, cavity loss and various exciton decay channels are included in the model. We found that in the presence of weak but finite exciton loss, the apparent lifetime of the lower polariton coincides with the out-of-cavity exciton lifetime and is independent of cavity-matter detuning. This is a simple explanation for the experimentally observed lifetimes for exciton polaritons and theoretically justifies the dark state reservoir hypothesis. Furthermore, if the upper polariton is initially populated, the system reaches the steady state very quickly, leading to single-exponential polariton relaxation. Starting from the lower polariton leads to a longer pre-steady-state time period, leading to double-exponential relaxation. Finally, we considered the effect of site orientational disorders and the exciton frequency disorderers. Under the collective limit, the effects of this disorder can be included in Fermi’s golden rule population dynamics without explicit sampling. For the exciton energy disorders, numerical calculations are needed. Our theoretical framework is applicable to interpret exciton–polariton experiments, especially related to the measured apparent lifetime of polaritons.

Chemical dynamics↗
DOE OSTI2024

Absolute Reactive Oxygen Species Densities in the Effluent of the COST Reference Source and Plasma-generated Atomic Oxygen Density Measurements in Liquid using TALIF

Reactive oxygen species are one of the most important species produced in a plasma. In particular the role of atomic oxygen and the hydroxyl radical OH has attracted many studies ranging from plasma medicine to plasma-driven biocatalysis and water treatment. Atomic oxygen and the hydroxyl radical OH are very reactive and react unselectively with many different substrates. The quantification of atomic oxygen and OH is challenging because of their reactivity and unknown reactions and rate constants.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Applications of Kinetic Methods in Thermal Analysis: A Review

Determining kinetic parameters such as activation energy (Ea), pre-exponential factor (A), rate constants, and reaction order (n) via thermal analysis techniques is important to material synthesis and fabrication, industrial production, biomedicals, energy storage, catalysis, etc. Various kinetic methods based on the non-isothermal techniques have been developed to obtain the kinetic parameters. In this review, we focus on introducing various kinetic methods in non-isothermal analysis and their applications in different reactions, including thermal decomposition, solid-state phase transformation, crystallization, thermal ignition, curing process, etc.

Zhang, Xin↗
DOE OSTI2021

Dynamics of Product Branching from Radical Intermediates in Elementary Combustion Reactions (Final Technical Report)

While the total rate constant for many elementary reactions is well-characterized, understanding the product branching in complex reactions presents a formidable challenge. To gain an incisive probe of such reactions, our experiments investigate the dynamics of the product channels that arise from transient radical intermediates along the bimolecular reaction coordinates. Our recent published work uses the methodology developed in my group in the last fifteen years, using both imaging and scattering apparatuses. The experiments generate a particular isomeric form of an unstable radical intermediate along a bimolecular reaction coordinate and study the branching between the ensuing product channels of the energized radical as a function of its internal rotational and vibrational energy under collision-less conditions. The experiments use a combination of: 1) measurement of product velocity and angular distributions in a crossed laser-molecular beam apparatus, with electron bombardment detection in my lab in Chicago or 2) with tunable vacuum ultraviolet photoionization detection at Taiwan's National Synchrotron Radiation Research Center (NSRRC), and 3) velocity map imaging using state-selective REMPI and single photon VUV ionization of radical intermediates and reaction products. We also employ tunable VUV photoionization detection in our imaging apparatus, using difference frequency four-wave mixing to produce photoionization light tunable from 8 to 10.8 eV. This year of support funded our work to photolytically generate the radical intermediate important in the OH + propene reaction when OH adds to an end carbon. Our other major result this past year was on the photoionization detection of vinoxy radicals. Our experiments at the NSRRC definitively show that vinoxy dissociatively ionizes to both CH 3 + and HCO + , and establishes the absolute partial photoionization cross section of vinoxy to CH 3 + at 10.5 and 11.44 eV. This should allow experimentalists detecting vinoxy radicals produced in bimolecular reactions such as OH + propene to get good branching fractions to that product channel. Finally we initiated studies of the photodissociation of 1-bromo-2-chloroethane and 1,1-bromochloroethane at 193 nm at the NSRRC.

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

Reactions of the ˙NO 3 radical with nuclear extraction ligands in alkane solution

Here, the reactivity of the nitrate radical (NO 3 ˙) with organophosphorus and amidic actinide and lanthanide complexing agents of interest to nuclear solvent extraction applications was measured, resulting in the first-ever reported bimolecular rate constants for this radicals’ reactions in dodecane solution. The order of reactivity for neutral organophosphorus compounds showed faster rate constants with increasing electron density on the phosphoryl phosphorus atom, indicating an increasing facility for electron abstraction reactions occurring in addition to H-atom abstraction from the ligand alkane chains. The only acidic organophosphorus compound investigated, HEH[EHP], showed low reactivity with the NO 3 ˙ radical, attributed to its dimerization in this non-polar solvent. Amide ligand reaction rates were faster than for organophosphorus molecules, suggesting more facile H-atom abstraction from carbonyl activated methylene and amyl groups. While all rate constants were slower than the diffusion-limited rate they were still rapid enough to result in significant oxidation of solvent extraction ligands in dodecane solution.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2020

Kinetic study of reaction C 2 H 5 + HO 2 in a photolysis reactor with time-resolved Faraday rotation spectroscopy

The rate constant and branching ratios of ethyl reaction with hydroperoxyl radical, C 2 H 5 + HO 2 (1), a key radical-radical reaction for intermediate temperature combustion chemistry, were measured in situ for the first time in a photolysis Herriott cell by using mid-IR Faraday rotation spectroscopy (FRS) and UV-IR direct absorption spectroscopy (DAS). The microsecond time-resolved diagnostic technique in this work enabled the direct rate measurements of the target reaction at 40 and 80 mbar and reduced the experimental uncertainty considerably. C 2 H 5 and HO 2 radicals were generated by the photolysis of (COCl) 2 /C 2 H 5 I/CH 3 OH/O 2 /He mixture at 266 nm. By direct measurements of the transient profiles of C 2 H 5 , HO 2 and OH concentrations, the overall rate constant for this reaction at 297 K was determined as k 1 (40 mbar) = (3.8 ± 0.8) × 10 –11 cm 3 molecule –1 s –1 and k 1 (80 mbar) = (4.1 ± 1.0) × 10 –11 cm 3 molecule –1 s –1 . As a result, the direct observation of hydroxyl radical (OH) indicated that OH formation channel was the major channel with a branching ratio of 0.8 ± 0.1.

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

The kinetics of the reactions of ground-state NH with N 2 O and implications for ammonia combustion

Quantum chemistry and canonical transition state theory have been applied to derive rate constants k for reactions of 3 Σ − imidogen with nitrous oxide. 3 channels were quantified. The fastest leads to NNH+NO and may be summarized as k = 2.8 × 10 7 T 1.68 exp.(−23.1 kcal mol −1 /RT) cm 3 mole −1 s −1 . Two much slower pathways lead to N 2 + 3 HNO and spin-forbidden 1 HNO. Modeling of literature data, from an NH 3 /O 2 flame and jet-stirred reactor measurements on dilute NH 3 /N 2 O mixtures, with the new kinetic information suggests that these channels are too slow to consume N 2 O significantly, and confirms that a widely used proposed rate constant for NH+N 2 O overestimates the reactivity.

Marshall, Paul [University of North Texas, Denton,↗
DOE OSTI2025

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↗
DOE OSTI2025

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

Experimental and kinetic modeling study of tetralin: A naphtheno-aromatic fuel for gasoline, jet and diesel surrogates

Distillate fuels contain significant proportions of naphtheno-aromatic components and tetralin is a suitable surrogate component to represent this molecular moiety. The presence of aromatic and naphthyl rings makes kinetic modeling of tetralin very challenging. Primary radicals formed during the oxidation of tetralin can be aryl, benzylic or paraffinic in nature. Using available information on reaction paths and rate constants of naphthenes and alkyl-aromatics, in this study a kinetic model of tetralin has been developed with emphasis on low-temperature chemistry and high-pressure conditions. Due to the lack of high-level quantum chemical calculations on reaction pathways of tetralin, analogous rates from ab-initio studies on benzylic and paraffinic radicals have been adopted here. Some modifications to the reaction rate rules are incorporated to account for the unique characteristics of tetralin's molecular structure. Important reaction channels have been identified using reaction path and brute force sensitivity analyses. In order to investigate the model performance at low temperatures, new experiments are carried out in a rapid compression machine on blends of tetralin and 3-methylpentane. Blending of low-reactivity tetralin with a high-reactivity alkane allowed the investigation of tetralin ignition at very low temperatures (665 – 856 K). The kinetic model developed in the current study is found to predict the current experiments and literature data adequately. The new model will aid in high-fidelity surrogate predictions at engine-relevant conditions.

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

Challenges in Kinetic modeling of ammonia pyrolysis

Ammonia pyrolysis reactions have implications for its ignition and oxidation in engines and gas turbines. In the present work, the chemistry of ammonia pyrolysis is investigated by kinetic modeling and theory. Rate constants for key reactions are carefully evaluated based on available experimental and theoretical results. The high pressure limit k1,∞ for NH 2 + H (+M) ⇄ NH 3 (+M) (R1) is calculated to be essentially the collision frequency, indicating that dissociation of ammonia in combustion processes will be at or close to the low pressure limit even at engine and gas turbine conditions. The chemical kinetic model is validated against reported shock tube measurements of NH 3 , NH 2 , and NH in ammonia pyrolysis. Predictions are in good agreement with observations for dilute conditions (≤ 0.5% NH 3 ), but the model appears to underpredict the NH 3 consumption rate at longer times in less dilute mixtures. At short reaction times, thermal dissociation of NH3, together with the NH 3 + H reaction, controls conversion. At longer times, secondary reactions involving NH 2 and NH become important due to their impact on the radical pool. Predictions become sensitive to formation and consumption of diazene (tHNNH and cHNNH). Several of the key steps in the ammonia pyrolysis mechanism are radical-radical reactions that are difficult to measure accurately and challenging to calculate theoretically, and a more comprehensive experimental characterization is desirable to support further model development.

10 SYNTHETIC FUELS↗
DOE OSTI2020

An extended methodology for automated calculations of non-Boltzmann kinetic sequences: H + C 2 H 2 + X and combustion impact

It is generally assumed in phenomenological kinetic models that bimolecular reactions only occur between species whose rovibrational energy follows a Boltzmann (thermal) distribution. That is, any complexes initially formed in non-Boltzmann distributions are assumed to be thermalized by energy-transferring collisions prior to bimolecular reactions. Given the high mole fractions of reactive species, X, in combustion environments, reactive collisions of the complexes with X often occur on the same timescale as energy-transferring collisions – yielding sequences proceeding through non-Boltzmann intermediates across multiple potential energy surfaces. Recent studies have shown that such non-Boltzmann kinetic sequences can have substantial impact on the global reactivity in combustion systems. Simulations of these non-Boltzmann reaction sequences, which can be described in phenomenological kinetic models via chemically termolecular reactions, require that rovibrational excitation from one potential energy surface be carried over to the next. Here this paper presents an extended theoretical and computational methodology that couples multiple master equations and derives rate constants for phenomenological reactions describing the conversion of thermal reactants to thermal products for use in phenomenological kinetic schemes. The methodology is then implemented using in-house scripts for non-Boltzmann sequences involving C2H$^*_3$ + X (with X = O 2 , H, and OH) where C 2 H$^*_3$ is formed via H + C 2 H 2 association – which were identified as having strong potential for influencing combustion predictions in a recent study. The results reveal that non-Boltzmann reaction sequences for X = O 2 (the primary focus of this paper) significantly alters the total conversion rate from H + C 2 H 2 to products and product branching fractions from those of thermal sequential pathways. Furthermore, the present results demonstrate that non-Boltzmann reaction sequences have significant impact – as high as an order of magnitude – on predicted ignition delay times. Similarly, they yield significantly different dependence of ignition delay times with temperature and O 2 mole fraction – yielding signatures that are likely observable experimentally.

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

Kinetics and Energetics of Electron Transfer to Dimer Radical Cations

In this work, spectra of the dimer cations naphthalene (Nap 2 •+ ) and ethene (Ethene 2 •+ ) were measured in liquid dichloromethane (DCM). The spectra peak at very different energies, 1.2 and 3.3 eV. In DCM dimerization stabilizes Nap 2 •+ by Δ$G_\text{d}°$(Nap 2 •+ ) = -218 meV relative to the monomer Nap •+ as determined from the dimerization equilibrium constant. Both dimers can transfer a positive charge to hole acceptor molecules, but for both the rate constants rise more gradually with reaction energetics than do many charge transfer reactions previously studied. A striking observation finds that the rate constant for hole transfer from the Nap 2 •+ dimer to phenanthrene is smaller by two decades than that from biphenyl•+ monomer to Nap, although both reactions have the same –Δ$\textit{G}°$=0.05 eV. A plausible interpretation for these observations is the presence of an energy of reorganization, λ(M 2 ), for the dimer that involves movement apart of the two partners in the dimer. While the dimerization equilibrium cannot be measured for Ethene 2 •+ the charge transfer data imply that both Δ$G_\text{d}°$(Ethene 2 •+ ) and λ(Ethene 2 •+ ), are considerably larger, perhaps by factors of 2-4 than for Nap 2 •+ .

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

Enhanced Menshutkin S N 2 Reactivity in Mesoporous Silica: The Influence of Surface Catalysis and Confinement

A significant enhancement in the Menshutkin S N 2 reaction between 1-methylimidazole (MeIm) and methyl thiocyanate (MeSCN) is observed when the reaction is confined in the nanoscale silica pores of MCM41 and SBA15. The experiments in the silica pores are conducted without surrounding bulk reaction mixture. The influences of temperature, pore radius, and surface chemistry on the kinetics of the confined reaction are analyzed with time-dependent infrared spectroscopy, molecular dynamics simulations, and ab initio calculations. The rate constant of the pseudo-first order reaction increases with decreasing pore size, and the activation energy is found to decrease by 5.6 kJ/mol in the smallest pore studied (2.8 nm) relative to the bulk reaction. The rate constant dependence on pore size is accurately described by a two-state model in which molecules within the 4.6 Å interfacial layer experience a 2.4-fold rate constant increase relative to those reacting at the bulk rate further away from the interface. The removal of polar silanol groups from the silica surface via passivation with trimethylsilyl chloride results in bulk-like kinetics despite a reduction in the pore diameter, demonstrating the role of silanols as catalytic sites. Here, electronic structure calculations of the energy profile on a model silica surface confirm that silanol groups, particularly those of the vicinal type, can reduce the activation energy and reaction endothermicity through the donation of hydrogen bonds to the reactant, transition state, and product complexes.

36 MATERIALS SCIENCE↗