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

Molecular-level insights into structure and dynamics in ionic liquids and polymer gel electrolytes

We report designing new electrolytes requires a better understanding of the correlation between their transport properties and their molecular structure. In this work, we present a detailed study of ionic liquids and polymer gel electrolytes probing their structure and dynamics by nuclear magnetic resonance (NMR) spectroscopy. In particular, ammonium- and phosphonium-based ionic liquids (ILs) are combined with different LiTFSI concentrations, and then with different ratios of poly(methylmethacrylate) polymer. The temperature dependence of self-diffusion coefficients of mobile species, D Li+ , D TFSI- and D P4441+ or D N4441+ , measured by pulsed field gradient (PFG) NMR spectroscopy obeys the Arrhenius equation. Diffusivity of [P4441][TFSI] ILs is found to be greater than those of the [N4441][TFSI] ILs in both liquid and gel electrolytes. Solid state NMR experiments including 13 C and 19 F MAS, 13 C{ 19 F}, 13 C{ 1 H} CPMAS probe the local structure and molecular-level interactions between ions and polymer in gel electrolytes, particularly for samples with higher PMMA content (≥25 wt%). Finally, the fast field cycling relaxometry has been used to unveil the rotational and translational dynamics of P4441 + or N4441 + and TFSI - by measuring 1 H and 19 F R 1 relaxation rate profiles at different temperatures. A comprehensive NMR analysis including relaxation studies at low magnetic field provides decisive new insights regarding the formation of ionic clusters and the interaction of ions with the polymer chain in the case of the gel electrolytes.

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

Crowded electrolytes containing redoxmers in different states of charge: Solution structure, properties, and fundamental limits on energy density

Nonaqueous redox flow batteries use liquid electrolytes containing redox-active organic molecules (redoxmers) as their energy storage medium. To maximize energy density, the redoxmer concentration needs to be maximized while maintaining low viscosity and high ionic conductivity. During charge, a redoxmer molecule pairs with an ion in the electrolyte while another ion migrates across the membrane to maintain electric neutrality. In a crowded electrolyte, this reconstitution changes physical and chemical properties of the solution. To explore these behaviors, a phenothiazine redoxmer fully miscible with acetonitrile was used, and electrochemical charge was mimicked by chemical oxidation. The solutions were examined using small-angle X-ray scattering, nuclear magnetic resonance, and conductometry and modeled using classical molecular dynamics. Overall, our study indicates that physical and structural properties of redoxmer solutions in both states of charge make it exceedingly difficult to increase the redoxmer concentrations over 2 M at any temperature without compromising dynamic properties of such solutions. The cause for this limitation is proximity to a gel-like regime in which fluidity, diffusivity, and ionic conductivity exponentially decrease with increasing concentration. This tendency is compounded by non-Arrhenius behavior of the electrolyte: a small increase in the concentration outruns gains in fluidity and conductivity at a higher temperature. Thus the properties of crowded electrolytes generally make it impossible to operate when gel-like behavior sets in. Pushing the redoxmer concentration to 2.5-3 M might be possible for small redoxmer molecules, but it would require the use of ionic liquid electrolytes at 340-360 K. (C) 2021 Elsevier B.V. All rights reserved.

25 ENERGY STORAGE↗
DOE OSTI2022

Electrochemical, thermodynamic, and physical properties of tetradecyltrihexylphosphonium ([P 6,6, 6,14 ] + ) and methyl-propyl piperidinium containing ionic liquids and their propylene carbonate solutions

