Quantum Chemistry–Based Force Field for Simulations of Energetic Dinitro Compounds

ABSTRACT A quantum chemistry–based force field for molecular dynamics simulations of energetic dinitro compounds has been developed, based on intermolecular binding energies, molecular geometries, molecular electrostatic potentials, and conformational energies obtained from quantum chemistry calculations on model compounds. Nonbonded parameters were determined by fitting experimental densities and heats of vaporizations of model compounds. Torsional parameters were parameterized to reproduce accurately the relative conformational energy minima and barriers in 2,2-dinitropropane, di-methoxy di-methyl ether, 2,2-dinitro-3-methoxypropane, and bis(2,2-dinitropropyl)formal. Molecular dynamics simulations using the developed force field accurately reproduce thermodynamic and transport properties of 1,1-dinitroethane, 2,2-dinitropropane, and a eutectic mixture of bis(2,2-dinitropropyl)formal and bis(2,2-dinitropropyl)acetal.