Efficient sampling using Macrocanonical Monte Carlo and density of states mapping

Abstract

In the context of Monte Carlo sampling for lattice models, the complexity of the energy landscape often leads to Markov chains being trapped in local optima, thereby increasing the correlation between samples and reducing sampling efficiency. This study proposes a Monte Carlo algorithm that effectively addresses the irregularities of the energy landscape through the introduction of the estimated density of states. This algorithm enhances the accuracy in the study of phase transitions and is not model-specific. Although our algorithm is primarily demonstrated on the two-dimensional square lattice model, the method is also applicable to a broader range of lattice and higher-dimensional models. Furthermore, the study develops a method for estimating the density of states of large systems based on that of smaller systems, enabling high-precision density of states estimation within specific energy intervals in large systems without sampling. For regions of lower precision, a re-weighting strategy is employed to adjust the density of states to enhance the precision further. This algorithm is not only significant within the field of lattice model sampling but may also inspire applications of the Monte Carlo method in other domains.

arXiv comment14 pages, 9 figures

References in corpus (25)

• An efficient, multiple range random walk algorithm to calculate the density of states
• Multicanonical Ensemble: A New Approach to Simulate First-order Phase Transitions
• Self-Learning Monte Carlo Method
• Fast Algorithm to Calculate Density of States
• Transition Matrix Monte Carlo Method
• Two phase transitions in the fully frustrated $XY$ model
• The Dynamic Exponent of the Two-Dimensional Ising Model and Monte Carlo Computation of the Sub-Dominant Eigenvalue of the Stochastic Matrix
• A generic, hierarchical framework for massively parallel Wang-Landau sampling
• Multi-GPU Accelerated Multi-Spin Monte Carlo Simulations of the 2D Ising Model
• Multicritical behavior in the fully frustrated XY model and related systems
• Scalable replica-exchange framework for Wang-Landau sampling
• Large-scale Monte Carlo simulation of two-dimensional classical XY model using multiple GPUs
• Finding Low-Temperature States with Parallel Tempering, Simulated Annealing and Simple Monte Carlo
• Kosterlitz-Thouless and Potts transitions in a generalized XY model
• Wang-Landau sampling: Improving accuracy
• Tensor network approach to the two-dimensional fully frustrated XY model and a chiral ordered phase
• Wang-Landau sampling: a criterion for halting the simulations
• Spin-chirality decoupling and critical properties of a two-dimensional fully frustrated XY model
• Control of accuracy in the Wang-Landau algorithm
• Finite versus zero-temperature hysteretic behavior of spin glasses: Experiment and theory
• Successive transitions and intermediate chiral phase in a superfuilid ^3He film
• An Efficient Simulation Protocol for Determining the Density of States: Combination of Replica-Exchange Wang-Landau Method and Multicanonical Replica-Exchange Method
• The Wang-Landau Algorithm as Stochastic Optimization and Its Acceleration
• Dynamically Optimized Wang-Landau Sampling with Adaptive Trial Moves and Modification Factors
• A New Monte Carlo Algorithm for Free Energy Calculation