The crystal structure of β-RDX—an elusive form of an explosive revealed

High-pressure structural studies and pressure-induced sensitisation of 3,4,5-trinitro-1H-pyrazole

Herein we report the first high-pressure study of the energetic material 3,4,5-trinitro-1H-pyrazole (3,4,5-TNP) using neutron powder diffraction and single-crystal X-ray diffraction. A new high-pressure phase, termed Form II, was first identified through a substantial change in the neutron powder diffraction patterns recorded over the range 4.6-5.3 GPa, and was characterised further by compression of a single crystal to 5.3 GPa in a diamond-anvil cell using X-ray diffraction. 3,4,5-TNP was found to be sensitive to initiation under pressure, as demonstrated by its unexpected and violent decomposition at elevated pressures in successive powder diffraction experiments. Initiation coincided with the sluggish phase transition from Form I to Form II. Using a vibrational up-pumping model, its increased sensitivity under pressure can be explained by pressure-induced mode hardening. These findings have potential implications for the safe handling of 3,4,5-TNP, on the basis that shock- or pressure-loading may lead to significantly increased sensitivity to initiation.

Ab Initio Crystal Structure Prediction of the Energetic Materials LLM-105, RDX, and HMX

Crystal structure prediction (CSP) is performed for the energetic materials (EMs) LLM-105 and α-RDX, as well as the α and β conformational polymorphs of 1,3,5,7-tetranitro-1,3,5,7-tetraazacyclooctane (HMX), using the genetic algorithm (GA) code, GAtor, and its associated random structure generator, Genarris. Genarris and GAtor successfully generate the experimental structures of all targets. GAtor's symmetric crossover scheme, where the space group symmetries of parent structures are treated as genes inherited by offspring, is found to be particularly effective. However, conducting several GA runs with different settings is still important for achieving diverse samplings of the potential energy surface. For LLM-105 and α-RDX, the experimental structure is ranked as the most stable, with all of the dispersion-inclusive density functional theory (DFT) methods used here. For HMX, the α form was persistently ranked as more stable than the β form, in contrast to experimental observations, even when correcting for vibrational contributions and thermal expansion. This may be attributed to insufficient accuracy of dispersion-inclusive DFT methods or to kinetic effects not considered here. In general, the ranking of some putative structures is found to be sensitive to the choice of the DFT functional and the dispersion method. For LLM-105, GAtor generates a putative structure with a layered packing motif, which is desirable thanks to its correlation with low sensitivity. Our results demonstrate that CSP is a useful tool for studying the ubiquitous polymorphism of EMs and shows promise of becoming an integral part of the EM development pipeline.

XDM-corrected hybrid DFT with numerical atomic orbitals predicts molecular crystal lattice energies with unprecedented accuracy

Molecular crystals are important for many applications, including energetic materials, organic semiconductors, and the development and commercialization of pharmaceuticals. The exchange-hole dipole moment (XDM) dispersion model has shown good performance in the calculation of relative and absolute lattice energies of molecular crystals, although it has traditionally been applied in combination with plane-wave/pseudopotential approaches. This has limited XDM to use with semilocal functional approximations, which suffer from delocalization error and poor quality conformational energies, and to systems with a few hundreds of atoms at most due to unfavorable scaling. In this work, we combine XDM with numerical atomic orbitals, which enable the efficient use of XDM-corrected hybrid functionals for molecular crystals. We test the new XDM-corrected functionals for their ability to predict the lattice energies of molecular crystals for the X23 set and 13 ice phases, the latter being a particularly stringent test. A composite approach using a XDM-corrected, 25% hybrid functional based on B86bPBE achieves a mean absolute error of 0.48 kcal mol^-1 per molecule for the X23 set and 0.19 kcal mol^-1 for the total lattice energies of the ice phases, compared to recent diffusion Monte-Carlo data. These results make the new XDM-corrected hybrids not only far more computationally efficient than previous XDM implementations, but also the most accurate density-functional methods for molecular crystal lattice energies to date.

Molecular Dynamics Study on the Reaction of RDX Molecule with Si Substrate

