Nonreactive molecular dynamics force field for crystalline hexahydro-1,3,5-trinitro-1,3,5 triazine

An empirical nonreactive force field has been developed for molecular dynamics (MD)/Monte Carlo simulation of the formation, diffusion, and agglomeration of point defects in the crystal lattice of the alpha modification of hexahydro-1,3,5-trinitro-1,3,5 triazine (RDX) using flexible molecules. Bond stretching and angle bending are represented by Morse and harmonic functions, and torsion by a truncated cosine series. Nonbonded interactions, both inter- and intramolecular, are described by Buckingham potentials separately parametrized. Intermolecular electrostatic interactions are treated via a Coulomb term coupled with a smooth 15.0 A cutoff radius. Parameters were taken in part from earlier published works and were determined partly by fitting to known molecular and crystal properties of RDX. In MD simulations at constant pressure and temperature, the model was able to stabilize and maintain the correct crystal structure, symmetry, and molecular conformation of alpha-RDX. Vibrational frequencies, lattice binding energy and dimensions, coefficients of thermal expansion, and several unusually short intermolecular distances are all reproduced in satisfactory agreement with experimental data.