next up previous contents
Next: 第9章 熱力学的進化 Up: 理論天体物理学特論I Previous: 第7章 Violent relaxation (続き)、2体緩和

第8章 2体緩和のタイムスケール、フォッカー・プランク方程式

前回は、実際にバックグラウンドの粒子も動いている場合について、2体緩和 によって粒子の速度、エネルギーがどう変化するかの期待値(のモーメント) を求めた。今日は、前回の結果を使って、まずいくつかの重要な場合について 2体緩和がどのように働くかを、数値実験と合わせて議論する。次に、上の結 果から分布関数の変化に対する方程式を導くことを考える。

8.1 2体緩和のタイムスケール

「タイムスケール」というのは、普通ある量 x の時間変化が

なる形で書ける時の T のことである。もちろん、一般にはTxを含む いろんなものの関数である。

そういった意味で考えやすいのは、温度(平均運動エネルギ)が違う2つの空間一様な 分布が重なりあっている時に、どのようにして2つが近付いていくかというの ものである。これは、もちろん時間が立てば熱平衡に近付くわけである。以下、 実際に計算してみる。

8.1.1 等分配のタイムスケール:理論

今、フィールドに質量 の粒子が一様に分布しており、テスト粒子として質量 のものがこれもまた一様に分布しているとする。さらに、どちらも速度 分布はマックスウェルで与えられるとする。ここでは等分配を考えるので、そ れぞれの粒子1個当たりのエネルギーを , と書く。今、テスト粒 子のエネルギー変化の平均を考えると、

と書けることになる。これに前回求めた を入れて実際に積分を 実行することができて、結果は

となる。

ここで、いくつかの極限的な場合を考えておくことは有益であろう。まず、 という状況を考えてみる。これはつまり非常 に重いものと軽いものが、同じような空間分布、速度分布で広がっている場合 である。この時は上の式で なので、

となる。なお、この時、変化率はテスト粒子の速度に分母の を通し てしか依存しないので、 の極限でエネルギー変化(減 速)のタイムスケールは一定値にいき、それは の時の値とそ れほど違わない。

次に、という状況を考えてみる。この時は上 の式を

となる。 ここで

を使った。 さらに、 などを使って書き直せば

を得る。つまり、速度が小さい極限では、一定の率でエネルギーをもらうわけ である。言いかえれば、温度が2倍になるタイムスケールというものは、温度 に比例して小さくなるともいえる。

さて、通常「2体緩和のタイムスケール」という時には何を指してい るかというと、この等分配のタイムスケールのことではないのが普通である。 が、時と場合によっていろんなものが出てくるが、まあ同じようなものである。 普通に使われるのは、

とするものである。 ここで はr.m.s. 速度である。になにか意味があるわけではなく、 こう定義したというだけである。

これは、ローカルな量で定義されていて、例えば系全体の緩和時間といったも のを考えるのにはちょっと不便なこともある。というわけで、いわゆる half-mass relaxation time というものを導入しておく。これは、 半径 の中に質量の 1/2があるとして、その中の密度は一様であるとし、 また前にやったようにビリアル定理からといった関係 を使えばでてくる。これは

となる。

ここで注意しないといけないことは、 はあくまでも球対称に近い系の half mass radius のあたりでの緩和時間であるに過ぎないということである。従って、 球状星団全体の緩和時間とか、あるいは楕円銀河、銀河団といったものには有 効な概念であるが、球対称から大きくずれた銀河とか、あるいは half mass radiusのずっと外側、ずっと内側では全く違ったものになっていることに注意 する必要がある。さらに、速度分布が非等方であるとか、回転がメインである とかでもまた話が全く変わってくる。このような場合、ローカルな緩 和時間、あるいはエネルギー変化自体の式に戻って考えないと、タイムスケー ルについて全く間違った推定をしてしまうことになる。

これがもっともクリティカルに効いて来るのは、実際の天体においてというよ りはむしろシミュレーションにおいてである。これについては後で具体例で議 論しよう。

8.1.2 等分配のタイムスケール:実験

さて、実際に数値実験でエネルギー等分配に近付く過程を見て、これまで理論 的に考えてきたもの、特に積分の上限がどうなっているかみようというわけだ が、これはそれほど簡単ではない。というのは、「空間分布が無限一様」とか 「等温」とかいう初期条件が設定できないためである。

一例として、 Farouki and Salpeter (1982 APJ 253, 512, 1994 APJ 427, 676)の実験を取り上げてみる。彼らは、実際に2種類の違う質量の粒子からなる、 初期に力学平衡にある系を考えた。初期には単位質量あたりの運動エネルギー が同じであるようにした(つまり、空間分布、速度分布ともに同じということ)。 したがって、重い方が熱平衡に比べて余計に運動エネルギーを持っていること になる。

正確にいうと、実際に力学平衡にある系を作ったわけではなく、適当に球とか 立方体のなかに一様に粒子をばらまいて、落ちつくまで待ってから使うという方法を とっている。速度は、全体としてビリアル平衡になるようにマックスウェル分 布を与えた。

さて、この系を進化させると何が起こるかをまず考えてみる。重い方の温度が 高いので、平均的には重い方から軽い方にエネルギーが流れるであろう。しか し、その結果温度が近付くであろうか?これはとりあえず何ともいえない。と いうのは、ポテンシャルが一様ではないのでエネルギーの変化は粒子の速度分 布だけでなく密度分布、すなわち粒子の位置と重力エネルギにも影響するから である。エネルギーの変化が温度にどう影響するかは、もともとの粒子の分布 にも依存するわけである。

さらにややこしいことに、 がなんであるかという問題もある。 通常、N体計算においては、2粒子間の重力ポテンシャル をそのまま使わないで、典型的には

