Spotlight #66: Fizyka materii skondensowanej — Ising, BEC, nadprzewodniki i szkła spinowe

Oto zaskakujący fakt: mniej więcej 40% wszystkich publikowanych badań z fizyki dotyczy fizyki materii skondensowanej. W tej dziedzinie przyznano więcej Nagród Nobla niż w jakiejkolwiek innej gałęzi fizyki. Mimo to sama nazwa jest niemal nieznana poza środowiskiem akademickim, ponieważ to, co badają fizycy materii skondensowanej, zwykle występuje pod innymi nazwami, gdy dociera do szerszej publiczności: nadprzewodnictwo, półprzewodniki, magnetyzm, nadciekłość, izolatory topologiczne. Symulacje w tym przeglądzie pozwalają zbadać pięć fundamentalnych zjawisk leżących u podstaw tej ogromnej i niezwykle ważnej dziedziny.

I. Model Isinga — jak porządek wyłania się z chaosu

Model Isinga jest zwodniczo prosty: umieść na siatce zbiór spinów, z których każdy może być skierowany w górę (+1) lub w dół (−1). Każdy spin oddziałuje jedynie ze swoimi najbliższymi sąsiadami. Mimo to ten model-zabawka oddaje istotę fizyki ferromagnetycznych przejść fazowych — spontaniczne pojawienie się makroskopowego porządku z mikroskopowej losowości, gdy temperatura spada poniżej wartości krytycznej.

Energia układu jest opisana hamiltonianem:

H = -J * sum_{} s_i * s_j - mu * B * sum_i s_i

  J    : stała sprzężenia (J > 0 ferromagnetyk, J < 0 antyferromagnetyk)
  s_i  : spin w węźle i (+1 lub -1)
   : suma po parach najbliższych sąsiadów
  mu   : moment magnetyczny
  B    : zewnętrzne pole magnetyczne

Temperatura krytyczna (2D sieć kwadratowa, dokładne rozwiązanie Onsagera):
  k_B * T_c = 2J / ln(1 + sqrt(2))  ≈ 2.269 J

Parametr porządku: namagnesowanie m = (1/N) * sum_i s_i
  m ≈ 0   dla T >> T_c  (nieuporządkowana faza paramagnetyczna)
  m ≈ ±1  dla T << T_c  (uporządkowana faza ferromagnetyczna)

Przejście fazowe przy T_c jest przejściem drugiego rodzaju: namagnesowanie rośnie w sposób ciągły od zera, gdy temperatura spada poniżej punktu krytycznego. W pobliżu T_c układ wykazuje zjawiska krytyczne — fluktuacje niezmiennicze skalowo, korelacje o charakterze potęgowym oraz rozbieżną podatność — opisywane przez uniwersalne wykładniki krytyczne, niezależne od szczegółów mikroskopowych. Dwuwymiarowy model Isinga jest jednym z bardzo nielicznych dokładnie rozwiązywalnych modeli w mechanice statystycznej, rozwiązanym przez Larsa Onsagera w 1944 roku w prawdziwym majstersztyku fizyki matematycznej.

Symulacja wykorzystuje algorytm Metropolisa-Hastingsa do próbkowania rozkładu Boltzmanna: proponowane jest losowe odwrócenie spinu, które zostaje przyjęte, jeśli obniża energię, a przyjęte z prawdopodobieństwem exp(-ΔE / k_B T), jeśli energię podwyższa. To podejście Monte Carlo naturalnie generuje fluktuacje termiczne, ściany domenowe i stopniowe zarodkowanie porządku magnetycznego w miarę przesuwania suwaka temperatury poniżej T_c.

II. Kondensacja Bosego-Einsteina — piąty stan materii

W 1924 roku Satyendra Nath Bose i Albert Einstein przewidzieli, że bozony — cząstki o spinie całkowitym — przy dostatecznie niskich temperaturach zapadną się wszystkie do jednego, najniższego energetycznie dostępnego im stanu kwantowego. W przeciwieństwie do fermionów (które podlegają zakazowi Pauliego i muszą zajmować różne stany), bozony chętnie tłoczą się w tym samym stanie. Rezultatem jest makroskopowy obiekt kwantowy: kondensat Bosego-Einsteina (BEC).