RDX is widely used in various explosion situations, and there are many studies on its detonation performance, safety, preparation, etc. Research on preparation of β-RDX is mainly conducted by experiments. In recent years, part of the research points to the use of substrate as a medium to produce β-RDX faster. Based on this guidance, our work aims to theoretically solve the physical and chemical processes that RDX may experience in the production process through numerical simulation. In this work, molecular dynamics simulation is set up for the interaction between RDX and a Si clean surface and a Si hydroxyl saturated surface separately, and a higher precision simulation is set up to verify the reliability of the results. NCI analysis is also used to guess the possible phase transition mechanism in the simulation results. In the simulation process, a 7 × 7 Si clean surface, a 3 × 3 Si clean surface, and a 7 × 7 Si-OH surface are set, and each surface adsorbs one α-RDX. The semiempirical Gfn1-xtb method is used for the 7 × 7 surface, and the DFT method is used for the 3 × 3 surface. The calculation results confirmed by high-precision results show that RDX molecules will react with the dangling bonds on the Si surface. Three conformations of RDX were found on the hydroxyl saturated surface of Si. The isosurface generated by the NCI method is used to analyze the reasons for the formation of these conformations.

Pressure and Polymorph Dependent Thermal Decomposition Mechanism of Molecular Materials: A Case of 1,3,5,-Trinitro-1,3,5,-triazine

1,3,5,-Trinitro-1,3,5,-triazine (RDX) serves as an important energetic material and is widely used as various solid propellants and explosives. Understanding the thermal decomposition behaviors of various polymorphs of RDX at high pressure and high temperature is significantly important for safe storage and handling. The present work reveals the early thermal decay mechanisms of two polymorphs (α- and ε-forms) of RDX at high pressure and high temperature by ReaxFF reactive molecular dynamic simulations and climbing image nudged elastic band (CI-NEB) static calculations. It is found that the thermal decomposition rate has positive and negative effects on the pressure for α- and ε-RDX, respectively. This difference originates from the difference of pressure effect on the intermolecular H transfer of the two polymorphs, as we confirm that the bimolecular H transfer rather than the NO₂ partition initiates the decay with a significantly lower energy barrier therein. This finding may be beneficial to understand the pressure and polymorph dependent effect on the decay of RDX and to develop a kinetic model for the combustion of solid RDX.

Low Temperature and High-Pressure Study of Bending L-Leucinium Hydrogen Maleate Crystals

The polymorphism of molecular crystals is a well-known phenomenon, resulting in modifications of physicochemical properties of solid phases. Low temperatures and high pressures are widely used to find phase transitions and quench new solid forms. In this study, L-Leucinium hydrogen maleate (LLHM), the first molecular crystal that preserves its anomalous plasticity at cryogenic temperatures, is studied at extreme conditions using Raman spectroscopy and optical microscopy. LLHM was cooled down to 11 K without any phase transition, while high pressure impact leads to perceptible changes in crystal structure in the interval of 0.0–1.35 GPa using pentane-isopentane media. Surprisingly, pressure transmitting media (PTM) play a significant role in the behavior of the LLHM system at extreme conditions—we did not find any phase change up to 3.05 GPa using paraffin as PTM. A phase transition of LLHM to amorphous form or solid–solid phase transition(s) that results in crystal fracture is reported at high pressures. LLHM stability at low temperatures suggests an alluring idea to prove LLHM preserves plasticity below 77 K.

Bottom-up coarse-grain modeling of plasticity and nanoscale shear bands in α-RDX

