HIGHER-ORDER KURAMOTO

Order param R0.00
Mean freq ⟨ω⟩0.00
dθᵢ/dt = ωᵢ + (K₁/N)Σⱼsin(θⱼ−θᵢ)
+ (K₂/N²)Σⱼₖsin(θⱼ+θₖ−2θᵢ)

R·eⁱψ = (1/N)Σₙ eⁱθₙ
Daido order param: Zₚ=⟨eⁱpθ⟩
Higher-order interactions enable explosive sync & chimera states