Atomistic-Scale Simulations of the Initial Chemical Events in the Thermal Initiation of Triacetonetriperoxide

Molecular Dynamics Simulation of Thermal Nonequilibrium and Chemical Reaction Processes in Hydrogen Combustion

Using reactive force field (ReaxFF) and molecular dynamics simulation, we investigate the combustion process of hydrogen-oxygen systems in initial thermal nonequilibrium states with different translational and rovibrational temperatures for oxygen. The system studied in this work contains 300 oxygen molecules and 700 hydrogen molecules with a density of 7 times the air density. For this system, the characteristic relaxation times of oxygen and hydrogen vibrational energies are 0.173 and 0.249 ns, respectively. 0.6% of hydrogen undergoes a chemical reaction with oxygen during the thermal nonequilibrium relaxation stage. For the distribution of translational energy and vibrational energy of oxygen in the thermal nonequilibrium state, the maximum mean error of the statistical distribution in the simulation and the Boltzmann distribution at temperature calculated from the average kinetic energy of molecules is about 2.25 × 10-5. At the same time, it was observed in the simulation that many-body interactions play a certain role in the combustion process. Furthermore, we compare the ignition time and temperature rise behavior of different combustion mechanisms and molecular dynamics simulations starting from the thermal equilibrium state. These results will provide meaningful references for the construction of thermal nonequilibrium combustion chemical reaction mechanisms.

Rapid and selective screening of organic peroxide explosives using acid-hydrolysis induced chemiluminescence

Organic peroxide explosives (OPEs) are unstable, non-military, contemporary security threats often found in improvised explosive devices. Chemiluminescence (CL) can be used to detect OPEs, via radical formation consisting of peroxide moieties (-O-O-) under acidic conditions. However, selectivity for specific OPEs is hampered by the ubiquitous background of H₂O₂. Herein, we report the differentiation of hexamethylene triperoxide diamine (HMTD), triacetone triperoxide (TATP), and methyl ethyl ketone peroxide (MEKP) by specific flow injection analysis-CL (FIA-CL) signal profiles, after H₂SO₄ treatment. The radical degradation pathway of each structure, and its corresponding FIA-CL profile, was explored using mass spectrometry to reveal the rapid loss of -O-O- from TATP and HMTD structures, while MEKP formed CL signal-sustaining oligomers, as opposed to the immediate attenuation of H₂O₂. The CL response for OPEs in an aqueous media, measured via the described FIA-CL method, enabled ultra-trace limits of detection down to 0.40 μM for MEKP, 0.43 μM for HMTD, and 0.40 μM for TATP (combined linear range 1-83 μM with 95% confidence limit, n = 12). Expanded uncertainties of measurement (UM) of MEKP = ±0.98, HMTD = ±1.03, and TATP = ±1.1 (UM included probabilities of false positive and false negative as well as standard deviations of % recoveries and limit of detections of OPEs). Direct aqueous sample introduction via FIA-CL thus offers the prospect of rapid and selective screening of OPEs in security-heightened settings (e.g., airports), averting false positives from more ubiquitous H₂O₂.

Ab initio neural network MD simulation of thermal decomposition of a high energy material CL-20/TNT

CL-20 (2,4,6,8,10,12-hexanitro-2,4,6,8,10,12-hexaazaisowurtzitane, also known as HNIW) is one of the most powerful energetic materials. However, its high sensitivity to environmental stimuli greatly reduces its safety and severely limits its application. In this work, ab initio based neural network potential (NNP) energy surfaces for both β-CL-20 and CL-20/TNT co-crystals were constructed. To accurately simulate the thermal decomposition processes of these two crystal systems, reactive molecular dynamics simulations based on the NNPs were performed. Many important intermediate species and their associated reaction paths during the decomposition had been identified in the simulations and the direct results on detonation temperatures of both systems were provided. The simulations also showed clearly that 2,4,6-trinitrotoluene (TNT) molecules in the co-crystal act as a buffer to slow down the chain reactions triggered by nitrogen dioxide and this effect is more significant at lower temperatures. Specifically, the addition of TNT molecules in the CL-20/TNT co-crystal introduces intermolecular hydrogen bonds between CL-20 and TNT molecules in the system, thereby increasing the thermal stability of the co-crystal. The current reactive molecular dynamics simulation is performed based on the NNP which helps in accelerating the speed of ab initio molecular dynamics (AIMD) simulation by more than 3 orders of magnitude while preserving the accuracy of density functional theory (DFT) calculations. This enabled us to perform longer-time simulations at more realistic temperatures that traditional AIMD methods cannot achieve. With the advantage of the NNP in its powerful fitting ability and transferability, the NNP-based MD simulation can be widely applied to energetic material systems.

