Rytm mózgu to miliony zegarów zgadzających się na wspólne tyknięcie
Elektroencefalogram (EEG) nie rejestruje pojedynczych neuronów – rejestruje sumowe pole elektryczne milionów z nich, a powód, dla którego ta suma pokazuje czyste rytmiczne fale zamiast szumu, to fakt, że duże populacje neuronów synchronizują swoje pobudzenia. Modelowanie dokładnie tego, jak i kiedy populacja wchodzi w synchronię, to zadanie modeli oscylatorów sprzężonych, a najczystszym z nich jest ten wprowadzony przez Yoshiki Kuramoto w 1975 roku.
W modelu Kuramoto każdy oscylator i ma własną naturalną częstotliwość ωᵢ (wybraną z rozkładu odzwierciedlającego różnorodność populacji) oraz fazę θᵢ, która postępuje zgodnie z tą częstotliwością, ale jest popychania przez każdy inny oscylator w populacji, ważona przez współczynnik sprzężenia K:
dθᵢ/dt = ωᵢ + (K/N) · Σⱼ sin(θⱼ − θᵢ) θᵢ = faza oscylatora i ωᵢ = jego naturalna częstotliwość K = współczynnik sprzężenia N = liczba oscylatorów parametr uporządkowania: r·e^(iψ) = (1/N) · Σⱼ e^(iθⱼ) (r ∈ [0,1] mierzy synchronizację) Poniżej krytycznego współczynnika sprzężenia Kc fazy oscylatorów pozostają rozproszone i parametr uporządkowania r waha się blisko zera – niepokój. Powyżej Kc grupa oscylatorów zamyka się razem, a r gwałtownie rośnie do 1: przejście fazowe synchronizacji, matematycznie takie samo nagłe przejście kolektywne jakby obserwowane w błyskach jednoczesnych świerszczy, oklaskujących tłumów wpadających w rytm i generatorach sieci energetycznej utrzymujących synchronizację.
dθᵢ/dt = ωᵢ + (K/N) · Σⱼ sin(θⱼ − θᵢ) θᵢ = phase of oscillator i ωᵢ = its natural frequency K = coupling strength N = number of oscillators order parameter: r·e^(iψ) = (1/N) · Σⱼ e^(iθⱼ) (r ∈ [0,1] measures synchrony)
Od abstrakcyjnych faz do nazwanych rytmów mózgowych
Rzeczywiste rytmy korowe konwencjonalnie dzielą się na pasma częstotliwości, z których każde jest związane z różnymi stanami poznawczymi: delta (poniżej 4 Hz, głęboki sen), teta (4–8 Hz, kodowanie pamięci i nawigacja), alfa (8–13 Hz, zrelaksowany czuwanie, szczególnie nad płatem wzrokowym z zamkniętymi oczami), beta (13–30 Hz, aktywna koncentracja i planowanie motoryczne) i gamma (30–100+ Hz, związane z łączeniem cech w pojedyncze postrzeganie oraz z przetwarzaniem sensorycznym). Sieć Kuramoto dostrojona do średniej częstotliwości naturalnej zbliżonej do, powiedzmy, 40 Hz i podwyższona powyżej progu synchronizacji, generuje koherentną oscylację, która w całości przypomina rytm gamma.
Podostymulacja i hamowanie ustalają prędkość rytmu
Biologicznie, rytmy korowe nie wyłaniają się z identycznych oscylatorów fazowych połączonych równomiernie – powstają w wyniku interakcji populacji neuronów pobudzających (E) i hamujących (I). Kanonicznym mechanizmem dla gamma jest PING (sieć piramidowo-interneuronowa gamma): komórki piramidalne pobudzone napędzają inhibitory interneurony, które z kolei przytłaczają sieć; gdy hamowanie zanika, podostymulacja buduje się ponownie aż do przerwania, a cykl powtarza się w okresie głównie ustalonym przez stały czas hamowania. To równowaga E-I jest przyczyną tego, że rytmy gamma są szybsze niż theta – szybki upadek synaptycznego hamowania w układzie E-I powoduje szybkie cyklowanie pętli, podczas gdy wolniejsze pętli z powtarzającym się pobudzeniem i hamowaniem zaangażowane w bardziej odległe obwody produkują wolniejsze rytmy theta i alfa.
Dlaczego synchroniczność to więcej niż ciekawość
Hipoteza komunikacji poprzez spójność postuluje, że zsynchronizowane fale gamma działają jako mechanizm ustalania czasu, otwierając krótkie okna, w czasie których dwie obszary mózgowe oscylujące w fazie mogą skutecznie wymieniać informacje, podczas gdy obszary poza fazą mówią do siebie bez sensu. Zaburzony równowaga E-I i nieprawidłowa synchroniczność stanowią również istotne klinicznie – padaczki to, na jednym poziomie, patologiczna nadmierna synchronizacja dużych populacji neuronów, a takie schorzenia jak psychoza są związane z obniżoną synchronizacją gamma podczas zadań wymagających łączenia cech sensorycznych.
Frequently asked questions
Co mierzy parametr uporządkowania Kuramoto?
Jest to pojedyncza liczba r pomiędzy 0 a 1 podsumowująca stopień, w jakim fazy wszystkich oscylatorów są ułożone w danym momencie – r bliskie 0 oznacza, że fazy są losowo rozrzucone po okręgu (nieuporządkowanie), a r bliskie 1 oznacza, że prawie wszystkie oscylatory dzielą tę samą fazę (pełne zsynchronizowanie).
Dlaczego fale gamma są szybsze niż fale theta?
Fale gamma powstają głównie dzięki szybkim pętlom sprzężenia zwrotnego excitatory-inhibitory (mechanizm PING), których okres jest ustalany przez szybkie opóźnienie synaptyczne, podczas gdy fale theta zależą od wolniejszych, bardziej rozproszonych pętli obwodów – różne stałe czasowe prowadzą do różnych charakterystycznych częstotliwości.
Czy sieć neuronów o różnych naturalnych częstotliwościach może zsynchronizować się?
Tak - to dokładnie to demonstruje model Kuramoto. O ile siła sprzężenia między oscylatorami przekracza krytyczny próg, zależny od rozrzutu ich naturalnych częstotliwości, duża część populacji zablokuje się w jednej wspólnej fazie, mimo że każdy indywidualny oscylator w izolacji obracałby się we własnym tempie.
Wypróbuj na żywo
Wszystko powyżej działa bezpośrednio w Twojej przeglądarce — otwórz Neural Oscillators i zmieniaj parametry podczas działania. Nic nie jest instalowane ani przesyłane na serwer, cały model działa w jednej karcie.
▶ Otwórz symulację Neural Oscillators