Головна Фізика та Механіка Топологічні Дефекти — Динаміка Вихорів

🌀 Топологічні Дефекти — Динаміка Вихорів

Досліджуйте вихрові дефекти в 2D XY-моделі. Спостерігайте за притяганням та анігіляцією дисклінацій у переході БКТ.

Фізика та Механіка2DСкладний60 FPS
topological-defects ↗ Відкрити окремо
DRAG · SCROLL · CLICK — керуйте прямо у вікні симуляції.

Про цю симуляцію

Ця симуляція моделює 2D XY-модель — класичні спіни фіксованої величини, розміщені на двовимірній ґратці, кожен з яких може вказувати в довільному напрямку в площині; енергія H = −J Σ cos(θᵢ−θⱼ) сприяє вирівнюванню сусідніх спінів. Топологічний дефект — це стабільна конфігурація впорядкованого поля, яку неможливо неперервно розкрутити до однорідного стану: у 2D XY-моделі вихори є точковими дефектами, навколо ядра яких спіни закручуються на ±2π. Число закручування є топологічним інваріантом — воно не може змінитися без утворення чи знищення пари дефект-антидефект.

🔬 Що показано

Симуляція показує вихори +1 та антивихори −1 у 2D XY-моделі, що притягуються один до одного через логарифмічний потенціал V ∝ −(J/π) ln(r/a), аналогічний до 2D-електростатики, де роль електричного заряду відіграє «вихровий заряд». Це притягання призводить до анігіляції пари вихор-антивихор: пара зближується, зливається і зникає, вивільняючи енергію у вигляді спінових хвиль. Нижче T_БКТ усі пари зв'язані цим притяганням.

🎮 Як користуватися

У симуляції можна спостерігати, як вихрові дефекти виникають і рухаються при зміні температури: нижче критичної температури переходу БКТ вихори існують лише у зв'язаних парах, а вище неї — вільно розповсюджуються по ґратці. Використовуйте кнопки «Скинути» та «На весь екран» на панелі інструментів, щоб перезапустити симуляцію або розглянути динаміку дисклінацій детальніше.

💡 Чи знали ви?

БКТ-перехід був першим прикладом топологічного фазового переходу — він не має звичайного параметра порядку. За це та пов'язані відкриття Костерліц, Таулесс і Голдейн отримали Нобелівську премію з фізики 2016 року. Фізика БКТ проявляється в тонких плівках гелію-4 (надплинний перехід), двовимірних надпровідниках, ультрахолодних атомних газах та масивах контактів Джозефсона.

Часті запитання

Що таке топологічний дефект?

Топологічний дефект — це стабільна конфігурація впорядкованого поля, яку неможливо неперервно розкрутити до однорідного стану. У 2D XY-моделі вихори є точковими дефектами, навколо ядра яких спіни закручуються на ±2π. Число закручування є топологічним інваріантом — воно не може змінитися без утворення чи знищення пари дефект-антидефект.

Що таке 2D XY-модель?

2D XY-модель складається з класичних спінів фіксованої величини, розміщених на двовимірній ґратці, кожен з яких може вільно вказувати в будь-якому напрямку в площині. Енергія H = −J Σ cos(θᵢ−θⱼ) сприяє вирівнюванню сусідніх спінів. На відміну від 2D моделі Ізінга, XY-модель має неперервну симетрію (обертання), що зумовлює якісно іншу фізику — зокрема перехід БКТ замість звичайного фазового переходу.

Що таке перехід Березинського–Костерліца–Таулесса (БКТ)?

Перехід БКТ — це топологічний фазовий перехід у двовимірних системах із неперервною симетрією. Нижче T_БКТ ≈ 0,89 J/k_B вихори існують лише у зв'язаних парах (вихор+антивихор); система має квазідалекий порядок з алгебраїчно спадаючими кореляціями. Вище T_БКТ теплові флуктуації розривають пари, створюючи вільні вихори, що руйнують квазіпорядок. Перехід БКТ не має звичайного параметра порядку — він став першим прикладом топологічного фазового переходу (Нобелівська премія 2016 року).

Що таке число закручування?

Число закручування підраховує, скільки разів напрямок спіну повертається на 2π при обході замкненого контуру. Вихор +1 має спіни, що закручуються на +360° навколо свого ядра; антивихор −1 закручується на −360°. Для контуру, що не охоплює жодного дефекту, число закручування дорівнює 0. Числа закручування додаються алгебраїчно: контур, що охоплює дефекти +1 і −1, дає закручування 0.

Чому пари вихор-антивихор притягуються?

Вихори протилежного знаку притягуються через логарифмічний потенціал V ∝ −(J/π) ln(r/a), аналогічний до 2D-електростатики, де роль електричного заряду відіграє вихровий заряд. Це притягання призводить до анігіляції пари вихор-антивихор: пара зближується, зливається і зникає, вивільняючи енергію у вигляді спінових хвиль. Нижче T_БКТ усі пари зв'язані цим притяганням.

Що відбувається при температурі переходу БКТ?

При T_БКТ енергетична вартість розриву пари вихорів стає порівнянною з виграшем ентропії від наявності вільних вихорів. Вище T_БКТ ентропія перемагає, і вихори розмножуються як вільні дефекти. Кореляційна довжина розходиться за характером істотної особливості при T→T_БКТ згори — якісно відмінно від звичайних фазових переходів другого роду.

Де перехід БКТ трапляється в реальних матеріалах?

Фізика БКТ проявляється в тонких плівках гелію-4 (надплинний перехід), двовимірних надпровідниках, ультрахолодних атомних газах у 2D і масивах контактів Джозефсона. Вона також описує плавлення двовимірних кристалів (теорія КТХНІ), де роль вихорів відіграють дислокації. Нобелівську премію з фізики 2016 року присуджено Костерліцу, Таулессу й Голдейну за ці та пов'язані відкриття.

Чим дисклінації відрізняються від дислокацій?

Дисклінації — це обертальні дефекти, при яких спінове поле має ненульове число закручування навколо ядра — саме це показує ця симуляція. Дислокації — це трансляційні дефекти в кристалічних ґратках, де ґратка має вектор Бюргерса (частковий зсув). Обидва типи є топологічними дефектами, але в різних параметрах порядку. При плавленні 2D-кристалів (теорія КТХНІ) спочатку розриваються дислокації, а потім дисклінації, що дає два окремі переходи.

Чому 2D XY-модель не має звичайного далекого порядку?

Теорема Мерміна–Вагнера доводить, що неперервні симетрії не можуть бути спонтанно порушені у 2D за скінченної температури (для короткодіючих взаємодій). Спінові хвилі (голдстоунівські моди) коштують лише логарифмічну енергію на великих довжинах хвиль — їхня ентропія завжди переважає енергетичну тенденцію до впорядкування. Нижче T_БКТ система натомість має квазідалекий порядок: кореляції ⟨cos(θᵢ−θⱼ)⟩ ∝ r^{-η(T)}, де η залежить від температури.

Як XY-модель моделюється методом Монте-Карло Метрополіса?

В алгоритмі Метрополіса один спін обирається випадково, і пропонується новий кут шляхом додавання рівномірного випадкового приросту. Обчислюється зміна енергії ΔE від зміни цього одного спіна (що стосується лише його чотирьох сусідів). Якщо ΔE < 0, рух приймається. Якщо ΔE > 0, він приймається з імовірністю exp(−ΔE/T), реалізуючи детальний баланс і забезпечуючи збіжність до розподілу Больцмана при температурі T.

Схожі симуляції