Zbigniew A. Dreger

Phase Diagram of Hexahydro-1,3,5-trinitro-1,3,5-triazine Crystals at High Pressures and Temperatures

Raman spectroscopy and optical imaging were used to determine the phase boundaries between various hexahydro-1,3,5-trinitro-1,3,5-triazine (RDX) polymorphs. Experiments were performed on single crystals at pressures up to 8.0 GPa and temperatures ranging from room temperature to 550 K. Several distinct pressure regions were found in the RDX response at elevated temperatures: (i) melting of alpha-RDX followed by decomposition, below 2.0 GPa, (ii) decomposition of alpha-RDX between 2.0 and 2.8 GPa, (iii) irreversible transformation of alpha- and gamma-RDX to epsilon-RDX between 2.8 and 6.0 GPa, and (iv) decomposition of gamma-RDX above 6.0 GPa. A triple point between the alpha-, gamma-, and epsilon-RDX was found at 3.7 GPa and 466 K. The alpha-gamma phase transition was confirmed to occur at the same pressure, approximately 3.7 GPa, regardless of temperature, in the range of 295-460 K. Furthermore, it was determined that epsilon-RDX (i) has limited chemical stability under the pressure and temperature conditions where it is produced and (ii) decomposes according to the autocatalytic rate law. The findings reported here have provided new information about the response of RDX crystals at high pressures and temperatures.

Shock Wave Induced Decomposition of RDX: Time-Resolved Spectroscopy

Time-resolved optical spectroscopy was used to examine chemical decomposition of RDX crystals shocked along the [111] orientation to peak stresses between 7 and 20 GPa. Shock-induced emission, produced by decomposition intermediates, was observed over a broad spectral range from 350 to 850 nm. A threshold in the emission response of RDX was found at about 10 GPa peak stress. Below this threshold, the emission spectrum remained unchanged during shock compression. Above 10 GPa, the emission spectrum changed with a long wavelength component dominating the spectrum. The long wavelength emission is attributed to the formation of NO2 radicals. Above the 10 GPa threshold, the spectrally integrated intensity increased significantly, suggesting the acceleration of chemical decomposition. This acceleration is attributed to bimolecular reactions between unreacted RDX and free radicals. These results provide a significant experimental foundation for further development of a decomposition mechanism for shocked RDX (following paper in this issue).

Raman Spectroscopy of High-Pressure-High-Temperature Polymorph of Hexahydro-1,3,5-trinitro-1,3,5-triazine (ε-RDX)

Raman spectroscopy was used to determine the vibrational structure and the stability of the high-pressure-high-temperature (HP-HT) polymorph of RDX after it had been quenched to room temperature. Although this polymorph has limited chemical stability under high pressure and temperature, we show that it is chemically and structurally stable from 0.6 GPa to at least 20 GPa at room temperature. Below 0.6 GPa, it readily converts to the alpha-polymorph. Pressure dependence of the vibrational structure of the HP-HT polymorph was measured and compared with the vibrational structures of other known RDX polymorphs: alpha, beta, and gamma. In contrast with previous suggestions, our data indicate that the HP-HT polymorph can have a different structure than the beta-polymorph. This finding supports the recent suggestion that the HP-HT polymorph should be given a separate designation, epsilon-RDX. Furthermore, symmetry correlation analyses of Raman spectra indicate that the HP-HT polymorph (epsilon-RDX) may assume the space group isomorphous with the C(2v)[C(1)(4)] point group and with molecules adopting the pseudo-AAA conformation.

Second-order elastic constants of pentaerythritol tetranitrate and cyclotrimethylene trinitramine using impulsive stimulated thermal scattering

Impulsive stimulated thermal scattering (ISTS) was used to determine the complete set of second-order elastic constants for pentaerythritol tetranitrate (PETN) and cyclotrimethylene trinitramine (RDX) single crystals. Despite the weak scattering efficiency of these materials, excellent signal quality was obtained by using an optical heterodyne detection approach. The elastic constants for PETN agree well with previous values obtained from ultrasonic velocity measurements. The elastic constants for RDX are consistent with previous values obtained from ultrasonic velocity measurements and from resonant ultrasound spectroscopy, but show significant differences with values obtained from Brillouin scattering data. The present results demonstrate that the ISTS method, with optical heterodyne detection, provides a useful and accurate approach for determining the elastic constants of energetic crystals.

Time-Resolved Spectroscopic Measurements of Shock-Wave Induced Decomposition in Cyclotrimethylene Trinitramine (RDX) Crystals: Anisotropic Response

Plate impact experiments on the (210), (100), and (111) planes were performed to examine the role of crystalline anisotropy on the shock-induced decomposition of cyclotrimethylenetrinitramine (RDX) crystals. Time-resolved emission spectroscopy was used to probe the decomposition of single crystals shocked to peak stresses ranging between 7 and 20 GPa. Emission produced by decomposition intermediates was analyzed in terms of induction time to emission, emission intensity, and the emission spectra shapes as a function of stress and time. Utilizing these features, we found that the shock-induced decomposition of RDX crystals exhibits considerable anisotropy. Crystals shocked on the (210) and (100) planes were more sensitive to decomposition than crystals shocked on the (111) plane. The possible sources of the observed anisotropy are discussed with regard to the inelastic deformation mechanisms of shocked RDX. Our results suggest that, despite the anisotropy observed for shock initiation, decomposition pathways for all three orientations are similar.

