Krauth 2024 HMC ECMC
From Werner KRAUTH
Revision as of 00:11, 26 November 2024 Werner (Talk | contribs) ← Previous diff |
Revision as of 00:11, 26 November 2024 Werner (Talk | contribs) Next diff → |
||
Line 1: | Line 1: | ||
'''W. Krauth''' '''''Hamiltonian Monte Carlo vs. event-chain Monte Carlo: | '''W. Krauth''' '''''Hamiltonian Monte Carlo vs. event-chain Monte Carlo: | ||
- | an appraisal of sampling strategies beyond the diffusive regime ''''' ''' (2024)''' | + | an appraisal of sampling strategies beyond the diffusive regime ''''' |
arXiv:2411.11690 (2024) | arXiv:2411.11690 (2024) | ||
Revision as of 00:11, 26 November 2024
'W. Krauth Hamiltonian Monte Carlo vs. event-chain Monte Carlo: an appraisal of sampling strategies beyond the diffusive regime arXiv:2411.11690 (2024)
Abstract We discuss Hamiltonian Monte Carlo (HMC) and event-chain Monte Carlo (ECMC) for the one-dimensional chain of particles with harmonic interactions and benchmark them against local reversible Metropolis algorithms. While HMC achieves considerable speedup with respect to local reversible Monte Carlo algorithms, its autocorrelation functions of global observables such as the structure factor have slower scaling with system size than for ECMC, a lifted non-reversible Markov chain. This can be traced to the dependence of ECMC on a parameter of the harmonic energy, the equilibrium distance, which drops out when energy differences or gradients are evaluated. We review the recent literature and provide pseudocodes and Python programs. We finally discuss related models and generalizations beyond one-dimensional particle systems.