next up previous
Next: 10 来週の予定 Up: 理論天体物理学特論I Previous: 8 Phase Mixing

9 Landau Damping

さて、モードがあって、世の中はその重ね合わせであるというのなら、それで 話はおしまいではないかとおもうのが人情だが、普通はそういう話にはなって いない。というのは、 Landau Damping という難しいものがあるということに なっているからである。これはどういう原理ででてくるかというと、要するに 複素数の $\omega$ があると信じて、そういう解を求めると、ちゃんとそうい うものが見つかるというものである。つまり、摂動の速度方向の分布 $f_a$を うまくとってやると指数関数的に減衰するモードがでてくる。ただし、注意し て欲しいのは、これは、 $f_a$に制限をつけないと出てこないということであ る。不連続な $f_a$を仮定すれば、時間のベキでしか減衰しないような解を構 成することも出来る。ただし、そういった解はモードの形、つまり位置、時刻 の指数関数の形に単純に書けるとは限らない。

逆にいえば、解が指数関数の形に書けると仮定すれば、ジーンズ波長より短け ればそれは指数関数的に減衰するわけである。

今、簡単のために重力がない一次元系を考える。この時 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)$に適当な形を仮定すれば、もっと速く減衰するものも作れる。

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

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



Jun Makino
2003/11/16