Модель Изинга в двумерной решетке
Модель Изинга является краеугольным камнем статистической механики, предоставляя упрощенную основу для понимания фазовых переходов и кооперативных явлений. Данный симулятор реализует двумерную модель Изинга на квадратной решетке. Каждый узел решетки содержит магнитную спиновую переменную, σᵢ, которая может находиться в одном из двух состояний: 'вверх' (+1) или 'вниз' (–1). Энергия системы описывается гамильтонианом: H = –J Σ_{<i,j>} σᵢ σⱼ – h Σᵢ σᵢ. Первая сумма берется по всем парам ближайших соседей и описывает обменное взаимодействие с силой J > 0, благоприятствующее выстраиванию спинов (ферромагнетизм). Вторая сумма представляет взаимодействие с внешним магнитным полем h. Симулятор эволюционирует систему, используя алгоритм Монте-Карло Метрополиса — стохастический процесс, который выбирает спиновые конфигурации с вероятностью, пропорциональной фактору Больцмана, exp(–H/k_B T), где k_B — постоянная Больцмана, а T — температура. Отслеживая среднюю намагниченность на спин, m = (1/N) Σᵢ σᵢ, и среднюю энергию на спин при изменении температуры, можно наблюдать резкий непрерывный фазовый переход при критической температуре T_c ≈ 2.27 J/k_B. Ниже T_c происходит спонтанное нарушение симметрии, и система приобретает ненулевую намагниченность (|m| > 0) даже при h=0. Выше T_c тепловые флуктуации разрушают дальний порядок, и m усредняется к нулю. Работа с этим симулятором показывает, как микроскопические взаимодействия (выравнивание соседей) приводят к макроскопическому порядку, иллюстрирует концепцию универсальности и предоставляет вычислительное окно в поведение реальных магнитных материалов, бинарных сплавов и других систем, проявляющих коллективный порядок.
Для кого: Студенты старших курсов и аспиранты физических, химических и материаловедческих специальностей, изучающие статистическую механику, вычислительную физику или фазовые переходы.
Ключевые понятия
- Модель Изинга
- Фазовый переход
- Критическая температура (T_c)
- Алгоритм Метрополиса
- Спонтанное нарушение симметрии
- Ферромагнетизм
- Моделирование Монте-Карло
- Распределение Больцмана
Как это работает
Классическая демонстрация равновесной статмеханики: локальные перевороты с вероятностью принятия exp(−ΔE / kT) связывают микроскопическую динамику с феноменологией фазового перехода.
Часто задаваемые вопросы
- Почему намагниченность иногда меняет направление во время моделирования при низкой температуре?
- В системе конечного размера с периодическими граничными условиями и без внешнего поля (h=0) энергетический барьер между двумя вырожденными основными состояниями (все спины вверх против всех спинов вниз) конечен. За достаточно длительное время моделирования редкие, но большие тепловые флуктуации могут преодолеть этот барьер, вызывая переворот всей системы. Это эффект конечного размера; в бесконечно большой термодинамической системе вероятность такого переворота была бы бесконечно малой, и нарушение симметрии было бы истинно спонтанным.
- Какие реальные системы фактически описывает двумерная модель Изинга?
- Хотя изначально она создавалась для моделирования ферромагнетизма, бинарные степени свободы модели Изинга делают её универсальной решеточной моделью. Она может описывать переход порядок-беспорядок в бинарных сплавах (например, медь-цинк), критическую точку жидкость-газ через модель решеточного газа и даже упрощенные модели нейронных сетей или динамики общественного мнения. Её наибольшая ценность — в качестве парадигмальной модели для понимания универсальных особенностей непрерывных фазовых переходов.
- Симуляция кажется 'скачкообразной' или зашумленной, особенно вблизи критической температуры. Это ошибка?
- Нет, это фундаментальная физическая особенность. Вблизи критической точки (T_c) корреляционная длина — типичный размер кластеров выровненных спинов — расходится. Система демонстрирует критические флуктуации на всех масштабах, что приводит к большим, медленным колебаниям макроскопических величин, таких как намагниченность и энергия. Эта возросшая дисперсия и медленная динамика являются характерными признаками критической точки и корректно воспроизводятся алгоритмом Монте-Карло.
- В чем основное упрощение модели Изинга по сравнению с реальным магнитом?
- Основное упрощение заключается в том, что реальные атомные спины являются квантовомеханическими векторами, которые могут быть направлены в любую сторону, а не строго 'вверх' или 'вниз'. Более точной классической моделью является модель Гейзенберга. Модель Изинга ограничивает спины одной осью, что математически упрощает задачу, сохраняя при этом основную физику фазового перехода, обусловленного конкуренцией между энергией взаимодействия и тепловым беспорядком.
Ещё из «Термодинамика»
Другие симуляторы в этой категории — или все 18.
Симулятор идеального газа
Движущиеся частицы в ящике. Наблюдайте закон PV=nRT в действии.
Gas Laws Interactive
Boyle's, Charles's, Gay-Lussac's laws with interactive piston.
Теплопередача
Теплопроводность, конвекция и излучение при наличии градиентов температуры.
Фазовая диаграмма
Диаграмма температура-давление с указанием фазовых переходов.
Двигатель Карно
Анимация цикла на PV-диаграмме с этапами цикла и расчётом КПД.
Распределение Максвелла–Больцмана
Гистограмма |v| из гауссовых компонент в сравнении с 3D плотностью вероятности скорости Максвелла; T и размер выборки.