1 ブラックホールのある系

1.1 中心ブラックホールのある星団の構造と進化

中心にブラックホール(というか、なんか重いもの)がある星団ではどの ような構造が見られることになるだろうか？これに理論的に答えたのは Bahcall and Wolf (1978) である。この頃には球状星団の中心にブラッ クホールがある可能性がかなり高いと考えられていたため、そのような方向の 研究が盛んであった。しかし、 球状星団の X 線源がほとんど Low-mass X-ray binary、つまり中性子星と小さな主系列星の連星であり、また必ずしも 星団中心にあるわけでもないということが 1980年代になって明らかになった ため、しばらくこの方向の研究は止まっていた。

それはともかく、Bahcall and Wolf はフォッカープランク方程式を数値的に 解くことで密度構造を決めたが、その結果は解析的に理解できることがわかっ ている。その考えは以下のようなものである。

中心部分の、ブラックホールの重力が支配的な領域を考え、また簡単のために 分布関数は等方的であるとする。 速度分散はポテンシャルで決まるので、ケ プラー速度になって速度は $v \propto r^{-1/2}$ になる。密度が $\rho$ であるとしよう。

中心に向かって温度があがっているので、熱は中心から外側に向かって流れる。 ここで、定常状態ならば熱流 $L$ が半径に依存しない。

大雑把にいうと、ある半径での熱流は、そこでの緩和時間くらいの間にその領 域の全エネルギーぐらいが流れ出すと考えることで見積もることができる。こ れはなんか根拠がない仮定であると思うかもしれないが、仮に密度が半径のべ きであるとすれば、 温度は元々半径のべきなので無次元量としての温度勾配 の大きさはどこでも同じになるため、この仮定は正しいことになる。

問題は、ではそういうべき乗の解はあるかどうかということだが、緩和時間は $t_r \sim v^3/\rho$ の程度、全エネルギーは $T = Mv^2 \sim \rho r^3 v^2$ の程度なので、 $T/t_r = 一定$ と置くことで

$$\rho \sim r^{-7/4}$$ (1)

という関係がでてくる。

こんな大雑把な計算でいいのかと思うであろうが、割合うまく数値計算の結果 を説明できている。

1.2 熱的進化以外の場合

ここまでは中心にブラックホールがある系について、周りにあるのが同一の質 量の星の集団であり、熱力学的な定常状態になれば密度が半径の $-7/4$ 乗の べき分布になることを導いた。これは理論的には美しいが、必ずしも非常に現 実的なケースとはいいがたい。以下では、より現実的と考えられるいくつかの 場合について、分布がどのようになるべきかを考えてみる。具体的には、以下 の 3 ケースを考える

1. 中心ブラックホールが断熱成長する場合

2. 力学的な時間スケールで「突然」中心ブラックホールができる場合

3. 質量分布がある系の熱力学的な進化

これらはそれぞれ、対応する現実的な系がある(かもしれない)と考えられる。

1.3 中心ブラックホールが断熱成長する場合

これは、例えばガス降着などでブラックホールが比較的ゆっくり成長する場合 に、周りの恒星集団の分布がどう変わるかという話である。ゆっくりといって も、力学的な時間スケールよりは十分に遅いが2体緩和の時間スケールよりは 速いものを考える。これは、銀河中心の巨大ブラックホールの場合にはありそ うな話である。

QSO や AGN のcentral engine は巨大ブラックホールへのガス降着であると考 えられているので、ガス降着が終わったあとの恒星系の分布は、この、ブラッ クホールが断熱成長した場合で与えられると考えられるであろう。この場合の 分布関数の変化を数値的および解析的に調べたのは Young (1980, ApJ, 242, 1232) である。以下、彼の論文の議論を要約しよう。

始めはブラックホールがなかったとして、分布関数が $f(E,J)$ であるとする。 考えないといけないことは以下の2つである。

• 最初に $(E,J)$ にいた星はどこに移動するか

• その結果 $f(E,J)$ はどう変化するか

というわけで、順番に考えていこう。

まず、$(E,J)$ にいた星がどこに移動するかであるが、 ポテンシャルは球 対称のままなのでその形が変化しても角運動量 $J$ は保存する。従って、 $E$ の変化だけを考える。 ポテンシャルの変化はゆっくりであるとしたので、 断熱不変量がある、具体的には、 radial action

$$I_R = \oint v_r dr = 2\int_{r_-}^{r_+} [2(E-\phi)-(J/r)^2]^{1/2}dr$$ (2)

が保存することになる。

初めに恒星系は有限サイズのコアを持っていたとしよう。この時、コアの十分 内側では、 ポテンシャルは中心密度を $\rho_0$ として

