Reactive molecular dynamics simulations of shock through a single crystal of pentaerythritol tetranitrate

Large-scale molecular dynamics simulations and the reactive force field ReaxFF were used to study shock-induced initiation in crystalline pentaerythritol tetranitrate (PETN). In the calculations, a PETN single crystal was impacted against a wall, driving a shockwave back through the crystal in the [100] direction. Two impact speeds (4 and 3 km/s) were used to compare strong and moderate shock behavior. The primary difference between the two shock strengths is the time required to exhibit the same qualitative behaviors with the lower impact speed lagging behind the faster impact speed. For both systems, the shock velocity exhibits an initial deceleration due to onset of endothermic reactions followed by acceleration due to the onset of exothermic reactions. At long times, the shock velocity reaches a steady value. After the initial deceleration period, peaks are observed in the profiles of the density and axial stress with the strongly shocked system having sharp peaks while the weakly shocked system developed broad peaks due to the slower shock velocity acceleration. The dominant initiation reactions in both systems lead to the formation of NO(2) with lesser quantities of NO(3) and formaldehyde also produced.