Strona główna Fizyka i Mechanika Defekty topologiczne — dynamika wirów

🌀 Defekty topologiczne — dynamika wirów

Poznaj defekty wirowe w dwuwymiarowym modelu XY. Obserwuj, jak dyslinacje +1/2 i −1/2 przyciągają się i anihilują w przejściu Bieriezinskiego-Kosterlitza-Thoulessa.

Fizyka i Mechanika2DZaawansowany60 FPS
topological-defects ↗ Otwórz osobno
Interfejs samej symulacji jest w języku angielskim.

Podobne symulacje

O tej symulacji

Ten symulator ewoluuje dwuwymiarowy model XY — klasyczne spiny o jednostkowej długości rozmieszczone na siatce, z których każdy może wskazywać dowolny kierunek w płaszczyźnie — zgodnie z dynamiką Monte Carlo Metropolisa w wybranej temperaturze. Wiry i antywiry to defekty topologiczne, w których kąt spinu nawija się o ±2π wokół punktu; poniżej temperatury przejścia Berezinskiego-Kosterlitza-Thoulessa (BKT) istnieją one wyłącznie jako ściśle związane pary, natomiast powyżej niej fluktuacje termiczne rozrywają pary i swobodne wiry zaczynają się mnożyć. Przeciąganie suwaka temperatury pozwala obserwować to rozpadanie się par na żywo, a przyciski Add Vortex / Add Antivortex umożliwiają ręczne umieszczenie defektów i obserwację, jak przeciwne ładunki spiralnie zbliżają się do siebie i anihilują.

🔬 Co przedstawia

Siatkę klasycznych spinów zabarwionych zgodnie z ich kątem w płaszczyźnie, na której wiry (nawinięcie +1/2) i antywiry (nawinięcie −1/2) pojawiają się jako małe wirujące wzory w polu kolorów. W niskiej temperaturze przetrwają tylko ściśle związane pary wir-antywir, połączone słabym przyciąganiem logarytmicznym; gdy temperatura wzrasta powyżej T_BKT ≈ 0,89 J/k_B, pary się rozpadają, a swobodne defekty rozpraszają się po siatce, niszcząc quasi-dalekozasięgowy porządek fazy niskotemperaturowej.

🎮 Jak korzystać

Użyj suwaka Temperature, by przesunąć układ poniżej lub powyżej przejścia BKT i obserwuj, jak pary wirów zachowują się inaczej po obu stronach. Naciśnij Play/Pause, by uruchomić ciągłą dynamikę Monte Carlo Metropolisa, lub Step, by przechodzić o jedno "przemiatanie" naraz i śledzić pojedyncze odwrócenia spinów. Kliknij Add Vortex (+1) i Add Antivortex (−1), by ręcznie umieścić parę defektów i obserwować ich przyciąganie, spiralne zbliżanie się i anihilację; Reset przywraca siatkę do stanu początkowego.

💡 Czy wiesz, że?

Nagroda Nobla z fizyki w 2016 roku została przyznana Kosterlitzowi, Thoulessowi i Haldane'owi między innymi za odkrycie, że przejście fazowe może zachodzić bez żadnego konwencjonalnego parametru porządku — przejście BKT pozostaje jednym z niewielu znanych przejść napędzanych wyłącznie rozpadem defektów topologicznych. Ta sama matematyka opisuje topnienie kryształów 2D, gdzie rolę wirów odgrywają dyslokacje.

Najczęściej zadawane pytania

Czym jest defekt topologiczny?

Defekt topologiczny to stabilna konfiguracja w uporządkowanym polu, której nie można w sposób ciągły rozwinąć do stanu jednorodnego. W dwuwymiarowym modelu XY wiry są punktowymi defektami, wokół których spin nawija się o ±2π. Liczba nawinięcia jest niezmiennikiem topologicznym — nie może się zmienić bez utworzenia lub zniszczenia pary defekt-antydefekt.

Czym jest dwuwymiarowy model XY?

Dwuwymiarowy model XY składa się z klasycznych spinów o stałej długości rozmieszczonych na siatce 2D, z których każdy może wskazywać dowolny kierunek w płaszczyźnie. Energia H = −J Σ cos(θᵢ−θⱼ) faworyzuje wyrównanie sąsiadów. W przeciwieństwie do dwuwymiarowego modelu Isinga, model XY ma symetrię ciągłą (obrotową), co prowadzi do jakościowo innej fizyki — w tym do przejścia BKT zamiast konwencjonalnego przejścia fazowego.

Czym jest przejście Berezinskiego-Kosterlitza-Thoulessa (BKT)?

