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.

Hamiltonian Adaptive Resolution Molecular Dynamics Simulation of Infrared Dielectric Functions of Liquids

Path Integral Molecular Dynamics within the Grand Canonical-like Adaptive Resolution Technique: Simulation of Liquid Water

Open Boundary Simulations of Proteins and Their Hydration Shells by Hamiltonian Adaptive Resolution Scheme

Modeling diffusive dynamics in adaptive resolution simulation of liquid water

Molecular dynamics simulation of dielectric properties of water

Adaptive Molecular Resolution Approach in Hamiltonian Form: An Asymptotic Analysis

Investigation of water-mediated intermolecular interactions with the adaptive resolution simulation technique

Molecular Dynamics Simulation for Dmf Aqueous Solution

ADCHα-I population analysis and constrained dipole moment density functional theory in force fields for molecular simulations

A new method to calculate broadband dielectric spectra of solvents from molecular dynamics simulations demonstrated with polarizable force fields

Liquid water simulation using hydrogen bond corrected SCAN and neural network potentials.

Molecular Simulation Of Lipids And Water: Atomistic, Coarse-Grained, And Mixed Resolutions

Isotope Effects in Liquid Water via Deep Potential Molecular Dynamics

An Efficient Method for Accurate Evaluation of the Site–site Direct Correlation Functions of Molecular Fluids

Ab Initio Simulation of Liquid Water without Artificial High Temperature

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

Application of Molecular Dynamics Simulations in Molecular Property Prediction II: Diffusion Coefficient.

Modeling Liquid Water by Climbing Up Jacob's Ladder in Density Functional Theory Facilitated by Using Deep Neural Network Potentials

Simple and Accurate One-Body Energy and Dipole Moment Surfaces for Water and Beyond

Liquid Li Structure and Dynamics: A Comparison Between OFDFT and Second Nearest‐neighbor Embedded‐atom Method