Химерні стани

1. Курамото та Баттогтох, 2002

Модель Курамото описує популяції зв'язаних фазових осциляторів і те, як вони переходять до повної синхронії зі зростанням сили зв'язку. Оригінальна модель Курамото використовувала глобальний (усі-з-усіма) зв'язок: кожен осцилятор однаково відчуває вплив кожного іншого.

У 2002 році Курамото та Баттогтох досліджували кільце ідентичних фазових осциляторів із нелокальним зв'язком — кожен осцилятор перебуває головним чином під впливом сусідів, причому вплив плавно спадає з відстанню вздовж кільця. Вони очікували або повної синхронії, або повної неузгодженості. Але чисельно виявили стійке співіснування: приблизно половина кільця злипалася у когерентний, синхронізований кластер, а інша половина залишалася неузгодженою й фактично хаотичною — нескінченно довго, з чіткою, стійкою межею між двома областями.

2. Нелокальний зв'язок: модель

Модель Курамото-Баттогтоха (KB) розміщує N ідентичних фазових осциляторів θᵢ на кільці, кожен зв'язаний з усіма іншими через ядро G, що спадає з відстанню:

Де:

• ω — власна частота, однакова для кожного осцилятора
• G(x) — ядро зв'язку, наприклад G(x) ∝ exp(−κ|x|), що спадає з відстанню вздовж кільця
• α — параметр фазового запізнення; ненульова α є необхідною умовою для появи химер взагалі

Критичним інгредієнтом є фазове запізнення α (іноді записується як π/2 − α). Без нього модель зводиться до звичайного притягуючого зв'язку, який завжди повністю синхронізується. При α, близькому до π/2, але меншому за нього, зв'язок стає майже (але не повністю) чисто дисипативним — це саме той режим, у якому з'являються химери.

3. Чому це дивно?

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

Це приклад спонтанного порушення симетрії в динамічній системі: керуючі рівняння мають більше симетрії, ніж їхні стійкі розв'язки. Розташування когерентного кластера саме по собі довільне — воно може опинитися будь-де на кільці залежно від початкових умов, але після формування розділ зберігається нескінченно довго.

4. Механізм: аргумент самоузгодженості середнього поля

У континуальній границі (N → ∞) Курамото та Баттогтох звели задачу до рівняння самоузгодженості для комплексного локального параметра порядку Z(x,t) = R(x,t)e^(iΨ(x,t)), де R вимірює локальну фазову когерентність у позиції x на кільці.

У когерентній області R ≈ 1 — осцилятори там обертаються зі спільною ефективною частотою, всі зчеплені між собою. У неузгодженій області R < 1, і локальні осцилятори дрейфують з частотами, розкиданими навколо середнього, ніколи не зчеплюючись. Рівняння самоузгодженості допускає стійкий розв'язок, у якому R(x) має саме такий двообласний профіль, із плавною, але швидкою перехідною зоною між ними — що відповідає чисельним спостереженням.

Важливо, що неузгоджені осцилятори не просто «неузгоджені» — вони відчувають ефективне середнє поле, створене когерентним кластером, і це середнє поле досить сильне, щоб не дати їм повністю розсинхронізуватися, але заслабке, щоб зчепити їх з когерентною групою. Вони перебувають у справді проміжному, частково зчепленому режимі.

5. Лабораторні підтвердження

Химерні стани залишалися суто теоретичними до 2012 року, коли дві незалежні групи отримали експериментальні докази:

• Хагерстром та ін. (2012) — кільце зв'язаних оптоелектронних (модульованих світлом) осциляторів, що демонструвало стійкі химери.
• Тінслі, Нкомо та Шовалтер (2012) — популяція хімічних осциляторів Білоусова-Жаботинського (див. також нашу статтю про реакцію-дифузію Білоусова-Жаботинського), зв'язаних через петлю зворотного зв'язку камера-світло, що показала співіснування синхронізованих і несинхронізованих підгруп.

Відтоді химери спостерігалися в механічних метрономах, хімічних осциляторах, лазерних масивах і симульованих нейронних мережах — підтверджуючи, що це явище не є артефактом однієї конкретної моделі, а стійкою рисою систем зв'язаних осциляторів загалом.

6. Значення: однопівкульний сон і не тільки

Химерні стани дають переконливу математичну метафору для однопівкульного повільнохвильового сну, що спостерігається у дельфінів, деяких птахів і тюленів: одна півкуля мозку показує синхронізовану повільнохвильову (сонну) активність, тоді як інша залишається активною та неузгодженою (притомною), дозволяючи тварині тримати одне око відкритим для хижаків, відпочиваючи іншою половиною мозку. Хоча біологічний механізм складніший за ідеалізоване кільце KB, якісна картина — ідентичні популяції нейронів, що розпадаються на когерентну й неузгоджену області — близько відповідає химерному феномену.

Химероподібну динаміку також пропонували як модель початку епілептичного нападу (локалізована гіперсинхронність на тлі нерегулярної активності) та подій розсинхронізації енергомереж, коли частина мережі зв'язаних генераторів випадає з такту решти.

7. Псевдокод

// N ідентичних фазових осциляторів на кільці, нелокальний зв'язок
const N = 200
const omega = 1.0       // ідентична власна частота
const K = 4            // сила зв'язку
const alpha = 1.45     // фазове запізнення, близько до PI/2
const kappa = 4        // швидкість спадання ядра

// Експоненційне ядро зв'язку, періодичне на кільці
function kernel(dx):
  return (kappa/2) * Math.exp(-kappa * Math.abs(dx))

theta = randomArray(N)  // майже ідентичні випадкові початкові фази

loop (кожен крок):
  dtheta = new Array(N)
  for i = 0 to N-1:
    sum = 0
    for j = 0 to N-1:
      dx = ringDistance(i, j, N)
      sum += kernel(dx) * Math.sin(theta[i] - theta[j] + alpha)
    dtheta[i] = omega - (K/N) * sum
  for i = 0 to N-1:
    theta[i] += dtheta[i] * dt

  renderPhases(theta)  // виникає когерентна дуга vs неузгоджена дуга