ブラックホール降着流・噴出流による フィードバック過程の 磁気流体シミュレーション

Slides:



Advertisements
Similar presentations
1 宇宙は何からできてくるか ? 理学部 物理 森川雅博 宇宙を満たす未知のエネルギー:暗黒エネル ギー 局在する見えない未知の物質:暗黒物質 銀河・星・ガス 何からできているか … 2006/7/25.
Advertisements

ジェット形成機構と粒子加速 機構 最近の話題 2012年9月6-7日 国立天文台 「巨大ブラックホールからの噴出流」研究会 高原文郎 ( 阪大理宇宙地球 )
硬 X 線で探るブラックホールと銀河の進化 深沢泰司(広大理) 最近の観測により、ブラックホールの形成と 銀河の進化(星生成)が密接に関係することが わかってきた。 ブラックホール観測の最も効率の良い硬 X 線で 銀河の進化を探ることを考える。 宇宙を構成する基本要素である銀河が、いつ どのように形成され、進化してきたか、は、宇宙の.
宇宙ジェット形成シミュレー ションの 可視化 宇宙物理学研究室 木村佳史 03S2015Z. 発表の流れ 1. 本研究の概要・目的・動機 2. モデルの仮定・設定と基礎方程式 3. シンクロトロン放射 1. 放射係数 2. 吸収係数 4. 輻射輸送方程式 5. 結果 6. まとめと今後の発展.
磁気流体降着円盤から噴出する磁気タワージェットの数値実験
YohkohからSolar-Bに向けての粒子加速
宇宙年齢10億年以前におけるSMBHの存在 遠方宇宙の観測で宇宙10億歳(z~6)未満で10億M⦿程度以上の活動銀河核中のSMBHの存在を確認 赤方偏移 z SMBH質量 [M⦿] URAS J ~2×109 M⦿ 宇宙7.5億歳(z~7)
星間物理学 講義3資料: 星間ガスの熱的安定性 星間ガスの力学的・熱的な不安定性についてまとめる。星形成や銀河形成を考える上での基礎。
高感度 VLBI 時代の QSO ターゲットを考えた
NGC 2043 銀河中 の 超光度X線源 (ULX) の スペクトル状態の遷移
(Fri) Astrophysics Laboratory MATSUO Kei
衝撃波によって星形成が誘発される場合に 原始星の進化が受ける影響
ブラックホール磁気圏研究会 ( 大阪市立大学)
加藤真理子1、藤本正樹2、井田茂1 1) 東京工業大学 2) JAXA/ISAS
原始惑星系円盤の形成と進化の理論 1. 導入:円盤の形成と進化とは? 2. 自己重力円盤の進化 3. 円盤内での固体物質の輸送
相対論的輻射流体力学における 速度依存変動エディントン因子 Velocity-Dependent Eddington Factor in Relativistic Photohydrodynamics 福江 純@大阪教育大学.
SHIGENOBU HIROSE AND JULIAN H. KROLIK
「Constraining the neutron star equation of state using XMM-Newton」
輻射優勢円盤のMHD数値実験 千葉大学宇宙物理学研究室 M2 松尾 圭 Thu.
GRS 等におけるジェット噴出と X 線強度変動の相関
ブラックホール降着・噴出流の 輻射磁気流体シミュレーション 大須賀 健 (国立天文台) -新型ジェットと母銀河への影響- 輻射加速 磁気収束
降着円盤磁気流体シミュレータの開発(2) 輻射磁気流体ソルバー
日本物理学会年次大会・総合パネル討論「現代プラズマ科学の 最前線:学際連携によるプラズマ理工学のさらなる展開」
Report from Tsukuba Group (From Galaxies to LSS)
降着円盤からの 相対論的輻射流体風 Relativistic Radiation Hydrodynamical Winds from Accretion Disks with Velocity-Dependent Eddington Factor Plane-Parallel Case 福江 純、秋月千鶴@大阪教育大学.
超巨大ブラックホール形成の鍵を握るAGNトーラスについて
ブラックホール降着流の大局的 3次元磁気流体数値実験
輻射圧駆動風の臨界点ついて 非相対論的領域 Radiatively Driven Spherical Wind Critical Points and Curves Nonrelativistic Regime 福江 純@大阪教育大学.
数値相対論の展望        柴田 大 (東大総合文化:1月から京大基研).
天文・宇宙分野2 松元亮治 「次世代スーパーコンピュータでせまる物質と宇宙の起源と構造」
S3: 恒星とブラックホール (上田、野上、加藤)
速度依存変動エディントン因子を用いた 相対論的鉛直輻射流 Velocity-Dependent Eddington Factor in Relativistic Photohydrodynamics Plane-Parallel Case 福江 純、秋月千鶴@大阪教育大学.
速度勾配依存 変動エディントン因子 Velocity-Gradient-Dependent Relativistic Variable Eddington Factor Plane-Parallel Case 福江 純@大阪教育大学.
ガンマ線バーストジェット内部における輻射輸送計算
Primordial Origin of Magnetic Fields in the Galaxy & Galaxies - Tight Link between GC and Cosmic B –  Y. Sofue1, M. Machida2, T. Kudoh3 (1. Kagoshima.
太陽・恒星フレアにおける輻射流体シミュレーション
明るいハードステートに対応する 光学的に薄い降着円盤モデル
論文紹介 Novae as a Class of Transient X-ray Sources K. Mukai, M
Fermi Bubble と銀河中心の巨大構造
ガンマ線バーストジェットの開き角 水田 晃(KEK) 井岡邦仁 (KEK) θj~CxΓ0-1 (C~1/5) ?
ブラックホール周辺の 磁場構造について 大阪市立大学 孝森 洋介 共同研究者 石原秀樹,木村匡志,中尾憲一(阪市大),柳哲文(京大基研)
「すざく」が NGC 4945 銀河中 に見付けた ブラックホール候補天体
銀河風による矮小銀河からの質量流出とダークマターハロー中心質量密度分布
磯部洋明 京都大学花山天文台 波動加熱勉強会 2004年2月23日
フォースフリーブラックホール 磁気圏 ~BH時空におけるGS方程式の解析~
基礎宇宙物理学 II 電磁流体力学入門 第1回 天体活動現象入門 2011年4月8日.
銀河・銀河系天文学 星間物理学 鹿児島大学宇宙コース 祖父江義明 .
S3: 恒星とブラックホール (上田、野上、加藤)
磁気回転不安定性によるブラックホール降着流の角運動量輸送機構の 解明
ON THE VERTICAL STRUCTURE OF RADIATION-DOMINATED ACCRETION DISKS N. J
Mitsuru Takeuchi and Shin Mineshige ApJ 486: ,1997 September 1
S5(理論宇宙物理学) 教 授 嶺重 慎 (ブラックホール)-4号館409 准教授 前田 啓一(超新星/物質循環)-4号館501
重要な効果 ブラックホールや中性子星(パルサー)磁気圏 銀河団スケールの加速(L×Bが大きい) 1020 eV以上
太陽フレアにおける Plasmoid-Induced-Reconnection の MHD シミュレーション
「すざく」であばく超光度X線源 (P4-7) rikne
新潟大学集中講義 ープラズマ物理学特論ー (天体電磁流体力学入門) 2004年1月19日ー1月21日
セイファート銀河中心核におけるAGNとスターバーストの結び付き
シミュレーションサマースクール課題 降着円盤とジェット
北大MMCセミナー 第62回 附属社会創造数学センター主催 Date: 2016年11月4日(金) 16:30~18:00
COE外国出張報告会 C0167 宇宙物理学教室 D2 木内 学 ascps
大阪市立大学 宇宙物理(重力)研究室 D2 孝森 洋介
ガンマ線バーストジェットの開き角 水田 晃(KEK) 井岡邦仁 (KEK) θj~CxΓ0-1 (C~1/5)
川島 朋尚 (国立天文台)、朝比奈 雄太 (国立天文台)、工藤祐己 (千葉大) supervised by 松本 洋介 (千葉大)
論文紹介07(2): ULXsの最近の論文から November 19, 2007 Tsunefumi Mizuno
スターバースト銀河NGC253の 電波スーパーバブルとX線放射の関係
大規模シミュレーションで見る宇宙初期から現在に至る星形成史の変遷
S5(理論宇宙物理学) 教 授 嶺重 慎 (ブラックホール)-4号館409 准教授 前田 啓一(超新星/物質循環)-4号館501
(Pop I+II連星起源と) 初代星連星起源 ロングガンマ線バースト
BH science for Astro-E2/HXD and NeXT mission
原始星からのX線発見と課題 (r-Ophの)T-Tauri星からX線放射とフレアーの発見
Presentation transcript:

ブラックホール降着流・噴出流による フィードバック過程の 磁気流体シミュレーション 松元亮治(千葉大理)

講演内容 共進化のマルチ連結シミュレーション 回転ガス円盤における角運動量輸送過程 降着円盤の大局的3次元磁気流体数値実験 円盤ダイナモとガス降着 降着円盤からのウィンド・ジェットの噴出 ジェット伝播シミュレーション 超臨界降着流のシミュレーション 輻射によるフィードバック

巨大ブラックホールと銀河の共進化のマルチ連結計算  巨大ブラックホールの起源 巨大ブラックホールへのフィーディング 銀河(バルジ) 銀河(10万光年)からブラック ホール(1光日)へのフィーディング 初代星形成 銀河(バルジ)形成 銀河核 宇宙論的進化 巨大ブラックホール 巨大ブラックホール形成 ジェットによる銀河への フィードバック 初代ブラックホール形成 電磁波放射による 銀河へのフィードバック 宇宙誕生後2-3億年 宇宙誕生後10-20億年 巨大ブラックホールからのフィードバック 4

流体系シミュレーションモデルと計算量 I 流体 演算量∝N3×Nstep 磁気流体 B 輻射流体 輻射磁気流体 磁気流体 N6×Nstep 計算領域を格子に分割して差分化 流体 3次元 ρ(t,x,y,z), v(t,x,y,z), P(t,x,y,z) 演算量∝N3×Nstep 磁気流体 輻射流束制限拡散(FLD)近似 +B(t,x,y,z) B 輻射輸送方程式を解く I 輻射流体 輻射磁気流体 磁気流体 N6×Nstep + I (t,x,y,z,n,q,f) +5年

ブラックホールへの質量供給

角運動量輸送問題 標準理論では粘性ストレス Trφ=αP と仮定 回転物質が落下するためには角運動量を失う必要がある。 降着のタイムスケール GAS SUPPLY 回転物質が落下するためには角運動量を失う必要がある。 標準理論では粘性ストレス Trφ=αP と仮定 降着のタイムスケール 観測と理論の比較から α ~ 0.01 – 0.1 流体乱流では α ~ O(0.001) too small ! H/r ~ 1 (RIAF)      0.1 (slim) 0.01 (thin) 降着円盤の理論モデルでは、長年、α粘性とよばれる現象論敵な粘性モデルが採用されてきました。これは、ストレステンソルの非対角成分が圧力に比例すると仮定し、その係数をαと置くというモデルです。観測との比較を通して、αの値は0.01から0.1程度と推定されてきましたがα粘性の起源は長年謎でした。 ガス温度に強く依存

差動回転円盤における磁気 回転不安定性 角運動量 重力 遠心力 R方向 Angular momentum 差動回転円盤における磁気  回転不安定性 角運動量 重力 遠心力 R方向 Angular momentum Balbus and Hawley (1991)

ブラックホール降着円盤の 大局的3次元MHDシミュレーション Machida and Matsumoto 2003

磁気エネルギーと角運動量輸送率の 時間変化 磁気エネルギー(対数) α= <-BrBφ/4πP0>  α~0.1 The megnetic energy increases with time and saturates when beta is aound 10. The ratio of the Maxwell stress to the gas pressure is about 0.1 in the innermost region. time TIME Machida and Matsumoto 2003

銀河ガス円盤の磁気流体シミュレーション (Nishikori et al. 2006) 重力ポテンシャル ダークマターを含む軸対称ポテンシャル( Miyamoto 1980) 初期条件 10kpcで密度最大になる角運動量一定のトーラス 方位角方向の弱い磁場 r=0.8kpc に吸収境界 赤道面対称性を仮定 Next, I would like to show the results of global 3D MHD simulations. We assume axisymmetric gravitational potential given by Miyamoto. In this potential, the effects of dark matter is included. Initial state is a constant angular momentum torus whose density maximum locates at 10Kpc. AT the initial state we assume weak toroidal magnetic fields. We include anomalous resistivity and Joule heating. We imposed absorbing BC at r=0.6Kpc. The right panel shows the density distribution and initial magnetic fields. We user 250*64*319 mesh points in cylindrical coordinate. 250*64*319 mesh

シミュレーション結果 2Gyr 3.5Gyr ρ+ B t = 3.8Gyr These panels show numerical results. Color shows density and curves show magnetic fields. As MRI grows, azimuthal magnetic field is amplified. The right panel shows the mean magnetic fields. The mean magnetic field is azimuthal but the strength of the turbulent magnetic field is comparable to the mean magnetic field.  ρ+ B t = 3.8Gyr

方位角方向磁場の逆転 10億年後

方位角方向磁場の逆転機構 +1 磁気回転不安定性が成長 +2 -1 パーカー不安定性による磁束流出 +2 浮上 -1

銀河中心への降着率 r=0.8kpcにおける降着率 =3.8Gyr This viewgraph shows the accretion rate to the central region. The accretion rate is 0.001 solar mass per year when the initial gas torus has mass about 10^8 solar mass. The accretion rate is small but high enough to support the activity of our galactic center. 初期にr=10kpcのトーラス状に n=1/cc のガスが分布している場合。より内側にガスが分布していれば降着率は増大。

銀河中心核ガス円盤の 大局的3次元MHDシミュレーション 2kpc Machida is also carrying out the simulation of the galactic central region. In this simulation, we imposed absorbing boundary condition at 200pc from the galactic center. The nuclear gas disk also becomes turbulent and mass accretion continues. This simulation produced magnetic loops inside 1kpc from the galactic center. Machida and Matsumoto (2006)

Global 3D MHD Simulation of the Galactic Center Gas Disk (Machida et al. 2008) Loop 銀河中心領域の模式図 磁気的加熱やコロナ形成によってH/rが増え、降着率が増大するのでは?

ブラックホール降着円盤とジェット 活動銀河中心核(AGN) ジェット 電波銀河M87中心核から噴出するジェット 1017cm VLA+ HALCA VLBA 43GHz Walker et al. 2007 活動銀河中心核(AGN) ジェット 電波銀河M87中心核から噴出するジェット  VLA(電波) 銀河中心核の活動性は超大質量ブラックホールへの質量降着によって駆動されていると考えられています。下図は降着円盤の模式図ですが、物質降着を可能にする角運動量輸送、円盤からのアウトフローの噴出等に磁場が本質的な寄与をしていることが明らかになってきました。さきほどの町田さんの講演でも紹介されましたように、私たちのグループでは降着円盤の時間発展を円盤全体を計算領域に含めたグローバルな磁気流体シミュレーションによって調べています。この講演ではSgrA*やこの図に示されたM87のように中心核が暗い低光度降着流の場合の計算結果について紹介します。 1017cm 連星系降着円盤(NASA) 銀河系内ジェット天体SS433

磁気流体ジェットの形成

降着円盤からのアウトフロー生成 磁力線と方位角方向磁場 vz=0.05c の面 This viewgraph shows the structurfe of the jet launchig region. Left panel shows the isosurface of vertical speed and solid curves in the right panel show magnetic field lines. Colors show azimuthal magnetic fields. We found that large-scale poloidal field is created and that outflows are accelerated along these field lines. 磁力線と方位角方向磁場 vz=0.05c の面

This viewgraph is a summary of our magnetic tower jet simulation This viewgraph is a summary of our magnetic tower jet simulation. When the central star is magnetically connected to the disk, well collimated magnetic tower develops. When the central object is a black hole and the disk has poloidal magnetic loop, the tower jet also appears.

ジェットの伝播 Bow Shock Jet-terminal Shock Internal Shock Cocoon 接触不連続面 Beam 密度の対数分布 等高線は圧力

Bφによるジェットの構造変化 密度のグラフ、等高線は圧力 Bφが加わりCocoonが細くなる ジェット先端の接触不連続面とBow Shockの間が薄くなる ジェットの伝播がBφを入れたほうが早い H1 MH-φ

銀河ジェットへの適用 コクーン 密度 温度 Gaibler et al. 2009

ジェット伝播シミュレーションの課題 非一様な密度分布中の磁気流体ジェットの伝播を3次元シミュレーションによって調べる。 磁気流体ジェットとの相互作用による星間ガスの加熱を調べる。 相対論的MHDシミュレーション 円盤から注入される磁場構造への依存性 円盤の質量流出率への依存性

降着円盤の輻射流体シミュレーション 降着率 α=0.01 r = 5rg Slim Advection Radiation RIAF SADM 次に降着率が高い場合の輻射流体シミュレーションについて紹介します。 表面密度 光学的に薄い 光学的に厚い Abramowicz et al. 1995

近傍銀河中で発見された 超光度X線源(ULX) M82 optical このシミュレーション結果は、近傍銀河中で発見された恒星質量ブラックホールのエディントン光度より明るい超光度X線源を、恒星質量ブラックホールへの超臨界降着によって説明できることを示唆しています。 Chandra X-ray Image X-ray NASA/CXC/SAO/PSU/CMU 恒星質量ブラックホールのエディントン光度を超えるX線源が存在

臨界降着率の100倍程度の降着率の場合のシミュレーション結果です。上図は温度、した図は密度分布を示します。回転軸方向に0.2c程度の速度でプラズマが加速されています。軸方向から見たときの光度はエディントン光度を超えています

解いている方程式はこれらの式で、輻射から働く力、逆コンプトン散乱による冷却などが含められています。輻射流速はFlux Limited Diffusion 近似を用いて計算されています。磁場はまだ考慮されていません。

輻射スペクトル 輻射流体シミュレーション結果から求めたスペクトル (Kawashima et al. 2010) ULXの輻射スペクトル (Gladstone  et al. 2009)

ジェット天体との関係 L > エディントン光度 ULX 光学的に厚いアウトフロー スペクトルはソフト エディントン光度を超えない    SS433? VLA(電波) この図に示すように、ブラックホールへの超臨界降着流によってULXとSS433のようなジェット天体を統一的に解釈することができます。すなわち、回転軸方向から見た場合はULX、軸と角度をなす方向から見た場合は光度は低く、スペクトルもソフトになり、ジェットを伴うSS433のような天体として観測されるというモデルです。 1017cm SS433から噴出するジェット

輻射磁気流体シミュレーション (Ohsuga,Mineshige,Mori,Kato 2009) 基礎方程式 Continuity Equation・・・・・・・ Equation of Motion・・・・・・・・ Gas Energy Equation・・・・・・・ Radiation Energy Equation・・ So far, I study accretion disks using RHD simulations. But magnetic fields are not solved and alpha-viscosity is assumed. Recently, many researchers perform MHD simulations. However, the radiation transfer is also important. Therefore, we carry out the RMHD simulations. In the RMHD simulations, multi-dimensional effects, radiation transfer, magnetic fields, and viscosity can be physically treated without assumption. RMHD simulations are most powerful tool to study accretion disks. Maxwell’s Equations・・・・・・ 輻射に関する項 磁場に関する項

次世代降着円盤シミュレータの開発 プラットフォーム : CANS 降着円盤への適用 並列計算機 向き最適化 シミュレーション結果解析 並列計算機 向き最適化 シミュレーション結果解析 シミュレーションエンジン 基本課題   ライブラリ リーマンソルバの改訂 相対論的磁気流体 ウェブページ 輻射磁気流体 可視化 降着円盤への適用 私たちのグループでは、輻射が重要になる領域まで含めた降着円盤のシミュレーションを実施するためのシミュレータを開発中です。以前に私たちが作成したコードのエンジン部にあたりリーマンソルバを改定するとともに、輻射磁気流体モジュールを加え、円盤の状態遷移やジェット形成に適用することを目指しています。 光学的に  薄 い円盤の 時間発展 ハード・ソフト 状態遷移 光学的に厚い円盤の進化 相対論的ジェットの形成

輻射によるフィードバック  Urry and Padvani 1995

RHDによるNLRの理論モデル r(pc)-数密度(1/cc)の関係。30-40pcまで100-1000 /ccのガスがまきあげられている。 Susa, Ohsuga, Wada 中心のBH(2e6Msun)+半径20pcの円盤 ionized 中心BHの輻射+円盤内のSN+流体 r(pc)-数密度(1/cc)の関係。30-40pcまで100-1000 /ccのガスがまきあげられている。 →NLR? Face-on ionized edge-on 緑:60deg 赤: cosin 青: no radiation Walsh et.al. 2008

まとめ ブラックホールへのガス供給:降着時間は円盤の厚さ(温度)に依存する。H/r=0.1程度なら磁気乱流による角運動量輸送により、円盤が百回転するタイムスケールで降着が進む。 超臨界降着によるブラックホール成長が可能。回転軸方向から観測したときの光度はエディントン光度を超える。 星間ガス中を伝播するジェットはkpcスケールで乱流を生成し、星間ガスを加熱する。 輻射磁気流体シミュレータを開発中。

END