Bifurcaciones de soluciones retorcidas en un límite continuo para el modelo de Kuramoto en gráficos de vecinos más cercanos

Información Básica

• ID del Artículo: 2510.22663
• Título: Bifurcaciones de soluciones retorcidas en un límite continuo para el modelo de Kuramoto en gráficos de vecinos más cercanos
• Autor: Kazuyuki Yagasaki (Departamento de Matemática Aplicada y Física, Escuela de Graduados en Informática, Universidad de Kyoto)
• Clasificación: math.DS (Sistemas Dinámicos)
• Fecha de Publicación: 5 de noviembre de 2025 (arXiv v2)
• Enlace del Artículo: https://arxiv.org/abs/2510.22663

Resumen

Este artículo investiga fenómenos de bifurcación de soluciones retorcidas en el límite continuo (LC) del modelo de Kuramoto de osciladores homogéneos definidos en gráficos de vecinos más cercanos. El estudio abarca gráficos densos deterministas, gráficos densos aleatorios y gráficos dispersos aleatorios, considerando los efectos del desfase de fase (phase-lag). Los autores aplican técnicas estándar de la teoría de sistemas dinámicos, específicamente la reducción de variedad central, para demostrar que el límite continuo experimenta bifurcación en valores de parámetros específicos: familias uniparamétricas de soluciones retorcidas se vuelven inestables, mientras que simultáneamente generan familias biparamétricas estables o inestables de soluciones retorcidas moduladas, cuya oscilación depende de la presencia del desfase de fase. Los resultados teóricos se verifican mediante simulaciones numéricas en gráficos densos deterministas, gráficos densos aleatorios y gráficos dispersos aleatorios.

Antecedentes de Investigación y Motivación

Definición del Problema

El modelo de Kuramoto es uno de los modelos matemáticos más representativos para estudiar fenómenos de sincronización en redes de osciladores acoplados. Este artículo se enfoca en el modelo de Kuramoto definido en el gráfico $G_n$:

$\frac{d}{dt}u_k^n(t) = \omega + \frac{1}{n^{\alpha_n}}\sum_{j=1}^n w_{kj}^n \sin(u_j^n(t) - u_k^n(t) + \sigma), \quad k \in [n]$

donde $u_k^n$ es la fase del nodo $k$, $\omega$ es la frecuencia natural, y $\sigma \in (-\frac{\pi}{2}, \frac{\pi}{2})$ es el parámetro de desfase de fase.

Importancia de la Investigación

1. Significado Teórico: Los estados retorcidos (twisted states) son una clase importante de soluciones estacionarias en el modelo de Kuramoto, con valor significativo para comprender la estructura del espacio de fases y patrones espaciales complejos (como estados quimera)
2. Valor Aplicado: Las redes de osciladores acoplados tienen aplicaciones amplias en física, química, biología, ciencias sociales e ingeniería
3. Desafío Matemático: El comportamiento de bifurcación de soluciones retorcidas en estructuras de gráficos de vecinos más cercanos aún no ha sido estudiado sistemáticamente

Limitaciones de Métodos Existentes

1. La investigación anterior se ha enfocado principalmente en análisis de estabilidad de soluciones retorcidas en el caso $\sigma=0$
2. Falta investigación sistemática para el caso $\sigma \neq 0$ y fenómenos de bifurcación de soluciones retorcidas
3. Garantías teóricas insuficientes para gráficos aleatorios

Motivación de la Investigación

Este artículo tiene como objetivo llenar estas brechas mediante análisis matemático riguroso para revelar los mecanismos de bifurcación de soluciones retorcidas en el modelo de Kuramoto en gráficos de vecinos más cercanos, estableciendo correspondencias precisas entre el modelo discreto y el límite continuo.

Contribuciones Principales

1. Marco de Análisis Teórico: Establece un marco teórico completo para la bifurcación de soluciones retorcidas en el límite continuo del modelo de Kuramoto en gráficos de vecinos más cercanos, utilizando técnicas rigurosas de reducción de variedad central
2. Teoremas de Bifurcación:
   • Demuestra que ocurre bifurcación tipo pitchfork cuando $\sigma=0$, generando familias biparamétricas estables o inestables de soluciones retorcidas moduladas estacionarias
   • Demuestra que ocurre bifurcación de Hopf cuando $\sigma \neq 0$, generando familias biparamétricas estables o inestables de soluciones retorcidas moduladas oscilantes
