Изинговая модель: фазовые переходы и метод Монте-Карло

Спін на сітці

Модель Ейзінга, запропонована Вільгельмом Ленцем у 1920 році та вирішена в одімені його учнем Ернстом Ейзінгом у 1925 році, розміщує дискретну змінну — спін, s = +1 або −1 — на кожній точці решітки. Кожен спін взаємодіє лише з найближчими сусідами, а загальна енергія конфігурації становить:

E = -J * ∑(s_i * s_j) для всіх сусідніх пар — J > 0 → сусіди віддають перевагу вирівнюванню (феромагнітний) J < 0 → сусіди віддають перевагу антивирівнюванню (антиферомагнітний) B → зовнішнє магнітне поле, упереджує спіни до +1 або -1 Це вся модель: немає квантової механіки, жодних безперервних полів, лише дискретні спіни та їх найближчі сусіди. Що робить її вражаючою полягає в тому, що цей надзвичайно спрощений набір правил достатньо для відтворення справжньої, чіткої фазової перемикання.

E = -J * sum_over_neighbour_pairs(s_i * s_j)  -  B * sum_over_all_sites(s_i)

J > 0  → neighbours prefer to align (ferromagnetic)
J < 0  → neighbours prefer to anti-align (antiferromagnetic)
B      → external magnetic field, biases spins toward +1 or -1

Метод Метаполіса

Точне перерахування кожної конфігурації спінів неможливе для будь-якого реального розміру решітки — наприклад, решітка 20x20 має 2^400 конфігурацій. Замість цього модель зращується за допомогою алгоритму Метаполіса-Гастінгса (1953): повторно пропонувати перевернути випадково обраний спін, обчислювати отриману зміну енергії ΔE та ймовірнісно приймати або відхиляти перетворення. Виберіть випадкове положення i ΔE = енергія_зміни_якщо_ми_перевернемо(i) // залежить лише від 4 сусідів i' Якщо ΔE <= 0: Перевернути(i) // завжди приймайте рух, який знижує енергію Інакше: Прийняти з ймовірністю exp(-ΔE / (k_B * T)) // іноді приймайте підйом енергії Повторити мільйони разів → решітка встановлюється в розподіл Больцмана за температурою T Це правило «один спін за раз» є важливим: воно задовольняє детальному балансу, математичну умову, яка гарантує збіжність моделювання до правильного теплового рівноважного розподілу замість іншого, нефізичного стаціонарного стану. Перевертання багатьох спінів з однієї знімка одночасно порушить цю гарантію.

pick a random site i
dE = energy_change_if_we_flip(i)     // only depends on i's 4 neighbours
if dE <= 0:
    flip(i)                          // always accept a move that lowers energy
else:
    accept with probability exp(-dE / (k_B * T))   // sometimes accept a rise

repeat millions of times → the lattice settles into the Boltzmann
distribution at temperature T

Критична температура та спонтанне магнітне упорідження

При високій температурі домінують термальні шуми: спіни випадково перевертаються, незалежно від сусідів, і середнє магнітне упорідження дорівнює нулю. При низькій температурі домінує взаємний вплив J: сусідні спіни вирівнюються у великі вирівняні області, і решітка розвиває ненульове спонтанне магнітне упорідження навіть без зовнішнього поля. Ларс Оґнсен розв’язав точну 2D модель Ising з квадратним лабіринтом у 1944 році та виявив різкий перехід за точної критичної температури, k_B T_c / J = 2 / ln(1 + √2) ≈ 2.269. Прямо при T_c довжина кореляції між спінами розбігається — ділянка вирівняних спінів може охопити весь ланцюг — що є відміткою справжньої безперервної (другої черги) фазової перемикання.

Чому фізики все ще цікавляться цією імітаційною моделлю

Реальна цінність моделі Ізінга полягає не в точному описі справжніх магнітів – вона є однією з небагатьох взаємодіючих статистичних моделей, які можна розв’язати аналітично, що робить її базовою відправною точкою для валідації будь-якого методу наближення (меандерна теорія, ренормована група, методи Монте-Карло). Її універсальна класова – конкретний спосіб масштабування фізичних величин поблизу критичної точки – виявляється описувати величезна кількість несхожих фізичних систем поблизу їх власних критичних точок, від двофазного розшарування бінарних рідин до деяких моделей нейронної активності, оскільки в умовах безперержної фазової перехід мікроскопічні деталі системи перестають мати значення, і лише її розмірність та симетрія визначають критичну поведінку.

Frequently asked questions

Спробуйте наживо

Усе, що вище, працює прямо у вашому браузері — відкрийте Ising Model і змінюйте параметри під час роботи. Нічого не встановлюється, нічого не завантажується на сервер, уся модель живе в одній вкладці.