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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Redox-Induced Structural Reorganization Dictates Kinetics of Cobalt(III) Hydride Formation via Proton-Coupled Electron Transfer

Two-electron, one-proton reactions of a family of [CoCp(dxpe)(NCCH 3 )] 2+ complexes (Cp = cyclopentadienyl, dxpe = 1,2-bis(di(aryl/alkyl)phosphino)ethane) form the corresponding hydride species [HCoCp(dxpe)] + through a stepwise proton-coupled electron transfer process. For three [CoCp(dxpe)(NCCH 3 )] 2+ complexes (dxpe = dppe (1,2-bis(diphenylphosphino)ethane), depe (1,2-bis(diethylphosphino)ethane), and dcpe (1,2-bis(dicyclohexylphosphino)ethane)), peak shift analysis was employed to quantify apparent proton transfer rate constants from cyclic voltammograms recorded with acids ranging 22 pK a units. The apparent proton transfer rate constants correlate with the strength of the proton source for weak acids, but these apparent proton transfer rate constants curiously plateau (k pl ) as the reaction becomes increasingly exergonic. The absolute apparent proton transfer rate constants across both these regions correlate with the steric bulk of the chelating diphosphine ligand, with bulkier ligands leading to slower kinetics (k plateau,depe = 3.5 x 10 7 M -1 s -1 , k plateau,dppe = 1.7 x 10 7 M -1 s -1 , k plateau,dcpe = 7.1 x 10 4 M -1 s -1 ). Mechanistic studies were conducted to identify the cause of the aberrant $k_PT^{app}$–ΔpK a trends. When the deuterated acids are employed, deuterium incorporation in the Cp ring is observed, indicating protonation of the CoCp(dxpe) species to form the corresponding hydride proceeds via initial ligand protonation. Digital simulations of cyclic voltammograms show ligand loss accompanying initial reduction gates subsequent PCET activity at higher driving forces. Together, these experiments reveal the details of the reaction mechanism: reduction of the Co(III) species is followed by dissociation of the bound acetonitrile ligand, subsequent reduction of the unligated Co(II) species to form a Co(I) species is followed by protonation, which occurs at the Cp ring, followed by tautomerization to generate the stable Co(III)-hydride product [HCoCp(dxpe)] + . Analysis as a function of chelating disphosphine ligand, solvent, and acid strength reveals that the ligand dissociation equilibrium is directly influenced by the steric bulk of the phosphine ligands and gates protonation, giving rise to the plateau of apparent proton transfer rate constant with strong acids. The complexity of the reaction mechanism underpinning hydride formation, encompassing dynamic behavior of the entire ligand set, highlights the critical need to understand elementary reaction steps in proton-coupled electron transfer reactions.

14 SOLAR ENERGY↗
DOE OSTI2020

An Experimental and Theoretical Study of the High Temperature Reactions of the Four Butyl Radical Isomers

The high temperature gas phase chemistry of the four butyl radical isomers (n-butyl,sec-butyl, iso-butyl, andtert-butyl) was investigated in a combined experimental and theoretical study. Organic nitrites were used as convenient and clean sources of each of the butyl radical isomers. Rate coefficients for dissociation of each nitrite were obtained experimentally and are at, or close to, the high pressure limit. Low pressure experiments were performed in a diaphragmless shock tube with laser schlieren densitometry at post-shock pressures of 65, 130, and 260 Torr and post-shock temperatures of 700-1000 K. Additional experiments were conducted with iso-butyl radicals at 805 K and 8.7 bar to elucidate changes in mechanism at higher pressures. These experiments were performed in a miniature shock tube with synchrotron-based photoionization mass spectrometry. The mass spectra confirmed that scission of the O-NO bond is the primary channel by which the precursors dissociate, but they also provided evidence of a minor channel (<7.7%) through HNO loss and formation of an aldehyde. These high pressure experiments were also used to determine the disproportionation/recombination ratio for iso-butyl radicals as 0.3. Reanalysis of the lower-temperature literature and the present data yielded rate constants for the disproportionation reaction, iso-butyl + iso-butyl = iso-butene + iso-butane. A chemical kinetics model was developed for the reactions of the butyl isomers that included new paths for highly energized adducts. These adducts are formed by the addition of H, CH 3 or C 2 H 5 to the butyl radicals. Accompanying theoretical investigations show that chemically activated pathways are competitive with stabilization of the adduct by collision under the conditions of the laser schlieren experiments. These calculations also show that at 10 bar and T< 1000 K stabilization is the only important reaction, but at higher temperatures, even at 10 bar, chemically activated product channels should also be considered. Branching fractions and rate coefficients are presented for these reactions. Furthermore, this study also highlights the importance of the radical structure for determining branching ratios for disproportionation and recombination of alkyl radicals, and these were facilitated by theoretical calculations of recombination rate coefficients for the four butyl radical isomers. Overall, the results reveal previously unknown features of butyl radical chemistry under conditions that are relevant to a wide range of applications and reaction mechanisms are presented that incorporate pressure dependent rate coefficients for the key steps.

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