3. Tratamiento Unificado: Maneja sistemáticamente tres casos: gráficos densos deterministas, gráficos densos aleatorios y gráficos dispersos aleatorios
4. Verificación Numérica: Verifica predicciones teóricas mediante simulaciones numéricas a gran escala, demostrando buena consistencia entre el modelo discreto y el límite continuo
5. Fórmulas Explícitas: Proporciona expresiones explícitas para puntos de bifurcación, soluciones de bifurcación y criterios de estabilidad

Explicación Detallada de Métodos

Definición de la Tarea

El objetivo de investigación es analizar el comportamiento de bifurcación de soluciones q-retorcidas $u(t,x) = 2\pi qx + \Omega t + \theta$ en la ecuación de límite continuo:

$\frac{\partial}{\partial t}u(t,x) = \omega + \int_I W(x,y)\sin(u(t,y) - u(t,x) + \sigma)dy, \quad x \in I=[0,1]$

donde la función de gráfico (graphon) es:

$W(x,y) = \begin{cases} p & \text{si } |x-y| \leq \kappa \text{ o } |x-y| \geq 1-\kappa \\ 0 & \text{en otro caso} \end{cases}$

El parámetro $\kappa \in (0, \frac{1}{2})$ actúa como parámetro de bifurcación.

Arquitectura del Modelo

1. Linealización y Problema de Valores Propios

Linealizando la ecuación de límite continuo cerca de la solución retorcida, se obtiene el operador lineal:

$\mathcal{L}\phi(x) = p\int_{x-\kappa}^{x+\kappa}\cos(2\pi q(y-x) + \sigma)(\phi(y) - \phi(x))dy$

Mediante expansión de Fourier se analizan los valores propios:

• $\phi(x) = 1$ corresponde a valor propio cero
• $\phi(x) = \cos 2\pi\ell x, \sin 2\pi\ell x$ corresponden a valores propios:
  • $\sigma=0$: $\lambda = p\chi_1(\kappa; \ell, q)$
  • $\sigma \neq 0$: $\lambda = p\chi_1(\kappa; \ell, q)\cos\sigma \mp ip\chi_2(\kappa; \ell, q)\sin\sigma$

donde la función clave es: $\chi_1(\kappa; \ell, q) = \begin{cases} \kappa + \frac{\sin 4\pi q\kappa}{4\pi q} - \frac{\sin 2\pi q\kappa}{\pi q} & \text{si } \ell=q \\ \frac{\sin 2\pi(\ell-q)\kappa}{2\pi(\ell-q)} + \frac{\sin 2\pi(\ell+q)\kappa}{2\pi(\ell+q)} - \frac{\sin 2\pi q\kappa}{\pi q} & \text{en otro caso} \end{cases}$

2. Reducción de Variedad Central

Cerca del punto de bifurcación $\kappa = \kappa_{1q}$ (satisfaciendo $\chi_1(\kappa_{1q}; 1, q) = 0$), se introduce el parámetro:

$\mu = p\bar{\chi}'_{1q}(\kappa - \kappa_{1q})$

Expandiendo la solución en serie de Fourier: $u(t,x) = 2\pi qx + \Omega t + \xi_0(t) + \sum_{j=1}^\infty(\xi_j(t)\cos 2\pi jx + \eta_j(t)\sin 2\pi jx)$

Mediante teoría de variedad central, el sistema se reduce a una variedad central tridimensional.

3. Análisis del Sistema Reducido

Caso $\sigma=0$: Introduciendo $r = \sqrt{\xi_1^2 + \eta_1^2}$, el sistema reducido es: $\dot{r} = \mu r - p\beta_0 r^3 + O(\sqrt{r^8 + \mu^4})$

donde $\beta_0 = \beta_1 + \frac{p\delta_1\rho_1}{\mu_2}$ es una constante determinada por parámetros del sistema.

