Reaction rates in shocked nitromethane from density functional tight binding simulations

The chemistry of energetic materials (EM) is characterized by rapid exothermic reactions that lead to dramatic increase of the pressure and temperature on the pico- to nanosecond timescales. Under these conditions, experiments have struggled to provide detailed insights into early and intermediate processes, and simulations have thus become a valuable tool to help interpret experiments and parameterize mesoscale models. We have performed molecular dynamics (MD) cook-off simulations of nitromethane under pressure with DFTB, a parameterized form of DFT that allows to simulate systems with hundreds of atoms, over hundreds of picoseconds, with explicit treatment of the electronic interactions and an accuracy close to that of DFT. We find drastically different times-to-explosion, even for the same initial T/P conditions, due to multiple complex and competitive chemical pathways. However, a simple effective reaction rate can be extracted, as long as multiple simulations are performed at each T/P to account for the stochastic component of detonation chemistry in EM. We built a two-step model for NM detonation that can be compared to experimental results and used in higher scale models.