next_inactive up previous
上へ: Stellar Dynamics Lecture Note

pt 称

重力多体系の進化

牧野淳一郎

1 ジーンズ不安定

ここまで、自己重力多体系の平衡形状を(ごく簡単な場合だけ)扱ってきたが、 ここからしばらくは平衡形状ではなく時間発 展について考えることにしたい。

時間発展といっても、もとの方程式が強い非線形性を(加速度の項に)持つの で、一般的な場合を解析的に扱うことはほとんど出来ない。そこで、ま ず、平衡状態から無限小だけずれている場合に対して線形化した発展を考える ということにする。さらに、話を簡単にするために、「無限一様」な平衡状態 とし、まず無衝突ボルツマン方程式ではなく流体の話を考える。

1.1 流体のジーンズ不安定

流体は、連続の式

\begin{displaymath}
{\partial \rho \over \partial t} + \nabla\cdot(\rho {\bf v}) = 0
\end{displaymath} (1)

オイラー方程式
\begin{displaymath}
{\partial {\bf v}\over \partial t} + ({\bf v}\cdot \nabla){\bf v}= - {1 \over
\rho} \nabla p - \nabla \Phi
\end{displaymath} (2)

ポアソン方程式
\begin{displaymath}
\nabla^2 \Phi = 4\pi G \rho
\end{displaymath} (3)

で記述される。さらに状態方程式がいる。これはいま圧力が密度だけの関数 で与えられるとする。(断熱でも等温でもなんでもいいが)

今、 $\rho, p, {\bf v}, \Phi$ をそれぞれ $\rho = \rho_0 +
\rho_1 $ という格好にして、添字 0 がつくものはもとの方程式の平衡解であ り、 $1$がつくものは小さい(二次以上の項を無視していい)として方程式を 書き直すと


\begin{displaymath}
{\partial \rho_1 \over \partial t} + \nabla\cdot(\rho_0 {\bf v}_1) + \nabla\cdot(\rho_1 {\bf v}_0) = 0
\end{displaymath} (4)


\begin{displaymath}
{\partial {\bf v}_1 \over \partial t} + ({\bf v}_0 \cdot \na...
..._0^2} \nabla p_0 - { 1\over \rho_0} \nabla p_1-
\nabla \Phi_1
\end{displaymath} (5)


\begin{displaymath}
\nabla^2 \Phi_1 = 4\pi G \rho_1
\end{displaymath} (6)


\begin{displaymath}
p_1 = \left({dp \over d\rho}\right)_0 \rho_1 = v_s^2 \rho_1
\end{displaymath} (7)

ここで $v_s$ は音速である。

既に述べたように、無限一様でいたるところ密度、圧力が等しいというのが平 衡で、速度も0だったとすると、上の2本は

\begin{displaymath}
{\partial \rho_1 \over \partial t} + \rho_1 \nabla\cdot(\rho_0 {\bf v}_1) = 0
\end{displaymath} (8)


\begin{displaymath}
{\partial {\bf v}_1 \over \partial t} = - { 1\over \rho_0} \nabla p_1-
\nabla \Phi_1
\end{displaymath} (9)

となる。下2本は見かけはかわらない。これを、 $\rho_1$ だけの式にすれば
\begin{displaymath}
{\partial^2 \rho_1 \over \partial t^2} - v_s^2 \nabla^2 \rho_1 - 4
\pi G \rho_0 \rho_1 = 0
\end{displaymath} (10)

さて、めでたく方程式が線形化されたので、分散関係を求めれば算数はおしま いだが、まずその前にどういう方程式かということを見ておこう。

最初の2項をみれば普通の波動方程式で、最後の項がポアソン方程式を通して でてくる重力の項である。したがって、波長が短い極限では普通の波動方程式 に近付く。これに対し、波長が長い極限では空間2階微分の項が効かなくなる ので、線形の常微分方程式になってしまう。

実際に分散関係を求めるために、解を

\begin{displaymath}
\rho_1 = Ce^{i({\bf k}\cdot {\bf x}- \omega t)}
\end{displaymath} (11)

として代入すれば
\begin{displaymath}
\omega^2 =v_s^2 k^2 - 4\pi G\rho_0
\end{displaymath} (12)

ということになる。したがって、
\begin{displaymath}
k_J^2 = {4\pi G \rho_0 \over v_s^2}
\end{displaymath} (13)