Caso $\sigma \neq 0$: En coordenadas polares $(\xi_1, \eta_1) = (r\cos\psi, r\sin\psi)$: $\dot{r} = \mu r - p\beta_\sigma r^3 + O(\sqrt{r^8 + \mu^4})$$\dot{\psi} = \nu_1 + O(\sqrt{r^2 + \mu^2})$

donde $\beta_\sigma$ depende de una expresión compleja que involucra $\sigma$.

Puntos de Innovación Técnica

1. Marco de Análisis Unificado: Primera demostración matemática completa de bifurcación de soluciones retorcidas en gráficos de vecinos más cercanos, cubriendo casos con y sin desfase de fase
2. Criterios de Bifurcación Precisos: Mediante análisis detallado de valores propios, proporciona ubicación exacta de puntos de bifurcación y criterios de discriminación de tipo de bifurcación
3. Correspondencia Discreto-Continua: Utilizando teoría de convergencia existente (Teoremas 2.2-2.7), establece correspondencia rigurosa entre comportamientos de bifurcación del modelo de Kuramoto discreto y límite continuo
4. Cálculo Explícito: Proporciona fórmulas de cálculo explícitas para todas las cantidades clave (puntos de bifurcación, soluciones de bifurcación, coeficientes de estabilidad)

Configuración Experimental

Configuración de Simulación Numérica

Este artículo realiza simulaciones numéricas del modelo de Kuramoto discreto, considerando tres tipos de gráficos:

1. Gráfico Denso Determinista (caso i): $p=1$, pesos $w_{kj}=1$ si y solo si se satisface la condición de vecino más cercano
2. Gráfico Denso Aleatorio No Dirigido (caso ii): Cuando se satisface la condición de vecino más cercano, se establece $w_{kj}=1$ con probabilidad $p$
3. Gráfico Disperso Aleatorio No Dirigido (caso iii): Cuando se satisface la condición de vecino más cercano, se establece $w_{kj}=1$ con probabilidad $n^{-\gamma}p$ ($\gamma=0.3$)

Configuración de Parámetros

• Escala del Sistema: Principalmente $n=1000$, con algunos casos usando $n=2000$
• Número de Torsión: $q \in \{1, 2, 3, 4\}$
• Desfase de Fase: $\sigma = 0$ o $\sigma = \pi/3$
• Frecuencia Natural: Establecida según la fórmula $\omega = -\frac{p\sin 2\pi q\kappa\sin\sigma}{\pi q}$ tal que $\Omega=0$
• Condiciones Iniciales: Seleccionadas aleatoriamente uniformemente cerca del estado retorcido $u_k^n = 2\pi qk/n$, en el rango $[-10^{-2}, 10^{-2}]$
• Solucionador Numérico: DOP853 (método Runge-Kutta de 8º orden)

Método de Análisis

1. Evolución Temporal: Registro de comportamiento a largo plazo ($t=1000$ o $t=2000$)
2. Análisis de Estado Estacionario: Estimación mediante mínimos cuadrados del estado retorcido más probable
3. Descomposición de Fourier: Extracción de amplitud y fase de términos modulados

Resultados Experimentales

Resultados Principales

1. Gráfico Denso Determinista (caso i)

Caso $\sigma=0$ (Figura 5):

• Para $q \in \{1,2,3,4\}$, cuando $\kappa$ es menor que el punto de bifurcación predicho teóricamente, el sistema permanece cerca del estado retorcido
• La Figura 7 muestra que los resultados numéricos coinciden casi perfectamente con el estado retorcido predicho teóricamente $u(t,x) = 2\pi qx + \bar{u}^n(t) - \pi q$

Caso $\sigma=\pi/3$ (Figura 6):

• $q=1, \kappa=0.31$: El estado retorcido permanece estable, con rotación lenta observada (efecto de tamaño finito)
• $q=2, \kappa=0.16$: Estado retorcido estable
• $q=2, \kappa=0.166$: Superando el punto de bifurcación, se observa solución retorcida modulada oscilante

Verificación de Soluciones de Bifurcación (Figuras 9-10):

• Figura 9: La desviación del estado retorcido coincide perfectamente con el término dominante $r(t)\sin(2\pi x + \psi(t))$ de la fórmula teórica (4.13)
• Figura 10(b): En $\kappa=0.166$ se observan oscilaciones pequeñas sostenidas, verificando la predicción de bifurcación de Hopf
• Figura 10(a): En $\kappa=0.16$ no hay oscilación, consistente con estabilidad antes de bifurcación

