Header

UZH-Logo

Maintenance Infos

Liquid water from first principles: investigation of different sampling approaches


Kuo, I F W; Mundy, C J; McGrath, M J; Siepmann, I J; VandeVondele, J; Sprik, M; Hutter, J; Chen, B; Klein, M L; Mohamed, F; Krack, M; Parrinello, M (2004). Liquid water from first principles: investigation of different sampling approaches. Journal of Physical Chemistry. B, 108(34):12990-12998.

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 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 Angstrom(2)/ps.}

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 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 Angstrom(2)/ps.}

Statistics

Citations

Dimensions.ai Metrics
320 citations in Web of Science®
326 citations in Scopus®
Google Scholar™

Altmetrics

Downloads

2 downloads since deposited on 27 Mar 2009
0 downloads since 12 months
Detailed statistics

Additional indexing

Item Type:Journal Article, refereed, original work
Communities & Collections:07 Faculty of Science > Department of Chemistry
Dewey Decimal Classification:540 Chemistry
Scopus Subject Areas:Physical Sciences > Physical and Theoretical Chemistry
Physical Sciences > Surfaces, Coatings and Films
Physical Sciences > Materials Chemistry
Language:English
Date:2004
Deposited On:27 Mar 2009 08:44
Last Modified:01 Dec 2023 02:43
Publisher:American Chemical Society
ISSN:1520-5207
OA Status:Closed
Publisher DOI:https://doi.org/10.1021/jp047788i