A practical method to avoid zero-point leak in molecular dynamics calculations: Application to the water dimer

G. Czakó, Alexey L. Kaledin, Joel M. Bowman

Research output: Contribution to journalArticle

58 Citations (Scopus)

Abstract

We report the implementation of a previously suggested method to constrain a molecular system to have mode-specific vibrational energy greater than or equal to the zero-point energy in quasiclassical trajectory calculations [J. M. Bowman, J. Chem. Phys. 91, 2859 (1989); W. H. Miller, J. Chem. Phys. 91, 2863 (1989)]. The implementation is made practical by using a technique described recently [G. Czakó and J. M. Bowman, J. Chem. Phys. 131, 244302 (2009)], where a normal-mode analysis is performed during the course of a trajectory and which gives only real-valued frequencies. The method is applied to the water dimer, where its effectiveness is shown by computing mode energies as a function of integration time. Radial distribution functions are also calculated using constrained quasiclassical and standard classical molecular dynamics at low temperature and at 300 K and compared to rigorous quantum path integral calculations.

Original languageEnglish
Article number164103
JournalThe Journal of Chemical Physics
Volume132
Issue number16
DOIs
Publication statusPublished - Apr 28 2010

Fingerprint

Dimers
Molecular dynamics
dimers
Trajectories
molecular dynamics
Water
trajectories
water
Distribution functions
zero point energy
radial distribution
distribution functions
energy
Temperature

ASJC Scopus subject areas

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

Cite this

A practical method to avoid zero-point leak in molecular dynamics calculations : Application to the water dimer. / Czakó, G.; Kaledin, Alexey L.; Bowman, Joel M.

In: The Journal of Chemical Physics, Vol. 132, No. 16, 164103, 28.04.2010.

Research output: Contribution to journalArticle

@article{af12efd0d9d847e091791a39bd545f3e,
title = "A practical method to avoid zero-point leak in molecular dynamics calculations: Application to the water dimer",
abstract = "We report the implementation of a previously suggested method to constrain a molecular system to have mode-specific vibrational energy greater than or equal to the zero-point energy in quasiclassical trajectory calculations [J. M. Bowman, J. Chem. Phys. 91, 2859 (1989); W. H. Miller, J. Chem. Phys. 91, 2863 (1989)]. The implementation is made practical by using a technique described recently [G. Czak{\'o} and J. M. Bowman, J. Chem. Phys. 131, 244302 (2009)], where a normal-mode analysis is performed during the course of a trajectory and which gives only real-valued frequencies. The method is applied to the water dimer, where its effectiveness is shown by computing mode energies as a function of integration time. Radial distribution functions are also calculated using constrained quasiclassical and standard classical molecular dynamics at low temperature and at 300 K and compared to rigorous quantum path integral calculations.",
author = "G. Czak{\'o} and Kaledin, {Alexey L.} and Bowman, {Joel M.}",
year = "2010",
month = "4",
day = "28",
doi = "10.1063/1.3417999",
language = "English",
volume = "132",
journal = "Journal of Chemical Physics",
issn = "0021-9606",
publisher = "American Institute of Physics Publising LLC",
number = "16",

}

TY - JOUR

T1 - A practical method to avoid zero-point leak in molecular dynamics calculations

T2 - Application to the water dimer

AU - Czakó, G.

AU - Kaledin, Alexey L.

AU - Bowman, Joel M.

PY - 2010/4/28

Y1 - 2010/4/28

N2 - We report the implementation of a previously suggested method to constrain a molecular system to have mode-specific vibrational energy greater than or equal to the zero-point energy in quasiclassical trajectory calculations [J. M. Bowman, J. Chem. Phys. 91, 2859 (1989); W. H. Miller, J. Chem. Phys. 91, 2863 (1989)]. The implementation is made practical by using a technique described recently [G. Czakó and J. M. Bowman, J. Chem. Phys. 131, 244302 (2009)], where a normal-mode analysis is performed during the course of a trajectory and which gives only real-valued frequencies. The method is applied to the water dimer, where its effectiveness is shown by computing mode energies as a function of integration time. Radial distribution functions are also calculated using constrained quasiclassical and standard classical molecular dynamics at low temperature and at 300 K and compared to rigorous quantum path integral calculations.

AB - We report the implementation of a previously suggested method to constrain a molecular system to have mode-specific vibrational energy greater than or equal to the zero-point energy in quasiclassical trajectory calculations [J. M. Bowman, J. Chem. Phys. 91, 2859 (1989); W. H. Miller, J. Chem. Phys. 91, 2863 (1989)]. The implementation is made practical by using a technique described recently [G. Czakó and J. M. Bowman, J. Chem. Phys. 131, 244302 (2009)], where a normal-mode analysis is performed during the course of a trajectory and which gives only real-valued frequencies. The method is applied to the water dimer, where its effectiveness is shown by computing mode energies as a function of integration time. Radial distribution functions are also calculated using constrained quasiclassical and standard classical molecular dynamics at low temperature and at 300 K and compared to rigorous quantum path integral calculations.

UR - http://www.scopus.com/inward/record.url?scp=77952412703&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=77952412703&partnerID=8YFLogxK

U2 - 10.1063/1.3417999

DO - 10.1063/1.3417999

M3 - Article

C2 - 20441254

AN - SCOPUS:77952412703

VL - 132

JO - Journal of Chemical Physics

JF - Journal of Chemical Physics

SN - 0021-9606

IS - 16

M1 - 164103

ER -