Rango de Parámetros Más Grande (Figura 11): En $\kappa=0.33$ (muy por encima del punto de bifurcación), se encuentran diferentes soluciones retorcidas moduladas, cuyo modo de Fourier de primer orden coincide bien con resultados numéricos.

2. Gráfico Denso Aleatorio (caso ii) y Gráfico Disperso Aleatorio (caso iii)

Evolución Temporal (Figuras 12-13):

• En ambos casos de gráficos aleatorios, el sistema mantiene el estado retorcido por debajo del punto de bifurcación
• Existen pequeñas fluctuaciones causadas por aleatoriedad
• Se observa fenómeno de rotación lenta cuando $\sigma=\pi/3$

Análisis de Estado Estacionario (Figuras 14-15):

• Los resultados numéricos son altamente consistentes con el estado retorcido predicho teóricamente
• La desviación causada por aleatoriedad es pequeña
• La Figura 4 muestra mapas de píxeles de la matriz de pesos para $n=1000, p=0.5$ (denso) y $n=2000, \gamma=0.3$ (disperso)

Resultados Cuantitativos

Constantes Clave (Tabla 1, $q \in \{1,2,3,4\}$):

Para $\sigma=0$, dado que $\bar{\chi}'_{1q}\beta_0 < 0$ (para $q \in \{1,2,3,4\}$), según el Teorema 4.2, la familia biparamétrica producida por bifurcación es inestable.

Caso $\sigma \neq 0$ (Figura 3):

• La dependencia de $\beta_\sigma$ respecto a $\sigma$ muestra que para $q=2$, en ciertos rangos de $\sigma$ se tiene $\beta_\sigma > 0$, produciendo soluciones estables
• Para $q \in \{1,3,4\}$, en todo el rango $\sigma \in (-\frac{\pi}{2}, \frac{\pi}{2})$ se tiene $\beta_\sigma < 0$, las soluciones de bifurcación son inestables

Hallazgos Experimentales

1. Consistencia Teoría-Numérica: Los resultados numéricos del modelo discreto son altamente consistentes con predicciones teóricas del límite continuo, incluso a escala moderada $n=1000$
2. Verificación de Tipo de Bifurcación:
   • $\sigma=0$: Se observan características de bifurcación tipo pitchfork
   • $\sigma \neq 0$: Se observan características de bifurcación de Hopf, con oscilaciones sostenidas
3. Robustez de Gráficos Aleatorios: El comportamiento de gráficos densos aleatorios y dispersos es similar al de gráficos deterministas, las predicciones teóricas siguen siendo válidas
4. Efectos de Tamaño Finito: Se observa rotación lenta cuando $\sigma \neq 0$, este es un efecto de tamaño finito, pero no afecta las características principales del comportamiento de bifurcación
5. Multiestabilidad: A valores más grandes de $\kappa$ se encuentran diferentes soluciones retorcidas moduladas, sugiriendo que el sistema puede poseer multiestabilidad rica

Trabajo Relacionado

Investigación del Modelo de Kuramoto

1. Trabajos Clásicos:
   • Kuramoto (1975, 1984): Propone el modelo
   • Strogatz et al. (2000): Investigación de sincronización
   • Acebrón et al. (2005): Revisión
2. Teoría de Límite Continuo:
   • Ermentrout (1985): Límite continuo temprano
   • Medvedev (2014-2019): Teoría de límite continuo para ecuaciones de calor no lineales en gráficos
   • Ihara & Yagasaki (2023): Redes de osciladores acoplados dependientes de multigráficos
3. Investigación de Estados Retorcidos:
   • Girnyk et al. (2012): Multiestabilidad de estados retorcidos en acoplamiento no local
   • Medvedev & Wright (2017): Estabilidad de estados retorcidos en modelo continuo de Kuramoto
   • Medvedev (2014): Osciladores de Kuramoto en redes de mundo pequeño

Teoría de Bifurcación

1. Reducción de Variedad Central: Guckenheimer & Holmes (1983), Haragus & Iooss (2011)
2. Sistemas de Dimensión Infinita: Aplicación a análisis de bifurcación de ecuaciones diferenciales parciales y ecuaciones diferenciales funcionales

