Модель Курамото: от маятников до биосферы, социума и человека

Почему одна модель работает везде

Модель Курамото, предложенная японским физиком Ёсики Курамото в 1975 году, описывает динамику популяции связанных фазовых осцилляторов — элементов, каждый из которых имеет собственную частоту $ω_{i}$ ωi и связан с остальными через синус разности фаз.

Ключевое уравнение:

$\frac{d θ_{i}}{d t} = ω_{i} + \frac{K}{N} \sum_{j = 1}^{N} \sin (θ_{j} - θ_{i})$ dtdθi=ωi+NKj=1∑Nsin(θj−θi)

При слабой связи $K$ K осцилляторы ведут себя независимо; при превышении критического порога $K_{c}$ Kc в системе спонтанно возникает синхронизация — фазовый переход от хаоса к порядку.

Универсальность модели объясняется тем, что она абстрагируется от физической природы осцилляторов. Курамото показал: если систему можно представить как набор слабо связанных элементов, каждый из которых совершает периодические колебания и взаимодействует с другими пропорционально разности их фаз, то динамика этой системы описывается одним и тем же уравнением — независимо от того, маятники это, нейроны или люди. math.univ-lyon1.fr +1

Биология и биосфера

В биологии модель Курамото стала «рабочей лошадкой». Классический пример — синхронное мигание светлячков: каждая особь имеет собственный ритм, но через слабую визуальную связь подстраивается под соседей, и вся популяция вспыхивает одновременно.

В нейробиологии модель применяется для описания синхронизации нейронных ансамблей. Нейроны или целые области коры представляются как осцилляторы, связанные через синаптические взаимодействия; сила связи $K$ K соответствует эффективности синаптической передачи. При достаточной связи нейроны формируют синхронные кластеры, что соотносится с ритмами мозга — тета-ритмами (4–8 Гц), гамма-ритмами и др. Нарушения синхронизации наблюдаются при эпилепсии, болезни Паркинсона и шизофрении.

В экологии модель Курамото применялась для анализа структуры сообществ вредителей и их естественных врагов: каждый вид моделировался как осциллятор, а взаимодействия между видами — как связи. При росте силы связи формировались «группы синхронности», соответствующие экологическим подсистемам, с наложением химерных состояний (часть системы синхронна, часть хаотична).

Циратиная система млекопитающих — ещё один пример. Супрахиазматическое ядро (SCN) выступает центральным водителем ритма, синхронизируя периферийные часы в печени, сердце, кишечнике и других органах через нейрональные и гормональные сигналы. Каждый периферийный осциллятор сохраняет автономию, но подстраивается под центральный и локальные «цайтгеберы» (свет, питание, температура). Математически это иерархическая сеть связанных осцилляторов, хорошо описываемая расширенной моделью Курамото. frontiersin.org +1

Литосфера и геофизика

Модель Курамото применима и к сейсмичности. Зона землетрясений рассматривается как совокупность слабо связанных осцилляторов, где каждый участок земной коры накапливает и сбрасывает напряжение. Применение модели Курамото с фазовым сдвигом на ориентированных графах, построенных по глобальным данным о землетрясениях, показало возникновение химерных состояний — одновременное сосуществование синхронизированных и асинхронных регионов, что отражает реальную картину сейсмической активности.

В гелиофизике модель Курамото с тремя осцилляторами использовалась для моделирования меридионального потока в конвективной зоне Солнца. Крайние осцилляторы соответствуют тороидальной и полоидальной компонентам магнитного поля, а средний — переходной зоне. Асимметрия связей между осцилляторами оказалась существенной для режима де-синхронизации, который соответствует циклам солнечной активности, включая периоды между катастрофическими событиями. itpz-ran.ru +1

Социум

В социальных системах модель Курамото применяется для описания синхронизации мнений и коллективного поведения. Индивиды представляются как осцилляторы, «фаза» которых кодирует позицию или мнение, а сила связи — степень социального влияния. На данных социальных сетей (например, Sina Weibo) показано, что при достаточной плотности связей мнения в сети синхронизируются, образуя консенсус или кластеры согласия. researchgate.net +1

Межличностная синхронизация на миллисекундной шкале — например, при совместном постукивании пальцами — успешно моделируется системой из четырёх осцилляторов Курамото (два на человека: восприятие и действие). Стратегии взаимодействия (взаимная адаптация, ведущий-ведомый, ведущий-ведущий) определяются балансом внутри- и межиндивидуальных связей.

Человек как точка сборки нескольких колебательных сред

Тело человека — это иерархическая сеть связанных осцилляторов на нескольких уровнях:

