Transition from patchlike to clusterlike inhomogeneity arising from hydrogen bonding in water

J Chem Phys. 2011 Feb 7;134(5):054121. doi: 10.1063/1.3552950.

Abstract

Assembling of water molecules via hydrogen bonding has been studied by molecular dynamics simulations using flexible potential model. The relationship between the number of H-bonds per molecule, n(HB), the size of H-bonded nets, k, and the size of patches of four-bonded molecules, k(4), has been examined for several thermodynamic states of water ranging from ambient to supercritical conditions. Two kinds of structural inhomogeneity have been found: the patchlike associated with the mean n(HB)> 2.0 and the clusterlike for n(HB)< 1.9. In compressed water up to ~473 K patches coexist with less ordered nets, both constituting the gel-like H-bonded network. The size of patches steeply decreases with the increasing temperature and the decreasing density of water. The inhomogeneity resulting from the presence of patches disappears above 473 K. This feature is associated with the rapid increase in the fraction of unbound molecules and with the breakage of the gel-like network into a variety of H-bonded clusters leading to the clusterlike structural inhomogeneity. In contrast to the patchlike inhomogeneity an increase in temperature and a decrease in density make this kind of inhomogeneity more pronounced. A degree of connectivity of H-bonds has been characterized by a parameter P(g) defined as the total fraction of molecules belonging to the H-bonded clusters of size k ≥ 5. The simulation-derived values of P(g) agree well with the predictions of the random bond theory giving the explicit expression for P(g) as a function of the mean n(HB). Going from ambient to supercritical conditions, we have found that the patchlike inhomogeneity is connected with the very slight reduction in P(g), whereas the clusterlike inhomogeneity generates a steep linear decrease of P(g) with the decreasing mean n(HB). The self-diffusion coefficient calculated for the thermodynamic states of water showing the clusterlike inhomogeneity has occurred to be inversely proportional to the density. We have also found that the clusterlike inhomogeneity is associated with the linear correlation between P(g) and the macroscopic properties of water: the static dielectric constant, the viscosity, and the density. The provided relationships allow one to estimate the degree of connectivity of hydrogen bonds from the measured macroscopic quantities.

MeSH terms

  • Hydrogen Bonding
  • Molecular Dynamics Simulation
  • Water / chemistry*

Substances

  • Water