Computationally inexpensive particle-based coarse-grained (CG) models are essential for use in molecular dynamics (MD) simulations of mesoscopically slow cooperative phenomena, such as plastic deformations in solids. Molecular crystals possessing complex symmetry present enormous practical challenges for particle-based coarse-graining at molecularly resolved scales, when each molecule is in a single-site representation, and beyond. Presently, there is no published pairwise non-bonded single-site CG potential that is able to predict the space group and structure of a molecular crystal. In this paper, we present a successful coarse-graining at a molecular level from first principles of an energetic crystal, hexahydro-1,3,5-trinitro-s-triazine (RDX) in the alpha phase, using the force-matching-based multiscale coarse-graining (MSCG/FM) approach. The new MSCG/FM model, which implements an optimal pair decomposition of the crystal Helmholtz free energy potential in molecular center-of-mass coordinates, was obtained by force-matching atomistic MD simulations of liquid, amorphous, and crystalline states and in a wide range of pressures (up to 20 GPa). The MSCG/FM potentials for different pressures underwent top-down optimization to fine-tune the mechanical and thermodynamic properties, followed by consolidation into a transferable density-dependent model referred to as RDX-TC-DD (RDX True-Crystal Density-Dependent). The RDX-TC-DD model predicts accurately the crystal structure of α-RDX at room conditions and reproduces the atomistic reference system under isothermal (300 K) hydrostatic compression up to 20 GPa, in particular, the Pbca symmetry of α-RDX in the elastic regime. The RDX-TC-DD model was then used to simulate the plastic response of uniaxially ([100]) compressed α-RDX resulting in nanoscale shear banding, a key mechanism for plastic deformation and defect-free detonation initiation proposed for many molecular crystalline explosives. Additionally, a comparative analysis of the effect of core-softening of the RDX-TC-DD potential and the degree of molecular rigidity in the all-atom treatment suggests a stress-induced short-range softening of the effective intermolecular interaction as a fundamental cause of plastic instability in α-RDX. The reported RDX-TC-DD model and overall workflow to develop it open up possibilities to perform high quality simulation studies of molecular energetic materials under thermal and mechanical stimuli, including extreme conditions.

Strategies for Achieving Balance between Detonation Performance and Crystal Stability of High-Energy-Density Materials

Performance-stability contradiction of high-energy-density materials (HEDMs) is a long-standing puzzle in the field of chemistry and material science. Bridging the gap that exists between detonation performance of new HEDMs and their stability remains a formidable challenge. Achieving optimal balance between the two contradictory factors is of a significant demand for deep-well oil and gas drilling, space exploration, and other civil and defense applications. Herein, supercomputers and latest quantitative computational strategies were employed and high-throughput quantum calculations were conducted for 67 reported HEDMs. Based on statistical analysis of large amounts of physico-chemical data, in-crystal interspecies interactions were identified to be the one that provokes the performance-stability contradiction of HEDMs. To design new HEDMs with both good detonation performance and high stability, the proposed systematic and comprehensive strategies must be satisfied, which could promote the development of crystal engineering of HEDMs to an era of theory-guided rational design of materials.

Modification of crystalline energetic salts through polymorphic transition: enhanced crystal density and energy performance

We presented a detailed investigation of polymorphic transition of energetic salts and explored a new path for modifying crystalline energetic salts.

Crystal structure evolution of an energetic compound dihydroxylammonium 5,5′-bistetrazole-1,1′-diolate induced by solvents

Recently, energetic ionic salts have become a research hotspot due to their attractive properties, such as high density, high heat of formation, and environmental friendliness. Dihydroxylammonium 5,5′-bistetrazole-1,1′-diolate (TKX-50) is a typical nitrogen-rich energetic ionic salt, which has broad application prospects. However, the research on the stability and crystal structure evolution of TKX-50 in different solvent systems is insufficient. Herein, we investigated the crystal structure transformations and searched for new solid forms of TKX-50 under different conditions via a solvent induction method. The phase composition of all screened samples was analyzed by powder or single-crystal X-ray diffraction. Three new solid forms of [NH₂(CH₃)₂⁺][BTO⁻], [NH₂(CH₃CH₂)₂⁺]₂[BTO²⁻], [NHOH(CH₃CCH₃)⁺][BTO⁻] H₂O were obtained from DMAC, DEF and AC/MT, respectively. Furthermore, the energetic properties were evaluated through EXPLO5.

Three new solvent compounds obtained by dihydroxylammonium 5,5′-bistetrazole-1,1′-diolate recrystallized from three solvents.

Crystal structure evolution of the high energetic compound carbonic dihydrazidinium bis[3-(5-nitroimino-1,2,4-triazolate)] induced by solvents

Four new solid forms of CBNT including CBNT·2H₂O, H₂BNT·2DMSO, H₂BNT·2H₂O and [NH₂(CH₃)₂⁺]₂[BNT²⁻]·2H₂O were successfully obtained from water, DMSO, BL/H₂O and DMF/H₂O through a solvent induction method.

Trends of the macroscopic behaviors of energetic compounds: insights from first-principles calculations

