Initial Decomposition of the Co-crystal of CL-20/TNT: Sensitivity Decrease under Shock Loading

Mechanical behaviors of CL-20 under an impact loading: A molecular dynamics study

Study on the dynamic process of CL-20 crystal under impact is critical for the safe utilization of energetic materials under extreme conditions. Herein, the mechanical and structural evolution of CL-20 under the impact of a diamond ball is investigated by using molecular dynamics simulation. The considerations are given to the effect of different impact velocity, impact direction and impact angle. It is found that a high impact velocity results in a large indentation depth and force, as well as more significant energy transition and the formation of a large number of molecular fragments. Moreover, CL-20 exhibits weak anisotropy along different impact directions due to the crystalline distribution anisotropy. Furthermore, the mechanical response of CL-20 is angle-dependent, which is caused by the discrepancy in local molecular re-arrangement. These results may enhance the understanding of the mechanical behavior of CL-20 and promote its wide applications.

Effects of Nanoparticle Size on the Thermal Decomposition Mechanisms of 3,5-Diamino-6-hydroxy-2-oxide-4-nitropyrimidone through ReaxFF Large-Scale Molecular Dynamics Simulations

ReaxFF-lg molecular dynamics method was employed to simulate the decomposition processes of IHEM-1 nanoparticles at high temperatures. The findings indicate that the initial decomposition paths of the nanoparticles with different sizes at varying temperatures are similar, where the bimolecular polymerization reaction occurred first. Particle size has little effect on the initial decomposition pathway, whereas there are differences in the numbers of the species during the decomposition and their evolution trends. The formation of the hydroxyl radicals is the dominant decomposition mechanism with the highest reaction frequency. The degradation rate of the IHEM-1 molecules gradually increases with the increasing temperature. The IHEM-1 nanoparticles with smaller sizes exhibit greater decomposition rate constants. The activation energies for the decomposition are lower than the reported experimental values of bulk explosives, which suggests a higher sensitivity.

Deep Potential Molecular Dynamics Study of Chapman-Jouguet Detonation Events of Energetic Materials

Detonation of energetic materials (EMs) is of great importance for military applications, while the understanding of detailed events and mechanisms for the detonation process is scarce. In this study, the first deep neural network potential NNP_Shock for molecular dynamics (MD) simulation of shock-induced detonation of EMs was generated based on a deep potential model, providing DFT accuracy but 106 times the computational efficiency. On this basis, we employ our deep potential to perform MD simulations of shock-induced detonation of high-performance EM material 2,4,6,8,10,12-hexanitro-2,4,6,8,10,12-hexaazaisowurtzitane (CL-20, C6H6N12O12) and obtain the theoretical Chapman-Jouguet (C-J) detonation velocities and pressures directly by multiscale shock technique (MSST) for the first time, which are in good agreement with experiment. In addition, the Hugoniot curves and initial reaction mechanisms were successfully obtained. Therefore, the NNP_Shock potential is competent in research of the detonation performance and shock sensitivity of CL-20, and the method can also be transplanted to studies of other EMs.

Structural evolution of CL-20/DNB cocrystals at high temperature: Phase transition and kinetics of thermal decomposition

As a typical new energetic material, CL-20/DNB cocrystals have been recognized as a promising explosive owing to their excellent comprehensive performance. The thermal decomposition behavior, structural evolution and dynamic process of CL-20/DNB cocrystals under high temperature were studied by means of thermogravimetric differential heating, X-ray diffraction, Raman spectroscopy to gain insight into the cocrystal materials. The study found that the decomposition of CL-20/DNB cocrystal is a heterogeneous process accompanied by the sublimation of DNB and structural change of CL-20. The phase transition of β → γ-CL-20 was observed at 120 °C. The kinetics of decomposition and the mechanism of micro structural evolution on CL-20/DNB cocrystals with heating were revealed. The primary NO⋯H hydrogen bonds of the cocrystal are broken, accompanied by the melting of DNB in the temperature range of 100-120 °C. Subsequently, the DNB single component decomposes completely, leading to lattice collapse of cocrystal; simultaneously, CL-20 undergoes a transition process from β phase to γ phase. Ultimately, γ-CL-20 gradually decomposes with increasing temperature. The activation energy of cocrystal is also obtained as 129 ± 10 kJ/mol. The understanding of cocrystal explosive was deepened and the further application was promoted.

Mechanism of the Impact-Sensitivity Reduction of Energetic CL-20/TNT Cocrystals: A Nonequilibrium Molecular Dynamics Study

