Development of All-Atom Empirical Potentials for Phloroglucinol (Phg) and Hexachlorocyclotriphosphazene (HCCP) Based on their Crystal Phase Structures

Two existing generic force fields have been augmented with partial charges and tuned in order to give intercompatible all-atom empirical potentials that can satisfactorily represent the known crystal phase structures of the organic phloroglucinol (Phg) (C₆H₆O₃) and inorganic hexachlorocyclotriphosphazene (HCCP) (N₃P₃Cl₆) molecules at several temperatures. It has been proposed that HCCP-Phg network polymers could act as efficient H₂ barrier layers in hydrogen storage tanks for cars. However, essential requirements for modeling such networks are adequate representations of both monomers in their pure dense phases using a common form of force field. Tests of their ability to maintain stable crystal structures have been made using classical molecular dynamics (MD) simulations on large 800-molecule supercells. The force fields have been optimized to match the densities calculated from the experimental unit cell dimensions at ambient conditions as well as the intermolecular potential energy, as estimated from experimental enthalpies of sublimation. For Phg, the crystal structure is stabilized by a network of hydrogen bonds and the Coulombic interactions contribute to over 55% of the total intermolecular potential energy. In contrast, the crystal structure of HCCP is intrinsically stabilized by the van der Waals terms. Both optimized force fields reproduce very well the orthorhombic symmetry of their respective crystals under constant-pressure <i>N</i><b><i>P</i></b><i>T</i> conditions. The model parameters tuned at ambient temperature also give reasonable agreement with crystallographic data at lower temperatures.