Improving the efficiency of extended ensemble simulations: The accelerated weight histogram method

We propose a method for efficient simulations in extended ensembles, useful, e.g., for the study of problems with complex energy landscapes and for free energy calculations. The main difficulty in such simulations is the estimation of the a priori unknown weight parameters needed to produce flat histograms. The method combines several complementary techniques, namely, a Gibbs sampler for the parameter moves, a reweighting procedure to optimize data use, and a Bayesian update allowing for systematic refinement of the free energy estimate. In a certain limit the scheme reduces to the 1/t algorithm of B.E. Belardinelli and V.D. Pereyra [Phys. Rev. E 75, 046701 (2007)]. The performance of the method is studied on the two-dimensional Ising model, where comparison with the exact free energy is possible, and on an Ising spin glass.

💡 Research Summary

The paper introduces the Accelerated Weight Histogram (AWH) method, a novel scheme for efficiently sampling extended ensembles, which are widely used to overcome free‑energy barriers in systems with rugged energy landscapes such as spin glasses, polymers, and complex statistical models. The central difficulty in extended‑ensemble simulations is the need to determine a set of weight parameters fₘ (the dimensionless free energies) that flatten the marginal distribution of the auxiliary parameter (e.g., temperature). Existing dynamic schemes such as Wang‑Landau and the 1/t algorithm adjust these weights on the fly, but they suffer from slow convergence, histogram resetting, or saturation of the modification factor.

AWH resolves these issues by combining three complementary ideas:

1. Gibbs‑sampler based parameter moves – For a given configuration x, the transition probability to any parameter value m′ is computed directly from the conditional distribution P(m′|x) ∝ exp