と書くと、

なお、一応念のために書いておくと、式(11)の形の解だけを 考えるのは任意の初期条件からの解が(連続性とかを仮定すれば)この形の解 の線形結合で表現できるからである。解の線形結合が解であるのは方程式が線 形だからであり、任意の解が表現できるのは要するにフーリエ変換が完全系を なすからである。

話を戻すと、定性的に見たときにわかったように、波長が短ければ普通の音波 として振舞うが、波長が $1/k_J$ より長いと時間の指数関数で進化すること になる。つまり、密度が上がり始めたらどんどんあがる(下がり始めたらどん どんさがる)ということになる。

ここで注意して欲しいのは、十分に波長が長いと必ず不安定になるということ である。これはつまり、重力があると無限に一様な状態というのは温度無限大 でない限り必ず不安定であるということである。

さて、 $k_J$ というのは直観的にはどういうものかというのを考えてみる。 $k_J$ は波数なので、その逆数に $2\pi$ を掛けると波長になる。これをジー ンズ波長 $\lambda_J$ といって、式にすると

\begin{displaymath}
\lambda_J= \sqrt{\frac{\pi}{G\rho_0}}v_s
\end{displaymath} (14)

である。

今、半径がジーンズ波長くらいの球を考えてみる。これの単位質量当りの運動 エネルギー(熱エネルギー)は、もちろん $v_s^2$ の程度である。これに対し て重力エネルギーは $GM_J/\lambda_J$ の程度、ここで $M_J$ はジーンズ質 量で半径 $\lambda_J$ の球の質量である。 $\lambda_J$ に上の式を入れて計算すると、結局 重力エネルギーが定数を別にして $v_s^2$ の程度になることがわかる。

つまり、熱エネルギーより重力エネルギーが大きくなるような長さの摂動は成 長するということになる。

2 恒星系でのジーンズ不安定

恒星系の場合を考えることにしよう。

2.1 CBEの線形化

何度も出てきたがもう一回式を書いておく。


\begin{displaymath}
{\partial f \over \partial t} + {\bf v}\cdot \nabla f - \nabla \Phi \cdot
{\partial f \over \partial {\bf v}} = 0,
\end{displaymath} (15)

ここで $f$は6次元位相空間での分布関数である。$\Phi$ は重力ポテンシャ ルであり、以下のポアソン方程式の解として与えられる。
\begin{displaymath}
\nabla ^2 \Phi = - 4\pi G \rho.
\end{displaymath} (16)

$G$ は重力定数であり、 $\rho$ は空間での質量密度
\begin{displaymath}
\rho = m\int d{\bf v}f,
\end{displaymath} (17)

である。 これを流体の時と同様に線形化して、その振舞いを調べる。分布関数を $f_0
+ f_1$、ポテンシャルを $\Phi_0 + \Phi_1$とし、添字0がつくほうは定常解で あるとして式を整理すれば
\begin{displaymath}
{\partial f_1 \over \partial t} + {\bf v}\cdot \nabla f_1
- ...
...\nabla \Phi_1 \cdot {\partial f_0 \over \partial {\bf v}}
= 0,
\end{displaymath} (18)


\begin{displaymath}
\nabla ^2 \Phi_1 = - 4\pi G \int f_1 d{\bf v}.
\end{displaymath} (19)

ということになる。これが線形化された無衝突ボルツマン方程式である。

2.2 ジーンズ不安定

上の線形化されたボルツマン方程式は、平衡解の回りならば何でも使えるが、 ここではもっとも簡単に解析できる場合として空間分布が一様な場合を考える。 これから $f_0$は速度だけの関数であり、 $\Phi_0$ は定数としていいことに なるので

\begin{displaymath}
{\partial f_1 \over \partial t} + {\bf v}\cdot \nabla f_1
- \nabla \Phi_1 \cdot {\partial f_0 \over \partial {\bf v}}
= 0,
\end{displaymath} (20)

と少し式が簡単になる。まず、流体の場合と同じような平面波型の解を考えて みよう。


$\displaystyle f_1$ $\textstyle =$ $\displaystyle f_a({\bf v}) \exp[i({\bf k}\cdot {\bf x}- \omega t)]$ (21)
$\displaystyle \Phi_1$ $\textstyle =$ $\displaystyle \Phi_a \exp[i({\bf k}\cdot {\bf x}- \omega t)]$ (22)

