Molecular dynamics simulations of HMX crystal polymorphs using a flexible molecule force field

Molecular dynamics simulations using a recently developed quantum chemistry-based atomistic force field [J. Phys. Chem. B 103 (1999) 3570] were performed in order to obtain unit cell parameters, coefficients of thermal expansion, and heats of sublimation for the three pure crystal polymorphs of octahydro-l,3,5,7-tetranitro-l,3,5,7-tetrazocine (HMX). The predictions for β-, α-, and δ-HMX showed good agreement with the available experimental data.