Determining the mechanical and decomposition properties of high energetic materials (α-RDX, β-HMX, and ε-CL-20) using a neural network potential

Molecular simulations of high energetic materials (HEMs) are limited by efficiency and accuracy. Recently, neural network potential (NNP) models have achieved molecular simulations of millions of atoms while maintaining the accuracy of density functional theory (DFT) levels. Herein, an NNP model covering typical HEMs containing C, H, N, and O elements is developed. The mechanical and decomposition properties of 1,3,5-trinitroperhydro-1,3,5-triazine (RDX), hexahydro-1,3,5-trinitro-1,3,5-triazine (HMX), and 2,4,6,8,10,12-hexanitrohexaazaisowurtzitane (CL-20) are determined by employing the molecular dynamics (MD) simulations based on the NNP model. The calculated results show that the mechanical properties of α-RDX, β-HMX, and ε-CL-20 agree with previous experiments and theoretical results, including cell parameters, equations of state, and elastic constants. In the thermal decomposition simulations, it is also found that the initial decomposition reactions of the three crystals are N-NO2 homolysis, corresponding radical intermediates formation, and NO2-induced reactions. This decomposition trajectory is mainly divided into two stages separating from the peak of NO2: pyrolysis and oxidation. Overall, the NNP model for C/H/N/O elements in this work is an alternative reactive force field for RDX, HMX, and CL-20 HEMs, and it opens up new potential for future kinetic study of nitramine explosives.

Temperature-dependent decomposition of the CL-20/MTNP cocrystal after phase separation

CL-20 (2,4,6,8,10,12-hexanitro-2,4,6,8,10,12-hexaazaisowurtzitane)-based cocrystals are attractive energetic cocrystals with a potential for high energy and low sensitivity, which account for nearly one-third of energetic cocrystals. The applications of cocrystal explosives require in-depth understanding of their thermal kinetics behaviors. Although thermal kinetics of the decomposition of CL-20-based cocrystals having no melting point have been studied, relevant research of CL-20-based cocrystals having a melting point, which are also the most frequently observed type, is still rare. In this study, the CL-20/MTNP (1-methyl-3,4,5-trinitropyrazole) cocrystal was chosen as a typical CL-20-based cocrystal having a melting point to investigate its thermal kinetics behavior. The thermal decomposition of CL-20/MTNP was identified to be a typical heterogeneous reaction with phase separation before decomposition. Due to the presence of intermolecular hydrogen bonds between CL-20 and molten MTNP after phase separation, the thermal decomposition behavior of CL-20/MTNP was strongly temperature-dependent. The complex decomposition reaction was separated into its three constituent pathways to simplify the kinetic analysis. On the basis of in-depth understanding of the decomposition process, the best functions of mechanism and kinetic parameters for each process of CL-20/MTNP decomposition were obtained using the model-fitting method. Finally, important thermal safety indicators, such as TMRad and SADT were simulated by combining the established kinetic models. This study provides further insights into the entire reaction process of the CL-20/MTNP cocrystal and would help in its better 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.

Molecular Dynamics Simulations of the Thermal Decomposition of RDX/HTPB Explosives

The addition of binders to energetic materials is known to complicate the thermal decomposition process of such materials. To assess this effect, the present work studied the thermal decomposition of cyclotrimethylene trinitramine (RDX)/hydroxy-terminated polybutadiene (HTPB) mixtures and of pure RDX over the temperature range of 2000-3500 K by combining the classical reaction and first-principles molecular dynamics methods. The incorporation of HTPB as a binder was found to significantly reduce the decomposition rate of RDX. At 3500 K, the decay rate constant of RDX in the RDX/HTPB system is 2.0141 × 10¹² s^-1, while it is 2.7723 × 10¹² s^-1 in the pure RDX system. However, the binder HTPB had little effect on the initial decomposition mechanism, which involved the rupture of N-NO₂ bonds to produce NO₂. The HTPB was predicted to undergo dehydrogenation and chain breaking. The free H resulting from these processes was predicted to react with low-molecular-weight intermediates generated by the RDX, resulting in greater equilibrium quantities of the final products H₂O and H₂ being obtained from the mixed system compared with pure RDX. HTPB-chain fragments were also found to combine with the primary RDX decomposition product NO₂ to inhibit the formation of N₂ and CO₂.

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.

Early thermal decay of energetic hydrogen- and nitro-free furoxan compounds: the case of DNTF and BTF

Exploration of the initial reactions of H-free and nitro-free energetic materials could enrich our understanding of the thermal decomposition mechanism of various energetic materials (EMs). In this work, two furoxan compounds, 3,4-dinitrofurazanfuroxan (DNTF) and benzotrifuroxan (BTF), were investigated to shed light on the decay mechanism of furoxan compounds based on the combination of self-consistent charge density functional tight binding and molecular dynamics simulations. The results show that DNTF and BTF decay via a unimolecular mechanism, and the transformation of the furoxan ring into a nitro group is suggested as a novel initial channel. Five initial steps of DNTF thermal decomposition are observed, including NO₂ loss and the N(O)-O bond cleavage of the central and peripheral rings. The bond cleavage of peripheral rings dominates the decay at low temperatures, while the central ring opening and C-NO₂ dissociation govern the high temperature decay. Besides, NO₂, CO and NO fragments are mainly yielded at high temperatures, while CO₃N₂ is dominant at low temperatures. The three-stage characteristic of the exothermic BTF decay is described under programmed heating conditions for the first time. Four initial steps of BTF thermal decomposition were identified, including furoxan ring opening reactions and the breakage of the 6-membered ring C-C bond. The cleavage of the N(O)-O bond is dominant in the initial step of BTF decomposition under different heating conditions, and the frequency increases with increasing temperature. In addition, the amounts of CON, ON and CO are higher at high temperatures, while C₂O₂N₂ shows an opposite trend. The findings of this work provide deep insights into the complicated sensitivity mechanism of EMs.

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.

