Warum fallen Spinkorrelationsfunktionen in Ising-Modellen exponentiell unterhalb der kritischen Temperatur ab?

Ich versuche, das 2D-Ising-Modell besser zu verstehen, insbesondere das Verhalten der Korrelationsfunktionen zwischen Abstandsspins R .

Ich habe eine Reihe von erklärenden Texten gefunden, die darauf hinzudeuten scheinen, dass dies sowohl über als auch unter der kritischen Temperatur liegt T C , die Korrelationsfunktion C ( R ) bei einer bestimmten Korrelationslänge exponentiell abfällt ξ , und dies hilft uns, typische Domänengrößen zu bestimmen. Bei T C , divergiert die Korrelationslänge gegen unendlich. Hier ist ein Bild von ( http://math.arizona.edu/~tgk/541/chap1.pdf ), um zu veranschaulichen, was ich meine ( β ist natürlich inverse Temperatur).

Geben Sie hier die Bildbeschreibung ein

  1. T > T C : Das macht für mich Sinn - benachbarte Spins sind praktisch unabhängig, so dass Domänen winzig sind und die Korrelationslänge mit zunehmender Temperatur gegen Null geht.

  2. T < T C : Das ergibt für mich keinen intuitiven Sinn - ich hatte den Eindruck, dass unten T C die Domänen waren groß genug, um sicherzustellen, dass eine spontane Magnetisierung beobachtet wurde. Das Obige weist jedoch darauf hin, dass die typische Korrelationslänge ( ξ ) mit steigender Temperatur gegen Null tendiert - und damit Domänen schrumpfen? Ich hätte gedacht:

    • Die Korrelationsfunktion C R sollte überhaupt nicht exponentiell abfallen, sondern konstant bleiben (wie in Kapitel 1, Seite 6 von Paul Fendleys Entwurf "Modern Statistical Mechanics" angegeben)
    • Wenn sie exponentiell abfallen, würden sie auf einen Wert > 0 abfallen (wie auf Seite 216 von Sethnas "Entropy, Order Parameters and Complexity" angegeben).
    • Wenn sie exponentiell und auf 0 abfallen, dann sicherlich als T 0 , die Korrelationslänge ξ sollte gegen unendlich tendieren, wie bei T = 0 wir wissen, dass alle Spins perfekt mit beliebig weit entfernten Spins korrelieren müssen?
  3. T = T C : Wenn die Korrelationslänge unendlich ist, würden wir sicherlich eine perfekte Korrelation und Gesamtmagnetisierung in der Region nahe der kritischen Temperatur sehen? Stattdessen nimmt bei den kritischen Temperaturreferenzen die Korrelationsfunktion die Form an C R R λ - die Intuition, wie dieser Sprung gemacht wird, ist mir nicht klar.

Offensichtlich habe ich grundlegend missverstanden, wie die Korrelationsfunktion und die kritische Länge funktionieren, wie sie sich auf Domänengrößen beziehen, oder beides. Ich würde mich sehr freuen, wenn jemand sagen könnte, wo ich falsch verstanden habe.

Als Referenz habe ich einen wirtschaftswissenschaftlichen, keinen physikalischen Hintergrund, muss aber die Intuition dieser Modelle für meine Doktorarbeit zur Meinungsdynamik verstehen. Mein primärer Referenztext ist "Dynamical Processes on Complex Systems" von Barrat et al.

Antworten (1)

Beachten Sie zunächst, dass, wie Sie sagen, die 2-Punkt-Funktion σ ich σ J nicht gegen Null tendiert J ich Wenn T < T C ; nämlich,

lim J ich σ ich σ J + = ( M ( T ) ) 2 ,
Wo M ( T ) = σ 0 + ist die spontane Magnetisierung (die + hochgestellt zeigt an, dass die Erwartung in der genommen wird + Zustand).

Also, wenn man sagt, dass Korrelationen exponentiell abfallen, wenn T < T C , spricht man eigentlich von den abgeschnittenen Korrelationsfunktionen. Die abgeschnittene 2-Punkt-Funktion ist definiert als

σ ich ; σ J + = ( σ ich σ ich + ) ( σ J σ J + ) + = σ ich σ J + σ ich + σ J + .
Sie misst, wie korreliert die Schwankungen an ich Und J Sind. (Probabilistisch gesehen ist dies einfach die Kovarianz zwischen den Zufallsvariablen σ ich Und σ J .) Die abgeschnittene 2-Punkt-Funktion konvergiert gegen 0 bei allen Temperaturen T . Dies ist eine allgemeine Tatsache, die in den reinen Phasen aller Modelle bei jeder Temperatur zutrifft.

Nun zu Ihrer Frage: Die abgeschnittenen Korrelationen fallen im 2D-Ising-Modell für alle tatsächlich exponentiell schnell ab T T C .

Wenn T < T C , sollten Sie sich das wie folgt vorstellen: Bei niedrigen Temperaturen nehmen die Spins typischerweise die gleichen Werte an (z + 1 im + Phase), mit nur seltenen Schwankungen. Eine nützliche Möglichkeit, diese Schwankungen zu betrachten, besteht darin, die Konfigurationen aus geometrischer Sicht zu betrachten: Zeichnen Sie ein Einheitslängensegment, das jedes Paar von nächstgelegenen Eckpunkten trennt, an denen die Spins entgegengesetzte Werte annehmen. Die Vereinigung dieser Segmente bildet die Peierls-Konturen der Konfiguration.

Geben Sie hier die Bildbeschreibung ein

Es ist leicht nachzuprüfen, dass die mit jeder solchen Kontur verbundenen Energiekosten proportional zu ihrer Länge sind. Beachten Sie auch, dass die Konturen eine vollständige Beschreibung der Konfiguration liefern (wenn Sie wissen, dass Sie sich in der + Phase).

Was hat das nun mit dem exponentiellen Abfall von Korrelationen zu tun? Wie oben erwähnt, die abgeschnittene 2-Punkt-Funktion σ ich ; σ J + misst, wie korreliert die Schwankungen an ich Und J Sind. Was ist das Ereignis, das zu einem gleichzeitigen Umdrehen beider Spin führen wird ich Und J ? Es ist nicht schwer, sich davon zu überzeugen, dass dies geschehen sollte, wenn eine große Kontur gleichzeitig die Eckpunkte umgibt ich Und J . Mit Hilfe von Korrelationsungleichungen kann man das nämlich überprüfen

0 σ ich ; σ J + P R Ö B ( es gibt eine Kontur, die beide umgibt  ich  Und  J ) ,
bei dem die + hochgestellt zeigt an, dass die Berechnung in erfolgt + Phase. Aber das Ereignis auf der rechten Seite hat eine Wahrscheinlichkeit, die in der Ferne exponentiell klein ist J ich .

All dies kann rigoros gemacht werden. Bei ausreichend niedrigen Temperaturen können Sie Cluster-Expansionstechniken verwenden (dies funktioniert in jeder Dimension D 2 ). Eine ausführliche Argumentation finden Sie beispielsweise in Theorem 5.27 in diesem Buch . Im Maß 2 , kann dies auch störungsfrei bewiesen werden, entweder durch explizite Berechnungen oder indem die Kosten einer großen Kontur mit der Oberflächenspannung im Modell in Beziehung gesetzt werden; siehe zum Beispiel dieses Papier .

Als letzte Bemerkung: Das obige Argument legt nahe, und Sie können dies rigoros machen (siehe zum Beispiel den oben erwähnten Beweis der Clustererweiterung), dass die Korrelationslänge gegen Null tendiert, wenn T 0 . Dies liegt daran, dass es äußerst unwahrscheinlich wird, dass eine Kontur zwei entfernte Eckpunkte umgibt, und Fluktuationen bei sehr niedrigen Temperaturen zu im Wesentlichen rein lokalen Ereignissen werden und somit (ungefähr) unabhängig an verschiedenen Eckpunkten auftreten.

Vielen Dank für eine sehr klare Antwort. Erstens denke ich, dass eine große Verwirrung entstanden ist, weil die < . > + Notation fehlte bei meiner Lektüre - als solche nahm ich an < σ ich > Null sein muss (durch Symmetrie) und implizit < σ ich σ J > ging auch gegen null. (Würde das Argument immer noch für freie Grenzen gelten?) Zweitens ist das Konturen-Argument sehr nett; Ich folge jedoch nicht ganz der letzten Zeile - sicherlich wird mit T -> 0 die Wahrscheinlichkeit einer Kontur, die zwei entfernte Spins umgibt, größer (wenn die Domänen wachsen), und daher sollte nach Ihrem anfänglichen Argument die Korrelationslänge zunehmen?
@Suhdohnimh: Ja, die abgeschnittene 2-Punkt-Funktion stimmt mit der Standard-2-Punkt-Funktion überein, wenn freie Randbedingungen verwendet werden, sodass sie nicht gegen 0 konvergiert. Tatsächlich tritt eine Konvergenz von Korrelationen gegen Null nur in den sogenannten Extremalzuständen auf (beim 2d-Ising-Modell bedeutet dies in den Phasen + und −). Nun zu Ihrer zweiten Frage: Die Wahrscheinlichkeit, eine bestimmte Kontur γ zu sehen, ist nach oben begrenzt durch e−2|γ|/T, wobei |γ| bezeichnet die Länge von γ. Beachten Sie, dass dies bei T↓0 gegen 0 konvergiert. [fortgesetzt werden...]
[...] Das ist leicht verständlich. Wenn T≪1 , nehmen fast alle Spins den Wert 1 an (im + Zustand), weil es so kostspielig wird, dass ein Spin nicht mit seinen Nachbarn übereinstimmt. Natürlich werden solche Schwankungen noch auftreten, aber sie werden sehr selten sein, und je größer sie sind, desto seltener werden sie. Tatsächlich ist diese Beobachtung der Kern des Peierls-Arguments, das verwendet wird, um die Existenz eines Phasenübergangs im Ising-Modell in den Dimensionen 2 und mehr zu beweisen. Ich empfehle dringend, dass Sie es sich ansehen (es ist nicht schwierig). [fortgesetzt werden...]
[...] Sie finden eine detaillierte Beschreibung des Peierls-Arguments in dem Buch, auf das ich mich in meiner Antwort beziehe, aber informellere finden Sie in vielen Lehrbüchern der statistischen Mechanik, wie denen von Huang, von Ma usw.
Nochmals vielen Dank, ich habe mehr gelesen und einige Ihrer Vorträge angeschaut, um ein besseres Verständnis zu bekommen. Ist diese Intuition richtig: Bei T>>Tc korrelieren die Spinfluktuationen aufgrund fehlender Ausrichtung zwischen lokalen Spins nicht. Wenn T -> Tc, verbessert sich die lokale Ausrichtung, und dies erhöht die Korrlänge. Bei T = 0 korrelieren die Schwankungen trivialerweise nicht, da keine Schwankungen auftreten. Bei 0<T<<Tc beeinflussen irgendwelche Schwankungen in einem gegebenen Spin weit entfernte Spins nicht, da weit entfernte Spins unter starkem Druck stehen, sich mit den unendlich anderen geordneten Spins global auszurichten. [TBC]
[...] Bei T->Tc nimmt der Druck zur globalen Ausrichtung ab und der Einfluss von Schwankungen bestimmter Spins nimmt zu, wodurch die Korrelationslänge verbessert wird. Bei T=Tc entspricht der globale Konformitätsdruck genau dem Einfluss bestimmter Schwankungen, und es wird kritisches Verhalten beobachtet. [TBC]
[...] Visuell sehen wir, was von High T -> Tc passiert, wir sehen Ordnungsbereiche im Gitter inmitten eines riesigen Feldes von Unordnung wachsen. Von Low T -> Tc sehen wir Bereiche der Unordnung ("Löcher") in einem Ordnungsfeld wachsen. Die Korrelationslänge misst die typische Größe von Domänen bzw. "Löchern". Bei T=Tc befinden wir uns genau in der Mitte zwischen den Feldern voller Ordnung und Unordnung.
Ja, deine Beschreibung ist richtig. Eine schöne alternative Möglichkeit, die Wirkung von Energie (die zu Ausrichtung führt) und Entropie (thermische Fluktuationen, die zu Unordnung führen) zu visualisieren, ist die zufällige Clusterdarstellung des Modells (es gibt eine kurze Diskussion im Buch und Verweise auf ausführlichere Behandlungen). [...]
[...] Bei letzterem wird eine Konfiguration in zwei Schritten erzeugt: Zunächst werden nach einem gut gewählten Wahrscheinlichkeitsmaß zufällig nächste Nachbarspins verbunden. Dann legt man denselben Spinwert fest ( + oder ) für alle Knoten, die verbunden sind. Der Wert für jeden Cluster wird unabhängig voneinander zufällig gewählt, mit Ausnahme des einzigartigen unendlichen Clusters (falls vorhanden), dessen Spins alle den Wert 1 annehmen müssen (in der + Zustand). Die unendliche Komponente induziert somit die globale Ordnung, während endliche Cluster die lokale Ordnung beschreiben. Es gibt genau dann einen unendlichen Cluster T < T C .
Eine wunderbare Erklärung, vielen Dank. PS Deine Vorlesungen sind großartig und deine Tafelschrift ist phänomenal.
@Suhdohnimh: Ich freue mich, dass ich helfen konnte, und vielen Dank für die netten Worte :) .