Przejście zachodzi wtedy, gdy termiczna długość fali de Broglie'a cząstek staje się porównywalna z odległością między cząstkami — mniej więcej wtedy, gdy funkcje falowe kwantowe zaczynają się nakładać:

Termiczna długość fali de Broglie'a:
  lambda_dB = h / sqrt(2 * pi * m * k_B * T)

Temperatura krytyczna (3D idealny gaz Bosego):
  T_c = (h^2 / (2 * pi * m * k_B)) * (n / zeta(3/2))^(2/3)
  zeta(3/2) ≈ 2.612  (funkcja dzeta Riemanna)
  n : gęstość liczbowa

Frakcja kondensatu poniżej T_c:
  N_0 / N = 1 - (T / T_c)^3

BEC jest opisywany przez pojedynczą makroskopową funkcję falową (parametr porządku):
  Psi(r, t) = sqrt(n_0) * exp(i * phi(r, t))
  n_0 : gęstość kondensatu
  phi : faza makroskopowa (spójna w całym kondensacie)

Kondensaty Bosego-Einsteina zostały po raz pierwszy stworzone w laboratorium w 1995 roku przez Erica Cornella i Carla Wiemana (pracujących z atomami rubidu) oraz przez Wolfganga Ketterlego (pracującego z atomami sodu), za co otrzymali Nagrodę Nobla w dziedzinie fizyki w 2001 roku. Wymagane temperatury — około 100 nanokelwinów, najzimniejsze miejsca w znanym wszechświecie — wymagają zaawansowanych technik chłodzenia laserowego i pułapkowania magnetycznego.

Najbardziej uderzającą właściwością BEC jest jego nadciekłość: poniżej temperatury kondensacji spójność kwantowa kondensatu tłumi dyssypację, a płyn płynie z zerową lepkością. Zamieszaj nadciekły płyn wystarczająco szybko, a wytworzy skwantowane wiry — linie wirowe, wokół których faza parametru porządku nawija się dokładnie o 2π, każda niosąc jeden kwant momentu pędu.

III. Złącza Josephsona — kwantowe tunelowanie w skali makroskopowej

Złącze Josephsona to jedno z najbardziej niezwykłych urządzeń kwantowych, jakie kiedykolwiek wymyślono: dwa nadprzewodniki oddzielone cienką barierą izolacyjną, przez którą tunelują kwantowo-mechanicznie pary Coopera (związane pary elektronów odpowiedzialne za nadprzewodnictwo). Całe złącze opisują zaledwie dwie liczby — różnica faz φ między makroskopowymi funkcjami falowymi dwóch nadprzewodników oraz prąd krytyczny I_c. Jednak z tych dwóch liczb wyłaniają się dwa fundamentalne zjawiska.

Efekt Josephsona DC (zerowe napięcie, niezerowy prąd):
  I = I_c * sin(phi)
  phi : różnica faz na złączu

Efekt Josephsona AC (niezerowe napięcie wytwarza oscylujący prąd):
  d(phi)/dt = (2e / hbar) * V
  f_J = (2e / h) * V  (częstość Josephsona)
  f_J ≈ 483.6 GHz / mV  (dokładna wartość, używana jako wzorzec napięcia)

Model RCSJ (złącze bocznikowane rezystancyjnie i pojemnościowo):
  C * d^2(phi)/dt^2 + (1/R) * d(phi)/dt + I_c * sin(phi) = I_bias

  Jest to matematycznie identyczne z tłumionym, wymuszonym wahadłem:
  m * d^2(theta)/dt^2 + b * d(theta)/dt + (mg/L) * sin(theta) = tau

Tłumienie prądu krytycznego w polu magnetycznym (wzór Fraunhofera):
  I_c(Phi) = I_c0 * |sin(pi * Phi / Phi_0) / (pi * Phi / Phi_0)|
  Phi_0 = h / 2e ≈ 2.07e-15 Wb  (kwant strumienia magnetycznego)

Efekt Josephsona DC jest niezwykły: nadprąd płynie przy zerowym przyłożonym napięciu, napędzany wyłącznie przez kwantową różnicę faz. Efekt Josephsona AC — w którym napięcie stałe wytwarza oscylujący nadprąd o precyzyjnie określonej częstotliwości — jest tak dokładny, że dziś definiuje wolt w jednostkach SI. Współczesne wzorce napięcia na całym świecie opierają się na macierzach złączy Josephsona.

