$\Delta$ as a Dynamical Order Parameter for Non-Abelian Synchronization:\ the Case of $SU(2)$

Resumen Técnico: $\Delta$ como un Parámetro de Orden Dinámico para la Sincronización No Abeliana

Planteamiento del Problema
La sincronización en redes de osciladores acoplados ha sido analizada tradicionalmente mediante fases escalares y el formalismo de la Función de Estabilidad Maestra (MSF, por sus siglas en inglés), que vincula la estabilidad con los autovalores de la matriz de acoplamiento. Si bien trabajos recientes han identificado obstrucciones topológicas a la sincronización en señales de dimensión impar (específicamente en complejos simpliciales como tetraedros), los intentos previos de resolver estas obstrucciones utilizando conexiones $SU(2)$ planas se vieron limitados. Específicamente, los modelos previos utilizaron operadores de Dirac hermíticos, los cuales no lograron establecer correctamente la división de las ramas de la MSF. Además, mientras que los modelos de Kuramoto no abelianos y las generalizaciones cuaterniónicas han demostrado que las fases no conmutativas enriquecen la dinámica colectiva, se carecía de una medida cuantitativa para la "división no abeliana" específica de la MSF. Este artículo aborda la necesidad de cuantificar la división entre las ramas de Lyapunov izquierda y derecha en redes de osciladores con conexiones de gauge no abelianas.

Metodología
El estudio se centra en una red de $N_0=4$ vértices y $N_1=6$ aristas dispuestas como un tetraedro, donde cada arista porta un vector de estado en el espacio cuaterniónico $\mathbb{H} \cong \mathbb{R}^4$. Los autores emplean una conexión de gauge $SU(2)$ definida por cuaterniones unitarios asignados a las incidencias de las aristas.

1. Marco Matemático: La dinámica se modela utilizando un oscilador de Stuart-Landau cuaterniónico. El acoplamiento está gobernado por un operador de Dirac. Para investigar los efectos no abelianos, los autores construyen un operador de Dirac no hermítico $\tilde{D} = D_0 + \epsilon S$, donde $D_0$ es el operador estándar hermítico y $S$ es una perturbación de tipo skew-Hermitiana. Crucialmente, $S$ se construye de tal manera que $S\mathbf{1}=0$, asegurando que el manifold de sincronización permanezca invariante.
2. Linealización y Desacoplamiento: La ecuación variacional se linealiza alrededor de una órbita periódica. Los autores demuestran que, debido a que el operador de acoplamiento actúa sobre los índices de la red mientras que el término de deriva actúa sobre los índices de estado, el sistema se desacopla basándose en los autovectores de $\tilde{D}$.
3. Definición de $\Delta$: Para un autovalor complejo $\lambda$ del operador de Dirac, la dinámica reducida se divide en dos canales independientes:
   • Canal Izquierdo: Gobernado por matrices de multiplicación por la izquierda $\chi_L(\lambda)$.
   • Canal Derecho: Gobernado por matrices de multiplicación por la derecha $\chi_R(\lambda)$.
     El parámetro de orden escalar $\Delta$ se define como la diferencia absoluta máxima entre los exponentes de Lyapunov máximos de estos dos canales: $\Delta = \|\Lambda_L(\nu) - \Lambda_R(\nu)\|_\infty$, donde $\nu = \sigma|\lambda|$.

Contribuciones Clave y Resultados Teóricos
El artículo proporciona pruebas rigurosas que establecen el comportamiento de este sistema:

• Teorema de Desacoplamiento (Teorema 3.1): Se demuestra que para cualquier conexión $SU(2)$ que satisfaga $\tilde{D}\mathbf{1}=0$, el operador variacional deja invariantes los subespacios izquierdo y derecho de $\tilde{D}$. No existe un acoplamiento lineal entre los canales izquierdo y derecho, lo que permite analizarlos de forma independiente.
• Límite Abeliano (Teorema 3.2): Se demuestra que para grupos de gauge abelianos (por ejemplo, $U(1)^n$ o $SO(2)$), la multiplicación por la izquierda y por la derecha coinciden. Consecuentemente, la función de división $\Delta(\nu)$ es idénticamente cero, confirmando que $\Delta$ es una medida específica de los efectos no abelianos.
• Conjetura de Monotonía (Conjetura 3.3): Los autores proponen que, para grupos de Lie compactos con embebimientos compatibles, $\Delta$ es monótona con respecto a la inclusión de subgrupos (por ejemplo, $SU(2) \subset SU(3) \subset SU(4)$), aunque esto se señala como un tema para trabajos futuros.

Resultados Numéricos
Las simulaciones numéricas utilizando un oscilador de Stuart-Landau cuaterniónico ($\mu=1$, $\sigma=0.8$) confirman el marco teórico:

• Caso Hermítico: Con $\epsilon=0$ (operador $D_0$ hermítico), el sistema se sincroniza exponencialmente y el error cae a la precisión de la máquina, confirmando la eliminación de la obstrucción topológica.
• Caso No Hermítico: Con una extensión no hermítica ($\epsilon=0.2$), el sistema exhibe una división no nula. La división absoluta máxima se calculó como $\Delta \approx 0.0679$.
• Robustez: La dependencia de $\Delta$ respecto a la fuerza de la perturbación $\epsilon$ muestra que la división persiste incluso para componentes no hermíticas arbitrariamente pequeñas (por ejemplo, $\Delta \approx 0.001$ para $\epsilon=10^{-3}$). La división crece monótonamente con $\epsilon$.
• Análisis Espectral: Las visualizaciones del espectro indican que los valores más grandes de $|\Delta|$ se concentran en regiones donde los autovalores poseen partes imaginarias considerables, y la magnitud de la división aumenta sistemáticamente con $\epsilon$.

Significancia y Reivindicaciones
El artículo establece a $\Delta$ como una medida práctica y cuantitativa de los efectos no abelianos en sistemas dinámicos acoplados. Al demostrar el desacoplamiento de las ramas de Lyapunov izquierda y derecha y al probar que $\Delta > 0$ específicamente para conexiones no hermíticas y no abelianas, el trabajo proporciona un método para distinguir la sincronización no abeliana de los casos abelianos. Los autores posicionan este trabajo como una "prueba de concepto" de $SU(2)$ en un único red, allanando el camino para un ranking general de los grupos de Lie basado en sus espectros de sincronización. Sugieren posibles aplicaciones futuras para $\Delta$ como herramienta de diagnóstico de simetrías de gauge ocultas y como parámetro de diseño en control multi-agente, pero limitan explícitamente sus afirmaciones actuales a la validación teórica y numérica sobre el modelo del tetraedro.