Hydrogen bond networks in water and methanol with varying interaction strengths

Metropolis Monte Carlo simulations of hydrogen-bonded liquids (water and methanol) were performed with the well tested effective pair potentials TIP5P and OPLS. The Coulomb contribution for the interaction potential was damped by a factor η varied from 1 to 0.49 for water and 1 to 0.15 for methanol. As a result, the networks formed by the hydrogen-bonded molecules presented interesting properties as a function of η, including small-world patterns and percolation transitions. These complex networks were analyzed by local (clustering coefficients, average degrees), semi-global (path lengths) and global (spectral densities) properties, and islands statistics. From these properties, small-world behavior was found for η in the range 0.60-0.75 for both liquids, interestingly independent of the molecular structure of the liquid. Phase transition behavior was observed for the average degrees and the clustering coefficient curves with critical values at 0.55 for water and 0.34 for methanol. Macroscopic properties such as mass density and vaporization enthalpy were also parametrically dependent on η and they presented phase transition behavior that coincides with the critical values obtained from the topological analysis. This is probably the first time that such phase transitions are observed for these quantities and shows a direct relation between macroscopic properties and topological features of hydrogen bond networks.