“… where D e σ , D e π , and D e ππ are the bond parameters for σ-bonds, π-bonds, and π-π bonds, respectively; p be 1 and p be 2 are the fitting parameters. ReaxFF [ 98 , 99 , 100 ] was initially developed for simulating combustion of hydrocarbon and reactions of energetic compounds, and the functional groups such as nitro group are well-parameterized, thus various energetic compounds such as RDX, HMX, CL-20, PETN, NM, TNT, TATB, and TATP were applied [ 115 , 116 , 117 , 118 , 119 , 120 , 121 , 122 , 123 , 124 , 125 , 126 , 127 , 128 , 129 , 130 , 131 , 132 , 133 , 134 , 135 , 136 ]. The ReaxFF-based MD simulations provide an efficient way to investigate the chemical reactions and other properties of various molecular crystals under particular conditions.…”