🌀 nonlinear sciences

Discrete-time Kuramoto model with phase lag: Linear stability analysis and onset of synchronization

本文通过推导非相干态的精确线性稳定性条件、解析确定洛伦兹频率分布下的同步阈值，并揭示了使其区别于连续时间对应模型的周期性与混沌态等独特非线性现象，对带有相位滞后的离散时间库拉莫托模型进行了研究。

原作者： Prashant M. Gade, Shamik Gupta

发布于 2026-09-22

📖 1 分钟阅读☕ 轻松阅读

原作者： Prashant M. Gade, Shamik Gupta

原始论文采用 CC BY 4.0 许可（http://creativecommons.org/licenses/by/4.0/）。 ✨ 这是对下方论文的AI生成解释。它不是由作者撰写或认可的。如需技术准确性，请参阅原始论文。阅读完整免责声明

在自然界中，从萤火虫的闪烁到心细胞的跳动，许多独立的个体往往会步调一致。这种被称为“同步”（synchronization）的现象，发生在许多具有各自自然节奏的独立振子在邻近个体的驱动下开始协同运动时。长期以来，科学家一直依赖于一个被称为库拉莫托模型（Kuramoto model）的数学框架来理解这种集体秩序是如何从混沌中涌现的。该模型将系统视为一个连续的时间流，其中振子的速度会平滑且逐渐地进行调整。它在解释大量耦合单元（如神经元或电网发电机）如何自发地锁定在单一的统一频率方面取得了巨大的成功。然而，许多现实世界的系统并不以平滑、不间断的流形式运行。相反，它们以离散的步骤运行，就像数字时钟向前跳动或计算机按离散间隔处理数据一样。这提出了一个根本性的问题：当时间本身被分解为独立的块而非连续流动时，同步的物理机制是否会发生变化？

研究人员 Prashant M. Gade 和 Shamik Gupta 致力于通过构建一个在离散时间步长下运行的著名模型来回答这个问题。他们没有观察振子如何平滑地漂移，而是设想了一个系统，其中每个单元的状态都是根据在每个时刻应用特定规则而同时更新的。他们在这一设置中引入了一个“相位滞后”（phase lag），它就像是一个延迟或一种在振子响应群体时的犹豫。在连续世界中，如果这种延迟变得过大，群体就会瓦解，同步也将变得不可能。该团队想要观察，当系统以跳跃而非流动的形式移动时，这一规则是否依然成立。通过分析这种逐步演化的数学过程，他们推导出了关于寻找特定相位的振子概率随时间变化的精确描述。这使他们能够准确计算出群体何时会失去其混沌的独立性并开始协同运动。

结果揭示了平滑世界与阶梯世界之间令人惊讶的差异。在连续模型中，一种特定的延迟（即响应正好落后半个周期）会完全阻止群体实现同步，无论成员之间的连接有多强。然而，在 Gade 和 Gupta 研究的离散时间版本中，即使在这些困难条件下，群体仍能实现同步。他们发现，同步的阈值（即群体锁定在一起的点）取决于延迟的大小，其方式与连续情况有着本质的不同。他们的计算表明，即使在振子之间存在排斥作用（即相互推开而非相互吸引）的情况下，离散系统仍然可以找到组织起来的方法。这表明，维持群体凝聚的机制足够强大，足以在那些会在平滑连续流中摧毁同步的条件下生存下来。

除了群体首次实现同步的时刻之外，离散系统的行为表现出了连续模型从未展现出的特性。一旦振子之间的耦合变得足够强，系统并不会仅仅稳定在一个稳态锁定状态。相反，它会开始表现出典型的非线性映射的复杂行为。研究人员观察到，群体的集体节奏可能会在两个截然不同的状态之间摆动，然后是三个，并最终变得混沌。在模拟中，他们观察到同步度会以一种规律的模式上升和下降，在不同的有序水平之间循环，最后破碎成不可预测的波动。这些发现表明，离定的时间版本模型支持着丰富多样的集体状态，包括周期性节奏和混沌，而这些在标准的连续时间描述中是被禁止的。

这项研究还挑战了一种广泛使用的数学简化方法，即 Ott-Antonsen 假设（Ott-Antonsen ansatz）。这种方法允许科学家将数百万个振子的复杂行为简化为一组更小、更易处理的方程，但它依赖于对系统如何演化的特定假设。研究人员发现，虽然这种简化方法在连续模型中表现良好，但在离散版本中却失效了。使这种简化有效的数学结构在时间以步长移动时并未得到保持。这意味着，适用于平滑流的优雅的低维描述不能直接应用于以离散跳跃演化的系统。团队通过运行包含数万个振子的规模庞大的计算机模拟，证实了他们的理论预测。这些模拟不仅符合他们计算出的同步阈值，还清晰地展示了其方程所预言的复杂多步节奏和混沌行为。

最终，这项工作证明了建模方式——无论是作为连续流还是作为一系列离散步骤——从根本上改变了同步的物理学。离散时间模型不仅仅是连续模型的数值近似，它是一个具有自身规则和可能性的截然不同的动力学范畴。系统在排斥力作用下仍能实现同步，并能表现出复杂的周期性和混沌状态，这表明离散时间动力学提供了更丰富的集体行为景观。通过确立同步的精确条件并揭示标准简化假设的失效，研究人员为理解如何在以步长运行的系统中（从数字网络到以离散间隔跳动的生物钟）产生秩序奠定了新的基础。

