The Ott–Antonsen ansatz shows that, for certain classes of distribution of the natural frequencies in systems of N globally coupled Kuramoto oscillators, the dynamics of the order parameter, in the limit N, evolves, under suitable initial conditions, in a manifold of low dimension. This is not possible when the frequency distribution, continued in the complex plane, has an essential singularity at infinity; this is the case, for example, of a Gaussian distribution. In this work, we propose a simple approximation scheme that allows one to extend also to this case the representation of the dynamics of the order parameter in a low dimensional manifold. Using the Gaussian frequency distribution as a working example, we compare the dynamical evolution of the order parameter of the system of oscillators, obtained by the numerical integration of the N equations of motion, with the analogous dynamics in the low dimensional manifold obtained with the application of the approximation scheme. The results confirm the validity of the approximation. The method could be employed for general frequency distributions, allowing the determination of the corresponding phase diagram of the oscillator system.

1.
E.
Ott
and
T. M.
Antonsen
, “
Low dimensional behavior of large systems of globally coupled oscillators
,”
Chaos
18
,
037113
(
2008
).
2.
A.
Pykovsky
,
M.
Rosenblum
, and
J.
Kurths
,
Synchronization: A Universal Concept in Nonlinear Sciences
(
Cambridge University Press
,
Cambridge
,
2001
).
3.
K.
Wiesenfled
,
P.
Colet
, and
S. H.
Strogatz
, “
Frequency locking in Josephson arrays: Connection with the Kuramoto model
,”
Phys. Rev. E
57
,
1563
(
1998
).
4.
J.
Buck
, “
Synchronous rhythmic flashing of fireflies. II
,”
Q. Rev. Biol.
63
,
265
(
1988
).
5.
A. T.
Winfree
,
The Geometry of Biological Time
(
Springer
,
New York
,
1980
).
6.
G.
Filatrella
,
A. H.
Nielsen
, and
N. F.
Pedersen
, “
Analysis of a power grid using a Kuramoto-like model
,”
Eur. Phys. J. B
61
,
485
(
2008
).
7.
M.
Bier
,
B. M.
Bakker
, and
H. V.
Westerhoff
, “
How yeast cells synchronize their glycolytic oscillations: A perturbation analytic treatment
,”
Biophys. J.
78
,
1087
(
2000
).
8.
S. H.
Strogatz
,
Sync: The Emerging Science of Spontaneous Order
(
Hyperion
,
New York
,
2003
).
9.
Y.
Kuramoto
, in International Symposium on Mathematical Problems in Theoretical Physics, Lecture Notes in Physics Vol. 39, edited by H. Arakai (Springer, New York, 1975).
10.
A.
Campa
,
T.
Dauxois
, and
S.
Ruffo
, “
Statistical mechanics and dynamics of solvable models with long-range interactions
,”
Phys. Rep.
480
,
57
(
2009
).
11.
If the parameters characterizing the distribution g(ω) are s, we expect that the dynamics depends on the ratio between these parameters and the coupling constant K, so that the phase diagram has dimension s.
12.
S. H.
Strogatz
, “
From Kuramoto to Crawford: Exploring the onset of synchronization in populations of coupled oscillators
,”
Physica D
143
,
1
(
2000
).
13.
S.
Gupta
,
A.
Campa
, and
S.
Ruffo
,
Statistical Physics of Synchronization
(
Springer
,
Cham
,
2018
).
14.
More rigorously, the continuity equation should be written for the N-particle distribution function, but the factorization approximation usually employed for systems with mean field interactions, when N, allows one to get an equation for the one-particle distribution function. For the passage from the N-particle distribution function to the one-particle distribution function, even in the more general context in which the system is subject to noise and the description has to be made with a Fokker–Planck equation, see, e.g., Ref. 13.
15.
It is not necessary to make also for g(ω) the assumption that |g(ω)|0 for Im(ω) since this is a consequence of the fact that it has no singularities in the lower half plane other than some poles at finite values; together with the physical requirement that g(ω) vanishes for ω real tending to ±, this implies that it has to vanish also for Im(ω).
16.
A linearized Vlasov equation would have, as a last term, something proportional to g(ω)dωf(ω,t), i.e., with the known function g(ω) outside the integral.
17.
We could not avoid to use the universally adopted symbol ω for the proper frequencies of the oscillators; thus, we had to use another symbol, ν, for the variable conjugated to the time t. This should not generate confusion.
18.
For the case M=1, the solution of the single equation, to which the system (20) reduces to, is known analytically also for the full nonlinear case,1 but here we consider only the linearized equation, since we want to compare with the result of the linear Eq. (22).
19.
E. A.
Martens
,
E.
Barreto
,
S. H.
Strogatz
,
E.
Ott
,
P.
So
, and
T. M.
Antonsen
, “
Exact results for the Kuramoto model with a bimodal frequency distribution
,”
Phys. Rev. E
79
,
026204
(
2009
).
20.
J. R.
Engelbrecht
and
R.
Mirollo
, “
Is the Ott-Antonsen manifold attracting?
,”
Phys. Rev. Res.
2
,
023057
(
2020
).
You do not currently have access to this content.