磁気リコネクション - 京都大学大学院理学研究科附属天文台
Download
Report
Transcript 磁気リコネクション - 京都大学大学院理学研究科附属天文台
コンピュータで探る
天体爆発現象
柴田一成
京大理学研究科・花山天文台
本講演のプラン
近年の観測が明らかにした驚くべき天体爆発現象を、最新の
スーパーコンピュータがどこまで明らかにしたか紹介
• はじめに
• 数値天体物理学(Numerical Astrophysics)入門
• 太陽フレア
– 原始星フレア
• 宇宙ジェット（活動銀河核、ブラックホール、原始星）
– （超新星とガンマ線バースト）
2001年４月６日 新入生フロンティアレクチャー
天体フレアと宇宙ジェット
柴田一成
京都大学理学研究科
花山天文台
はじめに
• １９世紀
永遠不変の静かな宇宙
↓
• ２０世紀前半 進化する宇宙
↓ 電波、Ｘ線、赤外線
• ２０世紀後半 活動する宇宙
宇宙は爆発だらけ！
星は進化する
• 形成
星間雲＝＞主系列星（太陽）
• 進化
主系列星（太陽）＝＞巨星
• 死 （軽い星：太陽） 惑星状星雲
＝＞白色矮星
（重い星） 超新星爆発
＝＞ブラックホール
中性子星
活動銀河核の発見
• １９６０年代ー 活動銀河が
続々と見つかる
クェーサー(QSO：準星）
電波銀河
セイファート銀河
＝＞中心核活動
きわめて遠方（宇宙初期）
宇宙最大の爆発
クェーサー ＝若い銀河の中心核
クェーサー
（３Ｃ２７３）
電波銀河（白鳥座A）
（距離＝５億光年、長さ＝３０万光年）
太陽フレア／ジェット
（京大飛騨天文台：Ｈα）
太陽プロミネンス噴出
（史上最大：１９４５年６月２８日）
私の青春の夢
宇宙最大の謎の活動銀河核（とそのジェット）を
生きているうちに解明したい！
しかし、本体の直接観測は遠い未来。
ならば身近な太陽面爆発（フレア）や他の類
似の天体フレア・宇宙ジェットにヒントを探し、
ＭＨＤプラズマ理論でせまる。（私の戦略）
天体
形成
（重力収縮）
天体
形成
降着円盤ー
ジェット
システムの
形成
恒星・降着円盤（銀河円盤）
内部におけるエネルギー発生による活動
日本のコンピュータは世界一！
地球シミュレータ(Earth Simulator)
http://www.hoise.com/primeur/02/articles/weekly/AE-PR-05-0259.html
Japanese 'Computenik'
Earth Simulator shatters
US supercomputer hegemony
Tokyo 20 April 2002 The Japanese Earth Simulator
is on-line and producing results that alarm the USA,
計算プロセッサのピーク性能8Gflops that considered itself as being leading in supercomputing
総プロセッサ数5120
technology. With over 35 Tflop/s, it five times outperforms
計算ノードのピーク性能64Gflops
the Asci White supercomputer that is leading the current
総計算ノード数640
TOP500 list. No doubt that position is for the Earth Simulator,
計算ノードの主記憶容量16GB
not only for the next list, but probably even for the coming
ピーク性能40Tflops
two years. In the New York Times, bench mark compiler Jack
主記憶容量10TB
Dongarra compares the event with the Sputnik, hence he
dubbed the Earth Simulator "Computenik".
数値天体物理学
Numerical Astrophysics
or Computational Astrophysics
• コンピュータを用いて、天体現象の物理学を
探求する学問
• コンピュータは、理論の望遠鏡
なぜ数値シミュレーションが
必要か？
• 基礎方程式（流体力学、電磁流体力学、一般
相対論、重力多体問題、、、）が非線形
• 天体爆発現象における物理量のダイナミック
レンジはきわめて広い＝＞非平衡開放系
流体方程式
（断熱、重力なし）
未知数５個： 密度（ρ）、速度ベクトル（ｖ）、圧力（ｐ）
方程式５個： 非線形連立偏微分方程式
質量保存
運動量保存
t
dv
( v) 0
p 0
dt
エネルギー保存
d
p
p v 0
dt ( 1)
d
v
ただし
dt
t
電磁流体（ＭＨＤ）方程式
（断熱、重力なし）
未知数８個： 密度（ρ）、速度ベクトル（ｖ）、圧力（ｐ）、磁場ベクトル（Ｂ）
方程式８個： 非線形連立偏微分方程式
質量保存
運動量保存
t
dv
エネルギー保存
誘導方程式
( v) 0
p
1
JB
c
dt
1 2
d
p
J
p v
dt ( 1)
ただし
B
c
c
rot ( v B
J)
J
t
4
rot B
さらに一般的な方程式は、
以下の式や項を含む
•
•
•
•
自己重力方程式（ポアソン方程式）
放射輸送方程式
非理想過程：放射冷却、熱伝導、粘性
ミクロ物理過程：宇宙線（非熱的粒子）、
原子・化学反応、、、
• 速度が光速に近づくと特殊相対論が必要
• ブラックホール近傍では一般相対論が必要
いずれにせよ、みな微分方程式
差分法
• 差分法： 微分を差分で近似。 有限の格子点上
で計算
t
n 1 n
t
• 粒子法： 超粒子の運動方程式を解く
ラグランジュ法。ＭＨＤには不向き。
（磁場を解く必要があるので、
格子点はどうしても必要）
参考文献
• 活動する宇宙（裳華房１９９９、柴田、福江、
松元、嶺重編）
• 数値流体力学ハンドブック（丸善２００２）
最終節「天体数値流体計算」（柴田）
• シミュレーション天文学最前線２００２
電子集録
http://th.nao.ac.jp/sim2002/proceedings/list.htm
詳しい文献リストについては、この中の私の原稿を参
照のこと
太陽フレア
Ｈα
太陽フレア
彩層
１万度
１９世紀中頃発見
黒点近傍で発生＝＞
磁気エネルギーが源
サイズ～（１－１０）万ｋｍ
全エネルギー
1029 - 1032erg
（水爆１０万ー１億個）
京大飛騨天文台
「ようこう」の見たコロナ
現在も観測中
軟Ｘ線望
遠鏡に
よる
（1 keV)
２００万度ー
数千万度
「ようこう」が見たフレア
軟Ｘ線(~1keV)
磁気リコネク
ション
磁気リコネクションとは？
磁力線
http://www-solar.mcs.st-andrews.ac.uk/~eric/
磁気エネルギーを短時間の内にプラズマの
運動エネルギーや熱エネルギーに変換
地球磁気圏でも発生
世界初の
熱伝導ＭＨＤリコネクション計算
（Yokoyama-Shibata 1998)
ようこう
軟Ｘ線望遠鏡
で観測された
太陽フレア
フレア温度の scaling law の
発見（Yokoyama and Shibata 1998, 2001）
8 p / B
T B
2
6/7
L
2/7
シミュレーションの効用
１） 定性的物理が良くわかる。教育的効果。
２） 天体物理学的モデリングを可能にする。観測と理
論の橋渡し。(e.g., Yokoyama and Shibata 1995)
３） 解析的手法では発見困難な物理の
発見の道具。シミュレーションは数値実験。
理論のカンニング（！）
例） scaling law (Yokoyama and Shibata 1998)
（世界初の熱伝導ＭＨＤリコネクションの
シミュレーションより発見）
太陽コロナにもジェット発見
（ようこう軟Ｘ線）
• 足元でマイク
ロ・フレア
• 長さ＝数万ー
数１０万ｋｍ
• 速度＝１０－
１０００ｋｍ／ｓ
太陽コロナ・ジェットのコンピュータ・
シミュレーション（横山・柴田１９９５）
ようこう軟Ｘ線望遠鏡
で観測された
太陽コロナＸ線ジェット
太陽コロナ・ジェットのコンピュータ・
シミュレーション（横山・柴田１９９５）
飛騨天文台
で観測された
太陽Hαジェット
コロナ質量放出（ＣＭＥ）
（ＳＯＨＯ／ＬＡＳＣＯ，
可視光／人工日食）
速度～１０－１０００ｋｍ／ｓ、質量～１０＾（１５）－１０＾（１６）ｇ
2D-MHD simulations of
emerging flux triggering mechanism for
flares/CMEs (Chen and Shibata 2000)
movie
原始星フレア
原始星フレア
• 超低温（１０Ｋ）の分子雲中の生まれたばかり
の星（原始星）で発生
• 時間スケールは太陽フレアとそれほど変わら
ない（数時間ー１０時間）
• 超高温、およそ１億度
（太陽フレアは、１－２千万度）
• 全エネルギーは、太陽フレアの１万倍以上
原始星フレア
（Ｘ線／あすか衛星：小山ら１９９５）
温度～
１億度
太陽フレアの
エネル
ギーの１
万倍以上
原始星フレア（ＹＬＷ１５Ａ）の
Ｘ線強度の時間変化
（あすか：Tsuboi et al. 2000）
原始星フレア
のモデル
（林、松元、
柴田 １９９６）
Long term evolution of dipole case
Uehara et al. 2003 in prep.
1 orbit ＝ 6.28
10.5AU
宇宙ジェット
宇宙ジェット
• 最初に発見された
宇宙ジェット
（１９１７）
• 楕円銀河Ｍ８７（電
波銀河）の中心核
から噴出
• 活動銀河核ジェット
の一種
活動銀河核ジェット
（電波銀河：Ｍ８７）
• 二つ目玉電波源
のエネルギー源
はジェット
• １光年以下のス
ケールから数１０
万年光年の大き
なスケールまで
同じ方向にジェッ
トが出ている
活動銀河核ジェットの根元に
降着円盤！
活動銀河核の正体は何か？
１．宇宙で最大級のエネルギー解放
＝＞大量の物質の静止エネルギー
の何割かを変換する必要あり
２．短時間の時間変動＝＞サイズは小
３．ジェットの速度は相対論的（光速に近い）
その正体は
超巨大ブラックホール（？）と
その周りを回る降着円盤
活動銀河核に
超巨大ブラックホールの証拠
野辺山宇宙電
波観測所の
電波天文学
者（中井、井
上、三好）が
発見
（１９９５）
ブラックホールと降着円盤
（コンピュータ・シミュレーション）
原始星ジェット
(ＨＨ１－２ ： 長さ約１光年）
（ハッブル宇宙望遠鏡：可視光観測）
電波銀河（白鳥座A）
（距離＝５億光年、長さ＝３０万光年）
原始星ジェットと降着円盤
宇宙ジェットの特徴（まとめ）
活動銀河核
中心天体
ジェットの
長さ
ジェットの
速度
脱出速度
近接連星系
原始星
超巨大ブラッ ブラックホー
クホール
ルまたは中
性子星
１００万光年 １０光年
原始星
c
0.3c - c
100km/s
c
0.3c - c
100km/s
１光年
宇宙ジェットは
いかにして発生したのか？
• 中心に星またはブラックホール
• そこにガスが降着＝＞降着円盤
• エネルギー源は重力エネルギー
• 重力エネルギーをいかにしてジェットの
運動エネルギーに変換するか？
宇宙ジェットの謎
• １光年以下から、１００万光年の大きなスケールま
でジェットが同じ方向に保たれているのはなぜか？
（＝＞回転）
• ジェットは、いかにして加速されるか？
• ジェットを、細長く絞っている力は何か？
• これらの謎を解くためのヒントが、
身近な太陽にある！！
太陽フレア／ジェットの
観測・理論からのヒント
• 直接のエネルギー源は磁気エネルギー（元
は核エネルギー）
• プラズマは磁気力によって加速される
• ヘリカル・フィラメント（磁力線）構造が普遍的
に出現
• 細長い構造は、磁気力によって容易に形成さ
れる
コロナはジェットだらけ！
（多くは、ヘリカル形状で回転
している）
宇宙ジェット（ハッブル）
SOHO/
CDS
SOHO/EIT LASCO
飛騨天文台／Ｈα
宇宙ジェットの電磁流体数値シミュレーション
降着円盤が
磁力線をねじ
ることにより
ヘリカル磁場
生成
磁力線
円盤
磁気力で加速
磁気張力で
ジェットを細く絞
る（ピンチ）
時間変化
ジェットの速度～ケプラー速度
Shibata & Uchida (1990、
1986)
ヘリカル原始星ジェット
宇宙ジェットの電磁流体モデル
（工藤、松元、柴田 ２００４）
Middle
Evidence Small
of collimation (Kudoh et al. 2004,
in prep)
Large
流体力学の新しい数値解法：CIP 法
（矢部教授（現東京工大）が１９９１年に開発）
接触不連続が精度良く解ける。
気体、液体、固体が同時に解ける
Numerical simulation of accretion disk
(Kudoh, Matsumoto, Shibata 2002, PASJ)
Magnetorotational Instability (Balbus and Hawley 1991)
leads to turbulence and reconnection
Similar to
ADAF
General
relativistic
MHD
simulation
(Koide,
Shibata,
Kudoh 1998,
1999,2000)
Maximum Lorentz factor
～２＝＞real limit ?
MHD jets from Kerr hole
magnetosphere
Koide, Meier, Kudoh,
Shibata (2000)
Koide et al. 2002 Science
むすび
• 近年の太陽観測とコンピュータシミュレーションの発展により、
太陽フレアのエネルギー解放は磁気リコネクションによると
いう説が、ほぼ確立した。原始星フレアも、磁気リコネクショ
ン説が有力。ただし、リコネクションの物理は未解決。
• 宇宙ジェット現象は、太陽面爆発現象と類似の電磁流体力
学的機構によって発生しているかもしれないという説が、活
発に議論されるようになってきた。ただし、まだ解決にはほど
遠い。
• 観測が発展するにつれ、宇宙はますます謎めいた爆発現象
に満ちていることが明らかになってきた。 例えば、ガンマ線
バースト。
ガンマ線バースト
• 発見以来２０年
以上、どこで起き
ているかも謎だっ
たが、最近ようや
く数１０億光年か
なたの銀河で発
生していることが
判明
• 正体はいまだ謎
（もしや電磁
ジェット？）
Central engine of gamma ray
bursts ?
Binary neutron
star merger model
Collapsar model
magnetar model
MHD simulation of collapsar
Mizuno et al. (2003) submitted to ApJ
general relativistic MHD simulation with
Schwarzchild black hole
Vjet ~ 0.2c
はじめに／おわりに
• １９世紀
永遠不変の静かな宇宙
↓
• ２０世紀前半 進化する宇宙
↓ 電波、Ｘ線、赤外線
• ２０世紀後半 活動する宇宙
宇宙は爆発だらけ！
↓ γ線、ニュートリノ、重力波、
スーパーコンピュータ
21世紀前半
？？する宇宙