Search papers, labs, and topics across Lattice.
The study investigates the reproducibility of ab initio molecular dynamics simulations of liquid water using the revPBE-D3 density functional, revealing significant discrepancies in previously published results for diffusion coefficient and density. By employing long-range machine-learning interatomic potentials trained on carefully converged DFT data, the authors achieve consensus across six independent codes, resolving these discrepancies. The corrected simulations demonstrate that prior studies overestimated density and underestimated the diffusion coefficient due to limitations in basis sets, pseudopotentials, and statistical sampling.
Widely used simulations of liquid water are surprisingly inconsistent, with errors exceeding 20% in key properties, but this work provides a robust, reproducible benchmark.
Liquid water is fundamentally important, and its accurate computer simulation has been the driving force for myriad methodological developments. Ab initio molecular dynamics with forces obtained from density functional theory (DFT) is now a standard tool widely used by researchers. However, we reveal that previous studies of liquid water using the same widely-used density functional (revPBE-D3) exhibit significant discrepancies with one another, varying by over 20% in the diffusion coefficient and 10% in the density, raising fundamental questions about reproducibility. By combining modern long-range machine-learning interatomic potentials that enable robust statistical sampling with carefully converged DFT training data, we resolve these discrepancies, achieving consensus across six diverse community codes. Our predictions differ markedly from previous literature: we show that most previous results overestimate the density and underestimate the diffusion coefficient of revPBE-D3 water due to basis set incompleteness and pseudopotential inconsistencies, coupled with limitations in statistical sampling (in some cases). These benchmark values provide a reliable reference for validating current and future implementations of DFT-based ab initio molecular dynamics. Reaching agreement establishes confidence and credibility and serves as a prerequisite for the systematic assessment of new density functionals and numerical approximations.
