Intermolecular Potential for the Hexahydro-1,3,5-trinitro-1,3,5-s-triazine Crystal (RDX): A Crystal Packing, Monte Carlo, and Molecular Dynamics Study

Abstract

We have developed an intermolecular potential that describes the structure of the α-form of the hexahydro-1,3,5-trinitro-1,3,5-s-triazine (RDX) crystal. The potential is composed of pairwise atom−atom (6-exp) Buckingham interactions and charge−charge interactions. The parameters of the Buckingham repulsion−dispersion terms have been determined through a combination of nonlinear least-squares fitting to observed crystal structures and lattice energies and trial-and-error adjustment. Crystal packing calculations were performed to determine the equilibrium crystallographic structure and lattice energy of the model. There are no significant differences in the geometrical structures and crystal energies resulting from minimization of the lattice energy with and without symmetry constraints. Further testing of the intermolecular potential has been done by performing symmetry-constrained isothermal−isobaric Monte Carlo simulations. The properties of the crystal (lattice dimensions, molecular orientation, and lat...