Search NASA⌕ Search

SEARCH · Search NASA

Results for “Arrhenius”

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 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Electron transfer in nonpolar media

The Marcus theory of electron transfer views fluctuating orientations of permanent dipoles as the nuclear mode bringing the donor and acceptor into the tunneling resonance. Electronic polarization of the solvent is excluded as the fast mode adiabatically following the electronic density. This view, valid for solids, does not apply to molecular liquids where molecular translations (density fluctuations) modulate the induction interaction of the donor-acceptor complex with the solvent. This mechanism of promoting radiation-less electronic transitions is considered here in the framework of the perturbation liquid-state theory. The reorganization energy of electron transfer in nonpolar solvents is nonzero and reaches the values of 0.1–0.3 eV for typical molecular sizes and solvents used in applications. The reorganization energy scales quadratically with the molecular polarizability of the solvent and decays as the inverse fifth power with the size of the donor and acceptor. The combination of the entropic character of the density fluctuations, driven by re-packing of molecular cores, with the short range of induction solute–solvent interactions leads to the violation of the fluctuation–dissipation theorem for the variance of the donor-acceptor energy gap. An explicit, approximately hyperbolic, dependence of the reorganization energy on temperature is predicted. It leads to a non-Arrhenius kinetic law for the rate constant of electron transfer in nonpolar liquid solvents.

14 SOLAR ENERGY↗
DOE OSTI2021

Mechanism and kinetics of the oxidation of 1,3-butadien-1-yl (n-C 4 H 5 ): a theoretical study

Here, ab initio CCSD(T)-F12/cc-pVTZ-f12//B3LYP/6-311G(d,p) calculations of the C 4 H 5 O 2 potential energy surface have been combined with Rice–Ramsperger–Kassel–Marcus Master Equation (RRKM-ME) calculations of temperature- and pressure-dependent rate constants and product branching ratios to unravel the mechanism and kinetics of the n-C 4 H 5 + O 2 reaction. The results indicate that the reaction is fast, with the total rate constant being in the range of 3.4–5.6 × 10 –11 cm 3 molecule –1 s –1 . The main products include 1-oxo-n-butadienyl + O and acrolein + HCO, with their cumulative yield exceeding 90% at temperatures above 1500 K. Two conformers of 1-oxo-n-butadienyl + O are formed via a simple mechanism of O 2 addition to the radical site of n-C 4 H 5 followed by the cleavage of the O–O bond proceeding via a van der Waals C 4 H 5 O∙∙∙O complex. Alternatively, the pathways leading to acrolein + HCO involve significant reorganization of the heavy-atom skeleton either via formal migration of one O atom to the opposite end of the molecule or its insertion into the C1–C2 bond. Not counting thermal stabilization of the initial peroxy adducts, which prevails at low temperatures and high pressures, all other products share a minor yield of under 5%. Rate constants for the significant reaction channels have been fitted to modified Arrhenius expressions and are proposed for kinetic modeling of the oxidation of aromatic molecules and 1,3-butadiene. As a secondary reaction, n-C 4 H 5 + O 2 can be a source for the formation of acrolein observed experimentally in oxidation of the phenyl radical at low combustion temperatures, whereas another significant (secondary) product of the C 6 H 5 + O 2 reaction, furan, could be formed through unimolecular decomposition of 1-oxo-n-butadienyl. Both the n-C 4 H 5 + O 2 reaction and unimolecular decomposition of its 1-oxo-n-butadienyl primary product are shown not to be a substantial source of ketene.

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

Reactivity of Zn + aq in high-temperature water radiolysis

For this work, reactivity of transients involving Zn + in high-temperature water radiolysis has been studied in the temperature range of 25–300 °C. The reduced monovalent zinc species were generated from an electron transfer process between the hydrated electron and Zn 2+ ions using pulse radiolysis. The Zn + species can subsequently be oxidized by the radiolytically-produced oxidizing species: ˙OH, H 2 O 2 and ˙H. We find that the absorption of monovalent zinc is very sensitive to the pH of the medium. An absorption maximum at 306–311 nm is most pronounced at pH 7 and the signal then decreases in acidic media where the reducing electrons are competitively captured by protons. At pH values higher than 7, hydroxo-forms of Zn 2+ are created and the maximum of the absorption signal begins to shift to the red spectral region. We find that the optical spectrum of Zn + aq cannot be fully explained in terms of a charge-transfer to solvent (CTTS) process, which was previously proposed. Reaction rates of most of the recombination reactions investigated follow the empirical Arrhenius relationship at temperatures up to 200 °C and have been determined at higher temperatures for the first time. A bimolecular disproportionation reaction of Zn + aq is not observed under the conditions investigated.

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