Study of the thermal decomposition mechanism of FOX-7 by molecular dynamics simulation and online photoionization mass spectrometry

The thermal decomposition mechanism of energetic materials is important for analyzing the combustion mechanisms of propellants and evaluating the safety of propellants during transport and storage. 1,1-Diamino-2,2-dinitroethylene (FOX-7) is an important insensitive energetic material that can be used as an oxidizer in propellants. However, the initial decomposition mechanism of FOX-7 is not clear to date. The ReaxFF molecular dynamics method is widely used in the investigation of the thermal decomposition mechanisms of energetic materials. Meanwhile, the combination of thermogravimetry with online photoionization time-of-flight mass spectrometry (TG-PI-TOF-MS) and online single-photon ionization time-of-flight mass spectrometry (SPI-TOF-MS) can reveal the decomposition products, which may be integrated with the results of the simulation. In this study, the primary thermal decomposition mechanism of 1,1-diamino-2,2-dinitroethylene (FOX-7) was studied by the ReaxFF molecular dynamics simulations and online photoionization mass spectrometry. The results of the molecular dynamics simulations showed that the primary decomposition step of FOX-7 is C–NO₂ cleavage; after this, C Created by potrace 1.16, written by Peter Selinger 2001-2019 ]]> O formation occurs via a three-membered ring transition state, followed by NO elimination. The remaining structure loses NH₂ and H, resulting in the formation of the NHCCO structure, which finally breaks down into HNC and CO. NH₂ reacts with an H atom to produce NH₃. A reversible intramolecular hydrogen transfer was also observed at 2500 K; however, it failed to dominate the decomposition reaction. During the decomposition of FOX-7, the major products are N₂, NH₃, CO₂, and H₂N₂ and the minor products are H₂O, HN₂, and H₂. The TG-PI-TOF-MS spectrum shows three signals, i.e., m/z = 18, 28, and 30, which can be assigned to H₂O, CO, and NO, respectively. Moreover, four signals at m/z = 72.72, 55.81, 45.79, and 29.88 corresponding to the products (NH₂)₂CCO, (NH₂)CCO, NO₂, and NO have been obtained in the SPI-TOF-MS spectrum. The experimental data obtained via online photoionization mass spectrometry further validated the results of the molecular dynamics simulations.

In this work, the primary thermal decomposition mechanism of 1,1-diamino-2,2-dinitroethylene (FOX-7) was studied by ReaxFF molecular dynamics simulations and online photoionization mass spectrometry.

Comparison study of carbon clusters formation during thermal decomposition of 1,3,5-triamino-2,4,6-trinitrobenzene and benzotrifuroxan: a ReaxFF based sequential molecular dynamics simulation

