N phase oscillators with natural frequencies ωi drawn from a distribution g(ω). All-to-all coupling:
dθi/dt = ωi + (K/N) Σj sin(θj − θi)
The order parameter r measures coherence:
r eiψ = (1/N) Σj eiθj
r ≈ 0 = incoherent, r ≈ 1 = fully phase-locked.
What I actually simulate: nearest-neighbour coupling on a ring, with additive white noise on each node:
dθi/dt = ωi + K[sin(θi+1 − θi) + sin(θi−1 − θi)] + σ ξi(t)
Where ξi(t) is unit white noise, indices mod N. The ωi are drawn from a Gaussian with mean 0 and spread Δω.
When K is too small to maintain lock, individual oscillators drift relative to the mean field. A phase slip is a 2π jump in the phase difference between an oscillator and the order parameter ψ.
For the homogeneous case (Δω = 0, σ = 0), the critical coupling for ring topology can be shown to be:
Kc = (π/2) Δωmax
where Δωmax is the maximum frequency difference between adjacent nodes. With noise, the boundary smears.
Personal research notes. Expect errors.