Przejście BKT to topologiczne przejście fazowe w dwuwymiarowych układach z symetrią ciągłą. Poniżej T_BKT ≈ 0,89 J/k_B wiry istnieją wyłącznie jako związane pary (wir+antywir); układ ma quasi-dalekozasięgowy porządek z algebraicznie zanikającymi korelacjami. Powyżej T_BKT fluktuacje termiczne rozrywają pary, tworząc swobodne wiry, które niszczą quasi-porządek. Przejście BKT nie ma konwencjonalnego parametru porządku — było pierwszym przykładem topologicznego przejścia fazowego (Nagroda Nobla 2016).

Czym jest liczba nawinięcia?

Liczba nawinięcia zlicza, ile razy kierunek spinu obraca się o 2π podczas okrążania zamkniętej pętli. Wir +1 ma spiny nawinięte o +360° wokół swojego rdzenia; antywir −1 nawija się o −360°. Dla pętli nieotaczającej żadnego defektu liczba nawinięcia wynosi 0. Liczby nawinięcia sumują się algebraicznie: otoczenie defektu +1 i −1 daje nawinięcie 0.

Dlaczego pary wir-antywir się przyciągają?

Wiry o przeciwnych znakach przyciągają się poprzez potencjał logarytmiczny V ∝ −(J/π) ln(r/a), analogiczny do elektrostatyki 2D, gdzie ładunek wiru odgrywa rolę ładunku elektrycznego. To przyciąganie prowadzi do anihilacji pary wir-antywir: para zbliża się, łączy i znika, uwalniając energię w postaci fal spinowych. Poniżej T_BKT wszystkie pary są związane tym przyciąganiem.

Co dzieje się w temperaturze przejścia BKT?

W T_BKT koszt energii swobodnej związany z rozerwaniem pary wirów staje się porównywalny z zyskiem entropii wynikającym z obecności swobodnych wirów. Powyżej T_BKT wygrywa entropia i wiry mnożą się jako swobodne defekty. Długość korelacji rozbiega się w sposób osobliwości istotnej, gdy T→T_BKT od góry — jakościowo inaczej niż w zwykłych przejściach drugiego rodzaju.

Gdzie w rzeczywistych materiałach występuje przejście BKT?

Fizyka BKT pojawia się w cienkich warstwach helu-4 (przejście nadciekłe), dwuwymiarowych nadprzewodnikach, ultrazimnych gazach atomowych w 2D oraz w tablicach złączy Josephsona. Opisuje również topnienie kryształów 2D (teoria KTHNY), gdzie rolę wirów odgrywają dyslokacje. Nagrodę Nobla 2016 przyznano Kosterlitzowi, Thoulessowi i Haldane'owi za te i pokrewne odkrycia.

Czym różnią się disklinacje od dyslokacji?

Disklinacje to defekty rotacyjne, w których pole spinowe ma niezerową liczbę nawinięcia wokół rdzenia — to właśnie pokazuje ta symulacja. Dyslokacje to defekty translacyjne w sieciach krystalicznych, gdzie sieć ma wektor Burgersa (częściowe przesunięcie). Oba typy są defektami topologicznymi, ale w różnych parametrach porządku. Podczas topnienia kryształu 2D (teoria KTHNY) najpierw rozpadają się dyslokacje, a potem disklinacje, co daje dwa odrębne przejścia.

Dlaczego dwuwymiarowy model XY nie ma konwencjonalnego porządku dalekozasięgowego?

Twierdzenie Mermina-Wagnera dowodzi, że symetrie ciągłe nie mogą zostać spontanicznie złamane w 2D w skończonej temperaturze (dla oddziaływań krótkozasięgowych). Fale spinowe (mody Goldstone'a) kosztują jedynie logarytmicznie mało energii przy dużych długościach fali — ich entropia zawsze przeważa nad energetyczną tendencją do porządkowania. Poniżej T_BKT układ ma za to quasi-dalekozasięgowy porządek: korelacje ⟨cos(θᵢ−θⱼ)⟩ ∝ r^{-η(T)}, gdzie η zależy od temperatury.

Jak symuluje się model XY metodą Monte Carlo Metropolisa?

W algorytmie Metropolisa losowo wybierany jest jeden spin, a nowy kąt proponowany jest przez dodanie jednostajnie losowego przyrostu. Obliczana jest zmiana energii ΔE wynikająca ze zmiany tego jednego spinu (dotycząca tylko jego czterech sąsiadów). Jeśli ΔE < 0, ruch jest akceptowany. Jeśli ΔE > 0, jest akceptowany z prawdopodobieństwem exp(−ΔE/T), co realizuje zasadę szczegółowej równowagi i zapewnia zbieżność do rozkładu Boltzmanna w temperaturze T.