Carbon rich clusters are usually found after the detonation of explosives, which greatly hinder their further decomposition into small molecules. A comparison study of thermal decomposition and clusters formation between 1,3,5-triamino-2,4,6-trinitrobenzene (TATB) and benzotrifuroxan (BTF) crystals was conducted to uncover the mechanisms behind their distinct differences in sensitivity and reaction violence, which has not been investigated in detail. The simulations of heating at 3500 K, then expansion and cooling were conducted through reactive molecular dynamics using the ReaxFF-lg force field. As a result, the initial low decay rate indicates that TATB is more stable than BTF under high temperatures, while once ignited it decays faster than BTF. Nevertheless, BTF decomposes more completely with a higher potential energy release, a greater amount of final products, and higher reaction frequencies, and shows higher reaction violence than TATB. More and heavier clusters occur in TATB crystals compared with those in BTF. Large clusters form during the heating process and then partly dissociate during expansion and cooling. A faster cooling rate facilitates larger clusters formation. Graphitic geometries as well as carbon rings and carbon chains are common in the stable clusters. Besides, further simulations show that a lower heating temperature facilitates larger clusters formation both in TATB and BTF. Our results are expected to deepen the insight into the mechanisms of carbon clusters formation and the different performances of TATB and BTF.

ReacNetGenerator: an automatic reaction network generator for reactive molecular dynamics simulations

Reactive molecular dynamics (MD) simulation makes it possible to study the reaction mechanism of complex reaction systems at the atomic level. However, the analysis of MD trajectories which contain thousands of species and reaction pathways has become a major obstacle to the application of reactive MD simulation in large-scale systems. Here, we report the development and application of the Reaction Network Generator (ReacNetGenerator) method. It can automatically extract the reaction network from the reaction trajectory without any predefined reaction coordinates and elementary reaction steps. Molecular species can be automatically identified from the cartesian coordinates of atoms and the hidden Markov model is used to filter the trajectory noises which makes the analysis process easier and more accurate. The ReacNetGenerator has been successfully used to analyze the reactive MD trajectories of the combustion of methane and 4-component surrogate fuel for rocket propellant 3 (RP-3), and it has great advantages in terms of efficiency and accuracy compared to traditional manual analysis.

Transforming the Accuracy and Numerical Stability of ReaxFF Reactive Force Fields

Molecular dynamics (MD) simulations provide an important link between theories and experiments. While ab initio methods can be prohibitively costly, the ReaxFF force field has facilitated in silico studies of chemical reactivity in complex, condensed-phase systems. However, the relatively poor energy conservation in ReaxFF MD has either limited the applicability to short time scales, in cases where energy propagation is important, or has required a continuous coupling of the system to a heat bath. In this study, we reveal the root cause of the unsatisfactory energy conservation, and offer a straightforward solution. The new scheme results in orders of magnitude improvement in energy conservation, numerical stability, and accuracy of ReaxFF force fields, compared to the previous state-of-the-art, at no additional cost. We anticipate that these improvements will open new avenues of research for more accurate reactive simulations in complex systems on long time scales.

A reactive force field molecular dynamics study of molecular nitrogen and water mixtures under high temperature and high pressure

Establishment of an equation of state (EOS) for detonation products is challenging through molecular modeling approaches. Molecular dynamics (MD) simulations with the reactive force fields (ReaxFF) were carried out in this work for molecular N₂ and H₂O and their mixtures, which are key components of detonation products, with high density and under high temperature, and compared with the experiment-based Becker-Kistiakowsky-Wilson (BKW) and the Lennard-Jones potential-based Tsien-Zhao (TZ) equations. Phase separations were noted for the mixtures, enhanced for denser systems under low temperature, and detracted for sparse systems under high temperature. ReaxFF-based MD simulations reproduce the BKW predictions for the systems with a wide range of densities, while the TZ predictions deviate at relatively low density. The deviations among the MD, BKW, and TZ results reduce with increasing temperature and decreasing density. Our calculations revealed the significant influence of phase separation on the predicted pressures, and verified that the MD/ReaxFF approach is potentially a feasible way to establish the EOS for detonation products.

ReaxFF MD Simulations of Peptide-Grafted Gold Nanoparticles