IV. Szkła spinowe — frustracja i zamrożony nieporządek

Co się stanie, gdy weźmiemy model Isinga i sprawimy, że stałe sprzężenia J_{ij} będą losowe — czasem ferromagnetyczne, czasem antyferromagnetyczne, dla każdego wiązania wybierane losowo? Rezultatem jest szkło spinowe: układ tak sfrustrowany, że nie potrafi znaleźć prostego uporządkowanego stanu podstawowego, i zamiast tego zamarza w jednej z astronomicznie wielu metastabilnych konfiguracji, wszystkie o niemal takiej samej energii.

Hamiltonian szkła spinowego:
  H = -sum_{} J_ij * s_i * s_j

Parametr porządku Edwardsa-Andersona (mierzy lokalne zamrożenie, nie porządek globalny):
  q_EA = lim_{t -> inf} (1/N) * sum_i 
  q_EA = 0  w fazie paramagnetycznej (spiny fluktuują swobodnie)
  q_EA > 0  w fazie szkła spinowego (spiny zamrożone, ale losowo zorientowane)

Frustracja: trójkąt wiązań antyferromagnetycznych (J < 0) jest sfrustrowany,
  ponieważ niemożliwe jest spełnienie wszystkich trzech wiązań jednocześnie.
  Na sieci trójkątnej z J < 0 wszędzie:
    każda konfiguracja pozostawia co najmniej jedno wiązanie niespełnione.

Krajobraz energii swobodnej:
  Eksponencjalnie wiele minimów lokalnych, oddzielonych barierami proporcjonalnymi do sqrt(N).
  Czas równoważenia ~ exp(alpha * N^sigma): rośnie szybciej niż jakikolwiek wielomian.

Szkła spinowe wykazują starzenie się — im starsza próbka, tym wolniej reaguje na zaburzenie — oraz efekty pamięci, dzięki którym szkło spinowe „pamięta” swoją historię termiczną. Te właściwości czynią szkła spinowe prototypem znacznie szerszej klasy problemów: każdego złożonego układu z wieloma konkurującymi ograniczeniami i eksponencjalnie postrzępionym krajobrazem energii. Zwijanie białek, sieci neuronowe, heurystyki optymalizacyjne, a nawet dynamikę społeczną analizowano za pomocą teorii szkła spinowego.

Model Sherringtona-Kirkpatricka (SK) rozszerza to na oddziaływania o nieograniczonym zasięgu i został rozwiązany przez Giorgia Parisiego w latach 80. za pomocą metody replik — jednego z najbardziej wymagających technicznie obliczeń w fizyce teoretycznej, za które Parisi otrzymał Nagrodę Nobla w dziedzinie fizyki w 2021 roku.

V. Topologiczne stany brzegowe — model SSH

Odkrycie, że stany kwantowe materii można klasyfikować nie tylko na podstawie łamania symetrii (jak teoria Landaua klasyfikuje ferromagnetyki i nadprzewodniki), ale też na podstawie topologii — globalnej właściwości funkcji falowej, której nie można zmienić poprzez gładkie deformacje — było jednym z przełomowych odkryć fizyki końca XX wieku, za które Thouless, Haldane i Kosterlitz otrzymali Nagrodę Nobla w 2016 roku. Najprostszym modelem oddającym fizykę topologiczną jest łańcuch Su-Schriffera-Heegera (SSH).

Hamiltonian SSH (1D łańcuch dimerów):
  H = -sum_i [ t1 * c_{i,A}^dag * c_{i,B} + t2 * c_{i,B}^dag * c_{i+1,A} ] + h.c.

  t1 : amplituda przeskoku wewnątrz komórki
  t2 : amplituda przeskoku między komórkami

Hamiltonian Blocha w przestrzeni pędów:
  H(k) = (t1 + t2 * cos(k)) * sigma_x + t2 * sin(k) * sigma_y

Liczba nawinięć (niezmiennik topologiczny):
  nu = (1 / 2pi) * integral_BZ d(phi(k))
  phi(k) = arg(t1 + t2 * exp(ik))

  nu = 0  jeśli |t1| > |t2|  (faza trywialna, brak stanów brzegowych)
  nu = 1  jeśli |t1| < |t2|  (faza topologiczna, stany brzegowe o zerowej energii)

