Cluster algorithms: Difference between revisions

From SklogWiki
Jump to navigation Jump to search
Line 18: Line 18:
depending on the product: <math> S_{i} S_{j} </math> (See [[Ising Models]] for details on the notation)
depending on the product: <math> S_{i} S_{j} </math> (See [[Ising Models]] for details on the notation)


* For pairs with <math> u_{ij}/k_B T < 0 </math> create a bond between the two spins with a given probability <math> p </math>
* For pairs of interacting sites (nearest neighbors) with <math> u_{ij}/k_B T < 0 </math> create a bond between the two spins with a given probability <math> p </math> using [[random numbers]])
 
: <math> p </math> will be chosen to be a function of <math> K </math>
* The bonds generated in the previous step are used to build up clusters of sites (spins).
 
* Build up the partition of the sites in different clusters
 
In each cluster all the spins will have the same state (either <math> S = 1 </math> or <math> S = -1 </math>)
 
* For each cluster, independtly, choose at random (with proability 1/2) to flip (invert the value of <math> S </math>) or not to flip, every spin on the cluster


== Wolf algorithm ==
== Wolf algorithm ==

Revision as of 19:31, 3 August 2007

WORKING ON THIS PAGE ...

Cluster algorithms in Monte Carlo Simulation.

These algorithms are mainly used in the simulation of Ising-like models. The essential feature is the use of collective motions of particles (spins) in a single Monte Carlo step.

An interesting property of some of these application is the fact that the percolation analysis of the clusters can be used to study phase transitions.

As an introductory example we will discuss the Swendsen-Wang technique (Ref 1) in the simulation of Ising Models.

Sketches of the Swendsen-Wang algorithm

  • Consider every pair interacting sites (spins)

In the current configuration the pair interaction can be either negative: of positive , depending on the product: (See Ising Models for details on the notation)

  • For pairs of interacting sites (nearest neighbors) with create a bond between the two spins with a given probability using random numbers)
will be chosen to be a function of
  • The bonds generated in the previous step are used to build up clusters of sites (spins).
  • Build up the partition of the sites in different clusters

In each cluster all the spins will have the same state (either or )

  • For each cluster, independtly, choose at random (with proability 1/2) to flip (invert the value of ) or not to flip, every spin on the cluster

Wolf algorithm

Invaded Cluter Algorithm

References

  1. Robert H. Swendsen and Jian-Sheng Wang, Nonuniversal critical dynamics in Monte Carlo simulations, Phys. Rev. Lett. 58, 86 - 88 (1987)