Phase locking and multistability in the topological Kuramoto model on cell complexes
Pith reviewed 2026-05-18 09:29 UTC · model grok-4.3
The pith
Phase-locked states in the topological Kuramoto model on cell complexes are organized by winding numbers on generalized independent cycles, with multistability arising only when boundaries contain at least five elements.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
The topological nonlinear Kirchhoff conditions characterize all phase-locked states of the topological Kuramoto model on cell complexes. These states are organized by winding numbers associated with generalized independent cycles. Using rings, Platonic solids, and regular simplices, the authors establish that boundaries must have at least five elements for multistability to arise and that independent winding numbers associated with lower- and higher-dimensional boundaries generate cascades of multistability across dimensions.
What carries the argument
Topological nonlinear Kirchhoff conditions that express phase-locked states through winding numbers on generalized independent cycles of the cell complex.
If this is right
- Phase locking on cell complexes is completely determined by the choice of independent cycles and their winding numbers.
- Multistability appears only when every relevant boundary contains at least five cells.
- Independent winding numbers on boundaries of successive dimensions produce chained multistable regimes.
- The framework applies uniformly to simplicial complexes and more general cell complexes.
- Collective dynamics become predictable from the combinatorial topology alone once the winding numbers are fixed.
Where Pith is reading between the lines
- The five-element threshold may limit multistability in biological or engineered networks whose higher-order cells are small, such as certain neural motifs or communication graphs.
- The same winding-number construction could be used to design oscillator networks with prescribed numbers of stable synchronization patterns by choosing the cell-complex topology in advance.
- Extending the model to weighted or directed cell complexes would test whether the multistability cascades survive when the underlying geometry is no longer regular.
Load-bearing premise
The topological nonlinear Kirchhoff conditions are assumed to exhaustively characterize all phase-locked states without additional dynamical constraints from the embedding or metric of the cell complex.
What would settle it
A numerical simulation or experiment on a cell complex whose boundary has only four elements that nevertheless exhibits multistable phase-locked states not predicted by the winding numbers would falsify the universal rule.
Figures
read the original abstract
Higher-order interactions fundamentally shape collective dynamics in oscillator networks. The topological Kuramoto model captures these effects by extending synchronization models to include interactions between cells of arbitrary dimension within simplicial and cell complexes. We introduce the topological nonlinear Kirchhoff conditions to characterize all phase-locked states of the topological Kuramoto model. These states are organized by winding numbers associated with generalized independent cycles, which quantify how phases wind around these cycles. Using rings, Platonic solids, and regular simplices as illustrative examples, we uncover a universal rule: boundaries must have at least five elements for multistability to arise. We further find that independent winding numbers associated with lower- and higher-dimensional boundaries generate cascades of multistability across dimensions. These results show how the topology and boundary structure of cell complexes influence phase locking and multistability, and provide a general framework for collective dynamics on cell complexes.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper claims to introduce the topological nonlinear Kirchhoff conditions that characterize all phase-locked states of the topological Kuramoto model on cell complexes. These states are organized by winding numbers associated with generalized independent cycles. Using examples of rings, Platonic solids, and regular simplices, the authors report a universal rule that boundaries must have at least five elements for multistability to arise, and that independent winding numbers on lower- and higher-dimensional boundaries generate cascades of multistability across dimensions.
Significance. If the central claims hold, the work is significant for supplying a topology-driven, parameter-free framework that links boundary structure directly to multistability in higher-order oscillator networks. The organization of equilibria by winding numbers on generalized cycles offers a systematic way to classify collective states on simplicial and cell complexes, extending classical Kuramoto theory beyond graphs.
major comments (2)
- [Introduction of topological nonlinear Kirchhoff conditions] The topological nonlinear Kirchhoff conditions are asserted to exhaustively characterize all phase-locked states. However, the underlying Kuramoto dynamics on a cell complex are generally weighted by edge/face lengths or embedding geometry; the manuscript provides no argument or test showing that these weights do not produce additional equilibria outside the reported winding-number organization. This assumption is load-bearing for the universal multistability rule.
- [Examples (rings, Platonic solids, simplices)] All concrete examples (rings, Platonic solids, regular simplices) employ uniform metrics. The claim that boundaries require at least five elements for multistability therefore rests on regular cases; the paper should examine at least one non-uniform or weighted complex to verify whether extra phase-locked states appear that violate the reported cross-dimensional cascades.
minor comments (2)
- [Abstract] The abstract could more explicitly separate the model definition from the new Kirchhoff conditions and the multistability rule.
- [Notation and definitions] Notation for generalized cycles and associated winding numbers would benefit from an explicit example equation showing their computation on a higher-dimensional boundary.
Simulated Author's Rebuttal
We thank the referee for the careful reading and constructive comments on our manuscript. We address each major comment point by point below, indicating the revisions we will implement.
read point-by-point responses
-
Referee: [Introduction of topological nonlinear Kirchhoff conditions] The topological nonlinear Kirchhoff conditions are asserted to exhaustively characterize all phase-locked states. However, the underlying Kuramoto dynamics on a cell complex are generally weighted by edge/face lengths or embedding geometry; the manuscript provides no argument or test showing that these weights do not produce additional equilibria outside the reported winding-number organization. This assumption is load-bearing for the universal multistability rule.
Authors: The topological nonlinear Kirchhoff conditions follow directly from setting the right-hand side of the topological Kuramoto equations to zero. Because the equations are written in terms of the coboundary operator, the resulting algebraic conditions on the phase differences are satisfied precisely when the net weighted sine contributions around every cycle vanish. The winding numbers label the distinct integer homology classes of these cycles; this labeling is a topological feature and therefore persists for arbitrary positive weights. We acknowledge that the current text does not spell out this generality explicitly nor supply a weighted numerical check. In the revised manuscript we will add a short general argument (in the main text or an appendix) showing that every equilibrium must obey the Kirchhoff conditions and is therefore organized by the winding numbers, together with one explicit numerical example on a weighted ring or simplicial complex confirming that no extraneous equilibria appear. revision: yes
-
Referee: [Examples (rings, Platonic solids, simplices)] All concrete examples (rings, Platonic solids, regular simplices) employ uniform metrics. The claim that boundaries require at least five elements for multistability therefore rests on regular cases; the paper should examine at least one non-uniform or weighted complex to verify whether extra phase-locked states appear that violate the reported cross-dimensional cascades.
Authors: We agree that the illustrative examples were chosen with uniform metrics for clarity. To test robustness we will insert a new subsection containing a non-uniform example (a ring with heterogeneous edge lengths, or a simplicial complex with perturbed face weights). In this example we will compute the admissible winding-number combinations, locate the corresponding equilibria, and verify that multistability still requires boundaries of size at least five and that the cross-dimensional cascades remain intact. Any minor quantitative shifts will be reported and discussed. revision: yes
Circularity Check
No circularity in derivation of multistability rule from topological conditions
full rationale
The paper introduces the topological nonlinear Kirchhoff conditions as a characterization of phase-locked states in the topological Kuramoto model on cell complexes, then derives the universal rule that boundaries require at least five elements for multistability along with cascades generated by independent winding numbers on lower- and higher-dimensional boundaries. This follows from direct analysis of the model equations applied to explicit examples including rings, Platonic solids, and regular simplices. No step reduces the claimed results to a fitted parameter, self-citation loop, or input by construction; the findings emerge from the topological structure without presupposing the multistability threshold or cascade behavior.
Axiom & Free-Parameter Ledger
axioms (1)
- domain assumption The topological Kuramoto model on cell complexes is well-defined by extending pairwise phase interactions to arbitrary-dimensional cells.
Lean theorems connected to this paper
-
IndisputableMonolith/Foundation/AlexanderDuality.leanalexander_duality_circle_linking unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
boundaries must have at least five elements for multistability to arise; independent winding numbers associated with lower- and higher-dimensional boundaries generate cascades of multistability across dimensions
What do these tags mean?
- matches
- The paper's claim is directly supported by a theorem in the formal canon.
- supports
- The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
- extends
- The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
- uses
- The paper appears to rely on the theorem as machinery.
- contradicts
- The paper's claim conflicts with a theorem or certificate in the canon.
- unclear
- Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.
Reference graph
Works this paper leans on
-
[1]
Parametrize the set of solutions candidates as(16) while taking into account the constraints|ψ[±]| ≤1
-
[2]
Choose a partitionS [n+1] =S • [n+1] ∪S ◦ [n+1] and S[n−1] =S • [n−1] ∪S ◦ [n−1]
-
[3]
Find the solutions of the set of equations(22)where z[±] is a vector of integers
-
[4]
We briefly discuss three important aspects of this al- gorithmic approach to phase locking
Check the stability from the JacobianJif necessary via Eq.(10). We briefly discuss three important aspects of this al- gorithmic approach to phase locking. First, phase locked states withS ◦ =∅are always linearly stable and as such are especially relevant for applications. States with S◦ ̸=∅are typically, but not always unstable. In fact, we will demonstr...
-
[5]
Hence, the set of solutions candidates (16) is given by ψ[+] =ψ [+] sp ,ψ [−] =ψ [−] sp +K −1
is spanned by the harmonic vector1. Hence, the set of solutions candidates (16) is given by ψ[+] =ψ [+] sp ,ψ [−] =ψ [−] sp +K −1
-
[6]
1ζ [−].(24) The phase conditions (22) reduce to 1⊤ f [n−1](ψ[−]) = 2πz [−].(25) We start by exploring the fully symmetric case when ω= 0, due to whichψ [+] = 0andψ [−] sp = 0. In Fig. 3 we plot the expression1 ⊤ f [n−1](ψ[−])/(2π)as a func- tion of the parameterζ [−], for each of the partitions S[0] =S •
-
[7]
Since there are six elements that can either be•or◦, we obtain 64 different partitions ofS[0]
∪S ◦ [0]. Since there are six elements that can either be•or◦, we obtain 64 different partitions ofS[0]. Due to symmetry, these 64 curves in Fig. 3 are grouped according to the 7 different multiset permutations (i.e. fixed size of both partitions). Phase locked states are found where the curves cross an integer value. All of the partitions, except those w...
-
[8]
= [5,6], withω= (0,0,0,0,+ω 0,−ω 0). We plot the phase conditions (25) divided by2πas a function ofζ[−] 6 1.00 0.75 0.50 0.25 0.00 0.25 0.50 0.75 1.00 [ ] 1.5 1.0 0.5 0.0 0.5 1.0 1.5 1 f0( [ ])/(2 ) 0 = 0.00 0 = 0.50 0 = 1.00 0 = 1.35 0 = 1.70 0 = 1.304 0 = 1.318 0 = 1.329 0.16 0.14 0.12 0.10 0.08 0.04 0.02 0.00 0.02 1.307 1.311 1.315 1.319 1.323 1.327 0 ...
-
[9]
=∅(solid lines) andS◦
-
[10]
Filled (empty) symbols mark stable (unstable) solutions
= [5,6](dashed lines) at selected values ofω0 in ω= (0,0,0,0, ω 0,−ω 0). Filled (empty) symbols mark stable (unstable) solutions. The number of intersections with inte- ger values decreases asω0 increases, i.e. phase locked states gradually disappear. The inset highlights a narrow interval of ω0 whereS ◦
-
[11]
= [5,6]crosses zero twice, producing one sta- ble and one unstable solution, whileS◦
-
[12]
Bottom row: Dependence of the corresponding phase- locked states onω0
=∅does not cross it. Bottom row: Dependence of the corresponding phase- locked states onω0. The stable solutionS◦
-
[13]
=∅(orange line) vanishes nearω 0 ≈1.309, giving rise to the stable solution S◦
-
[14]
Stable and unstable (dashed blue line) solutions ofS◦
= [5,6](solid blue line). Stable and unstable (dashed blue line) solutions ofS◦
-
[15]
= [5,6]annihilate in a saddle-node bifurcation nearω 0 ≈1.327. for chosenω0 values in Fig. 4. As before, we find that the former partitioning always gives rise to stable solutions, while the latter typically yields unstable ones. However, nearthebifurcationpoint, thephase-lockedstateS ◦
-
[16]
Increasing ω0 further, the stableS ◦
= [5,6]without losing stability. Increasing ω0 further, the stableS ◦
-
[17]
= [5,6]solution annihilates with another unstableS ◦
-
[18]
This scenario is shown in the bottom panel of Fig
= [5,6]solution in a saddle- node bifurcation. This scenario is shown in the bottom panel of Fig. 4. Finally, letusgeneralizeourresulttopolygonalfacesof arbitrary size. We consider only the all-normal partitions S•
-
[19]
=S [0] withω= 0. If the ring size iss=N [0] = N[1] ≥3, the number of stable phased-locked solutions 5 10 15 20 25 30 Ring size s 1 3 5 7 9 11 13 15Nstable Simulation 1 + 2 s 1 4 FIG. 5.Dependence of the number of stable phase-locked states on ring size.The number of stable solutionsN stable as a function of ring sizesin an all-normal partitioned ring. Sym...
-
[20]
Moreover, we once again setω=0 as our default study case
=S [0]. Moreover, we once again setω=0 as our default study case. The boundary matrices are given in [53]. Following algorithm 1 we first parametrize the set of solution candidates. We find that the kernels of the boundary operators are one-dimensional and spanned by the vectors1of appropriate dimension. Forω= 0, K[n] = diag(1 N[0])andK [n+1] = diag(1 N[2...
-
[21]
For definiteness, we once again setω= 0and all couplings to unity
=S [0], for all viablen values for simplexes up to 5 dimensions. For definiteness, we once again setω= 0and all couplings to unity. The explicit form of the boundary operators and the kernel matrices are given in [53]. Our results are summarized in Fig. 8. Note that whenever there is one stable solution, it refers to the trivial solution with winding numb...
-
[22]
A. Pikovsky, M. Rosenblum, and J. Kurths,Synchroniza- tion: A Universal Concept in Nonlinear Sciences(Cam- bridge University Press, 2001)
work page 2001
-
[23]
F. Dörfler and F. Bullo, Synchronization in complex net- works of phase oscillators: A survey, Automatica50, 1539 (2014)
work page 2014
-
[24]
Y. Kuramoto, Self-entrainment of a population of cou- pled non-linear oscillators, inInternational symposium on mathematical problems in theoretical physics: January 23–29, 1975, kyoto university, kyoto/Japan(Springer,
work page 1975
-
[25]
J. A. Acebrón, L. L. Bonilla, C. J. Pérez Vicente, F. Ri- tort, and R. Spigler, The kuramoto model: A simple paradigm for synchronization phenomena, Reviews of modern physics77, 137 (2005)
work page 2005
-
[26]
S.H.Strogatz,Fromkuramototocrawford: exploringthe onset of synchronization in populations of coupled oscil- lators, Physica D: Nonlinear Phenomena143, 1 (2000)
work page 2000
-
[27]
F. A. Rodrigues, T. K. D. Peron, P. Ji, and J. Kurths, The kuramoto model in complex networks, Physics Re- ports610, 1 (2016)
work page 2016
- [28]
-
[29]
R. Delabays, T. Coletta, and P. Jacquod, Multistability of phase-locking and topological winding numbers in lo- cally coupled kuramoto models on single-loop networks, Journal of Mathematical Physics57, 032701 (2016)
work page 2016
-
[30]
S. Jafarpour, E. Y. Huang, K. D. Smith, and F. Bullo, Flow and elastic networks on then-torus: Geometry, analysis, and computation, SIAM Review64, 59 (2022)
work page 2022
-
[31]
F. Hellmann, P. Schultz, P. Jaros, R. Levchenko, T. Kap- itaniak, J. Kurths, and Y. Maistrenko, Network-induced multistability through lossy coupling and exotic solitary states, Nature communications11, 592 (2020)
work page 2020
-
[32]
M. Breakspear, S. Heitmann, and A. Daffertshofer, Gen- erativemodelsofcorticaloscillations: neurobiologicalim- plications of the Kuramoto model, Frontiers in Human Neuroscience4, 190 (2010)
work page 2010
-
[33]
A. E. Sizemore, C. Giusti, A. Kahn, J. M. Vettel, R. F. Betzel, and D. S. Bassett, Cliques and cavities in the human connectome, Journal of Computational Neuro- science44, 115 (2018)
work page 2018
- [34]
-
[35]
M. Andjelković, B. Tadić, and R. Melnik, The topology of higher-order complexes associated with brain hubs in human connectomes, Scientific Reports10, 17320 (2020)
work page 2020
- [36]
-
[37]
E. Estrada and G. J. Ross, Centralities in simplicial complexes. applications to protein interaction networks, Journal of Theoretical Biology438, 46 (2018)
work page 2018
-
[38]
A. P. Millán, H. Sun, L. Giambagli, R. Muolo, T. Car- letti, J. J. Torres, F. Radicchi, J. Kurths, and G. Bian- coni, Topology shapes dynamics of higher-order net- works, Nature Physics21, 353 (2025)
work page 2025
-
[39]
F. Battiston, G. Cencetti, I. Iacopini, V. Latora, M. Lu- cas, A. Patania, J.-G. Young, and G. Petri, Networks beyond pairwise interactions: Structure and dynamics, Physics reports874, 1 (2020)
work page 2020
-
[40]
F. Battiston, E. Amico, A. Barrat, G. Bianconi, G. Fer- raz de Arruda, B. Franceschiello, I. Iacopini, S. Kéfi, V. Latora, Y. Moreno,et al., The physics of higher-order interactions in complex systems, Nature Physics17, 1093 (2021)
work page 2021
-
[41]
G. Bianconi,Higher-Order Networks: An Introduction to Simplicial Complexes(Cambridge University Press, 2021)
work page 2021
-
[42]
D. Mulder and G. Bianconi, Network geometry and com- 12 plexity, Journal of Statistical Physics173, 783 (2018)
work page 2018
-
[43]
C. Bick, T. Gross, H. A. Harrington, and M. T. Schaub, What are higher-order networks?, SIAM Review65, 686 (2023)
work page 2023
-
[44]
I. Iacopini, G. Petri, A. Barrat, and V. Latora, Simplicial models of social contagion, Nature Communications10, 2485 (2019)
work page 2019
-
[45]
A. Santoro, F. Battiston, M. Lucas, G. Petri, and E. Am- ico, Higher-order connectomics of human brain function reveals local topological signatures of task decoding, in- dividual identification, and behavior, Nature Communi- cations15, 10244 (2024)
work page 2024
- [46]
- [47]
- [48]
- [49]
-
[50]
Z. Wu, G. Menichetti, C. Rahmede, and G. Bianconi, Emergent complex network geometry, Scientific Reports 5, 10073 (2015)
work page 2015
-
[51]
M. Andjelković, B. Tadić, M. Mitrović Dankulov, M. Ra- jković, and R. Melnik, Topology of innovation spaces in theknowledgenetworksemergingthroughquestions-and- answers, PLOS ONE11, e0154655 (2016)
work page 2016
- [52]
- [53]
- [54]
-
[55]
A. P. Millán, J. J. Torres, and G. Bianconi, Explosive higher-orderkuramotodynamicsonsimplicialcomplexes, Physical Review Letters124, 218301 (2020)
work page 2020
-
[56]
R. Ghorbanchian, J. G. Restrepo, J. J. Torres, and G. Bianconi, Higher-order simplicial synchronization of coupled topological signals, Communications Physics4, 120 (2021)
work page 2021
-
[57]
M. Nurisso, A. Arnaudon, M. Lucas, R. L. Peach, P. Ex- pert, F. Vaccarino, and G. Petri, A unified framework for simplicial kuramoto models, Chaos34, 053118 (2024)
work page 2024
-
[58]
P. S. Skardal and A. Arenas, Higher order interactions in complex networks of phase oscillators promote abrupt synchronization switching, Communications Physics3, 218 (2020)
work page 2020
- [59]
-
[60]
T. Carletti, L. Giambagli, and G. Bianconi, Global topo- logical synchronization on simplicial and cell complexes, Physical review letters130, 187401 (2023)
work page 2023
-
[61]
P. S. Skardal and A. Arenas, Abrupt desynchronization and extensive multistability in globally coupled oscillator simplexes, Physical Review Letters122, 248301 (2019)
work page 2019
-
[62]
L. DeVille, Consensus on simplicial complexes: Re- sults on stability and synchronization, Chaos31, 023137 (2021)
work page 2021
-
[63]
L. V. Gambuzza, F. D. Patti, L. Gallo, S. Lepri, M. Ro- mance, R. Criado, M. Frasca, V. Latora, and S. Boc- caletti, Stability of synchronization in simplicial com- plexes, Nature Communications12, 1255 (2021)
work page 2021
-
[64]
F. Parastesh, M. Mehrabbeik, K. Rajagopal, S. Jafari, M. Perc, C. I. del Genio, and S. Boccaletti, Synchro- nization stability in simplicial complexes of near-identical systems, Physical Review Research7, 033039 (2025)
work page 2025
- [65]
-
[66]
T. M. Roddenberry, M. T. Schaub, and M. Hajij, Sig- nal processing on cell complexes, inIEEE International Conference on Acoustics, Speech and Signal Processing (ICASSP)(2022) pp. 8852–8856
work page 2022
-
[67]
L. Giambagli, L. Calmon, R. Muolo, T. Carletti, and G. Bianconi, Diffusion-driven instability of topological signals coupled by the dirac operator, Physical Review E106, 064314 (2022)
work page 2022
-
[68]
S. Barbarossa and S. Sardellitti, Topological signal pro- cessing over simplicial complexes, IEEE Transactions on Signal Processing68, 2992 (2020)
work page 2020
-
[69]
S. H. Strogatz,Nonlinear dynamics and chaos with stu- dent solutions manual: With applications to physics, bi- ology, chemistry, and engineering(CRC press, 2018)
work page 2018
- [70]
-
[71]
C. Hartmann, P. C. Böttcher, D. Gross, and D. Wit- thaut, Synchronized states of power grids and oscillator networks by convex optimization, PRX energy3, 043004 (2024)
work page 2024
-
[72]
I. Bačić, C. Hartmann, P. C. Böttcher, D. Gross, A. Be- nigni, and D. Witthaut, Existence of phase-cohesive solutions of the lossless power-flow equations, TechRxiv Preprint DOI:10.36227/techrxiv.174611978.89382822/ (2025)
-
[73]
Symmetry breaking in minimum dissipation networks
A. Parameswaran, I. Bačić, A. Benigni, and D. Witthaut, Symmetry breaking in minimum dissipation networks, arXiv preprint arXiv:2505.24818 (2025)
work page internal anchor Pith review Pith/arXiv arXiv 2025
-
[74]
I. Bačić, Topological kuramoto model on cell complexes, https://github.com/ibfzj/topological-kuramoto/ (2025)
work page 2025
- [75]
-
[76]
A. E. Motter, S. A. Myers, M. Anghel, and T. Nishikawa, Spontaneous synchrony in power-grid networks, Nature Physics9, 191 (2013)
work page 2013
-
[77]
R. Olfati-Saber, J. A. Fax, and R. M. Murray, Consen- sus and cooperation in networked multi-agent systems, Proceedings of the IEEE95, 215 (2007)
work page 2007
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.