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.

Classical density-functional theory for water

Improved association in a classical density functional theory for water

A molecular density functional theory to study solvation in water

A classical density functional theory for solvation across length scales

Hydration of Clays at the Molecular Scale: The Promising Perspective of Classical Density Functional Theory

Mass Density Fluctuations in Quantum and Classical descriptions of Liquid Water

Classical density functional theory to tackle solvation in molecular liquids

Consistent Density Functional Theory-Based Description of Ion Hydration Through Density-Corrected Many-Body Representations

The electronic structure of liquid water within density functional theory

First-principles energetics of water: a many-body analysis

A molecular density functional theory for associating fluids in 3D geometries

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

Solvation of complex surfaces via molecular density functional theory

Beyond the mean-field approximation for pair correlations in classical density functional theory: Reference inhomogeneous non-associating monomeric fluids for use with SAFT-VR Mie DFT

A doubly associated reference perturbation theory for water

Density functional theory (DFT) study of water autoionization in solvated clusters

Density Functionals with Quantum Chemical Accuracy: From Machine Learning to Molecular Dynamics

Density Functional Theory -- an introduction

Fine Tuning the Intermolecular Interactions of Water Clusters Using the Dispersion-Corrected Density Functional Theory

Optimization of an exchange-correlation density functional for water

Nonlocal density functional theory of water taking into account many-body dipole correlations: Binodal and surface tension of 'liquid-vapour' interface