next up previous
Next: 5 臨界点 Up: 理論天体物理学特論I Previous: 3 CBEの線形化

4 Jeans 不安定

上の線形化されたボルツマン方程式は、平衡解の回りならば何でも使えるが、 ここではもっとも簡単に解析できる場合として空間分布が一様な場合を考える。 これから $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$の場合と、不安定、すなわち振動 数が純虚数の場合を考えよう。



Jun Makino
2003/11/16