To assess the influence of these extemporaneous deformations on reactivity, we applied MB-SMD as four-body biasing potentials to all nitro groups in TATB molecules packed into its crystal structure with two sets of simulations: torsional deformations which rotate the nitro groups with respect to the carbon ring, and out-of-plane deformations which bends nitro groups upwards/inwards. To study the effect of these molecular strains, we performed a series of reactive MD simulations where we increase the temperature at various rates, varying stiffnesses and set angles of the external potential, until significant reactions are observed. Quite interestingly, the simulation results follow Kissinger statistics quite closely, even for highly deformed systems. The characteristic reaction time, td, defined as the point of maximum internal energy, follows the expected rate dependence equation, see S6 of the SI.