On the efficiency of pairwise Hamiltonian control to desynchronize the higher-order Kuramoto model

Technische Samenvatting: Over de efficiëntie van pairwise Hamiltonian-controle voor het desynchroniseren van het hogere-orde Kuramoto-model

Probleemstelling
Synchronisatie in gekoppelde oscillator-systemen is een alomtegenwoordig fenomeen met dubbelzinnige implicaties: het is gunstig in contexten zoals elektriciteitsnetten, maar nadelig in andere, zoals bij epileptische aanvallen. Terwijl traditionele controlesstrategieën vaak vertrouwen op pairwise-interacties (tweezijdige interacties), wijst recent bewijs erop dat veel natuurlijke en technische systemen hogere-orde (veel-lichamen) interacties vertonen. Deze interacties veranderen de systeemdynamica fundamenteel, met name door de lineaire stabiliteit van gesynchroniseerde toestanden te vergroten en tegelijkertijd hun attractiegebieden te verkleinen — een fenomeen dat wordt beschreven als "diepere maar kleinere" attractiegebieden.

Het centrale probleem dat wordt aangepakt, is hoe deze concurrerende effecten van hogere-orde interacties de efficiëntie van een minimaal invasieve pairwise-controlstrategie beïnvloeden die ontworpen is om dergelijke systemen te desynchroniseren. Specifiek onderzoeken de auteurs of de toegenomen lineaire stabiliteit (die desynchronisatie belemmert) of de verkleinde grootte van het attractiegebied (die desynchronisatie vergemakkelijkt) domineert wanneer een controleschema dat oorspronkelijk voor pairwise-netwerken is ontwikkeld, wordt toegepast op hogere-orde Kuramoto-modellen.

Methodologie
De studie maakt gebruik van een hoger-orde Kuramoto-model bestaande uit $N$ fase-oscillatoren met zowel pairwise ($K_1$) als drie-lichamen ($K_2$) interacties. Het systeem wordt beheerst door:
$$\dot{\theta}_i = \omega_i + \frac{K_1}{\langle k^{(1)} \rangle} \sum_{j} A^{(1)}_{ij} \sin(\theta_j - \theta_i) + \frac{K_2}{\langle k^{(2)} \rangle} \sum_{j,k} A^{(2)}_{ijk} \sin(\theta_j + \theta_k - 2\theta_i) + \mu \delta_i p_i$$
waarbij $\mu$ de controlsterkte is en $\delta_i$ aangeeft of knooppunt $i$ gecontroleerd wordt.

Het controlemechanisme is een minimaal invasieve pairwise Hamiltonian-controle, aangepast van Asllani et al. [35]. De controleterm $p_i$ is gedefinieerd als:
$$p_i = -\frac{1}{2K_1^2} \frac{M^2}{\langle k^{(1)} \rangle^2} R_M \tilde{R}_{M,i} \cos(\Psi_M - \tilde{\Psi}_{M,i})$$
Deze controle werkt alleen op een subset van $M$ gepinde knooppunten. Het meet het synchronisatieniveau van de gecontroleerde knooppunten en past een kracht toe die proportioneel is aan de synchronisatiemagnitudo, welke verdwijnt wanneer het systeem volledig gedesynchroniseerd is. Deze aanpak vermijdt de computationele complexiteit van het afleiden van hogere-orde controletermen.

De auteurs evalueren de controle-efficiëntie met behulp van twee ordeparameters:

1. Fasesynchronisatie ($\hat{R}$): De tijdgemiddelde Kuramoto-orde-parameter.
2. Frequentie-synchronisatie ($\hat{R}_{\dot{\theta}}$): Een maatstaf voor fasevergrendeling (frequentie-synchronisatie).

Numerieke experimenten werden uitgevoerd op twee topologieën:

• Hyperringen: Regelmatige structuren die 'twisted states' ondersteunen.
• Random Erdős-Rényi hypergrafen: Onregelmatige structuren zonder rotationele invariantie.

