PL 
EnglishDomyślny
Українська🧲 Sieć spinowa 2D — ferromagnetyk Heisenberga
Symuluj dwuwymiarową sieć spinów 1/2 metodą Monte Carlo Metropolisa. Obserwuj spontaniczną magnetyzację pojawiającą się poniżej temperatury Curie.
O sieci spinowej 2D
Model Isinga (uproszczenie modelu Heisenberga dla cząstek spin-1/2) umieszcza magnetyczny spin s = +1 lub s = −1 w każdym węźle dwuwymiarowej sieci kwadratowej. Spiny najbliższych sąsiadów oddziałują poprzez sprzężenie wymienne J. Hamiltonian ma postać H = −J·ΣSᵢSⱼ − μH·ΣSᵢ, gdzie pierwsza suma przebiega po wszystkich parach najbliższych sąsiadów, a druga reprezentuje sprzężenie z polem zewnętrznym H.
Dla sprzężenia ferromagnetycznego (J > 0) równoległe spiny obniżają energię, tworząc tendencję do uporządkowania. Poniżej temperatury krytycznej Tc = 2J/(k·ln(1+√2)) ≈ 2,27 J/k fluktuacje termiczne są zbyt słabe, by zaburzyć ten porządek, i układ rozwija samorzutną magnetyzację — makroskopowy moment magnetyczny w nieobecności jakiegokolwiek pola zewnętrznego. Powyżej Tc energia termiczna przewyższa oddziaływanie wymienne i układ jest paramagnetyczny (nieuporządkowany). To jest ferromagnetyczne przejście fazowe.
Symulacja wykorzystuje algorytm Monte Carlo Metropolisa: w każdym kroku czasowym wykonywanych jest N² losowych prób odwrócenia spinu. Każde odwrócenie jest akceptowane, jeśli obniża energię, lub odrzucane/akceptowane z prawdopodobieństwem Boltzmanna, jeśli podnosi energię. Poprawnie próbkuje to rozkład równowagi termicznej w temperaturze T. Obserwuj, jak sieć samorzutnie się porządkuje, gdy obniżasz T poniżej Tc, i traci porządek, gdy podnosisz je powyżej.
Najczęściej zadawane pytania
Czym jest model Heisenberga ferromagnetyzmu?
Model Heisenberga opisuje oddziaływania magnetyczne w krysztale. Każdy atom niesie spin (kwantowy moment magnetyczny). Sąsiadujące spiny oddziałują poprzez sprzężenie wymienne J: jeśli J > 0 (ferromagnetyczne), równoległe spiny obniżają energię. Hamiltonian ma postać H = −J·ΣSᵢ·Sⱼ. W 2D jest to wyidealizowany model cienkiej warstwy magnetycznej.
Czym jest algorytm Metropolisa?
Algorytm Metropolisa to metoda Monte Carlo do próbkowania konfiguracji równowagi termicznej. W każdym kroku: (1) wybierz losowy spin; (2) zaproponuj jego odwrócenie; (3) oblicz ΔE; (4) zaakceptuj zawsze, jeśli ΔE < 0; (5) jeśli ΔE > 0, zaakceptuj z prawdopodobieństwem exp(−ΔE/kT). Poprawnie próbkuje to rozkład Boltzmanna.
Czym jest temperatura Curie?
Temperatura Curie Tc to temperatura krytyczna, powyżej której ferromagnetyk traci samorzutną magnetyzację. Poniżej Tc dominują oddziaływania wymienne i spiny się porządkują. Powyżej Tc fluktuacje termiczne losowo ustawiają spiny. Dla dwuwymiarowego modelu Isinga Tc ≈ 2,27 J/k. Żelazo ma Tc = 1043 K, nikiel 631 K, kobalt 1388 K.
Dlaczego materiały stają się magnetyczne?
Magnetyzm wynika z kwantowo-mechanicznego oddziaływania wymiennego między elektronami. Gdy atomy mają częściowo wypełnione orbitale d lub f (metale przejściowe, takie jak Fe, Co, Ni), nakładanie się orbitali prowadzi do energii wymiennej preferującej równoległe ustawienie spinów. W żelazie to oddziaływanie utrzymuje się do 1043 K. Obserwowany makroskopowy magnetyzm jest zbiorowym efektem ~10²³ wyrównanych spinów elektronowych.
Jaka jest różnica między ferromagnetyzmem a antyferromagnetyzmem?
W ferromagnetyzmie (J > 0) sąsiadujące spiny preferują ustawienie równoległe, co prowadzi do samorzutnej magnetyzacji objętościowej. W antyferromagnetyzmie (J < 0) sąsiadujące spiny preferują ustawienie antyrównoległe, tworząc wzór szachownicy. Wypadkowa magnetyzacja antyferromagnetyka wynosi zero. Przykłady: Fe, Co, Ni są ferromagnetyczne; MnO, FeO, NiO są antyferromagnetyczne.
Czym jest domena magnetyczna?
Domena magnetyczna to obszar, w którym wszystkie spiny są wyrównane. Domeny powstają, ponieważ zmniejszają energię magnetostatyczną, mimo że kosztują energię wymienną na ściankach domenowych. Niezmagnesowany kawałek żelaza ma wiele wzajemnie znoszących się domen. Pole zewnętrzne wyrównuje domeny; twarde magnesy zatrzymują domeny po usunięciu pola.
Czym jest spin-1/2 i jak wiąże się z rzeczywistymi atomami?
Spin-1/2 oznacza wewnętrzny moment pędu ℏ/2, z dwiema możliwymi projekcjami: +ℏ/2 (góra) lub −ℏ/2 (dół). Elektrony są cząstkami spin-1/2. W modelu Isinga każdy węzeł ma s = +1 lub s = −1 reprezentujące te stany. Rzeczywiste atomy magnetyczne mogą mieć większe spiny (np. Mn²⁺ ma spin-5/2), ale modele spin-1/2 oddają istotną fizykę zjawiska.
Czym jest ciepło właściwe i dlaczego osiąga szczyt w temperaturze Curie?
Ciepło właściwe C = dE/dT mierzy pochłanianie energii na jednostkę wzrostu temperatury. W pobliżu przejścia fazowego występują duże fluktuacje energii. Powoduje to szczyt (rozbieżność w układach nieskończonych) ciepła właściwego w Tc. W rzeczywistych ferromagnetykach ten szczyt w kształcie λ jest obserwowalny w eksperymentach kalorymetrycznych.
Czym jest dwuwymiarowy model Isinga i czy jest dokładnie rozwiązywalny?
Dwuwymiarowy model Isinga został dokładnie rozwiązany przez Larsa Onsagera w 1944 roku. Dokładne Tc = 2J/(k·ln(1+√2)) oraz magnetyzacja M ∝ (1 − sinh⁻⁴(2J/kT))^(1/8) są znane analitycznie. Trójwymiarowy model Isinga nie jest dokładnie rozwiązywalny i wymaga metod Monte Carlo lub grupy renormalizacji.
Jak metoda Monte Carlo jest wykorzystywana w fizyce poza modelami spinowymi?
Metody Monte Carlo są szeroko stosowane w chromodynamice kwantowej na sieci, zwijaniu białek, chemii kwantowej, transporcie promieniowania jądrowego, wycenie opcji w finansach oraz modelowaniu klimatu. Wspólnym wątkiem jest losowe próbkowanie do obliczania całek wysokowymiarowych lub próbkowania ze złożonych rozkładów prawdopodobieństwa.
Symuluj dwuwymiarową sieć spin-1/2 metodą Monte Carlo Metropolisa. Obserwuj, jak samorzutna magnetyzacja tworzy się poniżej temperatury Curie.
3D · Three.js / WebGL renderer · 60 FPS target · runs fully client-side, no install