StartseiteChemie & MaterialienHeisenberg-Spin-Modell — Quantenmagnetismus

🧲 Heisenberg-Spin-Modell — Quantenmagnetismus

Metropolis-Monte-Carlo-Simulation des klassischen Heisenberg-Modells: 3D-Einheitsspins auf einem 2D-Gitter ordnen sich mit sinkender Temperatur zu sichtbaren magnetischen Domänen an.

Chemie & Materialien2DFortgeschritten60 FPS⚡ Plasma
heisenberg-spin ↗ Eigenständig öffnen

Funktionsweise

Jeder Platz auf dem Gitter trägt einen klassischen Spin: einen Einheitsvektor (Sx, Sy, Sz), der überall auf einer Kugel zeigen kann, nicht nur „auf“ oder „ab“ wie im einfacheren Ising-Modell. Benachbarte Spins wechselwirken über die Austauschenergie E = -J·Σ Sᵢ·Sⱼ, summiert über nächste-Nachbar-Paare mit periodischen Randbedingungen (das Gitter schließt sich an den Rändern). Positives J begünstigt parallele Ausrichtung (ferromagnetisch); negatives J begünstigt antiparallele Ausrichtung (antiferromagnetisch).

Die Simulation verwendet den Metropolis-Monte-Carlo-Algorithmus: jedes Animationsbild führt einen vollständigen Durchlauf aus (im Durchschnitt ein Update-Versuch pro Platz). Bei jedem Versuch wird ein Platz zufällig ausgewählt, eine neue zufällige Einheitsvektor-Orientierung vorgeschlagen, und die resultierende Energieänderung ΔE aus seinen vier Gitternachbarn berechnet. Züge, die die Energie senken, werden immer akzeptiert; Züge, die sie erhöhen, werden mit Wahrscheinlichkeit exp(-ΔE / kBT) akzeptiert. Über viele Durchläufe treibt dies das Gitter zum thermischen Gleichgewicht bei der gewählten Temperatur T, wobei es sich sichtbar bei niedrigem T zu großen farbigen geordneten Domänen umorganisiert und sich bei hohem T in rauschartige Unordnung auflöst.

Energie: E = -J·Σ<i,j> Sᵢ·Sⱼ - h·Σᵢ Sᵢz
Metropolis-Akzeptanz: P = min(1, exp(-ΔE / kBT))
Magnetisierung: |M| = |(1/N)·Σᵢ Sᵢ|, 0 ≤ |M| ≤ 1
Farbzuordnung: Farbton = atan2(Sy, Sx), Helligkeit ∝ Sz

Häufig gestellte Fragen

Was ist das Heisenberg-Spin-Modell?

Das klassische Heisenberg-Modell stellt jedes magnetische Moment als kontinuierlichen 3D-Einheitsvektor dar, statt als feste Richtung. Benachbarte Spins wechselwirken über eine Austauschenergie E = -J·Sᵢ·Sⱼ, und das System entwickelt sich hin zu dem Orientierungsmuster, das die Gesamtenergie bei einer gegebenen Temperatur minimiert.

Wie unterscheidet sich das Heisenberg-Modell vom Ising-Modell?

Das Ising-Modell beschränkt jeden Spin auf zwei diskrete Zustände, auf oder ab. Das Heisenberg-Modell erlaubt Spins, in jede Richtung auf der Einheitskugel zu zeigen, was kontinuierliche Rotationsfreiheit ergibt — eine realistischere Beschreibung vieler realer magnetischer Materialien.

Was macht der Metropolis-Monte-Carlo-Algorithmus?

Bei jedem Schritt wird ein zufälliger Platz ausgewählt und eine neue zufällige Spin-Orientierung vorgeschlagen. Senkt sie die Energie, wird sie immer akzeptiert; erhöht sie die Energie, wird sie mit Wahrscheinlichkeit exp(-ΔE / kBT) akzeptiert. Dies tastet die Boltzmann-Verteilung ab, sodass das Gitter das thermische Gleichgewicht bei Temperatur T widerspiegelt.

Warum ist die Akzeptanzwahrscheinlichkeit exp(-ΔE/kBT)?

Dies ist der Boltzmann-Faktor aus der statistischen Mechanik: Die relative Wahrscheinlichkeit eines Zustands mit Energie E ist proportional zu exp(-E/kBT). Ihn als Akzeptanzregel zu verwenden, garantiert detailliertes Gleichgewicht, sodass simulierte Konfigurationen zum korrekten thermischen Ensemble konvergieren.

