Як це працює
Кожен вузол ґратки містить класичний спін: одиничний вектор (Sx, Sy, Sz), що може вказувати будь-куди на сфері, а не лише "вгору" чи "вниз", як у простішій моделі Ізінга. Сусідні спіни взаємодіють через обмінну енергію E = -J·Σ Sᵢ·Sⱼ, підсумовану по парах найближчих сусідів з періодичними граничними умовами. Додатне J сприяє паралельному вирівнюванню (феромагнетизм); від'ємне J сприяє антипаралельному вирівнюванню (антиферомагнетизм).
Симуляція використовує алгоритм Монте-Карло Метрополіса: кожен кадр анімації виконує повний прохід (в середньому одна спроба оновлення на вузол). Для кожної спроби випадково обирається вузол, пропонується нова випадкова орієнтація одиничного вектора, і обчислюється зміна енергії ΔE від чотирьох сусідів по ґратці. Ходи, що знижують енергію, завжди приймаються; ходи, що підвищують її, приймаються з ймовірністю exp(-ΔE / kBT). За багато проходів це веде ґратку до термодинамічної рівноваги при обраній температурі T.
Прийняття Метрополіса: P = min(1, exp(-ΔE / kBT))
Намагніченість: |M| = |(1/N)·Σᵢ Sᵢ|, 0 ≤ |M| ≤ 1
Кольорове відображення: відтінок = atan2(Sy, Sx), яскравість ∝ Sz
Часті запитання
Що таке модель спінів Гейзенберга?
Класична модель Гейзенберга представляє кожен магнітний момент як неперервний 3D одиничний вектор, а не фіксований напрямок. Сусідні спіни взаємодіють через обмінну енергію E = -J·Sᵢ·Sⱼ.
Чим модель Гейзенберга відрізняється від моделі Ізінга?
Модель Ізінга обмежує кожен спін двома дискретними станами. Модель Гейзенберга дозволяє спінам вказувати в будь-якому напрямку на одиничній сфері — більш реалістичний опис багатьох реальних магнітних матеріалів.
Що робить алгоритм Монте-Карло Метрополіса?
На кожному кроці обирається випадковий вузол і пропонується нова випадкова орієнтація спіну. Якщо це знижує енергію, вона завжди приймається; якщо підвищує — приймається з ймовірністю exp(-ΔE / kBT).
Чому ймовірність прийняття дорівнює exp(-ΔE/kBT)?
Це фактор Больцмана зі статистичної механіки. Використання його як правила прийняття гарантує детальну рівновагу, тому змодельовані конфігурації збігаються до правильного термодинамічного ансамблю.
Що вимірює намагніченість |M|?
Намагніченість — векторне середнє всіх спінів, M = (1/N)·Σ Sᵢ. Її величина коливається від 0 (випадкові орієнтації) до 1 (всі спіни вирівняні) — найпростіший параметр порядку.
Чому 2D модель Гейзенберга не має справжнього фазового переходу при скінченній температурі?
Теорема Мерміна-Вагнера показує, що неперервні симетрії не можуть бути спонтанно порушені в 1D чи 2D короткодіючих системах при будь-якій T > 0. Ґратка все ще показує сильні кореляції ближнього порядку при низькій T.
Що представляє константа зв'язку J?
J — сила обмінної взаємодії між сусідніми моментами. Додатне J (феромагнітне) сприяє паралельним спінам; від'ємне J (антиферомагнітне) сприяє шаховому антипаралельному візерунку при низькій T.
Що робить керування зовнішнім полем?
Увімкнення поля додає доданок Зеємана -h·Sᵢz до енергії кожного спіну, сприяючи вирівнюванню вздовж z. Це імітує прикладене магнітне поле Bz.
Де класична модель Гейзенберга використовується в реальних дослідженнях?
Це стандартний полігон для феромагнетиків, антиферомагнетиків і фрустрованих магнетиків, що лежить в основі досліджень спінтроніки та квантового магнетизму.
Що таке прохід Монте-Карло (sweep)?
Один прохід намагається виконати одне оновлення Метрополіса для кожного вузла ґратки в середньому. Багато проходів дозволяють конфігурації релаксувати від випадкового старту до рівноваги.