next up previous
次へ: この文書について... 上へ: 理論天文学概論 戻る: 3 ビリアル定理

Subsections

4 ビリアル定理の応用

4.1 系の「比熱」

これまで何度か強調してきたように、無衝突ボルツマン方程式 CBE で記述さ れるシステムの力学平衡状態は熱平衡ではないし、自己重力系は熱平衡ではあ りえない。しかし、 King model のように近似的に熱平衡なものもあるので、 ここでエネルギーの出入りに対する系全体の応答、つまり比熱というものをちょっ と考えてみることにする。

ビリアル定理から $K = -E$ であったので、すぐにわかるように「エネルギー を奪うと運動エネルギーが増え、逆ならその逆になる」ということになってい る。つまり、見かけ上比熱が負になっている。

これは重力が効く系では普通のことで、例えば地球を回る人工衛星といったも のでも同様のことが起こっているわけである。

この、「見かけ上比熱が負」ということが、熱力学的不安定を通して構造形成 (自己組織化)が起きる基本的な理由である。

4.2 系の質量、 $M/L$

まず、系の「大きさ」についてどんなことがいえるかということをおさらいし ておく。$W$は全ポテンシャルであったので、定義により以下のように書ける:

\begin{displaymath}
W = -{1 \over 2} \sum_{i \ne j} {G m_i m_j \over \vert{\bf r}_i - {\bf r}_j\vert}
\end{displaymath} (25)

今、粒子の質量がすべて等しい場合(あるいは、分布関数が質量によらない場 合)を考えると、
\begin{displaymath}
r_v = \left< {1 \over \vert{\bf r}_i - {\bf r}_j\vert} \right>^{-1}
\end{displaymath} (26)

を導入して、
\begin{displaymath}
W = { -GM^2 \over 2r_v}
\end{displaymath} (27)

と書ける。この $r_v$ を普通ビリアル半径という。なお、 $r_g = 2 r_v$ の ことを gravitational radius 重力半径ということもあるが、「重力半径」と いうと時と場合によってはシュワルツシルト半径だったりすることもあるので 注意すること。とにかく、この $r_v$ を使って、$K = -W/2$を書き直せば
\begin{displaymath}
\vert v^2\vert = {GM \over 2r_v} = {GM \over r_g}
\end{displaymath} (28)

したがって、速度分散(系全体の平均)とビリアル半径がわかれば質量を決め られる。観測的には、球対称を仮定すれば、3軸方向の速度分散は等しいので、 視線方向の速度分散の3倍が3次元速度分散ということになる。問題はビリアル 半径のほうであるが、もちろん球対称(だけでは本当はだめで、さらに$M/L$ が半径によらずに一定であるという仮定も必要)を仮定すれば表面輝度分布か ら deprojection して計算できる。

大雑把な見積りでよければ、例えばhalf mass radius (質量の半分が入って いる半径)$r_h$ を適当に見積もって、それで $r_v$ の代わりにしてもそれほど大 きな誤差はない。典型的には

\begin{displaymath}
r_h \simeq 0.8 r_v,
\end{displaymath} (29)

程度である。具体例をあげれば、プラマーモデルでは
\begin{displaymath}
r_h = {3\pi \over 16 \sqrt{2^{2/3} - 1}} r_v = 0.77r_v,
\end{displaymath} (30)

ハーンキストモデルでは
\begin{displaymath}
r_h = {(1 + \sqrt2) \over 3}r_v,
\end{displaymath} (31)

であり、仮想的な密度一様の球というものを考えると
\begin{displaymath}
r_h = { 3\times 2^{2/3}\over 5} r_v = 0.95 r_v.
\end{displaymath} (32)

というようなことになる。実際に観測から直接決められるのは、多くの場合 $r_h$ ではなく有効半径 effective radius $r_e$ であり、これは質量ではな くて投影した輝度の半分が入っている半径である。一般に $r_e < r_h$ であ ることはいうまでもないが、まあ、それほど大きな差にはならない。

4.3 回転系の偏平度(あるいは楕円銀河の回転)

今、$z$軸まわりに対称な、回転している銀河を考えると、ポテンシャルテン ソルについては対称性から

\begin{displaymath}
W_{xx} = W_{yy}; \quad W_{ij} = 0 \quad (i\ne j)
\end{displaymath} (33)

で、運動エネルギについても同様にクロスタームは消えるので、結局テンサー ビリアル定理で残るのは
\begin{displaymath}
2 T_{xx} + \Pi_{xx} + W_{xx} = 0;\quad
2 T_{zz} + \Pi_{zz} + W_{zz} = 0
\end{displaymath} (34)

の2つだけである。この2つから、
\begin{displaymath}
{ 2 T_{xx} + \Pi_{xx} \over 2 T_{zz} + \Pi_{zz}} =
{ W_{xx} \over W_{zz}}
\end{displaymath} (35)

