A comparison of generalized hybrid Monte Carlo methods with and without momentum flip

Elena Akhmatskaya, Nawaf Bou-Rabee, Sebastian Reich

Research output: Contribution to journalArticlepeer-review

35 Scopus citations

Abstract

The generalized hybrid Monte Carlo (GHMC) method combines Metropolis corrected constant energy simulations with a partial random refreshment step in the particle momenta. The standard detailed balance condition requires that momenta are negated upon rejection of a molecular dynamics proposal step. The implication is a trajectory reversal upon rejection, which is undesirable when interpreting GHMC as thermostated molecular dynamics. We show that a modified detailed balance condition can be used to implement GHMC without momentum flips. The same modification can be applied to the generalized shadow hybrid Monte Carlo (GSHMC) method. Numerical results indicate that GHMC/GSHMC implementations with momentum flip display a favorable behavior in terms of sampling efficiency, i.e., the traditional GHMC/GSHMC implementations with momentum flip got the advantage of a higher acceptance rate and faster decorrelation of Monte Carlo samples. The difference is more pronounced for GHMC. We also numerically investigate the behavior of the GHMC method as a Langevin-type thermostat. We find that the GHMC method without momentum flip interferes less with the underlying stochastic molecular dynamics in terms of autocorrelation functions and it to be preferred over the GHMC method with momentum flip. The same finding applies to GSHMC.

Original languageEnglish (US)
Pages (from-to)2256-2265
Number of pages10
JournalJournal of Computational Physics
Volume228
Issue number6
DOIs
StatePublished - Apr 1 2009
Externally publishedYes

All Science Journal Classification (ASJC) codes

  • Numerical Analysis
  • Modeling and Simulation
  • Physics and Astronomy (miscellaneous)
  • Physics and Astronomy(all)
  • Computer Science Applications
  • Computational Mathematics
  • Applied Mathematics

Keywords

  • Langevin dynamics
  • Molecular dynamics
  • Monte Carlo methods
  • Sampling

Fingerprint

Dive into the research topics of 'A comparison of generalized hybrid Monte Carlo methods with and without momentum flip'. Together they form a unique fingerprint.

Cite this