On the efficiency of pairwise Hamiltonian control to desynchronize the higher-order Kuramoto model

技术摘要：关于利用成对哈密顿控制使高阶库拉莫模型去同步化的效率研究

问题陈述
耦合振荡器系统中的同步是一种普遍现象，具有双重影响：在同步电网等背景下是有益的，但在其他背景下（如癫痫发作）则是有害的。虽然传统的控制策略通常依赖于成对相互作用，但近期的证据表明，许多自然和工程系统表现出高阶（多体）相互作用。这些相互作用从根本上改变了系统动力学，特别是通过增加同步态的线性稳定性，同时缩小其吸引盆（basin of attraction）——这种现象被描述为“更深但更小”的吸引盆。

本研究解决的核心问题是：在应用最初为成对网络开发的控制方案来使此类系统去同步化时，这些相互竞争的效应（即增加的线性稳定性和缩小的吸引盆）如何影响成对控制策略的效率。具体而言，作者研究了当将该控制方案应用于高阶库拉мо模型时，是高阶相互作用引起的线性稳定性增加（阻碍去同步化）还是吸引盆的缩小（促进去同步化）占据了主导地位。

方法论
研究采用了一个包含 $N$ 个相位振荡器的更高阶库拉莫模型，其中包含成对（$K_1$）和三体（$K_2$）相互作用。系统的控制方程为：
$$\dot{\theta}_i = \omega_i + \frac{K_1}{\langle k^{(1)} \rangle} \sum_{j} A^{(1)}_{ij} \sin(\theta_j - \theta_i) + \frac{K_2}{\langle k^{(2)} \rangle} \sum_{j,k} A^{(2)}_{ijk} \sin(\theta_j + \theta_k - 2\theta_i) + \mu \delta_i p_i$$
其中 $\mu$ 是控制强度，$\delta_i$ 表示节点 $i$ 是否受控。

控制机制是一种极简侵入式成对哈密顿控制，改编自 Asllani 等人 [35]。控制项 $p_i$ 定义为：
$$p_i = -\frac{1}{2K_1^2} \frac{M^2}{\langle k^{(1)} \rangle^2} R_M \tilde{R}_{M,i} \cos(\Psi_M - \tilde{\Psi}_{M,i})$$
该控制仅作用于 $M$ 个被钉扎（pinned）节点的子集。它测量受控节点的同步水平，并施加一个与同步幅度成正比的力，当系统完全去同步时，该力趋于零。这种方法避免了推导高阶控制项的计算复杂性。

作者使用两种序参数来评估控制效率：

1. 相位同步 ($\hat{R}$)：时间平均库拉莫序参数。
2. 频率同步 ($\hat{R}_{\dot{\theta}}$)：相位锁定（频率同步）的度量。

数值实验在两种拓扑结构上进行：

• 超环（Hyperrings）：支持扭转态（twisted states）的规则结构。
• 随机 Erdős-Rényi 超图：缺乏旋转不变性的不规则结构。

研究系统地改变了高阶耦合强度 ($K_2$)、控制强度 ($\mu$)、受控节点数 ($M$) 以及系统状态与同步吸引盆中心的初始距离 ($\epsilon$)。

主要结果

1. 同步初始条件：
   当系统从接近完全同步的状态（吸引盆中心）开始时，增加高阶耦合强度 $K_2$ 会使去同步化变得更加困难。所需的临界控制强度 ($\mu_c$) 随 $K_2$ 单调增加。这表明在此状态下，高阶相互作用引起的线性稳定性增加超过了吸引盆缩小的影响。成对控制在 $\mu$ 足够大时仍然有效，但控制成本随 $K_2$ 而上升。

2. 扰动初始条件：
   当系统从远离同步状态的位置（较大的 $\epsilon$）开始时，出现了非单调行为：

   • 中等 $K_2$ 时：去同步化最困难；控制效率最低。
   • 较大 $K_2$ 时：去同步化变得更容易；控制效率提高。
     这种行为表明，对于较大的 $K_2$，吸引盆的缩小成为了主导因素，使得控制能够更容易地将轨迹推离吸引盆，尽管此时线性稳定性更高。

3. 部分控制：
   研究证实，即使仅控制一部分节点 ($M < N$)，也可以实现去同步化。所需的最小节点数 ($M_c$) 在同步初始条件下通常随 $K_2$ 增加，但在扰动初始条件下表现出非单调行为，这与观察到的控制强度的趋势相一致。

4. 频率 vs. 相位去同步：
   该控制能有效抑制相位同步。然而，频率同步（例如在超环中的扭转态或簇态）更为稳健，特别是在中等 $K_2$ 的情况下，这是由于结构对称性以及控制并未直接针对频率锁定。

意义与主张
本文声称提供了一个关于高阶相互作用与控制效率之间相互作用的清晰图景。其主要贡献在于证明了高阶相互作用的“更深但更小”效应会导致取决于初始状态的不同控制机制：

• 在同步状态附近，高阶相互作用由于增强了稳定性而成为控制的障碍。
• 在远离同步状态时，足够强的高阶相互作用可以通过缩小吸引盆来实际上促进控制。

作者断言，一种简单的、极简侵入式的成对控制策略足以使高阶库拉莫系统去同步化，前提是控制强度或受控节点数量足够。这一发现具有重要意义，因为它表明，对于旨在实现去同步化的过程，并不一定需要复杂的、计算量巨大的高阶控制项（如文献 [100] 中所提出的），从而提供了一条更具实际操作性的路径。研究结果强调了在设计具有高阶相互作用的复杂网络控制策略时，考虑线性稳定性与吸引盆几何形状的重要性。