Calculation of heat capacities of light and heavy water by path-integral molecular dynamics

Motoyuki Shiga, Wataru Shinoda

Research output: Contribution to journalArticlepeer-review

11 Citations (Scopus)


As an application of atomistic simulation methods to heat capacities, path-integral molecular dynamics has been used to calculate the constant-volume heat capacities of light and heavy water in the gas, liquid, and solid phases. While the classical simulation based on conventional molecular dynamics has estimated the heat capacities too high, the quantum simulation based on path-integral molecular dynamics has given reasonable results based on the simple point-charge/flexible potential model. The calculated heat capacities (divided by the Boltzmann constant) in the quantum simulation are 3.1 in the vapor H2 O at 300 K, 6.9 in the liquid H2 O at 300 K, and 4.1 in the ice Ih H2 O at 250 K, respectively, which are comparable to the experimental data of 3.04, 8.9, and 4.1, respectively. The quantum simulation also reproduces the isotope effect. The heat capacity in the liquid D2 O has been calculated to be 10% higher than that of H2 O, while it is 13% higher in the experiment. The results demonstrate that the path-integral simulation is a promising approach to quantitatively evaluate the heat capacities for molecular systems, taking account of quantum-mechanical vibrations as well as strongly anharmonic motions.

Original languageEnglish
Article number134502
JournalJournal of Chemical Physics
Issue number13
Publication statusPublished - Oct 3 2005
Externally publishedYes

ASJC Scopus subject areas

  • Physics and Astronomy(all)
  • Physical and Theoretical Chemistry


Dive into the research topics of 'Calculation of heat capacities of light and heavy water by path-integral molecular dynamics'. Together they form a unique fingerprint.

Cite this