Web Release Date: August 3,
Liquid Water from First Principles: Investigation of Different Sampling Approaches
and
Chemistry and Materials Science Directorate, Lawrence Livermore National Laboratory, Livermore, California 94550, Departments of Chemistry and of Chemical Engineering and Material Science, University of Minnesota, 207 Pleasant Street SE, Minneapolis, Minnesota 55455-0431, Department of Chemistry, University of Cambridge, Lensfield Road, Cambridge CB2 1EW, United Kingdom, Physical Chemistry Institute, University of Zurich, Winterthurerstrasse 190, CH-8057 Zurich, Switzerland, Department of Chemistry, Louisiana State University, Baton Rouge, Louisiana 70803-1804, Center for Molecular Modeling and Department of Chemistry, University of Pennsylvania, 231 South 34th Street, Philadelphia, Pennsylvania 19104-6323, and Laboratory of Physical Chemistry, ETH Zurich, USI-Campus, Via Giuseppe Buffi 13, CH-6904 Lugano, Switzerland
Received: May 21, 2004
In Final Form: July 6, 2004
Abstract:
A series of first principles molecular dynamics and Monte Carlo simulations were carried out for liquid water to investigate the reproducibility of different sampling approaches. These simulations include Car-Parrinello molecular dynamics simulations using the program CPMD with different values of the fictitious electron mass in the microcanonical and canonical ensembles, Born-Oppenheimer molecular dynamics using the programs CPMD and CP2K in the microcanonical ensemble, and Metropolis Monte Carlo using CP2K in the canonical ensemble. With the exception of one simulation for 128 water molecules, all other simulations were carried out for systems consisting of 64 molecules. Although the simulations yield somewhat fortuitous agreement in structural properties, analysis of other properties demonstrate that one should exercise caution when assuming the reproducibility of Car-Parrinello and Born-Oppenheimer molecular dynamics simulations for small system sizes in the microcanonical ensemble. In contrast, the molecular dynamics and Monte Carlo simulations in the canonical ensemble appear to be more reliable. Furthermore, in the case of canonical Car-Parrinello molecular dynamics simulations the application of Nosé-Hoover chain thermostats allows the use of larger fictitious electron masses. For the Becke-Lee-Yang-Parr exchange and correlation energy functionals and norm-conserving Troullier-Martins or Goedecker-Teter-Hutter pseudopotentials, these simulations at a fixed density of 1.0 g/cm3 and a temperature close to 315 K point to an overstructured liquid with a height of the first peak in the oxygen-oxygen radial distribution function of about 3.0, an underestimated value of the classical constant-volume heat capacity of about 70 J/(mol K), and an underestimated self-diffusion constant of about 0.04 Å2/ps.
Download the full text: PDF | HTML
