Potts-Modell

Potts-Modell

Version vom 9. März 2017, 10:24 Uhr von 79.201.131.111 (Diskussion) (Einordnung von MC-Simulationen und RNG-Theorie; genauerer Link zu DP Landau)
(Unterschied) ← Nächstältere Version | Aktuelle Version (Unterschied) | Nächstjüngere Version → (Unterschied)

Das Potts-Modell ist ein mathematisches Modell, welches das in der Statistischen Physik häufig verwendete Ising-Modell verallgemeinert. Auf einem Gitter befinden sich statt Spins mit nur zwei Zuständen, wie im Ising-Modell, Variablen mit q verschiedenen Zuständen. Im einfachsten Fall beschränkt sich die Wechselwirkung auf benachbarte Gitterplätze. Anwendung findet es unter anderem, außer in der Statistische Physik (insbesondere beim Studium von Phasenübergängen), auch in der Informatik (Signalverarbeitung) und Biologie (neuronale Netze). Das Modell wurde benannt nach Renfrey Potts, welcher das Modell 1951 in seiner Dissertation definierte. Einen Spezialfall behandelten vorher schon Julius Ashkin und Edward Teller 1943. Einen Überblick zu Geschichte und Analyse des Modells gibt ein Übersichtsartikel von Fa-Yueh Wu aus dem Jahr 1982.

Definition

Das Potts-Modell besteht aus einem d-dimensionalen Gittergraphen $ L(V,E) $, zum Beispiel einem zweidimensionalen Rechteckgitter, einer Menge von Knotenbelegungen $ S'^{|V|} $ (Knotenkonfigurationen) und einem Hamiltonoperator $ H' $ auf dieser Menge. Jeder Knoten wird mit einem Element aus der Menge

$ S'=\{\theta _{s}={\frac {2\pi s}{q}},0\leq s\leq q\}\quad q\in \{2,3,\ldots \} $

belegt. Diese können als Punkte auf dem 2-dimensionalen Einheitskreis interpretiert werden und sind die Richtungen, die die Spins auf den jeweiligen Gitterpunkten annehmen können. Der Hamiltonoperator ist im planaren Pottsmodell (auch Vektor-Potts-Modell oder Uhren-Modell, clock model) gegeben durch

$ H'(\sigma )=J_{1}\sum _{\langle i,j\rangle \in E}\cos(\theta _{s_{i}}-\theta _{s_{j}}),\quad \sigma =(\theta _{s_{1}},\theta _{s_{2}},\dots )\in S'^{|V|}. $

Summiert wird über alle benachbarten Knoten $ \langle i,j\rangle $ für $ i,j\in V $. Die Kopplungskonstante $ J_{1} $ beschreibt die Wechselwirkung zwischen Spins auf den benachbarten Knoten.

Alternativ zu dem gerade beschriebenen planaren Potts-Modell gibt es das Standard-Potts-Modell (oder einfach: Potts-Modell). Dabei werden die Knoten mit Elementen aus der Menge $ S=\{0,1,\dots ,q\} $ belegt. Der Hamiltonoperator ist gegeben durch

$ H(\sigma )=-J\sum _{\langle i,j\rangle \in E}\delta (s_{i},s_{j}),\quad \sigma =(s_{1},s_{2},\dots )\in S^{|V|}, $

wobei $ \delta $ das Kronecker-Delta ist. Das heißt, falls zwei benachbarte Knoten verschiedene Werte der Spins besitzen, verschwindet der entsprechende Summand. Das negative Vorzeichen von $ J $ ist eine Konvention, motiviert vom Ising-Modell. Das Standard-Potts-Modell ist ferromagnetisch für $ J>0 $ und antiferromagnetisch für $ J<0 $.

Die Unterscheidung von planarem und Standard-Potts-Modell stammt von Cyril Domb (1974)[1].

Verhältnis zu anderen statistischen Modellen

Allgemeine Version

Auf dem Gittergraphen $ L(V,E) $ mit der Menge der Knotenbelegungen $ \{1,2,\dots ,q-1\}^{|V|} $ kann eine allgemeinere Version des Potts-Modells definiert werden:

$ H(\sigma )=-\sum _{\langle i,j\rangle \in E}J_{ij}\delta (s_{i},s_{j})-\beta ^{-1}\sum _{i\in \{1,2,\dots \,|V|\}}h_{i}s_{i},\quad \sigma =(s_{1},s_{2},\dots )\in S^{|V|}. $

Im Unterschied zum ursprünglichen Modell variiert die Wechselwirkung zwischen den benachbarten Knoten. Außerdem kann ein äußeres Feld ergänzt werden:

$ -\beta ^{-1}\sum _{i\in \{1,2,\dots ,|V|\}}h_{i}s_{i} $

Hierbei ist wie üblich $ \beta ={\frac {1}{k_{b}T}} $.

