next up previous
次へ: 4 無衝突系におけるカオスの役割 上へ: 理論天体物理学特論I 戻る: 2 形成、相互作用の確率

3 カオス

そもそもカオスとは何でなぜそういうものを考えるかというところから話を始め る必要がある。

というわけで、定義から始めたいところだが、一応ここでは

ものをカオスということにする。

そういうものはいくらでもあるわけだが、まず、カオスで「ない」ものの例か ら考えていこう。

球対称ポテンシャルの中での粒子の運動はカオスではな い。これは、軌道が単純な動径方向の振動と角度方向の振動に分離できるから である。

これが、軸対称ポテンシャルになるともういきなり話がわからなくなる。 軸対称なので円筒座標 $(R,z)$ を考えると、ポテンシャル $U_g(R,z)$ の中 で $z$ 軸回りの角運動量 $L_z$ を持つ粒子の運動は、有効ポテンシャル

\begin{displaymath}
U(R,z) = U_g(R,z) +\frac{L_z^2}{2R^2}
\end{displaymath} (7)

の中での平面運動と同じである。従って、一般論として、軸対称ポテンシャル の中での運動を考えるには、単なる2次元平面$(x,y)$ の中での運動を考えれ ばよい。言い換えると、軸対称ポテンシャルの中での粒子の運動というのは、 極めて一般的な、2次元平面の中での運動と等価である。

Henon と Heiles (1964 AJ 69, 73) は、 $U$ として以下のようなものを考 えてみた。


\begin{displaymath}
H = \frac{1}{2}(x^2 + y^2) +x^2y - \frac{1}{3}y^3
\end{displaymath} (8)

何故この形を考えたかというのは、論文には以下のように書いてある

because: (1) it is analytically simple; this makes the computation of the trajectory easy; (2) at the same time, it is sufficiently complicated to give trajectories which are far from trivial, as will be seen below.

つまり、物理的あるいは天文学的になにかの意味があるというよりは、このポ テンシャルの中での軌道の振る舞いが妙であるということで選ばれたかなり人 工的な例ではある。

もちろん、最初の項は調和振動であり明確な意味がある。非線型の第2項の意 味はポテンシャルの等高線を書いて見ると分かる。

12cm \epsffile{HH1964fig1.ps}

このように、非線型項によりポテンシャルが三角形になっている。特に $U=1/6$ の等高線に注意すると、これは直線になっている。このことから、こ の三角形の頂点はポテンシャルの鞍点であり、力が 0 になることが分かる。

鞍点であるので、ポテンシャルが下がるのは三角形の内側を向いた方向とその 反対側の方向だけであり、それに直交する方向ではポテンシャルは増えている。 また、頂点以外では力は 0 でない。

軌道は例えばこんなふうである。これはエネルギー(ハミルトニアン)

\begin{displaymath}
H=\frac{v^2}{2} + U(x,y)
\end{displaymath} (9)

の値を $H=0.1$ とし、3つの違う初期条件から数値的に軌道を計算した結果である。

6cm \epsffile{henonwp0_0.ps}

6cm \epsffile{henonwp0_16.ps}

6cm \epsffile{henonwp0_45.ps}

これだけを見ていると別にどうということはないが、ここでいわゆるポアンカ レ断面というものを書いて見ると非常に面白いことになっていると分かる。ポ アンカレ断面とは、 この場合のような自由度が2のハミルトン力学系で、一方 の変数、例えば $x$ が 0 になる時に $y$, $v_y$ の値を書いたものである。

まず、 $H=1/12$ の場合を示す。

8cm \epsffile{HH1964fig4.ps}

ここで、つながって線になっているのは基本的には一つの軌道であり、分かれ た曲線は別の初期条件からの軌道に対応している。エネルギー保存を考えると、 $x=0$$v_x^2 \ge 0$なので、 $(y, v_y$) はある閉領域の中にくる。保存 量がエネルギーだけであれば、任意の初期条件から出発した軌道はこの領域の 全ての点をいつかは通ると考えられるが、 このようにそうならないというこ とは、なんだかわからない保存量のようなものがあるということを示している。 元々の軸対称ポテンシャルに戻ると軸回り角運動量は保存しているが、これは 有効ポテンシャルにして回転方向の速度を消去するのに使っているので、2次 元問題になった後では無関係である。

軸対称ポテンシャル内での運動にエネルギーと角運動量以外の保存量があるか どうかという問題は、「第3積分問題」と呼ばれる恒星系力学上の現在でも完 全には解決されているとはいいがたい問題の一つである。 とりあえずこの場 合には、第3積分が存在しているように見える。

但し、妙なのは、ポアンカレ断面上の軌跡が交点をもつように見えることであ る。実際には、この「交点」のところを精密に調べると、奇妙なことがわかる。 交点近くの十分に狭い領域では運動の積分が無くなって、ある領域内を軌道が 埋めているのである。

この領域はエネルギーを大きくするとどんどん大きくなる。次は $E=1/8$ の場合である。

8cm \epsffile{HH1964fig5.ps}

保存量を持つ、規則的な運動をする領域よりも、規則的ではない運動をする領 域のほうが広くなってしまっていることが分かる。少なくとも恒星系力学では、 このような場合に 規則的な軌道を regular orbit, そうでないものを irregular orbit または chaotic orbit という。

さらにエネルギーの上限に近い $-1/6$ をとったものが下である。

8cm \epsffile{HH1964fig6.ps}

キャプションには $1/6$ よりも大きいようなことが書いてあるが、気の迷い であろう。それはともかく、この場合にはほとんどの領域が規則的でない、つ まり、第3積分をもたないような軌道になっていることが分かる。

これはカオス的な運動の典型的な例である。

ここではカオスかどうかを、エネルギー以外の積分がある(ように見える)かど うかで判断したが、実質的に同じことになるもうひとつの定義がある。それは、 「十分に近い初期条件から出発した2つの軌道がどのように離れていくか」で ある。より厳密にはこれはリヤプノフ指数が正かどうかということになる。定 義の詳細はここでは省くが、大雑把にはリヤプノフ指数が正なら十分近い2つ の軌道が指数関数的に離れていく。このことは上の Henon-Heiles 系でも観 察され、周期的に見える領域では軌道は時間の1次で離れていくが、カオス的 な領域では指数関数的に離れる。


next up previous
次へ: 4 無衝突系におけるカオスの役割 上へ: 理論天体物理学特論I 戻る: 2 形成、相互作用の確率
Jun Makino 平成20年7月9日