Was misst die Magnetisierung |M|?

Die Magnetisierung ist der Vektormittelwert aller Spins, M = (1/N)·Σ Sᵢ. Ihr Betrag reicht von 0 (zufällige Orientierungen, kein Nettomoment) bis 1 (alle Spins perfekt ausgerichtet) — der einfachste Ordnungsparameter für magnetische Ordnung.

Warum hat das 2D-Heisenberg-Modell keinen echten Phasenübergang bei endlicher Temperatur?

Das Mermin-Wagner-Theorem zeigt, dass kontinuierliche Symmetrien in 1D- oder 2D-kurzreichweitigen Systemen bei keiner T > 0 spontan gebrochen werden können — langwellige Spinwellen-Fluktuationen zerstören echte Fernordnung. Das Gitter zeigt dennoch starke Nahordnungs-Korrelationen und große geordnete Domänen bei niedrigem T; verwandte planare (XY-artige) Modelle zeigen stattdessen einen topologischen Berezinskii-Kosterlitz-Thouless-(BKT)-Übergang statt eines konventionellen Phasenübergangs.

Was stellt die Kopplungskonstante J dar?

J ist die Stärke der Austauschwechselwirkung zwischen benachbarten Momenten, die quantenmechanisch aus der Überlappung von Elektronenwellenfunktionen und dem Pauli-Ausschlussprinzip entsteht. Positives J (ferromagnetisch) begünstigt parallele Spins; negatives J (antiferromagnetisch) begünstigt ein schachbrettartiges antiparalleles Muster bei niedrigem T.

Was bewirkt die Steuerung des externen Feldes?

Das Aktivieren des Feldes fügt der Energie jedes Spins einen Zeeman-Term -h·Sᵢz hinzu, der die Ausrichtung entlang z begünstigt. Dies imitiert ein angelegtes Magnetfeld Bz; selbst ein bescheidenes Feld zieht das Gitter zu einer Netto-Magnetisierung außerhalb der Ebene, sichtbar als gleichmäßig helleres oder dunkleres Werden der Leinwand.

Wo wird das klassische Heisenberg-Modell in der realen Forschung eingesetzt?

Es ist ein Standard-Testfeld für Ferromagnete, Antiferromagnete und frustrierte Magnete und liegt Simulationen zugrunde, die im Spintronik-Gerätedesign, in der Forschung zu magnetischen Speichermedien und in Studien zum Quantenmagnetismus verwendet werden, wo klassische Monte-Carlo-Ergebnisse quantenmechanische Vielteilchenberechnungen leiten.

Was ist ein Monte-Carlo-Durchlauf?

Ein Durchlauf versucht im Durchschnitt ein Metropolis-Update pro Gitterplatz (N×N Versuche für ein N×N-Gitter). Viele Durchläufe lassen die Konfiguration von einem zufälligen Start hin zum Gleichgewicht relaxieren; Akzeptanzrate und Magnetisierung stabilisieren sich typischerweise nach einigen Dutzend Durchläufen.

Über diese Simulation

Dieser Simulator führt eine Live-Metropolis-Monte-Carlo-Simulation des klassischen Heisenberg-Modells aus, bei der jeder Gitterplatz einen frei rotierenden 3D-Spin trägt statt eines einfachen Auf/Ab-Bits. Jedes Bild führt einen vollständigen Gitterdurchlauf vorgeschlagener Spin-Umkehrungen aus, akzeptiert oder abgelehnt nach der Boltzmann-Regel, sodass Sie einen echten thermischen Relaxationsprozess beobachten statt einer vorgefertigten Animation. Bei niedriger Temperatur organisiert sich die farbenfrohe Farbton-und-Helligkeit-Textur zu breiten geordneten Domänen; erhöhen Sie die Temperatur, und dasselbe Gitter löst sich in schimmerndes, unkorreliertes Rauschen auf — ein anschauliches, interaktives Bild davon, wie magnetische Ordnung mit thermischer Unordnung konkurriert.

🔬 Was es zeigt

Ein 2D-Gitter klassischer 3D-Spins, das sich unter Metropolis-Monte-Carlo-Dynamik entwickelt, gefärbt nach Farbton (Winkel in der Ebene) und Helligkeit (Komponente Sz außerhalb der Ebene), zusammen mit einem Live-Diagramm der Magnetisierung |M| über die letzten Durchläufe, sodass Sie beobachten können, wie sich das System in Echtzeit einpendelt oder schwankt.

🎮 Anwendung

