Direct MD simulations of terahertz absorption and 2D spectroscopy applied to explosive crystals

A direct molecular dynamics simulation of the THz spectrum of a molecular crystal is presented. A time-dependent electric field is added to a molecular dynamics simulation of a crystal slab. The absorption spectrum is composed from the energy dissipated calculated from a series of applied pulses characterized by a carrier frequency. The spectrum of crystalline cyclotrimethylenetrinitramine (RDX) and triacetone triperoxide (TATP) were simulated with the ReaxFF force field. The proposed direct method avoids the linear response and harmonic approximations. A multidimensional extension of the spectroscopy is suggested and simulated based on the nonlinear response to a single polarized pulse of radiation in the perpendicular polarization direction.