Effects of Defects on Thermal Decomposition of HMX via ReaxFF Molecular Dynamics Simulations

Abstract

Effects of molecular vacancies on the decomposition mechanisms and reaction dynamics of condensed-phase β-HMX at various temperatures were studied using ReaxFF molecular dynamics simulations. Results show that three primary initial decomposition mechanisms, namely, N-NO(2) bond dissociation, HONO elimination, and concerted ring fission, exist at both high and lower temperatures. The contribution of the three mechanisms to the initial decomposition of HMX is influenced by molecular vacancies, and the effects vary with temperature. At high temperature (2500 K), molecular vacancies remarkably promote N-N bond cleavage and concerted ring breaking but hinder HONO formation. N-N bond dissociation and HONO elimination are two primary competing reaction mechanisms, and the former is dominant in the initial decomposition. Concerted ring breaking of condensed-phase HMX is not favored at high temperature. At lower temperature (1500 K), the most preferential initial decomposition pathway is N-N bond dissociation followed by the formation of NO(3) (O migration), although all three mechanisms are promoted by molecular vacancies. The promotion effect on concerted ring breaking is considerable at lower temperature. Products resulting from concerted ring breaking appear in the defective system but not in the perfect crystal. The mechanism of HONO elimination is less important at lower temperature. We also estimated the reaction rate constant and activation barriers of initial decomposition with different vacancy concentrations. Molecular vacancies accelerate the decomposition of condensed-phase HMX by increasing the reaction rate constant and reducing activation barriers.

Compressive Shear Reactive Molecular Dynamics Studies Indicating That Cocrystals of TNT/CL-20 Decrease Sensitivity

“…Here the amount of NO 2 reached a maximum of about 9 ps and then started to decrease, indicating that some secondary reactions are already occurring in the CL-20 system. After this point one can expect, a quick consumption of NO 2 molecules with a faster production of stable reaction products, such as N 2 , H 2 O, and CO 2 , as shown in our previous cook-off simulations 26,32 .…”

ReaxFF reactive molecular dynamics on silicon pentaerythritol tetranitrate crystal validates the mechanism for the colossal sensitivity

“…This is because similar ReaxFF reactive force fields for other materials have been validated to predict accurately both the reactivity of bonds and mechanical properties of condensed phases. 16,[22][23][24][25][26][27][28][29][30][31] The studies of anisotropic sensitivity of PETN and HMX, 16,24 thermal decomposition of HMX (cyclotetramethylene tetranitramine), TATB (triamino-trinitrobenzene), and RDX, [25][26][27] shock dynamics of RDX and PBX (plastic bonded explosives), [28][29][30][31] and so forth using ReaxFF-RMD lead to the results in accordance with available experimental data, making it practical to describe chemical reactions occurring under various conditions during the large scale dynamical processes involving millions of atoms with currently available computational facilities.…”

Formation and growth mechanisms of natural metastable twin boundary in crystalline β-octahydro-1,3,5,7-tetranitro-1,3,5,7-tetrazocine: a computational study

“…8a and b). Zhou et al 3 found that the molecular vacancies accelerate the decomposition of HMX by increasing the reaction rate constant and reducing activation barriers. The VBM of (001)/[010] TB_2 is contributed by both the aNO 2 and eNO 2 groups, while the CBM is contributed dominantly by the aNO 2 groups.…”