Tetrahedrality, hydrogen bonding and the density anomaly of the central force water model. A Monte Carlo study

V. Ravnik,B. Hribar-Lee,O. Pizio,M. Luksic
DOI: https://doi.org/10.5488/CMP.24.33503
2021-10-09
Abstract:Monte Carlo computer simulations in the canonical and grand canonical statistical ensemble were used to explore the properties of the central force (CF1) water model. The intramolecular structure of the H$_2$O molecule is well reproduced by the model. Emphasis was made on hydrogen bonding, and on the tehrahedral, $q$, and translational, $\tau$, order parameters. An energetic definition of the hydrogen bond gives more consistent results for the average number of hydrogen bonds compared to the one-parameter distance criterion. At 300 K, an average value of 3.8 was obtained. The $q$ and $\tau$ metrics were used to elucidate the water-like anomalous behaviour of the CF1 model. The structural anomalies lead to the density anomaly, with a good agreement of the model's density with the experimental $\rho(T)$ trends. The chemical potential-density projection of the model's equation of state was explored. Vapour-liquid coexistence was observed at sufficiently low temperatures.
Soft Condensed Matter
What problem does this paper attempt to address?