ATR-SEIRAS Reveals Potential Inversion and Associated Electron Transfer Kinetics in the Reduction of Surface-Confined Anthraquinone

The detection of stable semiquinone radicals on an anthraquinone (AQ) layer chemically grafted to an electrode surface in aqueous electrolytes has been elucidated by using attenuated total reflection surface enhanced infrared absorption spectroscopy (ATR-SEIRAS). In very alkaline conditions (pH 13), the reduction of the AQ involves no proton transfer, but surface sensitive infrared spectroscopy reveals that the anthraquinone dianion forms a strong hydrogen bonding network with coadsorbed water, leading to irreversible features in the voltammetry. The potential dependence of the IR band assigned to the AQ radical is consistent with the enhanced hydrogen bonding network causing increased stabilization of the quinone radical and supports the predicted response of a system under mild potential inversion, whereby the formal potential for the reduction of the anthraquinone radical is positive of the reduction potential of the neutral AQ molecule. Time-resolved ATR-SEIRAS is used to measure the transient formation of the AQ •– radical, from which rate constant information can be extracted using the Butler–Volmer model involving two one-electron transfers without a direct disproportionation reaction. The potential dependence of the rate constants is consistent with the potential inversion and can be used to qualitatively simulate the measured cyclic voltammograms. In conclusion, the thermodynamic and kinetic analyses re-emphasize long established deficiencies associated with using one-electron reaction formalisms to characterize multi-electron systems.

37 - INORGANIC, ORGANIC, PHYSICAL AND ANALYTICAL C↗
DOE OSTI2024

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

Measurements of Radical Reactivity with an Imine, (CF 3 ) 2 CNH: Rate Constants for Chlorine Atoms and Hydroxyl Radicals and the Global Warming Potential

The rate constant kOH for the reaction of 1,1,1,3,3,3-hexafluoroprop-2-imine with OH radicals was measured relative to two reference compounds, CH 3 F and CH 3 CHF 2 , to be k OH = (4.2 ± 1.1) × 10 −14 cm 3 molecule −1 s −1 at 295 K. This implies an atmospheric lifetime with respect to consumption by OH of 0.75 years. Reaction with Cl atoms yielded k Cl = (7.9 ± 1.7) × 10 −16 cm 3 molecule −1 s −1 at 295 K, and reaction with O 3 has an upper limit of k O3 < 4 × 10 −23 cm 3 molecule −1 s −1 , so that the atmospheric consumption by Cl and O 3 is negligibly slow. Absolute infrared cross sections of the imine yield a radiative efficiency of 0.34 W m −2 ppb −1 , which is corrected to 0.23 W m −2 ppb −1 for the effects of atmospheric lifetime. The imine’s corresponding 100-year global warming potential is 64 ± 19. This value is an upper limit, given that heterogenous atmospheric removal paths, such as hydrolysis in water droplets, are not included.

climate metric↗
DOE OSTI2020

Applicability of the Delplot method for the determination of catalytic reaction sequences: Hydrodeoxygenation of $\gamma$-valerolactone on Ni 2 P/MCM-41

Contact-time measurements based on sites are used to investigate the reaction sequence of the hydrodeoxygenation of γ-valerolactone (GVL) on Ni 2 P/MCM-41. Here the study is carried out at 300 °C and 0.5 MPa with site times of 0.29–6.1 s. A reaction network is proposed in which pentanoic acid and 2-methyltetrahydrofuran (2-MTHF) are formed in parallel initially and pentanal, butane and pentane are formed subsequently. A first-rank Delplot analysis (selectivity) gives a good fit, but higher-rank Delplots do not. The lack of applicability is explained by the magnitude of the rate constants, in which the rate constant for one of the later steps (the consumption of pentanal) is substantially larger than those for the preceding steps. For this reason, pentanal demonstrated the characteristics of a primary product and consequently affected the discrimination of higher rank products. A simulation of a model sequence with assumed rate constants confirms the reason why the Delplot method does not work for this reaction. Since the rate constants cannot be determined a priori to be of close magnitude, the Delplot method cannot be generally used to ascertain reaction sequences.

