Modelo de Kuramoto

El modelo de Kuramoto (o modelo de Kuramoto-Daido) es un modelo matemático utilizado para describir la sincronización propuesto por Yoshiki Kuramoto (蔵本 由紀, Yoshiki Kuramoto?)

[4]​[5]​ Su formulación fue motivada por el comportamiento de sistemas de osciladores químicos y biológicos, y ha encontrado amplias aplicaciones en áreas como la neurociencia.

[6]​[7]​[8]​ Kuramoto se sorprendió bastante cuando vio que el comportamiento de algunos sistemas físicos, concretamente las matrices acopladas de uniones de Josephson, seguían su modelo.

[9]​ El modelo hace varias suposiciones, incluyendo que hay un acoplamiento débil, que los osciladores son idénticos o casi idénticos, y que las interacciones dependen sinusoidalmente de la diferencia de fase entre cada par de objetos.

En la versión más conocida del modelo de Kuramoto, se considera que cada uno de los osciladores tiene su propia frecuencia natural intrínseca

, y cada uno está acoplado por igual a todos los demás osciladores.

Sorprendentemente, este modelo totalmente no lineal puede resolverse exactamente cuando el número de osciladores N tiende a infinito, N→ ∞;[6]​ alternativamente, utilizando argumentos de autoconsistencia se pueden obtener soluciones de estado estacionario de r, el parámetro que describe el grado de sincronización del sistema.

[10]​ La forma más popular del modelo tiene las siguientes ecuaciones de gobierno:[11]​

, donde el sistema se compone de N osciladores de ciclo límite, con fases

y constante de acoplamiento K. Se puede añadir ruido al sistema mediante la introducción de un parámetro,

, que fluctúa en función del tiempo, generando la nueva ecuación:[11]​

Para resolver este modelo, definimos dos nuevos "parámetros de orden": r, que representa la coherencia de fase de la población de osciladores y ψ, que indica la fase media.

Así:[12]​ De esta manera, las ecuaciones de los osciladores ya no están explícitamente acopladas, sino que son los parámetros de orden los que gobiernan el comportamiento.

Se puede hacer una transformación más, a un marco de rotación en el que la media estadística de las fases sobre todos los osciladores es cero (es decir,

Consideremos ahora el caso cuando N tiende a infinito.

Tomemos la distribución normalizada de las frecuencias naturales intrínsecas como g(ω) y supongamos que la densidad de osciladores en una fase dada θ, con una frecuencia natural dada ω, en el tiempo t sigue la función

En este caso, los osciladores se distribuyen uniformemente a través de todas las fases posibles, y la población se encuentra en un estado estacionario (aunque los osciladores individuales siguen cambiando de fase de acuerdo con su ω intrínseco).

Cuando la fuerza de acoplamiento K es suficientemente fuerte, es posible una solución totalmente sincronizada.

En el estado totalmente sincronizado, todos los osciladores comparten una frecuencia común, aunque sus fases pueden ser diferentes.

Sin embargo, para valores intermedios de K, se produce un estado de sincronización parcial, un estado en el que solo algunos osciladores (los que están cerca de la frecuencia natural media del conjunto) se sincronizan; otros osciladores derivan de forma incoherente.

Matemáticamente, el estado tiene las ecuaciones:[14]​ para los osciladores bloqueados, y:[14]​ para osciladores a la deriva.

Este caso se produce cuando

Además del modelo original, que considera que todos los osciladores interaccionan igualmente entre ellos independientemente de su distancia, se han desarrollado variaciones que tienen en cuenta topologías de red compleja que admiten estados quimera y comportamientos locales diferenciados.

[15]​ En redes bidimensionales de Kuramoto con acoplamiento local difusivo, es habitual encontrar sincronía uniforme, ondas y espirales cuya estabilidad puede determinarse analíticamente utilizando los métodos del análisis de estabilidad de Turing.

[16]​ La topología sobre la que se estudia el modelo Kuramoto puede hacerse adaptativa mediante el uso de un modelo de aptitud, lo que mejora la sincronización y la percolación de forma autoorganizada.

[17]​ Algunos trabajos se han centrado en el modelo Kuramoto en redes en las que la fuerza de interconexión entre dos osciladores cualesquiera es arbitraria.

La dinámica de este modelo es la siguiente: donde

es un número real positivo no nulo si el oscilador

está conectado al oscilador

[18]​ Dado que el modelo de Kuramoto parece desempeñar un papel clave en la evaluación de los fenómenos de sincronización en el cerebro,[19]​ estas variaciones podrían allanar el camino para una comprensión más profunda de los fenómenos de sincronización neuronal.