First principles molecular dynamics simulation and thermal decomposition kinetics study of CL-20

2,4,6,8,10, 12-hexanitro-2,4,6,8,10, 12-hexazepane (CL-20) is a new energetic material with high performance and low sensitivity. In-depth study of the thermal decomposition mechanism of CL-20 is a necessary condition to improve its performance, ensure its safety, and optimize its application. On the basis of a large number of empirical force fields used in molecular dynamics simulation in the past, the machine learning augmented first-principles molecular dynamics method was used for the first time to simulate the thermal decomposition reaction of CL-20 at 2200 K, 2500 K, 2800 K, and 3000 K isothermal temperature. The main stable resulting compounds are N₂, CO₂, CO, H₂O, andH₂, where CO₂ and H₂O continue to decompose at higher temperatures. The initial decomposition pathways are denitration by N-N fracture, ring-opening by C-N bond fracture, and redox reaction involving NO₂ and CL-20. After ring opening, two main compounds, fused tricyclic pyrazine and azadicyclic, were formed, which were decomposed continuously to form monocyclic pyrazine and pyrazole ring structures. The most common fragments formed during decomposition are those containing two, three, four, and six carbons. The formation rule and quantity of main small molecule intermediates and resulting stable products under different simulated temperatures were analyzed.

Based on ab initio Bayesian active learning algorithm, efficient and accurate prediction of CL-20 is made using the dynamic machine learning function of Vienna Ab-Initio Simulation Package (VASP), which constructs the energy potential surface by learning a large number of data based on AIMD calculations. The result is a machine learning force field (MLFF). Then the molecular dynamics of CL-20 was simulated using the trained MLFF model. PAW pseudopotentials and generalized gradient approximation (GGA), namely, Perdew-Burke-Ernzerhof (PBE) functional, are used in the calculation. The plane wave truncation energy (ENCUT) is set to 550 eV, and using the Gaussian broadening, the thermal broadening size of the single-electron orbital is 0.05 eV. A van der Waals revision of the system with Grimme Version 3. The energy convergence accuracy (EDIFF) of electron self-consistent iteration is set to 1E-5 eV and 1E-6 eV, respectively. The two-step structure optimization is carried out using 1'1'1 k point grid and conjugate gradient method. The ENCUT was changed to 500 eV and EDIFF to 1E-5 eV, and NVT integration (ISIF = 2) of Langevin thermostat was used for machine learning force field training and AIMD simulation of the system.