A molecular dynamics simulation study of crystalline 1,3,5-triamino-2,4,6-trinitrobenzene as a function of pressure and temperature

Quantum chemistry-based dipole polarizable and nonpolarizable force fields have been developed for 1,3,5-triamino-2,4,6-trinitrobenzene (TATB). Molecular dynamics simulations of TATB crystals were performed for hydrostatic pressures up to 10 GPa at 300 K and for temperatures between 200 and 400 K at atmospheric pressure. The predicted heat of sublimation and room-temperature volumetric hydrostatic compression curve were found to be in good agreement with available experimental data. The hydrostatic compression curves for individual unit cell parameters were found to be in reasonable agreement with those data. The pressure- and temperature-dependent second-order isothermal elastic tensor was determined for temperatures between 200 and 400 K at normal pressure and for pressures up to 10 GPa on the 300 K isotherm. Simulations indicate considerable anisotropy in the mechanical response, with modest softening and significant stiffening of the crystal with increased temperature and pressure, respectively. For most properties the polarizable potential was found to yield better agreement with available experimental properties.