Synchronizacja Kuramoto: Matematyka stojąca za zachowaniem się motyli i rytmem serca

Model i parametr uporządkowania

W 1975 roku Yoshiki Kuramoto zaproponował N oscylatorów fazowych θᵢ(t), każdy z naturalną częstotliwością ωᵢ wybraną z rozkładu g(ω), połączonych ze sobą w sposób wszechnobazowy za pomocą funkcji sinusowej. Uporządkowanie jest mierzone przez parametr zespolony r·e^(iψ), który pozwala skrócić równanie do postaci średniej pola, w której każdy oscylator reaguje jedynie na zbiorowe pole, a nie indywidualnie na każdy z N-1 innych:

dθᵢ/dt = ωᵢ + (K/N)·Σⱼ sin(θⱼ − θᵢ) ← oryginalna postać N-ciałowa r·e^(iψ) = (1/N)·Σⱼ e^(iθⱼ) ← parametr uporządkowania dθᵢ/dt = ωᵢ + K·r·sin(ψ − θᵢ) ← równoważna postać średniego pola r ≈ 0 oznacza, że fazy rozrzucają się równomiernie po okręgu (nierozłączność); r ≈ 1 oznacza, że prawie każdy oscylator dzieli tę samą fazę (pełne uporządkowanie). To uproszczenie sprawia, że model Kuramoto jest dokładnie rozwiązywalny w granicach N → ∞.

dθᵢ/dt = ωᵢ + (K/N)·Σⱼ sin(θⱼ − θᵢ)     ← original N-body form
r·e^(iψ) = (1/N)·Σⱼ e^(iθⱼ)              ← order parameter
dθᵢ/dt = ωᵢ + K·r·sin(ψ − θᵢ)            ← equivalent mean-field form

Krytyczne sprzężenie i przejście fazowe

Poniżej krytycznego sprzężenia K_c, występuje tylko trywialne rozwiązanie r = 0 – oscylatory dryfują niezależnie. Powyżej K_c, pojawia się nienzętowy, zsynchronizowany odcinek rosnący w miarę jak r ∝ √(K − K_c): klasyczny przykład drugiego rzędu przejścia fazowego. Rozszerzając równanie samo-spójności dla małych r, otrzymujemy krytyczną wartość bezpośrednio w terminach gęstości szczytowej funkcji rozkładu częstotliwości g(0):

K_c = 2 / (π·g(0)) Funkcja Lorentza g(ω) = γ/[π(ω²+γ²)] → K_c = 2γ r ∝ √(K − K_c) dla K nieco powyżej K_c Powyżej K_c, oscylatory dzielą się na zablokowaną grupę (|ωᵢ| ≤ Kr, zsynchronizowaną z średnią fazą) i grupę dryfującą, która średnio daje zerowy wkład – mechanizm, który po raz pierwszy zauważył Arthur Winfree w jego modelu resetowania fazowego z 1967 roku, uproszczony przez Kuratowa do tej analitycznie wykonalnej postaci.

Kc = 2 / (π·g(0))
Lorentzian g(ω) = γ/[π(ω²+γ²)]  →  Kc = 2γ
r ∝ √(K − Kc)   for K slightly above Kc

Buzki, rozruszniki i neurony

Azjatyckie larwy gatunku Pteroptyx malaccae synchronizują tysiące błysków w ciągu 20 ms, wzajemnie się opóźniając lub przyspieszając fazowo – to przykład synchronizacji Kuramoto w natywnym środowisku. Sionoatrialny rozrusznik serca zawiera około 10 000 komórek pacemakera o lekko różniących się podstawowych rytmach, które synchronizują się za pomocą łączności szczelinowej; gdy współczynnik K_c spada poniżej wartości, obserwuje się migotanie przedsionków. Gamma oscylacje korowe (30–80 Hz) opierają się na synchronizacji szybkich interneuronów typu fast-spiking, a parametr r lokalnych potencjałów pola korelują z uwagą i pamięcią roboczą – a zaburzona synchronizacja jest związana ze schizofrenią i padaczkami. W około 20 000 neuronach zegara okołodziennego w jądrze nadjądrowe nawet prognozuje czas powrotu po jetlagu, skalując się jako 1/K.

Frequently asked questions

Wypróbuj na żywo

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