Ziehen Sie den Temperatur-Schieberegler, um den Ordnungsübergang zu durchqueren, kehren Sie das Vorzeichen von J um, um zwischen ferromagnetischer und antiferromagnetischer Kopplung zu wechseln, ändern Sie die Gittergröße für feinere oder gröbere Details, und schalten Sie das externe Feld um, um zu sehen, wie die Spins zur Ausrichtung entlang z gezogen werden.

💡 Wussten Sie schon?

Anders als das 3D-Heisenberg-Modell hat diese 2D-Version bei keiner Temperatur ungleich null echte magnetische Fernordnung (Mermin-Wagner-Theorem) — dennoch zeigt es bei niedrigem T bemerkenswerte Nahordnung und langsame, glasartig anmutende Domänendynamik, eine Physik, die eng mit dem Berezinskii-Kosterlitz-Thouless-Übergang in planaren Spinsystemen verwandt ist.

Häufig gestellte Fragen

Warum sieht das Gitter bei hoher Temperatur wie farbiges Rauschen aus?

Bei hohem T dominieren thermische Fluktuationen die Austauschkopplung J, sodass der Metropolis-Algorithmus fast jede vorgeschlagene Neuausrichtung akzeptiert. Jeder Spin zeigt schließlich in eine im Wesentlichen zufällige Richtung, unkorreliert mit seinen Nachbarn, was die Farbton/Helligkeit-Farbkarte als verrauschtes, unkorreliertes Sprenkelmuster darstellt.

Warum erscheinen bei niedriger Temperatur große farbige Domänen?

Bei niedrigem T dominiert die Austauschenergie, und falsch ausgerichtete Nachbarn kosten zu viel Energie, als dass die Metropolis-Regel sie leicht akzeptieren würde. Spins richten sich mit ihren Nachbarn in Flecken aus, die über aufeinanderfolgende Durchläufe wachsen und verschmelzen, was die breiten, sanft variierenden farbigen Bereiche erzeugt, die Sie sehen.

Was passiert, wenn ich J negativ setze?

Negatives J macht antiparallele benachbarte Spins energetisch günstiger als parallele. Bei niedriger Temperatur entwickelt das Gitter ein schachbrettartiges antiferromagnetisches Muster wechselnder Orientierungen statt gleichmäßiger ferromagnetischer Domänen.

Warum schwankt das |M|-Diagramm selbst bei niedrigem T noch?

Da es sich um ein endliches Gitter mit echtem thermischem Rauschen handelt, kommt die Magnetisierung nie perfekt zur Ruhe — sie schwankt um einen Gleichgewichtswert, dessen Größe von Temperatur und Gittergröße abhängt. Größere Gitter bei niedrigem T zeigen kleinere relative Fluktuationen, näher am idealisierten thermodynamischen Grenzfall.

Was sagt mir die Statistik zur Akzeptanzrate?

Die Akzeptanzrate ist der Anteil der während des letzten Durchlaufs akzeptierten vorgeschlagenen Spin-Neuausrichtungen. Sie ist typischerweise hoch bei hoher Temperatur (fast jeder Zug wird akzeptiert) und niedriger bei niedriger Temperatur (nur Züge, die die Energie kaum erhöhen oder senken, sind erfolgreich), was sie zu einer schnellen Diagnose dafür macht, wie „heiß“ oder „eingefroren“ das System gerade ist.

Warum Einheitsvektoren statt Winkel verwenden?

Die Darstellung jedes Spins als 3D-Einheitsvektor (Sx, Sy, Sz) macht das Skalarprodukt Sᵢ·Sⱼ — und damit die Austauschenergie — zu einer einfachen, numerisch stabilen Berechnung und lässt das Modell auf natürliche Weise die Neigung außerhalb der Ebene (hier als Helligkeit erfasst) zusammen mit der Rotation in der Ebene (als Farbton erfasst) einschließen, was ein einzelner Winkel in der Ebene nicht darstellen könnte.

⚙ Unter der Haube

Metropolis-Monte-Carlo-Simulation des klassischen Heisenberg-Modells: 3D-Einheitsspins auf einem 2D-Gitter ordnen sich mit sinkender Temperatur zu sichtbaren magnetischen Domänen an.

Heisenberg ModelMonte CarloMagnetismMetropolis AlgorithmCanvas 2D

2D · HTML5 Canvas 2D · Ziel 60 FPS · läuft vollständig clientseitig, keine Installation nötig

Was haben Sie gefunden?

Schritte zur Reproduktion hinzufügen (optional)