Using Car-Parrinello simulations and microscopic order descriptors to reveal two locally favored structures with distinct molecular dipole moments and dynamics in ambient liquid water

TL;DR

This study uses Car-Parrinello ab initio molecular dynamics to identify two distinct local water structures with different dipole moments and dynamics, revealing their coexistence under ambient conditions and implications for water's anomalies.

Contribution

It introduces a microscopic order parameter to distinguish two locally favored water structures and demonstrates their coexistence and properties under ambient conditions.

Findings

Two local water structures with distinct order and dipole moments identified.

High-{} molecules have a 6% higher dipole moment than Low-{}.

Low-{} molecules reorient faster and have less stable shells.

Abstract

Water is essential for life and technological applications, mainly for its unique thermodynamic and dynamic properties, often anomalous or counterintuitive. These anomalies result from the hydrogen-bonds fluctuations, as evidenced by studies for supercooled water. However, it is difficult to characterize these fluctuations under ambient conditions. Here, we fill this knowledge gap thanks to the Car-Parrinello ab initio molecular dynamics (MD) simulation technique. We calculate the local structural order parameter {\zeta}, quantifying the coordination shells separation, and find two locally-favored structures or states: High-{\zeta} and Low-{\zeta}. On average, High-{\zeta} molecules have a tetrahedral arrangement, with four hydrogen bonds, and the first and the second coordination shell well separated. The Low-{\zeta} molecules are less connected, partially merging the first and the second shells. The appearance of isosbestic points in the radial distribution functions and the collective density fluctuations at different length scales and timescales reveal that the two-state model, consistent with available experimental data for supercooled water, also holds under ambient conditions, as we confirm by analyzing the vibrational spectrum of both types of water molecules. Significant consequences of the structural differences between the two states are that High-{\zeta} molecules have a dipole moment 6 % higher than Low-{\zeta}. At the same time, Low-{\zeta} structures are more disordered and with more significant angular fluctuations. These differences are also reflected in the dynamics under ambient conditions. The Low-{\zeta} molecules decorrelate their reorientation faster than High-{\zeta} and merge their coordination shells within 0.2 ps, while the High-{\zeta} preserve the shell separation for longer times.