速度空間の方にも伝わっていく波とかいうのも考えられないわけではないが、 とりあえずそういうのは考えない。これらを上の線形化した式に入れれば
\begin{displaymath}
(-\omega + {\bf v}\cdot {\bf k}) f_a({\bf v}) - \Phi_a {\bf k}\cdot {\partial f_0
\over \partial {\bf v}} = 0
\end{displaymath} (23)


\begin{displaymath}
-k^2\Phi_a = 4\pi G\int f_ad {\bf v}
\end{displaymath} (24)

となる。これらから $f_a$ を消せば、 $\Phi_a$ も落ちて
\begin{displaymath}
1 + {4 \pi G \over k^2}\int {{\bf k}\cdot {\partial f_0 \ove...
...artial
{\bf v}} \over {\bf k}\cdot {\bf v}- \omega}d{\bf v}= 0
\end{displaymath} (25)

となって、$f_0$が与えられていれば $k$$\omega$ の関係、すなわち分散 関係を与える。

もっとも、これはちょっと困った式で、 ${\bf k}\cdot {\bf v}- \omega = 0 $が特 異点になっている。したがって、実数の振動数を考えるのはすこし厄介な話に なる。まず、臨界点、すなわち振動数が $0$の場合と、不安定、すなわち振動 数が純虚数の場合を考えよう。

3 臨界点

式が繁雑なので、いま、 波数ベクトルを$x$軸方向にとることにすると、分散 関係は

\begin{displaymath}
1 + {4 \pi G \over k^2}\int {k {\partial f_0 \over \partial
v_x} \over kv_x - \omega}d{\bf v}= 0
\end{displaymath} (26)

臨界安定で $\omega=0$とすれば、結局
\begin{displaymath}
k_J^2 = -4 \pi G \int { {\partial f_0 \over \partial
v_x} \over v_x }d{\bf v}
\end{displaymath} (27)

となる。 $f_0$$v_x=0$ で有限で微分可能なら積分は求まるので、これか ら $\omega=0$ となる波数 $k_J$は決まる。

流体の場合と対応をつけるために、速度分布 $f_0$をマックスウェル分布

\begin{displaymath}
f_0({\bf v}) = {\rho_0 \over (2\pi \sigma^2)^{3/2}} \exp\left({- v^2/2\over \sigma^2}\right)
\end{displaymath} (28)

にしてみる。 これは例によって
\begin{displaymath}
\int_{-\infty}^{\infty}e^{-{v^2 \over 2 \sigma^2}}dv = \sqrt{2\pi
\sigma^2}
\end{displaymath} (29)

を使って全部積分できて、
\begin{displaymath}
k_J^2 = {4 \pi G \rho_0 \over \sigma^2}
\end{displaymath} (30)

となる。これは、流体の場合と同じになっている。つまり、中立安定な波長 (ジーンズ波長)は、恒星系と流体で同じである。

4 不安定な場合

次に、振動数が純虚数の場合を考えてみる。この時は、 $\omega = i \gamma$ として元の式にマックスウェ ル分布を入れて整理すると

\begin{displaymath}
1 + {2\sqrt{2} \pi G\rho_0 \over k\sigma^3}\int
{v_x e^{-v_x^2/2\sigma^2} \over
kv_x - i \gamma}dv_x = 0
\end{displaymath} (31)

ここで、被積分関数の分母を実数にするために
\begin{displaymath}
{v_x (kv_x + i\gamma)e^{-v_x^2/2\sigma^2} \over
k^2v_x^2 + \gamma^2}dv_x = 0
\end{displaymath} (32)

と書き直すと、虚部は奇関数なので落ちる。実部は
\begin{displaymath}
\int_0^{\infty} {x^2 e^{-x^2} \over x^2 + \beta^2} dx =
{ 1 ...
...rt{\pi} - {1 \over 2} \pi\beta
e^{\beta^2}[1-{\rm erf}(\beta)]
\end{displaymath} (33)

なる関係を使って
\begin{displaymath}
k^2 = k_J^2\left\{ 1- {\sqrt{\pi}\gamma \over \sqrt{2}k\sigm...
...rf}\left({\gamma
\over \sqrt{2}k\sigma}\right)\right]\right\}
\end{displaymath} (34)

