[论文解读] Synchronization of Coupled Oscillators in a Local One-Dimensional Kuramoto Model

Jeremi K. Ochab, Paweł Góra|arXiv (Cornell University)|Aug 31, 2009

Nonlinear Dynamics and Pattern Formation参考文献 9被引用 18

一句话总结

本文研究了具有随机自然频率的有限一维环形局域耦合Kuramoto振子系统中的同步现象。通过解析与数值方法证明，多个稳定的同步解可共存，每种解由不同的绕数（winding number）表征，并推导出临界耦合强度 $k_c$ 的下界；在一维系统中，相位锁定并不等价于相位一致。

ABSTRACT

A modified Kuramoto model of synchronization in a finite discrete system of locally coupled oscillators is studied. The model consists of N oscillators with random natural frequencies arranged on a ring. It is shown analytically and numerically that finite-size systems may have many different synchronized stable solutions which are characterised by different values of the winding number. The lower bound for the critical coupling is given, as well as an algorithm for its exact calculation. It is shown that in general phase-locking does not lead to phase coherence in 1D.

研究动机与目标

• 研究有限尺寸的一维局域耦合振子系统是否能在无限一维系统不发生同步的情况下实现同步。
• 确定此类系统中多个稳定同步解出现的条件。
• 阐明绕数 $m$ 在决定解的稳定性与吸引盆大小中的作用。
• 推导临界耦合强度 $k_c$ 的解析下界，并提供其精确计算的算法。
• 澄清一维Kuramoto模型中相位锁定与相位一致的区别。

提出的方法

• 使用一阶常微分方程在环形结构上建立局域一维Kuramoto模型，实现最近邻耦合的振子相位描述。
• 推导线性系统 $\mathbf{A}\mathbf{s} = \frac{2}{k}\boldsymbol{\Delta}$，用于描述具有恒定相位差的稳态同步解。
• 引入自由参数 $p \in [-1,1]$ 参数化解，导致 $2^{N-1}$ 种潜在的相位差配置。
• 利用绕数 $m$ 通过方程 $\sum_{i=0}^{N-1} \arcsin\left(p + \frac{2}{k}\sum_{j=1}^{i} \Delta_j\right) = 2m\pi$ 对解进行分类，其中 $m \in [-\lfloor N/4\rfloor, \lfloor N/4\rfloor]$。
• 应用Gershgorin圆定理分析稳定性，根据 $|\phi_i - \phi_{i-1}| < \pi/2$（稳定）或 $> \pi/2$（不稳定）对解进行分类。
• 通过数值模拟关联绕数 $m$ 与特征值行为及吸引盆体积，假设更负的特征值对应更大的吸引盆。

实验结果

研究问题

• RQ1 有限一维局域耦合Kuramoto振子系统是否能在无限一维极限下不发生同步的情况下实现稳定同步？
• RQ2 此类系统中可共存的稳定同步解有多少种，其数量由什么决定？
• RQ3 同步解的绕数 $m$ 与其稳定性及吸引盆大小之间存在何种关系？
• RQ4 临界耦合强度 $k_c$ 如何依赖于系统尺寸与频率分布，是否可精确计算？
• RQ5 为何标准Kuramoto序参数量无法检测一维系统中的频率同步，需要何种替代度量？

主要发现

• 具有局域耦合的有限一维系统可支持多个稳定的同步解，每种解由不同的绕数 $m$ 表征。
• 稳定解的数量与系统尺寸 $N$ 成正比，解按 $m \in [-\lfloor N/4\rfloor, \lfloor N/4\rfloor]$ 分类。
• 临界耦合强度 $k_c$ 存在可通过所提算法精确计算的下界，且依赖于频率分布。
• 绕数绝对值较大的解因扰动矩阵特征值较不为负而被边缘化，如 $\lambda_i \sim -k\cos(2\pi m / N)$ 所示。
• 吸引盆与稳定性相关：特征值更负（更稳定）的解具有更大的吸引盆，数值证据表明在中等 $k$ 时 $m=0$ 占主导，而 $m=\pm1$ 在更高 $k$ 时出现。
• 在一维系统中，相位锁定不意味着相位一致；频率同步可在序参数不高时发生，表明其与平均场Kuramoto模型存在根本差异。