Застосування та синхронізація Куртових осциляторів |… |…

Модель та параметр порядку

У 1975 році Йошікі Курамото запропонував N фазових осциляторів θᵢ(t), кожен з власної природної частоти ωᵢ, обраної з розподілу g(ω), пов’язаних між собою через синусоїдну функцію. Синхронність вимірюється за допомогою комплексного параметра порядку r·e^(iψ), що дозволяє спрощувати рівняння до форми середнього поля, де кожен осцилятор реагує лише на колективне поле, а не на кожний з N-1 інших окремо:

dθᵢ/dt = ωᵢ + (K/N)·Σⱼ sin(θⱼ − θᵢ)

r·e^(iψ) = (1/N)·Σⱼ e^(iθⱼ)

(dθᵢ/dt = ωᵢ + K·r·sin(ψ − θᵢ))

Це спрощення робить модель Курамото точно розв’язуваною в межах 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

жива демонстрація · пов'язана симуляція● LIVE

Критичне з’єднання та фазова перехідність

Ниж критичного з’єднання K_c лише існує тривіальне рішення r = 0 — коли осцилятори змінюються незалежно. Вище K_c, з'являється нетривіальний синхронізований розгін, що росте як r ∝ √(K − K_c): класична фазова перехідність другого порядку. Розширення рівняння самоузгодженості для малих r дає критичне значення безпосередньо в термінах пікової щільності розподілу частот g(0):

K_c = 2 / (π·g(0)) Лорценовський g(ω) = γ/[π(ω²+γ²)] → K_c = 2γ r ∝ √(K − K_c) для K трохи вище K_c Вище K_c, осцилятори діляться на заблоковану групу (|ωᵢ| ≤ Kr, синхронізовану з середнім фазою) та групу, що змінюється, яка в середньому дає нульовий внесок — механізм, який вперше був підказаний моделлю ресета фази Артура Вінфри у 1967 році, яку Курамото спрощено до цієї аналітично обчислювальної форми.

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

Бджоли, пацієнти та нейрони

Східноазійські вогняні світлячки Pteroptyx malaccae синхронізують тисячі спалахів протягом 20 мс один за одним шляхом фазового просування або затримки до сусідніх – це стилі з’єднання синусоїдального типу Курамото у природі. Синусна вузлома (sinoatrial node) містить приблизно 10 000 кардіостимуляторних клітин із дещо різними внутрішніми частотами, які синхронізуються через з’єднання за отворами; коли з’єднання падає нижче K_c, результатом є аритмія. Гамма-генерація кору (30–80 Гц) залежить від синхронізації швидких інтернейронів, що виділяють імпульси, із параметром порядку r локальних потенціалів поля, який корелює з увагою та роботою пам’яті – а порушення синхронності пов’язане зі схильністю до шизофренії та епілепсії. Приблизно 20 000 нейронів годинника супречка (suprachiasmatic nucleus) навіть прогнозують час відновлення після перельотів, масштабування як 1/K.

Frequently asked questions

Що таке модель Курато́мо?

Модель Курато́мо описує N фазових осцилято́рів θᵢ(t) з природними частотами ωᵢ, підключених через dθᵢ/dt = ωᵢ + (K/N)Σⱼ sin(θⱼ − θᵢ). Кожен осцилято́р тяжіє до поточної фази інших; вище критичного значення коефіцієнту зв’язку популяція переходить від незв’язності до колективної синхронізації.

Що таке параметр порядку r у моделі Курато́мо?

Параметр порядку r·e^(iψ) = (1/N)Σⱼ e^(iθⱼ) вимірює синхронізацію: r ≈ 0 означає, що фази розподілені рівномірно (незв’язно), r ≈ 1 означає, що майже всі фази збігаються (повністю синхронізовано), а проміжні значення r описують часткову синхронізацію, коли деякі осцилято́ри заблоковані, а інші продовжують відхилятися.

Що таке критичний коефіцієнт зв’язку Kc?

Нижче Kc = 2/(π·g(0)), де g(0) є щільністю природних частот у середній точці популяції, популяція залишається незв’язною (r = 0). Вище Kc з'являється нетривіальне синхронізоване рішення, з r зростаючим як √(K − Kc) — класична друга-порядкова фазова перехід.

Спробуйте наживо

Усе, що вище, працює прямо у вашому браузері — відкрийте the simulation і змінюйте параметри під час роботи. Нічого не встановлюється, нічого не завантажується на сервер, уся модель живе в одній вкладці.

▶ Відкрити симуляцію the simulation