Understanding the structure-property relationships of energetic compounds is challenging. Herein, by including the experimental data, we systematically evaluated the microscopic characteristics of a series of transition metal carbohydrazide perchlorate (TMCP) complexes (MnCP, FeCP, CoCP, NiCP, ZnCP, and CdCP) by first-principles calculations. The calculated properties, i.e., lattice enthalpy, bulk modulus, and electronic structures, were correlated with their thermal decomposition temperatures and impact sensitivities, which indicated that the stability and sensitivity of the TMCP complexes were greatly changed through coordination with different metal ions. The trend was that a large lattice enthalpy indicated a better thermal stability. Complexes with a high impact sensitivity tended to have a smaller bulk modulus and pseudo-gap. The ultra-high impact sensitivity of FeCP may have been related to the unstable spin state with respect to the volume change in the lattice. The calculated bond order and bond dissociation energy did not fully reflect the impact and friction sensitivities in this study. In addition, the combination of crystal properties and local bond information may better describe the sensitivity trend for the TMCP energetic compounds. This analysis can be applied to other energetic compounds and may provide clues for the synthesis and assessment of novel energetic compounds.

Charge Distributions of Nitro Groups Within Organic Explosive Crystals: Effects on Sensitivity and Modeling

The charge distribution of NO₂ groups within the crystalline polymorphs of energetic materials strongly affects their explosive properties. We use the recently introduced basis-space iterated stockholder atom partitioning of high-quality charge distributions to examine the approximations that can be made in modeling polymorphs and their physical properties, using 1,3,5-trinitroperhydro-1,3,5-triazine, trinitrotoluene, 1-3-5-trinitrobenzene, and hexanitrobenzene as exemplars. The NO₂ charge distribution is strongly affected by the neighboring atoms, the rest of the molecules, and also significantly by the NO₂ torsion angle within the possible variations found in observed crystal structures. Thus, the proposed correlations between the molecular electrostatic properties, such as trigger-bond potential or maxima in the electrostatic potential, and impact sensitivity will be affected by the changes in conformation that occur on crystallization. We establish the relationship between the NO₂ torsion angle and the likelihood of occurrence in observed crystal structures, the conformational energy, and the charge and dipole magnitude on each atom, and how this varies with the neighboring groups. We examine the effect of analytically rotating the atomic multipole moments to model changes in torsion angle and establish that this is a viable approach for crystal structures but is not accurate enough to model the relative lattice energies. This establishes the basis of transferability of the NO₂ charge distribution for realistic nonempirical model intermolecular potentials for simulating energetic materials.

Van der Waals Density Functional Theory vdW-DFq for Semihard Materials

There are a large number of materials with mild stiffness, which are not as soft as tissues and not as strong as metals. These semihard materials include energetic materials, molecular crystals, layered materials, and van der Waals crystals. The integrity and mechanical stability are mainly determined by the interactions between instantaneously induced dipoles, the so called London dispersion force or van der Waals force. It is challenging to accurately model the structural and mechanical properties of these semihard materials in the frame of density functional theory where the non-local correlation functionals are not well known. Here, we propose a van der Waals density functional named vdW-DFq to accurately model the density and geometry of semihard materials. Using β -cyclotetramethylene tetranitramine as a prototype, we adjust the enhancement factor of the exchange energy functional with generalized gradient approximations. We find this method to be simple and robust over a wide tuning range when calibrating the functional on-demand with experimental data. With a calibrated value q = 1.05 , the proposed vdW-DFq method shows good performance in predicting the geometries of 11 common energetic material molecular crystals and three typical layered van der Waals crystals. This success could be attributed to the similar electronic charge density gradients, suggesting a wide use in modeling semihard materials. This method could be useful in developing non-empirical density functional theories for semihard and soft materials.

Boosting the performance of energetic materials through thermally-induced conformational transition

We presented an effective strategy to improve the performance of energetic materials through thermally-induced conformational transition.

A unified model of impact sensitivity of metal azides

A comprehensive theoretical study of 18 metal azides is reported.

Selective preparation of elusive and alternative single component polymorphic solid forms through multi-component crystallisation routes