という形に「ソフトニング」したものを使う。これは、数値計算上の困難を避 けるのが第一義的な理由である。つまり、純粋な ポテンシャルを使うと、 2つの粒子が無限に近付き得る。無限に近付くこと自体は問題ではないが、こ の時速度と重力エネルギーも無限大になる。従って、非常に高い精度で計算し ても、エネルギーの誤差が非常に大きくなってしまう。

もちろん、近付いたら解析解に切替える、あるいは摂動法にする、または座標 変換をして特異性を消すといった方法があるが、どれもかなり面倒である。そ こで、多くの数値実験で、ポテンシャルをすこし変えて原点で発散しないよう にする。この場合、インパクトパラメータが よりも小さい時には ほとんど曲がらなくなるので、速度変化の積分の下のほうを修正する必要が起 きる。大雑把にいって の寄与を無視すればよいことになる。

次のページの図に典型的な結果を示す。Farouki and Salpeter では、温度と かエネルギー変化を直接測定するのは断念し、その代わりに軽い粒子の half mass radius をとった。彼らは緩和時間として

つまり、半径の変化のタイムスケールをもって緩和時間であるとした。これは もちろん定数を別にすれば上で考えた定義と一致するものになっている。

この結果から、彼らは と粒子数 N について以下のような 関係を得た

上と下は の大きさが違い、上はの程度、下はの 程度である。ここでFは緩和時間を

と書くことにしたために出てくる定数である。26という定数は、前に出てきた を使って書き直したので、定数の が違う値になっ ているだけのものである。

彼らの主な主張は、「この振舞いは上のカットオフがシステムサイズであると いう理論と一致している」というものである。それはまあかなり確かといって もよいと思われる。が、 Fの値を推定すると 10 程度になり、かなり大き過 ぎるものになっている。これは、定義の違いもあるので評価はちょっと難しい。

もうちょっと良くあっている例を紹介しておこう。

これは Spurzem and Takahashi (1995, MNRAS 272, 772) によるものである。 ここでは、緩和時間といった定義のはっきりしないものを比べているのではな く、系の構造そのものの変化を、N体計算の結果とモーメント方程式から導 いた(球対称、等方の)フォッカー・プランク方程式の数値解とで比べている。 結果の一致は素晴らしいものである。

8.1.3 応用(問題)

さて、2体緩和についていろいろやってきたので、ちょっと手を動かして計算 してみるのもいいかと思う。以下に問題を出す。レポート課題の一つと思って もらってよい。

(2-1a)Singular isothermal な分布をしたダークマターハローのなかに、 自己重力的でない無限に薄い星のディスクがあるとする。ダークマター粒子の 質量が無限小ならばディスクは薄いままだが、数値計算では2体緩和のため厚 くなる。厚さが2倍になるタイムスケールを、厚さ、中心からの距離、その内 側のダークマターハローの質量、ハロー構成粒子の質量の関数として求めよ。 なお、結果は回転周期に対する比で表すこと。

(2-1b)上で求めたものから、それらしいパラメータを使って太陽がある 当たりでディスクが厚くなる様子の近似的な解析解を求めよ。ただし、時間は ハロー構成粒子の質量でスケールできることに注意せよ。

(2-1c)自己重力が無視できない薄いディスクについて、ディスク内の星 同士の2体緩和により厚さが2倍になるタイムスケールを上と同様に求めよ。 ただし、簡単のためダークハローがあって自己重力の効果は小さいとする。ま た、ダークハローに由来する2体緩和はないものとする。

(2-2a)SPH法とは、流体を大きさをもった粒子の重ね合わせで表現する 計算法である。無衝突系であるダークマターとホットガスが力学平衡にある状 態を、SPH法とN体を組み合わせて計算した場合、2体緩和の効果はどのよう に働くと考えられるか、例をあげて定量的に議論せよ。

(2-2b)さらに、適当な最近のSPHとN体を組み合わせて行なった計算 の論文を読んで、2体緩和の効果を評価し、それが結論に影響しているかどう かを考察せよ。

8.2 フォッカープランク近似

さて、数値計算がどうたらという話を別にすれば、2体緩和を考える理由はそ れにより系がどう進化するかを理解するということである。そのためには、粒 子の速度変化のモーメントの式から、分布関数の変化についての方程式を導く ことが有用であろう。

粒子の物理量の変化の1次と2次のモーメントから(高次の項の寄与を無視して) 分布関数に対する移流拡散方程式を導くことができる。この操作を通常フォッ カープランク近似といい、でてきた方程式をフォッ カープランク方程式という。以下、その導出を行なってみる。 これは、ボルツマン方程式の衝突項を求めることに対応する。

今、 が、速度 の粒子が、時間の間に速度変化 (の近傍の の範囲)を受ける確率であるとする。 すると、 たった後の分布関数は以下のように書ける

これは のところにやってくる可能性があるものすべてについての和をとっ ただけである。ここで、両辺を展開することを考える。左辺は

である。右辺は

ここで、

とおけることを使えば、微分と積分の順序を入れ換えて

これで、理屈の上では分布関数の変化が計算できるということになる。もちろ ん、実際にこれを解いて自己重力系の進化を調べるのは、必ずしも容易ではな い。その理由は、分布関数が6次元位相空間上で定義されること、それがジー ンズの定理を満たすように進化しなければいけないことである。具体的には、

というような困難があり、従来は、分布関数をエネルギーだけの関数と近似す る計算しか行なわれていなかった。1994年頃に、Takahashi が初めて の場合に信頼できる結果を得ることに成功した。しかし、これも、拡 散係数を求める時に、フィールドの分布は等方的であるという近似を行なって いる。

なお、モンテカルロ法など、もうちょっといい加減な方法ではある程度のこ とは出来ている。



Jun Makino
Mon Jun 1 23:17:40 JST 1998