On the pressure calculation for polarizable models in computer simulation

Péter T. Kiss, Andrs Baranyai

Research output: Contribution to journalArticle

8 Citations (Scopus)

Abstract

We present a short overview of pressure calculation in molecular dynamics or Monte Carlo simulations. The emphasis is given to polarizable models in order to resolve the controversy caused by the paper of M. J. Louwerse and E. J. Baerends Chem. Phys. Lett. 421, 138 (2006)10.1016/j.cplett.2006.01.087 about pressure calculation in systems with periodic boundaries. We systematically derive expressions for the pressure and show that despite the lack of explicit pairwise additivity, the pressure formula for polarizable models is identical with that of nonpolarizable ones. However, a strict condition for using this formula is that the induced dipole should be in perfect mechanical equilibrium prior to pressure calculation. The perfect convergence of induced dipoles ensures conservation of energy as well. We demonstrate using more cumbersome but exact methods that the derived expressions for the polarizable model of water provide correct numerical results. We also show that the inaccuracy caused by imperfect convergence of the induced dipoles correlates with the inaccuracy of the calculated pressure.

Original languageEnglish
Article number104109
JournalJournal of Chemical Physics
Volume136
Issue number10
DOIs
Publication statusPublished - Mar 14 2012

    Fingerprint

ASJC Scopus subject areas

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

Cite this