Hydrogen bonds in concreto and in computro: the sequel

Abstract

In the framework of our comparative research concerning hydrogen bonding in the crystalline and liquid phases we have carried out molecular dynamics (MD) simulations of liquid methanol. Six different rigid three site models are compared. Five of them had been reported in the literature and one (OM2) we developed by a fit to the experimental molar volume, heat of vaporization and neutron weighted radial distribution function. In general the agreement with experiment is satisfactory for the different models. None of the models has an explicit hydrogen bond potential, but five of the six models show a degree of hydrogen bonding comparable to experiments on liquid methanol. The analysis of the simulation hydrogen bonds indicates that there is a distinct preference of the O⋯O axis to lie in the acceptor lone pairs plane, but hardly any for the lone pair directions. Ab initio calculations and crystal structure statistics of OH⋯O hydrogen bonds agree with this observation. The O⋯O hydrogen bond length distributions are similar for most models. The crystal structures show a sharper O⋯O distribution. Explicit introduction of harmonic motion with a quite realistic root mean square amplitude of 0.08 Å to the thermally averaged crystal distribution results in a distribution comparable to OM2 although the maximum of the former is found at shorter distance. On the basis of the analysis of the static properties of all models we conclude that our OM2, Jorgenson's OPLS and Haughney, Ferrario and McDonald's HFM1 models are good candidates for simulations of liquid methanol under isothermal, isochoric conditions. Partly flexible and completely rigid OM2 are simulated at constant pressure and with fixed volume. The flexible simulations give essentially the same (correct) results under both conditions, which is not surprising because the flexible form was fitted under both conditions. Rigid OM2 has a similar potential energy but larger pressure in the isochoric case and larger energy and far larger volume in the isobaric case. Radial distribution functions and hydrogen bond geometries are very similar for all four cases. Only in the case of the osobaric rigid methanol does the volume expansion seem to be accompanied by a slight preference for tetrahedrality around the oxygen atom.