Here, the viscosity, conductivity, electrochemical window and related thermodynamic properties such as excess volume and dynamic viscosity deviation of two phosphonium ionic liquids were measured and calculated for “as-supplied” and dried neat liquids and also v/v mixtures of 99/1, 95/5, 90/10, 75/25, and 50/50 of the ionic liquids with propylene carbonate (PC). Tetradecyltrihexylphosphonium dicyanamide ([P 6,6,6,14 ] + dicyanamide), Cyphos 105, and tetradecyltrihexylphosphonium bis(trifluororomethanesulfonyl)imide (called bistriflimide) ([P 6,6,6,14 ] + bistriflimide), Cyphos 109, ionic liquids were studied. Generally speaking, there were slight differences in the measured properties of the wet and dried IL solutions which were reflected in differences if the calculated properties. The measured and calculated properties were compared with those collected and derived from a piperidinium-based ionic liquid, methyl-propyl piperidinium bistriflimide. Arrhenius plots for both the viscosity and conductivity were linear, with the linearity of the Litovitz plots being slightly higher. The viscosity and conductivity of the phosphonium solutions yielded Walden Plots that differed from Walden plots of piperidinium solutions and other ionic liquid solutions. For the phosphonium based ionic liquids, the ionicity was found to increase with increasing dilution in propylene carbonate. The electrochemical windows of the ionic liquids were determined as a function of the concentration of the ionic liquid present in the solution. The size of the window generally initially decreased and then increased with the addition of PC for the [P 6,6,6,14 ] + dicyanamide ionic liquid increasing from 3.4 V in the 95/5 v/v (volume Cyphos 105) to volume of PC solvent) to 3.7 in the 50/50 dried IL case. The electrochemical windows of the dried liquids were slightly larger 3.5 V for the 95/5 v/v and 4.2 V for the 50/50 v/v Cyphos 105/PC solutions. The electrochemical window was larger and the conductivity was higher for the piperidinium IL solutions, but since the ionicity increases with dilution, the phosphonium based IL solutions may prove favorable for processes such as electrochemical reduction.

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

Effects of helium cavity size and morphology on the strength of pure titanium

In this study, high-purity α-titanium was implanted with helium to observe cavity morphology – size, number density, and form (i.e., bubbles vs. voids) – effects on materials strength. Increasing implantation temperatures lead to an Arrhenius-type (i.e., exponential) increase in cavity size, a transition from spherical bubbles to faceted voids, and a marked increase in strength (~20–50%). Implanted samples deformed uniformly in contrast to nano-compression of unimplanted or alloyed titanium, but pillars with the largest faceted voids developed a type of bulging localized deformation. Affected voids were repeatedly sheared, becoming rough and shrunken. Dislocation-cavity interactions and the role of helium content within cavities are considered.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Catalytic resonance theory for the kinetic signatures of multiple wavelength photocatalysis

Multi-step surface chemistry promoted with one or two wavelengths of light exhibited distinct kinetic signatures indicative of the light-adsorbate interactions and the number of photon-sensitive elementary steps. In this work, kinetic Monte Carlo simulations identified the kinetic response of general surface mechanisms to variation in the per-site photon flux of one or two wavelengths of incident light. Photocatalytic rates were described via a non-dimensional photon flux, identifying multiple kinetic regimes with unique degrees of rate control. Under photon-controlled kinetic conditions, maximum quantum yield and turnover frequency were obtained under constant illumination, while pulsing of one or more light sources was shown to exhibit slower rates and less efficient photocatalytic promotion due to the inherent dynamic nature of light, which comprises a stream of photons. Simulations provided distinct kinetic signatures in catalytic rates and Arrhenius kinetics for specific light-surface-adsorbate interactions corresponding to photocatalytic promotion of specific steps in surface chemistry.

catalytic resonance↗
DOE OSTI2023

The P( 4 S) + NH( 3 Σ – ) and N( 4 S) + PH( 3 Σ – )reactions as sources of interstellar phosphorus nitride

Phosphorus nitride (PN) is believed to be one of the major reservoirs of phosphorus in the interstellar medium (ISM). For this reason, understanding which reactions produce PN in space and predicting their rate coefficients is important for modelling the relative abundances of P-bearing species and clarifying the role of phosphorus in astrochemistry. Here, in this work, we explore the potential energy surfaces of the $\textrm{P}(^4\textrm{S}) + \textrm{NH}(^3\Sigma^-)$ and $\textrm{N}(^4\textrm{S}) + \textrm{PH}(^3\Sigma^-)$ reactions and the formation of $\textrm{H}(^2\textrm{S}) + \textrm{PN}(^1\Sigma^+)$ through high accuracy ab initio calculations and the variable reaction coordinate transition state theory (VRC-TST). We found that both reactions proceed without an activation barrier and with similar rate coefficients that can be described by a modified Arrhenius equation ( $k(T)=\alpha\!\left( T/300 \right)^{\beta} \exp\!{(\!-\!\gamma/T)})$ with $\alpha=0.93\times 10^{-10}\rm cm^3\,s^{-1}$ , $\beta=-0.18$ and $\gamma=0.24\, \rm K$ for the $\textrm{P} + \textrm{NH} \longrightarrow \textrm{H} + \textrm{PN}$ reaction and $\alpha=0.88\times 10^{-10}\rm cm^3\,s^{-1}$ , $\beta=-0.18$ and $\gamma=1.01\, \rm K$ for the $\textrm{N} + \textrm{PH} \longrightarrow \textrm{H} + \textrm{PN}$ one. Both reactions are expected to be relevant for modelling PN abundances even in the cold environments of the ISM. Given the abundance of hydrogen in space, we have also predicted rate coefficients for the destruction of PN via H + PN collisions.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Ensembles of Metastable States Govern Heterogeneous Catalysis on Dynamic Interfaces

