Constrained molecular dynamics in the isothermal-isobaric ensemble and its adaptation for adiabatic free energy dynamics

T. Q. Yu, M. E. Tuckerman

Research output: Contribution to journalReview articlepeer-review

Abstract

The implementation of holonomic constraints within measure-preserving integrators for molecular dynamics simulations in the isothermal-isobaric ensemble is considered. We review the basic methodology of generating measure-preserving integrators for the microcanonical, canonical, and isothermal-isobaric ensembles and proceed to show how the standard SHAKE and RATTLE algorithms must be modified for the isothermal-isobaric ensemble. Comparison is made between constrained and unconstrained simulations employing multiple time scale integration techniques. Finally, we describe a temperature accelerated version of the isothermal-isobaric molecular dynamics approach, in which the cell matrix is adiabatically decoupled from the particles and maintained at a high temperature as a means of exploring polymorphism in molecular crystals. We demonstrate that constraints can be easily adapted for this new approach and, again, we compare the performace of this temperature-accelerated scheme with and without bond constraints.

Original languageEnglish (US)
Pages (from-to)183-209
Number of pages27
JournalEuropean Physical Journal: Special Topics
Volume200
Issue number1
DOIs
StatePublished - Nov 2011

ASJC Scopus subject areas

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

Fingerprint Dive into the research topics of 'Constrained molecular dynamics in the isothermal-isobaric ensemble and its adaptation for adiabatic free energy dynamics'. Together they form a unique fingerprint.

Cite this