Powrót, Który Fermi Nie Oczekiwano

Eksperyment komputerowy, który miał być nudny

W 1955 roku Enrico Fermi, John Pasta, Stanislaw Ulam i Mary Tsingou przeprowadzili jeden z pierwszych eksperymentów komputerowych w fizyce na MANIAC I: łańcuch mas połączonych sprężynami z niewielkim nieliniowym korektą do prawa Hooke'a, wprowadzony w pojedynczym trybie o niskiej częstotliwości. Oczekiwano, zgodnie z mechaniką statystyczną, że nieliniowość powoli rozprowadzi energię równomiernie na wszystkie tryby drzeniowe łańcucha – proces nazywany termalizacją lub rozdzielnością energii. Zamiast tego energia przemieszczała się pomiędzy kilkoma niskimi trybami i, po długim ale skończonym czasie, niemal w całości wróciła do pierwotnego stanu pojedynczego trybu. Ta bliska rezygnacja, całkowicie nieoczekiwana, jest teraz znana jako problem Fermi-Pasta-Ulam-Tsingou (FPUT) i wyjaśnienie tego zjawiska zajęło fizykom ponad dziesięć lat.

Łańcuch i jego nieliniowość

Łańcuch FPUT składa się z N mas połączonych ze sobą sprężynami w sąsiednich odstępach, gdzie siła pomiędzy sąsiednimi masami zawiera zwykły człon liniowy oraz niewielkie korekcje kwadratowe lub sześcioelementowe:

m * d²xᵢ/dt² = k(x₁₊₁ - 2xᵢ + xᵢ⌿₁) + α[(x₁₊₁ - xᵢ)² - (xᵢ - xᵢ⌿₁)²] (model α, kwadratowy)  lub z członem sześcioelementowym β[(x₁₊₁ - xᵢ)³ - (xᵢ - xᵢ⌿₁)]³ (model β)

Gdy α = β = 0, jest to po prostu liniowy łańcuch sprzężonych oscylatorów harmonicznych, których normalne tryby nigdy nie wymieniają energii w ogóle – liniowy łańcuch rozpoczęty z jednego trybu pozostaje w nim na zawsze. Tylko nieliniowy człon pozwala na sprzężenie i wymianę energii między trybami, co jest dokładnie tym, dlaczego oczekiwanie ostatecznego pełnego termiczno-równowagowego stanu wydawało się naturalne, a jego brak (przynajmniej w obserwowanym czasie) był tak zaskakujący.

m · d²x_i/dt² = k(x_{i+1} − 2x_i + x_{i−1}) + α[(x_{i+1}−x_i)² − (x_i−x_{i−1})²]   (α-model, quadratic)
      or with a cubic term β[(x_{i+1}−x_i)³ − (x_i−x_{i−1})³]                              (β-model)

Rozwiązanie solitonów

Wyjaśnienie, znalezione przez Zabuskiego i Kruskala w 1965 roku, mówiło o tym, że w granicy ciągłości równania łańcucha FPUT przybliżają równanie Korteweg-de Vries (KdV), znane z możliwości wspierania solitonów: stabilnych, lokalizowanych impulsów fal, które przechodzą przez siebie podczas kolizji i zachowują swój kształt oraz prędkość po tym wszystkim, co jest sygnaturą specjalnej, prawie integrabelnej struktury ukrytej pod pozornym chaosem nieliniowym. Ponieważ dynamika leżąca w jej rdzeniu jest zbliżona do systemu integrabelnego (czyli takiego z tak wieloma zachowanymi ilościami jak stopni swobody), trajektorie energii w przestrzeni modalnej pozostają ograniczone do quasi-okresowych orbit zamiast rozprzestrzeniać się ergodycznie — co dokładnie obserwował zespół Fermiego, zamiast rzeczywistej termalizacji.

Toda Lattice: Dokładnie Rozwiązywalny Brat

Około tego samego okresu, Morikazu Toda odkrył specyficzną siłę nieliniową, wykładniczą a nie wielomianową, dla której łańcuch jest dokładnie integrujący się — posiada on tyle samo pormanów co cząstek, a jego rozwiązania mogą być zapisane w zamkniętej postaci za pomocą rozwiązań solitonowych, które nigdy nie ulegają rozpadowi. Toda Lattice zajmuje miejsce obok FPUT: wielomianowy łańcuch FPUT nie jest dokładnie integrujący się, ale dla wystarczająco małej nieliniowości zachowuje się niemal tak, jakby to robił, przez wykładniczo długi czas, co stanowi współczesne rozumienie, dlaczego pojawia się rekurencja a nie natychmiastowe termiczne zrównoważenie.

Chaos w falach versus prawdziwa równomierność

Zwiększ siłę nieliniowości α lub β wystarczająco wysoko, albo poczekaj wystarczająco długo przy umiarkowanej sile, a łańcuch FPUT ostatecznie odchodzi od quasi-okresowych powtórzeń i ustala się w prawdziwej równomierności – zachowanie umowne związane z pojawieniem się dużych fal chaosu w przestrzeni faz amplitudy modów, gdzie efektywne wycentrowania Lyapunova ułamki stają się dodatnie. Przejście od uporządkowanego kapania energii do prawdziwego chaosu falowego, przy zmianie jednego parametru, jest jednym z najbardziej czystych demonstracji laboratoryjnych, jak deterministyczne układy nieliniowe przechodzą z zachowania bliskiego integracji w pełne termiczne zróżnicowanie.

Często zadawane pytania

Wypróbuj na żywo

Wszystko powyżej działa bezpośrednio w Twojej przeglądarce — otwórz Chaotic Waves i zmieniaj parametry podczas działania. Nic nie jest instalowane ani przesyłane na serwer, cały model działa w jednej karcie.