next_inactive up previous
上へ: Stellar Dynamics Lecture Note

pt 称

重力多体系の進化

牧野淳一郎

1 球対称モデル(続き)

1.1 King Model

等温モデルは、すでに述べたように熱平衡(エントロピーの変分が 0 )とい う重要な意味を持つ定常解ではあるが、なにしろ質量が無限大であり現実に存 在しないのでちょっと困るところがある。なにか適当な仮定を置くことで、 「おおむね等温モデルであり、なおかつ有限の大きさをもつ」というものを考 えることはできないだろうか?


\begin{displaymath}
f({\cal E}) = {\rho_1 \over (2\pi \sigma^2)^{3/2}} e^{{\cal ...
...\sigma^2)^{3/2}} \exp\left({\Psi -
v^2/2\over \sigma^2}\right)
\end{displaymath} (1)

上の分布関数で、質量が発散する理由は何かを思い出してみよう。その本質的 な理由は、分布関数がエネルギー無限大 ( ${\cal E}\rightarrow - \infty$)ま で0にならないことにある。

有限の質量のものが自己重力でまとまっているためには、すべての粒子のエネ ルギーが負でないといけないので、これでは自己重力系が表現出来ないのはあ る意味では当然のことといえる。

それならば、ある有限のエネルギー以上のものはないことにしてしまえばいい。 そのやり方にはいろいろあり得るが、とりあえず Lowered Maxwellian と呼ば れる以下のようなものを考える


\begin{displaymath}
f({\cal E}) = \cases{ {\rho_1 \over (2\pi \sigma^2)^{3/2}} (...
...E}
/\sigma^2}-1) & $({\cal E}> 0)$\cr
0 & $({\cal E}\le 0)$
}
\end{displaymath} (2)

これは${\cal E}= 0 $$f=0$ となるように、 1 を引いたというだけである。 これしか方法がないというわけではないが、これは扱いやすいこともあってもっ ともよく使われている。

これはいかにも人工的な感じがすると思うが、等温モデルから有限質量でなお かつ有限半径のモデルにするもっとも簡単な方法なので、まあ、そういうもの と思って欲しい。

例によって、まず速度空間で積分すれば