なんだかよくわからないが、まあ、 $k$$\omega$ ($\gamma$) の関係を与 えてはいる。 $k<k_J$なら、ある実数 $\gamma$(正負どちらでも)があって 上を満たす。ただし、 $k=k_J$ の場合と $k=0$ の極限を除いては、値は流体 の場合とは一致しない。

というわけで、波長がジーンズ波長よりも長いモードは流体と同様不安定で、 勝手に成長することになる。

5 van Kampen mode

さて、それでは、音波に当たるような振動数が純実数のモードというものはあ るのだろうか?これについては、 BT を含めて標準的な教科書でも若干混乱し た記述がなされていることがある。以下、 van Kampen の論文 (1955, Physica, 21, 949) に沿った理解を試みることにする。分散関係の一つ前の式 に戻ってみると、

\begin{displaymath}
(-\omega + {\bf v}\cdot {\bf k}) f_a({\bf v}) + {4 \pi G \ov...
...\cdot {\partial f_0
\over \partial {\bf v}}\int f_ad{\bf v}= 0
\end{displaymath} (35)

となっている。面倒なのでまた $y, z $方向についてはあらかじめ積分すると、


\begin{displaymath}
(-\omega + k v_x) g(v_x) = - {8 \pi^2 G v_x \over k} f_0(v_x)\int gdv_x = 0
\end{displaymath} (36)

ただし、$g$$f_a$$y, z $ 方向に積分したものである。また、 $f_0$は 等方的であるとし、さらに
\begin{displaymath}
\int\int v^{-1}{d f_0 \over dv} dv_ydv_z = -2\pi f_0(v_x)
\end{displaymath} (37)

なる関係を使った。これは、$v_x$一定の平面上で極座標に変換すればすぐに 出てくる。上の式は $g$ についての線形斉次な方程式で、
\begin{displaymath}
g(v_x) = -{8 \pi^2 G \over k}{v_x f_0(v_x) \over \omega - kv_x}
\end{displaymath} (38)

という解を持つ。ここでは、 $\int g dv_x = 1$となるように規格化した。 ただし、意味のある解であるためには
\begin{displaymath}
-{8 \pi^2 G \over k}\int {v_x f_0(v_x) \over \omega - kv_x}dv_x = 1
\end{displaymath} (39)

を満たすようになっている必要がある。これは、最初に出した分散関係と実は 同じ式である。

最初にもいったように、この式がよろしくないのは $v = \omega/k$で積分が 特異になるからである。意味のある解を求める一つの考え方は、 $g$を超関数 に拡張してしまうことである:

\begin{displaymath}
g(v_x) = -{8 \pi^2 G v_x f_0(v_x)\over k}\left[{\cal P}{ 1\over \omega
- kv_x} + \lambda\delta(\omega - kv_x)\right]
\end{displaymath} (40)

ここで ${\cal P}$は主値をとるということで、要するに積分の不定性を $\delta$関数のほうに押しつけてみたというだけである。実際に、 規格化を満たすという条件から $\lambda$の値を決めることができる。この 「モード」は、任意の $\omega$$k>k_J$なる$k$ のすべての組合せに対し て存在する。これを van Kampen mode という。

これがいったいどういうものかを少し考えてみよう。まず、 $\delta$関数の 分は、 $v_x = \omega/k$ のところにだけ値があるということを示している。 つまり、位相速度が摂動を受けていないもともとの速度と等しい、いいかえれ ば与えたものがそのままラグランジュ的に動いていくようなものである。これ は、いま重力がまったくない極限を考えれば、単に摂動がまわりと相互作用す ることなくそのまま動いていくというものであると考えられる。これにたいし て、もう一つの項は重力による応答を示していると考えていい。

このへんはちょっとお話しになって申し訳ない(詳細を知りたいひとは van Kampen の原論文に当たってほしい)が、 van Kampen モードは完備であるこ とがわかっている。つまり、任意の摂動をvan Kampen mode の組合せとして表 現できる。

6 Phase Mixing

さて、任意の摂動が、減衰しないモードの組合せとして表現できるのなら、そ れはなにか音波のような伝わっていく波になっているのであろうか?実はそん なことはない。これは、以下の簡単な例で示すことができる。1次元で、重力 がない系である有界な領域に摂動(overdensity)を与えたとする。で、簡単の ために周期境界で、右から出ていったものが左から入ってくるとしよう。もち ろん、ちゃんと適当なポテンシャルを考えるとか、無限一様な場合を考えると かしても本質は同じである。(無限一様な場合は、単に周期的な摂動を与えた というのに対応する)

12cm \epsffile{note5fig1.eps}

すると、時間がたつに従って摂動が引き延ばされていくということがわかる。 このために、例えば密度の変化といった量は時間がたつにしたがって減衰して いく。

つまり、 singular な van Kampen mode 自体は減衰しないが、それを重ね合 わせた有限の広がりを持つ摂動は減衰するように見える。これが通常 phase mixing と呼ばれるものである。

7 Phase mixing による波の減衰

Phase mixing でどう波が減衰するのかを考えてみよう。

今、簡単のために重力がない一次元系を考える。この時 van Kampen mode は 単なる $\delta$関数なので、初期の波数が $k$であるような摂動は

\begin{displaymath}
f_1(x,v,t) = g(v)\exp[ik(x - vt)]
\end{displaymath} (41)

という形をしている。これは、自明な解になっているとい うことは式を良く見ればわかる。

しかし、注意して欲しいのは、これは最初に安定性解析の時に仮定した $f_a(v)\exp[i(kx - \omega t)]$という形とは違うということである。これは、 振動数 $\omega$$v$そのものであり、速度空間のなかでの位置に依存する ためである。

このことをいいかえると、上のような「自然な」解があるにもかかわらず、 モード解析をするとvan Kampen mode のような singular なものが出てくるの は、モード解析の仮定として位相速度が粒子の速度に寄らないようなものを考 えたからであるともいえる。

さて、上の「自然な」解はどのように振舞うかをちょっと見てみよう。 密度は、

\begin{displaymath}
\rho_1(x,t) = e^{ikx} \int g(v)e^{ivt}dv
\end{displaymath} (42)

これから、例えば $g(v)$
\begin{displaymath}
g(v) = \cases{
1/2v_0 & $\vert v\vert<v_0$\cr
0 & otherwize}
\end{displaymath} (43)

みたいなものだと、 $1/t$で減衰するという解がでてくる。

このような減衰が起きるのは、初期条件が非常に特別なものであるためである ということに注意してほしい。つまり、速度ごとに波の位相速度が違うのに、 初期条件としてその空間位相がすべてそろったものを考えたわけである。そう すると、時間がたてば位相はずれていくので速度方向に積分して見た波の振幅 は小さくなってくことになる。これは、速度方向の「波数」に時間が生で入っ てくるためである。

なお、 $g(v)$に適当な形を仮定すれば、もっと速く減衰するものも作れる。

もちろん、無衝突ボルツマン方程式にしたがった進化は可逆過程である(エン トロピーを生成しない)ので、原理的には逆に振幅が大きくなるような初期条 件も存在していないといけない。実際、減衰していく解で、どこかで時間反転 すればそういう解が作れるわけである。

なお、今日の話で要領を得ないと思った人は、、、この辺は基本的で大事な問 題であるわりにはあまり良いテキストも論文もないので、なかなか難しい。

この文書について...

重力多体系の進化

この文書はLaTeX2HTML 翻訳プログラム Version 2002-2-1 (1.71)

Copyright © 1993, 1994, 1995, 1996, Nikos Drakos, Computer Based Learning Unit, University of Leeds,
Copyright © 1997, 1998, 1999, Ross Moore, Mathematics Department, Macquarie University, Sydney.

日本語化したもの( 2002-2-1 (1.71) JA patch-1.9 版)

Copyright © 1998, 1999, Kenshi Muto, Debian Project.

Copyright © 2000, Jun Nishii, Project Vine.

Copyright © 2001, 2002, Shige TAKENO, Niigata Inst.Tech.

Copyright © 2002, KOBAYASHI R. Taizo, Project Vine.

を用いて生成されました。

コマンド行は以下の通りでした。:
latex2html -nomath_parsing -local_icons -show_section_numbers -split 0 note4-e.tex.

翻訳は Jun Makino によって 平成25年7月16日 に実行されました。


next_inactive up previous
上へ: Stellar Dynamics Lecture Note
Jun Makino 平成25年7月16日