Connecting Hodge and Sakaguchi-Kuramoto through a mathematical framework for coupled oscillators on simplicial complexes

The paper addresses how to extend the Kuramoto model of coupled oscillators to higher-order structures (simplicial complexes) while incorporating frustration (phase-lag) in a way that is orientation independent for higher-order simplices. Classical Kuramoto models on graphs capture synchronization phenomena but are limited to pairwise interactions, and the Sakaguchi-Kuramoto variant introduces a phase-lag that leads to rich dynamics. Recent interest in higher-order interactions motivates generalizations to simplicial complexes, either by keeping node phases with higher-order couplings or by promoting phases to higher-order simplices (edges, faces, etc.). The authors aim to (i) formulate a weighted simplicial Kuramoto model with phases on k-simplices, (ii) introduce linear and nonlinear frustration terms in an orientation-invariant way, and (iii) study the resulting dynamics—particularly on edges—using Hodge-theoretic tools. The significance lies in connecting Hodge theory to frustrated synchronization, unveiling how topology (holes, orientations) and higher-order couplings shape synchronization, desynchronization, and complex dynamics.

The introduction surveys extensive work on Kuramoto dynamics on graphs, including time delays, signed and time-varying interactions, and stochasticity. The Sakaguchi-Kuramoto model, with a frustration (phase-lag) parameter, produces chimera states and remote synchronization and appears in diverse applications. Higher-order interactions are modeled via hypergraphs and simplicial complexes, prompting generalizations of dynamical systems, including Kuramoto models with node phases but higher-order couplings and models with phases on higher-order simplices. Prior work demonstrated explosive synchronization in edge dynamics projected to nodes/faces under specific nonlinear couplings and later local couplings between orders. These provide the foundation for the present study, which adds weights and both linear and nonlinear frustration in a rigorous, orientation-independent framework.

• Mathematical framework: The authors use discrete differential geometry on simplicial complexes. They define weighted coboundary operators N_k = B_k^{k+1} W_{k+1} and Hodge Laplacians L_k = N_k^† N_k + N_{k-1} N_{k-1}^†, with weight matrices W_k (diagonal, positive).
• Weighted simplicial Kuramoto model: For k-cochains θ^{(k)}, the unfrustrated dynamics are θ̇^{(k)} = −N_{k−1}^† sin(N_{k−1} θ^{(k)}) − N_k sin(N_k θ^{(k)}), yielding linear consensus θ̇^{(k)} ≈ −L_k θ^{(k)} near ker(L_k). Internal frequencies arise from components in ker(L_k) via a rotating frame.
• Frustration (orientation-invariant): To generalize Sakaguchi-Kuramoto, they introduce lift matrices V_k and projections to positive/negative entries to construct a frustration operator F^{(k)}(N_k): x ↦ N_k x + a_k. The frustrated model becomes θ̇^{(k)} = −F^{(k−1)}(N_{k−1})[sin(N_{k−1} θ^{(k)})] − (N_k V_{k+1}^T) sin(V_{k+1} N_k θ^{(k)} + a_{k+1}). Here a_k is a linear frustration (natural frequency term) and a_{k+1} a nonlinear frustration (phase-lag). This construction is invariant to (k+1)-simplex orientation.
• Focus on edge dynamics (k=1): Main working equation θ̇^{(1)} = −a_1 N_0 sin(N_0 θ^{(1)}) − (N_1 V_2^T) sin(V_2 N_1 θ^{(1)} + a_2). The model is invariant to face (2-simplex) orientation but depends on edge orientation.
• Hodge decomposition: Decompose θ^{(k)} into gradient (Im N_{k−1}), harmonic (ker L_k), and curl (Im N_k) components using orthonormal projectors P_grad, P_harm, P_curl. Track their time evolutions and slopes (drifts) at stationarity.
• Simplicial order parameter (SOP): Define a generalized order parameter R_k^{(s)}(θ) as the average cosine over lifted differences (normalized by counts/weights of adjacent simplices) that equals 1 when θ^{(k)} lies in the harmonic space, generalizing full synchronization to harmonic phase locking.
• Numerical experiments: Simulate edge oscillators scanning frustration parameters (a_1, a_2), under different edge orientation configurations on:
  1. Single triangle with one face (no hole), including fully consistent orientation and with one edge flipped.
  2. Small complex with one hole; test effects of flipping specific edge orientations.
  3. Larger Delaunay triangulation with two holes.
  • Measurements: SOP (mean and std), slopes of projections onto gradient/curl/harmonic subspaces, and largest Lyapunov exponent (mean and quartiles across edges). Explore effect of face weight w in a two-triangle complex (one full, one empty) by varying w ∈ [0,1].
• Code availability: Simulations reproducible via provided GitHub/PyPI/Zenodo resources.

