蔵本モデル
蔵本モデルは蔵本由紀によって提案された同期現象を記述する数学モデルである。特に、相互作用のある非線形振動子集団の振る舞いを記述するモデルである。このモデルは化学的、生物学的な非線形振動子系の振る舞いを示唆するものであり、幅広い応用が見られる。
このモデルの前提として、完全に独立（またはほぼ独立）した振動子に弱い相互作用がはたらくこと、そしてこの相互作用は二つの振動子間の位相差の正弦関数として与えられる、という仮定がある。
定義
最も知られた形式の蔵本モデルの場合、各々の振動子らは固有振動数[math]\omega_i[/math]を持ち、他の全ての振動子と等しく相互作用している、と考えられる。驚くべきことに、この非線形モデルは[math]N \to \infty[/math]の極限において上手く変形することで厳密に解くことができる。
最も知られた蔵本モデルの形式は次のような支配方程式に従う。
[math] \frac{\partial \theta_i}{\partial t} = \omega_i + \frac{K}{N} \sum_{j=1}^{N} \sin(\theta_j - \theta_i), \qquad i = 1 \ldots N[/math],
ここで、系はN個のリミットサイクル振動子から構成される。
また、系にノイズを加えることができる。この場合、方程式は書き換えられて、
[math] \frac{\partial \theta_i}{\partial t} = \omega_{i}+\zeta_{i}+\dfrac{K}{N}\sum_{j=1}^N\sin(\theta_{j}-\theta_{i}) [/math],
ここで、[math]\zeta_{i}[/math]は揺らぎを表し、時刻の関数である。ホワイトノイズを考えれば、
[math]
\langle\zeta_{i}(t)\rangle=0
[/math] ,
[math]
\langle\zeta_{i}(t)\zeta_{j}(t')\rangle=2D\delta_{ij}\delta(t-t')
[/math]
となる。ここで[math]D[/math]はノイズの強さを表す。
変形
蔵本モデルは次のようになる。 「秩序」パラメータrとψを次のように定義する。
[math]re^{i \psi} = \frac{1}{N} \sum_{j=1}^{N} e^{i \theta_j} [/math].
ここでr、ψは振動子集団の平均場の振幅、位相である。 この変形を適用することで（任意の N で）支配方程式は次のようになる。
[math] \frac{\partial \theta_i}{\partial t} = \omega_i + K r \sin(\psi-\theta_i) [/math].
こうして振動子の方程式はもはや陽的には結合されてはおらず、その代わりに秩序パラメータが振る舞いを決める。振動子集団の位相分布が均一であれば（少なくとも[math]N \to \infty[/math] の極限で）更に変形が行われて [math]\psi=0[/math]となり、支配方程式は次のようになる。
[math] \frac{\partial \theta_i}{\partial t} = \omega_i - K r \sin(\theta_i) [/math].
Nが大きい場合の極限
[math]N \to \infty[/math] の場合を考えよう。 固有振動数の分布がg（ω）（正規化されていると仮定する）で表されるとする。 時刻tでの位相θ、固有振動数ωにおいて、振動子の密度が [math]\rho(\theta, \omega, t)[/math]であるとする。 正規化の要請から次の式を満たす。
[math] \int_{-\infty}^{\infty} \rho(\theta, \omega, t) \, d \theta = 1. [/math]
振動子の密度の連続の式は次のようになる。
[math] \frac{\partial \rho}{\partial t} + \frac{\partial}{\partial \theta}[\rho v] = 0, [/math]
ここで、vは振動子のドリフト速度であり、[math]N \to \infty[/math] における支配方程式の変形から、
[math] \frac{\partial \rho}{\partial t} + \frac{\partial}{\partial \theta}[\rho \omega + \rho K r \sin(\psi-\theta)] = 0. [/math]
最後に、[math]N \to \infty[/math]での秩序パラメータの定義を書き直そう。 [math]\theta_i[/math]は（ωでの）アンサンブル平均で、和は積分で置き換えられるので次のようになる。
[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]
解
全ての振動子がランダムに動くインコヒーレントな状態の解は [math]\rho = 1/(2\pi)[/math]に対応する。 [math]r = 0[/math]の場合、振動子の間に全く相関は無い。 集団の振動子の位相分布が一様であれば、集団は静的に安定な状態である（けれども個々の振動子の位相は自らの固有のωに従って変化し続けている）。
Kが十分強いとき、完全に同期した解が実現する。 完全に同期した状態では、全ての振動子は、個々の位相は異なれども、共通の振動数をとる。
部分的に同期した場合の解は、 固有振動数の値が近い幾つかの振動子のみが同期し、他の振動子はばらばらに動く状態を引き起こす。 数学的には、同期した振動子は、
[math]\rho = \delta\left(\theta - \psi - \arcsin\left(\frac{\omega}{K r}\right)\right)[/math]
となり、ばらばらに動く振動子は、
[math]\rho = \frac{\rm{normalization \; constant}}{(\omega - K r \sin(\theta - \psi))}[/math]
となる。振動子は[math]|\omega| \lt K r [/math]の場合は同期でき、そうでない場合はばらばらな動きになる。
関連分野
脚注
- ↑ Xiao Fan Wang and Guanrong Chen (2003). “Complex Networks: Small-World, Scale-Free and Beyond”. IEEE CIRCUITS AND SYSTEMS MAGAZINE 3 (1): 16-19 . 2013年3月29日閲覧..
参考文献
- 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 [1]
- Strogatz, S. From Kuramoto to Crawford: exploring the onset of synchronization in populations of coupled oscillators, Physica D 143 (2000) 1–20