A transferable, simple, method for producing previously elusive and novel polymorphic forms of important active pharmaceutical ingredients (APIs; paracetamol (acetaminophen), piroxicam and piracetam) is demonstrated. Nitrogen heterocyclic co-molecules are employed to influence the self-assembly crystallisation process in a multi-component environment. Previously unknown solvates have also been synthesised by this method.

Polymorphic Phase Control of RDX-Based Explosives

The polymorphic phase of 1,3,5-trinitro-1,3,5-triazine (RDX) was examined as a function of mass loading, solvent, and sample deposition technique. When RDX was deposited at a high mass loading, the vibrational modes in the obtained Raman spectra were indicative of concomitant polymorphism as both the α-RDX and β-RDX phases were present. At low mass loadings, only β-RDX was observed regardless of solvent when using the drop cast crystallization method. However, α-RDX (the thermodynamically stable polymorphic phase observed with visible quantities of the explosive) was observed when RDX deposits were dry transferred. Observation of α-RDX was independent of the initial mass loading or the initial deposition solvent when using the dry transfer methodology. These data indicate that the use of the dry transfer preparation method can be used to successfully prepare RDX-based test articles with the α-RDX phase regardless of the solvent used to initially dissolve the RDX, the initial deposition technique, or the mass loading.

Effect of pressure gradient and new phases for 1,3,5-trinitrohexahydro-s-triazine (RDX) under high pressures

Herein, pressure-induced phase transitions of RDX up to 50 GPa were systematically studied under different compression conditions.

Exchange-Hole Dipole Dispersion Model for Accurate Energy Ranking in Molecular Crystal Structure Prediction II: Nonplanar Molecules

The crystal structure prediction (CSP) of a given compound from its molecular diagram is a fundamental challenge in computational chemistry with implications in relevant technological fields. A key component of CSP is the method to calculate the lattice energy of a crystal, which allows the ranking of candidate structures. This work is the second part of our investigation to assess the potential of the exchange-hole dipole moment (XDM) dispersion model for crystal structure prediction. In this article, we study the relatively large, nonplanar, mostly flexible molecules in the first five blind tests held by the Cambridge Crystallographic Data Centre. Four of the seven experimental structures are predicted as the energy minimum, and thermal effects are demonstrated to have a large impact on the ranking of at least another compound. As in the first part of this series, delocalization error affects the results for a single crystal (compound X), in this case by detrimentally overstabilizing the π-conjugated conformation of the monomer. Overall, B86bPBE-XDM correctly predicts 16 of the 21 compounds in the five blind tests, a result similar to the one obtained using the best CSP method available to date (dispersion-corrected PW91 by Neumann et al.). Perhaps more importantly, the systems for which B86bPBE-XDM fails to predict the experimental structure as the energy minimum are mostly the same as with Neumann's method, which suggests that similar difficulties (absence of vibrational free energy corrections, delocalization error,...) are not limited to B86bPBE-XDM but affect GGA-based DFT-methods in general. Our work confirms B86bPBE-XDM as an excellent option for crystal energy ranking in CSP and offers a guide to identify crystals (organic salts, conjugated flexible systems) where difficulties may appear.

Carbon dioxide binary crystals via the thermal decomposition of RDX at high pressure

A carbon dioxide and nitrous oxide solid solution has been captured in a diamond anvil cell following the thermal decomposition of 1,3,5-trinitroperhydro-1,3,5-triazine (RDX) at high temperatures and pressures. This is the first time a carbon dioxide binary solid has been observed at high pressure. This observation has stimulated low temperature crystallographic studies of this binary system using recently developed gas absorption apparatus and computational modelling.

Exchange-Hole Dipole Dispersion Model for Accurate Energy Ranking in Molecular Crystal Structure Prediction

Accurate energy ranking is a key facet to the problem of first-principles crystal-structure prediction (CSP) of molecular crystals. This work presents a systematic assessment of B86bPBE-XDM, a semilocal density functional combined with the exchange-hole dipole moment (XDM) dispersion model, for energy ranking using 14 compounds from the first five CSP blind tests. Specifically, the set of crystals studied comprises 11 rigid, planar compounds and 3 co-crystals. The experimental structure was correctly identified as the lowest in lattice energy for 12 of the 14 total crystals. One of the exceptions is 4-hydroxythiophene-2-carbonitrile, for which the experimental structure was correctly identified once a quasi-harmonic estimate of the vibrational free-energy contribution was included, evidencing the occasional importance of thermal corrections for accurate energy ranking. The other exception is an organic salt, where charge-transfer error (also called delocalization error) is expected to cause the base density functional to be unreliable. Provided the choice of base density functional is appropriate and an estimate of temperature effects is used, XDM-corrected density-functional theory is highly reliable for the energetic ranking of competing crystal structures.

