Nanoconfined water plays a crucial role in nanofluidics, biology, and cutting-edge technologies. The process of melting water monolayers and quasi-two-dimensional confined water involves, as an intermediate stage, the hexatic phase—a state that lies between solid and liquid and is characterized by quasi-long-range orientational order and short-range translational order. However, the influence of hydrogen bond (HB) cooperativity in this process has not been thoroughly investigated. This gap hampers our understanding of the phase behavior of confined water and limits the accuracy of our models. To address this, we extend the water model developed by Franzese and Stanley, which explicitly includes many-body interactions (MBIs) of HBs. We distinguish the contributions of three-body and five-body HB-MBIs. Our Monte Carlo calculations in the isobaric–isothermal ensemble produce a detailed pressure–temperature phase diagram, revealing polymorphism and polyhexaticity: low-density square ice and high-density triangular ice are separated from the liquid phase by distinct hexatic phases. Three-body interactions notably promote crystallization and can destabilize the low-density hexatic phase, while cooperative five-body interactions help restore it, thus modifying the thermodynamic landscape. These findings demonstrate that HB-MBIs are key to determining the phase behavior of confined water, influencing phenomena such as non-monotonic specific heat, maximum-density lines, and the accessibility of the liquid–liquid critical point. Beyond advancing theoretical understanding, these results have wide-ranging implications for nanofluidics, interfacial science, and applications in biology, food technology, and pharmaceutics, where controlling water under confinement is essential.
Vilanova et al. (Mon,) studied this question.