Innovación de Este Artículo

1. Primera Investigación Sistemática: Comportamiento de bifurcación de soluciones retorcidas en gráficos de vecinos más cercanos
2. Efecto de Desfase de Fase: Primera vez analizando bifurcación cuando $\sigma \neq 0$
3. Marco Unificado: Cubriendo tres tipos de gráficos: determinista, denso aleatorio, disperso aleatorio
4. Demostración Rigurosa: Proporciona demostración matemática completa, no solo observación numérica

Conclusiones y Discusión

Conclusiones Principales

1. Teoremas de Bifurcación:
   • Demuestra que soluciones q-retorcidas en límite continuo del modelo de Kuramoto en gráficos de vecinos más cercanos experimentan bifurcación en $\kappa=\kappa_{1q}$
   • $\sigma=0$: Bifurcación tipo pitchfork, produciendo soluciones moduladas estacionarias
   • $\sigma \neq 0$: Bifurcación de Hopf, produciendo soluciones moduladas oscilantes
2. Criterios de Estabilidad:
   • La estabilidad de soluciones de bifurcación está determinada por el signo de $\bar{\chi}'_{1q}\beta_0$ (o $\bar{\chi}'_{1q}\beta_\sigma$)
   • Para $q \in \{1,2,3,4\}$, cuando $\sigma=0$ las soluciones de bifurcación son inestables
   • Cuando $\sigma \neq 0$, $q=2$ puede producir soluciones de bifurcación estables
3. Correspondencia Discreto-Continua: La verificación numérica muestra que el comportamiento del modelo de Kuramoto discreto es altamente consistente con predicciones teóricas del límite continuo

Limitaciones

1. Rango de Parámetros: El análisis principal se enfoca en el caso $\ell=1$, aunque bifurcaciones con $\ell \geq 2$ también existen, las soluciones producidas son necesariamente inestables
2. Restricción de Estructura de Gráfico: Solo considera gráficos de vecinos más cercanos, estructuras de gráficos más generales (como redes de mundo pequeño) requieren investigación adicional
3. Términos No Lineales: La reducción de variedad central solo retiene términos hasta tercer orden, efectos de términos de orden superior no se discuten en detalle
4. Escala Numérica: Simulaciones numéricas principalmente realizadas a escala $n=1000$, verificación a escala más grande es limitada
5. Soluciones Oscilantes en Gráficos Aleatorios: No se observan claramente soluciones retorcidas moduladas oscilantes en gráficos aleatorios, puede requerir técnicas numéricas más refinadas

Direcciones Futuras

1. Control por Retroalimentación: Los autores reportarán resultados relacionados de control por retroalimentación en 32
2. Estructuras de Gráfico Más Generales: Extensión a redes de mundo pequeño, redes libres de escala, etc.
3. Osciladores Heterogéneos: Consideración de efectos de distribución de frecuencias naturales
4. Bifurcación Global: Investigación de estructura de bifurcación en rango de parámetros más grande
5. Estados Quimera: Exploración de relación entre soluciones retorcidas y estados quimera

Evaluación Profunda

Fortalezas

1. Rigor Matemático:
   • Proporciona demostración matemática completa, desde análisis de valores propios hasta reducción de variedad central
   • Todas las cantidades clave tienen expresiones explícitas
   • Utiliza teoría de convergencia existente para establecer correspondencia discreto-continua
2. Sistematicidad de Métodos:
   • Manejo unificado de tres tipos de gráficos
   • Consideración simultánea de casos con y sin desfase de fase
   • Proceso completo desde análisis de estabilidad lineal hasta bifurcación no lineal
3. Combinación Teoría-Numérica:
   • Predicciones teóricas altamente consistentes con resultados numéricos
   • Simulaciones numéricas cubren múltiples combinaciones de parámetros
   • Análisis de visualización detallado
4. Contribución Técnica:
   • Análisis detallado de patrones de Fourier
   • Implementación específica de reducción de variedad central
   • Construcción explícita de soluciones de bifurcación
5. Claridad de Escritura:
   • Estructura razonable, lógica clara
   • Apéndice proporciona derivaciones detalladas
   • Figuras abundantes, fácil de entender

