Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “phase kinetics”

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

Nickel hydroxide–nickel carbonate competitive growth on carbonate surfaces

The thermodynamic and kinetic factors controlling the competitive heterogeneous nucleation and growth of ubiquitous metal carbonate and hydroxide phases are poorly understood. In this work, calcite (CaCO 3 ) and magnesite (MgCO 3 ) powders were reacted with NiCl 2 (0–600 μM) for 7 days at 22 °C and 5 °C. The reacted powders were analyzed with X-ray photoelectron spectroscopy (XPS), scanning electron microscopy, and energy-dispersive X-ray spectroscopy to characterize the Ni surface precipitates formed. Evidence from these techniques pointed to the formation of mixed Ni carbonate-Ni hydroxide amorphous surface precipitates. On calcite, XPS detected primarily Ni(OH) 2 despite the initial solutions being more supersaturated with respect to gaspéite (NiCO 3 ) than to theophrastite (Ni(OH) 2 ) by a factor of 17–18. In contrast, NiCO 3 was the dominant component detected by XPS on magnesite in the same conditions. Decreasing the temperature had the effect of increasing the proportion of NiCO 3 to the detriment of Ni(OH) 2 . The experimental observations were consistent with low lattice/cation size mismatch favoring NiCO 3 nucleation and temperature most influencing Ni(OH) 2 nucleation. Comparison to previous work on Co-reacted powders indicated the differences in lattice/cation size mismatch and/or water exchange rate impacted the composition of the surface precipitates more than the relative thermodynamic stabilities of the competing minerals. This work explored the heterogeneous growth regime of Ni carbonate and hydroxide phases on carbonate surfaces and shed light on the factors that control the competition between surface precipitates when mineral surfaces are in contact with aqueous solutions supersaturated with respect to multiple mineral phases. In conclusion, these results contribute to geochemists’ efforts toward interpreting data from geochemical systems with elevated Ni concentrations, improving Ni environmental remediation and recovery strategies, and predicting the fate and transport of Ni in geochemical systems.

Heterogeneous nucleation↗
DOE OSTI2021

Predicting zeolites’ stability during the corrosion of nuclear waste immobilization glasses: Comparison with glass corrosion experiments

During the long-term corrosion of nuclear waste glasses under nuclear waste disposal conditions, the precipitation of zeolitic phases has been linked to a delayed acceleration in glass corrosion (known as “Stage III”). Hence, predicting the thermodynamic propensity for zeolites to form upon the dissolution of nuclear waste glasses is key to ensure their long-term performance. Here, we compile a unified, internally-consistent thermodynamic database “clay20” to estimate the stability of clay and feldspar phases relevant to nuclear waste immobilization glasses, including beidellite(Mg, Ca, Na, K), kaolinite, montmorillonite(Mg, Ca, Na, K), nontronite(Mg, Ca, Na, K), saponite(Ca, Na, K), and albite. Based on this, we report a geochemical modeling method allowing us to predict the stability of secondary phases (including zeolites, calcium–silicate–hydrate gels, and clays) upon the dissolution of nuclear waste immobilization glasses. We show that this approach offers a realistic description of the stability of the secondary phases forming during the dissolution of two archetypical model nuclear glasses (namely, the International Simple Glass, ISG, and WVUTh-203) under conditions relevant to nuclear waste disposal (T = 90°C, p = 1 bar) as a function of pH. We find that the formation of silica and clay secondary phases is thermodynamically favored at low pH (pH < 10), whereas, in contrast, zeolite (analcime) and calcium–silicate–hydrate phases are favored at high pH (pH > 10.5). This suggests that thermodynamics (i.e., not solely kinetics) plays a key role in determining the range of solution pH wherein stage III corrosion may occur, i.e., when zeolite formation is favored.

Zhen-Wu, Bi Yun↗
DOE OSTI2019

A Review on Lignin Liquefaction: Advanced Characterization of Structure and Microkinetic Modeling

