[论文解读] A toric deformation method for solving Kuramoto equations

Tianran Chen, Robert Davis|arXiv (Cornell University)|Oct 12, 2018

Nonlinear Dynamics and Pattern Formation参考文献 34被引用 8

一句话总结

本文提出了一种环形形变同伦方法，通过利用邻接多面体的正则单模三角剖分，高效计算柯尔莫哥洛夫模型在环形网络中的所有频率锁相配置。该方法避免了昂贵的混合体积计算，采用整数提升以保证数值稳定性，并精确追踪 $ N\binom{N-1}{\floor{(N-1)/2}} $ 条同伦路径，与已知解的上界完全一致。

ABSTRACT

The study of frequency synchronization configurations in Kuramoto models is a ubiquitous mathematical problem that has found applications in many seemingly independent fields. In this paper, we focus on networks whose underlying graph are cycle graphs. Based on the recent result on the upper bound of the frequency synchronization configurations in this context, we propose a toric deformation homotopy method for locating all frequency synchronization configurations with complexity that is linear in this upper bound. Loosely based on the polyhedral homotopy method, this homotopy induces a deformation of the set of the synchronization configurations into a series of toric varieties, yet our method has the distinct advantage of avoiding the costly step of computing mixed cells. We also explore the important consequences of this homotopy method in the context of direct acyclic decomposition of Kuramoto networks and tropical stable intersection points for Kuramoto equations.

研究动机与目标

• 开发一种高效的同伦延续方法，用于求解环形网络上的柯尔莫哥洛夫方程，其复杂度与孤立解的已知上界呈线性关系。
• 显式构造与环形图相关的邻接多面体的正则单模三角剖分，强化先前关于体积和面的结果。
• 改进柯尔莫哥洛夫网络的有向无环分解方案，使其可分解为具有唯一锁相解的原始子网络。
• 在热带代数几何中解释结果，表明由柯尔莫哥洛夫系统定义的热带超曲面恰好具有 $ N\binom{N-1}{\floor{(N-1)/2}} $ 个稳定自交点，每个重数为一。
• 建立环形网络中邻接多面体剖分、有向无环分解与热带稳定交点之间的三重等价关系。

提出的方法

• 该方法采用多面体同伦框架，但用基于邻接多面体 $ P_{C_N} $ 的正则单模三角剖分的环形形变替代标准的混合单元计算。
• 采用特殊赋值 $ \operatorname{val}(c_i^R) = 0 $ 和 $ \operatorname{val}(a_{ijk}^R) = 1 $ 或 $ 2 $，取决于奇偶性和边，该赋值诱导出具有稳定、重数为一的自交点的热带化。
• 初始系统由三角剖分导出，可在 $ O(N^2) $ 串行时间或 $ O(N) $ 并行时间内求解，避免了解二项式系统的需要。
• 同伦通过整数提升（取值于 $ \{0,1,2\} $）将柯尔莫哥洛夫系统形变为一组环道簇，确保数值稳定性，并避免随机提升带来的病态问题。
• 该算法精确追踪 $ N\binom{N-1}{\floor{(N-1)/2}} $ 条同伦路径，每条路径对应一个孤立解，由于三角剖分的正则性，不会出现路径追踪失败。
• 利用三角剖分改进有向无环分解，将网络分解为具有唯一复数锁相配置的原始子网络。

实验结果

研究问题

• RQ1 能否为环形网络上的柯尔莫哥洛夫方程构造一种同伦方法，使其以线性复杂度精确追踪已知孤立解的上界？
• RQ2 环形图的邻接多面体的正则单模三角剖分具有何种结构，其如何实现高效解的计算？
• RQ3 如何改进柯尔莫哥洛夫网络的有向无环分解，以实现分解为具有唯一解的原始子网络？
• RQ4 在特定赋值下，柯尔莫哥洛夫系统的热带解释是什么，其稳定交点如何与解的数量相关？
• RQ5 在环形网络的背景下，邻接多面体的正则剖分、有向无环分解与热带稳定交点之间是否存在三重等价关系？

主要发现

• 所提出的同伦方法以与上界 $ N\binom{N-1}{\floor{(N-1)/2}} $ 线性相关的复杂度，计算环形网络上柯尔莫哥洛夫方程的所有孤立解，与理论解数完全一致。
• 该方法避免了标准多面体同伦方法中昂贵的混合单元和混合体积计算，后者是主要性能瓶颈。
• 初始系统由三角剖分导出，可在 $ O(N^2) $ 时间内串行求解，或在 $ O(N) $ 时间内并行求解，显著提升了效率。
• 使用取值于 $ \{0,1,2\} $ 的整数提升确保了数值稳定性，消除了随机提升方案中常见的病态问题。
• 显式构造了邻接多面体 $ P_{C_N} $ 的正则单模三角剖分，为同伦与分解方法提供了几何基础。
• 在指定赋值下，柯尔莫哥洛夫系统的热带化恰好产生 $ N\binom{N-1}{\floor{(N-1)/2}} $ 个稳定自交点，每个重数为一，从几何上确认了解的数量。