On the optimal relaxation rate for the Metropolis algorithm in one dimension

Abstract

We study the relaxation of the Metropolis Monte Carlo algorithm corresponding to a single particle trapped in a one-dimensional confining potential, with even jump distributions that ensure that the dynamics verifies detailed balance. Previous work suggested that, for smooth jump distributions, the fastest relaxation rate is obtained as a result of the competition between diffusive and rejection-dominated dynamics. In this work, we show that a new regime comes into play for two-peaked jump distributions, where the relaxation dynamics is neither dominated by diffusion nor rejection: the eigenmodes adopt an oscillatory form, reminiscent of charge density waves (CDW) -- thus we term this new regime the CDW regime. Using a combination of numerical and analytical techniques, the parameter regions corresponding to diffusion, rejection, and CDW are characterised, as well as the transition lines between them -- i.e. a phase diagram is built. The optimal relaxation rate is located at the triple point of phase coexistence, where the transition lines (diffusive-rejection, diffusive-CDW, and CDW-rejection) intersect. Our theoretical framework is checked versus the numerical diagonalisation of the master equation. We also briefly discuss more sophisticated attempts at optimising the relaxation rate to equilibrium.

25 pages, 10 figures, 1 table

References in corpus (17)

• Shortcuts to adiabaticity: concepts, methods, and applications
• Intermittent search strategies
• Stochastic Resetting and Applications
• Anomalous cooling and heating - the Mpemba effect and its inverse
• When the Hotter Cools More Quickly: Mpemba Effect in Granular Fluids
• Exponentially faster cooling in a colloidal system
• The Mpemba index and anomalous relaxation
• Pre-Cooling Strategy Allows Exponentially Faster Heating
• Anomalous heating in a colloidal system
• Driving rapidly while remaining in control: classical shortcuts from Hamiltonian to stochastic dynamics
• Eigenvalue crossing as a phase transition in relaxation dynamics
• The energy level crossing behavior and quantum Fisher information in a quantum well with spin-orbit coupling
• Mpemba effect in a Langevin system: population statistics, metastability and other exact results
• Non-equilibrium memory effects: granular fluids and beyond
• Optimum Monte Carlo Simulations: Some Exact Results
• Optimal resetting strategies for search processes in heterogeneous environments
• Metropolis Monte Carlo sampling: convergence, localization transition and optimality