Thermal stability ranking of energetic crystals via a neural network potential-enabled MD simulation protocol

• Developed an improved MD protocol for predicting thermal stability of EMs. • Particle structures reduce T d error by up to 400 K than periodic structures. • Lower heating rates (0.001 K/ps) reduce T d error to as low as 80 K. • Achieved strong correlation (R 2 = 0.969) between MD results and experiments. • Proposed a correction model to bridge MD-predicted and experimental T d values. Assessing the thermal stability of energetic materials (EMs) remains challenging due to the limitations of traditional experimental and computational methods. This study develops an optimized molecular dynamics (MD) protocol based on a neural network potential (NNP) to enable reliable quantitative prediction of EM thermal stability. Key improvements include the use of nanoparticle models and reduced heating rates.​ Systematic investigations on RDX show that nanoparticle structures mitigate decomposition temperature ( T d ) overestimation compared to periodic models, with surface effects dominating over particle size. Lower heating rates (e.g., 0.001 K/ps) further reduce deviation, bringing RDX T d within 80 K of experimental values (vs. >400 K in conventional simulations). Kissinger analysis of the heating rate- T d relationship supports the feasibility of optimizing heating rates to align with experimental T d .​ Applied to eight representative EMs, the optimized protocol yields thermal stability rankings in excellent agreement with experiments (R 2 =0.96), outperforming traditional periodic models (R 2 =0.85). This work establishes a robust computational framework for EM thermal stability evaluation, particularly valuable in data-limited scenarios.