Higher-order interactions for controlling time-delayed Kuramoto model

Resumen Técnico: Interacciones de orden superior para el control del modelo de Kuramoto con retardo temporal

Planteamiento del problema
La sincronización colectiva de sistemas de osciladores acoplados es fundamental para diversas aplicaciones de ingeniería y biología, incluyendo la dinámica neuronal, las redes eléctricas y la robótica de enjambres. Sin embargo, las redes del mundo real a menudo involucran retardos temporales debido a las velocidades finitas de propagación de la señal. Estos retardos convierten al sistema en uno de dimensión infinita (descrito por ecuaciones diferenciales con retardo), lo que hace que su análisis y control sean significativamente desafiantes. Si bien los retardos temporales pueden inducir comportamientos complejos como la bistabilidad y la histéresis, el modelado matemático convencional a menudo los descuida o se basa en aproximaciones de desfase de primer orden (por ejemplo, los modelos de Sakaguchi–Kuramoto). Estas aproximaciones de pares convencionales no logran reproducir las características dinámicas distintivas de los sistemas con retardo, específicamente la bistabilidad entre los estados sincronizados e incoherentes.

Metodología
Los autores proponen un marco de control libre de retardo para el modelo de Kuramoto con retardo temporal, interpretando las interacciones de pares con retardo como interacciones efectivas de orden superior (no de pares). La metodología procede a través de los siguientes pasos:

1. Aproximación de orden superior: Partiendo de un sistema de $N$ osciladores de fase acoplados globalmente con interacciones con retardo temporal y una entrada de control, los autores aproximan los términos con retardo para una fuerza de acoplamiento $\epsilon$ pequeña. Esto produce un modelo de Kuramoto libre de retardo que contiene tanto términos de interacción de pares como términos de interacción de tres cuerpos (orden superior) efectivos, donde el retardo temporal $\tau$ aparece como desfases.
2. Reducción de dimensionalidad: Los autores aplican el ansatz de Ott–Antonsen (OA) a la función de densidad de probabilidad de las fases de los osciladores. Asumiendo una distribución Lorentziana de las frecuencias naturales (que se reduce a osciladores idénticos en el límite), derivan un sistema dinámico bidimensional para el parámetro de orden complejo $z(t) = R(t)e^{i\Psi(t)}$.
3. Promediación: Para aislar la dinámica lenta del grado de sincronización $R$, los autores emplean el método de promediación de segundo orden. Esto elimina la dinámica de fase rápida, resultando en una ecuación diferencial ordinaria (ODE) unidimensional reducida para el parámetro de orden $R$. Esta ecuación reducida retiene explícitamente términos de orden $O(\tau\epsilon^2)$, que capturan los efectos de las interacciones de orden superior.
4. Diseño de control: Se diseñan dos esquemas de control por retroalimentación de campo medio basados en la ecuación reducida:
   • Retroalimentación lineal: $u(t) = C R(t) \cos \Psi(t)$.
   • Retroalimentación no lineal: $u(t) = C R(t) (R(t) - E) \cos \Psi(t)$, introduciendo un parámetro adicional $E$ para apuntar a niveles de sincronización específicos.
5. Validación: Los diagramas de estabilidad derivados de la ecuación reducida de orden superior se comparan con aquellos de una aproximación de pares convencional (que descuida los términos $O(\tau\epsilon^2)$). Estas predicciones se validan mediante simulaciones numéricas de las ecuaciones diferenciales con retardo originales.

Contribuciones clave y resultados

• Precisión predictiva superior: Las simulaciones numéricas demuestran que la aproximación de orden superior predice la dinámica del sistema original con retardo de manera más precisa que la aproximación de pares convencional. Esto se cuantifica utilizando el Error Cuadrático Medio Global (RMSE), que es consistentemente menor para el modelo de orden superior a través de diferentes valores de retardo ($\tau$).
• Realización de la bistabilidad: El estudio reproduce con éxito la bistabilidad entre el estado incoherente ($R=0$) y el estado totalmente sincronizado ($R=1$) utilizando el control de retroalimentación lineal. Crucialmente, este régimen bistable es capturado únicamente por la ecuación reducida de orden superior; la aproximación de pares falla al predecir este fenómeno, mostrando solo regiones monoestables.
• Estabilización de estados intermedios: El esquema de control de retroalimentación no lineal permite la estabilización de niveles de sincronización intermedios arbitrarios ($0 < R < 1$). Mediante el ajuste de los parámetros $C$ y $E$, el sistema puede estabilizarse en un valor de referencia deseado $R_{ref}$.
• Regímenes bistables complejos: El marco revela un régimen bistable donde el estado totalmente sincronizado ($R=1$) coexiste con un estado sincronizado intermedio ($R=R_{ref}$). Esta estructura dinámica compleja es predicha con precisión por el modelo de orden superior, pero está ausente en la aproximación de pares.

Significado
El artículo sostiene que interpretar los retardos temporales como interacciones de orden superior de dimensión finita proporciona una metodología poderosa para analizar y controlar sistemas de osciladores acoplados con retardo. El marco propuesto ofrece un modelo reducido libre de retardo que retiene la precisión predictiva esencial del sistema original de dimensión infinita. Este enfoque facilita simulaciones numéricas más rápidas y permite el diseño analítico de estrategias de control por retroalimentación que pueden realizar diversos estados de sincronización —incluyendo la bistabilidad y la sincronización intermedia— que de otro modo serían difíciles de lograr o predecir utilizando modelos de pares estándar. Los autores sugieren que esta perspectiva podría extenderse a otras clases de redes dinámicas con retardo temporal donde los efectos del retardo complican el análisis y el control.