Decreasing the density, RDX decomposition is less affected by the neighbor molecule, which is close to the isolated molecule in the previous DFT study . In summary, the effect of the neighbor molecule should be taken into consideration when studying the reaction mechanism of RDX and other EM in the condensed phase in practice. The present work develops a NN model for RDX at the level of ab initio calculations. Our NN model considers both atomic energy and force information from high-level DFT calculations. We demonstrate the accuracy of our NNP outperforms the widely used ReaxFF model. In particular, the ReaxFF model predicts unphysical atomic force compared to DFT calculations. This issue makes it unsuitable for the exploration of the processes in the reaction dynamics. Our NN model also exhibits a great computational efficiency, which is 10000 times faster than the DFT method, and 27 times faster than the ReaxFF method, allowing the possibility of investigating the complex reaction process of bulk RDX crystal from an atomic perspective. Thus, this study opens new opportunities for complex reactive systems to build reaction kinetics models.