English

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

Statistical Mechanics 2025-02-03 v3 Mathematical Physics math.MP

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.

Keywords

Cite

@article{arxiv.2402.11267,
  title  = {On the optimal relaxation rate for the Metropolis algorithm in one dimension},
  author = {A. Patrón and A. D. Chepelianskii and A. Prados and E. Trizac},
  journal= {arXiv preprint arXiv:2402.11267},
  year   = {2025}
}

Comments

25 pages, 10 figures, 1 table

R2 v1 2026-06-28T14:51:46.558Z