Development of a continuous synthesis process for carbamazepine using validated in-line Raman spectroscopy and kinetic modelling for disturbance simulation

Mitigation of failure modes in the continuous synthesis (CS) of a drug substance (DS) has the potential to widen the adoption of continuous manufacturing (CM) technologies by the pharmaceutical industry. Here, this work demonstrates the development of a robust continuous process for the synthesis of carbamazepine (CBZ), an essential medicine as per the World Health Organization (WHO), facilitated by kinetic modelling and monitored by in-line Raman spectroscopy. Accurate kinetic modelling and the use of validated process analytical technology (PAT) models for quantitative measurement were found to play an important role in developing CS of drug substances. Kinetic data for the formation of CBZ from iminostilbene (ISB) were collected by batch reaction sampling and high-performance liquid chromatography (HPLC) analysis. A non-linear solver and iterative method was applied to determine two sets of Arrhenius parameters simultaneously for the reaction system by minimizing the standard error of the model fit. The start-up and dynamic equilibrium stages for the CS of CBZ using a continuous stirred tank reactor (CSTR) were modelled based on the batch kinetic data and employed to optimize conversion and simulate process disturbances. An in-line Raman spectroscopy method was successfully developed, validated, and integrated to determine the concentrations of CBZ and ISB within the operating range for the CS. The CS kinetic model was evaluated experimentally from startup to dynamic equilibrium over 10 residence times with monitoring by HPLC and in-line Raman spectroscopy. The developed kinetic model in tandem with in-line Raman spectroscopy successfully predicted disturbances due to changes in process variables and can serve as a useful tool in the future design of advanced process control strategies for the continuous synthesis of CBZ.

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

Impact of lanthanide ion complexation and temperature on the chemical reactivity of N , N , N ', N '-tetraoctyl diglycolamide (TODGA) with the dodecane radical cation

The impact of trivalent lanthanide ion complexation and temperature on the chemical reactivity of N,N,N',N'-tetraoctyl diglycolamide (TODGA) with the n-dodecane radical cation (RH˙+) has been measured by electron pulse radiolysis and evaluated by quantum mechanical calculations. Additionally, Arrhenius parameters were determined for the reaction of the non-complexed TODGA ligand with the RH˙ + from 10–40 °C, giving the activation energy (E a = 17.43 ± 1.64 kJ mol –1 ) and pre-exponential factor (A = (2.36 ± 0.05) × 10 13 M –1 s –1 ). The complexation of Nd(III), Gd(III), and Yb(III) ions by TODGA yielded [LnIII(TODGA)3(NO3)3] complexes that exhibited significantly increased reactivity (up to 9.3× faster) with the RH˙ + , relative to the non-complexed ligand: k([Ln III (TODGA) 3 (NO 3 ) 3 ] + RH˙ + ) = (8.99 ± 0.93) × 10 10 , (2.88 ± 0.40) × 10 10 , and (1.53 ± 0.34) × 10 10 M –1 s –1 , for Nd(III), Gd(III), and Yb(III) ions, respectively. The rate coefficient enhancement measured for these complexes exhibited a dependence on atomic number, decreasing as the lanthanide series was traversed. Preliminary reaction free energy calculations—based on a model [Ln III (TOGDA)] 3+ complex system—indicate that both electron/hole and proton transfer reactions are energetically unfavorable for complexed TODGA. Furthermore, complementary average local ionization energy calculations showed that the most reactive region of model N,N,N',N'-tetraethyl diglycolamide (TEDGA) complexes, [Ln III (TEGDA) 3 (NO 3 ) 3 ], toward electrophilic attack is for the coordinated nitrate (NO 3 – ) counter anions. Furthermore, it is possible that radical reactions with the complexed NO 3 – counter anions dominate the differences in rates seen for the [Ln III (TODGA) 3 (NO 3 ) 3 ] complexes, and are likely responsible for the reported radioprotection in the presence of TODGA complexes.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2025

Decoupling electrode kinetics to elucidate reaction mechanisms in alkaline water electrolysis

