Car-Parrinello and path integral molecular dynamics study of the intramolecular hydrogen bonds in the crystals of benzoylacetone and dideuterobenzoylacetone

The dynamics of the intramolecular short hydrogen bond in the molecular crystal of benzoylacetone and its deuterated analogue are investigated using ab initio molecular dynamics simulations. A study on intramolecular hydrogen bonding in 1-phenyl-1,3-butadione (I) and 1-deuteroxy-2-deutero-1-phenylbut-1-en-3-one (II) crystals has been carried out at 160 K and 300 K on the CPMD method level and at 300 K on the PIMD method level. The analysis of the two-dimensional free-energy landscape of reaction coordinate delta-parameter and R-O center dot center dot center dot O distances shows that the hydrogen (deuter) between the two oxygen atoms adopts a slightly asymmetrical position in the single potential well. When the nuclear quantum effects are taken into account, very large delocalization of the bridging proton is observed. These studies indicate that hydrogen bonds in the crystal of benzoylacetone have characteristic properties for the type of bonding model resonance-assisted hydrogen bonds (RAHB) without existing the equilibrium of the two tautomers. The infrared spectrum has been calculated, and a comparative vibrational analysis has been performed. The CPMD vibrational results appear to qualitatively agree with the experimental ones.