PL 
EnglishDomyślny
Українська〰️ Fermi-Pasta-Ulam Recurrence
A chain of nonlinear oscillators showing energy returning to mode 1 — the FPU recurrence paradox.
Co pokazuje symulacja
W 1955 roku Enrico Fermi, John Pasta i Stanisław Ulam (wraz z programistką Mary Tsingou) przeprowadzili jedną z pierwszych symulacji komputerowych w fizyce. Spodziewali się, że łańcuch słabo nieliniowych oscylatorów stopniowo osiągnie równowagę termiczną — energia równomiernie rozłoży się po wszystkich modach normalnych. Zamiast tego, po odwiedzeniu kilku wyższych modów, energia niemal całkowicie wróciła do pierwszego modu: rekurencja FPU. Zaprzeczyło to panującemu założeniu, że słaba nieliniowość sama w sobie wystarcza do termalizacji.
Jak korzystać z symulacji
Ustaw β (nieliniowość) na zero, by uzyskać idealnie liniowy łańcuch — mod 1 zachowuje całą swoją energię na zawsze. Zwiększ β, by obserwować, jak energia przelewa się do modów 2, 3, 4… a następnie wraca. Zwiększ energię początkową E₀ dla silniejszej nieliniowości i szybszych (ale mniej dokładnych) rekurencji. Śledź więcej modów, by zobaczyć pełne widmo. Żółta przerywana linia oznacza poziom ekwipartycji — wysokość, jaką osiągnąłby każdy słupek przy prawdziwej termalizacji. Kliknij Mod 2, by zacząć od pobudzenia modu 2.
Czy wiesz, że…
Paradoks FPU został rozwiązany w latach 60. XX wieku dzięki odkryciu solitonów: zlokalizowanych fal nieliniowych, które przechodzą przez siebie nawzajem bez rozpadu. W łańcuchu FPU energia koncentruje się w strukturach zbliżonych do solitonów, które okresowo zwracają ją do modu początkowego. Przy wystarczająco wysokiej energii układ w końcu się termalizuje — ale w skalach czasowych znacznie dłuższych niż pierwotnie oczekiwano. Problem zapoczątkował całą dziedzinę dynamiki nieliniowej i układów całkowalnych.
O rekurencji Fermiego-Pasty-Ulama
Ta symulacja całkuje model β FPU: łańcuch N=32 cząstek połączonych sprężynami o liniowej i sześciennej (β) sile przywracającej. Równania ruchu ẍⅰ = (xᵢ₋₁ − xᵢ)(1 + β(xᵢ₋₁ − xᵢ)²) + (xᵢ₋₁ − xᵢ)(1 + β(xᵢ₋₁ − xᵢ)²) są całkowane schematem Rungego-Kutty 4. rzędu dla dokładności. Energia jest początkowo ładowana do pojedynczego modu normalnego poprzez przesunięcie cząstek zgodnie z jego sinusoidalnym wzorem funkcji własnej. Energie modalne są obliczane poprzez transformatę sinusoidalną przemieszczeń cząstek.
Trzy wizualizacje aktualizują się jednocześnie: łańcuch cząstek pokazuje rzeczywiste przemieszczenia, kolorowane według prędkości; widmo słupkowe pokazuje, jak energia jest rozłożona wśród pierwszych N modów (żółta przerywana linia = poziom ekwipartycji); a wykresy historii czasowej śledzą energię każdego modu w ciągu ostatnich 300 próbek czasu. Gdy β = 0, tylko mod 1 ma energię i tam pozostaje. Przy β > 0 energia przecieka do modów 2, 3, 4… a następnie wraca — rekurencja. Rekurencję wykrywa się, gdy mod 1 niesie >70% całkowitej energii.
Najczęściej zadawane pytania
Czym jest rekurencja Fermiego-Pasty-Ulama?
To zaskakujące zjawisko, w którym energia wstrzyknięta do najniższego modu normalnego łańcucha oscylatorów nieliniowych rozprzestrzenia się do wyższych modów, ale następnie wraca niemal w całości do pierwszego modu, zamiast rozłożyć się równomiernie (termalizować), jak przewiduje mechanika statystyczna.
Co kontroluje parametr β?
β ustala siłę sześciennej nieliniowości w sile sprężyny. Przy β=0 łańcuch jest idealnie liniowy, a energia nigdy nie opuszcza modu 1. Gdy β rośnie, sprzężenie między modami rośnie, energia rozprzestrzenia się szybciej, rekurencje stają się mniej dokładne, a przy bardzo wysokim β układ może ostatecznie się termalizować w długich skalach czasowych.
Czym jest energia modalna i jak jest obliczana?
Mod normalny to zbiorowy wzorzec oscylacji całego łańcucha. Jego energia jest obliczana poprzez rzutowanie przemieszczeń i prędkości cząstek na sinusoidalną funkcję własną modu za pomocą dyskretnej transformaty sinusoidalnej, a następnie obliczenie energii kinetycznej i potencjalnej w tej współrzędnej modu.
Dlaczego wynik FPU zaskoczył fizyków?
Klasyczna mechanika statystyczna i twierdzenie o ekwipartycji przewidują, że słabo sprzężone układy powinny osiągnąć równowagę termiczną, z energią równo rozłożoną między wszystkie stopnie swobody. Fermi spodziewał się, że nastąpi to szybko. Zamiast tego energia wciąż wracała do pierwszego modu, sugerując, że słaba nieliniowość nie gwarantuje termalizacji w dostępnych skalach czasowych.
Czym jest żółta przerywana linia na widmie?
Oznacza ona poziom energii ekwipartycji — wysokość, jaką osiągnąłby każdy słupek, gdyby całkowita energia była równo rozłożona między wszystkie śledzone mody. Prawdziwa termalizacja doprowadziłaby wszystkie słupki do tej linii. Rekurencja FPU sprawia, że słupek 1 okresowo jest znacznie wyższy od poziomu ekwipartycji.
Jak solitony wiążą się z paradoksem FPU?
W 1965 roku Zabusky i Kruskal pokazali, że ciągła granica łańcucha FPU to równanie Kortewega-de Vriesa (KdV), które dopuszcza rozwiązania solitonowe. Te zlokalizowane fale poruszają się bez dyspersji i przechodzą przez siebie sprężyście. W dyskretnym łańcuchu struktury zbliżone do solitonów niosą energię, nie dzieląc się nią, co wyjaśnia rekurencję.
Co dzieje się przy bardzo wysokiej energii lub dużym β?
Przy silnej nieliniowości układ ostatecznie wymyka się z niemal całkowalnego reżimu solitonowego. Rekurencje stają się niepełne i nieregularne, a w bardzo długich czasach symulacji układ rzeczywiście zbliża się do ekwipartycji — prawdziwej termalizacji. Energię przejścia nazywa się progiem FPU lub progiem stochastyczności.
Czy symulacja jest numerycznie dokładna?
Stosowana jest metoda RK4 (Runge-Kutta czwartego rzędu) z małym krokiem czasowym, zachowującą całkowitą energię hamiltonianu z wysoką precyzją. Można to zweryfikować, obserwując odczyt Energii całkowitej — powinien pozostawać niemal stały przez całą symulację.
Dlaczego N=32 cząstki?
Oryginalne obliczenia FPU z 1955 roku wykorzystywały N=32 lub N=64. Przy 32 cząstkach wewnętrznych okres rekurencji jest wystarczająco krótki, by był widoczny w czasie rzeczywistym, a struktura modalna jest łatwa do wyświetlenia. Mniej cząstek daje mniej modów i grubsze widmo; więcej cząstek spowalnia obliczenia.
Czym jest model β FPU w porównaniu z modelem α?
Model α dodaje kwadratowy człon nieliniowy (αΔx²) do siły sprężyny; użyty tu model β dodaje człon sześcienny (βΔx³). Model β jest prostszy do analizy, ponieważ zachowuje symetrię odbicia x → −x łańcucha. Oba modele wykazują zjawisko rekurencji, ale z różnymi okresami rekurencji i wzorcami sprzężenia modów.
Jak to się ma do współczesnej fizyki i matematyki?
Problem FPU zapoczątkował badania nad układami całkowalnymi, chaosem hamiltonowskim i fizyką obliczeniową jako dyscypliną. Wpłynął na twierdzenie KAM (Kołmogorowa-Arnolda-Mosera) o trwałości orbit quasi-okresowych i połączył się z badaniami nad dyskretnymi oddechami, falami łotra (rogue waves) oraz dynamiką nieliniowych sieci, wciąż aktywnymi dziś.