The Nosé-Poincaré Method for Constant Temperature Molecular Dynamics

Stephen D. Bond*, Benedict J. Leimkuhler, Brian B. Laird

*Corresponding author for this work

Research output: Contribution to journalArticlepeer-review

Abstract / Description of output

We present a new extended phase space method for constant temperature (canonical ensemble) molecular dynamics. Our starting point is the Hamiltonian introduced by Nosé to generate trajectories corresponding to configurations in the canonical ensemble. Using a Poincaré time-transformation, we construct a Hamiltonian system with the correct intrinsic timescale and show that it generates trajectories in the canonical ensemble. Our approach corrects a serious deficiency of the standard change of variables (Nosé-Hoover dynamics), which yields a time-reversible system but simultaneously destroys the Hamiltonian structure. A symplectic discretization method is presented for solving the Nosé-Poincaré equations. The method is explicit and preserves the time-reversal symmetry. In numerical experiments, it is shown that the new method exhibits enhanced stability when the temperature fluctuation is large. Extensions are presented for Nosé chains, holonomic constraints, and rigid bodies.

Original languageEnglish
Pages (from-to)114-134
Number of pages21
JournalJournal of Computational Physics
Volume151
Issue number1
DOIs
Publication statusPublished - 1 May 1999

Fingerprint

Dive into the research topics of 'The Nosé-Poincaré Method for Constant Temperature Molecular Dynamics'. Together they form a unique fingerprint.

Cite this