プレスリリース:天の川銀河中心の超巨大ブラックホールに2013年夏に最接近し輝きを増すガス雲について

プレスリリース

天の川銀河中心の超巨大ブラックホールに2013年夏に最接近し輝きを増すガス雲について

2013年3月20日
日本天文学会

概要

fig0

ガス雲「G2」が超巨大ブラックホールSgr A*(白球)の強力な重力の影響で引き延ばされ、薄く潰されながら近傍を通過する。このときガス雲の輝きが急激に増す。この図は交差法による立体可視化可能にしたもの。クレジット:斎藤貴之(シミュレーション)、武田隆顕(可視化)

東京工業大学・斎藤貴之特任准教授らの研究グループは、3次元高分解能コンピュータシミュレーションを用いて、天の川銀河の中心にある超巨大ブラックホールSgr A*(サジタリウス エースター)に接近するガス雲「G2」が、今年夏~秋を中心に一年間明るく輝く可能性があることを明らかにしました。
シミュレーションによると、このガス雲はSgr A*の強力な重力の影響で、ガスの進行方向に延ばされるとともに、ガスの軌道平面に対して垂直方向に強力な圧縮を受けます。この圧縮によって加熱されたガス雲が明るく輝くと推測されます。このような圧縮によるエネルギー解放メカニズムは、本研究が世界で初めて示しました。
現在、世界中の望遠鏡が、このガス雲「G2」に何が起きるか、そしてSgr A*にどのような変化が訪れるかを明らかにするために観測を続けています。シミュレーションと観測データをつきあわせることで今後Sgr A*周辺環境の理解が大幅に進むものと期待されます。

研究グループ
斎藤貴之(東京工業大学・特任准教授/ELSI Affiliated scientist)、牧野淳一郎(東京工業大学・教授/ELSI Principal Investigator/理化学研究所計算科学研究機構チームリーダー)、朝木義晴(ISAS/JAXA・助教)、馬場淳一(東京工業大学・特任助教)、小麦真也(国立天文台・助教/合同アルマ観測所・CSVサイエンティスト)、三好真(国立天文台・助教)、長尾透(京都大学・特定准教授)、高橋真聡(愛知教育大学・教授)、武田隆顕(国立天文台・特任助教)、坪井昌人(ISAS/JAXA・教授)、若松謙一(岐阜大学・名誉教授)

※本記者発表は、日本天文学会2013年春季年会(3月20日~23日、埼玉大学)開催に先んじて行われました。

1.背景

現在ほぼすべての銀河はその中心に銀河質量の0.1%程度の超巨大ブラックホールを持つことが知られています。超巨大ブラックホールは一般に活発な活動を示し、銀河の形成と進化に重要な役割を果たしていると考えられています。我々の住む天の川銀河も例外ではなく、中心部に太陽のおよそ430万倍の質量を持つSgr A*(サジタリウス エースター)と呼ばれる超巨大ブラックホールが存在します。Sgr A*は過去の活動性の名残が報告されているものの、現在は不活性状態にあります。そのため、Sgr A*周辺環境はどうなっているかについて多くの研究者が注目しています。
2012年1月、ドイツマックスプランク研究所のGillessenらを中心とした研究グループが、Sgr A*に接近する地球の三倍の質量をもつガス雲「G2」の存在を報告しました(図1)。詳細な観測データを積み重ねることで、ガス雲「G2」がSgr A*の周りを運動する軌道にあり、2013年夏(その後の更新で秋と報告された)にSgr A*に最接近するということが明らかになりました(図2)。現在世界中の天文学者と望遠鏡がこの現象に注目しています。

fig1

図1 Gillessenらによって報告されたガス雲「G2」の位置。三つの色はそれぞれ異なる時刻の位置を表している。白い十が天の川銀河の超巨大ブラックホールSgr A*の位置を表している。S2はSgr A*の周りを運動している恒星。縦の白い棒はスケールを表しており、40ミリパーセク(mpc)、およそ8000天文単位(天文単位:地球と太陽の平均距離)。Gillessen et al. 2012 Nature 図1(c)より一部改変。

fig2

