Jest w tym pewna osobliwa satysfakcja — obserwować, jak ziarna piasku się osypują. Upuść jedno ziarno na stos piasku Abela i najczęściej nic dramatycznego się nie stanie. Ale czasem to pojedyncze ziarno wyzwala lawinę, która rozprzestrzenia się po całej siatce — kaskadę, której rozkład rozmiarów podlega prawu potęgowemu bez charakterystycznej skali. To samoorganizująca się krytyczność (SOC): system, który sam ustawia się na krawędzi chaosu bez żadnego zewnętrznego dostrajania. Zaproponowali ją Per Bak, Chao Tang i Kurt Wiesenfeld w 1987 roku i pozostaje ona jedną z najciszej rewolucyjnych idei fizyki statystycznej. Fala 105 wprowadza SOC do Twojej przeglądarki — wraz z 29 innymi symulacjami, równie mało prawdopodobnymi do spotkania na standardowym kursie licencjackim.
Fizyka na krawędzi porządku: SOC, szkło spinowe i złącza Josephsona
Stos piasku Abela i samoorganizująca się krytyczność
Model stosu piasku Abela żyje na dwuwymiarowej siatce całkowitoliczbowej. Każda komórka przechowuje pewną liczbę ziaren. Gdy komórka przekroczy próg — zwykle 3 dla siatki kwadratowej — osypuje się, rozdzielając po jednym ziarnie do każdego z czterech sąsiadów:
if grid[x][y] >= 4:
grid[x][y] -= 4
grid[x+1][y] += 1
grid[x-1][y] += 1
grid[x][y+1] += 1
grid[x][y-1] += 1
Model nazywa się abelowym, ponieważ ostateczny stan stabilny nie zależy od kolejności, w jakiej zachodzą osypania — nieoczywisty fakt algebraiczny, który nadaje systemowi bogatą strukturę grupową. Emergentne wzory fraktalne powstające przy zrzucaniu milionów ziaren na jedną centralną komórkę są uderzające: zagnieżdżone trójkąty, pierścienie w kształcie rombów i drobnoziarniste prążki interferencyjne wyglądające niemal krystalicznie, mimo że powstają z czysto lokalnych reguł. Nasza symulacja pozwala interaktywnie zasiewać stos i obserwować w czasie rzeczywistym, jak gromadzi się statystyka lawin.
Szkło spinowe: frustracja i nieporządek
Szkło spinowe to materiał magnetyczny, w którym konkurujące oddziaływania ferromagnetyczne
i antyferromagnetyczne tworzą frustrację — żadna konfiguracja spinów nie może
jednocześnie spełnić wszystkich wiązań. Model Edwardsa-Andersona odwzorowuje to na
siatce 2D, przypisując każdemu wiązaniu Jij losowy znak wybrany
z ±1. Energia konfiguracji wynosi:
E = -∑_⟨ij⟩ J_ij · s_i · s_j
gdzie suma przebiega po parach najbliższych sąsiadów. Minimalizacja tej energii w niskiej temperaturze jest notorycznie trudna — problem jest NP-zupełny w trzech wymiarach. Symulacja wykorzystuje metodę Monte Carlo Metropolisa-Hastingsa z powolnym harmonogramem wyżarzania, dzięki czemu można zaobserwować, jak system zamarza w jednym ze swoich licznych stanów metastabilnych. Szkła spinowe są koncepcyjnym przodkiem sieci Hopfielda i współczesnych krajobrazów energii w uczeniu maszynowym, więc to nie jest wyłącznie akademicka ciekawostka.
Złącze Josephsona: model RSJ
Gdy dwa nadprzewodniki oddziela cienka bariera izolacyjna, pary Coopera tunelują
kwantowo-mechanicznie — to efekt Josephsona. Model złącza rezystywnie bocznikowanego
(RSJ) traktuje złącze jako kondensator, rezystor i nieliniową indukcyjność połączone
równolegle. Dynamika cechowniczo niezmienniczej różnicy faz φ na złączu
podlega równaniu:
C · d²φ/dt² + (1/R) · dφ/dt + I_c · sin(φ) = I_bias
To równanie jest formalnie identyczne z równaniem wymuszonego, tłumionego wahadła — fakt,
który czyni je idealnym pedagogicznym pomostem między mechaniką klasyczną a transportem
kwantowym. Powyżej prądu krytycznego I_c złącze przełącza się ze stanu
zerowego napięcia w stan oscylujący, emitując promieniowanie mikrofalowe o częstotliwości
f = 2eV/h — podstawę wzorca napięcia Josephsona, najdokładniejszej realizacji
wolta w metrologii. Symulacja rysuje jednocześnie krzywą I-V i portret przestrzeni
fazowej.
Materia kwantowa i fizyka ultraniskich temperatur
Kondensat Bosego-Einsteina: równanie Grossa-Pitajewskiego
Gdy rozrzedzony gaz bozonów zostanie schłodzony do temperatury zaledwie nanokelwinów
powyżej zera absolutnego, makroskopowa frakcja atomów zapada się do tego samego kwantowego
stanu podstawowego. Zbiorowe zachowanie tego kondensatu opisuje pojedyncza zespolona
funkcja falowa ψ(r, t), rządzona równaniem Grossa-Pitajewskiego (GPE):
iℏ ∂ψ/∂t = [-ℏ²/(2m) ∇² + V(r) + g|ψ|²] ψ
Nieliniowy człon g|ψ|² koduje oddziaływania kontaktowe pola średniego
między atomami. GPE dopuszcza skwantowane wiry — defekty topologiczne, w których faza
ψ nawija się o 2π — będące znakiem rozpoznawczym
nadciekłości. Nasza symulacja 2D rozwiązuje GPE metodą operatora rozdzielonego z
krokowym całkowaniem w czasie na jednorodnej siatce, umożliwiając obserwację nukleacji
wirów, propagacji solitonów oraz prążków interferencyjnych pojawiających się, gdy dwa
kondensaty się nakładają. To ta sama fizyka, za którą Cornell, Ketterle i Wieman
otrzymali Nagrodę Nobla w 2001 roku — a działa ona w karcie Twojej przeglądarki.
Wszystkie prawdopodobieństwa funkcji falowej są renderowane w percepcyjnie jednorodnych mapach kolorów, gdzie faza jest zakodowana odcieniem, a gęstość — jasnością, co daje natychmiastowy wizualny odczyt stanu kwantowego bez żadnej wcześniejszej wiedzy o liczbach zespolonych.
Geometria, fraktale i psychoakustyka
Mozaika Penrose'a metodą inflacji
Mozaiki Penrose'a pokrywają płaszczyznę z pięciokrotną symetrią obrotową — wzorem
kategorycznie zabronionym przez klasyczne twierdzenie o ograniczeniu krystalograficznym.
Odkrył je Roger Penrose w 1974 roku i antycypowały one eksperymentalne odkrycie
kwazikryształów przez Dana Shechtmana dekadę później (Nagroda Nobla z chemii 2011).
Algorytm inflacji generuje mozaikę rekurencyjnie: zaczynając od małego fragmentu
kafelków typu latawiec i strzała, powtarzalnie dzieli każdy kafelek według ustalonych
reguł podstawienia, skalując siatkę na każdym kroku o złoty podział
φ = (1+√5)/2.
Symulacja wykonuje algorytm inflacji w przeglądarce i pozwala przybliżać samopodobną strukturę. Ponieważ mozaika jest nieokresowa, nigdy nie znajdziesz translacyjnej jednostki powtarzalnej — a mimo to nie jest przypadkowa. Ma dalekozasięgowy porządek orientacyjny, widoczny w jej obrazie dyfrakcyjnym, który pokazuje ostre piki Bragga ułożone z symetrią ikozaedryczną.
Uszczelka Apolloniusza: rekurencyjne upakowanie okręgów
Mając trzy wzajemnie styczne okręgi, Apoloniusz z Pergi pokazał już w 200 roku p.n.e., że istnieją dokładnie dwa okręgi styczne do wszystkich trzech. Jeśli powtarzalnie zastosujesz tę konstrukcję do każdej nowej trójki stycznych okręgów, otrzymasz uszczelkę Apolloniusza — fraktal o wymiarze Hausdorffa równym w przybliżeniu 1,3057. Krzywizny wszystkich okręgów w dowolnym całkowitoliczbowym upakowaniu Apolloniusza spełniają twierdzenie Kartezjusza o okręgach:
(k₁ + k₂ + k₃ + k₄)² = 2(k₁² + k₂² + k₃² + k₄²)
gdzie k = 1/r to krzywizna (odwrotność promienia). Jeśli cztery początkowe
krzywizny są liczbami całkowitymi, każdy okrąg w uszczelce ma całkowitą krzywiznę — teorio-
liczbowy cud, który łączy geometrię klasyczną z teorią form kwadratowych oraz grupą
apolloniuszową w arytmetyce. Symulacja renderuje tysiące okręgów w jednej klatce,
wykorzystując rekurencyjne zstępowanie z wczesnym odrzucaniem okręgów mniejszych niż
rozmiar piksela.
Ton Sheparda: nieskończone schody dźwięku
Ton Sheparda to iluzja słuchowa: dźwięk, który sprawia wrażenie, jakby wznosił się (lub opadał) w wysokości w nieskończoność, choć w rzeczywistości nigdy nie staje się wyższy ani niższy. Powstaje przez nałożenie fal sinusoidalnych w odstępach oktawowych, z dzwonową obwiednią amplitudy w przestrzeni logarytmicznej częstotliwości. Gdy każdy składnik pojawia się od dołu widma i zanika u góry, odbierana wysokość dźwięku spiralnie rośnie bez ograniczeń — akustyczny odpowiednik schodów Eschera.
Symulacja wizualizuje obwiednię widmową i chwilową fazę każdego składnika częściowego w czasie rzeczywistym, wykorzystując Web Audio API i wyświetlacz kaskadowy Canvas 2D. Możesz regulować tempo wznoszenia, szerokość pasma obwiedni oraz liczbę składników częściowych, co czyni mechanizm iluzji przejrzystym. Tony Sheparda pojawiają się w muzyce filmowej — zwłaszcza w filmach Christophera Nolana Dunkierk i Mroczny Rycerz Powstaje — tworząc podprogowe poczucie nieubłaganej eskalacji.
Uwaga o renderowaniu: wszystkie 30 symulacji fali 105 celuje w Canvas 2D, a nie WebGL. Była to świadoma decyzja: Canvas 2D działa wszędzie, w tym na tabletach o niskiej mocy i starszych telefonach z Androidem, a dla symulacji opartych na siatce i geometrii API 2D jest często szybsze w iteracji niż kod shaderów. Budżet wydajności wynosi 16 ms na klatkę; symulacje przekraczające go na urządzeniu średniej klasy automatycznie przechodzą na zmniejszoną rozdzielczość siatki.
Wypróbuj sam
Najlepszym sposobem na wyrobienie sobie intuicji dla tych modeli jest ich zepsucie. Spróbuj podkręcić prąd polaryzacji powyżej punktu krytycznego w złączu Josephsona i zobacz, jak krzywa I-V przeskakuje z zerowego napięcia do skończonego. Zasiej stos piasku asymetrycznie i sprawdź, czy granica fraktalna pozostaje symetryczna. Zamieszaj kondensat ruchomym potencjałem przeszkody i policz nukleujące się skwantowane wiry. Oto trzy symulacje, od których warto zacząć:
Pełny katalog fali 105 — wszystkie 30 symulacji — można przeglądać przez indeks kategorii. Każda strona symulacji zawiera przegląd matematyczny, kontrolki parametrów z podpisanymi jednostkami oraz odnośnik do literatury źródłowej, gdyby samo okno przeglądarki nie wystarczyło.
Co dalej
Fala 105 zamyka obecny sprint. Modele z tej fali leżą na przecięciu fizyki materii skondensowanej, teorii liczb i nauki o percepcji — dziedzin, które rzadko dzielą wspólny program nauczania, lecz mają ze sobą więcej wspólnej matematyki, niż zwykle przyznają ich praktycy. Równanie złącza Josephsona i wahadło to to samo równanie różniczkowe. Uszczelka Apolloniusza i teoria form automorficznych są, w ścisłym sensie, tą samą strukturą. Grupa stosu piasku i algebra abstrakcyjna są bliższymi krewnymi, niż spodziewałby się ktoś, kto uczył się każdej z nich na osobnym kursie.
Fala 106 jest już w planach. Docelowe tematy to systemy reakcji-dyfuzji na powierzchniach nieeuklidesowych, modele metapopulacyjne epidemiologii z rzeczywistymi danymi o mobilności oraz powrót do podwójnego wahadła z mapami cieplnymi wykładnika Lapunowa. Jeśli masz symulację, którą chciałbyś zobaczyć — albo model, który Twoim zdaniem jest niedostatecznie reprezentowany w edukacji naukowej opartej na przeglądarce — strona kontaktowa jest zawsze otwarta. Celem nie jest odtworzenie tego, co można znaleźć w podręczniku. Celem jest sprawić, by rzeczy, które podręczniki zostawiają jako ćwiczenia, naprawdę działały, na Twoich oczach, już teraz.