Kuramoto model

Kuramoto model

The Kuramoto model, first proposed by Yoshiki Kuramoto (蔵本 由紀 "Kuramoto Yoshiki"), is a mathematical model used to describe synchronization. More specifically, it is a model for the behavior of a large set of coupled oscillators. Its formulation was motivated by the behavior of systems of chemical and biological oscillators, and it has found widespread applications.Fact|date=March 2008

The model makes several assumptions, including that there is weak coupling, identical or nearly identical oscillators, and that interactions depend sinusoidally on the phase difference between the two objects.


In the most popular version of the Kuramoto model, each of the oscillators is considered to have its own intrinsic natural frequency omega_i, and each is coupled equally to all other oscillators. Surprisingly, this fully nonlinear model can be solved exactly, in the infinite-"N" limit, with a clever transformation and the application of self-consistency arguments.

The most popular form of the model has the following governing equations:: frac{partial heta_i}{partial t} = omega_i + frac{K}{N} sum_{j=1}^{N} sin( heta_j - heta_i), qquad i = 1 ldots N,where the system is composed of "N" limit-cycle oscillators.

Noise can be added to the system. In that case, the original equation is altered to::frac{partial heta_i}{partial t} = omega_{i}+zeta_{i}+dfrac{K}{N}sum_{j=1}^Nsin( heta_{j}- heta_{i}),where zeta_{i} is the fluctuation and a function of time. If we consider the noise to be white noise, then:langlezeta_{i}(t) angle=0 ,:langlezeta_{i}(t)zeta_{j}(t') angle=2Ddelta_{ij}delta(t-t')
with D denotes the strength of noise.


The transformation that allows this model to be solved exactly (at least in the "N" → ∞ limit) is as follows.Define the "order" parameters "r" and "ψ" as:re^{i psi} = frac{1}{N} sum_{j=1}^{N} e^{i heta_j} .Here "r" represents the phase-coherence of the population of oscillators, and "ψ" indicates the average phase. Applying this transformation, the governing equation becomes: frac{partial heta_i}{partial t} = omega_i + K r sin(psi- heta_i) .Thus the oscillators' equations are no longer explicitly coupled; instead the order parameters govern behavior. A further transformation is usually done, to a rotating frame in which the statistical average of phases over all oscillators is zero. That is, psi=0. Finally, the governing equation becomes: frac{partial heta_i}{partial t} = omega_i - K r sin( heta_i) .

Large "N" limit

Now consider the case as "N" tends to infinity. Take the distribution of intrinsic natural frequencies as "g"("ω") (assumed normalized). Then assume that the density of oscillators at a given phase "θ", with given natural frequency "ω", at time "t" is ho( heta, omega, t). Normalization requires that : int_{-pi}^{pi} ho( heta, omega, t) , d heta = 1.

The continuity equation for oscillator density will be : frac{partial ho}{partial t} + frac{partial}{partial heta} [ ho v] = 0, where "v" is the drift velocity of the oscillators given by taking the infinite-"N" limit in the transformed governing equation, i.e., : frac{partial ho}{partial t} + frac{partial}{partial heta} [ ho omega + ho K r sin(psi- heta)] = 0.

Finally, we must rewrite the definition of the order parameters for the continuum (infinite "N") limit. heta_i must be replaced by its ensemble average (over all ω) and the sum must be replaced by an integral, to give: r e^{i psi} = int_{-pi}^{pi} e^{i heta} int_{-infty}^{infty} ho( heta, omega, t) g(omega) , d omega , d heta.


The incoherent state with all oscillators drifting randomly corresponds to the solution ho = 1/(2pi). In that case r = 0, and there is no coherence among the oscillators. They are uniformly distributed across all possible phases, and the population is in a statistical steady-state (although individual oscillators continue to change phase in accordance with their intrinsic "ω").

When coupling "K" is sufficiently strong, a fully synchronized solution is possible. In the fully synchronized state, all the oscillators share a common frequency, although their phases are different.

A solution for the case of partial synchronization yields a state in which only some oscillators (those near the ensemble's mean natural frequency) synchronize; other oscillators drift incoherently. Mathematically, the state has: ho = deltaleft( heta - psi - arcsinleft(frac{omega}{K r} ight) ight) for locked oscillators, and: ho = frac{ m{normalization ; constant{(omega - K r sin( heta - psi))}for drifting oscillators. The cutoff occurs when |omega| < K r .


Acebrón, Juan A.; Bonilla, L. L.; Pérez Vicente, Conrad J.; Ritort, Félix; Spigler, Renato. The Kuramoto model: a simple paradigm for synchronization phenomena [http://scala.uc3m.es/publications_MANS/PDF/finalKura.pdf]

Strogatz, S. From Kuramoto to Crawford: exploring the onset of synchronization in populations of coupled oscillators, Physica D 143 (2000) 1–20

Wikimedia Foundation. 2010.

Look at other dictionaries:

  • Ising model — The Ising model, named after the physicist Ernst Ising, is a mathematical model in statistical mechanics. It has since been used to model diverse phenomena in which bits of information, interacting in pairs, produce collectiveeffects.Definition… …   Wikipedia

  • Yoshiki Kuramoto — is a Japanese physicist in the Nonlinear Dynamics group at Kyoto University who formulated the Kuramoto model …   Wikipedia

  • Neural oscillation — is rhythmic or repetitive neural activity in the central nervous system. Neural tissue can generate oscillatory activity in many ways, driven either by mechanisms localized within individual neurons or by interactions between neurons. In… …   Wikipedia

  • Smart grid — Public infrastructure …   Wikipedia

  • Scientific phenomena named after people — This is a list of scientific phenomena and concepts named after people (eponymous phenomena). For other lists of eponyms, see eponym. NOTOC A* Abderhalden ninhydrin reaction Emil Abderhalden * Abney effect, Abney s law of additivity William de… …   Wikipedia

  • Coupled map lattice — A coupled map lattice (CML) is a dynamical system that models the behavior of non linear systems (especially partial differential equations). They are predominantly used to qualitatively study the chaotic dynamics of spatially extended systems.… …   Wikipedia

  • List of mathematics articles (K) — NOTOC K K approximation of k hitting set K ary tree K core K edge connected graph K equivalence K factor error K finite K function K homology K means algorithm K medoids K minimum spanning tree K Poincaré algebra K Poincaré group K set (geometry) …   Wikipedia

  • Kurt Wiesenfeld — is an American physicist working primarily on non linear dynamics. His works primarily concern stochastic resonance, spontaneous synchronization of coupled oscillators, and non linear laser dynamics. Since 1987, he has been professor of physics… …   Wikipedia

  • Synchronization — Synchrony redirects here. For linguistic synchrony, see Synchronic analysis. For the X Files episode, see Synchrony (The X Files). For similarly named concepts, see Synchronicity (disambiguation). Not to be confused with data… …   Wikipedia

  • List of Battle Royale characters — The following is a list of characters that appear in the novel, manga, and film versions of Battle Royale. Contents 1 Primary characters 1.1 Shuya Nanahara 1.2 Noriko Nakagawa 1.3 …   Wikipedia

We are using cookies for the best presentation of our site. Continuing to use this site, you agree with this.