🔄 Synchronizacja Kuramoto

Model Kuramoto symuluje dużą populację sprzężonych oscylatorów, z których każdy ma nieco inną częstotliwość własną, a które dążą do synchronizacji, gdy wzajemne sprzężenie jest wystarczająco silne. Uchwytuje przejście fazowe od nieuporządkowanego (asynchronicznego) do uporządkowanego (synchronicznego) zachowania zbiorowego, będąc modelowym przykładem synchronizacji w układach złożonych.

Parametr porządku r mierzy stopień synchronizacji. Jest to wartość bezwzględna średniego wektora fazowego jednostkowego wszystkich oscylatorów: r = |średnia(exp(i·θⱼ))|. Gdy r = 0, oscylatory są równomiernie rozłożone wokół okręgu (niespójne). Gdy r = 1, wszystkie oscylatory mają tę samą fazę (w pełni zsynchronizowane). Przejście od r ≈ 0 do r > 0 następuje gwałtownie powyżej krytycznego sprzężenia Kc.

Synchronizacja napędzana sprzężeniem podobnym do modelu Kuramoto występuje między innymi w: synchronizacji błysków świetlików (południowoazjatyckie świetliki Pteroptyx wytwarzają idealnie zsynchronizowane błyski), komórkach rozrusznikowych serca w węźle zatokowo-przedsionkowym (miliony komórek bijących razem), oscylacjach neuronalnych w mózgu (oscylacje gamma koordynują przetwarzanie sensoryczne) oraz synchronizacji cykli miesiączkowych w grupach społecznych (choć dowody na to ostatnie zjawisko są dyskusyjne).

Generatory prądu podłączone do sieci elektroenergetycznej muszą wirować z tą samą częstotliwością (50 Hz w Europie, 60 Hz w Ameryce Północnej), aby uniknąć destrukcyjnej interferencji. Generatory są sprzężone za pośrednictwem wspólnych linii przesyłowych, a dynamikę ich synchronizacji fazowej opisuje wariant modelu Kuramoto. Utrata synchronizacji — gdy faza generatora się przesuwa — powoduje odłączenie go przez przekaźniki zabezpieczające, co może wywołać kaskadowe awarie sieci.

Sprzężenie krytyczne Kc to minimalna siła sprzężenia, przy której układ może spontanicznie się zsynchronizować. Dla N oscylatorów o częstotliwościach własnych rozłożonych zgodnie z rozkładem Lorentza o połówkowej szerokości γ, Kc = 2γ. Poniżej Kc stan niespójny r = 0 jest stabilny dla wszystkich warunków początkowych. Powyżej Kc synchronizacja rośnie jak r ~ √(1 − Kc/K) tuż powyżej przejścia, analogicznie do parametru porządku w pobliżu przejścia fazowego drugiego rodzaju.