$$\phi = \frac{2 \pi}{3}\rho_0 r^2$$ (3)

で与えられ、断熱不変量は

$$I_R = \oint v_r dr \propto \sqrt{\frac{3}{\rho_0 \pi}}E-J$$ (4)

(比例係数は無視)となる。 さて、ブラックホールが十分に成長した後、元々コアの中心近くにいた星はブ ラックホールのポテンシャルの深いところにいると考えよう。この時には、新 しいポテンシャルは $-M_{BH}/r$ であり、断熱不変量は

$$I_R^* = 2\pi(-J+M_BH\sqrt{-2E^*})$$ (5)

となる。これを解けば $E^*$ は求められる。但し、ここでは BH ができてか らの量に $*$ を付けて区別することにした。

さて、問題は、 $f^*(E^*, J)$ がどうなるかであり、これがわかれば密度分 布がわかる。少しややこしいのは、最初のエネルギーが同じであっても最終の エネルギーは角運動量 $J$ によって違うことで、このために算数が少し面倒 になる。

分布関数 $f$ ではなく、 $(E,J)$ 空間での分布関数 $N(E,J)$ を考えると、 ブラックホールの成長によってこれは滑らかに射影されるので

$$N^*(E^*,J)dE^*dJ = N(E,J)dEdJ$$ (6)

なる関係が成り立つ。$J$ は同じなので、これは

$$N^*(E^*,J)dE^* = N(E,J)dE$$ (7)

ということである。従って、 $dE/dE^*$ が計算出来ればいい。ここでは、 $I_R$ と $E$ の関係を使ってみよう。つまり、$I_R$ を $E$ で偏微分すると

$$\left.\frac{\partial I_R}{\partial E}\right\vert _J = P(E,J)$$ (8)

となることが知られている。但し、ここで $P(E,J)$ は半径方向の周期である。 (計算は簡単であるのでやってみること) 従って、

$$\left.\frac{dE^*}{dE}\right\vert _J = \frac{P(E,J)}{P(E^*,J)}$$ (9)

となってだいぶ目標に近づいてきた。

後は $f$ と $N$ の関係だが、これは単に

$$N(E,J) = 8\pi^2 Jf(E,J)P(E,J)$$ (10)

であるということが知られている。で、結局これらから何がわかるかというと、

$$f^*(E^*,J)=f(E,J)$$ (11)

ということである。と、これは Young の論文にそった議論だが、単にリウビ ルの定理からもこの場合に $f$ が保存するのは当然な気もする。まあ、それ はともかく、結局、初めにコアの中心近くにいたとすると $f$ は一定なので、 結局ブラックホールの近くではやはり $f^*$ が一定となる。この時は、速度 が $r^{-1/2}$ で上がるので、 $f$ を一定に保つためには $\rho \propto r^{-3/2}$ でないといけないことがわかる。

つまり、ブラックホールが恒星系の中心で断熱的に成長する場合には、ブラッ クホールの十分近くでは $\rho \propto r^{-3/2}$ のカスプができることに なる。

これは美しい理論であり、また重要な結果でもあるが、直接に天文学的な応用 があるかと言われると難しい。ブラックホールがあるという観測的な傍証があ るのはおもに巨大楕円銀河と近傍の円盤銀河であるが、どちらも中心スロープ が $-3/2$ とは遠く離れているからである。

具体的には、巨大楕円銀河ではスロープが非常に浅く、 $-0.2$ から $-1$ 程 度の範囲に分布する。これに対して、我々の銀河系や近傍の円盤銀河では、 $-2$ 的と考えられている(我々の銀河系でもブラックホールに本当に近い、距 離にして 1 pc くらいのところなのでなかなか良くわからないが)。

巨大楕円銀河でカスプが非常に浅いことは、ブラックホール、巨大楕円銀河の 起源を考える上では大きな問題である。もしも、 QSO が巨大楕円銀河の直接 の projenitor であって、 QSO のガスが無くなって静かになったものが巨大 楕円銀河であるとするなら、ブラックホールの周りの恒星の分布は $-3/2$乗 カスプになりそうなものだからである。

それ以前に、そもそも $-1$ よりも浅いカスプを作る方法はあるのかというの も問題である。

ここまでは、

• 2体緩和による進化での定常状態

• 断熱的なブラックホールの成長に対する無衝突恒星系の応答

を考えた。もうひとつの可能性として、力学的な時間スケールでブラックホー ルが形成されるか、あるいは系の中心以外のところから落ちてくることが考え られる。

1.4 銀河の合体とブラックホールの合体

