Emergence of higher-order interactions in systems of coupled Kuramoto oscillators with time delay

Resumen Técnico: Emergencia de Interacciones de Orden Superior en Sistemas de Osciladores de Kuramoto Acoplados con Retardo Temporal

Planteamiento del Problema
Los sistemas complejos del mundo real suelen exhibir interacciones que son inherentemente de orden superior (de muchos cuerpos), extendiéndose más allá de las redes tradicionales de pares. Si bien se sabe que las interacciones de orden superior enriquecen los comportos dinámicos como la sincronización, la histéresis y la bistabilidad, sus orígenes físicos son diversos. Concurrentemente, los retardos temporales en el acoplamiento son ubicuos en sistemas que van desde neuronas acopladas hasta redes eléctricas. El problema central abordado en este estudio es si las interacciones de pares con retardo temporal pueden reinterpretarse eficazmente como interacciones de orden superior. Específicamente, los autores investigan si la compleja dinámica colectiva derivada de los retardos temporales en osciladores de Kuramoto puede ser capturada por un modelo sin retardo que contenga interacciones explícitas de tres cuerpos.

Metodología
Los autores comienzan con un modelo de Kuramoto acoplado globalmente que presenta retardos temporales uniformes ($\tau$) en las interacciones de pares. Emplean un enfoque de expansión perturbativa, asumiendo una fuerza de acoplamiento ($\epsilon$) y un retardo temporal ($\tau$) moderadamente pequeños.

1. Derivación: Al expandir el término de acoplamiento de pares con retardo al segundo orden ($O(\epsilon^2)$), los autores derivan un sistema de ecuaciones diferenciales ordinarias sin retardo. El modelo derivado retiene los términos de pares originales pero introduce términos explícitos de interacción de tres cuerpos.
2. Simulación Numérica: Se comparan la dinámica del modelo original con retardo temporal, el modelo de orden superior derivado y un modelo de pares estándar (sin retardos ni términos de orden superior) mediante simulaciones numéricas para una gran población de osciladores idénticos ($N=300$). El grado de sincronización se cuantifica mediante el parámetro de orden de Kuramoto ($R$).
3. Análisis Analítico: Para comprender la estabilidad macroscópica, los autores aplican el ansatz de Ott–Antonsen (OA) al modelo de orden superior derivado. Asumen una distribución lorentziana de frecuencias naturales (aproximando frecuencias idénticas en el límite $\gamma \to 0$) para derivar una ecuación de forma cerrada para el parámetro de orden complejo. Esto permite un análisis de estabilidad lineal de los puntos fijos y de los ciclos límite, generando diagramas de fase para mapear regiones de incoherencia, sincronía total y bistabilidad.

Contribuciones Clave y Resultados

• Emergencia de Términos de Orden Superior: El estudio demuestra que la expansión del acoplamiento de pares con retardo temporal produce naturalmente un modelo de Kuramoto de orden superior. El término $O(\epsilon)$ corresponde al acoplamiento de pares estándar con un desfase, mientras que el término $O(\epsilon^2)$ corresponde a interacciones de tres cuerpos. En el caso de osciladores idénticos, la interacción de tres cuerpos toma la forma de un armónico $(2, -1, -1)$.
• Consistencia Cualitativa: Las simulaciones numéricas revelan que el modelo con retardo temporal y el modelo de tres cuerpos derivado exhiben transiciones de sincronización cualitativamente consistentes. Ambos modelos muestran transiciones de la incoherencia a la sincronía total y, crucialmente, regiones de bistabilidad donde coexisten estados coherentes e incoherentes. En contraste, un modelo de Kuramoto de pares estándar sin retardos ni términos de orden superior no logra reproducir esta bistabilidad bajo las mismas condiciones.
• Mecanismo de Bistabilidad: El análisis confirma que la bistabilidad observada en el sistema con retardo temporal es explicada por las interacciones efectivas de tres cuerpos en el modelo derivado. El análisis de estabilidad mediante el ansatz OA muestra que las estructuras de bifurcación (bifurcaciones de Hopf subcrítica y transcritica) del modelo de orden superior se alinean con las del modelo con retardo temporal cuando el producto $\epsilon\tau$ es pequeño.
• Limitaciones: Los autores señalan que surgen discrepancias cuantitativas entre los dos modelos cuando $\epsilon\tau$ es grande ($\epsilon\tau > 1$), lo cual es esperado dada la truncación de la expansión en el segundo orden.

Significancia y Reivindicaciones
El artículo sostiene que los retardos temporales en el acoplamiento de pares pueden reformularse eficazmente como interacciones de orden superior (específicamente de tres cuerpos). Este hallazgo proporciona una nueva perspectiva sobre cómo los retardos en el acoplamiento moldean la dinámica colectiva, sugiriendo que los comportamientos complejos inducidos por los retardos (como la bistabilidad) no requieren estrictamente la inclusión explícita de términos de retardo si se utiliza un marco de interacción de orden superior apropiado.

Los autores enfatizan que su objetivo no es proporcionar un esquema de aproximación de alta precisión para todos los sistemas con retardo temporal, sino más bien yuxtaponer los dos tipos de interacción para revelar su similitud cualitativa. Postulan que esta interpretación de orden superior efectiva podría ser adaptable a otras topologías de red (por ejemplo, redes de teselado regular o redes aleatorias) más allá del caso de acoplamiento global, facilitando potencialmente el análisis y el control de sistemas que involucran retardos temporales. Sin embargo, reconocen que para sistemas que exhiben dinámicas más complejas que la bistabilidad, pueden ser necesarios términos de orden superior más allá de los de tres cuerpos, y que un criterio sistemático para determinar el orden de truncación requerido sigue siendo un problema abierto.