The RDX decomposition process is investigated with the final NNP. The simulation is performed under the NVT ensemble with a 2x2x2 supercell. The equations of motion are integrated by the velocity Verlet method using periodic boundary conditions. A Nose-Hoover thermostat is applied with an equilibrium temperature of 1000, 1250, 1500, 1750, 2000, 2250, 2500 K and a dump parameter of 20 fs. A 100-ps MD simulation is performed with a time step of 0.l fs. Three parallel simulations are performed for each NVT MD simulation to ensure the statistical significance of the simulated results. The ReacNetGenerator is used to extract species and reactions from the MD trajectories.