$\displaystyle \rho$ $\textstyle =$ $\displaystyle {4\pi \rho_1 \over (2\pi \sigma^2)^{3/2}}\int_0^{\sqrt{2\Psi}}
\left[\exp\left({\Psi - v^2/2 \over \sigma^2}\right) -1
\right]v^2dv$  
  $\textstyle =$ $\displaystyle \rho_1 \left[e^{\Psi/\sigma^2} {\rm erf}\left({\sqrt{\Psi \over
\...
...\sqrt{4\Psi \over \pi \sigma^2}\left(1 + {2\Psi
\over 3 \sigma^2}\right)\right]$ (3)

ここで ${\rm erf}$は誤差関数で、積分が有限区間であるために出てくる。最後の項 は1引いている分の寄与である。ちなみに

\begin{displaymath}
{\rm erf}(z) = {2 \over \sqrt{\pi}}\int_0^z e^{-t^2}dt.
\end{displaymath} (4)

これでポテンシャルの関数として密度が求まったので、あとはポアソン方程式 に入れて数値的に解くだけである。ただし、 King model の場合境界条件につ いてすこしきちんと考える必要がある。

半径方向の分布は、中心から無限遠まで与えられるわけであるが、実は外側の 境界をどうとるべきかはちょっと自明ではないので、とりあえず中心から初期 値問題として解くことを考える。

初期条件としては、まず $d\Psi/dr = 0$ とする、すなわち、中心密度が有限 の解を考える。 $\Psi_0$ は任意に選べるので、これの値によっていろいろな 解がでてくる。

=10cm \epsffile{kings.ps}

これは実際に数値的に解いてみたものの例である。速く落ちるものから、 $\Psi_0$ が 1, 3, 6, 9, 12 と変えてみてある。

なお、横軸のスケールの $r_0$は、

\begin{displaymath}
r_0 = \sqrt{9\sigma^2 \over 4\pi G \rho_0}
\end{displaymath} (5)

として無次元化するのに使っている。これは、いわゆる「コア半径」というの とそこそこ一致するということになっている。通常、キングモデルのコア半径 というときにはこれをさす。観測的には、中心の表面輝度の 1/2 になるとこ ろとするのが普通である。

グラフからわかるように、有限の半径 $r_t$$\rho$は 0 になる。これは、 解いていったときに $\Psi$ が 0 になってしまうためである。この半径のこ とを King model の tidal radius 潮汐半径という。このモデルのばあい、 $\Psi$と本当のポテンシャル $\Phi$ の間に以下のような簡単な関係が成り立 つことに注意。

\begin{displaymath}
\Phi = -{GM \over r_t} - \Psi
\end{displaymath} (6)

ここで $M$は系の全質量である。

King Model は、球状星団のプロファイルのモデルとして非常によく使われて いる。なお、 $c = \log(r_t/r_0)$ のことを concentration parameter といっ て、観測データにキングモデルを合わせた論文では普通これがパラメータにな る。理論計算では $\Psi_0$ が使われるので、ちょっとややこしいことが多い。

1.2 多成分系の例

さて、今まで、単に無衝突ボルツマン方程式の球対称定常解というものの例を いろいろ見てきたわけだが、ちょっと飽きてきたので目先を変えてみよう。具 体的には、「球対称の恒星系のなかでのガスの分布」というものを考えてみる ことにする。というのは、現在でも観測データの解釈に使われている簡単なモ デルは、現在までにやった範囲で十分カバーできるものになっているからであ る。

多くの銀河団には高温ガスが存在しており、X線で観測可能である。

=10cm \epsffile{david1990.ps}

上は Einstein IPC による Hydra A 銀河団に(David et al. ApJ 1990, 356, 32) の像である。あと、 Chandra やすざくのとか見たことがある人も多いと思 うが、まあだいたいこんな感じにきれいに丸いのが普通である。 X 線ガスの 振舞いについての詳しい話は Sarazin の教科書かなにかを見てもらうとして、 とりあえず重要なのは、

ということである。この2つの仮定をつかっていろいろしらべていくことにし よう。もちろん、Chandra やすざくでもっと細かい構造があることはわかって きているが、大局的な理解は基本的に変わっていない。

まず、話を簡単にするために、

ということにする。これは、銀河団なんかだとちょっと怪しいが、まあ全然だ めというほど悪くはない。さらに、 とする。だんだん、ほんとかなあという気がしてくることと思うが、まあ、研 究の現場というのはそういうものです。

さて、これぐらいたくさん仮定すると、美しい理論を構築できることになる。 まず、ガスの速度分散と恒星系の速度分散が等しい場合というのを考えてみる。 ガスの密度$\rho_g$が従う方程式は、星の密度を $\rho_s$ として

\begin{displaymath}
{d \over dr} \left(r^2 {d \log \rho_{g} \over dr}\right) =
-{ 4\pi G m \over k_BT}\rho_s r^2
\end{displaymath} (7)

と書けるわけだが、$\rho_s$ 自体がこの方程式の解になっている。すると、 これは $\rho_g$の対数についての微分方程式なので、
\begin{displaymath}
\rho_g = C \rho_s
\end{displaymath} (8)

というのが解になっている。中心で微分が 0 とすれば他の解はない。つまり、 どこでもガスと星の密度比が等しいというのが平衡状態になる。この時は、も ちろん、ガスの自己重力を考慮しても平衡状態の分布の形は変わらな い。

では、温度(速度分散)が違う場合はどうなるかが問題である。ガスの温度 $T_g$と、恒星系に対応する「温度」 $T_s$の比を $T_g/T_s = 1/\beta$と置 くと、 ガスの密度が従う方程式は、

\begin{displaymath}
{d \over dr} \left(r^2 {d \log \rho_{g} \over dr}\right) =
-\beta{ 4\pi G m \over k_BT}\rho_s r^2
\end{displaymath} (9)

と書けるわけだが、この解が

\begin{displaymath}
\rho_g = C \rho_s^{\beta}
\end{displaymath} (10)

are? で与えられることは「容易にわかる」。すなわち、対数でグラフに書いて縦に $\beta $倍したものが、 解の形を表しているわけである。一般に、温度が高ければ等温解に比べて傾き が緩くなり、低ければより急になるわけである。

1.2.1 $\beta $モデル

と、ここまではちゃんと正しい話をしてきたが、実際に観測結果の解釈に使う 時にはさらに大胆な近似がなされる。以下、 Sarazin の X-ray emissions from clusters of galaxies 5.5節に従っていわゆる$\beta $モデルについて説 明し、その問題点について検討する。

式(10)はもちろん正しいし、有用なものではあるが、あま り使い回しがよくない。というのは、$\rho_s$のほうがそもそも数値的にしか 求まっていないからである。というわけで、$\rho_s$をなにか解析関数で置き 換えることを考える。さて、以下の質量分布

\begin{displaymath}
\rho = \rho_0 [1+ (r/r_c)^2]^{-3/2}
\end{displaymath} (11)

を考えてみると、これは $r<2.5r_c$程度ならば等温解の非常によい近似になっ ている。(遠くにいくとどんどんずれるが)これを恒星系の分布と思えば、ガ スの分布は
\begin{displaymath}
\rho = \rho_0 [1+ (r/r_c)^2]^{-3\beta/2}
\end{displaymath} (12)

と解析関数で書けることになって都合がよろしい。さらに、 X 線で観測する ということを考えると、実際に観測されるのは X 線の強度で、これは密度の2 乗を視線方向で積分したものと思っていい。つまり
\begin{displaymath}
I(r) = 2\int_0^{\infty} \rho(\sqrt{r^2+z^2})^2dz
\end{displaymath} (13)

となるわけだが、これは積分できて、
\begin{displaymath}
I(r) = C [1 + (r/r_c)^2]^{-3\beta + 1/2}
\end{displaymath} (14)

となる。定数はいろいろ出てくるけど省略。

左はうまくいく例である (Jones and Forman 1984)。右は先ほど見せた Hydra A

=6cm \epsffile{Jones1984a2256.ps} =6cm \epsffile{david1990fig2.ps}

これらはそれなりに典型的なものだが、もうちょっと最近の論文(Ikebe et al. 1997, ApJ 481, 660)をあげておくと

=8cm \epsffile{Ikebe1997HydraA.ps}
とまあ、こんな感じになる。

いくつかの観測結果をこのように $\beta $モデルでフィットした結果から、以 下のような問題があるということになっている。

これが問題であるかどうかについては、例えば以下のようなことも考慮する必 要があろう。

この文書について...

重力多体系の進化

この文書は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 note2-e.tex.

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


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