技术摘要：带有相位滞后的离散时间 Kuramoto 模型

问题陈述
本文研究了离散时间形式的 Kuramoto 模型，这是理解群体自发同步现象的一个范式框架。虽然连续时间版本已被广泛研究（特别是通过允许精确低维约化的 Ott–Antonsen (OA) 假设），但离散时间形式在解析研究方面仍较少被探索。作者探讨了一个基本问题：时间离散化本身是否会改变同步阈值和稳定性条件。具体而言，他们研究了一个在非线性映射下演化的全局耦合相位振子系统，其中包含一个相位滞后参数 ( $\beta$ )，旨在推导同步起始点的精确解析结果，并表征该阈值之后的非线性动力学。

研究方法
本研究聚焦于一个由具有洛伦兹分布自然频率 $\omega_i$ 的 $N$ 个全局耦合振子组成的系统。其动力学由以下非线性映射定义：
$\theta_i^{(n+1)} = \theta_i^{(n)} + \omega_i + K r_n \sin(\psi_n - \theta_i^{(n)} - \beta) \pmod{2\pi}$
其中 $K$ 是耦合常数， $\beta$ 是相位滞后， $(r_n, \psi_n)$ 是复序参数的幅度和相位。

在连续极限（ $N \to \infty$ ）下，作者使用单振子概率密度 $f(\theta, \omega, n)$ 来描述该系统。该密度的演化受与该非线性映射相关的 Frobenius-Perron (FP) 方程支配。研究方法如下：

线性稳定性分析： 作者对非相干态（ $r=0$ ）进行线性稳定性分析。他们将密度展开为傅里叶级数，并推导出第一傅里叶模 $\hat{f}_1(\omega, n)$ 的演化方程。
解析推导： 假设第一傅里叶模在下半复 $\omega$ 平面内是解析的，他们利用柯西留数定理处理序参数的表达式。这使得他们能够封闭方程组，并推导出序参数 $Z_n$ 的线性映射。
阈值计算： 通过寻找线性映射的不稳定条件（即 $|\lambda|=1$ ）来确定同步阈值。
数值验证： 通过对有限- $N$ 系统进行直接数值模拟，以验证解析阈值并探索超越该阈值后的非线性动力学（周期和混沌状态）。
OA 假设验证： 作者通过比较第二傅里叶模与第一模平方的关系，在离散时间机制下数值测试了 Ott–Antonsen 假设的有效性。

主要贡献与结果

精确同步阈值： 作者针对洛伦兹频率分布推导出了同步阈值 $K_c$ 的精确解析表达式：
$K_c = 2 \left( \sqrt{e^{2\gamma} - \sin^2 \beta} - \cos \beta \right)$
其中 $\gamma$ 是频率分布的半高全宽。这一结果与连续时间阈值（ $K_c = 2\Gamma/\cos\beta$ ）有着本质的区别。
排斥性耦合同步： 一个显著的发现是，即使在纯排斥耦合（ $\beta = \pi$ ）的情况下，同步依然存在。在连续时间 Kuramoto 模型中，当 $\beta \ge \pi/2$ 时是不可能实现同步的。然而，离散时间映射支持 $\beta = \pi$ 时的同步态，尽管这需要更大的耦合强度，但并未排除同步的可能性。
非线性现象： 在同步阈值之后，离散时间动力学表现出标准连续时间模型中所没有的丰富非线性行为。数值模拟显示：
- 周期轨道： 宏观序参数中出现了二周期和三周期轨道。
- 混沌： 对于足够大的耦合强度，系统表现出混沌特征，其回归映射展现出类似于马蹄映射（horseshoe maps）的拉伸和折叠结构。
Ott–Antonsen 假设的失效： 研究表明，断言所有傅里叶模满足 $\hat{f}_m = (\hat{f}_1)^m$ 的标准 OA 假设在离散时间动力学中并不保持。虽然该假设在阈值附近成立，但在进入同步机制深处后会出现显著偏差。这表明，离散时间系统无法像其连续时间对应物那样，被约化为一个控制序参数的简单一维流形。

意义与主张
本文声称，用离散时间非线性映射取代连续时间演化，从根本上改变了同步的物理特性。作者强调，离散时间同步构成了一个独特的动力学机制，其特征在于：

不同的稳定性条件： 映射的稳定性准则与流（flows）不同，这导致了修改后的同步阈值以及在排斥耦合下实现同步的可能性。
复杂的集体动力学： 与通常趋向于不动点的连续时间模型不同，离散时间模型支持集体倍周期、三周期以及混沌状态。
低维约化的局限性： 离散情形下 OA 假设的失效表明，在连续时间 Kuramoto 系统中有效的便捷低维描述并不能推广到离散时间映射中。

这项工作为系统研究离散时间同步提供了解析基础，并作为理解时间离散化如何影响耦合振子系统中集体现象的基准。

技术摘要：带有相位滞后的离散时间 Kuramoto 模型

类似论文