Die Wechselwirkungen, $ J_{ij} $, müssen nicht auf nächstbenachbarte Gitterplätze beschränkt werden. In verdünnten Potts-Modellen gibt es freie Gitterplätze (Gitter-Gas) oder auch Wechselwirkungen verschiedener Stärke. Durch geeignete Randbedingungen können interessante Effekte, wie z.B. Benetzung [2] oder Grenzflächenadsorption [3], induziert werden.

Das Ising-Modell als Spezialfall

Setzt man $ q=2 $, folgt aus dem Potts-Modell das Ising-Modell.

Das XY-Modell als Spezialfall

Für $ q\to \infty $ erhält man das XY-Modell, welches wiederum als Spezialfall des N-Vektor-Modells mit $ N=2 $ verstanden werden kann. Betrachtet man das planare Potts-Modell, so ist der Zustandsraum der Spins keine endliche Teilmenge des Einheitskreises, sondern der ganze 2-dimensionale Einheitskreis.

Ashkin-Teller-Modell

Das Ashkin-Teller-Modell ist das planare Potts-Modell mit q=4 Zuständen.

Sonstige

Es gibt auch Verbindungen zum Heisenberg-Modell, N-Vektor-Modell, (ice-rule) Vertex-Modellen und zur Perkolationstheorie (zuerst von P.W. Kasteleyn und C.M. Fortuin 1969 [4] für Bond-Perkolation, später auch für Site-Perkolation).

Das Kirchhoffsche Gesetz für Netzwerke aus linearen Widerständen ergibt sich als q=0 Grenzwert des Potts-Modells (Kasteleyn, Fortuin 1972).

Diskussion

Potts betrachtete das planare Modell und konnte ähnlich wie beim Ising-Modell mit Kramers-Wannier-Dualität den kritischen Punkt für das Rechteckgitter und q=2, 3, 4 bestimmen. Am Ende seiner Arbeit gab er den kritischen Punkt des Standart-Potts-Modells für alle q.

Das planare und das Standard-Modell sind identisch für q=2 (Ising-Modell) mit $ J=2J_{1} $ und q=3 mit $ J={\frac {3}{2}}J_{1} $. Außerdem ist das q=4 planare Modell für beliebige Gitter auf das q=2 Modell reduzierbar. Für $ q>4 $ gibt es dagegen keine offensichtlichen Bezüge zwischen dem planaren und Standard-Modell.

Auf einem zweidimensionalen Gitter hat das Potts-Modell mit $ J>0 $ einen Phasenübergang erster Ordnung für $ q>4 $ und ansonsten einen kontinuierlicher Phasenübergang (2. Ordnung) wie beim Isingmodell (q=2) (Rodney Baxter 1973, 1978). Baxter benutzte dabei die Identifizierung des zweidimensionalen Potts-Modells mit dem ice-rule Vertexmodell durch Temperley und Elliott Lieb (1971 für ein Gitter aus Quadraten).

Das eindimensionale Potts-Modell ist exakt lösbar (mit Hilfe der Transfer-Matrix-Methode) und ebenso das zweidimensionale Modell mit Wechselwirkungen zwischen benachbarten Gitterplätzen. Im Allgemeinen liefern insbesondere Monte-Carlo-Simulationen [5] und die Renormierungsgruppentheorie verläßliche Ergebnisse.

Potts-Maß

Mit der Hamilton-Funktion wie oben

$ H(\sigma )=-J\sum _{\langle i,j\rangle \in E}\delta (s_{i},s_{j}),\quad \sigma =(s_{1},s_{2},\dots )\in S^{|V|}. $

und der üblichen Definition der Zustandssumme $ Z $

$ Z(\beta )=\sum _{\sigma \in S^{|V|}}e^{-\beta H(\sigma )},\quad \beta \in ]0,\infty [. $

kann man das Potts-Maß $ \pi $ definieren, das als Wahrscheinlichkeitsmaß zu den Boltzmannverteilungen gehört:

$ \pi (\sigma )=Z(\beta )^{-1}e^{-\beta H(\sigma )},\quad \sigma \in S^{|V|}. $

Die Freie Energie ist wie üblich:

$ F(\beta )=-\beta ^{-1}\log(Z(\beta )),\quad \beta \in ]0,\infty [. $

Literatur

Einzelnachweise

  1. C.Domb, J. Phys. A, Band 7, 1974, S. 1335
  2. S.Dietrich, in Phase transitions and critical phenomena (Hrsg. C. Domb und J.L. Lebowitz), Band 12, 1988
  3. W.Selke, W.Pesch, Z. Phys. B, Band 47, 1982, S. 335
  4. P.W. Kasteleyn, C.M. Fortuin, J. Phys. Soc. Japan, Band 26 (Suppl.),1969, S.11
  5. D.P.Landau, K.Binder, A Guide to Monte Carlo Simulations in Statistical Physics, 2014