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 α 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 Å 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 α -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.
Bibliographical noteFunding Information:
We greatly appreciate the support of the Office of Naval Research (Contract No. N00014-99-1-0393) and Project Officer Dr. Judah M. Goldwasser. We are also grateful to the University of Minnesota Grant-in-Aid program and its Undergraduate Research Opportunities Program. One of the authors (S.B.) wants to thank Kevin J. Boyd for many enlightening discussions.