Inspired by the recent cocrystallization and theory of energetic materials, we theoretically investigated the intermolecular vibrational energy transfer process and the non-covalent intermolecular interactions between explosive compounds. The intermolecular interactions between 2,4,6-trinitrotoluene (TNT) and 2,4,6,8,10,12-hexanitro-2,4,6,8,10,12-hexaazaisowurtzitane (CL-20) and between 1,3,5,7-tetranitro-1,3,5,7-tetrazocane (HMX) and CL-20 were studied using calculated two-dimensional infrared (2D IR) spectra and the independent gradient model based on the Hirshfeld partition (IGMH) method, respectively. Based on the comparison of the theoretical infrared spectra and optimized geometries with experimental results, the theoretical models can effectively reproduce the experimental geometries. By analyzing cross-peaks in the 2D IR spectra of TNT/CL-20, the intermolecular vibrational energy transfer process between TNT and CL-20 was calculated, and the conclusion was made that the vibrational energy transfer process between CL-20 and TNTII (TNTIII) is relatively slower than between CL-20 and TNTI. As the vibration energy transfer is the bridge of the intermolecular interactions, the weak intermolecular interactions were visualized using the IGMH method, and the results demonstrate that the intermolecular non-covalent interactions of TNT/CL-20 include van der Waals (vdW) interactions and hydrogen bonds, while the intermolecular non-covalent interactions of HMX/CL-20 are mainly comprised of vdW interactions. Further, we determined that the intermolecular interaction can stabilize the trigger bond in TNT/CL-20 and HMX/CL-20 based on Mayer bond order density, and stronger intermolecular interactions generally indicate lower impact sensitivity of energetic materials. We believe that the results obtained in this work are important for a better understanding of the cocrystal mechanism and its application in the field of energetic materials.

Particle shapes significantly affect viscosity and flow behavior of energetic materials, and therefore affect their packability and processability. This study presents a computational geometry framework for automatically quantifying two-dimensional (2D) and three-dimensional (3D) particle shapes of energetic materials. A specimen by mixing three typical energetic materials including HMX (octahydro-1,3,5,7-tetranitro-1,3,5,7-tetrazocine), RDX (1,3,5-Trinitroperhydro-1,3,5-triazine) and AP (Ammonium Perchlorate) particles is used in this study. This specimen is scanned by high-resolution X-ray computed tomography (X-ray CT), yielding a volumetric image. An improved watershed analysis algorithm is used to process the volumetric image to identify individual 3D particles. The stereology sampling method is used to obtain 2D projections of 3D particles. Computational geometry techniques are developed by this study to analyze 2D particle projections and 3D particle geometries to compute seven commonly used shape descriptors, including convexity, circularity, intercept sphericity, area sphericity, diameter sphericity, circle ratio sphericity, and surface area sphericity. Results show that those different shape descriptors of energetic materials can be divided into three groups based on their numerical ranges. This study also evaluates the effectiveness and accuracy of 2D shape descriptors for quantifying the true 3D shapes. The inconsistent characterization results between 2D and 3D shape descriptors suggest that researchers should be cautious when using 2D images to characterize 3D particle shapes of energetic materials. The computational geometry framework and particle shape analysis results presented in this study can be potentially useful in numerical modeling, experimental analysis, and theoretical investigation for energetic materials.

The external electric field (E-field), which is an important stimulus, can change the decomposition mechanism and sensitivity of energetic materials. As a result, understanding the response of energetic materials to external E-fields is critical for their safe use. Motivated by recent experiments and theories, the two-dimensional infrared (2D IR) spectra of 3,4-bis (3-nitrofurazan-4-yl) furoxan (DNTF), which has a high energy, a low melting point, and comprehensive properties, were theoretically investigated. Cross-peaks were observed in 2D IR spectra under different E-fields, which demonstrated an intermolecular vibration energy transfer; the furazan ring vibration was found to play an important role in the analysis of vibration energy distribution and was extended over several DNTF molecules. Measurements of the non-covalent interactions, with the support of the 2D IR spectra, indicated that there were obvious non-covalent interactions among different DNTF molecules, which resulted from the conjugation of the furoxan ring and the furazan ring; the direction of the E-field also had a significant influence on the strength of the weak interactions. Furthermore, the calculation of the Laplacian bond order, which characterized the C-NO₂ bonds as trigger bonds, predicted that the E-fields could change the thermal decomposition process of DNTF while the positive E-field facilitates the breakdown of the C-NO₂ in DNTFⅣ molecules. Our work provides new insights into the relationship between the E-field and the intermolecular vibration energy transfer and decomposition mechanism of the DNTF system.

As a new type of energetic material, cocrystal explosives demonstrate many excellent properties, such as high energy density and low sensitivity, due to the interaction between the molecules of the two components. The known decomposition temperature is 235 °C for CL-20/HMX cocrystals at a faster heating rate. CL-20 molecules could separate from the cocrystal matrix and decompose at a higher temperature, much lower than the decomposition temperature. The current work provided deep insight into the isothermal structural evolution of CL-20/HMX cocrystals with slow roasting at 190 °C. We found that the initial decomposition originates from separating CL-20 molecules from the surface along the (010) plane of the cocrystals. The gas products, such as NO₂ and NO, escape from the largest exposed surface of the (010) plane and generates microbubbles and microholes. At the same time, the residual HMX molecules form δ-phase HMX crystals and shrink the volume by 72%. By increasing the time held at 190 °C, the decomposition of CL-20 molecules and recrystallization of the residual HMX molecules form a gully-like structure on the (010) plane of the CL-20/HMX cocrystal. After a long time at 190 °C, the CL-20 component completely decomposes, and all HMX molecules recrystallize in the δ-HMX form. The interaction between HMX and CL-20 molecules makes the decomposition rate of the CL-20/HMX cocrystal much slower than that of the CL-20 pure crystal with a similar decomposition activation energy during isothermal heating. This work can help to deeply understand the safety of CL-20/HMX cocrystal explosives at a temperature lower than the recognized decomposition temperature.