Heterogeneous catalysis is at the heart of chemical industry. Being able to tune and design efficient catalysts for processes of interest is of an utmost importance, and for this, the molecular-level understanding of heterogeneous catalysts is the first step, and indeed a prime focus of the modern catalysis research. For a long time, a single most thermodynamically stable structure of the catalytic interface attained in reaction conditions had been envisioned as the reactive phase. However, some catalytic interfaces continue to undergo structural dynamics in the steady state, triggered by high temperatures, pressures, and binding and changing reagents. Among particularly dynamic interfaces are such widely-used catalysts as crystalline and amorphous surfaced supporting (sub-)nano metallic clusters. Recently, it became clear that this dynamic fluxionality causes the supported clusters to populate many distinct structural and stoichiometric states in catalytic conditions. Hence, the catalytic interface should be viewed as an evolving statistical ensemble of many (not one) structures. As such, every member in the ensemble contributes to the properties of the catalyst differently, and in proportion to its probability to be populated. This new notion flips the established paradigm and calls for new theory, modeling approaches, operando measurements, and updated design strategies. The statistical ensemble nature of surface-supported sub-nano cluster catalysts can be exemplified by oxide-supported and adsorbate-covered Pt, Pd, Cu, CuPd clusters, catalytic toward oxidative and non-oxidative dehydrogenation. They have access to a variety of 3D and quasi-2D shapes. The compositions of their thermal ensembles are dependent on the cluster size, leading to size-specific catalytic activities and the famous “every atom counts” phenomenon. The support and adsorbates affect catalyst structures, and state of the reacting species causes the ensemble to change in every reaction intermediate. The most stable member of the ensemble dominates the thermodynamic properties of the corresponding intermediate, whereas the kinetics can be determined by more active but less populated metastable catalyst states, and that suggests that many earlier studies might have overlooked the actual active sites. Both effects depend on the relative time-scales of catalyst restructuring and reaction dynamics. The catalyst may routinely operate off-equilibrium. Ensemble phenomena lead to surprising exceptions from established rules of catalysis, such as scaling relations, and Arrhenius behavior. Catalyst deactivation is also an ensemble property, and its extent of mitigation can be predicted through the new paradigm. These findings were enabled by advances in theory, such as global optimization and subsequent utilization of multiple local minima, and pathways sampling, as well as operando catalyst characterization. The fact that the per-site and per-species resolution is needed for the description and predicting of catalyst properties gives theory the central role in catalysis research, as most experiments provide ensemble-average information and cannot detect the crucial minority species that may be responsible for the catalytic activity.

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

Solid-State, Single-Ion Conducting, Polymer Blend Electrolytes with Enhanced Li + Conduction, Electrochemical Stability, and Limiting Current Density