Lignin liquefaction microkinetics is a move toward a more first principles (i.e., ab initio) based understanding at the molecular level in reaction engineering. While the microkinetic modelling of reactions to obtain kinetic rate parameters of chemical reactions have been widely used in the field of gas phase combustion and heterogeneous catalysis, this approach has not been as thoroughly developed in the area of biomass thermochemical reactions (e.g., lignin pyrolysis, hydrothermal liquefaction). The difficulties in establishing the structure of complex heterogeneous materials, like lignin, is perhaps the main challenge in developing rational microkinetic descriptions of biomass thermochemical reactions. In this manuscript we review the current state of the art and the challenges to develop microkinetic models for lignin liquefaction technologies (e.g., pyrolysis, hydrothermal liquefaction, solvolysis). A general strategy for the development of microkinetic models for lignin liquefaction technologies is discussed. Here, the first hurdle is to obtain sufficiently rich experimental data of lignin underlying polymeric structure and methodologies to use this data to build realistic lignin structural representations. Some analytical techniques for lignin structural characterization and their associated data, as well as a correlation for calculating the degree of macromolecular lignin branching are discussed. The presence of small lignin oligomeric structures and the role of these structures in lignin pyrolysis is also addressed. The ways in which elementary deconstruction and repolymerization reactions occur within this structure to form a liquid intermediate, and how these deconstruction products continue to interact with each other until they are removed from the liquid intermediate is thoroughly discussed. Further, experimental work with model compounds and the effect of reaction parameters (e.g., temperature, pressure, vapor residence time) are reviewed.

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

Regeneration of Active Surface Alloys during Cyclic Oxidation and Reduction: Oxidation of H 2 on Pd/Ag(111)

The surface morphology and composition of a catalyst during excursions between oxidizing and reducing conditions can change substantially, especially in bimetallic alloys. Both thermodynamic and kinetic factors play a role in determining the properties of alloy surfaces where the active phase may be a metastable state. Previously, Ag oxide reduction was shown to be dramatically enhanced when Pd is on the surface; however, Pd is more stable when dissolved in Ag, raising the question as to whether a highly active Pd surface state will persist over multiple reaction cycles—a requirement for catalytic function. Experiments herein demonstrate that the enhanced chemical functionality due to the presence of Pd on the surface is retained, based on the enhanced rate of silver oxide reduction over multiple oxidation/reduction cycles for a Pd/Ag(111) model. Repeated oxidation and reduction promote PdAg alloying and reversible structural and compositional changes are detected using X-ray photoelectron spectroscopy. Furthermore, this study establishes that metastable phases can persist in reactive processes on surfaces, indicating their potential in heterogeneous catalysis.

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

Microhydration of the metastable N -Protomer of 4-Aminobenzoic acid by condensation at 80 K: H/D exchange without conversion to the more stable O- protomer

4-Aminobenzoic acid (4ABA) is a model scaffold for studying solvent-mediated proton transfer. Although protonation at the carboxylic group (O-protomer) is energetically favored in the gas phase, the N-protomer, where the proton remains on the amino group, can be kinetically trapped by electrospray ionization of 4ABA in an aprotic solvent such as acetonitrile. Here we report the formation of the hydrated deuterium isotopologues of the N-protomers, RND 3 + ·(H 2 O) n=1-3 , (R=C 6 H 4 COOD), which are generated by condensing water molecules onto the bare N-protomers in a liquid nitrogen cooled, radiofrequency octopole ion trap at 80 K. The product clusters are then transferred to a 20 K cryogenic ion trap where they are tagged with weakly bound D 2 molecules. The structures of these clusters are determined by analysis of their vibrational patterns obtained by resonant IR photodissociation. The resulting patterns confirm that the metastable N-protomer configuration remains intact even when warmed by sequential condensation of water molecules. Attachment of H 2 O molecules onto the RND 3 + head group also affords the opportunity to explore the possibility of H/D exchange between the acid scaffold and the proximal water network. The spectroscopic results establish that although the RND 3 + ·(H 2 O) n=1,2 clusters are formed without H/D exchange, the n = 3 cluster exhibits about 10% H/D exchange as evidenced by the appearance of the telltale HOD bands. Furthermore, the site of exchange on the acid is determined to be the acidic OH by the emergence of the OH stretching fundamental in the -COOH motif.

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

Schwinger-Keldysh effective action for hydrodynamics with approximate symmetries