High-energy low-sensitivity explosives are research objectives in the field of energetic materials, and the formation of cocrystals is an important method to improve the safety of explosives. However, the sensitivity reduction mechanism of cocrystal explosives is still unclear. In this study, CL-20/TNT, CL-20 and TNT crystals were taken as research objects. On the basis of the ReaxFF-lg reactive force field, the propagation process of the wave front in the crystals at different impact velocities was simulated. The molecular dynamics data were used to analyze the molecular structure changes and initial chemical reactions, and to explore the sensitivity reduction mechanism of the CL-20/TNT cocrystal. The results showed that the chemical reaction of the CL-20/TNT cocrystal, compared with the CL-20 single crystal, is different under different impact velocities. At an impact velocity of 2 km/s, polymerization and separation of the component molecules weakened the decomposition of CL-20. At an impact velocity of 3 km/s, the decay rates of CL-20 and TNT in the cocrystal decreased, and the intermediate products were enhanced, such as nitrogen oxides. At an impact velocity of 4 km/s, the cocrystal had little effect on the decay rates of the molecules and formation of CO₂, but it enhanced formation of N₂ and H₂O. This may explain the reason for the impact-sensitivity reduction of the CL-20/TNT cocrystal.

Molecular conformation evolutions of trans HNS under high pressure

The molecular conformation evolution of Hexanitrostilbene (HNS) under high pressure was systematically investigated using Raman and Fourier transform infrared (FTIR) spectroscopy. The vibration modes of HNS associated with C-H, nitro groups, CC and the ring have been analyzed and clarified in detail under ambient conditions. trans-HNS is symmetrically distributed about -CHCH-, and six nitro groups are symmetrically distributed under ambient conditions. Two molecular conformation changes of HNS were observed at 1.4 GPa and 5 GPa due to the variations of hydrogen-bond interaction between C-H (in the ring) and N-O and the distortion of trans olefin, respectively. The hydrogen-bond interaction between C-H (in the ring) and N-O strengthened at 1.4 GPa. It induced the degenerated symmetry of the nitro groups and the Raman changes of ν_as (NO₂), ν(CC), ν(C-C) and ν(C-H). In addition, the nonplanarity property of HNS and the sensitivity of trans olefin to pressure promoted the deformation of trans olefin, as well as the hydrogen bond interaction between C-H (in trans olefin) and N-O at about 5 GPa. When further loading pressure on HNS, the variations in the hydrogen-bond interaction between C-H and N-O restricted the vibrations of C-H, NO₂ and the ring. It blocked the nonradiative pathway and activated the strong fluorescent background in the Raman spectra as the pressure increased above 5.7 GPa. These current results reveal that there is no structural transformation and only conformational changes under high pressure for HNS.

CL-20/TNT decomposition under shock: cocrystalline versus amorphous

The cocrystallization strategy is considered to be an effective means to adjust the properties of explosives. Nevertheless, the underlying mechanism of the effect of the special cocrystal structure on the decomposition process is not clear enough. The present work compares the response processes of a CL-20/TNT cocrystal structure and an amorphous structure under shock waves with different velocities. The thermodynamic evolution, reactant decay, product formation, main initial reactions and cluster evolution are analyzed. As a result, we find that the amorphous structure is easier to compress than the cocrystal structure, achieving higher stress and temperature. These thermodynamic parameters have a strong correlation. For the amorphous structure, the chemical reaction of the system is more intense, the reactants decay faster, the products are more abundant, and the intermediate products can complete the transformation to stable products earlier. Furthermore, NO2 is the most important intermediate product, and its quantitative change can directly reflect the reaction process. The amorphous structure is more prone to decomposition reaction, and the cocrystal structure is more prone to polymerization reaction. The cluster size in the amorphous structure is smaller and more conducive to decomposition.

The role of electric field on decomposition of CL-20/HMX cocrystal: A reactive molecular dynamics study

Electric field can initiate decomposition or detonation of explosives, but underlying mechanism is unclear. Here, we performed ReaxFF molecular dynamics simulation for decomposition of a cocrystal, formed by 2,4,6,8,10,12-hexanitro-2,4,6,8,10,12-hexaazaisowurtzitane (CL-20) and 1,3,5,7-tetranitro-1,3,5,7-tetrazacyclooctane (HMX), solely induced by electric field. A new analytical method was proposed to obtain detailed decomposition mechanism. Results show that electric fields play important roles in decomposition of CL-20/HMX cocrystal, such as heating the system and causing the explosive to decompose. Strong constant field makes CL-20 molecules in the cocrystal decompose at significantly lower temperature, which greatly increases sensitivity. This is ascribed to the distinct decomposition mechanism that CN bond rupture dominates the initial step of CL-20's decomposition. Contrarily, oscillating field has a stronger heating effect but weaker influence on sensitivity. Moreover, HMX exhibits desensitizing effect in CL-20/HMX cocrystal under electric field. These results enhance our understanding of sensitivity mechanism beyond mechanical stimuli in explosives.