Odpowiedniość brzeg-objętość:
  Łańcuch w fazie topologicznej (nu = 1) z otwartymi warunkami brzegowymi
  ma dokładnie dwa stany o zerowej energii, zlokalizowane na obu końcach,
  niezależnie od długości łańcucha (dla dostatecznie dużego N).

Topologiczne przejście fazowe zachodzi przy t1 = t2: przerwa energetyczna w paśmie objętościowym zamyka się i ponownie otwiera, a liczba nawinięć przeskakuje z 0 na 1 (lub odwrotnie). Stany brzegowe są „chronione topologicznie” — nie można ich usunąć żadnym zaburzeniem, które zachowuje symetrie hamiltonianu i nie zamyka przerwy objętościowej. W rzeczywistych materiałach ta odporność oznacza, że przewodnictwo brzegowe w izolatorach topologicznych jest odporne na domieszki, nieporządek i niedoskonałości powierzchni.

Model SSH pierwotnie zaproponowano do opisu przewodzących polimerów (poliacetylenu) w 1979 roku, i pomógł wyjaśnić, dlaczego domieszkowany poliacetylen mógł przewodzić prąd elektryczny mimo posiadania nominalnej przerwy energetycznej. Dziś fizyka SSH pojawia się w zimnych atomach w sieciach optycznych, kryształach fotonicznych, akustycznych metamateriałach oraz stanach brzegowych dwuwymiarowych izolatorów topologicznych, takich jak studnie kwantowe HgTe — wszystkie opisywane przez ten sam niezmiennik topologiczny i tę samą odpowiedniość brzeg-objętość.

Wypróbuj sam

Pięć symulacji omówionych w tym przeglądzie obejmuje centralne koncepcje fizyki materii skondensowanej. Sugerujemy eksplorowanie ich w przedstawionej kolejności: model Isinga buduje intuicję dotyczącą przejść fazowych; symulacja BEC pokazuje, co się dzieje, gdy to statystyka kwantowa, a nie oddziaływania, napędza przejście fazowe; złącze Josephsona demonstruje spójność kwantową w skali makroskopowej; szkło spinowe ujawnia, co frustracja robi ze zdolnością układu do porządkowania się; a model SSH wprowadza perspektywę topologiczną, która przekształciła tę dziedzinę od lat 80.

• Model Isinga — magnetyczne przejście fazowe — przeciągnij temperaturę przez T_c i obserwuj zarodkowanie domen
• Kondensat Bosego-Einsteina — schłodź poniżej temperatury krytycznej i obserwuj wzrost frakcji kondensatu
• Złącze Josephsona — zwiększaj prąd polaryzacji ponad prąd krytyczny i obserwuj skok napięcia
• Szkło spinowe — porównaj wpływ szybkości chłodzenia na zamrożoną konfigurację
• Izolator topologiczny (SSH) — dostrajaj t1/t2 w poprzek granicy fazowej i obserwuj pojawianie się stanów brzegowych

Myśl na zakończenie

Fizykę materii skondensowanej czasem nazywa się „fizyką wszystkiego innego” — czyli wszystkiego, co nie jest cząstką elementarną ani gwiazdą. Zjawiska omówione w tym przeglądzie nie są ciekawostkami zamkniętymi w laboratoriach ultraniskich temperatur. Złącza Josephsona leżą u podstaw najczulszych detektorów pola magnetycznego, jakie kiedykolwiek zbudowano (SQUID-y), wzorca napięcia używanego przez każdy krajowy instytut metrologii oraz architektur kubitów rozwijanych na potrzeby obliczeń kwantowych. Izolatory topologiczne są badane jako platformy dla odpornych na błędy obliczeń kwantowych za pomocą fermionów Majorany. Matematyka szkła spinowego pojawia się w teorii głębokiego uczenia. Model Isinga jest wykorzystywany do modelowania wszystkiego — od zwijania białek po krachy na rynkach finansowych. Zrozumienie tych symulacji to nie tylko fizyka dla samej fizyki — to fizyka, która kształtuje technologię nadchodzącej dekady.