Functionalized gold nanoparticles have critical applications in biodetection with surface-enhanced Raman spectrum and drug delivery. In this study, reactive force field molecular dynamics simulations were performed to study gold nanoparticles, which are modified with different short-chain peptides consisting of amino acid residues of cysteine and glycine in different grafting densities in the aqueous environment. Our study showed slight facet-dependent peptide adsorption on a gold nanoparticle with the 3 nm core diameter. Peptide chains prefer to adsorb on the Au(111) facet compared to those on other facets of Au(100) and Au(110). In addition to the stable thiol interaction with gold nanoparticle surfaces, polarizable oxygen and nitrogen atoms show strong interactions with the gold surface and polarize the gold nanoparticle surfaces with an overall positive charge. Charges of gold atoms vary according to their contacts with peptide atoms and lattice positions. However, at the outmost peptide layer, the whole functionalized Au nanoparticles exhibit overall negative electrostatic potential due to the grafted peptides. Moreover, simulations show that thiol groups can be deprotonated and subsequently protons can be transferred to water molecules and carboxyl groups.

Structure-Property of Lithium-Sulfur Nanoparticles via Molecular Dynamics Simulation

Lithium-sulfur (Li-S) batteries offer higher energy densities than most reported lithium-ion batteries. However, our understanding of Li-S battery is still largely unknown at the level of the nanoscale. The structural properties of Li-S materials were investigated via molecular dynamics (MD) simulations using the ReaxFF force field. Several Li-S nanoparticles with different Li/S composition ratios (2:1 and 2:8) and various structures are studied. Our MD simulations show that among the four structures we constructed for Li₂S₈ nanoparticles, the core-shell structure is the most thermodynamically stable one during the charging (delithiation) process. In contrast to bulk crystal Li₂S, we find the presence of mixed lithium sulfide and polysulfide species are common features for these Li-S (Li₂S, Li₂S₈) nanoparticles. The complex distribution of these sulfide and polysulfide speciation are dictated by both stoichiometry and local atomic structures in the nanoparticle. These findings will provide insight into further development of functionalized lithium-sulfur cathodes.

Reactive Molecular Dynamics Simulations of the Thermal Decomposition Mechanism of 1,3,3-Trinitroazetidine

1,3,3-Trinitroazetidine (TNAZ) has a molecular formula of C₃ H₄ N₄ O₆ and the characteristics of low melting point, low impact sensitivity and good thermal stability. It is suitable for melt casting and pressed charges, and it has broad prospects for applications in low-sensitivity ammunition. In this study, the thermal decomposition of TNAZ crystals at high temperature was calculated by molecular dynamics simulation with the ReaxFF/lg reactive force field. The change in the potential energy of TNAZ, the formation of small-molecule products and clusters, and the initial reaction path of TNAZ were analysed. The kinetic parameters of different reaction stages in TNAZ thermal decomposition were obtained. The primary thermal decomposition reaction of TNAZ was found to be as follows: N-NO₂ and C-NO₂ bonds broke; a H atom on the quaternary ring was transferred to the nitro group; and the C-HNO₂ and N-HNO₂ bonds broke. The main decomposition products of TNAZ were thus NO₂ , NO, N₂ , H₂ O, CO₂ and HNO₂ , as well as macromolecular clusters. The size of the cluster structure was related to the reaction temperature, and the higher the temperature was, the smaller the cluster size was.

First-principles-based reaction kinetics from reactive molecular dynamics simulations: Application to hydrogen peroxide decomposition

This paper presents our vision of how to use in silico approaches to extract the reaction mechanisms and kinetic parameters for complex condensed-phase chemical processes that underlie important technologies ranging from combustion to chemical vapor deposition. The goal is to provide an analytic description of the detailed evolution of a complex chemical system from reactants through various intermediates to products, so that one could optimize the efficiency of the reactive processes to produce the desired products and avoid unwanted side products. We could start with quantum mechanics (QM) to ensure an accurate description; however, to obtain useful kinetics we need to average over ∼10-nm spatial scales for ∼1 ns, which is prohibitively impractical with QM. Instead, we use the reactive force field (ReaxFF) trained to fit QM to carry out the reactive molecular dynamics (RMD). We focus here on showing that it is practical to extract from such RMD the reaction mechanisms and kinetics information needed to describe the reactions analytically. This analytic description can then be used to incorporate the correct reaction chemistry from the QM/ReaxFF atomistic description into larger-scale simulations of ∼10 nm to micrometers to millimeters to meters using analytic approaches of computational fluid dynamics and/or continuum chemical dynamics. In the paper we lay out the strategy to extract the mechanisms and rate parameters automatically without the necessity of knowing any details of the chemistry. We consider this to be a proof of concept. We refer to the process as RMD2Kin (reactive molecular dynamics to kinetics) for the general approach and as ReaxMD2Kin (ReaxFF molecular dynamics to kinetics) for QM-ReaxFF-based reaction kinetics.