図2 ガス雲「G2」の観測された位置(赤丸)とそれから予想される軌道(赤実線、赤点線)。赤の点線は予測される今後の軌道を表している。青い点はSgr A*の近くを運動している恒星(S2)の観測された位置(青点線)と軌道(青実線)。オレンジの点のところにSgr A*が存在する。横軸、縦軸はそれぞれ赤経、赤緯を表している。Gillessen et al. 2012 Nature 図1(d)より、一部改変。

2.研究の詳細

fig3

図3 我々の計算セットアップ。ガス雲(ピンク)を観測によって得られたSgr A*(オレンジ)を回る楕円軌道上に置き、1995年からの状態の変化を追った。青実線が今回のシミュレーションで追った全期間に対応する。2013年夏-秋にSgr A*に最接近する。クレジット:斎藤他(2013)

もしガス雲「G2」が十分な広がりを持っていれば、太陽の430万倍の質量を持つ超巨大ブラックホールSgr A*の強力な重力下で運動の方向に引き延ばされる他に、運動面に垂直な方向に強力な圧縮を受けるはずです。そうすると、このときの圧縮のエネルギーを解放することでガス雲が明るく輝くことが予想されます。我々は、世界で初めて3次元かつ現実的なガスの振る舞いを考慮したコンピュータシミュレーションを行い、このガス雲「G2」が 2013年にSgr A*の極近傍を通過するときに、その強い重力の影響で何が起こるかを予測しました(従来の研究との比較は補足資料3を見てください)。シミュレーションでは観測から得られたSgr A*を回る軌道上にガス雲を置いてその運動を解きます(図3を参照)。現時点では、ガス雲の形状についてはっきりとわかっていないため、ここでは地球の3倍の質量を持つガス雲が一様の密度を持っており、その半径が125天文単位であると仮定しました。観測データに従い、ガス雲はSgr A*の周りを楕円軌道で運動するとします。また、これまでの観測事実からの示唆に従いSgr A*の周りには高温で薄いガスが広がっているとしました。
このシミュレーションは、HCPI戦略プログラム分野5「物質と宇宙の起源と構造」の研究の一環として東京工業大学の斎藤貴之(東京工業大学特任准教授/ELSI Affiliated scientist)が開発した宇宙流体シミュレーションコードASURAを用いて行いました。非常に大きな天体シミュレーションを高速に、そして正確に実行するために、独自の計算手法が多数組み込まれています(用語集参照)。シミュレーションには国立天文台天文シミュレーションプロジェクトのスーパーコンピュータCray XT4を用いました。

3.結果

fig4

図4 ガス雲の明るさの変化。横軸は西暦を表し縦軸はガス雲から一秒間に放出されるエネルギー(=明るさ)を表す。黒い線がシミュレーションから得られた明るさの変化。左側のパネル(a)は1995年から2030年まで、右側のパネル(b)は2011年夏から2015年夏までの期間のガス雲から放出されるエネルギーの量を表している。なお、太陽光度は4x1033erg/s。クレジット:斎藤他(2013)

図4はコンピュータシミュレーションで得られたガス雲の明るさの変化を表しています。この図より、ガス雲がSgr A*の極周辺を通過するときに非常に明るく輝くことがわかります。また、この輝く期間はおよそ一年間になります。最も明るい時期の明るさは、太陽の絶対等級(本来の明るさ)の50倍程度です。Sgr A*の周りにある高温ガスの量が異なる場合でも増光は高温ガスの影響を受けません(参考資料1を参照のこと)。このような増光の可能性は、1)三次元で、2)非常に高い分解能のシミュレーションを達成することによって、我々が今回初めて指摘することができました。

国立天文台4D2Uプロジェクトの武田隆顕(国立天文台特任助教)が開発した可視化ソフト Zindaiji3(詳細は用語集参照)を用いて三次元的な様子を表したのが図5およびアニメーションです。ガスから出てくる光によって色をつけています。暗いところは緑色に、たくさん光が出てくる明るいところはオレンジ、そして白っぽくなっています。ガス雲がSgr A*に近づくにつれて、その強力な重力の影響-潮汐力-により引き延ばされていることがわかります。このときガス雲は非常に薄くなります。特にSgr A*の近傍通過時にはシミュレーション開始時の1/100以下になっています。この効果もSgr A*の強力な重力による潮汐力の影響です。ガス雲がつぶされるにつれて、ガスの温度が上がりそのエネルギーを光で解放しているため、ガス雲が輝くようになります。図6では交差法をもちいることで、最接近のときのガス雲の姿を立体視することができます。補足資料2にはある時刻のガス雲の姿を異なる2方向からみた図があります。