We study hydrodynamic theories with approximate symmetries in the recently developed effective action approach on the Schwinger-Keldysh contour. We employ the method of spurious symmetry transformation for small explicit symmetry-breaking parameters to systematically constrain symmetry-breaking effects in the nonequilibrium effective action for hydrodynamics. We apply our method to the hydrodynamic theory of chiral symmetry in quantum chromodynamics at finite temperature and density and its explicit breaking by quark masses. We show that the spurious symmetry and the Kubo-Martin-Schwinger relation dictate that the Ward-Takahashi identity for the axial symmetry, i.e., the partial conservation of axial vector current (PCAC) relation, contains a relaxational term proportional to the axial chemical potential, whose kinetic coefficient is at least of the second order in the quark mass. In the phase where the chiral symmetry is spontaneously broken, and the pseudo-Nambu-Goldstone pions appear as hydrodynamic variables, this relaxation effect is subleading compared to the conventional pion mass term in the PCAC relation, which is of the first order in the quark mass. On the other hand, in the chiral symmetry restored phase, we show that our relaxation term, which is of the second order in the quark mass, becomes the leading contribution to the axial charge relaxation. Therefore, the leading axial charge relaxation mechanism is parametrically different in the quark mass across a chiral phase transition.

Goldstone bosons↗
DOE OSTI2022

Quantum Phases of Transition Metal Dichalcogenide Moiré Systems

Moiré systems provide a rich platform for studies of strong correlation physics. Recent experiments on heterobilayer transition metal dichalcogenide Moiré systems are exciting in that they manifest a relatively simple model system of an extended Hubbard model on a triangular lattice. Inspired by the prospect of the hetero-transition metal dichalcogenide Moiré system’s potential as a solid-state-based quantum simulator, we explore the extended Hubbard model on the triangular lattice using the density matrix renormalization group. Specifically, we explore the two-dimensional phase space spanned by the key tuning parameters in the extended Hubbard model, namely, the kinetic energy strength and the further-range Coulomb interaction strengths. We find competition between Fermi fluid, chiral spin liquid, spin density wave, and charge order. In particular, our finding of the optimal further-range interaction for the chiral correlation presents a tantalizing possibility.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Models of dislocation glide and strengthening mechanisms in bcc complex concentrated alloys

Abstract The mechanical response of complex concentrated alloys (CCAs) deviates from that of their pure and dilute counterparts due to the introduction of a combinatorially sized chemical concentration dimension. Compositional fluctuations constantly alter the energy landscape over which dislocations move, leading to line roughness and the appearance of defects such as kinks and jogs under stress and temperature conditions where they would ordinarily not exist in pure metals and dilute alloys. The presence of such chemical defects gives rise to atomic-level mechanisms that fundamentally change how CCAs deform plastically at meso- and macroscales. In this article, we provide a review of recent advances in modeling dislocation glide processes in CCAs, including atomistic simulations of dislocation glide using molecular dynamics, kinetic Monte Carlo simulations of edge and screw dislocation motion in refractory CCAs, and phase-field models of dislocation evolution over complex energy landscapes. We also discuss pathways to develop comprehensive simulation methodologies that connect an atomic-level description of the compositional complexity of CCAs with their mesoscopic dislocation-mediated plastic response with an eye toward improved design of CCA with superior mechanical response. Graphical abstract

36 MATERIALS SCIENCE↗
DOE OSTI2022

Developing A New NG Super-Absorbent Polymer (NG-SAP) for a Practical NG Storage System with Low Pressure, Ambient Temperature, and High Energy Density

This R&D project proposed a systematical study to develop a new class of NG super-absorbent polymer (NG-SAP) with the objective to achieve similar CNG volumetric storage level ~250 cm 3 -STP/cm 3 and fast sorption-desorption kinetics at ambient temperature but lower operational pressure (<100 bar). We suggested two consecutive research phases, including (i) a systematical investigation of various polymer systems to identify the suitable classes that can exhibit methane (C1) binding energy in the range of 20-25 kJ/mol and then (ii) focusing on the potential class polymers to identify the specific structure and morphology that can achieve the project objective.In first research phase (1st year), the major tasks included design and synthesis of several representative polymers with various structure attributes, examining their molecular structures, pore morphology, and BET surface area as the prepared samples, and then studying their C1 binding energy and sorption capacity.

03 NATURAL GAS↗
DOE OSTI2020

Lower length scale informed improvements to Bison U-Pu-Zr fuel swelling model

