A convex-geometric framework for fully phase-locked states in the finite Kuramoto model

1975 年提出的 Kuramoto 模型是解釋系統同步現象的核心。近期物理學者利用凸幾何架構，在頻率空間建構具 2n 個頂點的凸多面體，成功推導出有限數量振盪器達成完全相位鎖定所需最小耦合強度的數學上限公式。

1975 年 Kuramoto 模型的有限維度挑戰

在自然科學與工程領域中，互動振盪單元的同步現象無所不在，從生物的晝夜節律、神經元活動，到化學振盪與現代電力網格動態，都能見到其蹤影。為了以數學描述這些複雜的集體行為，物理學家在 1975 年提出了著名的 Kuramoto 模型。這個模型將每個獨立單元視為一個相位變數，並透過網路決定的耦合項來編碼彼此的交互作用，進而推導出豐富的巨觀行為，包含相位鎖定、部分同步與複雜的巨觀相變現象。

過去多數關於 Kuramoto 模型的理論研究，多半聚焦於熱力學極限（thermodynamic limit，即假設系統內的振盪器數量趨近於無限大），透過平均場理論來分析非相干狀態的穩定性喪失以及巨觀同調性的出現。然而，針對有限數量（finite-size）振盪器構成的系統，要找出其完全相位鎖定平衡點的存在條件與穩定性，本質上是一個高度非線性的固定點（fixed-point）數學難題。

當我們在高維度的環面上探討這組一階微分方程式時，系統的演化取決於每個振盪器的自然頻率、當前相位以及一個關鍵的實數參數——耦合強度 K。這項研究的核心目標，就是要精確找出在異質自然頻率下，系統能夠維持穩定完全相位鎖定狀態所需的最小耦合強度 $K_\ell$。

消除零特徵值的 n 維共轉座標系降維轉換

在原始的方程式中，只要將所有振盪器的相位加上一個相同的常數角度，系統依然會處於平衡狀態。這種平移不變性導致穩定的解並非孤立存在，反映在數學上，就是系統的 Jacobian 矩陣永遠具有一個特徵值為零，這讓穩定性分析變得異常棘手。

為了消除這個因均勻相位偏移產生的簡併現象，研究團隊引入了一個數學技巧：將視角切換到共轉座標系（co-rotating frame，一種藉由相對視角消除整體平移的參考系）。藉由將變數替換為相對相位，系統成功將變數數量從 $n+1$ 減少到 $n$ 個。在這個降維後的旋轉框架中，平衡解變成了孤立點，且數學論證證實這類平衡解的數量是有限的。

在這個簡化的 n 維系統中，當且僅當其 Jacobian 矩陣是負定矩陣（negative definite，指所有特徵值皆小於零的矩陣）時，我們才稱該相位鎖定平衡為穩定狀態。論文中的 Theorem A 明確指出，當頻率向量不為零時，系統會呈現嚴格的二分法：當耦合強度 K 夠小，系統不存在任何平衡解；而當 K 夠大時，必然存在一個穩定的平衡解，且系統的同步序參數 R 會隨著 K 嚴格遞增，最終逼近完美同步的數值 1。

將 Jacobian 穩定區域映射為凸集幾何

要界定系統何時開始出現穩定解，必須精確找出臨界耦合強度 $K_\ell$。研究團隊採取了高度幾何化的視角，首先在相位空間中定義出一個區域 $\Omega$，這個區域包含了所有能讓 Jacobian 矩陣特徵值皆小於零的相對相位。過去的動力系統研究已經證明，$\Omega$ 是一個開放且單連通的集合。

接下來是本研究最核心的數學轉換：透過 Kuramoto 向量場的映射，將相位空間的穩定區域 $\Omega$ 轉移到頻率空間，形成一個新的集合。先前的文獻已經證明，這個映像集合必然是一個凸集（convex set）。

這個幾何結構直接決定了有限系統是否能夠達成同步。如果在給定頻率和耦合強度的情況下，重新縮放的頻率向量落在這個凸集內部，系統就存在穩定的完全相位鎖定狀態；反之則無法穩定。因此，臨界耦合 $K_\ell$ 的幾何意義，就等於是從原點出發的射線，與該凸集邊界發生第一次交會點的參數倒數。這項轉換將原本複雜的動力系統穩定性問題，優雅地變成純粹的幾何相交難題。

透過 2n 個頂點建構多面體逼近耦合邊界

儘管確立了幾何意義，但該凸集的確切邊界通常難以透過解析方法直接求得。為了提供一個具備實踐價值的解答，研究團隊決定利用凸幾何的拓樸特性，建構一個多面體（polytope，由平坦邊界構成的幾何實體）從內部逼近這個凸集，進而推導出臨界強度的數學上限。

團隊首先在穩定區域 $\Omega$ 的邊界上，透過分析手法精確找出了 2n 個被標記的極值點。接著，他們將這些點映射到頻率空間，並以這些點為實體頂點，建構出一個包含 $2^n$ 個面的高維凸多面體 $\mathcal{P}_n$。

由於多面體 $\mathcal{P}n$ 完全被包覆在原本的映像凸集內部，當一條射線向外延伸時，必定會先穿透多面體 $\mathcal{P}_n$ 的邊界，然後才碰到真實凸集的邊界。透過找出射線與多面體邊界交會點，並取其參數倒數 $K_b$，就能保證獲得一個嚴格的上限值 $K_b \geq K\ell$。這個方法不僅提供了外圍逼近的途徑，更讓高維超平面方程的係數可以被精確求解。

Theorem B 提出的顯式解析臨界耦合上限

基於上述的凸多面體幾何建構，研究在 Theorem B 中給出了一個能直接套用於一般頻率向量的顯式解析上限（closed-form upper bound）。這意味著工程師或研究人員在處理有限網路時，不必再依賴繁雜且耗時的微分方程式數值模擬，只需將已知參數代入公式即可進行評估。

根據 Theorem B 的嚴謹推導，將自然頻率透過特定比例轉換為基底座標，並依照大於零或小於零將索引分組，最終能計算出一個明確的 $K_b$ 數學方程式。更值得注意的是，如果在系統設計時加入一個額外假設：即所有降維後的相對頻率總和為零，這個上限公式將會大幅度簡化為純粹的絕對值總和除以 $n+1$。

雖然研究者也坦承，這個透過多面體逼近得出的上限，在純數值量化層面上並非隨時絕對緊密，但如果頻率向量剛好與多面體的頂點對齊，這個邊界值就會是完全精確的。這項推導不僅賦予了 Kuramoto 模型一個具有高度操作性的分析工具，也徹底揭示了維持穩定完全相位鎖定解背後隱藏的優美幾何結構。

透過將動力系統的非線性固定點難題映射為頻率空間的凸集幾何相交問題，這項框架成功為有限振盪系統的同步臨界點，建立了無需模擬即可求得的精確數學解析上限。