Strona główna Chemia Heisenberg Spin Model — Quantum Magnetism

🧲 Heisenberg Spin Model — Quantum Magnetism

Metropolis Monte Carlo simulation of the classical Heisenberg model: 3D unit spins on a 2D lattice order into visible magnetic domains as temperature drops.

Chemia2DZaawansowany60 FPS
heisenberg-spin ↗ Otwórz osobno
Interfejs samej symulacji jest w języku angielskim.

Podobne symulacje

Jak to działa

Każde miejsce na sieci przechowuje spin klasyczny: wektor jednostkowy (Sx, Sy, Sz), który może wskazywać dowolny kierunek na sferze, a nie tylko „góra” lub „dół” jak w prostszym modelu Isinga. Sąsiednie spiny oddziałują poprzez energię wymiany E = -J·Σ Sᵢ·Sⱼ, sumowaną po parach najbliższych sąsiadów z periodycznymi warunkami brzegowymi (sieć zawija się na krawędziach). Dodatnie J faworyzuje ustawienie równoległe (ferromagnetyczne); ujemne J faworyzuje ustawienie antyrównoległe (antyferromagnetyczne).

Symulacja wykorzystuje algorytm Monte Carlo Metropolisa: każda klatka animacji wykonuje pełne przemiatanie (średnio jedna próba aktualizacji na miejsce). Przy każdej próbie losowo wybierane jest miejsce, proponowana jest nowa losowa orientacja wektora jednostkowego, a wynikowa zmiana energii ΔE jest obliczana na podstawie jego czterech sąsiadów sieciowych. Ruchy obniżające energię są zawsze akceptowane; ruchy podnoszące energię są akceptowane z prawdopodobieństwem exp(-ΔE / kBT). Po wielu przemiataniach prowadzi to sieć do równowagi termicznej w wybranej temperaturze T, widocznie reorganizując się w duże, kolorowe uporządkowane domeny przy niskim T i rozpuszczając się w nieuporządkowany szum przy wysokim T.

Energia: E = -J·Σ<i,j> Sᵢ·Sⱼ - h·Σᵢ Sᵢz
Prawdopodobieństwo akceptacji Metropolisa: P = min(1, exp(-ΔE / kBT))
Magnetyzacja: |M| = |(1/N)·Σᵢ Sᵢ|, 0 ≤ |M| ≤ 1
Mapowanie kolorów: barwa = atan2(Sy, Sx), jasność ∝ Sz

Najczęściej zadawane pytania

Czym jest model spinowy Heisenberga?

Klasyczny model Heisenberga reprezentuje każdy moment magnetyczny jako ciągły trójwymiarowy wektor jednostkowy, a nie ustalony kierunek. Sąsiednie spiny oddziałują poprzez energię wymiany E = -J·Sᵢ·Sⱼ, a układ ewoluuje w kierunku wzoru orientacji minimalizującego całkowitą energię w danej temperaturze.

Czym model Heisenberga różni się od modelu Isinga?

Model Isinga ogranicza każdy spin do dwóch dyskretnych stanów, góra lub dół. Model Heisenberga pozwala spinom wskazywać dowolny kierunek na sferze jednostkowej, dając ciągłą swobodę obrotową — bardziej realistyczny opis wielu prawdziwych materiałów magnetycznych.

Co robi algorytm Monte Carlo Metropolisa?

W każdym kroku wybierane jest losowe miejsce i proponowana jest nowa losowa orientacja spinu. Jeśli obniża energię, jest zawsze akceptowana; jeśli podnosi energię, jest akceptowana z prawdopodobieństwem exp(-ΔE / kBT). Próbkuje to rozkład Boltzmanna, więc sieć odzwierciedla równowagę termiczną w temperaturze T.

Dlaczego prawdopodobieństwo akceptacji to exp(-ΔE/kBT)?

To czynnik Boltzmanna z mechaniki statystycznej: względne prawdopodobieństwo stanu o energii E jest proporcjonalne do exp(-E/kBT). Użycie go jako reguły akceptacji gwarantuje szczegółową równowagę, więc symulowane konfiguracje zbiegają do poprawnego zespołu termicznego.

Co mierzy magnetyzacja |M|?

