Macroscale, mesoscale, and ab initio models of reacting shock physics are based, in their most general forms, on rate law descriptions of the chemical processes of interest. Reacting molecular dynamics simulations, by contrast, typically employ potential functions to model chemical reactions. An alternative reacting molecular dynamics formulation, employing nonholonomic Hamiltonian methods, models bonding-debonding as a kinetic process. Simulation results using this method are compared with experiment, for energy release and detonation products in HMX. The molecular dynamics simulations may be used to develop a macroscale, adiabatic model of the detonation chemistry.
Benchmarking Quasiclassical Mapping Hamiltonian Methods for Simulating Electronically Nonadiabatic Molecular Dynamics