となる。今、子午面環流みたいなのはないとすれば、 $T_{zz} = 0$ で、
\begin{displaymath}
2 T_{xx} = {1 \over 2} \int \rho v_{\phi}^2 d^3{\bf x}= {1 \over 2} M
v_0^2
\end{displaymath} (36)

(これは $v_0$、つまり平均の回転速度の定義と思って下さい)。同様に
\begin{displaymath}
\Pi_{xx} = M\sigma_0^2
\end{displaymath} (37)

これは視線方向(横から見ているとして)のランダム速度の平均 (r.m.s)であ る。あと、 $z$方向についてはなにかいえないかというと、 $x$方向との速度 分散の異方性パラメータ $\delta$ を導入して
\begin{displaymath}
\Pi_{zz} = (1-\delta)\Pi{xx}
\end{displaymath} (38)

と書いておくことにする(これは、球対称な系の速度分散の異方性 とはまったく違ったものであることに注意)。これらをまとめると、
\begin{displaymath}
{v_0^2 \over \sigma_0^2} = 2(1-\delta){W_{xx} \over W_{zz}} - 2
\end{displaymath} (39)

ということになる。

さて、 $W_{xx}/W_{zz}$と、銀河の「形」の関係についてであるが、一般に密 度が以下のように書ければ

\begin{displaymath}
\rho = \rho(m); \quad m^2 = a_1^2\sum_{i=1}^3 {x_i^2 \over a_i^2}
\end{displaymath} (40)

言い換えれば、密度が軸比一定の楕円体の表面上で一定なら、 $W_{ii}/W_{jj}$$a$だけの関数で $\rho(m)$の関数形によらないというこ とがわかっている。実際に書くとかなり面倒な形をしているので、具体的な形 は BT をみてもらうことにして、結果だけを使うと、要するに速度の非等方性 $\delta$ と偏平度 $\varepsilon 1 - b/a$を決めれば $v_0 \over \sigma_0$ が決まるということになる。どんな感じかを BT の図から示すと
10cm \epsffile{btfig4-5.ps}
10cm \epsffile{btfig4-6.ps}

まあこんな感じになる。上の図は、 $\varepsilon$の関数として $v_0 \over \sigma_0$ がどう変化するかを、 $\delta$ のいくつかの値について示したも のである。破線は、真横でない方から見たらどういうふうになるかを、いくつ かの場合について示したものである。下の図は、実際の観測結果と $\delta =
0$ の線を重ねて書いてみたものである。暗い楕円銀河は $\delta =
0$ の線 に近いところに集まっているのに対し、明るいものは可能な範囲全体にあり、 どちらかというと下の方(回転が小さく、おそらく異方性が大きいほう)に集 まっているということが見てとれる。このような構造の違いは、それぞれの銀 河がどうやってできたかということになにか重要な手がかりを与えているはず である。

とはいうものの、どうやってできたかということに対する標準的な理解という ものがあるわけでは必ずしもない。ここでは、せっかくなので一つの考え方と して、合体説(楕円銀河のうち、特に大きいものは銀河同士の合体でできた) をとった場合にはどのような説明が可能かということを見てみる。これ以外の 解釈がないというわけではない。

=10cm \epsffile{Okumura1991fig7.ps}

これは Okumura et al. (1991, PASJ 43, 781) によるシミュレーショ ンの結果と上のグラフを重ねたものである。シミュレーションはすべて2つの 等方的で回転を持たないプラマーモデルを2つ放物軌道でぶつけたもので、番 号順に最初の軌道の近点距離が大きくなる( 4 で $r_p/r_h = 2.6$である)。

大雑把にいうと、番号の大きいものほど軌道角運動量が大きいので、 $v_0/\sigma_0$が大きくなる。もっとも、ある程度より大きくなっていないこ とに注意してほしい。これは要するに合体するためにはかなりの軌道角運動量 を捨てないといけないので、残った系が持つ角運動量には上限があるためであ る。

合体後の銀河は特別な場合(1,7)を除いて3軸不等なので、その3方向からみ てグラフ上に点をとり、それを結んで3角形を書いている。大雑把にいって、 適当な方向からみたときにはこの3角形の中にあると思っていい(実際にある 方向から見た時にどう見えるかというのはもうちょっとややこしいが、いまは そのあたりの議論は省く)。

シミュレーションの結果からわかることは、仮に楕円銀河が合体で出来たとす れば、暗いものは初期に近点距離の大きな軌道(おそらく楕円軌道)をもち、 明るいものは近点距離の小さな軌道から合体したものが比較的多いということ になろう。


next up previous
次へ: この文書について... 上へ: 理論天文学概論 戻る: 3 ビリアル定理
Jun Makino 平成20年6月11日