Magnetyzacja to wektorowa średnia wszystkich spinów, M = (1/N)·Σ Sᵢ. Jej wielkość mieści się w zakresie od 0 (losowe orientacje, brak momentu netto) do 1 (wszystkie spiny idealnie ustawione) — najprostszy parametr porządku dla uporządkowania magnetycznego.

Dlaczego dwuwymiarowy model Heisenberga nie ma prawdziwego przejścia fazowego w skończonej temperaturze?

Twierdzenie Mermina-Wagnera pokazuje, że symetrie ciągłe nie mogą zostać spontanicznie złamane w jedno- lub dwuwymiarowych układach krótkozasięgowych przy żadnym T > 0 — fluktuacje fal spinowych o długiej długości fali niszczą prawdziwy porządek dalekiego zasięgu. Sieć nadal wykazuje silne korelacje krótkiego zasięgu i duże uporządkowane domeny przy niskim T; pokrewne modele płaskie (typu XY) wykazują topologiczne przejście Berezinskiego-Kosterlitza-Thoulessa (BKT) zamiast konwencjonalnego przejścia.

Co reprezentuje stała sprzężenia J?

J to siła oddziaływania wymiennego między sąsiednimi momentami, wynikająca mechanicznie kwantowo z nakładania się funkcji falowych elektronów i zasady wykluczania Pauliego. Dodatnie J (ferromagnetyczne) faworyzuje spiny równoległe; ujemne J (antyferromagnetyczne) faworyzuje wzór szachownicy antyrównoległej przy niskim T.

Co robi kontrola pola zewnętrznego?

Włączenie pola dodaje wyraz Zeemana -h·Sᵢz do energii każdego spinu, faworyzując ustawienie wzdłuż z. Naśladuje to przyłożone pole magnetyczne Bz; nawet umiarkowane pole ciągnie sieć w stronę netto magnetyzacji poza płaszczyzną, widocznej jako jednolite rozjaśnienie lub przyciemnienie płótna.

Gdzie klasyczny model Heisenberga jest wykorzystywany w prawdziwych badaniach?

To standardowy poligon testowy dla ferromagnetyków, antyferromagnetyków i magnetyków sfrustrowanych, leżący u podstaw symulacji wykorzystywanych w projektowaniu urządzeń spintronicznych, badaniach nośników zapisu magnetycznego oraz badaniach magnetyzmu kwantowego, gdzie klasyczne Monte Carlo prowadzi obliczenia kwantowe wielu ciał.

Czym jest przemiatanie Monte Carlo?

Jedno przemiatanie próbuje średnio jedną aktualizację Metropolisa na miejsce sieci (N×N prób dla sieci N×N). Wiele przemiatań pozwala konfiguracji zrelaksować się od losowego startu w stronę równowagi; szybkość akceptacji i magnetyzacja zwykle stabilizują się po kilkudziesięciu przemiataniach.

O tej symulacji

Ten symulator uruchamia żywą symulację Monte Carlo Metropolisa klasycznego modelu Heisenberga, gdzie każde miejsce sieci przechowuje swobodnie obracający się trójwymiarowy spin zamiast prostego bitu góra/dół. Każda klatka wykonuje pełne przemiatanie sieci z proponowanymi zmianami spinów, akceptowanymi lub odrzucanymi zgodnie z regułą Boltzmanna, więc obserwujesz prawdziwy proces relaksacji termicznej, a nie wcześniej przygotowaną animację. W niskiej temperaturze kolorowa tekstura barwy i jasności organizuje się w szerokie uporządkowane domeny; podkręć temperaturę, a ta sama sieć rozpuszcza się w migoczący, nieskorelowany szum — żywy, interaktywny obraz tego, jak porządek magnetyczny rywalizuje z nieporządkiem termicznym.

🔬 Co pokazuje

Dwuwymiarową sieć klasycznych trójwymiarowych spinów ewoluujących pod dynamiką Monte Carlo Metropolisa, kolorowanych barwą (kąt w płaszczyźnie) i jasnością (składowa poza płaszczyzną Sz), wraz z żywym wykresem magnetyzacji |M| w ostatnich przemiataniach, dzięki czemu można obserwować, jak układ się stabilizuje lub fluktuuje w czasie rzeczywistym.

🎮 Jak korzystać