42 ENGINEERING↗
DOE OSTI2023

Impact of Pycnonuclear Fusion Uncertainties on the Cooling of Accreting Neutron Star Crusts

Abstract The observation of X-rays during quiescence from transiently accreting neutron stars provides unique clues about the nature of dense matter. This, however, requires extensive modeling of the crusts and matching the results to observations. The pycnonuclear fusion reaction rates implemented in these models are theoretically calculated by extending phenomenological expressions and have large uncertainties spanning many orders of magnitude. We present the first sensitivity studies of these pycnonuclear fusion reactions in realistic network calculations. We also couple the reaction network with the thermal evolution code dStar to further study their impact on the neutron star cooling curves in quiescence. Varying the pycnonuclear fusion reaction rates alters the depth at which nuclear heat is deposited although the total heating remains constant. The enhancement of the pycnonuclear fusion reaction rates leads to an overall shallower deposition of nuclear heat. The impurity factors are also altered depending on the type of ashes deposited on the crust. These total changes correspond to a variation of up to 9 eV in the modeled cooling curves. While this is not sufficient to explain the shallow heat source, it is comparable to the observational uncertainties and can still be important for modeling the neutron star crust.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Theoretical kinetics predictions for NH 2 + HO 2

Recent modeling studies of NH 3 oxidation, which are motivated by the prospective role of ammonia as a zero-carbon fuel, have indicated significant discrepancies between existing literature mechanisms. In this study high level theoretical kinetics predictions have been obtained for the reaction of NH 2 with HO 2 , which has previously been highlighted as an important reaction with high sensitivity and high uncertainty. The potential energy surface is explored with coupled cluster calculations including large basis sets and high-level corrections to yield high accuracy (~0.2 kcal/mol) estimates of the stationary point energies. Variational transition state theory is used to predict the microcanonical rate constants, which are then incorporated in master equation treatments of the temperature and pressure dependent kinetics. For the radical-radical channels, the microcanonical rates are obtained from variable reaction coordinate transition state theory implementing directly evaluated multireference electronic energies. The analysis yields predictions for the total rate constant as well as the branching to the NH 3 + O 2 , H 2 NO + OH, and HNO + H 2 O channels. Rate constants are also reported for the H 2 NO + OH reaction as they arise naturally from the analysis. The rate constant and branching fraction determined in this work for the NH 2 + HO 2 reaction deviate significantly from values used in most previous modeling studies. The fact that the main product channel is chain terminating, rather than propagating, has strong implications for modeling NH 3 ignition and oxidation, in particular at intermediate temperatures and elevated pressure.

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

Impact of cement composition, brine concentration, diffusion rate, reaction rate and boundary condition on self-sealing predictions for cement-CO 2 systems

Geological CO 2 storage (GCS) plays an important role in curbing CO 2 emissions by reducing the carbon footprint of difficult to decarbonize operations and achieve negative CO 2 emissions through activities like Bioenergy with Carbon Capture and Storage (BECCS) and Direct Air Carbon Capture and Storage (DACCS). Leakage of CO 2 through wells is an important concern when it comes to deployment of large-scale GCS. Existing wells in sites that are otherwise suitable for GCS can act as conduits for stored CO 2 to escape the reservoir. There is broad consensus that the main risk of leakage through wellbores is via fractures/damaged pathways. The results from several studies evaluating the permeability evolution of cement fractures in wells upon leakage of CO 2 agree that smaller fracture apertures, slower brine velocities and higher brine residence times promote self-sealing of fractures by mineral precipitation. Quantitatively, however, the differences in sealing conditions are significant and are typically attributed to differences in experimental conditions or model assumptions. Here we examine the sensitivity of our model, describing CO 2 leakage through wellbores, to cement composition, brine concentration, diffusion rates, and reaction rates. We also evaluate the impact of the boundary condition to allow comparisons between observations from experiments performed at constant flow rate and model predictions made at constant pressure conditions. Our results show that diffusion and reactions rates have the most impact on the self-sealing criteria for cement-CO 2 systems. In addition, conditions associated with self-sealing of fractures at constant flow rate require longer fractures, smaller fracture apertures, and slower velocities than under constant pressure.