fig5

図5 アニメーションより抜粋。左上の数字が西暦を表している。ガス雲がSgr A*(中心の白丸)に接近する過程で引き延ばされ、そして薄くなっていることがわかる。白丸のサイズは300 Rs、ここでRsはSgr A*の事象の地平面(シュバルツシルト半径)、と便宜上大きく取ってあり、本来のSgr A*の事象の地平面よりずっと大きいことに注意。アニメーションはウエブページよりダウンロードできる。クレジット:斎藤貴之(シミュレーション)、武田隆顕(可視化)

fig6

図6 ガス雲がSgr A*を通過する時(2013年春)の姿。交差法による立体視が可能。交差法とは、右目と左目でそれぞれ左側のパネルと右側のパネルを見るようにすることで立体視する方法のこと。ガス雲が非常に薄くなって通過する様子がわかる。クレジット:斎藤貴之(シミュレーション)、武田隆顕(可視化)

4.まとめと今後

本研究では、天の川銀河の超巨大ブラックホールSgr A*に接近しつつある三地球質量ガス雲「G2」の運動をコンピュータシミュレーションし、とくに最接近する2013年にガス雲「G2」に何が起きるかを調べました。我々の3次元シミュレーションによると、ガス雲「G2」はSgr A*の強力な重力の影響で軌道面に伸ばされるとともに軌道面鉛直方向に圧縮され、その圧縮のエネルギーを光として放出するために明るく輝くことがわかりました。もちろんこれは現在までに知られているガス雲「G2」についての情報をもとに行われたシミュレーションであり、あくまで可能性です。現在、Sgr A*にガス雲「G2」が接近することでどのような変化が起きるかをとらえるために、すばる望遠鏡、ALMA望遠鏡、チャンドラ望遠鏡など地上と衛星軌道にある世界中の望遠鏡が Sgr A*周辺に注目しています。今後の観測によりこのガス雲に実際に何が起きたかが明らかになり、本研究の確からしさも明らかになるものと思われます。
ガス雲「G2」についての研究は、天文学の研究としては極めてまれな「答え合わせのできる研究」になっています。その意味でも非常に興味深い天体です。このガス雲の研究により、超巨大ブラックホールSgr A*周辺環境や、その活動性のメカニズムについての理解が進むものと期待されます。

補足資料、用語集など

fig7

図7 Sgr A*の周りにある高温ガスの量を変えた場合の明るさの変化を表示。黒い線は図7と同じもの。赤い線はSgr A*の周りにある高温ガスの質量が1/10の場合の明るさの変化、青い線は高温ガスがない場合の明るさの変化。2013年の増光はどれも同程度。クレジット:斎藤他 (2013)

補足資料1:Sgr A* 周りの高温ガスの影響
Sgr A*の周りにどの程度の高温ガスがどのように分布しているかは、実はまだはっきりしていません。そこで我々はその影響についても調べました。図7は、Sgr A*の周りを取り巻く高温ガスの量を変えたときのガス雲の光度の変化の様子を表しています。黒い線(Run 1)が図4に書かれていた線に対応します(高さが異なっていることに注意)。Run 2は高温ガスの質量を1/10にしたときの振る舞い、また、Run 3は高温ガスを無くしたときの振る舞いになっています。ガス雲が Sgr A*に最接近したときの明るさは三つのモデルで大差がありません。したがって、圧縮による増光には、高温ガスはほとんど影響しないということがわかりました。我々の見積もりによると、むしろガス雲がどの程度の質量を持っているか、またどのぐらいの広がりを持っているかが鍵となります。

fig8

図8 ガス雲の形状の変化の様子。特にSgr A*に最接近する時期(2006年春-2013年夏)を表示。パネル(a,b,c,d)は軌道面上でのガス雲の形状。パネル(a’,b’,c’,d)は、ガス雲を軌道面垂直方向にカットし、高さ方向の形状の変化を表示した。色はガスのおのおのの場所における密度に対応している。それぞれのパネルの左上に書かれている数字が時期(西暦)を表している。また、それぞれのパネルの右上に書かれている白い横棒はスケールを表している。単位は天文単位(au)。パネル(e)には、ガス雲の軌道(黒い楕円)、ガス雲の初期位置(青六角形)、またパネル(a,b,c,d)に対応する領域(赤四角)が書かれている。アスタリスク(*)がSgr A*の位置を表している。クレジット:斎藤他(2013)