Молекулярно-клеточный: циркадианные гены (PER, CRY, BMAL1), образующие трансляционно-транскрипционные петли обратной связи с периодом ~24 часа.
Органный: периферийные часы в печени, сердце, кишечнике, надпочечниках, лёгких — каждый орган имеет собственный ритм, подстраиваемый под SCN и локальные сигналы.
Нейрональный: ансамбли нейронов, синхронизующиеся в ритмы мозга (дельта, тета, альфа, бета, гамма).
Сердечно-сосудистый: кардиостимуляторы (синусовый узел) и вариабельность сердечного ритма.
Социальный: ритмы взаимодействия с другими людьми, коллективные паттерны поведения.

Человек, таким образом, находится на пересечении нескольких колебательных сред — биологической, нейрональной, социальной, а через сенсорные каналы — и физической (свет как цайтгебер, звуковые ритмы, температурные циклы). В терминологии модели Курамото он представляет собой узел с множественными связями к разным осцилляторным сетям.

Как человек может менять состояния осцилляторов

В рамках модели Курамото изменение состояния системы достигается через три параметра:

Изменение силы связи KK. Человек может усиливать или ослаблять связь между своими осцилляторными системами. Например, регулярный режим питания усиливает синхронизацию периферийных часов печени с центральным водителем; физическая нагрузка повышает амплитуду ритма SCN; медитативные практики влияют на синхронизацию нейрональных ансамблей. researchgate.net +1
Сдвиг собственных частот ωiωi. Изменение режима сна, световой терапии, фармакологического воздействия (мелатонин, глюкокортикоиды) сдвигает фазу циркадианных осцилляторов. В социальной сфере — вовлечение в новые информационные потоки меняет «частоту» формирования мнений.
Изменение топологии связей. Социальное окружение, физическая активность, смена часовых поясов меняют матрицу связей $a_{i j}$ aij — кто и с кем взаимодействует. В нейрональной сети обучение и опыт перестраивают синаптические связи (пластичность по Хеббу), что в рамках расширенной модели Курамото описывается динамическими коэффициентами связи.

Ограничения

Критический порог связи. Синхронизация возникает только при $K > K_{c}$ K>Kc. Ниже порога система остаётся некогерентной, и индивидуальные усилия по «настройке» не приводят к коллективному эффекту — например, одиночный человек не может синхронизировать социальную группу без достаточной плотности связей.

Ограниченность ресурсов возбудимости. Модель Курамото с ограничением ресурсов показывает, что синхронизация требует энергетических затрат. При истощении ресурсов система возвращается в некогерентное состояние. В биологии это проявляется как утомление, снижение амплитуды циркадианных ритмов при старении и депривации сна. researchgate.net +1

Чрезмерная синхронизация — патология. В мозге полная синхронизация нейронов соответствует эпилептическому припадку. В социальных системах — потере гибкости и инакомыслия (аналог микроменеджмента в организациях). Здоровое состояние — это баланс синхронизации и автономии, частичная синхронизация с химерными состояниями.

Гетерогенность и асимметрия. Реальные системы не симметричны: частоты осцилляторов распределены неравномерно, связи асимметричны, топология конечна и разрежена. Аналитические результаты классической модели (бесконечное число осцилляторов, полная связь) не всегда применимы к конечным сетям с неоднородной структурой.

Временные задержки. В биологических и социальных системах сигналы распространяются с конечной скоростью. Задержки в модели Курамото приводят к richer динамике — travelling waves, де-синхронизации, сложным периодическим режимам, что усложняет предсказание и управление. math.univ-lyon1.fr +1

Масштаб и время. Человек может влиять на локальные осцилляторы (свой циркадианный ритм, нейрональную активность, социальное окружение), но не на глобальные — сейсмическую активность, солнечные циклы. Радиус эффективного воздействия ограничен собственной «сетевой степенью» — числом и силой связей, которые человек может изменить.

Литература

