Feedback-optimized parallel tempering Monte Carlo

Abstract

We introduce an algorithm to systematically improve the efficiency of parallel tempering Monte Carlo simulations by optimizing the simulated temperature set. Our approach is closely related to a recently introduced adaptive algorithm that optimizes the simulated statistical ensemble in generalized broad-histogram Monte Carlo simulations. Conventionally, a temperature set is chosen in such a way that the acceptance rates for replica swaps between adjacent temperatures are independent of the temperature and large enough to ensure frequent swaps. In this paper, we show that by choosing the temperatures with a modified version of the optimized ensemble feedback method we can minimize the round-trip times between the lowest and highest temperatures which effectively increases the efficiency of the parallel tempering algorithm. In particular, the density of temperatures in the optimized temperature set increases at the "bottlenecks'' of the simulation, such as phase transitions. In turn, the acceptance rates are now temperature dependent in the optimized temperature ensemble. We illustrate the feedback-optimized parallel tempering algorithm by studying the two-dimensional Ising ferromagnet and the two-dimensional fully-frustrated Ising model, and briefly discuss possible feedback schemes for systems that require configurational averages, such as spin glasses.

12 pages, 14 figures

References in corpus (4)

• Optimizing the ensemble for equilibration in broad-histogram Monte Carlo simulations
• Performance Limitations of Flat Histogram Methods and Optimality of Wang-Landau Sampling
• Dynamics of the Wang-Landau algorithm and complexity of rare events for the three-dimensional bimodal Ising spin glass
• Analytical study of tunneling times in flat histogram Monte Carlo

Cited by in corpus (11)

• The ALPS project release 1.3: open source software for strongly correlated systems
• Generalized Ensemble and Tempering Simulations: A Unified View
• Make life simple: unleash the full power of the parallel tempering algorithm
• Using the average spectrum method to extract dynamics from quantum Monte Carlo simulations
• Phase diagram of the Ising square lattice with competing interactions
• Propagation and Ghosts in the Classical Kagome Antiferromagnet
• On Dynamics and Optimal Number of Replicas in Parallel Tempering Simulations
• Critical behavior of the random-anisotropy model in the strong-anisotropy limit
• First-order phase transitions: A study through the parallel tempering method
• Optimized Broad-Histogram Ensembles for the Simulation of Quantum Systems
• Optimized Folding Simulations of Protein A