Insuficiencias

1. Escenarios de Aplicación:
   • Principalmente análisis teórico, discusión limitada de contexto de aplicación práctica
   • Significado físico y valor aplicado de soluciones retorcidas pueden elaborarse más
2. Complejidad Computacional:
   • El cálculo de coeficientes de bifurcación involucra integrales complejas, puede ser difícil para funciones de gráfico generales
   • Detalles de implementación numérica (como criterios de convergencia) no se explican en detalle
3. Análisis de Gráficos Aleatorios:
   • El análisis teórico de gráficos aleatorios depende principalmente de resultados existentes
   • Efectos de aleatoriedad en comportamiento de bifurcación no se discuten profundamente
4. Comportamiento Global:
   • Se enfoca principalmente en bifurcación local
   • Comportamiento lejos del punto de bifurcación y dominio de atracción global no se estudian sistemáticamente
5. Análisis Comparativo:
   • Comparación insuficiente con gráficos completos u otras estructuras de gráficos
   • Particularidad de gráficos de vecinos más cercanos no se enfatiza suficientemente

Influencia

1. Contribución Teórica:
   • Llena vacío en investigación de bifurcación de soluciones retorcidas en gráficos de vecinos más cercanos
   • Proporciona marco de análisis generalizable a otras estructuras de gráficos
   • Enriquece teoría de bifurcación del modelo de Kuramoto
2. Valor Metodológico:
   • Demuestra aplicación efectiva de reducción de variedad central en dinámica de redes
   • Establece puente riguroso entre modelos discretos y límite continuo
   • Proporciona referencia para análisis de otras redes de osciladores acoplados
3. Valor Práctico:
   • Proporciona orientación teórica para diseño de sistemas como redes neuronales y redes eléctricas
   • Predicción de ubicación de punto de bifurcación ayuda a controlar comportamiento del sistema
   • Criterios de estabilidad pueden usarse para selección de parámetros
4. Reproducibilidad:
   • Derivación teórica detallada, fácil de verificar
   • Parámetros de experimento numérico claros
   • Aunque código no es público, detalles de implementación son claros

Escenarios Aplicables

1. Investigación Teórica:
   • Análisis de bifurcación de sistemas dinámicos en redes
   • Modelado matemático de fenómenos de sincronización
   • Aplicación de teoría de límite de gráficos
2. Sistemas Prácticos:
   • Patrones de ritmo en redes neuronales
   • Sincronización de fase en sistemas eléctricos
   • Comportamiento colectivo de osciladores biológicos
   • Dinámica de opinión en redes sociales
3. Enseñanza:
   • Casos de teoría de bifurcación de sistemas dinámicos
   • Ejemplos de aplicación de técnica de reducción de variedad central
   • Intersección de ciencia de redes y sistemas dinámicos

Referencias

Este artículo cita 32 referencias, principalmente incluyendo:

1. Modelo de Kuramoto: Kuramoto (1975, 1984), Acebrón et al. (2005)
2. Límite Continuo: Serie de trabajos de Medvedev (2014-2019), Ihara & Yagasaki (2023)
3. Teoría de Bifurcación: Guckenheimer & Holmes (1983), Haragus & Iooss (2011)
4. Estados Retorcidos: Girnyk et al. (2012), Medvedev & Wright (2017)
5. Trabajos Anteriores del Autor: Yagasaki (2025a, 2025b)

Evaluación General: Este es un artículo de alta calidad que realiza investigación sistemática y rigurosa de bifurcación de soluciones retorcidas del modelo de Kuramoto en gráficos de vecinos más cercanos. Las principales fortalezas radican en la completitud del análisis matemático, la integración cercana de teoría y numérica, y el tratamiento unificado de gráficos deterministas y aleatorios. El artículo no solo demuestra la existencia de bifurcación, sino que también proporciona expresiones explícitas para puntos de bifurcación, soluciones de bifurcación y criterios de estabilidad. Las simulaciones numéricas apoyan fuertemente las predicciones teóricas. Este trabajo proporciona base teórica importante para comprender comportamiento colectivo de osciladores en redes, con contribuciones significativas tanto a sistemas dinámicos como a ciencia de redes.