A quantum-based molecular dynamics study of the ICM-102/HNO3 host-guest reaction at high temperatures

The contradiction between energy and safety of explosives is better balanced by the host-guest inclusion strategy. Understanding the reaction mechanism of the host-guest explosive is necessary. To deeply analyze the role of the small guest molecules in the host-guest system, a quantum-based molecular dynamics method was used to calculate the initial decomposition reaction of the new host-guest explosive ICM-102/HNO3 against the pure ICM-102 at several high temperatures. The incorporation of HNO3 had no significant influence on the initial decomposition step of ICM-102. Conversely, the earliest intramolecular hydrogen transfer reaction is delayed partly because the H and O atoms of HNO3 connect with the O and H atoms of ICM-102, respectively. As the reaction proceeds, guest molecules get heavily involved in the reaction and increase the reaction rate. The generation rate and quantity of the small oxidizing molecules in the final product were increased significantly in the ICM-102/HNO3 system. These mechanisms revealed that HNO3 molecules inhibit the early stages of the initial decomposition of ICM-102 to some extent, and play an important role in accelerating the decomposition subsequently.

Decomposition mechanism scenarios of CL-20 co-crystals revealed by ReaxFF molecular dynamics: similarities and differences

Understanding the similarities and differences of decomposition mechanisms of CL-20 and its cocrystals is of great interest for practical applications of CL-20 cocrystals. The responses of CL-20 cocrystals to thermal stimulus were investigated using ReaxFF molecular dynamics simulations of two representative cocrystals, CL-20/HMX and CL-20/TNT, under adiabatic conditions and comparing to the baseline system of pure CL-20. The comprehensive chemical details were revealed with the aid of the unique code of VARxMD. The three CL-20-involved reactive systems all exhibit a distinct three-stage character during adiabatic decomposition when using the double peaks of the major intermediate NO₂ amount as the boundary. By taking advantage of the three-stage classification, a clear scenario for the similar stimulus-response of the CL-20 cocrystals can be elucidated for the dominant primary decomposition of CL-20 in stage I and the transition of favored chemical mechanisms from the generation of intermediates/radicals in stage II into their consumption to form stable products in stage III. The similar chemical behaviors are rooted in the dominance of CL-20 chemistry in the initial response of its cocrystals to thermal stimulus. The prolonged reaction zone uncovers the slowed decomposition kinetics of CL-20/HMX and CL-20/TNT, which is associated with the altered kinetics of CL-20 decomposition specifically by N-NO₂ bond scission and CL-20 skeleton decay. The retarded CL-20 decomposition in its cocrystals consequently results in more moderate self-heating and less violent exothermic reactions that agrees with the experimental observations of improved stability and damaged detonation performance of CL-20 cocrystals, particularly for CL-20/TNT. The results obtained in this work suggest that ReaxFF MD simulations can provide useful insight for the modulated chemical properties of varied CL-20 cocrystals.

Design, preparation, characterization and formation mechanism of a novel kinetic CL-20-based cocrystal

2,4,6,8,10,12-Hexanitrohexaazaisowurtzitane (CL-20)-based cocrystals have gained increasing attention as a means of obtaining insensitive high explosives. However, the design of ideal candidates for these cocrystals remains difficult. This work compares the crystal energies of the CL-20-dinitrobenzene (DNB) and CL-20-2,4,6-trinitrotoluene (TNT) cocrystals with those of the respective pure coformers. The results indicate that the cocrystal formation is driven by the differences in the energies of the cocrystals and the coformers. Furthermore, analysis via Hirshfeld surfaces and two-dimensional fingerprint plots confirms that the O...O, O...H, O...N and C...O interactions were the main force for stabilizing the CL-20-based cocrystal structure. Based on these findings, a novel energetic-energetic cocrystal of CL-20-2,4,6-trinitrophenol (TNP) was designed and prepared by means of a rapid method for solvent removal. The crystal structure was investigated via powder X-ray diffraction methods, solid-state nuclear magnetic resonance spectroscopy and Fourier transform infrared spectroscopy. The results revealed that the O-H...O hydrogen bonding interaction between the phenolic hydroxyl group of TNP and nitro groups of CL-20, as well as nitro...π, nitro...nitro and O_NO2...π(N)_NO2 interactions, based on the benzene ring and nitro groups, are the main interactions occurring in the cocrystal.