補足資料2:ガス雲の変形、二つの方向から
図8は、とくに最接近する途中のガス雲の状態の変化の様子を表したものです。軌道面と軌道面垂直方向からみることで、実際のガス雲の立体構造が把握できます。
パネル (a,b,c,d) は、軌道面でのガス雲の形状の時間変化を表しています。また、パネル(e) がそれぞれのパネルの軌道面上での位置を表しています。Sgr A*の強力な重力でよりガス雲が引き延ばされていることがわかります。パネル (a,b,c,d) にある白い線の方向にガス雲をカットし、軌道面垂直方向の形状の変化を見たのがパネル (a’,b’,c’,d’) です。パネル(a’,b’) と パネル (c’,d’) では倍率が4倍違うことに注意してください。
これらのパネルを見るとガス雲は徐々に高さ方向に潰されて薄くなっていくことがわかります。これもSgr A*からうける潮汐力の効果です。特にSgr A*の近傍を通過する際にはシミュレーション開始時の1/100以下程度まで薄くなります(パネル c’, d’)。

補足資料3:他のコンピュータシミュレーションによる研究との比較
ガス雲「G2」について、これまでに行われてきたコンピュータシミュレーションによる計算結果と我々の結果の比較をします。最初にガス雲「G2」の状態の変化についてのシミュレーションを行ったのは、Schartmann et al. ApJ Vol. 755, Article id:155 (2012)です。彼らは、軌道面内だけに限定した2次元のガス雲のシミュレーションを多数行い、ガス雲の変形の様子を調べました。したがって、このシミュレーションでは、ガス雲が軌道面垂直方向に潰される過程がそもそも含まれていません。また、簡単のために放射冷却過程(ガスが持つエネルギーを電磁波で放出して失う過程)について扱っていないため、そもそもどの程度輝くかの予測ができないという問題があります。
3次元のシミュレーションを最初に行ったのは、Annios et al. ApJ Vol. 759, article id 132 (2012)でした。しかし、3次元になった分、シミュレーションの空間分解能が低く、最小でも4天文単位までしか表現できていません。また、彼らも放射冷却過程を解いていません。彼らのシミュレーションを見ると、ガスが十分につぶれることができていません。空間分解能が足りないこと、ガスのエネルギーが失われないこと、この2つが今回我々が指摘した増光現象を示唆できなかった原因だと考えられます。
我々のシミュレーションでは、(1)3次元で非常に高い空間分解能(0.1天文単位まで分解)をもち、(2)ガスの放射冷却過程を考慮したことにより、超巨大ブラックホールの強力な重力の影響でガス雲「G2」が圧縮され、それにより明るく輝くという一連の変化を指摘することに成功しました。そしてこれは、従来のコンピュータシミュレーションでは表現できていなかった現象です。このような高い分解能を持ち放射冷却過程を考慮したシミュレーションは、計算量が増えるために非常に大きな困難を伴います。このシミュレーションの達成には、斎藤が開発した並列シミュレーションコードASURAの高い性能が遺憾なく発揮されました。
我々は、コンピュータシミュレーションの他に数式を用いてこの増光がもっともらしい明るさになっているということを確認してあります(詳細割愛します)。

単位について
太陽質量:およそ2×1030kg
地球質量:およそ6×1024kg
天文単位(au):およそ1.5×108km
パーセク(pc):およそ3.1×1013km
※詳しくは理科年表(丸善出版)などを参考にしてください。

用語集
超巨大ブラックホール:恒星の質量を遙かに上回る巨大(大質量)なブラックホールのこと。典型的な質量は数100万太陽質量から1億太陽質量程度。