Shock wave induced decomposition of RDX: quantum chemistry calculations.

Quantum chemical calculations on single molecules were performed to provide insight into the decomposition mechanism of shocked RDX. These calculations complement time-resolved spectroscopy measurements on shock wave compressed RDX crystals (previous paper, this issue). It is proposed that unimolecular decomposition is the primary pathway for RDX decomposition in its early stages and at stresses lower than approximately 10 GPa. This decomposition leads to the generation of broadband emission from 350 to 850 nm. Chemiluminescence from (2)B1 and (2)B2 excited states of NO2 radicals is associated with a major portion of the experimentally observed emission spectrum (>400 nm). The remaining portion (<400 nm) of the emission spectrum primarily results from excited HONO intermediates. It is proposed that for stresses higher than 10 GPa, bimolecular reactions between radical decomposition products and unreacted RDX molecules become the dominant pathway. This radical assisted homolysis pathway is cyclic and leads to the acceleration of decomposition, with increased production of low energy NO2 radicals. These radicals produce emission that is stronger in the long wavelength portion of the spectrum. Finally, a comprehensive chemical decomposition mechanism is put forward that is consistent with the experimental observations of shock-induced emission in RDX crystals.

Decomposition of γ-cyclotrimethylene trinitramine (γ-RDX): relevance for shock wave initiation.

To elucidate the reactive behavior of RDX crystals at pressures and temperatures relevant to shock wave initiation, Raman spectroscopy and optical imaging were used to determine the pressure-temperature (P-T) stability and the decomposition of γ-RDX, the high pressure phase of RDX. Experiments were performed on single crystals in a diamond anvil cell at pressures from 6 to 12 GPa and at temperatures up to 600 K. Evidence for the direct decomposition of γ-RDX above 6 GPa, without the involvement of other phases, is provided. The upper limit of the P-T locus for the γ-RDX thermal decomposition was determined. A refined P-T phase diagram of RDX is presented that includes the current findings for γ-RDX. The static compression results are used to gain key insight into the shock initiation of RDX, including a determination of the RDX phase at decomposition and understanding the role of pressure and temperature in accelerating shock induced decomposition. This study has established the important role that γ-RDX plays in decomposition of RDX under static and shock compression conditions; thus theoretical modeling of RDX decomposition at high pressures and temperatures needs to incorporate the γ-phase response.

Density Functional Theory Calculations of Pressure Effects on the Structure and Vibrations of 1,1-Diamino-2,2-dinitroethene (FOX-7)

Pressure effects on the Raman vibrations of an energetic crystal FOX-7 (1, 1-diamino-2, 2-dinitroethene) were examined using density functional theory (DFT) calculations. High accuracy calculations were performed with a periodic plane-wave DFT method using norm-conserving pseudopotentials. Different exchange-correlation functionals were examined for their applicability in describing the structural and vibrational experimental data. It is shown that the PBE functional with an empirical dispersion correction by Grimme, PBE-D method, reproduces best the molecular geometry, unit cell parameters, and vibrational frequencies. Assignments of intramolecular Raman active vibrations are provided. The calculated pressure dependence of Raman shifts for the intramolecular and lattice modes were found to be in good agreement with the experimental data; in particular, the calculations predicted correctly a decrease of frequencies for the NH2 stretching modes with pressure. Also, in accord with experiments, the calculations indicated some instances of modes mixing/coupling with increasing pressure. This work demonstrates that the dispersion-corrected PBE functional can account for the structural and vibrational properties of FOX-7 crystal at ambient and high pressures.

High-pressure effects in pyrene crystals: vibrational spectroscopy.

The response of pyrene crystals to high pressure was examined using Raman and FTIR spectroscopies. Raman spectra of external and internal modes were measured up to 11 GPa. Changes in the external modes were observed at approximately 0.3 GPa, indicating the onset of a phase transition. We demonstrated that at this pressure pyrene I (P2(1)/a, 4 mol/unit cell) transforms to pyrene III (P2(1)/a, 2 mol/unit cell). Further increase of pressure produced a gradual broadening of the internal modes and an increase of fluorescence background, indicating the formation of another phase above 2.0 GPa. Irreversible chemical changes were observed upon gradual compression to 40 GPa. FTIR spectroscopy of the recovered product indicated a transformation of pyrene into an amorphous hydrogenated carbon (a-C:H) structure.

Density Functional Theory Calculations of Pressure Effects on the Vibrational Structure of α-RDX

Pressure effects on the vibrational structure of alpha-RDX were examined using density functional theory (DFT) up to 4 GPa. The calculated vibrational frequencies at ambient conditions are in better agreement with experimental data than are previous single molecule calculations. The calculations showed the following pressure-induced changes: (i) larger shifts for lattice modes and for internal modes associated with the CH(2) and NO(2) groups as compared to the pressure shifts for modes associated with the triazine ring, (ii) enhancement of mixing between different vibrations, for example, between NN stretching and CH(2) scissor, wagging, twisting vibrations, and (iii) increase in mixing between translational lattice vibrations and the NO(2) wagging vibrations, reducing the distinction between internal and lattice modes. The calculated volume and lattice constants at ambient pressure are larger than the experimental values, due to the inability of the present density functional approach to correctly account for van der Waals forces. Consequently, the pressure-induced frequency shifts of many modes deviate substantially from experimental data for pressures below 1 GPa. With increasing pressure, both the lattice constants and the frequency shifts agree more closely with experimental values.