数理解析研究所講究録 1244 巻 2002 年 47-55

47-

カオス場における間欠性と対称性
等方 2 次元反応拡散系を用いた解析-

京大院情報 本河知明 (Tomoaki Honkawa)藤坂博一 (Hirokazu Fujisaka)(Syuji Miyazaki)
Graduate School of Informatics, Kyoto Univ.
宮崎修次

概要
モードが縮退しているような対称性 (等方性) を持つ 2 次元空間において, 一様カオス場
力坏安定化する際, どのような現象が観測されるかについて調べた. 一様カオス場力坏安
定化した直後, その系の不変多様体近傍間を飛び移る現象が現れると予想していたが, 数1 次元的な振る舞いをすることが値実験の結果, その予想に反して対称性が破れており,分かった. またその破れた対称性は, 熱雑音下, またはシステムサイズを大きくしていく
ことで回復し, その際, 不変多様体近傍間を間欠的に飛び移る現象が観測された.

1

序論

振動子を複数個結合させた結合振動子系は, 単体の振動子では見られない多様な振る舞いを示すことが知られている. この結合振動子系で現れる現象の一つに, 同期・非同期転移がある. このような現象は, 振動子を多数結合させた系, 例えぱ, 空間自由度を持つ系でも観測されることが, スピン波不安定の実験などで知られている [1].
本研究では, モードが縮退しているような対称性 (等方性) を持つ 2 次元空間において,一様カオス場力坏安定化する際, どのような現象が観測されるかについて調べた.
今回我々が扱ったのは, 等方 2 次元反応拡散系である.

(1)$\ovalbox{\tt\small REJECT}(r, t)=F(X, t)+\hat{D}\nabla^{2}X$

この系は特解として $\hat{D}\nabla^{2}X=0$ となる一様解 (同期解) を持ち, 一様解軌道が $tarrow\infty$ で占$=0$
める不変多様体 $M$ が存在する. 2 次元系の場合はこの他にも,$\hat{D}\frac{\partial^{2}}{\partial x}\mathrm{z}X\neq 0,$$D^{2} \frac{\partial}{\partial y}\tau X$

を満たすと, その逆でを満たす特解が占める不変多様体の二つの不変多様体が存在する.
特解が占める不変多様体
方向にも 方向にも非一様な微小摂動を加えたとき, 系の安定性は
さて, この系にどうなるであろうか ? 系のコントロールパラメータを変化させていくと, 一様状態から非$\hat{D}\frac{\partial^{2}}{\partial x}\mathrm{Z}X=0,\hat{D}^{2}\frac{\partial}{\partial y}\tau X\neq 0$$oV_{e}$$V_{y}$$x$$y$

一様状態への転移が起こる. つまり, 不変多様体 $M$ が不安定化する転移点が存在する.
やも同時に不安定化するか否かである. も我々の興味は, この転移点において方向のモードと 方向のモードが共に不安定化し,し同時に不安定化するのであれば,近傍間を交互に飛び移る (パターンの向きが 方向と 方向とで交互に変化す$V_{l}$$x$$V_{\mathfrak{B}},$$V_{y}$$V_{y}$

$y$

$x$

$y$

48
る) ような現象が観測されることが予想される. もしこのように不安定モードが縮退して存在するならば, 縮退のない 1 次元系とは異なる振舞を示すはずである.

2 章で空間的一様な同期カオス場の線形安定性につぃて説明し, 次に 3 章
一様カオス場力坏安定化した直後, モードが縮退してぃる 2 次元系ではどのような振舞をするの力 , 1 次元系との比較を行う.
そこでまず,

で,$\mathrm{a}$2-

空間的一様な同期カオス場の線形安定性

空間自由度を持つカオス場の例として, 振幅 , 周波数
反応拡散系を考える.
$a$$\dot{X}(r,t)$$\Omega$の一様な周期外場を加えた,$=$$F(X,\mathrm{Y})+a\mathrm{c}\mathrm{o}\mathrm{e}\Omega t+D\nabla^{2}X$$=$$G(X,\mathrm{Y})+D\nabla^{2}\mathrm{Y}$

