Dynamical traps in Wang-Landau sampling of continuous systems: Mechanism and solution
arXiv:1508.01888 · doi:10.1103/PhysRevE.92.023306
Abstract
We study the mechanism behind dynamical trappings experienced during Wang-Landau sampling of continuous systems reported by several authors. Trapping is caused by the random walker coming close to a local energy extremum, although the mechanism is different from that of critical slowing down encountered in conventional molecular dynamics or Monte Carlo simulations. When trapped, the random walker misses entire or even several stages of Wang-Landau modification factor reduction, leading to inadequate sampling of configuration space and a rough density-of-states even though the modification factor has been reduced to very small values. Trapping is dependent on specific systems, the choice of energy bins, and Monte Carlo step size, making it highly unpredictable. A general, simple, and effective solution is proposed where the configurations of multiple parallel Wang-Landau trajectories are inter-swapped to prevent trapping. We also explain why swapping frees the random walker from such traps. The efficacy of the proposed algorithm is demonstrated.
References in corpus (12)
- Determining the density of states for classical statistical models: A random walk algorithm to produce a flat histogram
- On the Wang-Landau Method for Off-Lattice Simulations in the "Uniform" Ensemble
- Understanding and Improving the Wang-Landau Algorithm
- Wang-Landau algorithm for continuous models and joint density of states
- Wang-Landau Algorithm: a Theoretical Analysis of the Saturation of the Error
- A generic, hierarchical framework for massively parallel Wang-Landau sampling
- Convergence and Refinement of the Wang-Landau Algorithm
- Wang-Landau Monte Carlo simulation of isotropic-nematic transition in liquid crystals
- Density of states of a binary Lennard-Jones Glass
- Wang-Landau molecular dynamics technique to search for low-energy conformational space of proteins
- Wang-Landau study of the 3D Ising model with bond disorder
- Scaling and self-averaging in the three-dimensional random-field Ising model