next up previous
Next: 2 恒星系でのジーンズ不安定とランダウ減衰 Up: 理論天体物理学特論I Previous: 理論天体物理学特論I

Subsections


1 ジーンズ不安定

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

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

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$ の程度になることがわかる。

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


next up previous
Next: 2 恒星系でのジーンズ不安定とランダウ減衰 Up: 理論天体物理学特論I Previous: 理論天体物理学特論I
Jun Makino
2003/11/16