2005年12月27日 基礎物理学研究所研究会「高エネルギー天体物理学の最前線」

宇宙ジェットの形成：MHDモデルか否か？
高原 vs 柴田 論争(2001) revisited

柴田一成 京大理花山天文台（工藤哲洋 国立天文台）

講演内容• 宇宙ジェットの形成： ＭＨＤモデルか否か？
高原ｖｓ柴田（２００１年3月@宇宙研） revisited– ＭＨＤジェット駆動には，大局的磁場の存在が必要か？
– エネルギー解放率は十分か？
– ローレンツ因子の大きなジェットを作るにはどうすれば良いか？
– 岡本論争について: ＭＨＤジェットのコリメーション– Ｅ（粒子）＞＞Ｅ（磁場）という観測をどう考えるか？
– 今後の課題: 加速と構造, collimation, 観測• その後の進展（２００１－２００５）• 残された問題点

2001年３月２９日 「ディスクとジェット」研究会（宇宙研）

宇宙ジェットの形成：
ＭＨＤモデルか否か？
高原ｖｓ柴田柴田一成、工藤哲洋

２．宇宙ジェットの理論モデル
原始星の場合宇宙ジェットの特徴（まとめ）活動銀河核中心天体

ジェットの長さ

近接連星系

超巨大ブラッ ブラックホー
クホール
ルまたは中性子星
１００万光年 １０光年

星形成領域
原始星

１光年

ジェットの速度

光速

0.3 – １光速 100km/s

脱出速度

光速

0.3 – １光速 100km/s

宇宙ジェットの磁気流体モデル
磁場と回転によってアウトフローを加速するモデル

ガス圧加速モデル 音速<<ジェットの速度
磁力線

放射圧加速モデル 放射の運動量<<ジェットの運動量
磁気力で加速するモデル（磁気流体モデル）

ジェット
降着円盤
原始星、ブラックホール

磁気遠心力
磁気圧

Blandford, Lovelace, Uchida-Shibata, Pudritz-Norman,Shu,,,1

宇宙ジェットの電磁流体数値シミュレーション

宇宙ジェットの電磁流体モデル（工藤、松元、柴田 ２０００）

Shibata & Uchida (1986)
降着円盤が磁力線をねじることによりヘリカル磁場生成

磁力線

円盤

磁気力で加速磁気張力でジェットを細く絞る（ピンチ）時間変化
ジェットの速度～ケプラー速度（円盤の回転速度）～中心天体近傍では、その脱出速度

ＭＨＤジェット駆動には，大局的磁場の存在が必要か？

高原： 活動銀河核や連星系では大局的磁場は弱い（？）
からＭＨＤジェットはできない
工藤・柴田：
１）大局的磁場は弱くても、ＭＨＤジェットはできる！
２）局所的磁場でも、ＭＨＤジェットはできる！

磁場に対する依存性 (from 2.5D-MHD
simulation by Kudoh et al. 1998)
質量放出率

定常理論（Kudoh & Shibata 1997）

ジェットの速度／Ｖｋ

1/ 6
∝ Emg

max
1/ 2
∝ Emg

2

V Emg =  A  ＝磁気エネルギー／重力エネルギー
 VK 2

V Emg =  A  ∝ B 2p
 VK 2磁気遠心力 ｖｓ 磁気圧

Jet velocityV jet ≈ Bϕ /( 4πρ )1/ 21/ 3 Ω2 B p2 r 4 ≈ & M

Michel (1969) velocity

1/ 3

V 
1/ 3
≈  A  Vk ∝ B p
 Cs 

原始星フレアのＭＨＤmodel

強い弱い

直線的ぎりぎり巻き

ρCs r

ρ Cs r 2

2

1/ 3

Ｖｋ：ケプラー速度（円盤回転速度）

VA = B p /(4πρ )1/ 2

局所的磁場でも、ジェットは加速できる
1. 降着円盤にダイポール磁場
2. 中心天体にダイポール磁場• ポロイダル磁場• 降着円盤近傍の磁力線形状• 質量流出量• 最終速度 (V∞ )• 適用範囲 V2 Vk  A 
 CsVk Bp
