The Kuramoto model is a mathematical model of synchronisation in a population of coupled oscillators. Each oscillator has a natural frequency drawn from a distribution and a phase that evolves over time; oscillators are coupled so that each one is attracted toward the average phase of the population. When the coupling strength exceeds a critical threshold, the population transitions from incoherent (all oscillating independently) to partially or fully synchronised (all oscillating together), a non-equilibrium phase transition analogous to ferromagnetism.
The model equations are dθᵢ/dt = ωᵢ + (K/N) Σⱼ sin(θⱼ - θᵢ), where θᵢ is the phase of oscillator i, ωᵢ is its natural frequency, K is the coupling constant, and N is the number of oscillators. The order parameter r = (1/N)|Σⱼ exp(iθⱼ)| measures global coherence (r = 0 is incoherent, r = 1 is fully synchronised). For a Lorentzian frequency distribution, Kuramoto derived an exact solution for the critical coupling Kc = 2γ, where γ is the half-width of the distribution.
The Kuramoto model captures synchronisation phenomena across biology, physics, and engineering. Fireflies synchronise their flashing using similar mutual coupling. Cardiac pacemaker cells synchronise to produce coherent heartbeats. Superconducting Josephson junctions (used in quantum computing hardware and sensitive magnetometers) synchronise their oscillations when coupled. Power grid generators must synchronise their AC frequencies (50 or 60 Hz) to maintain grid stability — desynchronisation causes blackouts.
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.