58 GEOSCIENCES↗
DOE OSTI2022

An experimental and modeling study on auto-ignition kinetics of ammonia/methanol mixtures at intermediate temperature and high pressure

A rapid compression machine (RCM) has been applied to measure the ignition delay times of NH 3 /CH 3 OH mixtures covering pressures of 20 and 40 bar, equivalence ratios of 0.5, 1.0 and 2.0, and temperatures between 845 and 1100 K. Here the measurements show that the NH 3 /CH 3 OH mixtures become more reactive with increasing methanol addition. Addition of merely 1% (molar basis) of CH 3 OH to NH 3 lowers the ignition temperature around 100 K at 40 bar in comparison to pure NH 3 . The ignition delay is a complex function of fuel mixture and stoichiometry. For the 1% CH 3 OH mixture, the leaner mixtures are more reactive, while the reverse trend is found for mixtures with 5%, 20% and pure CH 3 OH. Analysis of the pressure profiles shows three distinct ignition modes for NH 3/ CH 3 OH mixtures, facilitated by the pre-ignition heat release from NH 3 consumption. A detailed mechanism for ignition of NH 3 /CH 3 OH fuel blends has been developed, capable of reproducing the ignition behavior of mixtures with reasonable accuracy. A subset for amine / methanol interactions was established, with rate constants for the key reaction between NH 2 and CH 3 OH calculated from ab initio theory. A sensitivity analysis indicates that the critical reactions during the auto-ignition process vary with the CH 3 OH mole fraction in the fuel. The ammonia chemistry, namely NH 2 + NO, NH2 + NO 2 and NH 2 + HO 2 , is dominant for the mixture with 1% CH 3 OH, while the reactions related to CH 3 OH and H 2 O 2 are more important for the 20% CH 3 OH mixture. The interaction between ammonia and methanol shows a more prominent effect on auto-ignition for mixtures with 5% CH 3 OH in fuel as compared to those with 1% and 20% CH 3 OH. According to the modeling results, methanol addition is found to enrich the O/H radical pool, consuming ammonia and promoting auto-ignition through different reaction pathways.

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

Protein Electric Fields Enable Faster and Longer-Lasting Covalent Inhibition of β-Lactamases

The widespread design of covalent drugs has focused on crafting reactive groups of proper electrophilicity and positioning towards targeted amino-acid nucleophiles. In this work, we found that environmental electric fields projected onto a reactive chemical bond, an overlooked design element, play essential roles in the covalent inhibition of TEM-1 beta-lactamase by avibactam. Using the vibrational Stark effect, the magnitudes of the electric fields that are exerted by TEM active sites onto avibactam’s reactive C=O were measured and demonstrate an electrostatic gating effect that promotes bond formation yet relatively suppresses the reverse dissociation. These results suggest new principles of covalent drug design and off-target site prediction. Unlike shape and electrostatic complementary which address binding constants, electrostatic catalysis drives reaction rates, essential for covalent inhibition, and deepens our understanding of chemical reactivity, selectivity, and stability in complex systems.

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

Phenalenyl growth reactions and implications for prenucleation chemistry of aromatics in flames

The energetics and kinetics of phenalene and phenalenyl growth reactions were studied theoretically. Rate constants of phenalene and phenalenyl H-abstraction and C 2 H 2 addition to the formed radicals were evaluated through quantum-chemical and rate-theory calculations. The obtained values, assigned to all π radicals, were tested in deterministic and kinetic Monte Carlo simulations of aromatics growth under conditions of laminar premixed flames. Kekulé and non-Kekulé structures of the polycyclic aromatic hydrocarbons (PAHs) evolving in the stochastic simulations were identified by on-the-fly constrained optimization. The numerical results demonstrated an increased PAH growth and qualitatively reproduced experimental observations of Homann and co-workers of non-decaying PAH concentrations with nearly equal abundances of even and odd carbon-atom PAHs. The analysis revealed that the PAH growth proceeds via alternating and sterically diverse acetylene and methyl HACA additions. Furthermore, the rapid and diverse spreading in the PAH population supports a nucleation model as PAH dimerization, assisted by the non-equilibrium phenomena, forming planar aromatics first and then transitioning to the PAH–PAH stacking with size.

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

Helium metastable density determination in the COST reference source by absolutely calibrated optical emission spectroscopy

