Synchronization: The Kuramoto Model

Synchronization is one of nature's most universal phenomena: fireflies flash in unison, neurons oscillate together during memory consolidation, Huygens' clocks entrained on the same beam. Kuramoto's model gives an exact analytical description of this transition.

• Neuroscience: synchronization of hippocampal oscillators at gamma frequency (40 Hz) underlies memory consolidation; pathological oversynchronization is the mechanism of epileptic seizures
• Power grids: generator synchronization at 50/60 Hz is safety-critical; desynchronization after a cascade failure follows the reverse Kuramoto transition
• Cardiology: the sinoatrial node contains ~10000 pacemaker cells maintaining synchronous rhythm through a noisy Kuramoto model
• Laser physics: mode locking in a laser array is a direct Kuramoto application, used for coherent power combining in fiber laser arrays

Предварительные знания

• Phase oscillators
• Bifurcation theory
• Phase transitions (basics)

• Previous lesson: ds-27

Kuramoto model: coupled oscillators

Christiaan Huygens in 1665 first described synchronization: two pendulum clocks on the same beam synchronized within 30 minutes. In 1975, Yoshiki Kuramoto proposed an exactly solvable model of N coupled oscillators. Neural networks in the brain, pacemaker cells (N ~ 10000), and power grids show the same second-order phase transition.

Order parameter and mean field

Analogy with ferromagnetism: r is the analog of magnetization M, K_c is the analog of Curie temperature T_c. The second-order transition with r ~ (K - K_c)^{1/2} matches the mean-field exponent M ~ (T_c - T)^{1/2} in Landau theory.

Synchronization as a phase transition. Picture N metronomes on a wobbly table. At K=0 (rigid table) each ticks in its own rhythm. As K increases (table gets softer) metronomes start sensing each other through vibration. Past K_c a critical mass locks together, creating a collective field that pulls in more and more stragglers.

Phase transition and critical coupling

The Kuramoto transition is an emergent phenomenon: no single oscillator is synchronized in isolation, but the collective system undergoes an ordered phase transition. This is the mechanism behind gamma-band neural oscillations and power-grid frequency coordination.

Extensions of the Kuramoto model: networks with arbitrary topology (K replaced by coupling matrix A_ij/k_i), time delay (theta'_i = omega_i + K/N*sum sin(theta_j(t-tau) - theta_i)), adaptive plasticity (A_ij evolves with synchrony). On scale-free networks K_c -> 0 because <k^2>/<k> diverges.

Connections to other areas

The Kuramoto model bridges dynamical systems theory, statistical physics, and network science.

• Complex Networks and Graph Dynamics — Kuramoto on graphs is a central model of network synchronization
• Synchronization and Coupled Oscillators — Introductory chapter on synchronization in coupled-oscillator systems
• Bifurcation Theory — Onset of synchronization in Kuramoto is an order-disorder Hopf bifurcation

Итоги

• Kuramoto model: theta'_i = omega_i + K/N*sum sin(theta_j - theta_i); order parameter r = |sum exp(i*theta_j)|/N
• Phase transition: r = 0 for K < K_c (incoherence), r > 0 for K > K_c - a second-order transition
• Critical coupling: K_c = 2/(pi*g(0)), inversely proportional to frequency density at zero
• Self-consistency: oscillators split into phase-locked (|omega| < K*r) and drifting (|omega| >= K*r)
• Lorentzian distribution: exact solution r = sqrt(1 - K_c/K) for K > K_c = 2*gamma
• Ferromagnet analogy: r = magnetization, K_c = Curie temperature, exponent 1/2 - mean-field universality class

How does critical coupling K_c depend on the natural frequency distribution g(omega)?

K_c = 2/(pi*g(0)) depends only on the density of oscillators with zero detuning. Narrow distribution (large g(0)) means more oscillators near zero frequency, easier to synchronize. For Lorentzian g(0) = 1/(pi*gamma), giving K_c = 2*gamma.