UK 
EnglishТипова
Polski🧲 2D Ґраткова Модель Спінів — Феромагнетик Гайзенберга
Симулюйте 2D ґратку спінів методом Монте-Карло Метрополіс. Спостерігайте спонтанну намагніченість нижче температури Кюрі.
Про 2D Ґраткову Модель Спінів
Модель Ізінга (спрощення моделі Гайзенберга для частинок спін-1/2) розміщує магнітний спін s = +1 або s = −1 у кожному вузлі 2D квадратної ґратки. Спіни найближчих сусідів взаємодіють через обмінну сталу J. Гамільтоніан — H = −J·ΣSᵢSⱼ − μH·ΣSᵢ, де перша сума береться по всіх парах найближчих сусідів, а друга представляє зв'язок із зовнішнім полем H.
Для феромагнітного зв'язку (J > 0) паралельні спіни знижують енергію, створюючи тенденцію до вирівнювання. Нижче критичної температури Tc = 2J/(k·ln(1+√2)) ≈ 2.27 J/k теплові флуктуації недостатні, щоб порушити цей порядок, і система розвиває спонтанну намагніченість — макроскопічний магнітний момент за відсутності будь-якого зовнішнього поля. Вище Tc теплова енергія долає обмінну взаємодію, і система стає парамагнітною (невпорядкованою). Це феромагнітний фазовий перехід.
Симуляція використовує алгоритм Монте-Карло Метрополіса: на кожному кроці часу здійснюється N² випадкових спроб перевернути спін. Кожен переворот приймається, якщо він знижує енергію, або приймається/відхиляється з ймовірністю Больцмана, якщо він підвищує енергію. Це коректно відбирає розподіл теплової рівноваги при температурі T. Спостерігайте, як ґратка спонтанно впорядковується при зниженні T нижче Tc, і розупорядковується при підвищенні вище неї.
Поширені запитання
Що таке модель Гайзенберга феромагнетизму?
Модель Гайзенберга описує магнітні взаємодії в кристалі. Кожен атом несе спін (квантовий магнітний момент). Сусідні спіни взаємодіють через обмінну сталу J: якщо J > 0 (феромагнетик), паралельні спіни знижують енергію. Гамільтоніан — H = −J·ΣSᵢ·Sⱼ. У 2D це ідеалізована модель тонкої магнітної плівки.
Що таке алгоритм Метрополіса?
Алгоритм Метрополіса — метод Монте-Карло для вибірки конфігурацій теплової рівноваги. На кожному кроці: (1) обирається випадковий спін; (2) пропонується перевернути його; (3) обчислюється ΔE; (4) переворот завжди приймається, якщо ΔE < 0; (5) якщо ΔE > 0, приймається з ймовірністю exp(−ΔE/kT). Це коректно відбирає розподіл Больцмана.
Що таке температура Кюрі?
Температура Кюрі Tc — критична температура, вище якої феромагнетик втрачає свою спонтанну намагніченість. Нижче Tc домінують обмінні взаємодії, і спіни вирівнюються. Вище Tc теплові флуктуації рандомізують спіни. Для 2D моделі Ізінга Tc ≈ 2.27 J/k. Залізо має Tc = 1043 K, нікель — 631 K, кобальт — 1388 K.
Чому матеріали стають магнітними?
Магнетизм виникає через квантово-механічну обмінну взаємодію між електронами. Коли атоми мають частково заповнені d- або f-орбіталі (перехідні метали, такі як Fe, Co, Ni), перекриття орбіталей приводить до обмінної енергії, що надає перевагу паралельному вирівнюванню спінів. У залізі цей обмін зберігається аж до 1043 K. Макроскопічний магнетизм, який ми спостерігаємо, — колективний результат ~10²³ вирівняних електронних спінів.
У чому різниця між феромагнетизмом та антиферомагнетизмом?
При феромагнетизмі (J > 0) сусідні спіни надають перевагу паралельності, що приводить до спонтанної об'ємної намагніченості. При антиферомагнетизмі (J < 0) сусідні спіни надають перевагу антипаралельності, утворюючи шаховий візерунок. Результуюча намагніченість антиферомагнетика дорівнює нулю. Приклади: Fe, Co, Ni феромагнітні; MnO, FeO, NiO антиферомагнітні.
Що таке магнітний домен?
Магнітний домен — область, у якій усі спіни вирівняні. Домени утворюються, оскільки вони знижують магнітостатичну енергію, хоча й коштують обмінної енергії на доменних стінках. Ненамагнічений шматок заліза має багато взаємоскасовуючих доменів. Зовнішнє поле вирівнює домени; тверді магніти утримують домени після зняття поля.
Що таке спін-1/2 і як він пов'язаний з реальними атомами?
Спін-1/2 означає власний момент імпульсу ħ/2 із двома можливими проєкціями: +ħ/2 (вгору) або −ħ/2 (вниз). Електрони — частинки спін-1/2. У моделі Ізінга кожен вузол має s = +1 або s = −1, що представляють ці стани. Реальні магнітні атоми можуть мати більші спіни (наприклад, Mn²⁺ має спін-5/2), але моделі спін-1/2 відображають основну фізику.
Що таке питома теплоємність і чому вона піковує при температурі Кюрі?
Питома теплоємність C = dE/dT вимірює поглинання енергії на приріст температури. Поблизу фазового переходу відбуваються великі флуктуації енергії. Це спричиняє пік (розбіжність у нескінченних системах) питомої теплоємності при Tc. У реальних феромагнетиках цей λ-подібний пік спостережуваний у калориметричних експериментах.
Що таке 2D модель Ізінга і чи вона точно розв'язна?
2D модель Ізінга була точно розв'язана Ларсом Онсагером у 1944 році. Точні Tc = 2J/(k·ln(1+√2)) і намагніченість M ∝ (1 − sinh⁻⁴(2J/kT))^(1/8) відомі аналітично. 3D модель Ізінга не є точно розв'язною і потребує методів Монте-Карло або ренормгрупи.