Classical and quantum-mechanical studies of crystalline FOX-7 (1,1-diamino-2,2-dinitroethylene)

First principles molecular orbital and plane-wave ab initio calculations ha ve been used to investigate the structural and vibrational properties of th e highly efficient low sensitive explosive 1,1-diamino-2,2 dinitroethylene (FOX-7) in both the gas and solid phases. The ab initio molecular orbital c alculations performed at second-order (MP2) and fourth-order (MP4) Moller-P lesset levels and using density-functional theory (DFT) methods with B3LYP functional indicate that in the gas phase FOX-7 is the most stable conforme r relative to its cis-1,2 and trans-1,2 isomers. The calculated MP2 and DFT structures for the FOX-7 molecule agree well with the experimental X-ray c onfiguration but with twists of the nitro and amino groups much larger than in the solid phase. The calculated fundamental vibrational frequencies at the DFT level generally compare well. with the MP2 results. The IR spectra were computed for the three isomers. The structural properties of the FOX-7 crystal have been studied by a plane-wave DFT method. These calculations w ere done with periodic boundary conditions in all three directions. The opt imization of the crystal structure has been done with full relaxation of th e atomic positions and of the lattice parameters under P2(1)/n symmetry. Th e predicted crystal structure is in good agreement with X-ray data. We have developed an intermolecular potential to describe the structure of the FOX -7 crystal in the approximation of rigid molecules. This potential is compo sed of pairwise exp-6 Buckingham terms and Coulombic interactions. Crystal- packing calculations without symmetry constraints performed with the propos ed potential accurately reproduce the main crystallographic features and yi eld very good agreement with the estimated lattice energy. This intermolecu lar potential was further tested in isothermal-isobaric molecular dynamics simulations at atmospheric pressure and in the temperature range of 4.2-450 K. It is found that the increase of temperature does not significantly cha nge the orientations of the molecules inside the unit cell. The thermal exp ansion coefficients calculated for the model indicate anisotropic behavior with the largest expansion along the b crystallographic direction.