Przeciągnij suwak temperatury, by przekroczyć próg uporządkowania, odwróć znak J, by przełączyć się między sprzężeniem ferromagnetycznym a antyferromagnetycznym, zmień rozmiar sieci dla drobniejszego lub grubszego szczegółu, i przełącz pole zewnętrzne, by zobaczyć spiny ciągnięte w stronę ustawienia wzdłuż z.

💡 Czy wiesz, że...?

W przeciwieństwie do trójwymiarowego modelu Heisenberga, ta dwuwymiarowa wersja nie ma prawdziwego porządku magnetycznego dalekiego zasięgu w żadnej niezerowej temperaturze (twierdzenie Mermina-Wagnera) — a mimo to wciąż wykazuje wyraźny porządek krótkiego zasięgu i powolną, szklistą dynamikę domen przy niskim T, fizykę ściśle związaną z przejściem Berezinskiego-Kosterlitza-Thoulessa w płaskich układach spinowych.

Najczęściej zadawane pytania

Dlaczego sieć wygląda jak kolorowy szum przy wysokiej temperaturze?

Przy wysokim T fluktuacje termiczne dominują nad sprzężeniem wymiennym J, więc algorytm Metropolisa akceptuje niemal każdą proponowaną reorientację. Każdy spin kończy wskazując zasadniczo losowy kierunek nieskorelowany z sąsiadami, co mapa kolorów barwy/jasności renderuje jako zaszumiony, nieskorelowany deseń.

Dlaczego przy niskiej temperaturze pojawiają się duże kolorowe domeny?

Przy niskim T energia wymiany dominuje, a źle ustawieni sąsiedzi kosztują zbyt dużo energii, by reguła Metropolisa łatwo je zaakceptowała. Spiny ustawiają się z sąsiadami w plamach, które rosną i łączą się w kolejnych przemiataniach, dając szerokie, płynnie zmieniające się kolorowe obszary widoczne na obrazie.

Co się dzieje, gdy ustawię J na wartość ujemną?

Ujemne J sprawia, że antyrównoległe sąsiednie spiny są energetycznie korzystniejsze niż równoległe. Przy niskiej temperaturze sieć rozwija wzór antyferromagnetyczny przypominający szachownicę, z naprzemiennymi orientacjami, zamiast jednolitych domen ferromagnetycznych.

Dlaczego wykres |M| wciąż fluktuuje nawet przy niskim T?

Ponieważ jest to skończona sieć z prawdziwym szumem termicznym, magnetyzacja nigdy się doskonale nie ustabilizuje — fluktuuje wokół wartości równowagowej, której wielkość zależy od temperatury i rozmiaru sieci. Większe sieci przy niskim T wykazują mniejsze względne fluktuacje, bliższe wyidealizowanej granicy termodynamicznej.

Co mówi mi statystyka szybkości akceptacji?

Szybkość akceptacji to frakcja proponowanych reorientacji spinów zaakceptowanych podczas ostatniego przemiatania. Jest zwykle wysoka przy wysokiej temperaturze (niemal każdy ruch jest akceptowany) i niższa przy niskiej temperaturze (tylko ruchy nieznacznie podnoszące energię lub ją obniżające mają szansę powodzenia), co czyni ją szybką diagnostyką tego, jak „gorący” lub „zamrożony” jest obecnie układ.

Dlaczego używać wektorów jednostkowych zamiast kątów?

Reprezentowanie każdego spinu jako trójwymiarowego wektora jednostkowego (Sx, Sy, Sz) czyni iloczyn skalarny Sᵢ·Sⱼ — a więc i energię wymiany — prostym, numerycznie stabilnym obliczeniem, i pozwala modelowi naturalnie uwzględnić przechylenie poza płaszczyznę (uchwycone tutaj jako jasność) obok obrotu w płaszczyźnie (uchwyconego jako barwa), czego pojedynczy kąt w płaszczyźnie nie mógłby reprezentować.

⚙ Pod maską

Symulacja Monte Carlo Metropolisa klasycznego modelu Heisenberga: trójwymiarowe spiny jednostkowe na dwuwymiarowej sieci porządkują się w widoczne domeny magnetyczne wraz ze spadkiem temperatury.

Heisenberg ModelMonte CarloMagnetismMetropolis AlgorithmCanvas 2D

2D · HTML5 Canvas 2D · cel 60 FPS · działa w całości po stronie klienta, bez instalacji

Co znalazłeś?

Dodaj kroki odtworzenia (opcjonalnie)