The development of solid-state polymer electrolytes with high lithium conductivity is crucial to improve lithium-ion battery performance and ameliorate the safety challenges associated with current solvent-based electrolytes. Unfortunately, sluggish polymer segmental dynamics are known to constrain conductivity enhancements in solid-state polymer electrolyte systems, limiting overall performance. In this work, a glassy single-ion-conducting polymer, poly[lithium sulfonyl(trifluoromethane sulfonyl)imide methacrylate] (PLiMTFSI), was blended with a flexible polymer, poly(oligo-oxyethylene methyl ether methacrylate) (POEM), and the impact of PLiMTFSI molecular weight and ion concentration on the thermal and ion-conducting behavior of blend electrolytes was investigated. High ionic conductivities approaching 1 × 10 -2 S/cm at 150 °C were realized in this polymer blend electrolyte system as a result of decoupling Li + transport from polymer segmental dynamics. The decoupled ion transport was attributed to the packing frustration of the glassy PLiMTFSI – sufficient percolating free volume was generated to produce effective ion diffusion pathways. This decoupling was tunable as the ion transport could be altered from being closely coupled to the polymer segmental dynamics (Vogel–Tammann–Fulcher-like) to hopping (Arrhenius-like) by increasing the PLiMTFSI molecular weight and ion concentration. Moreover, the immobilized TFSI anion resulted in high Li + selectivity (Li + transference number = 0.9), high electrochemical stability (up to 4.7 V against Li + / Li), and limiting current density of 1.8 mA/cm 2 (electrolyte thickness = 0.05 cm). These features suggest that this single-ion-conducting, polymer blend electrolyte might be a promising alternative to a benchmark system – salt-doped polyethylene oxide. Moreover, the above characteristics can support the battery operation at higher voltages using energy-dense Li metal anodes, with faster charging rates and enhanced energy/power densities. Altogether, the results suggest that polymer chain packing frustration can be exploited to overcome the constraints of slow polymer segmental relaxations to achieve rapid and highly selective ion transport and enhanced performance in solid-state polymer electrolytes.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Using Data-Science Approaches to Unravel Insights for Enhanced Transport of Lithium Ions in Single-Ion Conducting Polymer Electrolytes

Solid polymer electrolytes have yet to achieve the desired ionic conductivity (>1 mS/cm) near room temperature required for many applications. This target implies the need to reduce the effective energy barriers for ion transport in polymer electrolytes to around 20 kJ/mol. In this work, we combine information extracted from existing experimental results with theoretical calculations to provide insights into ion transport in single-ion conductors (SICs) with a focus on lithium ion SICs. Through the analysis of temperature-dependent ionic conductivity data obtained from the literature, we evaluate different methods of extracting energy barriers for lithium transport. The traditional Arrhenius fit to the temperature-dependent ionic conductivity data indicates that the Meyer–Neldel rule holds for SICs. However, the values of the fitting parameters remain unphysical. Our modified approach based on recent work (Macromolecules 2023, 56, 15, 6051), which incorporates a fixed pre-exponential factor, reveals that the energy barriers exhibit temperature dependence over a wide range of temperatures. Using this approach, we identify anions leading to the energy barriers <30 kJ/mol, which include trifluoromethane sulfonimide (TFSI), fluoromethane sulfonimide (FSI), and boron-based organic anions. In our efforts to design the next generation of anions, which can exhibit the energy barriers <20 kJ/mol, we have performed density functional theory (DFT) based calculations to connect the chemical structures of boron-based anions via the binding energy of cation (lithium)-anion pairs with the experimentally derived effective energy barriers for ion hopping. Not only have we identified a correlation between the binding energy and the energy barriers, but we also propose a strategy to design new boron-based anions by using the correlation. This combined approach involving experiments and theoretical calculations is capable of facilitating the identification of promising new anions, which can exhibit ionic conductivity >1 mS/cm near room temperature, thereby expediting the development of novel superionic single-ion conducting polymer electrolytes.

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

Expanding Configurational Complexity through Dipole Dilution in Pseudohalide Argyrodite Ion Conductors

The advantageous properties of (pseudo)halide argyrodite ion conductors of the formula Li 6 PS 5 X (X = Cl – , Br – , I – , CN – ) have motivated extensive studies of their structure-transport relationships, particularly as they pertain to the role of atomic site disorder. The argyrodite structure can accommodate additional configurational complexity to promote ion transport via extended three-anion site mixing and the potential for orientational disorder of molecular anions. In this work, we explore a ternary anion system including the cyanide anion, expanding site disorder and introducing dipolar orientations as an additional degree of freedom. We prepared the series Li 6 PS 5 (CN) 1–x Br x , in which the dipolar cyanide anions are systematically diluted with bromide. We find that anion disorder, as determined by synchrotron and neutron diffraction and quantified by configurational entropy (S config ), is correlated with lowered activation barriers and increased lithium ion conductivity. We propose that S config describes the electrostatic heterogeneity of the Li environments, flattening the energetic landscape for ion transport. While anion substitution strongly impacts the activation barrier for transport, the temperature-independent Arrhenius prefactor does not follow the same trend. Through heat-capacity measurements of attempt frequency and deconvolution of terms within the prefactor, we rationalize the apparent decoupling of activation energy and prefactor to strong cyanide-lithium interactions that increase the entropy of migration. Together, these results expand the structure–property relationships in the argyrodite family to encompass multiple facets of disorder and the subsequent impact on lithium ion transport.

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

