[論文レビュー] 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

ひとこと要約

本稿は、ランダムな自然周波数を持つ局所的結合による1次元環状のカイラウト・オシレーター系における同期現象を研究する。解析的および数値的に、複数の安定な同期解が共存可能であり、それぞれが異なる巻き数（winding number）を持つことを示し、臨界結合強度 $k_c$ の下界を導出する。1次元系では、位相ロックが位相整合性（phase coherence）を意味しないことにも注意を向ける。

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.

研究の動機と目的

• 有限サイズの1次元局所的結合オシレーター系が、無限大の1次元系では同期しないにもかかわらず、同期を示すかどうかを調査すること。
• このような系において、複数の安定な同期解が出現する条件を特定すること。
• 解の安定性と吸引域の大きさを決定づける巻き数 $m$ の役割を同定すること。
• 臨界結合強度 $k_c$ の解析的下界を導出し、その正確な計算のためのアルゴリズムを提示すること。
• 1次元カイラウトモデルにおける位相ロックと位相整合性の違いを明確にすること。

提案手法

• 局所的1次元カイラウトモデルを、リング構造上の最近接結合を用いて、オシレーター位相の1階常微分方程式系として定式化する。
• 定常的同期解（定常的位相差を有する）を記述する線形系 $\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]$ と定義される。
• ゲルシュゴリンの円定理を用いて安定性を解析し、$|\phi_i - \phi_{i-1}| < \pi/2$（安定）または $> \pi/2$（不安定）の基準に基づいて解を分類する。
• 数値的シミュレーションを用いて、巻き数 $m$ と固有値の挙動、吸引域体積の関係を相関づける。より負の固有値を有する解は、より大きな吸引域を有すると仮定する。

実験結果

リサーチクエスチョン

• RQ1 有限1次元局所的結合カイラウト・オシレーター系は、無限大1次元系では同期しないにもかかわらず、安定な同期を示すことができるか？
• RQ2 このような系に共存可能な異なる安定な同期解はいくつ存在するか？その数は何かによって決定されるか？
• RQ3 同期解の巻き数 $m$ とその安定性、吸引域の大きさの関係は何か？
• RQ4 臨界結合強度 $k_c$ は系のサイズや周波数分布にどのように依存するか？正確に計算可能か？
• RQ5 標準的なカイラウト順序パラメータがなぜ1次元系では周波数同期を検出できないのか？どのような代替指標が必要か？

主な発見

• 局所的結合を有する有限1次元系は、それぞれが異なる巻き数 $m$ を持つ複数の安定な同期解を支持可能である。
• 安定解の数は系のサイズ $N$ に比例し、$m \in [-\lfloor N/4\rfloor, \lfloor N/4\rfloor]$ により分類される。
• 臨界結合強度 $k_c$ には、提示されたアルゴリズムを用いて正確に計算可能な下界が存在し、周波数分布に依存する。
• $m$ の絶対値が大きい解は、摂動行列の固有値がより小さくなる（$\lambda_i \sim -k\cos(2\pi m / N)$ により）ため、劣位となる。
• 吸引域は安定性と相関する：より負の固有値（より安定）を持つ解はより大きな吸引域を有する。数値的証拠により、中程度の $k$ では $m=0$ が支配的であり、$k$ が高くなると $m=\pm1$ が現れる。
• 1次元系では、位相ロックが位相整合性を意味しない。周波数同期は高い順序パラメータがなくても発生可能であり、これは平均場カイラウトモデルとは根本的に異なる点である。