力学的な時間スケールで巨大ブラックホールを成長させるには、宇宙初期のゆ らぎから重力不安定で一気に作るという方法もありえる。しかし、現在の標準的な インフレーション宇宙モデルでは最初に重力不安定からコラプスする質量スケー ルはかなり小さく、その中のバリオン質量はもっと小さいので巨大ブラックホー ルを一気に作るのは無理がある。

現実的な階層的な構造形成モデルを考えると、大きな銀河はより小さな銀河が 合体することで形成されたということになる。この時に、銀河中心にあるブラッ クホールには何が起きるだろうか？、また、銀河の中心部の構造はどうなるの だろうか？

銀河同士の合体は、前に議論した「激しい緩和」の典型であり、十分に緩和が 進む前に構造が固まってしまう。ブラックホールが中心にある場合、これは、 合体してもブラックホールは初めから結構中心近くに行くということである。

まず、合体前の銀河がほぼ自己重力的でブラックホール質量よりも大きなコアを持っ ていた場合を考える。この場合には、コアは合体の時に大きく構造が変化する が、ブラックホールはコアの中ないしは非常に近くにいるであろう。ブラック ホールは回りの星よりも圧倒的に重いので、力学的摩擦によりコア中心に沈む。 2つブラックホールがあればそれらは連星になり、しばらくは回りの星をはね 飛ばして進化する。

合体前の銀河がコアをもたず、 $\rho \propto r^{-2}$ ないしはそれより急な ハローから、ブラックホールの重力が支配的な領域に滑らかにつながった構造 をしている場合も、合体中にブラックホールはほぼ中心まで沈み、2つあれば 連星ブラックホールになって近くの星をはね飛ばすことになる。

今、ブラックホール以外は無衝突系である極限的な場合を考えると、ブラック ホール連星は合体してできた銀河の中心にいって、重心運動の速度は 0 になっ て止まっている。ブラックホール連星の軌道長半径程度まで近づいた星は基本 的にははね飛ばされて無くなるので、ある程度よりもエネルギーが低い粒子は 存在しないことになる。細かいことをいうとはね飛ばされる条件はエネルギー ではなく近点でのブラックホールからの距離だが、とりあえずエネルギーに下 限がある場合を考える。

この時には、簡単な議論から中心に $\rho \propto r^{-1/2}$ のカスプがで きることがわかる。以下、 Nakano and Makino (1999, ApJL 525, 77) に従っ てこのことを示す。

分布関数がエネルギー$E$だけの関数だ とすると、ある半径 $r$ での密度は

$$\rho(r)=4\pi \int^{0}_{\phi(r)} f({E}) \sqrt{2\left[%% {E}-\phi(r) \right]} d{ E},$$ (12)

で与えられる。 これは速度空間での積分を、積分変数を $E$ にとって書 き換えただけである。ここで、 $f(E)$ はあるエネルギー $E_0$ よりも深い ところ($E<E_0$) では 0 であると仮定しよう。

ブラックホールに十分に近いところでは、ポテンシャルが十分に深いとすると、 積分区間の下限は $\phi(r)$ だが、そこに来る前に $f(E)$ は 0 になるので 上の式は

$$\rho(r)=4\pi \int^{0}_{E_0} f({E}) \sqrt{2\left[%% {E}-\phi(r) \right]} d{ E},$$ (13)

と書き換えられる。さらに、根号の中で、 $\vert\phi(r)\vert \gg \vert E\vert$ とすれば、 $E/\phi(r)$で展開できて、

$$\rho(r) = 4\pi \int^{0}_{\cal E_{\rm0}} f(E) \sqrt{2\left[%% E-\phi(r)\right]} d{ E}$$

$$= 4 \sqrt{-2\pi \phi(r)} \int^{0}_{ E_{\rm0}} f( E) \left[%% 1-\fra... ...}{\phi(r)}+O\left({\left[\frac{{ E}}{\phi(r)}\right]}^{2} \right) \right] d{ E}$$

$$\sim 4 \sqrt{-2\pi \phi(r)} \int^{0}_{ E_{\rm0}} f( E) d{ E}$$ (14)

最後の変形では結局展開の 0 次の項だけを残した。ここで、ポテンシャルが 十分に深い、ブラックホールが支配的な領域を考えていたので、

$$\phi(r) \sim - G M_{\rm BH}/r$$ (15)

であるから、結局

$$\rho(r) \propto \sqrt{-\phi(r)} \sim \sqrt{\frac{G M_{\rm BH}}{r}}.$$ (16)

となってブラックホールに十分に近いところでは $\rho \propto r^{-1/2}$ のカスプができることがわかる。

この結果は、本質的には

1. 分布関数が $f(E)$ で書ける(等方的である)
2. $E\rightarrow -\infty$ で f(E) が十分速く 0 になる

という 2 つの仮定だけによっているので、割合に一般性が高いと考えられる であろう。