De studie varieerde systematisch de hogere-orde koppelsterkte ($K_2$), de controlsterkte ($\mu$), het aantal gecontroleerde knooppunten ($M$) en de initiële afstand ($\epsilon$) van de systeemtoestand tot het centrum van het synchronisatie-attractiegebied.

Belangrijkste Resultaten

1. Gesynchroniseerde Initiële Condities:
   Wanneer het systeem start nabij de volledig gesynchroniseerde toestand (het centrum van het attractiegebied), maakt het verhogen van de hogere-orde koppelsterkte $K_2$ desynchronisatie moeilijker. De kritieke controlsterkte ($\mu_c$) die vereist is om desynchronisatie te bereiken, neemt monotoon toe met $K_2$. Dit geeft aan dat in dit regime de toegenomen lineaire stabiliteit veroorzaakt door hogere-orde interacties domineert over het verkleinen van het attractiegebied. De pairwise-controle blijft effectief als $\mu$ voldoende groot is, maar de kosten van de controle stijgen met $K_2$.

2. Geperturbeerde Initiële Condities:
   Wanneer het systeem verder van de gesynchroniseerde toestand start (grotere $\epsilon$), treedt een niet-monotoon gedrag op.

   • Intermediaire $K_2$: Desynchronisatie is het moeilijkst; de controle is het minst efficiënt.
   • Grote $K_2$: Desynchronisatie wordt gemakkelijker; de controle-efficiëntie verbetert.
     Dit gedrag suggereert dat voor grote $K_2$, het verkleinen van het attractiegebied de dominante factor wordt, waardoor de controle trajecten gemakkelijker uit het gebied kan duwen ondanks de hogere lineaire stabiliteit.

3. Partiële Controle:
   De studie bevestigt dat desynchronisatie kan worden bereikt, zelfs wanneer slechts een fractie van de knooppunten ($M < N$) wordt gecontroleerd. Het minimale aantal vereiste knooppunten ($M_c$) neemt over het algemeen toe met $K_2$ voor gesynchroniseerde initiële condities, maar vertoont niet-monotoon gedrag voor geperturbeerde initiële condities, wat de trends volgt die bij de controlsterkte zijn waargenomen.

4. Frequentie versus Fase Desynchronisatie:
   De controle onderdrukt fasesynchronisatie effectief. Echter, frequentie-synchronisatie (bijv. in 'twisted' of 'cluster' states) is robuuster, met name voor gemiddelde $K_2$ op hyperringen, vanwege structurele symmetrieën en het feit dat de controle niet direct op frequentievergrendeling is gericht.

Betekenis en Claims
Het artikel claimt een duidelijker beeld te geven van de wisselwerking tussen hogere-orde interacties en controle-efficiëntie. De primaire bijdrage is het aantonen dat het "diepere maar kleinere" effect van hogere-orde interacties leidt tot verschillende controle-regimes afhankelijk van de initiële toestand:

• Nabij de gesynchroniseerde toestand fungeren hogere-orde interacties als een obstakel voor controle vanwege de verhoogde stabiliteit.
• Ver van de gesynchroniseerde toestand kunnen voldoende sterke hogere-orde interacties de controle juist faciliteren door het attractiegebied te verkleinen.

De auteurs stellen dat een eenvoudige, minimaal invasieve pairwise controlestrategie voldoende is om hogere-orde Kuramoto-systemen te desynchroniseren, mits de controle-intensiteit of het aantal gecontroleerde knooppunten adequaat is. Deze bevinding is significant omdat het suggereert dat complexe, computationeel dure hogere-orde controletermen (zoals voorgesteld in Ref. [100]) niet strikt noodzakelijk zijn voor desynchronisatie, wat een meer praktische implementatieroute biedt. De resultaten benadrukken het belang van het overwegen van zowel lineaire stabiliteit als de geometrie van het attractiegebied bij het ontwerpen van controlestrategieën voor complexe netwerken met hogere-orde interacties.