Density, structure and dynamics of water: the effect of Van der Waals interactions

TL;DR

This study investigates how incorporating van der Waals interactions in density functional theory improves the accuracy of simulating water's density and diffusion, revealing structural changes and the importance of exchange functional choice.

Contribution

It demonstrates that using the van der Waals density functional significantly enhances the simulation of water's properties compared to traditional GGA functionals.

Findings

vdW-DF improves water density and diffusivity predictions

Increases non-H-bonded interstitial sites in water structure

Weakens hydrogen bond network, affecting coordination shells

Abstract

It is known that ab initio molecular dynamics (AIMD) simulations of liquid water, based on the generalized gradient approximation (GGA) to density functional theory (DFT), yield structural and diffusive properties in reasonable agreement with experiment only if artificially high temperatures are used in the simulations. The equilibrium density, at normal conditions, of DFT water has been recently shown by Schmidt et al. [J. Phys. chem. B, 113, 11959 (2009)] to be underestimated by different GGA functionals for exchange and correlation, and corrected by the addition of interatomic pair potentials to describe van derWaals (vdW) interactions. In this contribution we present a DFTAIMD study of liquid water using several GGA functionals as well as the van der Waals density functional (vdW-DF) of Dion et al. [Phys. Rev. Lett. 92, 246401(2004)]. As expected, we find that the density of water is grossly underestimated by GGA functionals. When a vdW-DF is used, the density improves drastically and the experimental diffusivity is reproduced without the need of thermal corrections. We analyze the origin of the density differences between all the functionals. We show that the vdW-DF increases the population of non-H-bonded interstitial sites, at distances between the first and second coordination shells. However, it excessively weakens the H-bond network, collapsing the second coordination shell. This structural problem is partially associated to the choice of GGA exchange in the vdW-DF. We show that a different choice for the exchange functional is enough to achieve an overall improvement both in structure and diffusivity.