Understanding metastable phase transformation during crystallization of RDX, HMX and CL-20: experimental and DFT studies

Multiphase growth during crystallization severely affects deliverable output of explosive materials. Appearance and incomplete transformation of metastable phases are a major source of polymorphic impurities. This article presents a methodical and molecular level understanding of the metastable phase transformation mechanism during crystallization of cyclic nitramine explosives, viz. RDX, HMX and CL-20. Instantaneous reverse precipitation yielded metastable γ-HMX and β-CL-20 which undergo solution mediated transformation to the respective thermodynamic forms, β-HMX and ε-CL-20, following 'Ostwald's rule of stages'. However, no metastable phase, anticipated as β-RDX, was evidenced during precipitation of RDX, which rather directly yielded the thermodynamically stable α-phase. The γ→β-HMX and β→ε-CL-20 transformations took 20 and 60 minutes respectively, whereas formation of α-RDX was instantaneous. Density functional calculations were employed to identify the possible transition state conformations and to obtain activation barriers for transformations at wB97XD/6-311++G(d,p)(IEFPCM)//B3LYP/6-311G(d,p) level of theory. The computed activation barriers and lattice energies responsible for transformation of RDX, HMX and CL-20 metastable phases to thermodynamic ones conspicuously supported the experimentally observed order of phase stability. This precise result facilitated an understanding of the occurrence of a relatively more sensitive and less dense β-CL-20 phase in TNT based melt-cast explosive compositions, a persistent and critical problem unanswered in the literature. The crystalline material recovered from such compositions revealed a mixture of β- and ε-CL-20. However, similar compositions of RDX and HMX never showed any metastable phase. The relatively long stability with the highest activation barrier is believed to restrict complete β→ε-CL-20 transformation during processing. Therefore a method is suggested to overcome this issue.

New phase of ammonium nitrate: A monoclinic distortion of AN-IV

A new phase of ammonium nitrate (AN) is found using first principles evolutionary crystal structure search. It is this polymorph that is associated with the phase transition to previously unidentified phase, which was detected in experiment at 17 GPa upon appearance of the two extra peaks in Raman spectrum. The new phase has a monoclinic unit cell in the P21/m space group symmetry (AN-P21/m) and is similar to the known phase IV of AN (AN-IV) except the ammonium molecules are oriented differently relative to the nitrate molecules. The calculated free energy of AN-P21/m is found to be lower than AN-IV at pressures above 10.83 GPa. The equation of state of both AN-P21/m and AN-IV phases (volume vs hydrostatic pressure at room temperature) has been obtained within the quasi-harmonic approximation. The calculated Raman spectrum of both AN-P21/m and AN-IV as a function of pressure is in a good agreement with experiment. The energetic competitiveness of AN-IV and AN-P21/m at ambient conditions suggests a possibility of the phase transition in a small pressure-temperature range near ambient pressure and temperature.

A new polymorph of metacetamol

The existence of a new polymorph of metacetamol together with its properties are reported for the first time.

DFT study of structural, electronic, and absorption properties of crystalline β-RDX under pressures

The structural, electronic, and absorption properties of crystalline hexahydro-1,3,5-trinitro-1,3,5-triazine (RDX) under hydrostatic pressure of 0–50 GPa were studied using density functional theory. As the pressure increases, the lattice constants and the unit cell volume decrease gradually. The compressibility of RDX crystal is anisotropic. The different changing trends of bond lengths and bond angles under pressures show that both the ring opening and the N–NO2 cleavage are possible to trigger the RDX decomposition. The interactions between electrons, especially for the valence electrons, are strengthened under the influence of pressure. The absorption spectra of RDX display more and stronger bands in the fundamental absorption region at high pressure.

Mechanisms and risk assessments on the N-nitration of N-acetylhexahydro-s-triazines: understanding the preparation of RDX (2)