Due to renewed interest in metallic fuels from the U-(Pu)-Zr material system, significant improvements have recently been made to BISON's capability to simulate metallic fuel, particularly in the area of swelling. In this report, lower length scale simulations that have been conducted to inform the engineering-scale models in BISON during FY20 are described. For the high-temperature regions of the fuel that are dominated by the $\gamma$ phase of U-(Pu)-Zr, two new models for swelling have implemented in BISON. Phase-field simulations were conducted to calculate parameters for these and other BISON models that control when gaseous swelling ceases and fission gas release begins. For the low-temperature regions of the fuel dominated by the phase of uranium, microstructurally resolved simulations of pore growth and crystal deformation were performed to understand the effect of irradiation-generated crystal shape changes. Conformally meshed microstructures with pores were simulated with a crystal plasticity model and with a burnup eigenstrain model. Plastic deformation will not cause pore size to increase; burnup eigenstrain (irradiation-induced crystal shape changes) will increase the volume of existing porosity and cause it to become anisotropic in shape. Additional model development work for early-stage fission gas cluster/bubble formation in the $\gamma$ phase were performed, by coupling with defect clustering within cluster dynamics framework. Sensitivity analysis were conducted to identify the critical thermo-kinetic parameters needed for rigorous assessment.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

Kinetics and Reaction Mechanisms of Acetic Acid Hydrodeoxygenation over Pt and Pt–Mo Catalysts

In this work, kinetic measurements for silica-supported Pt and Pt-Mo catalysts were collected in vapor-phase acetic acid hydrodeoxygenation by varying the hydrogen partial pressure between 18 and 72 kPa and acetic acid partial pressure between 7 and 18 kPa at 423-473 K. At all testing conditions, the Pt-Mo catalyst was more active and selective. In addition, the apparent activation energy for Pt-Mo of 76±3 kJ/mol was lower than that of 84±4 kJ/mol for Pt. The apparent reaction orders were also different. The order in hydrogen of 0.8±0.1 for Pt-Mo changed to zero at higher hydrogen partial pressures while that for Pt remained constant at 0.6±0.1. A near-zero order in acetic acid for Pt-Mo changed to -2.1 at higher acetic acid pressures while that for Pt remained constant at -2.9±0.3. These differences in reaction kinetics as well as in selectivity trends with changes in temperature and feed composition indicated a change in the reaction mechanism for Pt-Mo. The catalysts were characterized with hydrogen temperature programmed desorption, oxygen temperature programmed oxidation and transmission electron microscopy with energy-dispersed X-ray spectroscopy elemental mapping. Mo was present in the form of subnanometer-size clusters on the surface of Pt nanoparticles. Both Pt and Pt-Mo catalysts were stable under the reaction conditions for 10 h, and the size and structure of Pt and Pt-Mo particles remained mostly unchanged, without coke accumulation. Density functional theory calculations show that surface acetate is not a major reaction intermediate on both Pt and Pt-Mo and, instead, C-OH bond splitting with the formation of acetyl is the first reaction step in hydrodeoxygenation. The activation energy for this step is dramatically lower on Pt-Mo. Furthermore, the activity of acetyl on Pt-Mo is different. Unlike on pure Pt, the reaction of C-O bond splitting becomes exothermic with a lower activation energy on Pt-Mo, and the reaction of C-C bond splitting, in contrast, becomes endothermic with a higher activation energy, explaining the experimentally observed higher activities and selectivities of Pt-Mo. In addition, the calculations demonstrate that a pair of Pt-Mo surface atoms acts as a single active site where Mo serves as a preferential binding anchor for O atoms. The presence of Mo atoms changes the structure and energy of adsorbed and reacting surface species.

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

Multiple freezing–melting pathways of high-density ice through ice XXI phase at room temperature

Various metastable ice phases and their complicated transition pathways have been found by pressurization at low temperatures at which slow kinetics and high metastability are easily achieved. By contrast, such diversity is less expected at room or elevated temperatures. Here, using a combination of a dynamic diamond anvil cell and X-ray free electron laser techniques, we demonstrate that supercompressed water transforms into ice VI through multiple freezing–melting pathways at room temperature, hidden within the pressure region of ice VI. These multiple transition pathways occur via a metastable ice (more specifically, ice XXI with body-centred tetragonal structure ($I\bar{4}2d$)) discovered in this study and a metastable ice VII that exists within the pressure range of ice VI. We find that supercompressed water structurally evolves from high-density water to very-high-density water, causing multiple transition pathways. These findings provide an insight to find more metastable ice phases and their transition pathways at elevated temperatures.