,
(2)

$\{$

$\dot{\mathrm{Y}}(r,t)$は反応物質の濃度場を表し, $F$ およひ $G$ は空間的に一様に起こる反応項である. $D$
は正の拡散係数であり, 反応物質に依らず, また方向にも依らない等方なものとする札$X,\mathrm{Y}$

この系は周期的境界条件の下で,$\dot{X}^{0}(r,t)$,$=$$F(X^{0},\mathrm{Y}^{0})+a\mathrm{c}\mathrm{o}\mathrm{e}\Omega t$$=$$G(X^{0}, \mathrm{Y}^{0})$

(3)

$\{$

$\dot{\mathrm{Y}}^{0}(r,t)$を満たす, 空間的に一様な特解を持っ. この一様解の最大リアプノフ指数は正とし, 一様解はカオス的に振る舞っているものとする.
$X^{0}(t),\mathrm{Y}^{0}(t)$$\lambda_{0}$

この一様解に非一様な微小摂動を加え, この摂動に対する線形安定性$X(r,t)=X^{0}(t)+\tilde{X}(r,
t),$
を考える.
を式 (2) に代入しフーリ
工変換を行なうと,のフーリエ係数は,$\tilde{X}(r, t),\tilde{\mathrm{Y}}(r,t)$$\mathrm{Y}(r, t)=\mathrm{Y}^{0}(t)+\tilde{\mathrm{Y}}(r, t)$$\tilde{X}(r, t),\tilde{\mathrm{Y}}(r, t)$$(\tilde{X_{k}.}.(t)\tilde{\mathrm{Y}}_{k}(t))=(\begin{array}{ll}F_{X}(t)-D|k|^{2} F_{\mathrm{Y}}(t)G_{X}(t) G_{\mathrm{Y}}(t)-D|k|^{2}\end{array})$$(\begin{array}{ll} \tilde{X}_{r}(t)-\sim \tilde{\mathrm{Y}}_{r}(t)\end{array})$(4).

に従う. 但し,ある.
さて, 波数 $k$ のモードの時間発展の漸近形は,$F_{X}(t)=\partial F(X, \mathrm{Y})/\partial X|_{X=X^{\mathrm{O}}(t),\mathrm{Y}=\mathrm{Y}^{0}(t)}$を表し, 他の記述もこれと同様で$\tilde{X}_{k}(t),\tilde{\mathrm{Y}}_{k}(t)\propto\exp[\Lambda_{k}(t)t]$のように表すことができる.

$\Lambda_{k}(t)$は $k$ モードの局所拡大率で,$\lambda_{k}\equiv\lim_{tarrow\infty}\frac{1}{t}\int_{0}^{t}\Lambda_{k}(s)ds=\lambda_{0}-Dk^{2}$1

一般的には拡散係数は R 応物質 |\breve \acute 依 ff$\mathrm{b}$,(5)
これを長時間平均した,(6)

また異方性のあるような系も存在するが, ここでは解析を簡単にするため, そういう依存性は持たせていない. パターン形成などの研究の際は, 依存性を持たせるのが普$\backslash .\mathrm{B}^{[]}7’ \mathrm{A}\mathrm{X}$

49

$*$

$k$

.

図 1: $k$ モードリアプノフ指数

を $k$ モードリアプノフ指数と言う.

この量の正負によって,する一様解の線形安定性を判断することができる.

つまり,

$k$

モードの微小な揺らぎに対

この量が負ならすのモードは安定であり, 正なら不安定である.
が正である場合を考えているので, $k=0$ 以外にも,
いま一様解のリアプノフ指数となる臨界波数 k。が存在し,ある波数領域のモードカ坏安定化しうる.
$\lambda_{0}$$\lambda_{k_{c}}=0$$)(7)$$k_{c}=\sqrt{\frac{\lambda_{0}}{D}}$と求まる. k>k。のモードは $\lambda_{k}<0$ となり安定であるが, k<k。のモードは $\lambda_{k}>0$ なので不安定化し, その波数に対応したサイズの空間構造が出現する. 言い換えると, シス

テムサイズ $L$ を大きくしていくと, 空間的一様解が安定から不安定へ転移する臨界シス
テ\Delta サイズ$L_{c}= \frac{2\pi}{k_{c}}=2\pi\sqrt{\frac{D}{\lambda_{0}}}$(8)

このことを示したのが図 1 である. L<L。では不安定モードは $k=0$ だけであるが, L>L。では, $k=0$ 以外にも不安定モードが存在し, その波数に対応したサが存在する.

イズの空間構造が現れる.
なお, 以下でモードを記述する際は,$k= \frac{2\pi}{L}n=\frac{2\pi}{L}(n_{x},.n_{y})$

のように,

$n$

,

を以って表すことにする. また,$n_{x},$$k$$n_{y}=0,$ $\pm 1,$ $\pm 2,$ $\cdots$

モードリアプノフ指数についても(9)$\lambda_{k}\equiv$\lambda (nx’n、のように表すことにする.

以下の章で行なう数値実験では, 反応項としてブラッセレータモデル$F(X, \mathrm{Y})=A+X^{2}\mathrm{Y}-(B+1)X$,$G(X, \mathrm{Y})=BX-X^{2}\mathrm{Y}$(10)を使用する [2, 8, 9]. パラメータは, 濃度場 $A=0.4,$ $B=1.2$ , 周期外場の振幅 $a=0.12$ ,周期外場の振動数 $\Omega=0.9$ , 拡散係数 $D=0.1$ と固定する. また, コントロールパラメー

50
夕としてシステムサイズ $L$ を採用する. このとき最大リアプノフ指数は0347 とと定まる. 3 章およひ 4 章ではなり, これから臨界システムサイズがL。
の範囲で計算し,章ではを考え, 1068
5
L\sim 2L。
$\lambda_{0}\ovalbox{\tt\small REJECT} \mathrm{Q}$$L_{c}\ovalbox{\tt\small REJECT} 1067$$L\ovalbox{\tt\small REJECT}$$(\ovalbox{\tt\small REJECT} 1 \omega 1L_{c})\ovalbox{\tt\small REJECT} L\ovalbox{\tt\small REJECT} 1121(\ovalbox{\tt\small REJECT} 105L_{c})$について計算する.

不安定化直後における 1 次元系と 2 次元系の比較

3

前章で行なった一様状態の線形安定性の議論より, 等方 2 次元反応拡散系における $k$

モードリアプノフ指数$\lambda k$

は

(11)

$\lambda_{k}=\lambda_{0}-D[(\frac{2\pi n_{x}}{L})^{2}+(\frac{2\pi n_{y}}{L})^{2}]$となる. つまり, この系が持っている等方性により, 一様解が不安定化する方向が $x$ 方向と$y$方向の二つが存在し, 不安定モードは縮退して存在すると予想される. 例えばシステ

ムサイズが L\sim >L。だとすると, $k=0$ 以外に $(\pm 1,0)$ モードと $(0, \pm 1)$ モードカ坏安定化すると予想される. そうであるなら, 不安定化する方向が一方向しかない 1 次元系とは異なる振舞が見られるはずである.
さて,2 次元系の数値実験について, 一様状態から

$x$

方向と

$y$

方向のどちらにどの程度外れているかを見るために以下の時系列を考える.
$l_{x}(t)$$\equiv$$\sqrt{|\tilde{X}_{(1,0)}(t)|^{2}+|\tilde{\mathrm{Y}}_{(1,0)}(t)|^{2}}$$l_{y}(t)$$\equiv$$\sqrt{|\tilde{X}_{(0,1)}(t)|^{2}+|\tilde{\mathrm{Y}}_{(0,1)}(t)|^{2}}$

,

(12)

(13)

また, 直接 $l_{x}(t),$ $l_{y}(t)$ の時系列を見る以外に, 次のような物理量$\theta(t)=\arctan\frac{l_{y}(t)}{l_{x}(t)}$

,

(14)

$(0\leq\theta\leq\pi/2)$

方向にパターンができて 方向に一様なときは $\theta(t)=0$ , 逆に 方向にパを定義する.
ターンができて 方向に一様なときは $\theta(t)=\pi/2$ を示すよう定義してあり, この量を見
$x$

$y$

$y$

$x$

ることで $x$ 方向と$y$方向のいずれの方向にパターンができているかを知ることができる.
この結果は予想と異なり, 過渡

図 2 は $L=10.70(=1.0028L_{c})$ における時系列である.

状態を経た後にどちらか一方の不変多様体に吸引され, 不変多様体

近傍間を飛ひ移るという現象は観測されないことを示している. つまり, 一様状態が不安定化した直後
では, 等方性は破れており, 1 次元的な振舞しか示さない.
さて, この 1 次元的な振舞を適当な統計量を用いて確認してみよう. 2 次元系では 方
向と 方向のいずれに軌道が落ち込むか分からないので, 空間的一様状態からのずれを表す変数 $l(t)$ を次のように定義する.
$aV_{e},$$aV_{e}$$x$$y$$l(t)=\sqrt{|l_{x}(t)|^{2}+|l_{y}(t)|^{2}}$

.

(15)

51

$\backslash$

ごさ$\iota\pi$$\mathrm{t}\Pi$(a) 時系列$\theta(t)$.

(b) 時系列$l_{x}(t),$ $l_{y}(t)$.

近傍の異なる初期状態から時間発展させた軌道, (b) はそのうちの一つのはとの時系列をプロットしてある. (b) は過渡状態後,軌道についてのみ不安定化せず,のみが間欠的に強く不安定化している. $(L=10.70)$
図 2:(a) は$oV_{e}$$l_{x}(t)$$l_{y}(t)$$l_{y}(t)$$l_{x}(t)$また,1 次元系については(16)$l(t)=\sqrt{|\tilde{X}_{1}(t)|^{2}+|\tilde{\mathrm{Y}}_{1}(t)|^{2}}$とする.

一般に結合振動子系などにおいて, カオス特解である同期状態力坏安定化する際, オン
オフ間欠性と呼ばれる現象が観測されることが知られている [3, 4, 5, 6, 7]. この現象は,少数自由度系だけでなく, 大自由度系においても 1 次元系については観測されることが報告されている [2]. もし 2 次元系でも 1 次元的な振舞をしているのであれば, 1 次元系と同じ統計性を示すはずである.
オンオフ間欠性に特徴的な統計則として, (1) バーストの強度 (振幅戸の分布 $P(l),$ $(2)$の三つがよくラミナー状態の継続時間 の分布 $Q(\tau)$ , (3) $l(t)$ のパワースペクトル調べられている. 我々も 1 次元系と 2 次元系について, これら統計則の比較を行なった.
図 3 に示したように, 1 次元系と 2 次元系とで統計則は一致した. これは上でも述べたよ
次元的な振舞をしているか
うに, 2 次元系においても, 一様状態力坏安定化. し’ た直$I(\omega)$$\tau$$\acute{\{}\mathrm{g}\mathrm{B}\mathrm{V}\overline{\mathrm{i}}$らである.

本章をまとめると, 一様解の線形安定性の議論では, 一様解力坏安定化した直後は $(\pm 1,0)$
モードと $(0,\pm 1)$ モードが同時に不安定化すると思われていたが, 実は等方性が破れており,

どちら力\vdash 方のモードが不安定化したら, もう一方は不安定化せず,

1 次元的なパターン

しかできないことが分かった.

この破れた対称性はどのようにしたら回復するのであろうか ? そこで, 4 章ではシステムサイズは本章と同じままにして熱雑音を加えた場合, 5 章では熱雑音は加えずに,
システ\Delta . サイズを大きくしていった場合について考えることにする.
では,

52

嘗

$6\mathrm{e}$(a) $P(l)$(b) $Q(\tau)$(c)$I(\omega)$

図 3: オンオフ間欠性の特徴的な統計性について,

1 次元系と 2 次元系の比較. いずれの統計性も 1 次元系と 2 次元系とで一致しており, 2 次元系においても 1 次元的な振舞をしている.

4

熱雑音の影響前章で, 一様状態力坏安定化しオンオフ間欠性が出現するのと同時に, 対称性 (等方性)が破れることが数値実験により見出された. しかし, 現実世界の現象においても対称性の破れが観測されるとは限らないであろう. 何故なら, 数値実験を行なったモデルは理想化
された決定論的方程式であるのに対し, 現実世界では熱雑音の影響を受けると考えられるからである. 果して熱雑音下でも対称性は破れているのか ? それとも対称性を保持しているのか ?

ここでは, 熱雑音としてガウス白色ノイズ $\xi x(r,t),$$\xi_{\mathrm{Y}}(r, t)$を加えた$\dot{X}(r,t)$$=$$F(X,\mathrm{Y})+a\mathrm{c}\mathrm{o}\mathrm{e}\Omega t+D\nabla^{2}X+\xi \mathrm{x}(r,t)$$\dot{\mathrm{Y}}(r, t)$$=$$G(X, \mathrm{Y})+D\nabla^{2}\mathrm{Y}+\xi_{\mathrm{Y}}(r,t)$

(17)

$\{$

を考える [11]. 相加ノイズ

$\xi$

は
$\langle\xi_{l}(r,t)\rangle=0$,$\langle\xi\iota(r,t)\xi_{l’}(r’,t’)\rangle=2\gamma\delta_{l,l’}\delta(r-r’)\delta(t-t’)$の性質を持つとする.

$\gamma$(18)$l=\{X,\mathrm{Y}\}$

,

$l,$

(19)

$l’=\{X,\mathrm{Y}\}$はノイズの強度を表す. また, 境界条件は周期的境界条件としてある.
図 4 は, 相加ノイズ下でのの時系列である. このとき, 不変多様体
近傍間を間欠的に飛ひ移る現象が観測され, 対称性が回復している. っまり, 熱雑音下では対称性は破れない.
$\theta(t)$$aV_{e},V_{y}$53(a) $L=10.68,\gamma=10^{-8}$(c)(b) $L=11.00,\gamma=10^{-8}$

図 4: 相加ノイズを加えた時の時系列$\theta(t)$$L=\mathrm{I}\mathrm{L}\mathrm{O}\mathrm{O}$, $\gamma=10^{-12}$. ノイズを加えると, 対称性が回復する.

L。から離れるに従い, また, ノイズの強度$\gamma$が小さくなるに従い,$oV_{e},$$V_{y}$$L$が近傍に滞在する時間が長くなっている.

システ\Delta サイズと対称性

5

前章では熱雑音を加えることで対称性が回復することを見たが, では, 熱雑音なしで,システムサイズを大きくしていくことにより対称性は回復しないのだろうか ?1 次元的なパターンではなく 2 次元的なパターンが現れるには $(\pm 1,\pm 1)$ モードの不安定化が必要だと考えられる. そこで, システムサイズ $L$ を変化させたときの分岐図を図 5 に示す. 図 5より, 線形安定性の議論から予想される L\sim >\psi L。の領域では周期窓が現われているが,さらにシステムサイズを大きくしていくと,

1 次元的な振舞から 2 次元的な振舞への転移が見られることが分かる.

$\frac{8}{\approx}$(a) $L/L_{c}$ :$l_{(1,0)}$(b) $L/L_{c}$ :(c) $L/L_{c}$ :$l_{(0,1)}$$l_{(1,1)}$

によって, 1 次元的な周期運
2 次元系の 090L。<L<2.10L。における分岐図.
動が現れる領域, 1 次元的な非周期運動の領域, 2 次元的な周期運動の領域, 2 次元的な非周期運動の領域が見える. L/L。を大きい方から 005 刻みずつプロット.

図 5:

そこで次に,$L$

2 次元的なカオス運動が現われる L\sim >1.95L。について時系列を調べた.

どちらの不変多様体近傍に滞在しているかを示す$\theta(t)$の時系列と,$(1, 1)$

モードの時系列

54を示したのが図 6 である.

これより, 不変多様体近傍間を遷移するタイミングと同時に,(強度) 力状きくなっていることが分かる. このことは, パターンの向きの変化には $(1,1)$ モードの不安定化が関係していることを示唆している.

$(1,1)$ モードの振幅

このことを確認するために, 図 7 を見てみよう. これより, $\sqrt{l_{(1,0)}(t)^{2}+l0,1)((t)^{2}}$ の振幅が大きくなってもは 0 付近しか取らず, パターンの向きをあまり変えないが, $l_{(1,1)}(t)$$\delta\theta(t)$の振幅力状きくなると,$\delta\theta(t)$ターンの向きを変えるには,(a)も大きく動くことが可能であることが分かる. っまり,
モードが重要である.

パ$(\pm 1,\pm 1)$$\theta(t)$図7:$l_{(1,1)}(t)$$\sqrt{l_{(1,0)}(t)^{2}+l_{(0,1)}(t)^{2}}$:$\delta\theta(t)$:$\delta\theta(t)$と. $(1,0)$ モードや $(0,1)$ モードの強度が大きくてもはあまり大きく変化
しないが, $(1,1)$ モードカ状きくなるとは$\delta\theta$$\delta\theta$

大きく変化しうる. $\delta\theta(t)=\theta(t)-\theta(t-T)$ .

(b) $l_{(1,1)}(t)$図 6:(a)$\theta(t)$およひ(b) $(1,1)$ モードの時系列.
$(L=21.50=2.01L_{c})$6結論
モードが縮退している等方 2 次元系において, カオス的な一様状態が不安定化した直後
1 次元的な振舞をしており, オンオフ間欠性が観測されるでは, 対称性 (等方性) が破れ,

55
ことが分かった.
この破れた対称性は,熱雑音がないときは対称性が破れているようなパラメータ領域であっても, 熱雑音を加えることによって対称性が回復した.
また, 熱雑音がなくても, コントロールパラメータであるシステムサイズを大きくする
ことによっても対称性が回復することを確かめた. この対称性の回復には $(\pm 1,\pm 1)$ モードの不安定化が重要であることが分かった.

参考文献[1] A. Krawiecki and A. Sukiennicki: Acta Phys. Pol. 88(1995), 269. ;F. T. R\"odelsperger, A. enys and H. Benner: Phys. Rev. Lett. 75(1995), 2594. ;J. Becker, F. R\"odekperger, Th. Weyrauch, H. Benner, W. Just and A. enys: Phys.
$\check{\mathrm{C}}$$\check{\mathrm{C}}$Rev. $\mathrm{E}59(1999)$ , 1622.
[2] H. Fujisaka, K. Ouchi, H. Hata, B. Masaoka and S. Miyazaki: Physica $\mathrm{D}114(1998)$ ,237.
[3] H. Ehjisffia and T. Yamada: Prog. Theor. Phys. 74(1985), 918; 75(1986), 1087.
[4] 藤坂博–: 応用数理$\mathrm{V}\mathrm{o}\mathrm{l}.9$

No 1, 28 (岩波書店, 1999)[5] E. Ott and J. C. Sommerer: Phys. Lett. A188(1994), 39.
[6] Y. -C. Lai and C. Grebogi: Phys. Rev.

$\mathrm{E},$$52$(1995), R3313.

[7] M. Zhan and G. Hu: Phys. Rev. $\mathrm{E}62(2000),$ $375$ .
[8] T. Kai and K. Tomita: Prog. Theor. Phys. Suppl. 64(1978), 280. ;K. Tomita and T. Kai: Prog. Theor. Phys. 61(1979), 54.

[9] T. Yamada and H. Hhjisaka: Prog. Theor. Phys. 70(1983), 1240.
[10] M. Abramowitz and I. A. Stegan (Eds.): Handbook of Mathematical Functions,Chap. 25, Dover Publications, New York, 1972.

[11] A. C.

$\mathrm{B}\mathrm{r}\mathrm{a}\acute{\mathrm{n}}\mathrm{k}\mathrm{a}$and D. M. Heyes: Phys. Rev. $\mathrm{E}58(1998)$ , 2611.