ECMC 2021 Hoellmer
From Werner KRAUTH
Dipole rotations in non-reversible Markov Chains
- Philipp Höllmer, Bethe Center for Theoretical Physics, University of Bonn, Nussallee 12, 53115 Bonn, Germany
- A. C. Maggs, CNRS UMR7083, ESPCI Paris, PSL Research University, 10 rue Vauquelin, 75005 Paris, France
- Werner Krauth, Laboratoire de Physique de l’Ecole normale supérieure, ENS, Université PSL, CNRS, Sorbonne Université, Université Paris-Diderot, Sorbonne Paris Cité, Paris, France
Abstract: In recent years, several algorithms belonging to the family of event-chain Monte Carlo (ECMC) have been proposed. They share the concept of particle and displacement lifting where an active particle is displaced on a non-interacting trajectory in a lifted direction. The individual algorithms differ, however, in the update scheme of the lifting variables in the events that are required to construct a non-reversible continuous time Markov chain. We study straight ECMC with direction sweeps, reflective ECMC, Newtonian ECMC, and forward ECMC for many simplified two-dimensional extended flexible dipoles that resemble water molecules in the context of molecular simulation models. Here, we point out possible pitfalls of the straight and reflective versions. We show that the dynamics of the polarization of many dipoles are analogous to a simple Gaussian random walk and a path integral. The polarization's integrated autocorrelation times show that straight ECMC, which was the superior variant for the problem of two-dimensional hard disks, rotates the dipoles slowly in comparison to the other methods. In comparison to a reversible Metropolis algorithm, the optimal ECMC algorithm yields a speedup that increases for lower densities up to a factor of $60$, which thus motivates its application to systems of long-range interacting extended molecules at the core of the JeLLyFysh project.
Slides
Recording
Further material
