|
|
| (One intermediate revision by one other user not shown) |
| Line 1: |
Line 1: |
| 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 process|biological]] oscillators, and it has found widespread applications such as in [[Neural oscillation#Mathematical description|neuroscience]]. | | The title of the author is Numbers. He utilized to be unemployed but now he is a pc operator but his marketing never comes. South Dakota is her beginning location but she needs to transfer because of her family. Doing ceramics is what her family members and her appreciate.<br><br>My web-site :: home std test ([http://raybana.com/chat/pg/groups/174883/solid-advice-for-dealing-with-a-candida/ visit the up coming article]) |
| | |
| The model makes several assumptions, including that there is weak coupling, that the oscillators are identical or nearly identical, and that interactions depend sinusoidally on the phase difference between each pair of objects.
| |
| | |
| == Definition ==
| |
| | |
| [[File:KuramotoModel.ogv|thumb|240px| right|Simulation of Kuramoto model showing neural synchronization and oscillations in the mean field]]
| |
| | |
| In the most popular version of the Kuramoto model, each of the oscillators is considered to have its own intrinsic [[natural frequency]] <math>\omega_i</math>, 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:
| |
| :<math> \frac{d \theta_i}{d t} = \omega_i + \frac{K}{N} \sum_{j=1}^{N} \sin(\theta_j - \theta_i), \qquad i = 1 \ldots N</math>,
| |
| where the system is composed of ''N'' limit-cycle oscillators.
| |
| | |
| <br />Noise can be added to the system. In that case, the original equation is altered to:
| |
| :<math>
| |
| \frac{d \theta_i}{d t} = \omega_{i}+\zeta_{i}+\dfrac{K}{N}\sum_{j=1}^N\sin(\theta_{j}-\theta_{i})
| |
| </math>,
| |
| where <math>\zeta_{i}</math> is the fluctuation and a function of time. If we consider the noise to be white noise, then
| |
| :<math>
| |
| \langle\zeta_{i}(t)\rangle=0
| |
| </math> ,
| |
| :<math>
| |
| \langle\zeta_{i}(t)\zeta_{j}(t')\rangle=2D\delta_{ij}\delta(t-t')
| |
| </math>
| |
| <br />with <math>D</math> denoting the strength of noise. | |
| | |
| == Transformation ==
| |
| 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
| |
| :<math>re^{i \psi} = \frac{1}{N} \sum_{j=1}^{N} e^{i \theta_j} </math>.
| |
| Here ''r'' represents the phase-[[coherence (physics)|coherence]] of the population of oscillators, and ''ψ'' indicates the average phase. Applying this transformation, the governing equation becomes
| |
| :<math> \frac{d \theta_i}{d t} = \omega_i + K r \sin(\psi-\theta_i) </math>.
| |
| 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, <math>\psi=0</math>. Finally, the governing equation becomes
| |
| :<math> \frac{d \theta_i}{d t} = \omega_i - K r \sin(\theta_i) </math>. | |
| | |
| == Large ''N'' limit ==
| |
| Now consider the case as ''N'' tends to infinity. Take the distribution of intrinsic natural frequencies as ''g''(''ω'') (assumed [[normalization constant|normalized]]). Then assume that the density of oscillators at a given phase ''θ'', with given natural frequency ''ω'', at time ''t'' is <math>\rho(\theta, \omega, t)</math>. Normalization requires that
| |
| :<math> \int_{-\pi}^{\pi} \rho(\theta, \omega, t) \, d \theta = 1. </math> | |
| | |
| The [[continuity equation]] for oscillator density will be
| |
| :<math> \frac{\partial \rho}{\partial t} + \frac{\partial}{\partial \theta}[\rho v] = 0, </math> | |
| where ''v'' is the drift velocity of the oscillators given by taking the infinite-''N'' limit in the transformed governing equation, i.e.,
| |
| :<math> \frac{\partial \rho}{\partial t} + \frac{\partial}{\partial \theta}[\rho \omega + \rho K r \sin(\psi-\theta)] = 0. </math>
| |
| | |
| Finally, we must rewrite the definition of the order parameters for the continuum (infinite ''N'') limit. <math>\theta_i</math> must be replaced by its ensemble average (over all ω) and the sum must be replaced by an integral, to give
| |
| :<math>
| |
| r e^{i \psi} = \int_{-\pi}^{\pi} e^{i \theta} \int_{-\infty}^{\infty} \rho(\theta, \omega, t) g(\omega) \, d \omega \, d \theta.
| |
| </math>
| |
| | |
| == Solutions ==
| |
| The [[Coherence (physics)|incoherent]] state with all oscillators drifting randomly corresponds to the solution <math>\rho = 1/(2\pi)</math>. In that case <math>r = 0</math>, 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
| |
| :<math>\rho = \delta\left(\theta - \psi - \arcsin\left(\frac{\omega}{K r}\right)\right)</math>
| |
| for locked oscillators, and | |
| :<math>\rho = \frac{\rm{normalization \; constant}}{(\omega - K r \sin(\theta - \psi))}</math>
| |
| for drifting oscillators. The cutoff occurs when <math>|\omega| < K r </math>.
| |
| | |
| == Variations on the models ==
| |
| There are two types of variations of models based on the original model presented above, one who deals with changing the topological structure of the model, the second are more related to models who are inspired by kuramoto model but dont have the same functional form.
| |
| | |
| === Variations on topology ===
| |
| Beside the original model, which has an all-to-all topology, a sufficiently dense [[Complex Network|complex network]] like topology is amenable to the mean-field treatment used in the solution of the original model (see Transformation and Large N Limit above for more info).
| |
| One also may ask for the behavior of models in which there are intrinsically local , like one dimensional topologies which the chain and the ring are prototipical examples. In such topologies, in which the coupling is not scalable according to <math>\frac{1}{N}</math>, it's not possible to apply the canonical mean-field approach, so one must relies upon case-by-case analysis, making use of simmetries whenever it is possible, which may give basis for abstraction of general principles of solutions.
| |
| | |
| == References ==
| |
| *{{citation
| |
| |last = Kuramoto | first = Yoshiki
| |
| |title = Lecture Notes in Physics, International Symposium on Mathematical Problems in Theoretical Physics
| |
| |volume = 39
| |
| | publisher = Springer-Verlag, New York
| |
| | page = 42
| |
| |year = 1975
| |
| |editor= H. Araki
| |
| }}.
| |
| *{{citation
| |
| |last = Kuramoto | first = Yoshiki
| |
| |title = Chemical Oscillations, Waves, and Turbulence
| |
| |publisher = Springer-Verlag, New York
| |
| |page = 164
| |
| |year = 1984
| |
| }}.
| |
| *{{citation
| |
| | last = Strogatz | first = S. | author-link = Steven Strogatz
| |
| | doi = 10.1016/S0167-2789(00)00094-4
| |
| | issue = 1–4
| |
| | journal = Physica D
| |
| | pages = 1–20
| |
| | title = From Kuramoto to Crawford: exploring the onset of synchronization in populations of coupled oscillators
| |
| | volume = 143
| |
| | year = 2000|bibcode = 2000PhyD..143....1S }}.
| |
| *{{citation
| |
| | last1 = Acebrón | first1 = Juan A.
| |
| | last2 = Bonilla | first2 = L. L.
| |
| | last3 = Vicente | first3 = Pérez
| |
| | last4 = Conrad | first4 = J.
| |
| | last5 = Ritort | first5 = Félix
| |
| | last6 = Spigler | first6 = Renato
| |
| | journal = Reviews of Modern Physics
| |
| | pages = 137–185
| |
| | title = The Kuramoto model: a simple paradigm for synchronization phenomena
| |
| | url = http://scala.uc3m.es/publications_MANS/PDF/finalKura.pdf
| |
| | volume = 77
| |
| | year = 2005|bibcode = 2005RvMP...77..137A |doi = 10.1103/RevModPhys.77.137 }}.
| |
| *{{citation
| |
| | last1 = Cumin | first1 = D.
| |
| | last2 = Unsworth | first2 = C. P.
| |
| | doi = 10.1016/j.physd.2006.12.004
| |
| | issue = 2
| |
| | journal = Physica D
| |
| | pages = 181–196
| |
| | title = Generalising the Kuromoto model for the study of neuronal synchronisation in the brain
| |
| | volume = 226
| |
| | year = 2007|bibcode = 2007PhyD..226..181C }}.
| |
| | |
| [[Category:Exactly solvable models]]
| |
| [[Category:Lattice models]]
| |
| [[Category:Partial differential equations]]
| |