Application of Q2MM to predictions in stereoselective synthesis

Quantum-Guided Molecular Mechanics (Q2MM) can be used to derive transition state force fields (TSFFs) that allow the fast and accurate predictions of stereoselectivity for a wide range of catalytic enantioselective reactions. The basic ideas behind the derivation of TSFFs using Q2MM are discussed and the steps involved in obtaining a TSFF using the Q2MM code, publically available at github.com/q2mm, are shown. The applicability for a range of reactions, including several non-standard applications of Q2MM, is demonstrated. Future developments of the method are also discussed.

Predicted detonation properties at the Chapman-Jouguet state for proposed energetic materials (MTO and MTO3N) from combined ReaxFF and quantum mechanics reactive dynamics

The development of new energetic materials (EMs) with improved detonation performance but low sensitivity and environmental impact is of considerable importance for applications in civilian and military fields. Often new designs are difficult to synthesize so predictions of performance in advance is most valuable. Examples include MTO (2,4,6-triamino-1,3,5-triazine-1,3,5-trioxide) and MTO3N (2,4,6-trinitro-1,3,5-triazine-1,3,5-trioxide) suggested by Klapötke as candidate EMs but not yet successfully synthesized. We propose and apply to these materials a new approach, RxMD(cQM), in which ReaxFF Reactive Molecular Dynamics (RxMD) is first used to predict the reaction products and thermochemical properties at the Chapman Jouguet (CJ) state for which the system is fully reacted and at chemical equilibrium. Quantum mechanics dynamics (QMD) is then applied to refine the pressure of the ReaxFF predicted CJ state to predict a more accurate final CJ point, leading to a very practical calculation that includes accurate long range vdW interactions needed for accurate pressure. For MTO, this RxMD(cQM) method predicts a detonation pressure of P_CJ = 40.5 GPa and a detonation velocity of D_CJ = 8.8 km s^-1, while for MTO3N it predicts P_CJ = 39.9 GPa and D_CJ = 8.4 km s^-1, making them comparable to HMX (P_CJ = 39.5 GPa, D_CJ = 9.1 km s^-1) and worth synthesizing. This first-principles-based RxMD(cQM) methodology provides an excellent compromise between computational cost and accuracy including the formation of clusters that burn too slowly, providing a practical mean of assessing detonation performances for novel candidate EMs. This RxMD(cQM) method that links first principles atomistic molecular dynamics simulations with macroscopic properties to promote in silico design of new EMs should also be of general applicability to materials synthesis and processing.

Atomic-Level Simulation Study of n-Hexane Pyrolysis on Silicon Carbide Surfaces

Ethylene production plays a key role in the petrochemical industry. The severe operation conditions of ethylene thermal cracking, such as high-temperature and coke-formation, pose challenges for the development of new corrosion-resistant and coking-resistant materials for ethylene reactor radiant coils tubes (RCTs). We investigated the performance of ceramic materials such as silicon carbide (SiC) in severe pyrolysis conditions by using reactive force field molecular dynamics (ReaxFF MD) simulation method. Our results indicate that β-SiC surface remains fully stable at 1500 K, whereas increased temperature results in melted interface. At 2500 K, fully grown cross-linked-graphene-like polycyclic aromatic hydrocarbon coking structure on SiC surfaces was observed. Such coking was particularly severe in the carbon-side of the surface slab. The coking structures were mainly derived from surface atoms at the initial 3.0 ns, as a result of the loss of interfacial hydroxyl layer and further hydrothermal corrosion. The SiC substrate surface enhances the ethylene cracking rate and also leads to different intermediate-state compounds. Our fundamental research will have significant and broad impact on both petrochemical industry and academic research in materials science, petrochemistry, and combustion chemistry.