とはいえ、上の仮定は結構強い仮定なので、実際に成り立つか？という問題は ある。まず、数値計算ではどんな風かというのを見てみよう。この辺は信用で きる数値計算結果はあまりないので、手前味噌になるが Makino and Ebisuzaki 1996 (ApJ 465, 527) を例にしよう。

図に、中心ブラックホールを持つ銀河(左)ともたない銀河(右)について、合体 したものを種にしてまた合体させるということを繰り返した時に密度プロファ イルがどう進化するかを示す。合体する毎に半径はほぼ2倍になるが、この図では横軸をスケールして重なるようにしている。

ブラックホールがある場合だと、合体を繰り返した時に中心に浅い密度カスプ が形成され、その領域の銀河の半径に比べた相対的な大きさはほぼ一定になっ ていることがわかる。この部分の傾きはほぼ $-1/2$ であり、 Nakano and Makino の理屈の通りになっているといえる。

まあ、論文の出版年から想像されるように、 Makino and Ebisuzaki は「結果 はこうなった」と書いているだけでなぜそうなるのかという説明は Nakano and Makino までちゃんとはなかったわけである。こういうは本当はよろしく ないという考え方もある。というのは、特に定性的な結果については、数値シ ミュレーションが正しいか間違っているかということの判断は、解析的な解釈 がつかない限り一般には非常に難しいからである。 NFW プロファイルの話の 時に議論したように、ユニバーサルな結果と思ったものが、ユニバーサルに計 算間違いしていたせいということもある。

ブラックホールがある銀河の合体の、この数値シミュレーションについては、 やった人は計算の信頼性についてはそれなりの自信があったわけだが、 Nakano and Makino の結果がでて初めて確信を持ってスロープが $-1/2$ にな るといえるようになった。

観測はどんな風かというと、想像されるように数値シミュレーションほど綺麗 になにかが決まるわけではない。が、まあ、傾向はある。

この図は Nakano and Makino からとったものであるが、 HST による Virgo Cluster の楕円銀河等の観測 (Gebhardt, K., et al., 1996, AJ, 112, 105, Faber, S. M., et al., 1997, AJ, 114, 1771) から、中心 (0.1'') の体積輝度密度の傾きと、回転速度と速 度分散の比をプロットしたものである。 シンボルの違いは絶対等級で、塗りつぶしてあるほうが明るい。

明るい楕円銀河はスロープが $-1$ よりも浅く、回転が小さいのに対して、暗 い楕円銀河がスロープが $-2$ 近くであり回転が卓越しているという傾向がか なりはっきりあることがわかる。

3 次元に密度を deprojection する時の誤差等の見 積もりに問題がある可能性がまだなくはないようだが、定性的に明るい銀河が 浅いスロープを持つという結果自体は確実なものといえるであろう。

1.5 質量分布のある恒星系の熱力学的な進化

最後に、「現在の球状星団の中心にブラックホールがあるとすればどんなふう に見えるはずか」という問題についてちょっと考えておこう。

2回前に話をしたように、現在 core collapse をしている球状星団では中心に ブラックホールが形成されるということはありそうにない。 が、球状星団が 生まれた直後の大質量星がある時には暴走的な合体から中質量ブラックホール が形成された可能性はある。

その後現在まで進化したら球状星団はどのように見えるか？というのがここで の問題である。単純に考えると、前回話をしたようにブラックホールの近くで の密度分布は $-7/4$ 乗のカスプになると考えられるが、これは等質量の星だ けを考えた場合である。ブラックホールの存在を考えない、 gravothermal catastrophe の場合には、シミュレーションでは中性子星や重い白色矮星は理 論通りのほぼ $-2.3$ 乗のカスプを作るが、見える星はもっと軽いためにもっ と浅いカスプになるというのは既に見た通りである。

ブラックホールがある場合にも同様に、中性子星や重い白色矮星はほぼ理論通りの カスプを作るが、見える星はずっと浅いというのが最近のシミュレーションの 結果(Baumgardt et al. 2004) である。

上の図は3次元的な密度分布を星の質量毎に書いたものである。初期条件はか なり深いキングモデルの中心にブラックホールをおいたものである。

ブラックホールの近くではどの質量でもカスプになっているが、そのスロー プは質量が小さいと小さくなることがわかる。明るい星 (turn-off mass の辺 り)の分布を2次元に投影したのが下の図で、基本的にフラットなコアを持つ 「普通の」キングモデルでフィットできる球状星団のように見える。

というわけで、この辺はまだ 100% 信用してもらっても困るが、球状星団の中心 にブラックホールがあるとすれば、その球状星団は比較的大きくフラットなコ アを持つように見えるということになる。