Free energy of liquid water on the basis of quasichemical theory and ab initio molecular dynamics

We use ab initio molecular dynamics as a basis for quasichemical theory evaluation of the free energy of water near conventional liquid thermodynamic states. The Perdew-Wang-91 (PW91), Perdew-Burke-Ernzerhof (PBE), and revised PBE (rPBE) functionals are employed. The oxygen radial density distribution using the rPBE functional is in reasonable agreement with current experiments, whereas the PW91 and PBE functionals predict a more structured oxygen radial density distribution. The diffusion coefficient with the rPBE functional is in reasonable accord with experiments. Using a maximum entropy procedure, we obtain [Formula presented] from the coordination number distribution [Formula presented] for oxygen atoms having n neighbors. Likewise, we obtain [Formula presented] from [Formula presented] the probability of observing cavities of specified radius containing n water molecules. The probability [Formula presented] is a measure of the local chemical interactions and is central to the quasichemical theory of solutions. The probability [Formula presented] central to the theory of liquids, is a measure of the free energy required to open cavities of defined sizes in the solvent. Using these values and a reasonable model for electrostatic and dispersion effects, the hydration free energy of water in water at 314 K is calculated to be [Formula presented] with the rPBE functional, in encouraging agreement with the experimental value of [Formula presented].