Study of thermal conductivity of ice clusters after impact deposition on the silica surfaces using the ReaxFF reactive force field

During aircraft or spacecraft missions, ice accumulates on different parts of their surface elements. An important parameter affecting the ability to remove this ice from the surface is the heat transfer characteristics of the accumulated ice. The ice heat transfer is related to the process of ice formation and its density and internal structure. In this study we investigate the effects of the ice and silica structure and the ice cluster attachment mechanism to the silica surface on the thermal conductivity (TC) of the attached ice cluster using the ReaxFF reactive force field. The purpose of this study is to investigate the thermal transport in amorphous and crystalline ice after high-velocity deposition on the silica surfaces. A dual thermostat method has been applied for the calculation of TC values. The validity of this method has been verified by comparing the calculated values of TC for crystal and amorphous ice with available experimental values. Our calculations show that the TC values of both crystal and amorphous ice drop after deposition on the silica surfaces. This decrease in the TC is more significant for the ice deposition on suboxide silica surfaces. Furthermore, crystal ice shows higher TC values than amorphous ice after accumulation. However, when crystal ice impacts on the silica surface at 1 km s(-1) impact speed, the crystalline shape of the ice cluster is lost to a considerable level and the TC values obtained for the ice clusters in such cases are closer to amorphous ice TC values. We observed a decrease in the TC values when ionic species are added inside the ice clusters.

Differentiation among peroxide explosives with an optoelectronic nose

Forensic identification of batches of homemade explosives (HME) poses a difficult analytical challenge. Differentiation among peroxide explosives is reported herein using a colorimetric sensor array and handheld scanner with a field-appropriate sampling protocol. Clear discrimination was demonstrated among twelve peroxide samples prepared from different reagents, with a classification accuracy >98%.

Anisotropic shock sensitivity in a single crystal δ-cyclotetramethylene tetranitramine: a reactive molecular dynamics study

The anisotropic shock sensitivity in a single crystal δ-cyclotetramethylene tetranitramine (δ-HMX) was investigated using the compress-shear reactive dynamics (CS-RD) computational protocol. Significant anisotropies in the thermo-mechanical and chemical responses were found by measuring the shear stress, energy, temperature, and chemical reactions during the dynamical process for the shock directions perpendicular to the (100), (010), (001), (110), (101), (011), and (111) planes. We predict that δ-HMX is sensitive for the shocks perpendicular to the (111), (011), (110), and (101) planes, which is intermediate to the (100) and (010) plane and is insensitive to the (001) plane. The internal energy accumulated within the duration of the surmounting shear stress barrier is a useful criterion to distinguish the sensitive directions from the less sensitive ones. The molecular origin of the anisotropic sensitivity is suggested to be the intermolecular steric arrangements across a slip plane induced by shock compression. The shear deformation induced by the shock along the sensitive direction encounters strong intermolecular contacts and has small intermolecular free space for geometry relaxation when the molecules collide, leading to high shear stress barriers and energy accumulation, which benefits the temperature increase and initial chemical bond breaking that trigger further reactions.

A reactive molecular dynamics study on the anisotropic sensitivity in single crystal α-cyclotetramethylene tetranitramine

Anisotropic sensitivity is related to the different intermolecular steric arrangements across the slip plane induced by shock compression along various orientations.

Reactive molecular dynamics simulation on the disintegration of Kapton, POSS polyimide, amorphous silica, and teflon during atomic oxygen impact using the ReaxFF reactive force-field method

