Higher-order interactions for controlling time-delayed Kuramoto model

Sintesi Tecnica: Interazioni di ordine superiore per il controllo del modello di Kuramoto con ritardo temporale

Definizione del Problema
La sincronizzazione collettiva dei sistemi di oscillatori accoppiati è fondamentale per varie applicazioni ingegneristiche e biologiche, tra cui la dinamica neuronale, le reti elettriche e la robotica di sciami. Tuttavia, le reti del mondo reale spesso comportano ritardi temporali dovuti alle velocità finite di propagazione del segnale. Questi ritardi rendono il sistema a dimensione infinita (descritto da equazioni differenziali a ritardo), rendendo l'analisi e il controllo significativamente complessi. Sebbene i ritardi temporali possano indurre comportamenti complessi come la bistabilità e l'isteresi, la modellazione matematica convenzionale spesso li trascura o si affida ad approssimazioni di primo ordine del ritardo di fase (ad esempio, i modelli di Sakaguchi–Kuramoto). Queste approssimazioni pairwise convenzionali non riescono a riprodurre le caratteristiche dinamiche tipiche dei sistemi ritardati, specificamente la bistabilità tra stati sincronizzati e incoerenti.

Metodologia
Gli autori propongono un framework di controllo privo di ritardo per il modello di Kuramoto con ritardo temporale, interpretando le interazioni pairwise con ritardo come interazioni effettive di ordine superiore (non pairwise). La metodologia procede attraverso i seguenti passaggi:

1. Approssimazione di ordine superiore: Partendo da un sistema di $N$ oscillatori di fase accoppiati globalmente con interazioni con ritardo temporale e un input di controllo, gli autori approssimano i termini ritardati per una piccola forza di accoppiamento $\epsilon$. Ciò produce un modello di Kuramoto privo di ritardo che contiene sia termini di interazione pairwise che termini di interazione a tre corpi (di ordine superiore) effettivi, dove il ritardo temporale $\tau$ appare come ritardi di fase.
2. Riduzione della dimensionalità: Gli autori applicano l'ansatz di Ott–Antonsen (OA) alla funzione di densità di probabilità delle fasi degli oscillatori. Assumendo una distribuzione Lorentziana delle frequenze naturali (che si riduce a oscillatori identici nel limite), derivano un sistema dinamico bidimensionale per l'ordine parametrico complesso $z(t) = R(t)e^{i\Psi(t)}$.
3. Media (Averaging): Per isolare la dinamica lenta del grado di sincronizzazione $R$, gli autori utilizzano il metodo di media del secondo ordine. Questo elimina la dinamica di fase veloce, risultando in un'equazione differenziale ordinaria (ODE) unidimensionale ridotta per l'ordine parametrico $R$. Questa equazione ridotta mantiene esplicitamente i termini di ordine $O(\tau\epsilon^2)$, che catturano gli effetti delle interazioni di ordine superiore.
4. Progettazione del Controllo: Due schemi di controllo feedback mean-field sono progettati basandosi sull'equazione ridotta:
   • Feedback Lineare: $u(t) = C R(t) \cos \Psi(t)$.
   • Feedback Nonlineare: $u(t) = C R(t) (R(t) - E) \cos \Psi(t)$, introducendo un parametro aggiuntivo $E$ per mirare a specifici livelli di sincronizzazione.
5. Validazione: I diagrammi di stabilità derivati dall'equazione ridotta di ordine superiore sono confrontati con quelli di un'approssimazione pairwise convenzionale (che trascura i termini $O(\tau\epsilon^2)$). Queste previsioni sono validate attraverso simulazioni numeriche del sistema originale di equazioni differenziali a ritardo temporale.

Contributi Chiave e Risultati

• Accuratezza Predittiva Superiore: Le simulazioni numeriche dimostrano che l'approssimazione di ordine superiore predice la dinamica del sistema ritardato originale con maggiore accuratezza rispetto all'approssimazione pairwise convenzionale. Ciò è quantificato utilizzando l'Errore Quadratico Medio (RMSE) globale, che è costantemente inferiore per il modello di ordine superiore attraverso diversi valori di ritardo ($\tau$).
• Realizzazione della Bistabilità: Lo studio riesce a riprodurre la bistabilità tra lo stato incoerente ($R=0$) e lo stato completamente sincronizzato ($R=1$) utilizzando il controllo feedback lineare. Fondamentalmente, questa regione bistabile è catturata solo dall'equazione ridotta di ordine superiore; l'approssimazione pairwise fallisce nel prevedere questo fenomeno, mostrando solo regioni monostabili.
• Stabilizzazione di Stati Intermedi: Lo schema di controllo feedback nonlineare consente di stabilizzare arbitrari livelli di sincronizzazione intermedi ($0 < R < 1$). Regolando i parametri $C$ ed $E$, il sistema può essere stabilizzato a un valore di riferimento desiderato $R_{ref}$.
• Regimi Bistabili Complessi: Il framework rivela un regime bistabile in cui lo stato completamente sincronizzato ($R=1$) coesiste con uno stato sincronizzato intermedio ($R=R_{ref}$). Questa complessa struttura dinamica è predetta accuratamente dal modello di ordine superiore, ma è assente nell'approssimazione pairwise.

Significatività
L'articolo sostiene che interpretare i ritardi temporali come interazioni di ordine superiore a dimensione finita fornisca una metodologia potente per analizzare e controllare sistemi di oscillatori accoppiati con ritardo. Il framework proposto offre un modello ridotto privo di ritardo che mantiene l'essenziale accuratezza predittiva del sistema originale a dimensione infinita. Questo approccio facilita simulazioni numeriche più veloci e consente la progettazione analitica di strategie di controllo feedback che possono realizzare diversi stati di sincronizzazione — inclusa la bistabilità e la sincronizzazione intermedia — che sarebbero altrimenti difficili da ottenere o prevedere utilizzando i modelli standard pairwise. Gli autori suggeriscono che questa prospettiva potrebbe essere estesa ad altre classi di reti dinamiche con ritardo temporale in cui gli effetti del ritardo complicano l'analisi e il controllo.