Molecular dynamics simulations of the melting of 1,3,3-trinitroazetidine

Physical properties of condensed-phase 1,3,3-trinitroazetidine (TNAZ) have been computed with molecular dynamics (MD) and a nonreactive, fully flexible force field formulated by combining the intramolecular interactions obtained from the Generalized AMBER Force Field and the rigid-molecule force field developed by Sorescu-Rice-Thompson [J. Phys. Chem. B 1997, 101, 798] (AMBER-SRT). The results are compared with MD calculations, using the AMBER force field. The predicted densities of crystalline TNAZ from both force fields are about 10% lower than the experimental value. The calculated thermodynamic melting point at 1 atm from the AMBER-SRT force field is 390 K, in good agreement with the measured value of 374 K, while the AMBER force field predicts a thermodynamic melting point of 462 K. The lattice parameters and the molecular and crystal structures calculated with the AMBER-SRT force field are in excellent agreement with experiment. Simulations with the AMBER-SRT force field were also used to generate the isotherm of TNAZ up to 4 GPa and the bulk modulus and its pressure derivative.