Atomic oxygen (AO) is the most abundant element in the low Earth orbit (LEO). It is the result of the dissociation of molecular oxygen by ultraviolet radiation from the sun. In the LEO, it collides with the materials used on spacecraft surfaces and causes degradation of these materials. The degradation of the materials on the surface of spacecrafts at LEO has been a significant problem for a long time. Kapton polyimide, polyhedral oligomeric silsesquioxane (POSS), silica, and Teflon are the materials extensively used in spacecraft industry, and like many other materials used in spacecraft industry, AO collision degradation is an important issue in their applications on spacecrafts. To investigate the surface chemistry of these materials in exposure to space AO, a computational chemical evaluation of the Kapton polyimide, POSS, amorphous silica, and Teflon was performed in separate simulations under similar conditions. For performing these simulations, the ReaxFF reactive force-field program was used, which provides the computational speed required to perform molecular dynamics (MD) simulations on system sizes sufficiently large to describe the full chemistry of the reactions. Using these simulations, the effects of AO impact on different materials and the role of impact energies, the content of material, and temperature of material on the behavior of the materials are studied. The ReaxFF results indicate that Kapton is less resistant than Teflon toward AO damage. These results are in good agreement with experiment. These simulations indicate that the amorphous silica shows the highest stability among these materials before the start of the highly exothermic silicon oxidation. We have verified that adding silicon to the bulk of the Kapton structure enhances the stability of the Kapton against AO impact. Our canonical MD simulations demonstrate that an increase in the heat transfer in materials during AO impact can provide a considerable decrease in the disintegration of the material. This effect is especially relevant in silica AO collision. Considerable experimental efforts have been undertaken to minimize such AO-based degradations. As our simulations demonstrate, ReaxFF can provide a cost-effective screening tool for future material optimization.

Development of a ReaxFF reactive force field for ammonium nitrate and application to shock compression and thermal decomposition

We have developed a new ReaxFF reactive force field parametrization for ammonium nitrate. Starting with an existing nitramine/TATB ReaxFF parametrization, we optimized it to reproduce electronic structure calculations for dissociation barriers, heats of formation, and crystal structure properties of ammonium nitrate phases. We have used it to predict the isothermal pressure-volume curve and the unreacted principal Hugoniot states. The predicted isothermal pressure-volume curve for phase IV solid ammonium nitrate agreed with electronic structure calculations and experimental data within 10% error for the considered range of compression. The predicted unreacted principal Hugoniot states were approximately 17% stiffer than experimental measurements. We then simulated thermal decomposition during heating to 2500 K. Thermal decomposition pathways agreed with experimental findings.

The identification of chlorinated acetones in analyses of aged triacetone triperoxide (TATP)

The organic peroxide explosive triacetone triperoxide (TATP) is regularly encountered by law enforcement agents in various stages of its production, storage, or usage. In a previous study, it has been demonstrated that isolated, rigorously purified, TATP may degrade to form a series of chlorinated acetones when directly treated with excess concentrated hydrochloric acid. The current study extends this work to examine whether this phenomenon may be measured during the more feasible scenario of aging of rudimentarily purified TATP contaminated with trace reaction mixture. It was demonstrated that solid-phase microextraction gas chromatography/mass spectrometry analyses of aged TATP that was synthesized utilizing hydrochloric acid catalyst may identify the presence of the degradation products chloroacetone and 1,1-dichloroacetone. Upon aging of TATP synthesized utilizing either sulfuric or nitric acid catalyst, no acid specific degradation products could be identified. These findings may be exploited by forensic chemists in the analyses of TATP samples.

Molecular dynamics simulations of weak detonations

Detonation of a three-dimensional reactive nonisotropic molecular crystal is modeled using molecular dynamics simulations. The detonation process is initiated by an impulse, followed by the creation of a stable fast reactive shock wave. The terminal shock velocity is independent of the initiation conditions. Further analysis shows supersonic propagation decoupled from the dynamics of the decomposed material left behind the shock front. The dependence of the shock velocity on crystal nonlinear compressibility resembles solitary behavior. These properties categorize the phenomena as a weak detonation. The dependence of the detonation wave on microscopic potential parameters was investigated. An increase in detonation velocity with the reaction exothermicity reaching a saturation value is observed. In all other respects the model crystal exhibits typical properties of a molecular crystal.