Wang-Landau method: Difference between revisions

From SklogWiki
Jump to navigation Jump to search
m (more references)
Line 45: Line 45:
== Extensions ==
== Extensions ==
The Wang-Landau method has inspired a number of simulation algorithms that
The Wang-Landau method has inspired a number of simulation algorithms that
use the same strategy in different contexts.
use the same strategy in different contexts. For example:


* [[Inverse Monte Carlo|Inverse Monte Carlo]] methods
* [[Inverse Monte Carlo|Inverse Monte Carlo]] methods

Revision as of 11:36, 10 July 2008

The Wang-Landau method was proposed by F. Wang and D. P. Landau (Ref. 1-2) to compute the density of states, , of Potts models; where is the number of microstates of the system having energy .

Sketches of the method

The Wang-Landau method, in its original version, is a simulation technique designed to achieve a uniform sampling of the energies of the system in a given range. In a standard Metropolis Monte Carlo in the canonical ensemble the probability of a given microstate, , is given by:

;

whereas for the Wang-Landau procedure one can write:

 ;

where is a function of the energy. changes during the simulation in order produce a predefined distribution of energies (usually a uniform distribution); this is done by modifying the values of to reduce the probability of the energies that have been already visited, i.e. If the current configuration has energy , is uptdated as:

 ;

where it has been considered that the system has discrete values of the energy (as happens in Potts Models), and .

Such a simple scheme is continued until the shape of the energy distribution approaches the one predefined. Notice that this simulation scheme does not produce an equilibrium procedure, since it does not fulfil detailed balance. To overcome this problem, the Wang-Landau procedure consists in the repetition of the scheme sketched above along several stages. In each subsequent stage the perturbation parameter is reduced. So, for the last stages the function hardly changes and the simulation results of these last stages can be considered as a good description of the actual equilibrium system, therefore:

;

where , and is the Kronecker Delta.

If the probability distribution of energies is nearly unifom: ; then

Extensions

The Wang-Landau method has inspired a number of simulation algorithms that use the same strategy in different contexts. For example:

  • Computation of the phase equibria of fluids
  • Control of polydispersity by chemical potential tunning

References

  1. Fugao Wang and D. P. Landau "Efficient, Multiple-Range Random Walk Algorithm to Calculate the Density of States", Phys. Rev. Lett. 86, 2050 - 2053 (2001)
  2. Fugao Wang and D. P. Landau "Determining the density of states for classical statistical models: A random walk algorithm to produce a flat histogram", Physical Review E 64 056101 (2001)
  3. D. P. Landau, Shan-Ho Tsai, and M. Exler "A new approach to Monte Carlo simulations in statistical physics: Wang-Landau sampling", American Journal of Physics 72 pp. 1294-1302 (2004)
  4. Georg Ganzenmüller and Philip J. Camp "Applications of Wang-Landau sampling to determine phase equilibria in complex fluids", Journal of Chemical Physics 127 154504 (2007)
  5. R. E. Belardinelli and V. D. Pereyra "Wang-Landau algorithm: A theoretical analysis of the saturation of the error", Journal of Chemical Physics 127 184105 (2007)
  6. R. E. Belardinelli and V. D. Pereyra "Fast algorithm to calculate density of states", Physical Review E 75 046701 (2007)