Quick Search:

uzh logo
Browse by:

Zurich Open Repository and Archive

Maintenance: Tuesday, July the 26th 2016, 07:00-10:00

ZORA's new graphical user interface will be relaunched (For further infos watch out slideshow ZORA: Neues Look & Feel). There will be short interrupts on ZORA Service between 07:00am and 10:00 am. Please be patient.

Permanent URL to this publication: http://dx.doi.org/10.5167/uzh-3169

McGrath, M J; Siepmann, J I; Kuo, I F W; Mundy, C J; VandeVondele, J; Hutter, J; Mohamed, F; Krack, M (2005). Isobaric-isothermal Monte Carlo simulations from first principles: Application to liquid water at ambient conditions. ChemPhysChem, 6(9):1894-1901.

[img] PDF - Registered users only
View at publisher


A series of first-principles Monte Carlo simulations in the isobaric-isothermal ensemble were carried out for liquid water at ambient conditions (T = 298 K and p = 1 atm). The Becke-Lee-Yang-Parr (BLYP) exchange and correlation energy functionals and norm-conserving Goedecker-Teter-Hutter (GTH) pseudopotentials were employed with the CP2K simulation package to examine systems consisting of 64 water molecules. The fluctuations in the system volume encountered in simulations in the isobaric-isothermal ensemble require a reconsideration of the suitability of the typical charge-density cutoff and the regular grid-generation method previously used for the computation of the electrostatic energy in first-principles simulations in the microcanonical or canonical ensembles. In particular, it is noted that a much higher cutoff is needed and that the most computationally efficient method of creating grids con result in poor simulations. Analysis of the simulation trajectories using a very large charge-density cutoff at 1200 Ry and four different grid-generation methods point to a significantly underestimated liquid density of about 0.8 g cm(-3) resulting in a somewhat understructured liquid (with a value of about 2.7 for the height of the first peak in the oxygen-oxygen radial distribution function) for BLYP-GTH water at ambient conditions. In addition, a simulation using a charge-density cutoff at 280 Ry yields a higher density of 0.9 g cm(-3) showing the sensitivity of the simulation outcome to this parameter.


61 citations in Web of Science®
60 citations in Scopus®
Google Scholar™



1 download since deposited on 26 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
Deposited On:26 Mar 2009 09:34
Last Modified:05 Apr 2016 12:26
Publisher:Wiley VCH
Publisher DOI:10.1002/cphc.200400580

Users (please log in): suggest update or correction for this item

Repository Staff Only: item control page