Sgr A*:天の川銀河の中心部分に存在する、超巨大ブラックホール。天の川銀河中心部が天球面上で射手座の方角にあることからこの名をつけられた。A*なのは、まず銀河中心方向の明るい電波源を順に番号をつけたときに一番明るかったということでAがついたこと、そしてさらに高い分解能で見たときに点源のように見えたことから*をつけたことによる。現在、Sgr A*の周りを運動する星の詳細な軌道が調べられており、それよりSgr A*の質量がおよそ430万太陽質量であると推定されている。

fig9

図9 潮汐力による引き延ばしと圧縮の効果。青丸で表現された天体1には、天体2(Stellite)によって引き起こされた潮汐力により、天体2がある方向に引き延ばされ、また、天体1-2を結ぶ線に対して垂直方向に圧縮を受ける。出典:wikipedia Tidal fieldより。

潮汐力:潮汐力は、広がりを持った天体1に別な天体2からの重力が働いているときに生じる2次的な力である。天体1では、天体2がある方向に引き延ばされ、それ以外の方向には圧縮される力が働く(図8を参照のこと)。地球上では、月および太陽の重力の影響で潮の満ち引きを引き起こす。

シミュレーションコード ASURA:東京工業大学の斎藤貴之特任准教授/ELSI Affiliated scientistが開発した、銀河形成シミュレーションのための並列計算コード。N体/SPH法と呼ばれる計算法を採用している。N体法は重力計算を質点の間に働くニュートン重力として計算する方法である。この重力計算は一般に非常に負荷が高い。ASURAでは、GRAPEと呼ばれる重力計算専用機を使って重力計算できるように設計されている。本研究ではこのGRAPEのソフトウェアエミュレータであるPhantom-GRAPEを用いている。SPH法とは、流れを構成する要素(流体粒子)の重ね合わせで密度などの物理量を表現する流体計算法である。本研究で用いたSPH法は、斎藤と牧野が開発した、従来のSPH法で指摘されていた問題点を克服した新しいタイプのSPH法である。

可視化ソフト Zindaiji:国立天文台4D2Uプロジェクトの武田隆顕特任助教が開発した、粒子系シミュレーション専用可視化ソフト。市販のソフトでは難しい大量の粒子情報を元に可視化を行うために独自開発された。立体視、ドーム投影用データを作成できる。これまでZindaijiを用いて数多くの天文シミュレーションの可視化が行われており、その成果は国立天文台4D2Uプロジェクトのドームを代表に日本国内のプラネタリウムや、世界中のプラネタリウムプログラムにも利用されている。開発ページはこちら。また、これまでの成果を見ることができる。なお、Zindaijiの命名は国立天文台のある三鷹市のそばにある深大寺にちなむ。

HPCI戦略プログラム分野5「物質と宇宙の起源と構造」:スーパーコンピュータ「京」に代表される国内のスーパーコンピュータを高度利用して研究を進めるハイパフォーマンスコンピューティングインフラ(HPCI)の戦略プログラムの一つ。ビッグバンに始まる宇宙の歴史の中での素粒子から原子核、星・銀河形成に至る物質路宇宙の起源と構造を統一的に理解することを目指している。ウエブページ

関連情報
ウエブ上の資料(デジタル画像ファイル、アニメーションファイル) http://v1.jmlab.jp/~saitoh/G2/Press/
アーカイブサーバに登録されている論文原稿 http://arxiv.org/abs/1212.0349
関連ウエブページ http://www.phyas.aichi-edu.ac.jp/~takahasi/Project_Horizon/index.html

日本天文学会2013年春季年会における関連講演
企画セッション「超巨大ブラックホールの起源」:B07a “三地球質量ガス雲 G2 の近点通過によるフレア現象について” 斎藤貴之(東京工業大学)本研究
企画セッション「超巨大ブラックホールの起源」:B06a “銀河系中心SgrA*ブラックホールの周辺環境” 高橋真聡(愛知教育大学)
企画セッション「超巨大ブラックホールの起源」:B08a “短基線VLBIによるSagittarius A*の2013年事象の観測” 坪井昌人(宇宙航空研究開発機構)

問い合わせ先

斎藤貴之 東京工業大学・特任准教授/ELSI Affiliated scientist
TEL 03-5734-2243 FAX 03-5734-3538
E-mail: saitoh(-at-)geo.titech.ac.jp
URL http://v1.jmlab.jp/~saitoh/G2/Press/

カテゴリー: プレスリリース, 新着情報   パーマリンク