🔄 Kuramoto Synchronization

The Kuramoto model simulates a large population of coupled oscillators, each with a slightly different natural frequency, that tend to synchronise when mutual coupling is strong enough. It captures the phase transition from disordered (asynchronous) to ordered (synchronous) collective behaviour, serving as a paradigmatic model for synchronisation in complex systems.

The order parameter r measures the degree of synchronisation. It is the magnitude of the average unit phasor across all oscillators: r = |mean(exp(i*theta_j))|. When r = 0, oscillators are uniformly spread around the circle (incoherent). When r = 1, all oscillators have the same phase (fully synchronised). The transition from r ≈ 0 to r > 0 occurs sharply above the critical coupling Kc.

Synchronisation driven by Kuramoto-like coupling occurs in: firefly flash synchronisation (Southeast Asian Pteroptyx fireflies produce perfectly synchronous flashing), cardiac pacemaker cells in the sinoatrial node (millions of cells beating together), neural oscillations in the brain (gamma oscillations coordinate sensory processing), and menstrual cycle synchronisation in social groups (though evidence for the last is debated).

Power generators connected to the electricity grid must all spin at the same frequency (50 Hz in Europe, 60 Hz in North America) to avoid destructive interference. Generators are coupled through shared transmission lines, and the dynamics of their phase synchronisation is described by a variant of the Kuramoto model. Loss of synchronisation — when a generator's phase slips — causes protective relays to disconnect it, potentially triggering cascading grid failures.

The critical coupling Kc is the minimum coupling strength at which the system can spontaneously synchronise. For N oscillators with Lorentzian-distributed natural frequencies of half-width gamma, Kc = 2*gamma. Below Kc, the incoherent state r = 0 is stable for all initial conditions. Above Kc, synchronisation grows as r ~ sqrt(1 - Kc/K) just above the transition, analogous to an order parameter near a second-order phase transition.