Reactive molecular dynamics simulation of the high-temperature pyrolysis of 2,2′,2′′,4,4′,4′′,6,6′,6′′-nonanitro-1,1′:3′,1′′-terphenyl (NONA)

2,2′,2′′,4,4′,4′′,6,6′,6′′-Nonanitro-1,1′:3′,1′′-terphenyl (NONA) is currently recognized as an excellent heat-resistant explosive. To improve the atomistic understanding of the thermal decomposition paths of NONA, we performed a series of reactive force field (ReaxFF) molecular dynamics simulations under extreme conditions of temperature and pressure. The results show that two distinct initial decomposition mechanisms are the homolytic cleavage of the C–NO(2) bond and nitro–nitrite (NO(2) → ONO) isomerization followed by NO fission. Bimolecular and fused ring compounds are found in the subsequent decomposition of NONA. The product identification analysis under finite time steps showed that the gaseous products are CO(2), N(2), and H(2)O. The amount of CO(2) is energetically more favorable for the system at high temperature or low density. The carbon-containing clusters are a favorable growth pathway at low temperatures, and this process was further demonstrated by the analysis of diffusion coefficients. The increase of the crystal density accelerates the decomposition of NONA judged by the analysis of reaction kinetic parameters and activation barriers. In the endothermic and exothermic stages, a 20% increase in NONA density increases the activation energies by 3.24 and 0.48 kcal mol(−1), respectively. The values of activation energies (49.34–49.82 kcal mol(−1)) agree with the experimental data in the initial decomposition stage.