[论文解读] Efficient moment-based approach to the simulation of infinitely many heterogeneous phase oscillators

Iván León, Diego Pazó|arXiv (Cornell University)|Jan 31, 2023

Nonlinear Dynamics and Pattern Formation参考文献 43被引用 1

一句话总结

本文提出了一种基于矩的数值方法，用于高效模拟无限多个异质相位振子，通过相位的傅里叶-埃米特谱分解对振子密度进行表示。该方法将无限维动力学简化为有限维常微分方程组，从而能够准确计算集体态、分岔现象以及李雅普诺夫指数——尤其在混沌区域中，直接模拟因有限尺寸效应和计算成本过高而失效。

ABSTRACT

The dynamics of ensembles of phase oscillators are usually described considering their infinite-size limit. In practice, however, this limit is fully accessible only if the Ott-Antonsen theory can be applied, and the heterogeneity is distributed following a rational function. In this work, we demonstrate the usefulness of a moment-based scheme to reproduce the dynamics of infinitely many oscillators. Our analysis is particularized for Gaussian heterogeneities, leading to a Fourier-Hermite decomposition of the oscillator density. The Fourier-Hermite moments obey a set of hierarchical ordinary differential equations. As a preliminary experiment, the effects of truncating the moment system and implementing different closures are tested in the analytically solvable Kuramoto model. The moment-based approach proves to be much more efficient than the direct simulation of a large oscillator ensemble. The convenience of the moment-based approach is exploited in two illustrative examples: (i) the Kuramoto model with bimodal frequency distribution, and (ii) the "enlarged Kuramoto model" (endowed with nonpairwise interactions). In both systems, we obtain new results inaccessible through direct numerical integration of populations.

研究动机与目标

• 为克服大规模但有限振子群体直接数值模拟的局限性，后者受有限尺寸涨落和高计算成本的影响。
• 为研究具有连续异质性的全局耦合相位振子的热力学极限，开发一种数值高效的直接模拟替代方法。
• 实现对不稳定集体态和集体李雅普诺夫指数的计算，这些量在标准直接积分方法中无法获得。
• 将基于矩的方法的应用范围扩展至奥特-安东森框架之外，特别是针对非有理频率分布（如高斯分布）的情形。
• 为研究传统模拟方法失效的系统中的复杂集体现象（包括分岔和混沌动力学）提供一个稳健的框架。

提出的方法

• 将振子密度分解为相位方向的傅里叶模和自然频率方向的埃米特多项式基，得到傅里叶-埃米特谱表示。
• 推导出描述振子密度在热力学极限下演化的矩的常微分方程（ODE）层级系统。
• 对矩系统应用截断与闭合近似，其中线性闭合方法在测试案例中表现有效且可靠。
• 利用数值延拓技术访问不稳定集体态（如部分同步解），无论其稳定性如何。
• 通过矩系统计算集体李雅普诺夫指数，避免了有限模拟中正指数的O(N)缩放问题。
• 将该方法应用于两个模型系统：双模库朗托模型和具有非成对相互作用的扩展库朗托模型。

实验结果

研究问题

• RQ1 基于矩的方法是否能有效模拟超越奥特-安东森框架的无限多个异质相位振子的热力学极限？
• RQ2 对于具有高斯频率分布的系统，不同闭合方案（如线性闭合）在截断矩层级时的准确性和稳定性如何？
• RQ3 基于矩的方法是否能够访问直接模拟有限群体无法获得的不稳定集体态和分岔？
• RQ4 在热力学极限下，最大李雅普诺夫指数的收敛行为如何？其渐近值能否通过矩方法可靠估计？
• RQ5 在具有非成对相互作用的复杂系统（如扩展库朗托模型）中，基于矩的方法表现如何？

主要发现

• 基于傅里叶-埃米特分解的矩方法能够准确且高效地模拟具有高斯异质性的相位振子集合的无限尺寸极限。
• 矩层级的线性闭合提供了可靠且计算高效的近似，在包括解析可解的库朗托模型在内的多个测试案例中得到验证。
• 该方法成功揭示了双模库朗托模型中此前未被注意的漂移-叉分岔线，从而完善了已知的相图。
• 通过矩系统实现的数值延拓可访问不稳定的部分同步态，从而实现对分岔结构的完整表征。
• 在集体混沌中，最大李雅普诺夫指数的渐近值估计为 λ ≈ 1.26×10⁻⁴，其收敛行为满足 Λ(N)−λ ∝ N⁻⁰.⁶⁶，幂律行为通过矩方法模拟得到证实。
• 该方法揭示了扩展库朗托模型中 σ → 0 的奇异极限，此时系统表现出慢-快动力学，且慢相持续时间发散，与异质性趋于零时的解析预期一致。