Alkaline water electrolysis (AWE) presents key advantages, including reduced material costs, enhanced operational stability, and compatibility with non-precious metal catalysts, positioning it as a scalable route for hydrogen production. In this study, we introduce a minimally invasive single-cell configuration incorporating a reference electrode via diaphragm extension to form an internal ion channel. This setup, combined with an interfaced potentiostat and auxiliary electrometer, enables real-time, independent monitoring of anode and cathode behavior, offering high-resolution electrochemical diagnostics. While it is well established that the hydrogen evolution reaction (HER) exhibits sluggish kinetics in alkaline media, our study reveals that this limitation persists even in practical AWE systems where nickel-based substrates are used as electrodes. This observation is supported by both experimental data and voltage breakdown modeling. Arrhenius-type analysis reveals that localized electric fields induced by catalysts shift the reaction kinetics from classical Butler–Volmer behavior toward a Marcus-like regime, where interfacial molecular dynamics and bimolecular charge transfer dominate. We propose a semi-empirical model and a surficial reaction mechanism to describe these dynamics. This work underscores the critical need for cathode innovation and provides a rational framework for designing advanced catalysts and electrode architectures to optimize AWE performance.

08 HYDROGEN↗
DOE OSTI2025

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

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

Density Functional Theory↗
DOE OSTI2020

Polysulfone shock compressed above the decomposition threshold: Velocimetry and modeling of two-wave structures

Polysulfone was shock loaded to pressures of 14.6–26.2 GPa in a series of gas gun-driven plate-impact experiments measuring material response with embedded electro-magnetic particle velocity gauges and optical velocimetry. The embedded electro-magnetic particle velocity gauges did not show a distinct two-wave structure but did show rounding that suggested a reaction but not a distinct separation of the reactants and product waves. In contrast, the transmission experiments fielded with optical velocimetry, with product pressures ranging from 21.0 to 26.2 GPa, showed well-defined two-wave structures due to shock-driven chemical decomposition of the polymer to products at a higher density. Distinct two-wave structures have not previously been published in literature during polymer compression; here, we observed these two-wave structures at both the polymer/lithium fluoride and polymer/polymethylmethacrylate interfaces. Hydrodynamic simulations of the experiments were performed using a newly constructed SESAME equation of state (EOS) for the polymer, a thermochemical EOS for the decomposition products, and an Arrhenius reaction rate model for the kinetics of the reaction. Finally, simulation results also demonstrated two-wave structures but were unable to quantitatively reproduce either the embedded gauge or optical velocimetry data.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Hydrodynamic simulations of shock-driven chemistry in polyimide

Here, the development for equations of state both for polyimide and its reaction products is presented along with hydrodynamic simulations linking the two EOS through an Arrhenius rate law. The equations of state compare favorably to available data, and the hydrodynamic simulations are able to qualitatively reproduce many of the features seen in both legacy shock data and more modern embedded gauge data. However, quantitative agreement is not achieved, suggesting that new rate laws are needed to fully capture the effect of the reaction over a wide range of conditions.

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

Numerical simulations of the Los Alamos gapstick experiment

The gapstick is a high explosive (HE) sensitivity test recently developed at Los Alamos. The experiment design was motivated by the traditional gap test and consists of a series of HE and inert pellets in a rate stick configuration. The inert pellets are made increasingly longer so that eventually the attenuated shock is unable to initiate the next HE pellet. While the gapstick is a mechanically simple experiment it poses several challenges for numerical simulation. In particular, accurately modeling the HE initiation and detonation phenomena requires a reactive burn model with sufficient mesh resolution to capture the reaction scales. In this work, the Scaled Uniform Reactive Front (SURF) and Arrhenius Wescott-Stewart-Davis (AWSD) burn models are used for simulations of a PBX 9501 (95 wt% HMX, 5% binder) gapstick. Furthermore, the material model for the inert pellet material, Polyvinylidene Fluoride (PVDF), has a direct influence on the shock propagation and a new equation of state (EOS) calibration for PVDF is developed using available Hugoniot data. Despite some challenges, the simulations are able to reasonably predict shock transit velocities and detonation failure in the gapstick.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Experimental observations of shock-wave-induced bubble collapse and hot-spot formation in nitromethane liquid explosive

