Molecular dynamics simulations of shock waves in oriented nitromethane single crystals: plane-specific effects.

Molecular dynamics simulations of supported shock waves (shock pressure P(s) ∼ 15 GPa) propagating along the [110], [011], [101], and [111] directions in crystalline nitromethane initially at T = 200 K were performed using the nonreactive Sorescu-Rice-Thompson force field [D. C. Sorescu, B. M. Rice, and D. L. Thompson, J. Phys. Chem. B 104, 8406 (2000)]. These simulations, combined with those from a preceding study of shocks propagating along [100], [010], and [001] directions in nitromethane for similar conditions of temperature and shock pressure [L. He, T. D. Sewell, and D. L. Thompson, J. Chem. Phys. 134, 124506 (2011)], have been used to study the post-shock relaxation phenomena. Shocks along [010] and [101] lead to a crystal-crystal structure transformation. Shocks propagating along [011], [110], [111], [100], and [001] exhibit plane-specific disordering, which was characterized by calculating as functions of time the 1D mean square displacement (MSD), 2D radial distribution function (RDF), and 2D orientation order parameter P(2)(θ) in orthogonal planes mutually perpendicular to the shock plane; and by calculating as functions of distance behind the shock front the Cartesian components of intermolecular, intramolecular, and total kinetic energies. The 2D RDF results show that the structural disordering for shocks along [100], [110], and [111] is strongly plane-specific; whereas for shocks along [001] and [011], the loss of crystal structural order is almost equivalent in the orthogonal planes perpendicular to the shock plane. Based on the entire set of simulations, there is a trend for the most extensive disordering to occur in the (010) and (110) planes, less extensive disordering to occur in the (100) plane, and essentially no disordering to occur in the (001) plane. The 2D P(2)(θ) and 1D MSD profiles show, respectively, that the orientational and translational disordering is plane-specific, which results in the plane-specific structural disordering observed in the 2D RDF. By contrast, the kinetic energy partitioning and redistribution do not exhibit plane specificity, as shown by the similarity of spatial profiles of the Cartesian components of the intermolecular, intramolecular, and total kinetic energies in orthogonal planes perpendicular to the shock plane.