Mechanism of the improvement of the energy of host-guest explosives by incorporation of small guest molecules: HNO3 and H2O2 promoted C-N bond cleavage of the ring of ICM-102

Host–guest materials exhibit great potential applications as an insensitive high-energy–density explosive and low characteristic signal solid propellant. To investigate the mechanism of the improvement of the energy of host–guest explosives by guest molecules, ReaxFF-lg reactive molecular dynamics simulations were performed to calculate the thermal decomposition reactions of the host–guest explosives systems ICM-102/HNO₃, ICM-102/H₂O₂, and pure ICM-102 under different constant high temperatures and different heating rates. Incorporation of guest molecules significantly increased the energy level of the host–guest system. However, the initial reaction path of the ICM-102 molecule was not changed by the guest molecules. The guest molecules did not initially participate in the host molecule reaction. After a period of time, the H₂O₂ and HNO₃ guest molecules promoted cleavage of the C–N bond of the ICM-102 ring. Stronger oxidation and higher oxygen content resulted in the guest molecules more obviously accelerating destruction of the ICM-102 ring structure. The guest molecules accelerated the initial endothermic reaction of ICM-102, but they played a more important role in the intermediate exothermic reaction stage: incorporation of guest molecules (HNO₃ and H₂O₂) greatly improved the heat release and exothermic reaction rate. Although the energies of the host–guest systems were clearly improved by incorporation of guest molecules, the guest molecules had little effect on the thermal stabilities of the systems.

Chemical reactions of a CL-20 crystal under heat and shock determined by ReaxFF reactive molecular dynamics simulations

Studying the chemical reactions of hexanitrohexaazaisowurtzitane (CL-20) under heat and shock is helpful to understand its sensitivity and shock initiation mechanism. In this work, several molecular dynamics simulations were performed under three different conditions: high temperature, high temperature and pressure, and shock. The formation and breakage of chemical bonds, changes of bond lengths, and initial reactions were analysed. It was found that the main small-molecule product of CL-20 during initial decomposition under the three different conditions was always NO2, but the generation pathways were different. At high temperatures, NO2 was generated by the direct cleavage of N-NO2 bonds. In contrast, high pressure and shock promoted the transfer of O atoms to N atoms connected to NO2, leading to the breakage of N-NO2 bonds. Almost all NO2 originated from the transfer of O atoms under the shock conditions.

Polymerization Effects on the Decomposition of a Pyrazolo-Triazine at high Temperatures and Pressures

4-amino-3-aminopyrazole-8-trinitropyrazolo-[5, 1-c] [1, 2, 4]triazine (PTX, C5H2N8O6) has good detonation performance, thermal stability and low mechanical sensitivity, which endow it with good development prospects in insensitive ammunition applications. To study the effects of polymerization on the decomposition of PTX, the reaction processes of PTX at different conditions were simulated by quantum chemistry and molecular dynamics methods. In this paper, the effects of polymerization on the decomposition of PTX were studied in terms of species information, reaction path of PTX, bond formation and bond cleavage, evolution of small molecules and clusters, and kinetic parameters at different stages. The results show that under the high-temperature and high-pressure conditions, the initial reaction path of unimolecular PTX in the thermal decomposition is mainly the cleavage of C-NO2 bonds. At the same time, there are many polymerization reactions in thermal decomposition process, which may greatly affect the reaction rate and path. The higher the degree of polymerization, the larger equilibrium value of potential energy, the less energy release of thermal decomposition. Compared with the activation energy of other explosives, the activation energy of PTX is higher than that of β-HMX and lower than that of TNT.

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.

Thermal stability mechanism via energy absorption by chemical bonds bending and stretching in free space and the interlayer reaction of layered molecular structure explosives

Layered molecular structure explosives have the characteristic of great thermal stability.

A flexible-molecule force field to model and study hexanitrohexaazaisowurtzitane (CL-20) – polymorphism under extreme conditions

The quantum-chemistry based force field (FF) developed for HMX by Smith and Bharadwaj (SB) [G. D. Smith and R. K. Bharadwaj, J. Phys. Chem. B, 1999, 103(18), 3570–3575] is transferred to another nitramine of different stoichiometry: hexanitrohexaazaisowurtzitane (CL-20 or HNIW). The modification of a single parameter alongside a very small number of add-ons related to carbon–carbon bonds, angles and dihedrals lead to two SB FF variants denoted SB-CL20 and SB-CL20 + CCNN. These flexible-molecule FFs should inherit the predictive capabilities of SB FF. For this purpose, we perform Molecular Dynamics simulations at ambient temperature and selected pressures. The modeled structures of the various CL-20 polymorphs are consistent with experimental data. Focusing on the ε-polymorph, we determine an equation of state which consolidates the general trend underpinned by most published results, and we confirm the increasing stiffness of the crystal under pressures up to 90 GPa. Moreover, we link some subtle pressure-induced changes of the elastic and structural properties to the flexibility and mobility of well-identified nitro groups. Finally, the simulations of the γ ↔ ζ phase transition suggest different multiple-step direct and reverse thermodynamic paths.

The quantum-chemistry based force field developed by Smith and Bharadwaj is transferred to hexanitrohexaazaisowurtzitane (CL20), revealing pressure-induced alterations of ε-CL20.