Preprint Article Version 1 Preserved in Portico This version is not peer-reviewed

Comparing Markov Chain Samplers for Molecular Simulation

Version 1 : Received: 6 September 2017 / Approved: 7 September 2017 / Online: 7 September 2017 (03:43:35 CEST)

A peer-reviewed article of this Preprint also exists.

Skeel, R.D.; Fang, Y. Comparing Markov Chain Samplers for Molecular Simulation. Entropy 2017, 19, 561. Skeel, R.D.; Fang, Y. Comparing Markov Chain Samplers for Molecular Simulation. Entropy 2017, 19, 561.

Abstract

Markov chain Monte Carlo sampling propagators, including numerical integrators for stochastic dynamics, are central to the calculation of thermodynamic quantities and determination of structure for molecular systems. Efficiency is paramount, and to a great extent, this is determined by the integrated autocorrelation time (IAcT). This quantity varies depending on the observable that is being estimated. It is suggested that it is the maximum of the IAcT over all observables that is the relevant metric. Reviewed here is a method for estimating this quantity. For reversible propagators (which are those that satisfy detailed balance), the maximum IAcT is determined by the spectral gap in the forward transfer operator, but for irreversible propagators, the maximum IAcT can be far less than or greater than what might be inferred from the spectral gap. This is consistent with recent theoretical results (not to mention past practical experience) suggesting that irreversible propagators generally perform better if not much better than reversible ones. Typical irreversible propagators involve a parameter controlling the mix of ballistic and diffusive movement. To gain insight into the effect of the damping parameter for Langevin dynamics, its optimal value is obtained here for a multidimensional quadratic potential energy function.

Keywords

Markov chain Monte Carlo; stochastic dynamics integrators; decorrelation time; integrated autocorrelation time

Subject

Physical Sciences, Atomic and Molecular Physics

Comments (0)

We encourage comments and feedback from a broad range of readers. See criteria for comments and our Diversity statement.

Leave a public comment
Send a private comment to the author(s)
* All users must log in before leaving a comment
Views 0
Downloads 0
Comments 0
Metrics 0


×
Alerts
Notify me about updates to this article or when a peer-reviewed version is published.
We use cookies on our website to ensure you get the best experience.
Read more about our cookies here.