Helium metastable densities in the COST Reference Microplasma Jet are estimated for a variety of He/N2 admixtures and dissipated powers by applying a collisional-radiative model to absolutely calibrated optical emission spectroscopy measurements. This is accomplished by delineating the excitation mechanisms that result in the N2(C–B) and N2+(B–X) emission bands, the latter of which is strongly coupled to the presence of helium metastables. A number of other plasma parameters are established and discussed for each operating condition including the electron energy distribution function, reduced electric field, rate constants, and electron density. With these parameters, the reaction rates for the primary ionization pathways are also calculated, emphasizing the importance of helium metastables for discharge sustainment. Good agreement with the existing literature is found for most plasma parameters and for helium metastable densities, in particular. A clear [N2]−1 relationship between the nitrogen concentration and density of helium metastables is demonstrated, as has been identified in previous studies in analogous atmospheric pressure plasma jets. This validates the efficacy of this optical technique for determining helium metastable densities and establishes it as a viable, and in many cases, more accessible alternative to other means of quantifying helium metastables in low-temperature plasmas.

Physics↗
DOE OSTI2020

Subsurface Nitrogen Dissociation Kinetics in Lithium Metal from Metadynamics

The dissociation of molecular nitrogen in lithium is of interest for several promising technologies, such as the catalytic synthesis of ammonia in ambient or mild conditions. In this work we simulate nitrogen dissociation in the lithium BCC (110) surface at ambient and elevated temperatures using density functional theory (DFT) metadynamics simulations. The rate constants at temperatures of 300, 400, and 500 K are calculated by statistical analysis of the reaction time distributions from the accelerated simulations. This approach finds and estimates rate constants for transition pathways out of the initial state; the required input is the stable initial state and a reasonable choice of collective variable. A single collective variable is used in this case: the N–N distance. The results are robust to changes in metadynamics parameters, and the reaction time distributions follow the expected exponential distribution. We show that the metadynamics-derived rate constants are in agreement with results from the conventional harmonic approximation approach using a climbing image nudged elastic band (NEB) transition state search. The reaction barriers from metadynamics and the NEB/harmonic approximation agree to within 0.02–0.04 eV at all temperatures studied. This paper demonstrates that the harmonic approximation provides an accurate description of the rate constants for nitrogen dissociation in lithium metal, even at temperatures near or above the melting point of lithium, lending credence to previous and future theoretical studies using this approximation. Moreover, this work demonstrates a step toward the automated exploration and discovery of reaction mechanisms and associated rate constants for elementary surface-catalyzed reactions using DFT-based metadynamics.

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

High-Entropy Ceramic Coatings: Transformative New Materials for Environmentally-Compatible Thin-Film Insulators against High-T Molten Salts

An emerging new class of promising materials, High-Entropy Ceramics (HECs) or Compositionally-Complex Ceramics (CCCs), have been successfully fabricated in the last three years. This project evaluates HECs/CCCs in the concentrated solar power (CSP) context. In this project, we have identified a couple of simple, and key descriptors governing thermal conductivity in HECs. Size disorder among the cations was found to have a strong influence on the thermal conductivity across multiple systems and throughout a vast compositional space. Since identifying this, we have optimized previous compositions to achieve low thermal conductivities of $k$~$1-1.4Wm{-1}K{-1}$. We have further reduced our selection of promising compositions based on project guidelines and applied a theory that can potentially predict reactivity with molten salts to identify a small set of HECs that are suitable for further molten salt heat transfer media (HTM) stability testing. Many compositions reveal poor resistance to molten chloride salt after being exposed for 24 h compared to the reference structural alloy, Haynes 230. However, after 100 h, a few HECs/CCCs show strong parabolic rate constant effects and provide greater resistance to chloride salt compared to Haynes 230 alloy, which shows more linear kinetics. Postmortem characterization of the microstructure reveals the parabolic rate constant originates from a dense reaction layer form at the interface. This feature is absent in Haynes 230, which only shows grain boundary infiltration. The HECs/CCCs were also put in contact with the metal alloy under Ar for 800 h at 800°C and show negligible interaction. Lastly, parabolic rate constant HECs/CCCs also have low thermal conductivity $(<2 Wm^{-1} K^{-1})$ at 800°C. These specimens satisfy all project requirements and demonstrate their promising application in next-generation CSP.

14 SOLAR ENERGY↗
DOE OSTI2023

