Hard-disk dipoles and non-reversible Markov chains
arXiv:2111.11943 · doi:10.1063/5.0080101
Abstract
We benchmark event-chain Monte Carlo (ECMC) algorithms for tethered hard-disk dipoles in two dimensions in view of application of ECMC to water models in molecular simulation. We characterize the rotation dynamics of dipoles through the integrated autocorrelation times of the polarization. The non-reversible straight, reflective, forward, and Newtonian ECMC algorithms are all event-driven, and they differ only in their update rules at event times. They realize considerable speedups with respect to the local reversible Metropolis algorithm. We also find significant speed differences among the ECMC variants. Newtonian ECMC appears particularly well-suited for overcoming the dynamical arrest that has plagued straight ECMC for three-dimensional dipolar models with Coulomb interactions.
9 pages, 7 figures
References in corpus (7)
- Event-chain Monte Carlo: foundations, applications, and prospects
- Efficient Equilibration of Hard Spheres with Newtonian Event Chains
- Event-Chain Monte-Carlo Simulations of Dense Soft Matter Systems
- Chaining of hard disks in nematic needles: particle-based simulation of colloidal interactions in liquid crystals
- Anomalous diffusion analysis of the lifting events in the event-chain Monte Carlo for the classical XY models
- Sparse hard-disk packings and local Markov chains
- Direction-sweep Markov chains
Cited by in corpus (8)
- Accelerating equilibrium spin-glass simulations using quantum annealers via generative deep learning
- Molecular simulation from modern statistics: Continuous-time, continuous-space, exact
- Sparse hard-disk packings and local Markov chains
- Direction-sweep Markov chains
- Sampling algorithms in statistical physics: a guide for statistics and machine learning
- Microscopic Mechanisms of Diffusion Dynamics: A Comparative Efficiency Study of Event-Chain Monte Carlo Variants in Dense Hard Disk Systems
- Comparison of the clock, stochastic cutoff, and Tomita Monte Carlo methods in simulating the dipolar triangular lattice at criticality
- Necessary and sufficient symmetries in Event-Chain Monte Carlo with generalized flows and Application to hard dimers