Preprint Article Version 1 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.

Journal reference: Entropy 2017, 19, 561
DOI: 10.3390/e19100561

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.

Subject Areas

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

Readers' Comments and Ratings (0)

Leave a public comment
Send a private comment to the author(s)
Rate this article
Views 0
Downloads 0
Comments 0
Metrics 0
Leave a public comment

×
Alerts
Notify me about updates to this article or when a peer-reviewed version is published.