Although the N-nitration by nitric acid is widely used to synthesize nitramines in biological, medical, and explosive industries, little is known about the microscopic behavior when the nitrated substrates are tertiary amines. Hexahydro-1,3,5-triacetyl-s-triazine (TRAT) nitrated into hexahydro-1,3,5-trinitro-s-triazine (RDX) was theoretically investigated at the MP2/cc-PVDZ level. An O-to-N transnitration mechanism was put forward for the N-nitration of N-acetyl tertiary amines, including the formation of diverse complexes R'N(COCH3)RNO2(+) and deacetylate. The electron transfer results in the complex formation, and the acetyl-to-nitro electrophilic displacement leads to deacetylate. Presumably, the carbonyl groups (C═O) in N-acetyl tertiary amines serve as the hinged joint in the electron transfer. Three successive N-nitrations transform TRAT into RDX; their electron transfers are strongly exothermic by -21.1, -19.5, and -15.4 kcal/mol relative to TRAT + 3NO2(+), repectively, and their electrophilic displacements possess low activation Gibbs free energies of 9.0, 6.8, and 7.5 kcal/mol relative to the σ-complexes 6, 11, and 14, respectively. The rate constants of the single electron transfer (SET) and the acetyl-to-nitro displacement were estimated roughly by Marcus and transition-state (TS) theories, respectively, indicating that they are both fast with the strong exothermicity. The available experimental phenomena were well interpreted by the computational results.

Vibrational spectra of an RDX film over an aluminum substrate from molecular dynamics simulations and density functional theory

We report calculated vibrational spectra in the range of 0-3,500 cm(-1) of RDX (hexahydro-1,3,5-trinitro-1,3,5-triazine) molecules adsorbed on a model aluminum surface. A molecular film was modeled using two approaches: (1) density functional theory (DFT) was used to optimize a single RDX molecule interacting with its periodic images, and (2) a group of nine molecules extracted from the crystal structure was deposited on the surface and interacted with its periodic images via molecular dynamics (MD) simulations. In both cases, the molecule was initialized in the AAA conformer geometry having the three nitro groups in axial positions, and kept that conformation in the DFT examination, but some molecules were found to change to the AAE conformer (two nitro groups in axial and one in equatorial position) in the MD analysis. The vibrational spectra obtained from both methods are similar to each other, except in the regions where collective RDX intermolecular interactions (captured by MD simulations) are important, and compare fairly well with experimental findings.

High pressure-high temperature decomposition of γ-cyclotrimethylene trinitramine

Decomposition of γ-cyclotrimethylene trinitramine (γ-RDX) under high pressure-high temperature conditions was examined to elucidate the reactive behavior of RDX crystals. Vibrational spectroscopy measurements were obtained for single crystals in a diamond anvil cell (DAC) at pressures from 6 to 12 GPa and temperatures up to 600 K. Global decomposition rates, activation energies, and activation volumes at several pressures and temperatures below the P-T locus for the γ-RDX decomposition were obtained. Similar to ε-RDX, but in contrast to α-RDX, we found that pressure decelerates the decomposition of γ-RDX. The decomposition deceleration with pressure in the γ-phase can be attributed to pressure-inhibiting bond homolysis step(s). The main decomposition species were identified as N(2)O, CO(2), and H(2)O, in accord with the species reported for the α-phase decomposition at high pressures. This work complements previous studies on RDX at HP-HT conditions and provides comprehensive results on the reactive behavior of γ-RDX; the γ-phase plays a key role in RDX decomposition at P-T conditions relevant to shock wave initiation.

Comparison of visible and near-infrared Raman cross-sections of explosives in solution and in the solid state

Raman cross-sections of explosives in solution and in the solid state have been measured using visible and near-infrared excitation via secondary calibration. These measurements are valuable for both fundamental scientific purposes and applications in the standoff detection of explosives. The explosive compounds RDX, HMX, TNT, 2,4-DNT, 2,6-DNT, and ammonium nitrate were measured using discrete excitation wavelengths ranging from 532 nm to 785 nm. A comparison of the spectral features and cross-sections between the solid state and solution was performed. Comparison is also made to cross-sections measured with deep ultraviolet excitation.