Random cluster model


In statistical mechanics, probability theory, graph theory, etc. the random cluster model is a random graph that generalizes and unifies the Ising model, Potts model, and percolation model. It is used to study random combinatorial structures, electrical networks, etc. It is also referred to as the RC model or sometimes the FK representation after its founders Cees Fortuin and Piet Kasteleyn. The random cluster model has a critical limit, described by a conformal field theory.

Definition

Let be a graph, and be a bond configuration on the graph that maps each edge to a value of either 0 or 1. We say that a bond is closed on edge if, and open if. If we let be the set of open bonds, then an open cluster or FK cluster is any connected component in union the set of vertices. Note that an open cluster can be a single vertex.
Suppose an edge is open independently with probability and closed otherwise, then this is just the standard Bernoulli percolation process. The probability measure of a configuration is given as
The RC model is a generalization of percolation, where each cluster is weighted by a factor of. Given a configuration, we let be the number of open clusters, or alternatively the number of connected components formed by the open bonds. Then for any, the probability measure of a configuration is given as
Z is the partition function, or the sum over the unnormalized weights of all configurations,
The partition function of the RC model is a specialization of the Tutte polynomial, which itself is a specialization of the multivariate Tutte polynomial.

Special values of ''q''

The parameter of the random cluster model can take arbitrary complex values. This includes the following special cases:
  • : linear resistance networks.
  • : negatively-correlated percolation.
  • : Bernoulli percolation, with.
  • : the Ising model.
  • : -state Potts model.

    Edwards-Sokal representation

The Edwards-Sokal representation of the Potts model is named after Robert G. Edwards and Alan D. Sokal. It provides a unified representation of the Potts and random cluster models in terms of a joint distribution of spin and bond configurations.
Let be a graph, with the number of vertices being and the number of edges being. We denote a spin configuration as and a bond configuration as. The joint measure of is given as
where is the uniform measure, is the product measure with density, and is an appropriate normalizing constant. Importantly, the indicator function of the set
enforces the constraint that a bond can only be open on an edge if the adjacent spins are of the same state, also known as the SW rule.
The statistics of the Potts spins can be recovered from the cluster statistics, thanks to the following features of the ES representation:
  • The marginal measure of the spins is the Boltzmann measure of the q-state Potts model at inverse temperature.
  • The marginal measure of the bonds is the random-cluster measure with parameters q and p.
  • The conditional measure of the spin represents a uniformly random assignment of spin states that are constant on each connected component of the bond arrangement.
  • The conditional measure of the bonds represents a percolation process on the subgraph of formed by the edges where adjacent spins are aligned.
  • In the case of the Ising model, the probability that two vertices are in the same connected component of the bond arrangement equals the two-point correlation function of spins, written.

    Frustration

There are several complications of the ES representation once frustration is present in the spin model. In particular, there is no longer a correspondence between the spin statistics and the cluster statistics, and the correlation length of the RC model will be greater than the correlation length of the spin model. This is the reason behind the inefficiency of the SW algorithm for simulating frustrated systems.

Two-dimensional case

If the underlying graph is a planar graph, there is a duality between the random cluster models on and on the dual graph. At the level of the partition function, the duality reads
On a self-dual graph such as the square lattice, a phase transition can only occur at the self-dual coupling.
The random cluster model on a planar graph can be reformulated as a loop model on the corresponding medial graph. For a configuration of the random cluster model, the corresponding loop configuration is the set of self-avoiding loops that separate the clusters from the dual clusters. In the transfer matrix approach, the loop model is written in terms of a Temperley-Lieb algebra with the parameter. In two dimensions, the random cluster model is therefore closely related to the O model, which is also a loop model.
In two dimensions, the critical random cluster model is described by a conformal field theory with the central charge
Known exact results include the conformal dimensions of the fields that detect whether a point belongs to an FK cluster or a spin cluster. In terms of Kac indices, these conformal dimensions are respectively and, corresponding to the fractal dimensions and of the clusters.

History and applications

RC models were introduced in 1969 by Fortuin and Kasteleyn, mainly to solve combinatorial problems. After their founders, it is sometimes referred to as FK models. In 1971 they used it to obtain the FKG inequality. Post 1987, interest in the model and applications in statistical physics reignited. It became the inspiration for the Swendsen–Wang algorithm describing the time-evolution of Potts models. Michael Aizenman and coauthors used it to study the phase boundaries in 1D Ising and Potts models.