Search

The accurate representation of multidimensional potential energy surfaces is a necessary requirement for realistic computer simulations of molecular systems. The continued increase in computer power accompanied by advances in correlated electronic structure methods nowadays enables routine calculations of accurate interaction energies for small systems, which can then be used as references for the development of analytical potential energy functions (PEFs) rigorously derived from many-body (MB) expansions. Building on the accuracy of the MB-pol many-body PEF, we investigate here the performance of permutationally invariant polynomials (PIPs), neural networks, and Gaussian approximation potentials (GAPs) in representing water two-body and three-body interaction energies, denoting the resulting potentials PIP-MB-pol, Behler-Parrinello neural network-MB-pol, and GAP-MB-pol, respectively. Our analysis shows that all three analytical representations exhibit similar levels of accuracy in reproducing both two-body and three-body reference data as well as interaction energies of small water clusters obtained from calculations carried out at the coupled cluster level of theory, the current gold standard for chemical accuracy. These results demonstrate the synergy between interatomic potentials formulated in terms of a many-body expansion, such as MB-pol, that are physically sound and transferable, and machine-learning techniques that provide a flexible framework to approximate the short-range interaction energy terms. [ABSTRACT FROM AUTHOR]

We carried out extensive calculations of liquid water at different temperatures and pressures using the BK3 model suggested recently [P. T. Kiss and A. Baranyai, J. Chem. Phys. 138, 204507 (2013)]. In particular, we were interested in undercooled regions to observe the propensity of water to form tetrahedral coordination of closest neighbors around a central molecule. We compared the found tetrahedral order with the number of hydrogen bonds and with the partial pair correlation functions unfolded as distributions of the closest, the second closest, etc. neighbors. We found that contrary to the number of hydrogen bonds, tetrahedrality changes substantially with state variables. Not only the number of tetrahedral arrangements increases with lowering the pressure, the density, and the temperature but the domain size of connecting tetrahedral structures as well. The difference in tetrahedrality is very pronounced between the two sides of the Widom line and even more so between the low density amorphous (LDA) and high density amorphous (HDA) phases. We observed that in liquid water and in HDA, the 5th water molecule, contrary to ice and LDA, is positioned between the first and the second coordination shell. We found no convincing evidence of structural heterogeneity or regions referring to structural transition. [ABSTRACT FROM AUTHOR]