Editing Potts model

{{short description|Model in statistical mechanics generalizing the Ising model}}
In [[statistical mechanics]], the '''Potts model''', a generalization of the [[Ising model]], is a model of interacting [[Spin (physics)|spins]] on a [[crystalline lattice]].<ref>{{Cite journal |last=Wu |first=F. Y. |date=1982-01-01 |title=The Potts model |url=https://link.aps.org/doi/10.1103/RevModPhys.54.235 |journal=Reviews of Modern Physics |volume=54 |issue=1 |pages=235–268 |doi=10.1103/RevModPhys.54.235|bibcode=1982RvMP...54..235W }}</ref> By studying the Potts model, one may gain insight into the behaviour of [[ferromagnet]]s and certain other phenomena of [[solid-state physics]].  The strength of the Potts model is not so much that it models these physical systems well; it is rather that the one-dimensional case is [[exactly solvable]], and that it has a rich mathematical formulation that has been studied extensively.

The model is named after [[Renfrey Potts]], who described the model near the end of his 1951 Ph.D. thesis.<ref>{{Cite journal |last1=Potts |first1=R. B. | authorlink1=Renfrey Potts |date=January 1952 |title=Some generalized order-disorder transformations |url=https://www.cambridge.org/core/journals/mathematical-proceedings-of-the-cambridge-philosophical-society/article/abs/some-generalized-orderdisorder-transformations/5FD50240095F40BD123171E5F76CDBE0 |journal=[[Mathematical Proceedings of the Cambridge Philosophical Society]] |language=en |volume=48 |issue=1 |pages=106–109 |doi=10.1017/S0305004100027419 |bibcode=1952PCPS...48..106P |s2cid=122689941 |issn=1469-8064}}</ref> The model was related to the "planar Potts" or "[[Z N model|clock model]]", which was suggested to him by his advisor, [[Cyril Domb]]. The four-state Potts model is sometimes known as the '''Ashkin–Teller model''',<ref>{{Cite journal |last1=Ashkin |first1=J. |last2=Teller |first2=E. |date=1943-09-01 |title=Statistics of Two-Dimensional Lattices with Four Components |url=https://link.aps.org/doi/10.1103/PhysRev.64.178 |journal=Physical Review |volume=64 |issue=5–6 |pages=178–184 |doi=10.1103/PhysRev.64.178|bibcode=1943PhRv...64..178A }}</ref> after [[Julius Ashkin]] and [[Edward Teller]], who considered an equivalent model in 1943.

The Potts model is related to, and generalized by, several other models, including the [[XY model]], the [[Heisenberg model (classical)|Heisenberg model]] and the [[N-vector model]]. The infinite-range Potts model is known as the [[Kac model]]. When the spins are taken to interact in a [[non-abelian group|non-Abelian]] manner, the model is related to the [[flux tube model]], which is used to discuss [[Color confinement|confinement]] in [[quantum chromodynamics]]. Generalizations of the Potts model have also been used to model [[grain growth]] in metals, [[coarsening]] in [[foam]]s, and statistical properties of [[Protein structure prediction|proteins]].<ref name=":1" /> A further generalization of these methods by [[James Glazier]] and [[Francois Graner]], known as the [[cellular Potts model]],<ref>{{Cite journal |last1=Graner |first1=François |last2=Glazier |first2=James A. |date=1992-09-28 |title=Simulation of biological cell sorting using a two-dimensional extended Potts model |url=https://link.aps.org/doi/10.1103/PhysRevLett.69.2013 |journal=Physical Review Letters |volume=69 |issue=13 |pages=2013–2016 |doi=10.1103/PhysRevLett.69.2013|pmid=10046374 |bibcode=1992PhRvL..69.2013G }}</ref> has been used to simulate static and kinetic phenomena in foam and biological [[morphogenesis]].

== Definition ==

=== Vector Potts model ===
The Potts model consists of ''spins'' that are placed on a [[lattice (group)|lattice]]; the lattice is usually taken to be a two-dimensional rectangular [[Euclidean space|Euclidean]] lattice, but is often generalized to other dimensions and lattice structures.

Originally, Domb suggested that the spin takes one of <math>q</math> possible values {{Citation needed|date=May 2022}}, distributed uniformly about the [[circle]], at angles
: <math>\theta_s = \frac{2\pi s}{q},</math>
where <math>s = 0, 1, ..., q-1</math> and that the interaction [[Hamiltonian mechanics|Hamiltonian]] is given by
: <math>H_c = J_c\sum_{\langle i, j \rangle} \cos \left( \theta_{s_i} - \theta_{s_j} \right)</math>
with the sum running over the nearest neighbor pairs <math>\langle i,j \rangle</math> over all lattice sites, and <math>J_c</math> is a coupling constant, determining the interaction strength. This model is now known as the '''vector Potts model''' or the '''clock model'''. Potts provided the location in two dimensions of the phase transition for <math>q = 3,4</math>. In the limit <math>q \to \infty</math>, this becomes the [[XY model]].

=== Standard Potts model ===
What is now known as the standard '''Potts model''' was suggested by Potts in the course of his study of the model above and is defined by a simpler Hamiltonian:
: <math>H_p = -J_p \sum_{(i,j)}\delta(s_i,s_j) \,</math>
where <math>\delta(s_i, s_j)</math> is the [[Kronecker delta]], which equals one whenever <math>s_i = s_j</math> and zero otherwise.

The <math>q=2</math> standard Potts model is equivalent to the [[Ising model]] and the 2-state vector Potts model, with <math>J_p = -2J_c</math>. The <math>q=3</math> standard Potts model is equivalent to the three-state vector Potts model, with <math>J_p = -\frac{3}{2}J_c</math>.

=== Generalized Potts model ===
A generalization of the Potts model is often used in statistical inference and biophysics, particularly for modelling proteins through [[direct coupling analysis]].<ref name=":1">{{Cite journal |last1=Shimagaki |first1=Kai |last2=Weigt |first2=Martin |date=2019-09-19 |title=Selection of sequence motifs and generative Hopfield-Potts models for protein families |url=https://link.aps.org/doi/10.1103/PhysRevE.100.032128 |journal=Physical Review E |volume=100 |issue=3 |pages=032128 |arxiv=1905.11848 |bibcode=2019PhRvE.100c2128S |doi=10.1103/PhysRevE.100.032128 |pmid=31639992 |s2cid=167217593}}</ref><ref>{{Cite journal |last1=Mehta |first1=Pankaj |last2=Bukov |first2=Marin |last3=Wang |first3=Ching-Hao |last4=Day |first4=Alexandre G. R. |last5=Richardson |first5=Clint |last6=Fisher |first6=Charles K. |last7=Schwab |first7=David J. |date=2019-05-30 |title=A high-bias, low-variance introduction to Machine Learning for physicists |journal=Physics Reports |volume=810 |pages=1–124 |arxiv=1803.08823 |bibcode=2019PhR...810....1M |doi=10.1016/j.physrep.2019.03.001 |issn=0370-1573 |pmc=6688775 |pmid=31404441}}</ref> This generalized Potts model consists of 'spins' that each may take on <math>q</math> states: <math>s_i \in \{1,\dots,q\}</math> (with no particular ordering). The Hamiltonian is,
: <math>
H = \sum_{i < j} J_{ij}(s_i,s_j) + \sum_i h_i(s_i),
</math>
where <math>J_{ij}(k,k')</math> is the energetic cost of spin <math>i</math> being in state <math>k</math> while spin <math>j</math> is in state <math>k'</math>, and <math>h_i(k)</math> is the energetic cost of spin <math>i</math> being in state <math>k</math>. Note: <math>J_{ij}(k,k') = J_{ji}(k',k)</math>. This model resembles the [[Sherrington-Kirkpatrick model]] in that couplings can be heterogeneous and non-local. There is no explicit lattice structure in this model.

== Physical properties ==
=== Phase transitions ===
Despite its simplicity as a model of a physical system, the Potts model is useful as a model system for the study of [[phase transition]]s. For example, for the standard ferromagnetic Potts model in <math>2d</math>, a phase transition exists for all real values <math>q \geq 1</math>,<ref name=":0">{{Cite journal |last1=Beffara |first1=Vincent |last2=Duminil-Copin |first2=Hugo |date=2012-08-01 |title=The self-dual point of the two-dimensional random-cluster model is critical for q ≥ 1 |journal=Probability Theory and Related Fields |language=en |volume=153 |issue=3 |pages=511–542 |doi=10.1007/s00440-011-0353-8 |s2cid=55391558 |issn=1432-2064|doi-access=free }}</ref> with the critical point at <math>\beta J = \log(1 + \sqrt{q})</math>. The phase transition is continuous (second order) for <math>1 \leq q \leq 4</math> <ref>{{Cite journal |last1=Duminil-Copin |first1=Hugo |last2=Sidoravicius |first2=Vladas |last3=Tassion |first3=Vincent |date=2017-01-01 |title=Continuity of the Phase Transition for Planar Random-Cluster and Potts Models with <math>{1 \le q \le 4}</math> |url=https://doi.org/10.1007/s00220-016-2759-8 |journal=Communications in Mathematical Physics |language=en |volume=349 |issue=1 |pages=47–107 |doi=10.1007/s00220-016-2759-8 |arxiv=1505.04159 |s2cid=119153736 |issn=1432-0916}}</ref> and discontinuous (first order) for <math>q > 4</math>.<ref>{{cite journal
| last1=Duminil-Copin | first1=Hugo | authorlink1=Hugo Duminil-Copin
| last2=Gagnebin | first2=Maxime
| last3=Harel | first3=Matan
| last4=Manolescu | first4=Ioan
| last5=Tassion | first5=Vincent
| date=2021
| title=Discontinuity of the phase transition for the planar random-cluster and Potts models with <math>q>4</math>
| arxiv=1611.09877
| journal=Annales Scientifiques de l'École Normale Supérieure
| volume=54
| issue=6
| pages=1363-1413
| doi=10.24033/asens.2485}}</ref>

For the clock model, there is evidence that the corresponding phase transitions are infinite order [[BKT transition]]s,<ref name="lyxt19" /> and a continuous phase transition is observed when <math>q \leq 4</math>.<ref name="lyxt19" /> Further use is found through the model's relation to [[percolation theory|percolation]] problems and the [[Tutte polynomial|Tutte]] and [[chromatic polynomial]]s found in combinatorics. For integer values of <math>q \geq 3</math>, the model displays the phenomenon of 'interfacial adsorption' <ref>{{Cite journal |last1=Selke |first1=Walter |last2=Huse |first2=David A. |date=1983-06-01 |title=Interfacial adsorption in planar potts models |url=https://doi.org/10.1007/BF01304093 |journal=Zeitschrift für Physik B: Condensed Matter |language=en |volume=50 |issue=2 |pages=113–116 |doi=10.1007/BF01304093 |bibcode=1983ZPhyB..50..113S |s2cid=121502987 |issn=1431-584X}}</ref> with intriguing critical [[wetting]] properties when fixing opposite boundaries in two different states {{Clarify|reason=is this for the standard potts or clock?|date=May 2022}}.

=== Relation with the random cluster model ===

The Potts model has a close relation to the Fortuin-[[Pieter Kasteleyn|Kasteleyn]] [[random cluster model]], another model in [[statistical mechanics]]. Understanding this relationship has helped develop efficient [[Markov chain Monte Carlo]] methods for numerical exploration of the model at small <math>q</math>, and led to the rigorous proof of the critical temperature of the model.<ref name=":0" />

At the level of the partition function <math>Z_p = \sum_{\{s_i\}} e^{-H_p}</math>, the relation amounts to transforming the sum over spin configurations <math>\{s_i\}</math> into a sum over edge configurations <math>\omega=\Big\{(i,j)\Big|s_i=s_j\Big\}</math> i.e. sets of nearest neighbor pairs of the same color. The transformation is done using the identity<ref>{{cite book |last=Sokal |first=Alan D. |title=Surveys in Combinatorics 2005 |chapter=The multivariate Tutte polynomial (alias Potts model) for graphs and matroids |year=2005 |arxiv=math/0503607 |pages=173–226 |doi=10.1017/CBO9780511734885.009|isbn=9780521615235 |s2cid=17904893 }}</ref> 
: <math>
e^{J_p\delta(s_i,s_j)} = 1 + v \delta(s_i,s_j) \qquad \text{ with } \qquad v = e^{J_p}-1 \ . 
</math> 
This leads to rewriting the partition function as 
: <math>
Z_p = \sum_\omega v^{\#\text{edges}(\omega)} q^{\#\text{clusters}(\omega)}
</math>
where the ''' FK clusters''' are the connected components of the union of closed segments <math>\cup_{(i,j)\in\omega}[i,j]</math>. This is proportional to the partition function of the random cluster model with the open edge probability <math>p=\frac{v}{1+v}=1-e^{-J_p}</math>. An advantage of the random cluster formulation is that <math>q</math> can be an arbitrary complex number, rather than a natural integer.

Alternatively, instead of FK clusters, the model can be formulated in terms of '''spin clusters''', using the identity 
: <math>
e^{J_p\delta(s_i,s_j)} = (1 - \delta(s_i,s_j)) + e^{J_p} \delta(s_i,s_j)\ . 
</math>
A spin cluster is the union of neighbouring FK clusters with the same color: two neighbouring spin clusters have different colors, while two neighbouring FK clusters are colored independently.

== Measure-theoretic description ==
The one dimensional Potts model may be expressed in terms of a [[subshift of finite type]], and thus gains access to all of the mathematical techniques associated with this formalism.  In particular, it can be solved exactly using the techniques of [[transfer operator]]s.  (However, [[Ernst Ising]] used combinatorial methods to solve the [[Ising model]], which is the "ancestor" of the Potts model, in his 1924 PhD thesis). This section develops the mathematical formalism, based on [[measure theory]], behind this solution.

While the example below is developed for the one-dimensional case, many of the arguments, and almost all of the notation, generalizes easily to any number of dimensions. Some of the formalism is also broad enough to handle related models, such as the [[XY model]], the [[Heisenberg model (classical)|Heisenberg model]] and the [[N-vector model]].

=== Topology of the space of states ===
Let ''Q'' = {1, ..., ''q''} be a finite set of symbols, and let
: <math>Q^\mathbf{Z}=\{ s=(\ldots,s_{-1},s_0,s_1,\ldots) :  s_k \in Q \; \forall k \in \mathbf{Z} \}</math>
be the set of all bi-infinite strings of values from the set ''Q''. This set is called a [[full shift]]. For defining the Potts model, either this whole space, or a certain subset of it, a [[subshift of finite type]], may be used. Shifts get this name because there exists a natural operator on this space, the [[shift operator]] τ : ''Q''<sup>'''Z'''</sup> → ''Q''<sup>'''Z'''</sup>, acting as
: <math>\tau (s)_k = s_{k+1}</math>

This set has a natural [[product topology]]; the [[base (topology)|base]] for this topology are the [[cylinder set]]s
: <math>C_m[\xi_0, \ldots, \xi_k]= \{s \in Q^\mathbf{Z} : s_m = \xi_0, \ldots ,s_{m+k} = \xi_k \}</math>
that is, the set of all possible strings where ''k''+1 spins match up exactly to a given, specific set of values ξ<sub>0</sub>, ..., ξ<sub>''k''</sub>. Explicit representations for the cylinder sets can be gotten by noting that the string of values corresponds to a [[p-adic number|''q''-adic number]], however the natural topology of the q-adic numbers is finer than the above product topology.

=== Interaction energy ===
The interaction between the spins is then given by a [[continuous function (topology)|continuous function]] ''V'' : ''Q''<sup>'''Z'''</sup> → '''R''' on this topology. ''Any'' continuous function will do; for example
: <math>V(s) = -J\delta(s_0,s_1)</math>
will be seen to describe the interaction between nearest neighbors. Of course, different functions give different interactions; so a function of ''s''<sub>0</sub>,  ''s''<sub>1</sub> and ''s''<sub>2</sub> will describe a next-nearest neighbor interaction.  A function ''V'' gives interaction energy between a set of spins; it is ''not'' the Hamiltonian, but is used to build it. The argument to the function ''V'' is an element ''s'' ∈ ''Q''<sup>'''Z'''</sup>, that is, an infinite string of spins. In the above example, the function ''V'' just picked out two spins out of the infinite string: the values ''s''<sub>0</sub> and ''s''<sub>1</sub>.  In general, the function ''V'' may depend on some or all of the spins; currently, only those that depend on a finite number are exactly solvable.

Define the function ''H<sub>n</sub>'' : ''Q''<sup>'''Z'''</sup> → '''R''' as
: <math>H_n(s)= \sum_{k=0}^n V(\tau^k s)</math>

This function can be seen to consist of two parts: the self-energy of a configuration [''s''<sub>0</sub>, ''s''<sub>1</sub>, ..., ''s<sub>n</sub>''] of spins, plus the interaction energy of this set and all the other spins in the lattice. The {{nowrap|''n'' → ∞}} limit of this function is the Hamiltonian of the system; for finite ''n'', these are sometimes called the '''finite state Hamiltonians'''.

=== Partition function and measure ===
The corresponding finite-state [[partition function (statistical mechanics)|partition function]] is given by
: <math>Z_n(V) = \sum_{s_0,\ldots,s_n \in Q} \exp(-\beta H_n(C_0[s_0,s_1,\ldots,s_n]))</math>
with ''C''<sub>0</sub> being the cylinder sets defined above. Here, β = 1/''kT'', where ''k'' is the [[Boltzmann constant]], and ''T'' is the [[temperature]]. It is very common in mathematical treatments to set β = 1, as it is easily regained by rescaling the interaction energy. This partition function is written as a function of the interaction ''V'' to emphasize that it is only a function of the interaction, and not of any specific configuration of spins.  The partition function, together with the Hamiltonian, are used to define a [[measure (mathematics)|measure]] on the Borel σ-algebra in the following way: The measure of a cylinder set, i.e. an element of the base, is given by
: <math>\mu (C_k[s_0,s_1,\ldots,s_n]) =  \frac{1}{Z_n(V)}  \exp(-\beta H_n (C_k[s_0,s_1,\ldots,s_n]))</math>

One can then extend by countable additivity to the full σ-algebra. This measure is a [[probability measure]]; it gives the likelihood of a given configuration occurring in the [[Configuration space (physics)|configuration space]] ''Q''<sup>'''Z'''</sup>. By endowing the configuration space with a probability measure built from a Hamiltonian in this way, the configuration space turns into a [[canonical ensemble]].

Most thermodynamic properties can be expressed directly in terms of the partition function. Thus, for example, the [[Helmholtz free energy]] is given by
: <math>A_n(V)=-kT \log Z_n(V)</math>

Another important related quantity is the [[topological pressure]], defined as
: <math>P(V) = \lim_{n\to\infty} \frac{1}{n} \log Z_n(V)</math>
which will show up as the logarithm of the leading eigenvalue of the [[transfer operator]] of the solution.

=== Free field solution ===
The simplest model is the model where there is no interaction at all, and so ''V'' = ''c'' and ''H<sub>n</sub>'' = ''c'' (with ''c'' constant and independent of any spin configuration).  The partition function becomes
: <math>Z_n(c) = e^{-c\beta} \sum_{s_0,\ldots,s_n \in Q} 1</math>

If all states are allowed, that is, the underlying set of states is given by a [[full shift]], then the sum may be trivially evaluated as
: <math>Z_n(c) = e^{-c\beta} q^{n+1}</math>

If neighboring spins are only allowed in certain specific configurations, then the state space is given by a [[subshift of finite type]]. The partition function may then be written as
: <math>Z_n(c) = e^{-c\beta} |\mbox{Fix}\, \tau^n| =  e^{-c\beta} \mbox{Tr} A^n</math>
where card is the [[cardinality]] or count of a set, and Fix is the set of [[Fixed point (mathematics)|fixed points]] of the iterated shift function:
: <math>\mbox{Fix}\, \tau^n =  \{ s \in Q^\mathbf{Z} : \tau^n s = s \}</math>

The ''q'' × ''q'' matrix ''A'' is the [[adjacency matrix]] specifying which neighboring spin values are allowed.

=== Interacting model ===
The simplest case of the interacting model is the [[Ising model]], where the spin can only take on one of two values, ''s<sub>n</sub>'' ∈ {{mset|−1, 1}} and only nearest neighbor spins interact. The interaction potential is given by
: <math>V(\sigma) = -J_p s_0 s_1\,</math>

This potential can be captured in a 2 × 2 matrix with matrix elements
: <math>M_{\sigma \sigma'} = \exp \left( \beta J_p \sigma \sigma' \right)</math>
with the index σ, σ′ ∈ {−1, 1}.  The partition function is then given by
: <math>Z_n(V) = \mbox{Tr}\, M^n</math>

The general solution for an arbitrary number of spins, and an arbitrary finite-range interaction, is given by the same general form. In this case, the precise expression for the matrix ''M'' is a bit more complex.

The goal of solving a model such as the Potts model is to give an exact [[closed-form expression]] for the partition function and an expression for the [[Gibbs state]]s or [[equilibrium state]]s in the limit of ''n'' → ∞, the [[thermodynamic limit]].

== Applications ==
=== Signal and image processing ===
The Potts model has applications in signal reconstruction. Assume that we are given noisy observation of a piecewise constant signal ''g'' in '''R'''<sup>''n''</sup>. To recover ''g'' from the noisy observation vector ''f'' in '''R'''<sup>''n''</sup>, one seeks a minimizer of  the corresponding inverse problem, the ''L<sup>p</sup>''-Potts functional ''P''<sub>γ</sub>(''u''), which is defined by
: <math> P_\gamma(u) = \gamma \| \nabla u \|_0 + \| u-f\|_p^p = \gamma \# \{ i : u_i \neq u_{i+1} \} + \sum_{i=1}^n  |u_i - f_i|^p</math>

The jump penalty <math>\| \nabla u \|_0</math> forces piecewise constant solutions and the data term <math>\| u-f\|_p^p</math> couples the minimizing candidate ''u'' to the data ''f''. The parameter γ > 0 controls the tradeoff between regularity and [[data fidelity]]. There are fast algorithms for the exact minimization of the ''L''<sup>1</sup> and the ''L''<sup>2</sup>-Potts functional.<ref>{{Cite journal |last1=Friedrich |first1=F. |last2=Kempe |first2=A. |last3=Liebscher |first3=V. |last4=Winkler |first4=G. |date=2008 |title=Complexity Penalized M-Estimation: Fast Computation |jstor=27594299 |journal=Journal of Computational and Graphical Statistics |volume=17 |issue=1 |pages=201–224 |doi=10.1198/106186008X285591 |s2cid=117951377 |issn=1061-8600}}</ref>

In image processing, the Potts functional is related to the segmentation problem.<ref>{{Cite journal |last1=Krähenbühl |first1=Philipp |last2=Koltun |first2=Vladlen |date=2011 |title=Efficient Inference in Fully Connected CRFs with Gaussian Edge Potentials |url=https://proceedings.neurips.cc/paper/2011/hash/beda24c1e1b46055dff2c39c98fd6fc1-Abstract.html |journal=Advances in Neural Information Processing Systems |publisher=Curran Associates, Inc. |volume=24|arxiv=1210.5644 }}</ref> However, in two dimensions the problem is NP-hard.<ref>{{Cite journal |last1=Boykov |first1=Y. |last2=Veksler |first2=O. |last3=Zabih |first3=R. |date=November 2001 |title=Fast approximate energy minimization via graph cuts |url=https://ieeexplore.ieee.org/document/969114 |journal=IEEE Transactions on Pattern Analysis and Machine Intelligence |volume=23 |issue=11 |pages=1222–1239 |doi=10.1109/34.969114 |issn=1939-3539}}</ref>

== See also ==

* [[Random cluster model]]
* [[Critical three-state Potts model]]
* [[Chiral Potts model]]
* [[Square-lattice Ising model]]
* [[Minimal model (physics)|Minimal models]]
* [[Z N model]]
* [[Cellular Potts model]]

== References ==
{{Reflist|refs=
<ref name="lyxt19">{{cite journal | last1=Li | first1=Zi-Qian | last2=Yang | first2=Li-Ping | last3=Xie | first3=Z. Y. | last4=Tu | first4=Hong-Hao | last5=Liao | first5=Hai-Jun | last6=Xiang | first6=T. | title=Critical properties of the two-dimensional <math>q</math>-state clock model | journal=[[Physical Review E]] | year=2020 | volume=101 | issue=6 | page=060105 | doi=10.1103/PhysRevE.101.060105 | pmid=32688489 | arxiv=1912.11416v3 | bibcode=2020PhRvE.101f0105L | s2cid=209460838 }}</ref>
}}

== External links ==
* {{cite web |title=Code for efficiently computing Tutte, Chromatic and Flow Polynomials |first1=Gary |last1=Haggard |first2=David J. |last2=Pearce |first3=Gordon |last3=Royle |url=http://www.mcs.vuw.ac.nz/~djp/tutte/ }}

{{Statistical mechanics topics}}
{{Stochastic processes}}

[[Category:Spin models]]
[[Category:Exactly solvable models]]
[[Category:Statistical mechanics]]
[[Category:Lattice models]]