More converged, less accurate? Reassessing standard choices for ab initio water using machine learning potentials

Accurately simulating the properties of liquid water remains a central challenge in molecular simulations. In this work, we use machine learning potentials to investigate how the convergence settings of electronic structure calculations impact the predicted structural and dynamical properties of simulated water and ice. We evaluate the true performance of several reference methods in classical and path-integral molecular dynamics. When we compare a popular, computationally pragmatic revPBE0-D3 setup against a highly converged one, our results reveal that its widely reported experimental agreement degrades. Applying the same highly converged settings to the $\mathrm{\omega}$B97X-rV functional, we find an improved agreement with experimental results. MP2 with a triple-$\zeta$ basis set commonly used for liquid water shows poor performance, which is indicative of insufficient convergence. These findings underscore the need for fully converged reference calculations when evaluating the fundamental accuracy of electronic structure methods and developing reliable models for aqueous systems.