Michel Kapfer Krauth 2013
From Werner KRAUTH
Revision as of 13:36, 1 April 2014 Werner (Talk | contribs) ← Previous diff |
Revision as of 20:17, 3 May 2014 Werner (Talk | contribs) Next diff → |
||
Line 8: | Line 8: | ||
[http://arxiv.org/pdf/1309.7748 Electronic version (from arXiv, original version)] | [http://arxiv.org/pdf/1309.7748 Electronic version (from arXiv, original version)] | ||
+ | |||
+ | [http://scitation.aip.org/content/aip/journal/jcp/140/5/10.1063/1.4863991 Published version (Journal of Chemical Physics - Subscription needed) |
Revision as of 20:17, 3 May 2014
M. Michel, S. C. Kapfer, W. Krauth Generalized event-chain Monte Carlo: Constructing rejection-free global-balance algorithms from infinitesimal steps Journal of Chemical Physics 140 54116 (2014)
Paper
Abstract In this article, we introduce a factorized Metropolis filter and the concept of infinitesimal Monte Carlo moves to design a rejection-free Markov-chain Monte Carlo algorithm for interacting particle systems that breaks detailed balance yet satisfies global balance. This event-driven algorithm generalizes the recent hard-sphere event-chain Monte Carlo method without introducing any discretizations in time or in space. We demonstrate considerable speed-ups of this method with respect to the classic local Metropolis algorithm. The new algorithm generates a continuum of samples of the stationary probability density. This allows us to derive an exact formula for the pressure that is obtained as a byproduct of the simulation without any additional computations.
Electronic version (from arXiv, original version)
[http://scitation.aip.org/content/aip/journal/jcp/140/5/10.1063/1.4863991 Published version (Journal of Chemical Physics - Subscription needed)