The integrated DL_POLY/DL_FIELD/DL_ANALYSER software platform for molecular dynamics simulations for exploration of the synthonic interactions in saturated benzoic acid/hexane solutions

Three separately developed software Molecular Dynamics packages at Daresbury Laboratory, namely DL_FIELD (DL_F), DL_POLY and DL_ANALYSER, have been integrated to form an efficient computational infrastructure to investigate the detailed solution chemistry of saturated benzoic acid in hexane solutions. These software capabilities are demonstrated, in combination with the Synthonic Engineering tools and density functional theory (DFT) calculations, to assess the extent that the solute-solute intermolecular synthonic interactions in solution mirrors the synthons in the crystal structure. The results show that the majority of the COOH groups are forming OH … O H-bonds, which are a combination of classic OH … O homo-dimers and three membered H-bonding clusters. The formation of pi-pi stacking interactions is observed, but in far fewer numbers than observed for the OH … O interactions. The DFT simulations of the IR spectra of the multiple benzoic acid aggregates extracted from the MD trajectories provides further in-depth analysis of previously published IR data, by matching simulated peaks to the experimental peaks, hence identifying the exact bonding modes that are responsible for such peaks. This study demonstrates the value of a multi-scale and multi-technique approach to exploring the molecular transition pathway from solution to crystal structure.