B 1/ 3

V Vk  A  Cs Emg ,c < Emg < 1Emg < Emg ,cE mg ≡ 磁気エネルギー／重力エネルギー

E mg ,c ≡ ρ s / ρ d

ダイポール磁場の場合(Kudoh et al. 2001)

２成分
Polar jet
Equatorial wind
（ただし、非定常過程？）

ローレンツ因子～１０は磁場加速説では困難？

（Hayashi, Shibata,Matsumoto 1996)3

ローレンツ因子の大きなジェットを作るにはどうすれば良いか？

一般相
対論的ＭＨＤジェット• パルサー風と良く似た状況設定を作る(Koide,Shibata,Kudoh
1998、1999,2000,Aoki et al.)

– 急激に広がる FLOW TUBE（磁場）– 非常に希薄なコロナ（ジェット）＝＞ダイポール磁場
電子・陽電子プラズマ

最大のローレンツ因子～２＝＞real limit ?

ローレンツ因子

Flow tube の断面積

A ∝ r − ( 2+α )Flow tube の断面積

A ∝ r − ( 2+α )ＭＨＤジェットのコリメーション
岡本論争について• 岡本（天文月報 ２０００年３月）「宇宙ジェットーあるパラダイムの終焉」
宇宙ジェットは磁場によってコリメートできるか？
「定常・軸対称ＭＨＤアウトフローはセルフコリメイトしない」

4

岡本説の論拠

これまでのコリメーションに関する研究
2次元軸対称定常解

桜井（1985）

天体の近く

遠方

θ
Z

アウトフローは２成分（polar jet
Equatorial wind)rrR

ダイポール磁場の場合(Kudoh et al. 2001)

２成分
Polar jet
Equatorial wind
（ただし、非定常過程？）

ＪＸＢ力

岡本説への反論（天文月報別冊）• 桜井：ダイポール磁場も、十分時間がたつとスプリットモノポール磁場になる（ただし、その証明はまだない）• その他の反論（新田、内田、工藤、他）

SmallMiddle
コリメーションの証拠(Kudohet al. 2001)

Ｅ（粒子）＞＞Ｅ（磁場）という観測をどう考えるか？

Large5

ジェット中の磁場は弱いので磁場仮説は困難？

高原：観測によればＡＧＮジェット
では、粒子エネルギーが磁気エネルギーより桁違いに大きい。
工藤、柴田：それでもジェットは磁気的に加速できる

観測事実：AGN jet 中の粒子エネ
ルギー＞＞磁気エネルギー• 解決策• １） 磁気チューブの断面積が急激に広がっている (Kudoh et al. 1998)• ２） 磁気リコネクションにより、磁気エネル
ギー＝＞粒子エネルギー(Blandford 2000)

MHDアウトフローの性質
Michel(1969)
放射状の磁場を仮定すると、磁場のエネルギーが、十分運動エネルギーに変換されずに無限遠方に達する。
（fast point が無限遠方。）
Begelman & Li (1994)
放射状よりも開いた磁束管の中を流れるアウトフローは、すべての磁場のエネルギーを運動エネルギーに変換できる。

放射状より開いた磁力線
放射状

アウトフローのコリメーションとも関係している。

磁力管の開き具合は、アウトフローのコリメーションをきちんと求めないとわからない。しかし、それを解くのは難しい。

Takahashi & Shibata (1998)の研究をAGNジェットに応用

そこで、磁力管の開き具合を適当な関数で与えて、磁力管に沿った流れのみを求めるという研究方針をとる。

降着円盤からのアウトフローr ≈ 1012 rNS ≈ 109 rLC