We studied the collapse of individual helium gas bubbles in the homogeneous explosive nitromethane (NM) to investigate effects of hot-spot formation on the detonation process. A bubble was injected into a NM sample, and a shock wave from an explosive detonator compressed the bubble, creating a localized hot spot. We measured shock and detonation wave speeds with optical velocimetry, and we used a high-speed camera to image the shock propagation and the pre- and post-bubble collapse processes. An infrared camera image showed the residual radiance temperature distribution after the bubble collapse, and an optical fiber pyrometer measured the time-resolved thermal radiance. We measured the optical spectra of light emitted from detonating NM without a bubble and from a collapsing bubble in shocked, undetonated NM. We estimated temperatures of the detonation fronts and of the hot spots formed by bubble collapse. To study the incipient detonation process, we performed all bubble collapse experiments at pressures below the threshold for creating a sustained detonation. Where the bubble collapsed, we observed an opaque, thermally emissive region believed to be chemical reaction products. Chemical reactions in NM can be produced with lower shock pressures (~1 GPa) when a helium bubble is present than without a bubble (~10 GPa). We used hydrodynamic modeling to predict shock wave propagation, extent of chemical reaction, and subsequent temperature rise from the collapsing bubble. Simulations using a temperature-dependent Arrhenius burn model gave much better results than reactive burn models that depend only on pressure and density.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Effect of thickness on metal-to-semiconductor transition in 2-dimensional TiN thin films

Titanium nitride thin films have been grown on c-plane sapphire substrates using a pulsed laser deposition technique in the thickness range of 6–45 nm. X-ray diffraction (XRD) analysis has demonstrated TiN (111) as the preferred orientation of growth on the sapphire substrates. The XRD measurements have also indicated that orientational alignment between the TiN and the sapphire improved with an increase in the TiN film thickness. A change in the resistivity behavior of the TiN thin films from metallic to semiconducting has been observed as the TiN film thickness is reduced below 15 nm. Analyzing and fitting of TiN films’ conductivity data have shown that while the Arrhenius law governs their conductivity in the temperature range of 300–350 K, conductivity values of the films follow the variable range hopping mechanism below 300 K.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Probabilistic effects of porosity and chemical kinetics on the shock initiation of an octahydro-1,3,5,7-tetranitro-1,3,5,7-tetrazocine (HMX) based explosive

In this work, we investigate the effects of porosity and chemical reactivity on the shock initiation of an HMX (octahydro-1,3,5,7-tetranitro-1,3,5,7-tetrazocine)-based polymer bonded explosive due to thin flyer impact using two-dimensional mesoscale simulations. Microstructural samples are computationally generated, and the calculated shock response is compared to experiments. The sample dimensions are 0.1×1.05 mm 2 , and each sample accounts for grain scale heterogeneities, including randomized grain morphology and pore spacing. Simulations are performed using the Lawrence Livermore National Laboratory (LLNL) arbitrary Lagrangian–Eulerian multi-physics code ALE3D. Chemistry and equation of states are handled with LLNL's thermochemical code Cheetah. We employ a modified Arrhenius-based chemical kinetics model that scales the temperature exponent fitting parameter with flyer shock conditions. This new model is calibrated to the mean response of experimental shock initiation thresholds for three separate Kapton flyer thicknesses. A simplified shock initiation criterion is proposed, which classifies the macroscale initiation response based on cross-sectional temperature profiles of the shocked sample. This criterion is used to further investigate the probabilistic response of the sample to variations in the mean porosity and chemical kinetics rate, and we find an inverse correlation between the shock pressure and the confidence interval of the experimental threshold velocity. Furthermore, the methodology presented in this work for quantifying the probabilistic nature of shock initiation is material agnostic. These mesoscale studies may be used to inform continuum scale models and expedite the development of future, tailored energetic formulations.

36 MATERIALS SCIENCE↗
DOE OSTI2021

A reactive flow model for the 3,3'-diamino-4,4'-azoxyfurazan based plastic bonded explosive (PBX 9701)

Here, this paper describes the calibration of the Arrhenius Wescott–Stewart–Davis (AWSD) reactive flow model for the recently developed high explosive PBX 9701, which consists of 97% 3,3'-diamino-4,4'-azoxyfurazan (DAAF) and 3% FK-800 binder by weight. DAAF-based explosive formulations have several desirable qualities as they are relatively insensitive to non-shock insults but have higher performance than triaminotrinitrobenzene based formulations. Equations of state for the explosive reactants and products are calibrated using a combination of existing and new experimental data and theoretical calculations. The AWSD rate law calibration utilizes both one-dimensional shock-to-detonation and multidimensional rate stick experiments to capture the shock initiation and propagating detonation regimes. Validation of the calibrated model is demonstrated through comparison with recent gas-gun experiments.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗
DOE OSTI2022

Effect of proton irradiation temperature on persistent photoconductivity in zinc oxide metal-semiconductor-metal ultraviolet photodetectors

