TY - JOUR
T1 - Liquid water from first principles
T2 - Investigation of different sampling approaches
AU - Kuo, I. Feng W
AU - Mundy, Christopher J.
AU - McGrath, Matthew J.
AU - Siepmann, J. Ilja
AU - VandeVondele, Joost
AU - Sprik, Michiel
AU - Hutter, Jürg
AU - Chen, Bin
AU - Klein, Michael L.
AU - Mohamed, Fawzi
AU - Krack, Matthias
AU - Parrinello, Michele
PY - 2004/8/26
Y1 - 2004/8/26
N2 - 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 Nose-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/cm 3 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.
AB - 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 Nose-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/cm 3 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.
UR - http://www.scopus.com/inward/record.url?scp=4444372137&partnerID=8YFLogxK
UR - http://www.scopus.com/inward/citedby.url?scp=4444372137&partnerID=8YFLogxK
U2 - 10.1021/jp047788i
DO - 10.1021/jp047788i
M3 - Article
AN - SCOPUS:4444372137
VL - 108
SP - 12990
EP - 12998
JO - Journal of Physical Chemistry B
JF - Journal of Physical Chemistry B
SN - 1520-6106
IS - 34
ER -