Determining efficient temperature sets for the simulated tempering method

Abstract

In statistical physics, the efficiency of tempering approaches strongly depends on ingredients such as the number of replicas $R$, reliable determination of weight factors and the set of used temperatures, ${\mathcal{T}}_R=\{T_1,T_2,\dots ,T_R\}$. For the simulated tempering (SP) in particular -- useful due to its generality and conceptual simplicity -- the latter aspect (closely related to the actual $R$) may be a key issue in problems displaying metastability and trapping in certain regions of the phase space. To determine ${\mathcal{T}}_R$'s leading to accurate thermodynamics estimates and still trying to minimize the simulation computational time, here it is considered a fixed exchange frequency scheme for the ST. From the temperature of interest $T_1$, successive $T$'s are chosen so that the exchange frequency between any adjacent pair $T_r$ and $T_{r+1}$ has a same value $f$. By varying the $f$'s and analyzing the ${\mathcal{T}}_R$'s through relatively inexpensive tests (e.g., time decay toward the steady regime), an optimal situation in which the simulations visit much faster and more uniformly the relevant portions of the phase space is determined. As illustrations, the proposal is applied to three lattice models, BEG, Bell-Lavis, and Potts, in the hard case of extreme first-order phase transitions, always giving very good results, even for $R=3$. Also, comparisons with other protocols (constant entropy and arithmetic progression) to choose the set ${\mathcal{T}}_R$ are undertaken. The fixed exchange frequency method is found to be consistently superior, specially for small $R$'s. Finally, distinct instances where the prescription could be helpful (in second-order transitions and for the parallel tempering approach) are briefly discussed.

10 pages, 14 figures

References in corpus (18)

• An efficient, multiple range random walk algorithm to calculate the density of states
• Exchange Monte Carlo Method and Application to Spin Glass Simulations
• Simulated Tempering: A New Monte Carlo Scheme
• Multicanonical Ensemble: A New Approach to Simulate First-order Phase Transitions
• Parallel Tempering: Theory, Applications, and New Perspectives
• The incomplete beta function law for parallel tempering sampling of classical canonical systems
• Make life simple: unleash the full power of the parallel tempering algorithm
• Liquid polymorphism, order-disorder transitions and anomalous behavior: a Monte Carlo study of the Bell-Lavis model for water
• Comparing parallel and simulated tempering enhanced sampling algorithms at phase transition regimes
• General approach for studying first-order phase transitions at low temperatures
• Diffusion anomaly and dynamic transitions in the Bell-Lavis water model
• Comparing different protocols of temperature selection in the parallel tempering method
• Exploiting a semi-analytic approach to study first order phase transitions
• Obtaining pressure versus concentration phase diagrams in spin systems from Monte Carlo simulations
• Simulated Tempering and Magnetizing: An Application of Two-Dimensional Simulated Tempering to Two-Dimensional Ising Model and Its Crossover
• A comparative study of selected parallel tempering methods
• A simple protocol for the probability weights of the simulated tempering algorithm: applications to first-order phase transitions
• Hydration and anomalous solubility of the Bell-Lavis model as solvent

Cited by in corpus (2)

• Irreversible simulated tempering
• Efficient simulated tempering with approximated weights: Applications to first-order phase transitions