The electrical and structural characteristics of 50-nm-thick zinc oxide (ZnO) metal-semiconductor-metal ultraviolet (UV) photodetectors subjected to proton irradiation at different temperatures are reported and compared. The devices were irradiated with 200 keV protons to a fluence of 1016 cm−2. Examination of the x-ray diffraction (XRD) rocking curves indicates a preferred (100) orientation prior to irradiation, with decrease in crystal quality afterward. Additionally, peak shifts in XRD and Raman spectra of the control sample relative to well-known theoretical positions are indicative of tensile strain in the as-deposited ZnO films. Shifts toward theoretical unstrained positions are observed in the irradiated films, which indicates partial relaxation. Raman spectra also indicate increase in oxygen vacancies (VO) and zinc interstitial defects (Zni) compared to the control sample. Additionally, transient photocurrent measurements performed on each sample at different temperatures showed up to 2× increase in photocurrent decay time constants for irradiated samples vs the control. This persistent photoconductive behavior is linked to the activation of electron and hole traps near the surface, and to the desorption and reabsorption of O2 molecules on the ZnO surface under the influence of UV light. Using an Arrhenius model, trap activation energies were extracted and, by comparing with known energies from the literature, the dominant defects contributing to persistent photoconductivity for each irradiation condition were identified. The persistence of differences in photocurrent transients between different samples months after irradiation indicates that the defects introduced by the suppression of thermally activated dynamic annealing processes have a long-term deleterious effect on device performance.

Heuser, Thomas A. (ORCID:0000000225642400)↗
DOE OSTI2022

Effect of void positioning on the detonation sensitivity of a heterogeneous energetic material

We show although it is well-established that voids profoundly influence the initiation and reaction behaviors of heterogeneous energetic materials such as polymer-bonded explosives (PBX) and propellants, there has been little study of how void location in different constituents in the microstructures of such materials affect the macroscale behavior. Here, we use three-dimensional (3D) mesoscale simulations to study how void placement within the reactive grains versus the polymer binder influences the shock-to-detonation transition (SDT) in a polymer-bonded explosive. The material studied here has a microstructure comprised of 75% PETN (pentaerythritol tetranitrate) grains and 25% HTPB (hydroxyl-terminated polybutadiene) polymer binder by volume. Porosities up to 10% in the form of spherical voids distributed in both the grains and polymer are considered. An Arrhenius reactive burn relation is used to model the chemical kinetics of the PETN grains under shock loading, thereby resolving the heterogeneous detonation behavior of the PBX. The influence of void location on the shock initiation sensitivity of the material is quantitatively ranked by comparing the predicted run distance to detonation (RDD) for each sample. The analysis includes inherent quantification of uncertainties arising from the stochastic variations in the microstructure morphologies and void distributions by using statistically equivalent microstructure sample sets (SEMSS), leading to probabilistic formulations for the RDD as a function of shock pressure. The calculations reveal that the location of voids in the composite microstructure significantly affects the RDD. Specifically, voids exclusively within the grains cause the PBX to be more sensitive (having shorter RDD) than voids in the polymer binder. Unique probabilistic relationships are derived to map the probability of observing RDD for each void location material case, allowing for prediction of initiation behavior anywhere in the shock pressure – RDD space. These findings agree with trends reported in the literature.

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

Experimental and computational studies of the kinetics of the reaction of hydrogen peroxide with the amidogen radical

In this work, the pulsed-laser photolysis/laser-induced fluorescence method is used to study the kinetics of the reaction of NH 2 with H 2 O 2 to yield a second-order rate constant of (2.42 ± 0.55) × 10 –14 cm 3 molecule –1 s –1 at 412 K in 10–22 mbar of Ar bath gas. There are no prior measurements for comparison. To check this value and enable reliable extrapolation to other temperatures, we also compute thermal rate constants for this process over the temperature range 298–3000 K via multi-structural canonical variational transition-state theory with small-curvature multidimensional tunneling (MS-CVT/SCT). The CVT/SCT rate constants are derived using a dual-level direct dynamics approach utilizing single-point CCSD(T)-F12b/cc-pVQZ-F12 energies—corrected for core-valence and scalar relativistic effects—and M06-2X/MG3S geometries, gradients, and Hessians—for all stationary and non-stationary points along the reaction path. The multistructural method with torsional anharmonicity, based on a coupled torsional potential, is then employed to calculate correction factors for the rate constants, accounting for the comprehensive effects of torsional anharmonicity on the kinetics of this reaction system. The final MS-CVT/SCT rate constants are found to be in good agreement with our measurements and can be expressed in modified Arrhenius form as 2.13 × 10 –15 (T/298 K) 4.02 exp(–513 K/T) cm 3 molecule –1 s –1 over the temperature range of 298–3000 K.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