Neutral gas pressure dependence of ion–ion mutual neutralization rate constants using Landau–Zener theory coupled with trajectory simulations

In this computational study, we describe a self-consistent trajectory simulation approach to capture the effect of neutral gas pressure on ion–ion mutual neutralization (MN) reactions. The electron transfer probability estimated using Landau–Zener (LZ) transition state theory is incorporated into classical trajectory simulations to elicit predictions of MN cross sections in vacuum and rate constants at finite neutral gas pressures. Electronic structure calculations with multireference configuration interaction and large correlation consistent basis sets are used to derive inputs to the LZ theory. The key advance of our trajectory simulation approach is the inclusion of the effect of ion-neutral interactions on MN using a Langevin representation of the effect of background gas on ion transport. For H + – H – and Li + – H(D) – , our approach quantitatively agrees with measured speed-dependent cross sections for up to ~10 5 m/s. For the ion pair Ne + – Cl – , our predictions of the MN rate constant at ~1 Torr are a factor of ~2 to 3 higher than the experimentally measured value. Similarly, for Xe + – F – in the pressure range of ~20 000–80 000 Pa, our predictions of the MN rate constant are ~20% lower but are in excellent qualitative agreement with experimental data. Here, the paradigm of using trajectory simulations to self-consistently capture the effect of gas pressure on MN reactions advanced here provides avenues for the inclusion of additional nonclassical effects in future work.

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

Isothermal Microcalorimetric Analysis of Li/CF x Battery Discharge

Lithium/carbon monofluoride (Li/CF x ) batteries generate substantial amounts of heat during discharge in part due to the large deviation of the loaded voltage from the measured and thermodynamically predicted open circuit voltage. Here, in this study, we further analyze this system by estimating the equilibrium voltage (V eq ) and its temperature dependence (dV eq /dT) over the entire discharge range. Based on these results, the ohmic and entropic heat contributions to the overall heat flow are calculated from experimental data. The ohmic heat flow is consistent with an electrochemical impedance spectroscopy model and electrode density measurements, indicating decreasing porosity of the electrode during discharge. Entropic heating increases during discharge as the cell reaction becomes increasingly entropically unfavorable. The total energy dissipated by the cell (electrical and thermal) remains similar over the entire discharge, evidence for a two-phase reaction with a constant, and rate-independent, ΔG rxn of about -465 kJ mol -1 . However, near the end of discharge (>90% DOD), ΔG rxn changes significantly indicating possible secondary reactions.

25 ENERGY STORAGE↗
DOE OSTI2023

Automated identification and calculation of prompt effects in kinetic mechanisms using statistical models

The kinetics of prompt dissociation involves rovibrationally excited species (generally formed by exothermic reactions) which may dissociate or isomerize prior to thermalization via collisions with the bath gas. Treating such rovibrationally excited species (so-called "hot" species) with standard kinetic phenomenology may result in incorrect macroscopic representation of their reactivity. Here this work presents the first fully automated methodology for the calculation of prompt effects of a chosen species in a kinetic mechanism, including (i) reaction selection; (ii) theoretical calculation of rate constants and prompt branching fractions; and (iii) final rate constant fitting. The energy partition between hot fragments is estimated using a variety of statistical models, including a new physically sound microcanonical statistical model based on the rovibrational density of states of the fragments. The methodology is validated against literature data for the prompt dissociations of HCO and C 3 H 7 radicals. The microcanonical statistical model is in better agreement with trajectory simulations for larger species and is thus applicable for practical systems that typically involve large molecules, for which direct dynamics calculations are impractical. The automated workflow is applied to the evaluation of the effects of prompt dissociation for two isomeric radicals C 4 H 7 1-3 (1-methylallyl) and C 4 H 7 1-4 (3-buten-1-yl). Twelve H-atom abstraction reactions are selected and the corresponding rate constants are computed with first principles theory. The microcanonical statistical model predicts that prompt dissociations of C 4 H 7 1-3 and C 4 H 7 1-4 are already significant at 1000K, resulting in differences of up to an order of magnitude at 2000K with respect to the phenomenological thermal rate constants. To illustrate the effects of prompt dissociation on simulations of experimental data, the calculated prompt rate constants are implemented in both CRECK and C3MechV3.3 kinetic mechanisms. Simulations of experimental flame data illustrate the noticeable impact of prompt dissociation kinetics on the high-temperature combustion reactivity of C 4 H 8 -1 and C 4 H 8 -2.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