[论文解读] Synchronization in a stochastic Hebbian network of phase oscillators

Alexander Isakov, L. Mahadevan|arXiv (Cornell University)|Apr 8, 2014

Nonlinear Dynamics and Pattern Formation参考文献 2被引用 4

一句话总结

该论文提出了一种相位振子的随机赫布网络，其中耦合强度根据噪声驱动的动力学缓慢演化，即使在振子完全相同的情况下，也能实现同相与反相状态之间的转换。关键结果表明，同步在多种拓扑结构中均具有鲁棒性，相变由噪声强度和耦合演化速率决定，且通过有效序参量动力学的最小模型可解释观测到的翻转行为。

ABSTRACT

We investigate a generalized Kuramoto phase-oscillator model with Hebb-like couplings that evolve according to a stochastic differential equation on various topologies. Numerical simulations show that even with identical oscillators, there is a regime in the nearest-neighbor coupling topologies and a complex network topology where oscillators move between an in phase and anti-phase state. Phase diagrams show the transition probabilities as a function of the noise strength and rate of evolution of network coupling. A minimal theoretical model allows us to understand these transitions.

研究动机与目标

• 建立一个耦合强度基于相位差演化、受神经系统赫布可塑性启发的网络同步模型。
• 研究噪声和缓慢耦合动力学如何影响具有相同振子的网络中的同步转换。
• 理解在最近邻和复杂网络拓扑中，间歇性同相与反相状态的涌现机制。
• 构建一个最小理论模型，解释观测到的同步与反同步状态之间的切换行为。
• 探索同步在不同网络结构（包括随机无标度网络）中的鲁棒性。

提出的方法

• 通过引入时间尺度分离，对 Kuramoto 模型进行改进：快速相位动力学与缓慢、随机的耦合演化，后者由类似赫布的规则控制。
• 对相位动力学（公式 2）和耦合演化（公式 3）均使用随机微分方程，引入噪声项以增强现实性。
• 提出一种耦合演化规则：当振子同步时增强连接，反同步时减弱连接，以模拟赫布可塑性。
• 采用序参量 $O_1$、$O_2$ 和 $O_K$ 量化相位与耦合同步，支持相图分析。
• 在环形、二维晶格和优先连接网络上进行数值模拟，以研究切换行为与稳定性。
• 通过近似耦合动力学的总和，推导出 $O_1(t)$ 的简化平均场模型，将 $O_K$ 与 $O_1$ 联系起来，解释序参量的指数衰减。

实验结果

研究问题

• RQ1 噪声强度和耦合演化速率如何影响相同振子网络中同相与反相同步之间的转换？
• RQ2 网络拓扑结构（环形、二维晶格、无标度网络）在促进或抑制同步与反同步状态之间的间歇性切换中起什么作用？
• RQ3 能否通过一个最小理论模型解释观测到的切换动力学以及序参量 $O_1$ 的指数衰减？
• RQ4 同步在多大程度上对拓扑异质性具有鲁棒性，特别是在最小度和平均度不同的无标度网络中？
• RQ5 耦合符号 $\text{sign}(K_{ij})$ 对预测邻近振子之间实际相位关系的准确性如何？

主要发现

• 即使在完全相同的振子下，由于噪声和缓慢耦合演化的作用，该系统在环形和二维最近邻拓扑中仍表现出同相与反相状态之间的间歇性切换。
• 相图显示，切换频率强烈依赖于噪声强度 $\sigma$ 和耦合演化速率 $\epsilon$，高噪声会导致快速失同步。
• 耦合强度的序参量 $O_K$ 密切跟踪 $O_1$，表明耦合动力学是宏观相位同步的可靠预测指标。
• 在无标度网络中，系统对拓扑异质性表现出鲁棒性：高连接度节点因数量稀少，对全局行为影响甚微；而提高最小度会抑制切换行为。
• 当 $l=3$（最小度为 3）时，即使在 $\sigma=0.30$ 的情况下也未观察到一致的翻转，表明在高连通性约束下切换行为被完全抑制。
• 理论模型预测 $O_1(t)$ 指数衰减至稳态，衰减时间尺度由 $\epsilon$ 决定，与无翻转区域的数值结果一致。