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

Sylke Boyd, Matthew Gravelle, Peter Politzer

    Research output: Contribution to journalArticlepeer-review

    38 Scopus citations


    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.

    Original languageEnglish (US)
    Article number104508
    JournalJournal of Chemical Physics
    Issue number10
    StatePublished - 2006

    Bibliographical note

    Funding 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.

    Fingerprint Dive into the research topics of 'Nonreactive molecular dynamics force field for crystalline hexahydro-1,3,5-trinitro-1,3,5 triazine'. Together they form a unique fingerprint.

    Cite this