Anomalies and Local Structure of Liquid Water from Boiling to the Supercooled Regime as Predicted by the Many-Body MB-pol Model

Thomas E. Gartner III, Kelly M. Hunter, Eleftherios Lambros, Alessandro Caruso, Marc Riera, Gregory R. Medders, Athanassios Z. Panagiotopoulos, Pablo G. Debenedetti and Francesco Paesani, J. Phys. Chem. Lett. 2022, 13, XXX, 3652–3658 (2022)

For the last 50 years, researchers have sought molecular models that can accurately reproduce water’s microscopic structure and thermophysical properties across broad ranges of its complex phase diagram. Herein, molecular dynamics simulations with the many-body MB-pol model are performed to monitor the thermodynamic response functions and local structure of liquid water from the boiling point down to deeply supercooled temperatures at ambient pressure. The isothermal compressibility and isobaric heat capacity show maxima near 223 K, in excellent agreement with recent experiments, and the liquid density exhibits a minimum at ~208 K. A local tetrahedral arrangement, where each water molecule accepts and donates two hydrogen bonds, is found to be the most probable hydrogen-bonding topology at all temperatures. This work suggests that MB-pol may provide predictive capability for studies of liquid water’s physical properties across broad ranges of thermodynamic states, including the so-called water’s “no man’s land” which is difficult to probe experimentally.