2D Спин Латівця: Вирощування Магніту з Монетки

Магніт – мільйони монетних подвоєнь, що узгоджуються один з одним

Феромагніт, такий як залізо, мікроскопічно є решіткою атомних магнітних моментусів — спінів. Кожен спін віддає перевагу спрямуватись у тому ж напрямку, що й його сусіди, оскільки це знижує локальну обмінну енергію. При високій температурі термічне тремтіння випадково змінює напрямок спінів швидше, ніж зв’язок може їх вирівняти, і матеріал не показує жодного магнітного моменту. Охолодіть його нижче критичної температури, і зв’язок перемагає: сусідні спіни починають узгоджуватися, ростуть домени з вирівняними спінами, і весь зразок може розвинути магнітний момент без зовнішнього поля. Це спонтанне магнетоутворення, і це один із найчистіших прикладів макроскопічного порядку, що виникає за рахунок чисто локальних взаємодій.

Енергія

Енергія решітки спінів з найближчим ферромагнітним зв’язком визначається як гамільтоніан Гейзенберга:

H = -J ∑⟨i,j⟩ Sᵢ · Sⱼ J = константа обміну (J > 0 сприяє паралельному, "ферромагнітному" вирівнюванню) ⟨i,j⟩ = сума по кожній парі найближчих сусідів на решітці, порахована один раз Sᵢ = спін на сайті i — вектор одиничного значення для спінів Гейзенберга, ±1 для спінів Ізинга Цей гамільтоніан навмисно простий: взаємодіють лише сусіди, а сила зв’язку однакова всюди. Ця простота є ключовою – це мінімальна модель, яка відтворює справжню ферромагнітну фазову перехід, і вона розв’язується точно в деяких випадках (2D модель Ізинга, наприклад, Ларсом Онсогером у 1944 році), що робить її базовою критерієм для перевірки кожної симуляції та наближення.

H = -J · Σ⟨i,j⟩ Sᵢ · Sⱼ

J     = exchange coupling constant (J > 0 favours parallel, "ferromagnetic" alignment)
⟨i,j⟩ = sum over each pair of NEAREST NEIGHBOURS on the lattice, counted once
Sᵢ    = spin at site i — a unit vector for Heisenberg spins, ±1 for Ising spins

Алгоритм Метапопування

Точне перерахування кожної конфігурації спінів, щоб знайти теплові середні, неможливе для будь-якої решітки більшого розміру, ніж кілька сайтів — простір станів зростає експоненціально. Алгоритм Метапопування (1953) замість цього будує випадковий рух через простір конфігурацій, частота відвідування якого в довгостроковій перспективі відповідає розподілу Больцмана, не обчислюючи функцію розділу:

повторювати для кожного кроку Монте-Карло: вибрати випадковий сайт i запропонувати перевернути (або змінити орієнтацію) спін Sᵢ обчислити ΔE = E_new - E_old (враховуються лише 4 найближчих сусіди — це дешево!) якщо ΔE <= 0: прийняти перевертання // завжди приймати покращення інакше: прийняти з ймовірністю exp(-ΔE / kT) // іноді приймати гірші стани інакше: відхилити, зберегти Sᵢ незмінним Ключовий момент полягає в тому, що перевертання локальне: оскільки змінюються лише 4 (або 6, у 3D) найближчі зв’язки між спінами, ΔE є кількома додаваннями незалежно від розміру решітки, отже один повний «прохід» по всіх сайтах дешевий. Запустіть достатньо проходів і ансамбль відвіданих конфігурацій задовольняє детальному балансу — правило прийняття/відхилення побудоване таким чином, що при рівновазі переходи між будь-якими двома станами відбуваються з такими швидкістями, як узгоджено з розподілом Больцмана при температурі T. Це обґрунтовує те, що середні значення за часом над конфігураціями, отриманими в результаті моделювання, є тепловими середніми.

repeat for each Monte Carlo step:
  pick a random site i
  propose flipping (or reorienting) spin Sᵢ
  compute ΔE = E_new - E_old   (only 4 neighbours matter — cheap!)
  if ΔE <= 0:
    accept the flip                         // always take improvements
  else:
    accept with probability exp(-ΔE / kT)   // sometimes take worse states
  else: reject, keep Sᵢ unchanged

Спостереження за переходом фази

Відстежуйте параметр порядку — середнє магнітне моменту на спін, m = (1/N)·Σ Sᵢ — під час повільного зниження температури. Вище критичної температури Tc, m коливається близько нуля: спіни практично не корельовані за дуже короткої відстані. Нижче Tc, m стабілізується біля ±1: майже кожен спін погоджується. Поруч із Tc самій, решітка показує домени будь-якого розміру одночасно, хаотично коливається, і довжина кореляції (наскільки поширюється вплив спіну) розбігається — відмітна риса безперервної (другої черги) фазової переходу, і причина того, що критичні явища біля Tc дотримуються універсального масштабування за законом степеня, яке поділяється на різних фізичних системах.

Практичні нотатки щодо моделювання

Ближче до Tc, окреме перемикання спінів методом Метапунктування страждає від сповільненого охолодження: великі області потребують дуже багато часу для перемикання одного спіну за раз, тому балансування та декореляція також стають повільними саме там, де ви найбільше хочете отримати хорошу статистику. Алгоритми кластерів, такі як Волфф і Швенсен-Ванг, вирішують цю проблему шляхом одномородного перемикання цілих корельованих кластерів спінів одним кроком, що значно скорочує ефективний час релаксації поблизу критичної точки. Для інтерактивної демонстрації на моніторів у реальному часі звичайне перемикання окремих спінів зазвичай достатньо — формування та розчинення доменів відбувається точно так, як це було призначено, навіть без прискорення на рівні кластерів.

Frequently asked questions

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

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