Lee, Yun-Hee [Korea Research Institute of Standard↗
DOE OSTI2024

Microsphere LiMn 0.6 Fe 0.4 PO 4 /C cathode with unique rod-like secondary architecture for high energy lithium ion batteries

LiMn x Fe 1-x PO 4 /C is considered a promising next-generation cathode material with significant commercial potential, inheriting the safety of LiFePO 4 while offering higher energy densities. However, the extremely low conductivity and the Jahn-Teller effect induced by Mn 3+ limit its practical capacity and rate performance. Effective modifications can be achieved through particle nanonization and uniform carbon coating. Here, in this work, we synthesized microspherical LiMn 0.6 Fe 0.4 PO 4 /C cathode materials using a hydrothermal method combined with spray drying carbon coating. The cathode material exhibits a microsphere structure composed of aggregated nanorods with a uniform 3 nm carbon coating, showing good dispersibility, small specific surface area and high tap density. In-situ diffraction analysis showed that expanding the single-phase solid solution region during (de)lithiation can reduce the energy barrier for electron transport, improve the kinetics of the (dis)charge process, and enhance both cycling and rate performance. The initial capacity at 0.1C can reach 155 mAh/g, and the capacity remains at 133.5 mAh/g with a retention rate of 97.1 % after 300 cycles. The synergistic effect of particle nanonization and uniform carbon coating endows the LiMn x Fe 1-x PO 4 /C material with excellent electrochemical performance.

25 ENERGY STORAGE↗
DOE OSTI2025

Resolution requirements for numerical modeling of neutrino quantum kinetics

Neutrino quantum kinetics is a rapidly evolving field in computational astrophysics, with a primary focus on collective neutrino oscillations in core-collapse supernovae and postmerger phases of binary neutron star mergers. In recent years, there has been considerable debate concerning resolution dependence in numerical simulations. In this paper, we conduct a comprehensive resolution study in both angular- and spatial directions by using two independent schemes of quantum kinetic neutrino transport: finite volume and pseudospectral methods. We complement our discussion by linear stability analysis including inhomogeneous modes. Our result suggests that decreasing spatial resolutions underestimates the growth of flavor instability, and then leads to wrong asymptotic states of flavor conversions, which potentially has a critical impact on astrophysical consequences. We further delve into numerical results of low resolution simulations, that reveals the underlying mechanism responsible for numerical artifacts caused by insufficient resolutions. Finally, this study settles the debate on requirements of resolutions and serves as a guideline for numerical modeling of quantum kinetic neutrino transport.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Investigation of kinetic inductance and microwave loss in thin-film TaC x N 1−x superconducting resonators

The promise of tantalum for realizing superconducting quantum devices has generated interest in its compound films, particularly nitrides. Among these, cubic-phase tantalum carbonitride (TaC x N 1−x ) offers reduced susceptibility to oxidation and a high critical temperature, yet its microwave properties remain largely unexplored. Here, in this study, we investigate plasma-enhanced atomic layer deposition of cubic-phase TaC x N 1−x thin films for superconducting microwave circuits. Structural and transport measurements reveal nanocrystalline morphology with sub-10 nm grains and superconductivity in the dirty limit. Coplanar waveguide resonators exhibit moderately high kinetic inductance (16.8 pH/sq for 30 nm films) with potential for enhancement through dimensional scaling. The films also support high internal quality factors exceeding 10 5 at 50 mK in the single-photon regime, comparable to granular aluminum. Loss analysis identifies the two-level systems as the dominant limiting mechanism, with potential for further reduction through interface engineering. These results establish atomic layer deposited TaC x N 1−x as a promising material for scalable, low-loss, high-inductance superconducting circuits.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Resistive wall mode stability and resonant field amplification in MAST high beta plasma

The n=1 (n is the toroidal mode number) resistive wall mode (RWM) instability and the resonant field amplification (RFA) due to a stable RWM response are numerically investigated for a MAST high-pressure plasma scenario, utilizing the MARS-F (Liu Y.Q. et al 2000 Phys. Plasmas 7 3681) and MARS-K (Liu Y.Q. et al 2008 Phys. Plasmas 15 112503) codes, for the purpose of understanding the potential RWM behavior in MAST-U plasmas. Ideal kink stability analysis identifies a target plasma, with parameters similar to that of the reference equilibrium reconstructed from a MAST high-pressure discharge, that accesses the RWM regime. The unstable n=1 RWM for the target equilibrium is subject to strong damping by the plasma toroidal flow and/or the drift kinetic effects from thermal particles. As a result, the mode is found to be stable under the experimental flow conditions, or even without flow stabilization if drift kinetic stabilization is included. The stability prediction is robust against variation of the assumed resistive wall minor radius. Active magneto-hydrodynamic (MHD) spectroscopy modeling, using the magnetic coils designed for controlling the edge localized modes (ELMs) in MAST-U as the antenna, shows strong resonant field amplification due to a stable RWM response in the target plasma. Maximal amplification, of 6-7 times larger than the applied vacuum field within the plasma, is obtained assuming a coil phasing of 120-180 degrees between the upper and lower rows of the ELM control coils. The MHD-kinetic hybrid model predicts generally higher RFA than the fluid model, in particular near the low-field and high-field sides of the plasma boundary surface. Furthermore, the MARS-F/K modeling thus shows that the RWM in MAST-U is subject to strong flow and kinetic stabilization, but can nevertheless be detected by active MHD spectroscopy.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Computational simulation of asymmetric phase transformation in cracked Li₇La₃Zr₂O₁₂: Variant selection and chemo-mechanical implications

Coupling between microcracks and phase transformation in ion-conducting ceramics can jointly affect mechanical responses and ion transport. In this work we investigate the cubic-to-tetragonal phase transformation in Li₇La₃Zr₂O₁₂ in the presence of a microcrack under hydrostatic loading and quantify its implications for crack-tip stress concentration and Li-ion transport using phase-field and molecular dynamics simulations. The phase transformation exhibits a strong asymmetry between hydrostatic tension and compression. Under tension, the crack edge nucleates one tetragonal variant that amplifies the crack-tip stress intensity and promotes crack opening. Under compression, the crack tip nucleates a different tetragonal variant that enhances the stress-induced crack-closure tendency. Effective Li diffusivity analysis shows faster transport degradation under compression due to accelerated transformation kinetics, exposing a trade-off between mechanical stability and ionic conductivity. These results highlight the intertwined nature of cracking, phase transformation, and ionic transport in ion-conducting oxides and provide mechanistic insights into chemo-mechanical degradation of solid electrolytes.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Experimental Investigation of the Thermal Decomposition Pathways and Kinetics of TATB by Isotopic Substitution

Real-time measurements of the product gases arising from the thermal decomposition of triamino-trinitro benzene (TATB), its deuterated analogue, and plastically bonded TATB (LX-17) are presented in this study. Gas-phase decomposition products are identified by IR absorption spectroscopy. The frequency shifts in rovibrational spectra due to isotopic substitution and the change in rate of formation of decomposition products due to the kinetic-isotope-effect (KIE) help elucidate the decomposition pathways. The formation of H 2 O precedes other molecules (e.g., HCN, HNCO) during decomposition. After the concentrations of HCN and HNCO molecules reach a peak, their amounts gradually decrease. The concentrations of the other decomposition products (e.g., NH 3 and CO 2 ) rapidly rise after an induction period, which is attributed to the presence of autocatalytic reactions. The trends of chemical evolution are similar for all the samples, but their kinetic behaviors are different. This indicates the rates of consistent pathways are changed during thermal decomposition. Additionally, the kinetics of deuterated TATB decomposition is slower than that of unsubstituted TATB due to the KIE (kH/kD~1.41). The rate of LX-17 decomposition is slightly lower than unsubstituted TATB (k TATB /k LX-17 ~1.15). The KIE is more pronounced during the early stage of decomposition, which is attributed to the first steps of TATB decomposition involving water formation (i.e., H vs D transfer). The KIE slows down the formation of all gases, including those lacking hydrogen (e.g., CO 2 ). These results suggest the TATB thermal decomposition mechanism might involve a series of pathways rather than a set of independent and parallel reactions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