Isoconversional Kinetic Analysis of Oxygen Desorption for Sr 0.75 Ca 0.25 FeO 3-δ Perovskite in Dry and Steam-Based Environments

The desorption kinetics of Sr 0.75 Ca 0.25 FeO 3-δ perovskite material are examined in both dry and steam-based environments using redox gaseous products of a laboratory-scale fixed bed. First, desorption kinetics associated with the dry environment are proposed based on Friedman’s isoconversional method. It is found from the reconstructed reaction model that the desorption kinetics of this perovskite are controlled by a three-step mechanism, where hypothetically, the phase-boundary reaction rapidly prevails first, followed by diffusion and/or nucleation, and finally, unimolecular decay or random nucleation-controlled reaction. Subsequent analyses of the Arrhenius parameters indicated that both the early and later desorption stages, respectively, controlled by the phase boundary and unimolecular/random nucleation reactions, are associated with low energetic demand. In contrast, a high energetic penalty is required for the diffusion/nucleation reaction. Further, a satisfactory a priori verification of the proposed kinetics suggested that the three-step reaction mechanism reasonably describes the desorption of this perovskite in a dry environment. With the presence of steam in the desorption environment, oxygen production is substantially inhibited, but the hypothetical three-step mechanism is still found to control the desorption. Meanwhile, the primary effects of the steam on these mechanisms include the widening of the conversion extent associated with the diffusion/nucleation reaction, along with a higher energetic demand compared to the dry environment.

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

Microwave-Assisted Plastic Upcycling: Dynamic Data Reconciliation, Parameter Estimation, and Kinetic Modeling

Microwave (MW)-assisted catalytic pyrolysis offers a promising pathway for efficient plastic upcycling. This work develops an integrated modeling framework combining dynamic data reconciliation, a temperature-dependent rate model, and a yield model to represent the time-varying production rate of components in MW-assisted LDPE pyrolysis conducted in a batch reactor. An Arrhenius-type rate model with a temperature-dependent reaction order is developed. A biexponential correlation is proposed for the yield of gaseous products that enables to capture the evolving product formation behavior during conversion. In the yield correlation, one term is used to represent the initial increase in yield, reflecting the rapid formation of intermediate or primary products at the early stages of the reaction when a larger fraction of the reactant remains available. As conversion progresses, the influence of this term gradually diminishes. The other term accounts for the subsequent decrease in the predicted yield, representing secondary reactions such as further cracking or coke formation that reduce the concentration of certain products at higher conversion. The model is found to accurately represent reconciled experimental flow rate profiles from an in-house MW-assisted catalytic batch reactor for major products, including ethylene, ethane, 1-butene, and benzene, across 250−350 °C. Ethylene remains the dominant product but decreases from about 41.95% at 250 °C to 30.14% at 350 °C, while heavier products increase significantly, with 1-butene rising to nearly 8.37% and benzene reaching 2.17% at intermediate temperatures. The model shows that the ethylene production rate can be maximized at around 270 °C. The models developed in this work can be utilized for process optimization, reactor design and scale-up of microwave-assisted plastic conversion technologies, and economic analysis.

