WOLFRAM|DEMONSTRATIONS PROJECT

Kuramoto Model for Phase Locking of Coupled Oscillators

​
coupling parameter
K = 0.05
K =
K
C
= 0.1
K = 0.25
K = 0.5
running time
0
This Demonstration shows the collective synchronization of
N
coupled oscillators according to the Kuramoto model [1], which is a modified version of the Winfree model [2] for the population of biological oscillators.
Many systems in biology, physics and engineering show the phenomenon of collective synchronization, such as heart pacemaker cells, flashing fireflies, crickets chirping in unison, laser arrays and microwave oscillators.
The Kuramoto model is defined through the following set of
N
time-dependent coupled differential equations

θ
i
=
ω
i
+
K
N
N
∑
j=i
sin
θ
j
-
θ
i
,i=1,…,N
where
θ
i
(t)
is the time-dependent phase of the
th
i
oscillator and
ω
i
is its natural frequency chosen from a probability density distribution (in this Demonstration, a Gaussian). The starting values of the
θ
i
are chosen from a uniform density distribution. The parameter
K
is the coupling strength between each pair
(i,j)
of oscillators. While
K
is below a certain critical value
K
C
, the oscillators behave completely incoherently: you can see the points giving the
θ
i
(t)
values spreading at random around the unit circle. When
K
exceeds the critical value, the oscillator population locks in a partially synchronized group running around the unit circle at the same frequency. By further increasing
K
, more and more oscillators are added into the synchronized group locked in phase.
This phase locking can thus be interpreted as a temporal phase transition, where the
K
-dependent order parameter gives the amount of synchronization in the oscillator population as
N
and time tend to infinity.
The degree of synchronization is given by the order parameter
1
N

N
∑
k=1
iθ
k
e

, which tends to 1 as
K>
K
C
for increasing
N
and time, as you can check by running this Demonstration.
To optimize graphical performance, the solutions for the set of
θ
i
(t)
are precomputed in the initialization code for a given set of coupling parameters
K
. In order to further explore the dependence of synchronization onset on model parameters (number of oscillators, coupling parameter, initial conditions), change the corresponding variables in the initialization code.