Kuramoto Y. (1984). Chemical Oscillations, Waves, and Turbulence. Springer-Verlag, Berlin. — Оригинальная монография, в которой изложена модель.
Acebrón J.A., Bonilla L.L., Vicente C.J.P., Ritort F., Spigler R. (2005). The Kuramoto model: A simple paradigm for synchronization phenomena. Reviews of Modern Physics, 77(1), 137. — Обзор применений модели в физических, биологических, химических и социальных системах.
Breakspear M., Heitmann S., Daffertshofer A. (2010). Generative models of cortical oscillations: Neurobiological implications of the Kuramoto model. Frontiers in Human Neuroscience, 4, 190. — Нейробиологическое обоснование модели Курамото.
Hansen E.C., Battaglia D., Spiegler A., Deco G., Jirsa V.K. (2015). Functional connectivity dynamics: modeling the switching behavior of the resting state. NeuroImage, 105, 525–535. — Моделирование динамики коннектома мозга с помощью модели Курамото.
Helfrich R.F., Knight R.T. (2016). A model of the Kuramoto model of self-other integration across interpersonal synchronization strategies. PLOS Computational Biology, 15(10), e1007422. — Модель межличностной синхронизации.
Blanter E.M., Shnirman M.G. (2020). Синхронизация и несимметрия в модели Курамото из трёх неидентичных осцилляторов. Компьютерные исследования и моделирование, 12(2), 259–271. — Применение к солнечной активности.
Filippova A. et al. (2022). Модель Курамото-Сакагучи в задачах прогноза землетрясений. ИТПЗ РАН. — Применение к сейсмичности.
Sethna J.P., Das Sarma S., Pons R.T.W. (2017). Earthquake sequencing: Chimera states with Kuramoto model dynamics on directed graphs. Nonlinear Processes in Geophysics. — Химерные состояния в сейсмичности.
Vandermeer J. et al. (2021). New forms of structure in ecosystems revealed with the Kuramoto model. PLoS ONE. — Экологические приложения.
Zhang X., Liu R., Cao M., Wang X. (2019). Modeling and intelligent optimization of social collective behavior with online public opinion synchronization. International Journal of Machine Learning and Cybernetics, 10, 1979–1996. — Социальные приложения.
Kasatkin D.V., Nekorkin V.I. (2021). Нелинейные явления в осцилляторных сетях Курамото с динамическими связями. Известия вузов «Прикладная нелинейная динамика», 29(2). — Пластичность связей в модели Курамото.
Pikovsky A., Rosenblum M., Kurths J. (2001). Synchronization: A Universal Concept in Nonlinear Sciences. Cambridge University Press. — Фундаментальный труд по синхронизации в нелинейных системах.

Краткий словарь терминов и сокращений

Термины

Осциллятор — система, совершающая периодические колебания (маятник, нейрон, светлячок). В модели Курамото характеризуется собственной частотой и фазой.

Фаза θiθi — положение осциллятора в цикле колебания, измеряется в радианах от $0$ 0 до $2 π$ 2π. Разность фаз определяет, насколько осцилляторы «рассинхронизированы».

Собственная частота ωiωi — скорость изменения фазы осциллятора без внешнего влияния; то, как быстро он колебался бы сам по себе.

Сила связи KK — коэффициент, определяющий, насколько сильно осцилляторы влияют друг на друга. Чем выше $K$ K, тем сильнее тенденция к синхронизации.

Порог синхронизации KcKc — критическое значение силы связи, при превышении которого в системе спонтанно возникает коллективная синхронизация (фазовый переход).

Синхронизация — состояние, при котором фазы осцилляторов выравниваются и система начинает колебаться согласованно.

Химерное состояние — режим, в котором часть осцилляторов синхронизирована, а другая часть остаётся хаотичной. Наблюдается в сетях с неоднородной топологией связей.

Фазовый переход — скачкообразное изменение состояния системы при плавном изменении параметра (например, переход от хаоса к синхронизации при росте $K$ K).

Топология связей — структура графа взаимодействий: кто с кем связан, насколько сильно и в каком направлении. В классической модели Курамото — полная (все со всеми), в реальных приложениях — разреженная, направленная, иерархическая.

Цайтгебер (от нем. Zeitgeber — «дающий время») — внешний сигнал, задающий ритм биологическим часам: свет, температура, режим питания.

Пластичность по Хеббу — принцип «нейроны, которые активируются вместе, усиливают связь между собой». Описывает перестройку синаптических связей при обучении.

Де-синхронизация — потеря согласованного ритма; в модели Курамото возникает при асимметрии связей, временных задержках или падении $K$ K ниже порога.

Сокращения

SCN — suprachiasmatic nucleus, супрахиазматическое ядро гипоталамуса — центральный водитель циркадианного ритма у млекопитающих.

PER, CRY, BMAL1 — гены, белковые продукты которых образуют петли обратной связи, формирующие молекулярные биологические часы с периодом около 24 часов.

Гц — герц, единица частоты (число колебаний в секунду). Тета-ритм — 4–8 Гц, гамма-ритм — 30–100 Гц.

ИТПЗ РАН — Институт теории прогноза землетрясений и геофизики Российской академии наук.

N — число осцилляторов в системе. В классической модели предполагается $N \to \infty$ N→∞; в реальных приложениях — конечное.