Takahashi & Shibata (1998)
流線に沿ってr 2 B p ∝ r −α (α = 0が放射状）

カニ星雲

回転アウトフローパルサー風と仮定し、赤道面の流れを解きパルサー風の問題に応用した。

α ≈ 0.4 であれば観測を説明

ブラックホール

衝撃波
運動エネルギー＞磁場のエネルギー
1： 0.0033

降着円盤•一般相対論的MHD方程式（定常、軸対称、Cold）•流れは一次元。赤道面。
•ブラックホールの最内安定軌道からアウトフローが流れている。
磁力線の角速度はそこでのケプラー角速度。
•最終速度はγ＝１０。

r 2 B p ∝ r −α を仮定。
ブラックホール半径の100倍から1000倍の所で、磁場のエネルギーが運動エネルギーの10分の1以下になるための条件を調べる。

6

結果１

結果２ブレーザーの磁場ブレーザーの条件を満たすためには、

議論

α > 0.6Bogovalov & Tsinganos (1999)

α > 0.6 は現実的にありうる値か？

非相対論

単純なMHDアウトフローの理論では現実的ではないと考えられている。

相対論 γ = 300

磁力管の開き具合表すパラメータαは、アウトフローのコリメーションと関係している。
ログスケールでコリメーション放射状より開いた磁力線
放射状Sakurai (1985)

非相対論

これまでのコリメーションの研究では、磁力線はもっとゆっくり、ログスケールでコリメーションするという結果が得られている。

相対論的MHDアウトフローはさらにコリメーションしにくい
磁場のエネルギーが運動エネルギーに変換されにくい。
Li (1996), Bogovalov (2001)

今後の課題: 加速と構造

今後の課題: コリメーション• ジェットの内部構造（ノット）の起源• ジェットの３次元安定性（＝＞磁気リコネク
ション＝＞粒子加速）• 相対論的ジェットの場合：
ローレンツ因子＞１０の高速ジェットはいかにして形成されるか？
＝＞ガンマ線バーストへの応用• 電子陽電子磁気的ジェットの物理？

• 回転星ダイポール磁場では、定常状態でどれだけの磁束が極方向にコリメートするか？
• 定常軸対称ＭＨＤアウトフローの数値解（今までのところ、全空間をセルフコンシステントに解いたのは、桜井(1985)のみ）＝＞数理物理の超難問

7

今後の課題： 観測• ＭＨＤジェットの観測的証拠をさがす– ジェットの回転とヘリカル磁場– ジェットの根本の降着円盤の磁場測定• ジェットのコリメーションのスケールを観測的に解明• ジェットとフレアの関係の検証

その後の進展
２００１－２００５

MHD model of GRBs
(Spruit et al. 2001 A&A 369, 694 )リコネクションを利用したエネルギー変換モデル

Flare/CME model ofGRBs driven by magneticreconnection(Aoki, Yashiro, Shibata 200４)
Only large flares canproduce large plasmoids
that have enough energy toescape frommagnetosphere of thecentral engine of a GRB(cf. Negoro and Mineshige2002)

一般相対論的MHDシミュレーション

8

MHD jets from Kerr holemagnetosphere

Vmax=0.86c(Lorentz factor 2.0)

Koide （２００４）

τS=rS/c(Unit of time) Bφ2 log 2  ρc 
z/ rS
Kerr
blackhole

Koide, Meier, Kudoh,
Shibata (2000)

Koide et al. 2002 SciencecR/rS

ErgosphereLines: Magnetic fieldsurfaces
Arrows: Velocity ofplasmaRelativistic Outflow driven by
Magnetic Field from Ergospheret = 10.7τSKerr black hole

ジェットの3次元安定性
PlasmaErgosphereMagnetic field lines

Magnetic field
flux tube acts aspropeller screw!

Koide （２００４）

3D stability• MHD jets are linearly unstable for helical KH
instability, but nonlinearly stable(Ouyed and Pudritz 2002)• MHD jets show non-axisymmetric structure
because of unstable character of accretion disk.
Neverthelss, the basic characteristics of jet
velocity and its parameter dependence are thesame as those in 2.5D (Kato. S et al., Kigure etal. 2005).

93D non-axisymmetric structure of a jet
produced by non-axisymmeric MHD
instability of the disk (Kigure-Shibata 2005)

MHDジェットの観測的証拠

Model R6

Spinning jets have finally been
observed by HST (Bacciotii et al. 2002)

Modeling Faraday rotation measureof AGN jet (Kigure et al. 2004)

0.5” = 110 AU
Optical jet fromDG TauCase of Collapsing cloud/starCase of collapsingmagnetized rotating star

• Supernova/collapsar– Lebranc-Wilson,Symbalicity– Proga et al.,Mizuno et al.

• Star formation– Tomisaka10

MHD simulation of collapsar
Mizuno et al. (2004) ApJ

中心星磁気圏と円盤の相互作用によるジェット形成

general relativistic MHD simulation withSchwarzchild black hole
Vjet ~ 0.2c

MHD model of protostellar jets as anextention of Hayashi et al (1996) model(Uehara et al. 2005)

Long term evolution of dipole caseUehara et al. 2005. to be submitted

1 orbit ＝ 6.28

Knots
Are formed
By sausage
Instability ?

Hayashi et al (1996)
A jet consists of two component:
reconnection outflow anddisk wind

Jet velocity in the case of dipole field(Uehara et al. 2005, submitted)
Maximum Velocity of Jet （Vz,max）Vz,max）
∝Bp(∝Bp1/3

This isdue to reconnectionrelated-accelerationof coronal outflow

for dipole field
for uniform field)
Vz,max / Vko

局所磁場をもつ円盤から噴出するジェット∝Bp∝Bp211

KatoY, Mineshige, Shibata (2004)
3D sim. (ApJ)

Case of Initially Localized weak field(Kudoh, Matsumoto, Shibata 2002, PASJ)
Magnetorotational Instability (Balbus and Hawley 1991)
leads to turbulence and reconnection
Magnetic Buoyancy(Parker instability) carries
magnetic flux into corona
at a rate of 0.1 – 0.3 VA
Similar to solar interior-corona
If external gas pressure is low,
Strong jet is formed

This toroidal field dominated jet is launched by magnetic pressure (similar toShibata and Uchida 1985, Turner et al. 1999, Kudoh et al. 2002), and is also
Similar to “magnetic tower” of Lynden-Bell (1996)

Long term evolution –nonsteady, intermittent ejection of jetswith many reconnection events
density

磁気回転不安定性の効果と長期時間変動•

Why do jets and disks can never reachsteady state ?
Because Magnetorotational Instability is so
powerful (Balbus and Hawley 1991)

Sato et al. 2005 I brahim et al. 2005
in prep

Magneto-rotational instability• Balbus-Hawley (1991)• [Chandrasekhar (1961), Velikhov (1959)]• Explains viscosity of accretion disks

Kudoh et al
2002 PASJ

12

Mass Accretion Rate vs Mass Ejection Rate

Power spectrum of mass ejection ratefrom MHD simulations(Sato et al, Ibrahim et al. .2005, to be submitted)

101
100
10−1

power

Mass Ejection Rate

Mass Accretion Rate

Observations of Time variabilityof X-rays from BH accretion disk
(Miyamoto, Kitamoto, Negoro)102

See poster

10−2
10−3
10−4
10−5

frequency
time

コリメーション
Tsinganos（２００３）

13

今後の課題: 加速と構造

今後の課題: コリメーション• ジェットの内部構造（ノット）の起源• ジェットの３次元安定性（＝＞磁気リコネクション＝＞粒子加速）• 相対論的ジェットの場合：
ローレンツ因子＞１０の高速ジェットはいかにして形成されるか？
＝＞ガンマ線バーストへの応用• 電子陽電子磁気的ジェットの物理？

• 回転星ダイポール磁場では、定常状態でどれだけの磁束が極方向にコリメートするか？
• 定常軸対称ＭＨＤアウトフローの数値解（今までのところ、全空間をセルフコンシステントに解いたのは、桜井(1985)のみ）＝＞数理物理の超難問

今後の課題： 観測• ＭＨＤジェットの観測的証拠をさがす– ジェットの回転とヘリカル磁場– ジェットの根本の降着円盤の磁場測定• ジェットのコリメーションのスケールを観測的に解明• ジェットとフレアの関係の検証

14