• Orientation-invariant frustration: The proposed lifting/projection framework yields a simplicial Sakaguchi-Kuramoto model whose nonlinear frustration term is independent of (k+1)-simplex orientation while enabling linear and nonlinear frustrations (a_k, a_{k+1}).
• Edge dynamics structure and coupling: With frustration, Hodge subspaces (gradient, harmonic, curl) become hierarchically coupled; the gradient dynamics can influence curl via the nonlinear term (Eqs. 32–33). In absence of nonlinear frustration (a_2=0), subspaces decouple.
• Single-face triangle (no hole): Despite zero-dimensional harmonic space, full synchronization (R=1) can still be achieved for nonzero frustrations (region in (a_1,a_2)-space). Edge orientation matters: a flipped edge produces rich regimes including non-constant gradient and curl components, Lissajous-like trajectories when curl is constant, and sharp transitions (in gradient drift) as parameters vary.
• Complex with a hole: Identify two critical nonlinear frustration values α2g (onset of non-constant gradient) and α2c (onset of non-constant curl), with α2c ≥ α2g. Empirically, curl becomes non-constant only when gradient is also non-constant. Edge orientation strongly alters the presence and order of transitions; some configurations exhibit re-synchronization (phase re-locking) at higher a_2.
• Larger complex (two holes): More complex dynamics with no re-synchronization at large a_2; standard deviation of SOP increases, and largest Lyapunov exponent becomes positive as soon as α2 > α2g and grows further for larger a_2, indicating chaos. Distinct chaotic regimes are suggested by changes in projection slopes and Lyapunov trends.
• Effect of weights: Varying face weight w modulates the influence of nonlinear frustration. For w→0 (face vanishes), dynamics lie in ker(L) (no nonlinear frustration). Increasing w produces regimes with non-vanishing gradient and, for larger w, non-vanishing curl projections. At w=1, behavior resembles the simple triangle case. Weights at different orders modulate coupling strengths and can control the impact of nonlinear frustration.
• General phenomena: Partial loss of synchronization aligned with Hodge subspaces; emergence of simplicial phase re-locking at high frustration in small, symmetric complexes; strong dependence on edge orientation and topology (holes) for synchronization transitions and complexity of dynamics.

The study connects Hodge-theoretic structure with the Sakaguchi-Kuramoto phase-lag in a rigorous higher-order framework, clarifying how orientation-invariant nonlinear frustration can be defined on simplicial complexes. By decomposing dynamics into Hodge subspaces, the authors show that frustration introduces hierarchical coupling between gradient and curl components, leading to rich transitions from phase-locked (harmonic) synchronization to non-stationary and chaotic regimes. The observed thresholds (α2g, α2c) and the rule that curl dynamics become non-stationary only when gradient dynamics do so, reveal structural constraints imposed by topology and the lifted coupling. Edge orientation and the presence/location of holes critically modulate synchronization loss, re-locking, and chaos. These insights suggest new control and design strategies for higher-order oscillator systems by tuning orientations, weights, and topology to achieve desired synchronization patterns or to avoid chaotic regimes.

The paper extends the simplicial Kuramoto model to include arbitrary weights and both linear (a_k) and nonlinear (a_{k+1}) frustration, yielding a simplicial Sakaguchi-Kuramoto model that is invariant to (k+1)-simplex orientation. It generalizes synchronization to harmonic phase locking via a simplicial order parameter and leverages Hodge decomposition to dissect dynamics. The work demonstrates diverse behaviors—from constant and phase-locked to chaotic dynamics—whose emergence depends strongly on edge orientation, topology (holes), and weights. Future research directions include: systematic exploration of weights and their constraints, comparison of consensus vs diffusion formulations, analysis of the interplay between linear and nonlinear frustrations, and a broader topological characterization (e.g., role of Betti numbers, hole localization, and symmetry) for predicting transitions between non-stationary and chaotic regimes, potentially enabling control over partial or cluster synchronizations.

• Analytical thresholds: Exact analytical determination of α2g and α2c is difficult and depends on both topology and orientation; results are largely numerical and example-based.
• Scope of examples: Detailed studies focus on small and one larger 2D complexes; generality across broader classes of complexes (dimensions, geometries) is not fully established.
• Weights: Effects of weights are explored phenomenologically in a simple example; a systematic analytical and computational study of weight choices and inter-order relations is left for future work.
• Interplay of frustrations: The combined effects of linear and nonlinear frustrations are not comprehensively explored.
• Formulation choice: The study emphasizes the consensus formulation; diffusion formulation differences in weighted settings are noted but not investigated.
• Orientation dependence: While the nonlinear term is orientation-invariant for faces, edge orientation still affects dynamics, complicating general predictions.