Deviations from global coupling in adaptive oscillator networks: a mean-field theory for the variance of coupling weights

技術要約：適応型振動子ネットワークにおけるグローバル結合からの逸脱

問題提起
複雑な物理的および生物学的システムは、ノードのダイナミクスとエッジの重みが共進化する適応型ネットワークとしてモデル化されることが多い。しかし、このようなシステムの平均場簡約化は通常、平均的な結合強度のみを追跡するため、結合が実質的に均質であるレジームと、構造化された連結性が創発するレジームを区別できない。Ott-Antonsen（OA）多様体簡約（一様な結合重み $A_{ij} = 1$ を仮定した、グローバル結合された蔵本振子の研究など）のような既存の理論は、重みがノードのダイナミクスに基づいて進化する適応型ネットワークにおいて、高次の統計モーメント（具体的には結合重みの分散）が重要になり、システムがOA多様体から逸脱する可能性を見落としている。本研究が扱う中心的な問題は、適応型ネットワークがいつグローバル結合システムのように振る舞い、いつ複雑で構造化された結合パターンが創発するのかという条件を特定することであり、また、これら適応型重みの分散をどのように解析的に記述するかである。

手法
著者らは、ローレンツ分布に従う固有振動数 $\omega_i$ を持つ、非同一の蔵本位相振子のネットワークを研究している。ネットワークのダイナミクスは以下によって支配される：
$$\dot{\theta}_i = \omega_i + \frac{K}{N} \sum_{j=1}^N A_{ij}(t) \sin(\theta_j - \theta_i)$$
結合重み $A_{ij}(t)$ は、以下の適応則に従って進化する：
$$\dot{A}_{ij} = \mu G(\theta_j - \theta_i) + \gamma(1 - A_{ij})$$
ここで、$G$ は位相差の任意の関数であり、$\mu$ は適応率、$\gamma$ は減衰項である。本研究では、$G$ についての2つの特定のケース、すなわち対称型（$G_s(\phi) = \cos \phi$）と反対称型（$G_a(\phi) = \sin \phi$）を検討する。

手法は以下の2段階で進行する：

1. 一次平均場解析： 著者らは、平均結合重み $\bar{A}$ と位相コヒーレンス $R$ に関する方程式を導出する。これらの方程式の不動点と分岐構造を解析し、微視的なシステムによる数値シミュレーションと比較する。
2. 二次モーメント閉鎖： グローバル結合からの逸脱に対処するため、著者らは結合重みの分散 $V_A = \text{Var}(A_{ij})$ に関する閉形式の平均場方程式を導出する。これには以下が含まれる：
   • 共分散 $C_A = \text{Cov}(A_{ij}, G_{ij})$ の定義。
   • 適応駆動 $G_{ij}$ を、静的な時間平均成分（$\Gamma_{ij}$）と変動成分（$g_{ij}$）への分解。
   • 高次モーメント（例えば $S = \langle |c|^2 \rangle$）を、位相コヒーレンス $R$ と位相ロックした振子の割合を用いて表現するためのOA多様体仮定の利用。
   • 平均重み $\bar{A}$、位相コヒーレンス $R$、重みの分散 $V_A$、および静的成分と変動成分の共分散に関する連立微分方程式の導出。

主な貢献

1. 分散方程式の導出： 本論文は、一様な初期条件から出発して、適応型結合重みの分散に対する平均場方程式を与える、新しい二次モーメント閉鎖を提示している。これにより、結合の不均質性が支配的になる時期を解析的に予測することが可能となる。
2. 非単調依存性： 著者らは、振子の不均質性（$\Delta$）と相対的な重みの分散（$V_A/\bar{A}^2$）との間の、位相コヒーレンスによって媒介された非線形で非単調な関係を導出した。
3. 対称性に依存するダイナミクス： 本研究は、対称型および反対称型の適応則の影響を解析的に区別している：
   • 対称的適応： 位相コヒーレンスが平均結合重みを駆動するフィードバックループを生じさせ、非適応型または反対称型システムには存在しない双安定レジームを作り出す。これは、コヒーレントな振子の強力に結合されたコアの創発を促進する。
   • 反対称的適応： 平均結合重みは一定（$\bar{A} \to 1$）に保たれるが、分散は著しく増大する。この規則は、コヒーレントなコア内での反対称結合をもたらし、同期状態を不安定化させ、カオスのような複雑なダイナミクスを誘発する可能性がある。

結果

• 分岐構造： 対称的適応の場合、位相コヒーレンスと平均結合重みの相互作用が分岐図を変化させ、フォールド分岐と、初期条件に応じてシステムが同期状態または非同期状態のいずれかに存在できる双安定領域を導入する。これは、ヒステリシスを示す数値シミュレーションによって確認されている。
• 分散と不均質性： 導出された方程式は、相対的な重みの分散が特定の振子の不均質性のレベルでピークに達することを予測する。対称的なケースでは、グローバル結合からの逸脱は、周波数分布の裾の部分にある振子（脱同期する振子）に由来する。反対称的なケースでは、偏差は分布の中心部（非ゼロの位相ラグを伴って同期する振子）に由来し、OA多様体からのより強い逸脱を引き起こす。
• 平均場理論の精度： 平均場方程式は、対称的なケースにおいて、不均質性、コヒーレンス、および重みの分散の関係を正確に予測する。しかし、反対称的なケースで高い適応率（$\mu > \gamma$）を持つ場合、平均場の分散の予測は、同期が発生する低不均質性のレジームにおいて微視的なシミュレーションから乖離する。これは、高次の効果または非OA多様体ダイナミクスが重要になっていることを示唆している。

意義
著者らは、本研究が適応型ネットワークの統計物理学を理解するための決定的な解析ツールを提供すると主張している。適応型ネットワークがグローバル結合システムとして振る舞うレジーム（平均重みが支配的なモーメントである場合）と、構造化された結合パターンが形成されるレジーム（分散が支配的な場合）を明確にすることで、導出された方程式は複雑な時空間パターンの創発を予測する枠組みを提供する。本研究の結果は、ノードの不均質性と適応則の対称性が、ネットワークの挙動を制御する重要なパラメータであることを示唆しており、電力網から神経集団に至るまでのシステムの安定性を理解することへの示唆を含んでいる。本論文は、標準的なOA簡約の限界を超え、高次の統計モーメントを考慮した適応型システムのための低次元方程式を導出するためのステップとして位置づけられている。