Abstract:The accurate representation of the structural and dynamical properties of water is essential for simulating the unique behavior of this ubiquitous solvent. Here we assess the current status of describing liquid water using ab initio molecular dynamics, with a special focus on the performance of all the later generation Minnesota functionals. Findings are contextualized within the current knowledge on DFT for describing bulk water under ambient conditions and compared to experimental data. We find that, contrary to the prevalent idea that local and semilocal functionals overstructure water and underestimate dynamical properties, M06-L, revM06-L, and M11-L understructure water, while MN12-L and MN15-L overdistance water molecules due to weak cohesive effects. This can be attributed to a weakening of the hydrogen bond network, which leads to dynamical fingerprints that are over fast. While most of the hybrid Minnesota functionals (M06, M08-HX, M08-SO, M11, MN12-SX, and MN15) also yield understructured water, their dynamical properties generally improve over their semilocal counterparts. It emerges that exact exchange is a crucial component for accurately describing hydrogen bonds, which ultimately leads to corrections in both the dynamical and structural properties. However, an excessive amount of exact exchange strengthens hydrogen bonds and causes overstructuring and slow dynamics (M06-HF). As a compromise, M06-2X is the best performing Minnesota functional for water, and its D3 corrected variant shows very good structural agreement. From previous studies considering nuclear quantum effects (NQEs), the hybrid revPBE0-D3, and the rung-5 RPA (RPA@PBE) have been identified as the only two approximations that closely agree with experiments. Our results suggest that the M06-2X(-D3) functionals have the potential to further improve the reproduction of experimental properties when incorporating NQEs through path integral approaches. This work provides further proof that accurate modeling of water interactions requires the inclusion of both exact exchange and balanced (non-local) correlation, highlighting the need for higher rungs on Jacob's ladder to achieve predictive simulations of complex biological systems in aqueous environments.

Short-Range Structural Transformations in Water at High Pressures

Vibrational Features of Water at the Low-Density/High-Density Liquid Structural Transformations

High Pressure-Temperature Brillouin Study of Liquid Water: Evidence of the Structural Transition from Low-Density Water to High-Density Water

Water at Positive and Negative Pressures

Dynamics of Water Clusters Confined in Ionic Liquid at an Elevated Pressure

Water at negative pressure: Nuclear quantum effects

In Situ Brillouin Scattering Study of Water in High Pressure and High Temperature Conditions

X-ray Diffraction and Molecular Dynamics Study of Medium-range Order in Ambient and Hot Water

Effect of pressure on the anomalous response functions of a confined water monolayer at low temperature

Characterization of the Hydrogen-Bond Network in High-Pressure Water by Deep Potential Molecular Dynamics

Unraveling thermodynamic anomalies of water: A molecular simulation approach to probe the two-state theory with atomistic and coarse-grained water models

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

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

Path-integral simulation of ice VII: Pressure and temperature effects

Structural and dynamic properties of SPC/E water

Evolution of the structure of amorphous ice - from low-density amorphous (LDA) through high-density amorphous (HDA) to very high-density amorphous (VHDA) ice

Ephemeral ice-like local environments in classical rigid models of liquid water

Effect of pressure on the phase behavior and structure of water confined between nanoscale hydrophobic and hydrophilic plates

Structure and dynamics of liquid water from ab initio simulations: Adding Minnesota density functionals to Jacob's ladder

Experimental study of the polyamorphism of water. I. The isobaric transitions from amorphous ices to LDA at 4 MPa