Spotlight #56 – Chaos, Ekologia i Neuronauka

Teoria chaosu to nauka o systemach silnie wrażliwych na warunki początkowe — niewielkie różnice w stanie startowym rosną wykładniczo, czyniąc długoterminowe przewidywanie niemożliwym. A jednak systemy te nie są losowe: podążają za dokładnymi deterministycznymi równaniami, a ich trajektorie są ograniczone do geometrycznie pięknego obiektu: dziwnego atraktora.

Efekt motyla

Edward Lorenz odkrył chaos deterministyczny w 1963 roku, uruchamiając uproszczony model konwekcji atmosferycznej. Zauważył, że ponowne uruchomienie symulacji z wydruku (zaokrąglonego do 3 miejsc po przecinku zamiast 6) po krótkim czasie dawało zupełnie inną prognozę. Trzy równania, które to kodują, to:

ẋ = σ(y − x)
ẏ = x(ρ − z) − y
ż = xy − βz

Przy σ=10, ρ=28, β=8/3 trajektoria kreśli słynnego dwupłatowego motyla, nigdy się nie powtarzając, wirując wokół dwóch niestabilnych punktów stałych w nieskończoność. Największy wykładnik Lapunowa wynosi λ≈0,905, co oznacza, że bliskie trajektorie rozchodzą się o czynnik e w każdej jednostce czasu.

Geometria fraktalna

Atraktor Lorenza ma wymiar Hausdorffa około 2,06 — więcej niż powierzchnia, ale mniej niż bryła. Ten ułamkowy wymiar jest geometrycznym podpisem chaosu: atraktor jest nieskończenie sfałdowany, ale zajmuje zerową objętość. Atraktor Thomasa przy b≈0,19 pokazuje uderzającą trójkrotnie symetryczną strukturę, najlepiej badaną poprzez przeciąganie i oglądanie z różnych kątów.

Odkryj w bibliotece:

🌀 Dziwne atraktory 🌿 Diagram bifurkacyjny 🎲 Gra w chaos 🦋 System Lorenza

Ekologia populacyjna bada, jak gatunki wchodzą w interakcje, konkurują i współewoluują w ekosystemach. Równania Lotki-Volterry (niezależnie wyprowadzone przez Alfreda Lotkę w 1925 roku i Vito Volterrę w 1926 roku) ujmują istotę dynamiki drapieżnik-ofiara:

dH/dt = rH − aHL        (ofiara: rosnie w tempie r, zjadana w tempie a na drapieznika)
dL/dt = baHL − mL      (drapieznicy: rozmnazaja sie kosztem ofiary, gina w tempie m)

Rozwiązaniem są okresowe oscylacje: populacja ofiar osiąga szczyt, drapieżnicy pozostają w tyle i też osiągają szczyt, ofiary załamują się, drapieżnicy głodują i również się załamują, ofiary się odbudowują. Okres i amplituda zależą od warunków początkowych — w przeciwieństwie do modelu Kuramoto, nie ma tu ustalonego atraktora.

Dlaczego modele agentowe mają znaczenie

Wersja równania różniczkowego zakłada nieskończone, idealnie wymieszane populacje. Prawdziwe ekosystemy są skończone i przestrzenne. Symulacja wilk-owca pokazuje, że struktura przestrzenna tworzy refugia — łaty, w których ofiary przetrwają fale wilków — co może ustabilizować to, co w przeciwnym razie byłoby rozbieżnymi oscylacjami lub kaskadami wymierania. Szum demograficzny stochastyczny (indywidualne losowe zgony i narodziny) również daje jakościowo inne zachowanie niż gładkie ODE.

Wypróbuj ten eksperyment

1. Zacznij od domyślnych parametrów i obserwuj 2-3 cykle populacyjne.
2. Zwiększ zysk pokarmowy wilków do 70. Obserwuj, jak wilki dominują, a owce wymierają.
3. Użyj przycisku Więcej owiec, by zresetować do warunku początkowego z przewagą owiec, i porównaj przebiegi przejściowe.
4. Zmniejsz odrastanie trawy do 1. Czy ekosystem może się utrzymać?

Odkryj w bibliotece:

🐺 Drapieżnictwo wilk-owca 🐦 Stado Boids 🪸 Rafa koralowa 🦠 Kolonia bakterii

Mózg jest pełen oscylujących neuronów: oscylacje gamma (40 Hz) leżą u podstaw uwagi i świadomości, rytmy theta (4–8 Hz) są powiązane z pamięcią i nawigacją, fale alfa (8–12 Hz) odzwierciedlają synchronię korową w stanie spoczynku. Jak populacja neuronów o różnorodnych częstotliwościach wewnętrznych osiąga zbiorową oscylację? Model Kuramoto daje matematyczną odpowiedź: przejście fazowe.

Od pojedynczych oscylatorów do zbiorowego rytmu

W symulacji każdy oscylator reprezentuje neuron (lub populację neuronalną) z własnym naturalnym tempem wyładowań ωᵢ. Wyraz sprzężenia (K/N)Σsin(φⱼ−φᵢ) działa jak pobudzenie synaptyczne: neurony starają się wyładowywać w fazie z sąsiadami. Poniżej K_c sieć jest asynchroniczna — każdy neuron wyładowuje się we własnym tempie, a średnia populacji się znosi. Powyżej K_c pojawia się makroskopowa oscylacja, której nie ma w programie żadnego pojedynczego neuronu.

Przejawy biologiczne

• Zegar dobowy — neurony w jądrze nadskrzyżowaniowym (SCN) synchronizują się, tworząc 24-godzinny zegar biologiczny. Jet lag to zdarzenie desynchronizacji.
• Arytmia serca — komórki rozrusznikowe tracą synchronię podczas migotania przedsionków; defibrylacja to wymuszona resynchronizacja.
• Choroba Parkinsona — patologiczna nadmierna synchronia w obwodach dopaminowych (blokowanie pasma beta) jest obecnie celem głębokiej stymulacji mózgu.
• Świetliki — Pteroptyx malaccae synchronizują błyskanie w całych drzewach, makroskopowe zjawisko Kuramoto widoczne gołym okiem.

Odkryj w bibliotece:

🔄 Synchronizacja Kuramoto 🧠 Oscylacje fal mózgowych ❤️ Potencjał czynnościowy serca 🕐 Rytm dobowy