Damahe, Harish [West Virginia Univ., Morgantown, W↗
DOE OSTI2025

La 3 CuTe 5 : A Narrow-Gap Semiconductor with Indirect Gap and Dual-Regime Thermally Activated Transport

Metal-chalcogenide systems remain a long-standing research topic because of their structural diversity and potential to host emergent phenomena. Here, we report a new compound, La 3 CuTe 5 , synthesized from the halide-flux method. Single-crystal X-ray diffraction studies indicate the structure to be unique among reported ones. Here, the compound crystallizes in a novel structure type adopting the orthorhombic space group Pnma and a unit cell of a = 24.3947(14) Å, b = 4.4232(2) Å, and c = 10.2142(5) Å. The tetrahedral [CuTe 4 ] building blocks form chains along [010] by corner sharing and link [LaTe 7 ] and [LaTe 8 ] polyhedra via edge sharing, resulting in a three-dimensional bulk structure. Thermal analysis results indicate that the material remains stable with a temperature up to 950 °C and decomposable at 1400 °C. First-principles calculations reveal an indirect electronic band gap and flat valence bands dominated by Te p and Cu d states. Optical absorption measurements yield a band gap of ∼0.65 eV, consistent with semiconducting behavior observed in transport measurements. Fittings to the temperature-dependent resistivity reveal two thermally activated regimes associated with Arrhenius-type conduction and three-dimensional variable range hopping, respectively.

Chalcogenides↗
DOE OSTI2021

Transformation of an Embedded Five-Membered Ring in Polycyclic Aromatic Hydrocarbons via the Hydrogen-Abstraction–Acetylene-Addition Mechanism: A Theoretical Study

Five-membered rings are constituents of many polycyclic aromatic hydrocarbons (PAHs) and their presence on the edges of large PAHs has been repeatedly observed experimentally. However, modern kinetic combustion models often do not consider the growth of PAHs through the transformation of the five-membered rings. Here, in connection with the above, we carried out a theoretical study of the mechanism of Hydrogen-Abstraction-aCetylene-Addition (HACA) transformation of an embedded five-membered ring on the armchair PAH edge to a six-membered ring, considering cyclopenta[d,e,f]phenanthrene (4,5-methylenephenanthrene) as a prototype system for this process. The potential energy surface for the reactions of cyclopenta[d,e,f]phenanthrenyl radicals produced by direct H abstractions from cyclopenta[d,e,f]phenanthrene with acetylene has been compiled at the G3(MP2,CC)//B3LYP/6-311G(d,p) level of theory including zero-point vibrational energy corrections. The computed energies and molecular parameters were then used to solve the Rice-Ramsperger-Kassel-Marcus Master Equation (RRKM-ME) in order to calculate reaction rate at various pressures and temperatures, which were fitted to modified Arrhenius equation for further kinetic modeling. The results show that the HACA transformation of the embedded five-membered ring to a six-membered ring is possible, albeit slow. The most viable reaction mechanism involves the R2 + C 2 H 2 reaction, where the acetylene molecules adds to a σ-radical in the six-membered ring adjacent to the five-membered ring via a low entrance barrier. The predominant product of R2 + C 2 H 2 is predicted to be 3 ethynyl 4H cyclopenta[def]phenanthrene Pr5 via immediate H elimination from the initial addition complex. Next, Pr5 undergoes H-assisted isomerization to 4aH-pentaleno[4,3,2,1-cdef]phenanthrene Pr4, and the latter adds an H atom eventually forming the 1-pyrenylmethyl radical Pr3: R2 + C 2 H 2 ⇌ 3 ethynyl 4H cyclopenta[def]phenanthrene (Pr5) + H ⇌ 4aH-pentaleno[4,3,2,1-cdef]phenanthrene (Pr4) + H Pr5 + H ⇌ Pr4 + H Pr4 + H → 1-pyrenylmethyl (Pr3) This HACA sequence may be competitive with the methyl radical addition to the R1 radical formed by H abstraction from the CH 2 group in the five-membered ring of cyclopenta[d,e,f]phenanthrene, which provides a pathway to pyrene following two H atom losses. Relative contributions of the two mechanisms of the five- to six-membered ring transformation would strongly depend on the branching ratios of the R1 and R2 radicals produced by the H abstractions and the available concentration of C 2 H 2 vs. CH 3 and hence differ in different flames.

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

The Role of Methylaryl Radicals in the Growth of Polycyclic Aromatic Hydrocarbons: The Formation of Five-Membered Rings

The regions of the C 13 H 11 potential energy surface (PES) related to the unimolecular isomerization and decomposition of the 1-methylbiphenylyl radical and accessed by the 1-/2-methylnaphthyl + C 2 H 2 reactions have been explored by ab initio G3(MP2,CC)//B3LYP/6-311G(d,p) calculations. The kinetics of these reactions relevant to the growth of polycyclic aromatic hydrocarbons (PAH) under high-temperature conditions in circumstellar envelopes and in combustion flames has been studied employing the RRKM-Master Equation approach. The unimolecular reaction of 1-methylbiphenylyl proceeding via a five-membered ring closure followed by H elimination is predicted to be very fast, on a submicrosecond scale above 1000 K and to result in the formation of an embedded five-membered ring in the 9H-fluorene product. The 1-/2-methylnaphthyl + C 2 H 2 reaction mechanism involves acetylene addition to the radical on the methylene group followed by a six- or five-membered ring closure and aromatization via an H atom loss. Despite of the complexity of the C 13 H 11 PES, these straightforward pathways are dominant in the high-temperature regime (above ~1000 K), with the prevailing products being phenalene, with a significant contribution of 1H-cyclopenta(a)naphthalene, for 1-methylnaphthyl + C 2 H 2 , and 1H-cyclopenta(b)naphthalene and 3H-cyclopenta(a)naphthalene, for 2-methylnaphthyl + C 2 H 2 . The methylnaphthyl reactions with acetylene represent a clean source of the three-ring PAHs, but they are relatively slow owing to the high entrance barriers of ~10 kcal/mol, with the rate constants of about an order of magnitude lower as compared to those for naphthyl + allene and σ-aryl + C 2 H 2 . Here, the 1-methylnaphthyl + C 2 H 2 and 2-methylnaphthyl + C 2 H 2 reactions represent prototypes for PAH growth by an extra six- and five-membered ring on a zigzag edge or a corner of PAH and the generated modified Arrhenius expressions are recommended for kinetic modeling of PAH expansion by the mechanism of acetylene addition to methylaryl radicals.

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

High-Temperature Reaction Kinetics of the e aq – and HO 2 • Radicals with Iron(II) Ions in Aqueous Solutions

Pulsed electron radiolysis was used to determine the chemical reaction kinetics and Arrhenius parameters for iron(II) reactions in aqueous solutions under irradiation. The second-order Fe 2+ reactions with the hydrated electron (e aq – ) and the perhydroxyl radical (HO 2 • ), arising from water radiolysis, were measured to high temperatures using custom-built flow-through cells with a multichannel optical detection system. The reaction with the HO 2 • radical was found to proceed via the formation of a metal-ion adduct species, Fe 2+ –HO 2 • . Additionally, the adduct’s molar extinction coefficient and its first-order decay rate coefficients are reported.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2023

Bimolecular Reaction of Methyl-Ethyl-Substituted Criegee Intermediate with SO 2

Methyl-ethyl-substituted Criegee intermediate (MECI) is a four-carbon carbonyl oxide that is formed in the ozonolysis of some asymmetric alkenes. MECI is structurally similar to the isoprene-derived methyl vinyl ketone oxide (MVK-oxide) but lacks resonance stabilization, making it a promising candidate to help us unravel the effects of size, structure, and resonance stabilization that influence the reactivity of atmospherically important, highly functionalized Criegee intermediates. We present experimental and theoretical results from the first bimolecular study of MECI in its reaction with SO 2 , a reaction that shows significant sensitivity to the Criegee intermediate structure. Using multiplexed photoionization mass spectrometry, we obtain a rate coefficient of (1.3 ± 0.3) × 10 –10 cm 3 s –1 (95% confidence limits, 298 K, 10 Torr) and demonstrate the formation of SO 3 under our experimental conditions. Through high-level theory, we explore the effect of Criegee intermediate structure on the minimum energy pathways for their reactions with SO 2 and obtain modified Arrhenius fits to our predictions for the reaction of both syn and anti conformers of MECI with SO 2 (k syn = 4.42 × 10 11 T –7.80 exp(-1401/T) cm 3 s –1 and k anti = 1.26 × 10 11 T –7.55 exp(-1397/T) cm 3 s –1 ). In conclusion, our experimental and theoretical rate coefficients (which are in reasonable agreement at 298 K) show that the reaction of MECI with SO 2 is significantly faster than MVK-oxide + SO 2 , demonstrating the substantial effect of resonance stabilization on Criegee intermediate reactivity.

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

High-Pressure Rate Rules for Ether Alkylperoxy Radical Isomerization

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

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