学際大規模情報基盤共同利用・共同研究拠点 平成 28 年度共同研究 最終報告書 2017 年 5 月

jh160007-NAJ

次世代降着円盤シミュレータの開発
廣瀬重信（国立研究開発法人海洋研究開発機構）
降着円盤は、重力天体の周りに形成される回転ガス円盤であり、降着ガスの重力エネル
ギーを磁気エネルギーや輻射エネルギーに変換するエンジンとして、様々な天体活動現
象を駆動している。本研究課題では、この降着円盤におけるエネルギー変換過程を第一
原理から理解するために、角運動量輸送を担う乱流や、輻射輸送、一般相対論的効果な
ど、多様な物理を統合的に扱うことのできる「次世代降着円盤シミュレータ」の開発を目指している。

1. 共同研究に関する情報

現行の「降着円盤シミュレータ」では、磁気流(1) 共同研究を実施した拠点名

体力学シミュレーションによって磁気乱流を取り東京大学扱うことで、現象論的な角運動量パラメータを導(2) 共同研究分野入することなく、降着円盤の時間発展を調べるこ■

超大規模数値計算系応用分野とを可能にしているが、その対象は非相対論的な(3) 参加研究者の役割分担

現象に限られていた。また、輻射輸送や熱伝導に■ 全体統括（廣瀬重信・松元亮治）

ついては（その時間スケールが、多くの場合、降■ 一般相対論的輻射磁気流体力学（大須賀着円盤の力学時間よりも極端に短くなるため）時健・高橋博之・朝比奈雄太）

間陰解法で解く必要があるが、その計算効率の向■ 輻射スペクトル（川島朋尚・町田真美）

上が課題となっている。そこで、本課題では、現■ 自己重力ポアソンソルバー/陰解法（廣瀬行の「降着円盤シミュレータ」を改訂し、その適重信・松本倫明・簑島敬・花輪知幸・中用可能範囲を拡大するとともに、計算性能を向上

島研吾・伊田明弘・松本正晴・河合直聡）

させた「次世代降着円盤シミュレータ」を開発す■ コード整備（松本洋介・工藤祐己）

ることを目的とする。
（意義）

2.

研究の目的と意義

降着円盤は、近接連星系や、活動銀河中心核、（目的）

星形成領域等で観測される激しい時間変動や X 線

降着円盤は、重力天体の周りに形成される回転

放射、ジェット噴出などを伴う様々な活動現象をガス円盤であり、降着ガスの重力エネルギーを磁駆動している。また、太陽系をはじめとする惑星

気エネルギーや輻射エネルギーに変換するエンジ系形成の母体となる原始惑星系円盤も降着円盤のンとして、様々な天体活動現象を駆動している。

一種である。したがって、降着円盤の構造や安定

我々のグループでは、この過程を第一原理からシ性、時間変動、輻射スペクトル、ジェット噴出機

ミュレーションするために、磁気乱流による角運構等を解明することは、理論天体物理学の中心テ

動量輸送や、物質と輻射の相互作用（運動量交換ーマの一つである。

とエネルギー交換）、高温領域と低温領域間での熱本課題を通して、輻射輸送、熱伝導、自己重力、伝導、星形成や活動銀河核においては自己重力、

一般相対論効果などを取り入れた「次世代降着円

ブラックホール近傍においては一般相対論的効果、 盤シミュレータ」が構築されることにより、様々など、多様な物理を統合的に扱うことのできる「降な重力天体の周囲の降着円盤シミュレーションに着円盤シミュレータ」の開発を目指している。

対応できるほか、降着円盤のダイナミクスと無撞
1

学際大規模情報基盤共同利用・共同研究拠点 平成 28 年度共同研究 最終報告書 2017 年 5 月

着な輻射場が得られることになり、シミュレーシ

4. 前年度までに得られた研究成果の概要

ョン結果を観測結果と直接比較することが可能にA）一般相対論的輻射磁気流体力学コードの開発なる。実際、こういったシミュレーション結果の我々はこれまで特殊相対論的輻射磁気流体コー「観測的可視化」は、本格稼働を始めているアルドを用いて超臨界（Eddington 降着率を超える）

マ望遠鏡など、精密化する観測を活かすためにも降着円盤からのアウトフロー加速について研究を必須である。

行ってきた。しかし、この計算ではブラックホー

現在、流体力学計算と組み合わせるために用い

ルによる重力を、ニュートン近似を用いて扱ってられている、角度モーメントを用いた輻射輸送解きた。従ってブラックホール周りの重力による光法は、クロージャ問題という本質的な問題を抱え

線の曲がりや、ブラックホールの回転といった効ている。これに対し、本課題の開発項目の一つで果を無視していた。また、計算は軸対称を仮定しある、輻射輸送方程式の直接解法はこの問題を解ていたため、円盤内のダイナモ機構を追う事が出決する最終手段であるが、基本的に（空間３次元・来ないといった課題が残っていた。

光線方向２次元・振動数１次元の）６次元計算とこれらの問題を解決すべく我々は３次元一般相なり、その計算量は膨大になる。
「京」コンピュー

対論的輻射磁気流体コードを構築した。基礎方程タを用いても本格的な計算は難しいが、ポスト「京」 式は以下である。
においては、むしろ本命の計算課題となることが𝜕!

−𝑔𝜌𝑢 ! + 𝜕!

!"

+ 𝜕!

−𝑔𝑇! =!

!
−𝑔𝑇!! Γ!"+ −𝑔𝐺!

!"+ 𝜕!

−𝑔𝑅! =!

!
−𝑔𝑅!! Γ!"− −𝑔𝐺!

予想される。したがって、本課題を通してそのプ

𝜕!

−𝑔𝑇!

ロトタイプを開発しておくことは、ポスト「京」

𝜕!

−𝑔𝑅!
𝜕!

に向けて、重要なステップとなる。

−𝑔𝜌𝑢 ! = 0−𝑔𝐵 ! + 𝜕! [ −𝑔(𝑏 ! 𝑢 ! − 𝑏 ! 𝑢 ! )] = 0

上から順に質量保存、磁気流体・輻射のエネル!
ギー運動量保存式、誘導方程式である。𝑔!" , Γ!"は

3. 当拠点公募型共同研究として実施した意義
降着円盤のシミュレーションでは、典型的に回メトリックの行列式とクリストッフェル記号で重転周期の数百倍の時間発展を追う必要があり、そ力やブラックホールの回転の効果を表す。

のための時間ステップは数百万ステップに及ぶ。

輻射場は輻射強度を運動量空間で積分したモーそのため、上記の多様な物理プロセスを取り入れ

メント量 によって定式化され、M1 クロージャー

た、本格的な３次元プロダクトランには「京」コを用いることで方程式を閉じる。輻射とガスは吸ンピュータでの計算が必須となる。そこで、東京

収・散乱過程により輻射４元力 を通じてエネルギ

大学情報基盤センターの FX10 を用いて、「京」のー・運動量のやり取りを行う。従ってこれらの式利用に向けたコードの開発・最適化およびテストを解くことでガスと輻射を無矛盾に決定する事が計算を行う。

出来る。本計算ではガスと輻射の相互作用としてまた、上述のように、輻射輸送や熱伝導を陰解

電子散乱と自由—自由吸収を考慮した。

法で解く場合や、自己重力を考慮する場合は、行

このコードを用いて計算を行った結果、降着円

列反転を高速に行うことが課題である。現在は、

盤内縁半径が従来の理論的考察よりも外側にある

ガウスザイデル法＋マルチグリッド法の組み合わこと、ブラックホールが高速回転している場合にせを適用しているが、必ずしも効率が良いとは言

円盤内縁の温度が無回転の場合よりも高温になるえない。そこで、東京大学情報基盤センターの中ことなど、一般相対論的効果が重要となる、新し島研吾氏のグループと協力し、マルチグリッド法

い結果が得られている。

の代わりに、クリロフ部分空間法を適用することによりコードの高速化が可能か検証を行う。

B)角度モーメントを用いない輻射輸送コードの開発

A)の一般相対論的輻射磁気流体シミュレーショ
2

学際大規模情報基盤共同利用・共同研究拠点 平成 28 年度共同研究 最終報告書 2017 年 5 月

ンにより得られたブラックホール超臨界降着流の輻射スペクトルを計算するために、モンテカルロ
法に基づく輻射スペクトル計算コード開発を進めている。前年度には 5 次および 8 次の embedded
Runge-Kutta 法（adaptive stepsize control アル
ゴリズムにより、積分のステップ幅を変化させな
がら計算する Runge-Kutta 法） を用いた測地線
方程式ソルバーを作成した。さらに、磁場の方向を考慮したシンクロトロン放射・吸収の計算モジ

図 1 中性子星(左)とブラックホール(右)へのガス降着の様ュールを作成した。

子 。カ ラー は ガ ス密 度 ( 左) と 輻射 エ ネ ルギ ー 密度 (右 ) 、 線はガス(左)と輻射(右)の流線、白線は光 球面を示す。

C）行列反転コード（熱伝導・輻射拡散などの陰解法

射磁気流体コードを用いるが、ガスと輻射の相互や自己重力ソルバ）の高効率化作用として電子散乱と自由—自由吸収の他に熱的昨年度は、東京大学情報基盤センターの中島研Compton 散乱も考慮した。

吾氏らによって開発された差分法向けに最適化さ

図 1 にその結果を示す。中性子星近傍でのガスれた行列解法ライブラリを、熱伝導問題に適用す

降着率は Eddington 降着率の 100 倍程度であり、るための基本的なコード開発を行った。一方で、中性子星へのガス降着が可能であることがわかる。

熱伝導のタイムスケールが降着円盤の力学時間にしかし中心天体近傍を見ると中性子星とブラック

比べてそれほど短くない場合には、Super-time

ホールで結果が異なる。ブラックホールの場合、

stepping 法がむしろ有効であることを確かめた。

ガスや輻射エネルギーは中心天体へと吸い込まれるため、ブラックホールへの流入口とアウトフロ

5. 今年度の研究成果の詳細

ーによる放出口の場所がはっきりと異なる。一方A）一般相対論的輻射磁気流体力学コードの開発

で中性子星の場合はガスやエネルギーを吸い込む

昨年度に引き続き、実用問題を通して、コードことができないため、中性子星近傍にガスと質量の適用範囲拡大を進めている。以下に、本年度のが溜まっているのがわかる。この貯められたエネ

科学的成果をまとめる。

ルギーによって非常に強いアウトフローを作るこ

昨年度までの計算では超臨界降着円盤の中心天とがわかった。200 重力半径でのパワーを見積も体としてブラックホールを想定していた。しかし2014 年に、これまでブラックホール超臨界降着円盤が中心エンジンであると考えられてきた超高光度 X 線源(ULX)の１つが中性子星である(ULX パルサー)ことがわかった。つまり、ブラックホールと異なり硬い表面を持つ中性子星でも大量のガス降着が可能であることが示された。
そこで、今年度前期では、我々は中性子星への超臨界降着が本当に可能か、そして降着に付随するアウトフローや光度はブラックホールと異なる

図 2 双極磁場を持つ中 性子星への超臨界降着の様子。カラーはのかを調べた。この目的のため、まずは簡単のた

密度(左)と輻射エネルギー (右)、矢印は速度場(左)と輻射フラッ

めに中性子星の磁場は十分に弱く、無視できると

クス(右)、青線は磁力線を 示す。

仮定した。計算は昨年度と同様、一般相対論的輻3

学際大規模情報基盤共同利用・共同研究拠点 平成 28 年度共同研究 最終報告書 2017 年 5 月

ると、ブラックホールでは輻射による光度が最も大きく、Eddington 光度を超えるが、運動エネル
ギーのパワーは Eddington 光度を下回る。一方で中性子星では運動エネルギーによるパワーが最も大きく、その次に輻射光度が大きい。このように同程度の降着率でもアウトフローの組成はブラッ
クホールと中性子星で異なることがわかった。さらに全アウトフロー強度(輻射＋運動エネルギー)
も中性子星の方が大きい。したがって観測される
ULX のある程度の割合は中性子星で説明することが可能であることがわかった。
本年度後期では、強磁場中性子星を想定し、中性子星の持つ双極磁場と円盤の相互作用を調べる

図 3 中性子星 への柱状の超臨界 降着の様子。(a)質量密度

ために同様の２次元計算を行った。前期の研究と分布（ 色）とガス速度の向き (矢印 )。(b)輻射エネルギー密の大きな違いは中性子星が 1010Gauss の双極磁場度分布（色）と輻 射フラックスの向き(矢印)。

を持つこと、また、強磁場中性子星を考えている

s-1 となり、これは観測される ULX パルサーと同程ため、シンクロトロン放射を考慮したことである。

度であることがわかった。

この結果を図 2 に示す。

また、本年度前期では、上記の研究で示されるよ

中性子星より遠方では双極磁場のエネルギーはうな中性子星極冠へと落下するガスの輻射流体シ

非常に小さいため、弱磁場中性子星の場合同様、

ミュレーションを実施している(図 3)。その結果、超臨界降着円盤が形成される。しかし中性子星近輻射は柱状の超臨界降着流の側面から抜け出し、傍（R=3R*, R*=10km は中性子星半径）では円盤が途

Eddington 光度のおよそ 100 倍に到達できること切れていることがわかる。これは双極磁場によるを明らかにした。これはパルスを示す ULX の 1 つ

磁気圧が卓越することで、円盤による降着が遮らである M82 X-2 の光度をよく再現している。このれるためである。従って円盤が途切れる半径=3R*

輻射流体シミュレーションによって「中性子星へ

は円盤の輻射圧と双極磁場の磁気圧が釣り合う半の柱状の超臨界降着が ULX の莫大な光度とパルス径で決まる。この半径より内側では磁場が卓越す

放射の両方を説明可能である」ことを示し、プレる磁気圏が形成される。

スリリースを行った。

磁気圏より遠方では磁気回転不安定性によって以上では、降着ガスは柱状の領域内を自由に動角運動量が外側へ輸送され、その結果ガス降着がくことが可能であると想定してが、さらに磁場が起こる。しかし磁気圏半径付近では磁場が強く、強い状況では完全に双極磁場に沿った方向にしか磁気回転不安定性は働かない。しかしこの半径でガスの運動が許されないと考えられる。そこで、は双極磁場による磁気トルクが非常に効率的に働

本年度後期では、ガスの運動方向を制限した状況

くために角運動量は内側へと輸送され、ガス降着での輻射流体シミュレーションへとコードの改良が起こる。その結果、磁気圏内ではガスは磁力線を行った。プレリミナリーな結果は得られつつあに沿って中性子星の磁極に向かって降着し、中性り、柱状であった構造が大幅に変わることがわか子星表面に降着柱を形成する。また、双極磁場にった。ガスは柱の外壁に沿って降着するものの、よって引き抜かれた角運動量は最終的に中性子星

柱内部のガスが吹き飛ばされるため、中の抜けたへと渡され、中性子星のスピンアップに繋がる。

筒状の構造になることがわかった。この場合でも本計算で見積もられたスピンアップ率は 10-11 s
4

学際大規模情報基盤共同利用・共同研究拠点 平成 28 年度共同研究 最終報告書 2017 年 5 月

光度はエディントン光度を超え、ULX パルサーの効と認められないという結果になった。

モデルとなり得ることがわかった。

今年度後期では、

C）行列反転コード（熱伝導・輻射拡散などの陰解法

n

BiCGSTAB 法の代わりに GCR 法を用いる

n

BiCGSTAB 法の前処理に GCR 法を用いるという方法も試みたが、いずれも、外側の非線形や自己重力ソルバ）の高効率化

本年度は、特に輻射輸送問題の陰解法に対して、問題（ニュートン法）の反復がある程度進むと、これまでのフルマルチグリッド法の代わりに、前内側の線形問題（GCR 法または BiCGSTAB 法）が収処理付き BiCGSTAB 法の適用を試みた。具体的には、 束しなくなる、という結果になった。
・ガスエネルギー式（非線形項）
・輻射エネルギー式（非線形項＋線形拡散項）6.

今年度の進捗状況と今後の展望の二本の方程式を連立させて、時間陰的に解く。

A）一般相対論的輻射磁気流体力学コードの開発

つまり、非線形方程式を解くニュートン・ラプソ

2014 年に中性子星への超臨界降着現象が観測

ン法の反復ループの内側で、線形問題（拡散方程的に指摘され、2016 年には同様の現象が他の２つ式の陰解法）を解くことになる。

の天体でも観測された。これら ULX パルサーとこの線形問題を解くのに、これまでフルマルチ呼ばれる天体はコンパクト天体の観測分野で非常グリッド法のサイクルを１〜２回決め打ちで繰りにホットな話題となっている。しかしその一方で、返し、その結果を外側のニュートン・ラプソン法

理論は進んでおらず、特にシミュレーション研究の反復ループに返していた。したがって、必ずしはその難しさから全く研究がなされてこなかった。

も線形問題を厳密に解いていたわけではなかった

今年度に行った双極磁場を持つ中性子星への超臨が、結果的に、（最終的に解きたい）非線形問題の界降着の一般相対論的輻射磁気流体シミュレーシ

ニュートン・ラプソン法のループは数回〜数十回

ョンは理論研究のマイルストーンとなる研究と考の反復で収束していた。

えている。

今年度前期では、このフルマルチグリッド法をその一方でこの研究では回転軸対称を仮定してそのまま、前処理付き BiCGSTAB 法に置き換えた。

いること、磁場強度が観測から示唆されている値

BiCGSTAB 法では前処理が重要であるが、まず手始よりも 2-3 桁小さい、という問題がある。前者のめに、簡単のため、SOR 法あるいはガウスザイデ

仮定は円盤降着の場合は比較的よい近似となりう

ル法（決め打ちで 20 回反復）を前処理に用いるこるが、３次元効果を考慮した場合、磁気圏半径でとにした。

非軸対称不安定性が成長するため、今後詳細な研このように作成したコードを用いて、具体的な究が必要である。また磁場強度についてもより広問題として、原始惑星系円盤における輻射輸送問

いパラメータ空間で計算し、観測との比較が必要

題を解いてみた結果（テスト計算のため手元のパとなる。

ソコンで 1 コアを用いたもの）を以下の表に示す。
このテスト結果からは、特に今回の問題（ニューB）角度モーメントを用いない輻射輸送コードの開発

トン・ラプソン法の反復ループ内の近似解を求め観測との比較においては、一般相対論的輻射輸る問題）に対しては、前処理付き BiCGSTAB 法は有送計算コードの開発を進めている。歪んだ時空中線形計算 L2 ノルム

BiCGSTAB 法反復数

Newton 法反復数

計算時間(秒)

Full MG

0.1-0.5（FMG サイクル 1 回） N/A

15

15

BiCGSTAB

0.000180-100

9

400

0.5

10-40

18

164

5

学際大規模情報基盤共同利用・共同研究拠点 平成 28 年度共同研究 最終報告書 2017 年 5 月

での光の伝搬および電子と光子の散乱のテスト計までニュートン・ラプソン法の反復ループ内で有算についてはすでに成功している。今後は実用化

効な近似解を求める手段に過ぎない、という点がに向けたチューニングを行う計画である。

ポイントであると思われる。今後の展望としては、

加えて、これまで採用してきた近似法である

今年度は完成に至らず試すことができなかったが、M1-closure 法を変動エディントン因子法へと発マルチグリッド法を前処理とした BiCGSTAB 法、お展させることで、より精緻な輻射流体シミュレー

よび、ニュートン・クリロフ法の有効性をぜひ確

ションを行うことも必要である。この新たなコー

かめたいと考えている。

ドの開発も進める計画である。

7.

研究成果リストC）行列反転コード（熱伝導・輻射拡散などの陰解法(1) 学術論文

や自己重力ソルバ）の高効率化-

Kawashima, T., Matsumoto, Y., Matsumoto,

本項目では、昨年度まで用いてきた、フルマル

R.,‘A Possible Time-Delayed Brightening

チグリッド法の代わりに、前処理付き BiCGSTAB 法

of the Sgr A* Accretion Flow after theを用いる方法を試すことができた。これまで、本Pericenter

項目の進捗の障壁となったのは、既存の行列反転

Publications of the Astronomical Society ofライブラリを用いるためには、線形拡散方程式のJapan, in press

拡散係数から行列要素データを作成しなければな-

PassageofG2

Cloud’,

Ogawa, T., Mineshige, S., Kawashima T.,

らない点である。特に、ここで対象としている問Ohsuga,

題はデカルト座標を用いているが、そのうち x 方hydrodynamic

向にずれ周期境界を用いているため、行列が複雑

super-Eddington accretor as a model for

な非対称行列となり、かつ、その「ずれ」の量がultra-luminous sources’, Publications of

毎時間ステップ異なるため、行列要素データの作

the Astronomical Society of Japan, in press

成が非常に煩雑になる。さらに、ガウス・ザイデ-

K.,Hashizume,

K.,‘Radiation

simulationsofaNomura, M., Ohsuga, K., ‘Line-driven disc

ル法を用いてきたユーザー側が、行列要素データ

wind model for ultrafast outflows in activeの取扱いに慣れていないという点も、実際上大きgalactic nuclei - scaling with luminosity’,

な問題であった。

2017,

今年度は、東京大学情報基盤センターの河合直Monthly

Notices

of

the

Royal

Astronomical Society, 465, 2873

聡氏の監修のもと、BiCGSTAB 法のアルゴリズムで-Takahashi, H.R., Ohsuga, K.,’ Formation ofは、行列が常に「行列とベクトルの積」で現れる

Overheated Regions and Truncated Disks

ことに着目し、その部分の計算を既存のガウス・aroundザイデル法のコード（の再利用）で実行するプロ

General

グラムを作成した。実際、こうすることで、ユー

Radiation-magnetohydrodynamics

ザー側からは、行列を意識せずに BiCGSTAB 法を用

Simulations’,2016,

いることが出来る他、ずれ周期境界を含め、複雑

Journal, 826, 23

な境界条件にも柔軟に対応できるというメリット-がある。
ただし、前節で示した通り、テスト計算の結果

Black

Holes:

Three-dimensional
Relativistic
The

AstrophysicalKawashima T., Mineshige S., Ohsuga K.,
Ogawa

T.,

model

of

‘A

radiation-hydrodynamics

accretionforpulsars’,

2016,

では、今回作成した、前処理付き BiCGSTAB 法の有ultra-luminous

効性は認められなかった。その理由は今後の検討

Publications of the Astronomical Society of

課題であるが、今回の問題では、行列反転があく

Japan, 68, 83
6

X-ray

columns

学際大規模情報基盤共同利用・共同研究拠点 平成 28 年度共同研究 最終報告書 2017 年 5 月

-

Dwarf

nova

outbursts

with

March 3-4, 2017 (招待講演)

magnetorotational turbulence, M. S. B.

--

Ultra-Luminous X-ray Source 研究会 - 多波(4): 3710-3726.

長で探る降着系の統一描像、口頭発表,

Convective quenching of field reversals in

Kanagawa (Japan), March 3-4, 2017 (招待講accretion disc dynamos, Matthew S. B.

演)
-

for

ULXs

and

Hydrodynamic

S. Hirose, Mon Not R Astron Soc (2016) 462

ULX-pulsars’,

Kawashima T., ‘A possible time-delayed

Salvesen, Shigenobu Hirose, Mon Not R

brightening of the Sgr A* accretion flow

Astron Soc (2017) 467 (3): 2625-2635.

after

The contribution of bulk Comptonization to

Challenges of AGN Jets 、 口 頭 発 表 , Tokyo

the soft X-ray excess in AGN, J. Kaufman,(Japan), January 17-20, 2017 (招待講演)-(2017) 467 (2): 1734-1747.
Gravito-turbulence

pericenterin

irradiatedpassageofG2’,

Ohsuga,K. ‘Theory of ultra-luminous X-Ray
sources’,Breaking

Super-Eddington

theAccretion

Limits:onCompact

protoplanetary discs, Shigenobu Hirose,

Objects, Sardinia (Italy), September 19-23,

Ji-Ming Shi, Mon Not R Astron Soc stx824.

2016 (招待講演)

DOI: https://doi.org/10.1093/mnras/stx824-‘Radiation

ModelsO. M. Blaes, S. Hirose, Mon Not R Astron Soc-

T.,

Coleman, I. Kotko, O. Blaes, J.-P. Lasota,

Coleman, Evan Yerger, Omer Blaes, Greg-

Kawashima-

Takahashi,H.R.‘Radiation

Matsumoto, T., Machida, M. N, Inutsuka,

Magnetohydrodynamic

S-I., 'Circumstellar disks and outflows in

Accretion Disks’, East Asian Numericalturbulent molecular cloud cores: possible

Astrophysics

formationmechanism

October 24-28, 2016 (招待講演)

systems',

The

for

Astrophysicalmisaligned
Journal,-839:69 (14pp), 2017 April 10

Meeting,

Simulations
Beijingof(China),Takahashi H.R. ‘Accretion and outflows
observed in GRRMHD Simulations’, ULXs andtheir environments, Strasbourg (France),(2) 国際会議プロシーディングス

June 13-16 (口頭発表)

N/A-

Takahashi

H.R"Relativistic

Magnetohydrodynamic

Radiation

Simulationsof(3) 会議発表(口頭，ポスター等)

Accretion Disks",International Congress on

国際会議発表Plasma Physics, Takao (Taiwan), June 27--July 1, 2016 (招待講演)

Takahashi H. R., ‘General Relativistic
Radiation Magnetohydrodynamic Simulations--"Thermal equilibrium curves obtained byof Disk Accretion onto a Black Hole or aradiation

Neutron Star’, The US-Japan Workshop andHirose, in Simulations and Modelling of

School on Magnetic Reconnection, Ehime

Relativistic MHD Accretion Discs, Oxford(Japan), March, 19-23, 2017 (招待講演)

University, July 10-15, 2016（招待講演）

Takahashi H. R., ‘Black Hole and NeutronStar

Accretion

Disk’,-

MHD

simulations”,"Gravito-turbulencein

Shigenobu

irradiatedUltra-Luminous

protoplanetary disks”, Shigenobu Hirose,

X-ray Source 研究会 - 多波長で探る降着系

in “Confronting MHD Theories of Accretionの統一描像、口頭発表, Kanagawa (Japan),

Disks with Observations”, Kavli Institute7

学際大規模情報基盤共同利用・共同研究拠点 平成 28 年度共同研究 最終報告書 2017 年 5 月

for Theoretical Physics, Jan 9 - Mar 30,-

simulations

2017（招待講演）-

ジウム第 29 回 理論懇シンポジウム「重力がin Astrophysical Objects', The US-Japan

織りなす宇宙の諸階層」、東北大学, 2016 年 12

Workshop

月 20-22 日(ポスター発表)

andSchoolMatsuyamaon

Magnetic(Japan),

March

-

in

研究推進連絡会」第 4 回ワークショップ、東Accretion Disks', 6th East-Asia School and

京大学天文学教育研究センター, 2016 年 12 月

Workshop

９-10 日(招待講演)

Matsumoto,R.,'Plasmaon

Dynamics

Laboratory,

Space,-

ofール研究推進連絡会」第 4 回ワークショップ、Protostellar Binary and Multiple Systems

東京大学天文学教育研究センター, 2016 年 12with

月９-10 日(口頭発表)

Matsumoto,

T.,

AMR

'Theoretical

Simulations',

ModelsEast-Asian-

よる降着円盤とアウトフロー研究’, 高エネ

講演）ル ギ ー 宇 宙 物 理 学 研 究 会 2016, 青 山 大 学 ,Matsumoto,
with

T.,

AMR

'Theoretical

Simulations',

Modelsof2016 年 11 月 30 日-12 月 2 日(招待講演)
高 橋 博 之 ‘Super Critical Accretion Onto
Weakly Magnetized Neutron Star’, Workshop

Monterey (USA), June 9 2016 （招待講演）

on Theoretical Astronomy 2016, 富山大学,

Matsumoto, T., 'A dynamical model of the2016 年 11 月 14 日-16 日(口頭発表)

with

the

ASTRONUM-

2016,

heliosphere

adaptivemesh-高橋博之 ‘降着円盤の一般相対論的輻射磁

refinement', JpGU2016, Makuhari (Japan),

気流体シミュレーション’, CfCA ユーザーズ

May 23 2016 （口頭発表）

ミーティング, 国立天文台, 2016 年 11 月 29

国内会議発表日-30 日(口頭発表)

高橋博之 ‘双極子磁場を持つ中性子星への超-高橋博之 ‘降着円盤の輻射磁気流体シミュ

臨界降着の一般相対論的輻射磁気流体シミュ

レーション’, 日本物理学会 2016 秋季年会、

レーション’, 日本天文学会春季年会, 九州大

金沢大学、9 月 13 日-15 日(招待講演)

学, 9 月 15 日-19 日(口頭発表)-

高橋博之‘輻射磁気流体シミュレーションにBeijing (China), October 24, 2016 （招待Protostellar Binary and Multiple Systems-

川島朋尚 ‘中性子星極冠への超臨界降着の輻
射流体シミュレーション’,「超巨大ブラックホ

Numerical Astrophysics Meetings (EANAM7),-

大須賀健 ‘超臨界降着による巨大ブラックホールの成長について’,「超巨大ブラックホール

July 11-July 16, 2016 (招待講演)-

accretion

Amplification and Magnetic Energy Release'Magnetic

Astrophysical Plasmas, Tsukuba (Japan),-super-criticalcolumns onto neutron stars’, 理論懇シンポR.,

19-March 23, 2017 (招待講演)

-

of

Field

Matsumoto,

Reconnection,-‘Radiation-hydrodynamic

川 島 朋 尚-高橋博之’弱磁場中性子星への超臨界降着の川島朋尚 ‘強磁場中性子星の超臨界降着柱の一般相対論的輻射磁気流体シミュレーショ

構造’, 日本天文学会春季年会, 九州大学, 9 月

ン’,流体力学会年会 2016, 名古屋工業大学、

15 日-19 日(口頭発表)

9 月 26 日-28 日(口頭発表)

高橋博之 ‘ブラックホール・中性子星降着円盤-

朝比奈雄太, 野村真理子, 大須賀健, ‘活動銀とアウトフロー’, 理論懇シンポジウム第 29 回

河核アウトフローによる星間ガスへのフィー理論懇シンポジウム「重力が織りなす宇宙のドバックの磁気流体計算’, 日本天文学会

諸階層」、東北大学, 2016 年 12 月 20-22 日 (招2016 年秋季年会, 愛媛大学, 2016 年 9 月 15

待講演)日(口頭発表)8

学際大規模情報基盤共同利用・共同研究拠点 平成 28 年度共同研究 最終報告書 2017 年 5 月-

朝比奈雄太, 野村真理子, 大須賀健, ‘宇宙ジ
ェットによる非一様な星間空間へのフィード
バックの磁気流体シミュレーション’, 日本流
体力学会年会 2016, 名古屋工業大学, 2016 年
9 月 28 日(口頭発表)-

川島朋尚, 嶺重慎, 大須賀健, 小川拓未
‘強磁場中性子星極冠への超臨界降着柱の輻
射流体シミュレーション：ULX パルサーへの応
用’, 日本天文学会 2016 年秋季年会, 愛媛大学, 2016 年 9 月 16 日(口頭発表)-

川島朋尚, 嶺重慎, 大須賀健, 小川拓未
‘強磁場中性子星への超臨界降着柱の輻射流
体シミュレーション’, 日本流体力学会年会
2016, 名古屋工業大学, 2016 年 9 月 28 日(口頭発表)-"中心星の照射を受ける原始惑星系円盤における重力乱流", 廣瀬重信, Ji-Ming Shi, 日本天文学会秋季年会, 愛媛大学, 2016 年 9 月 16 日
（口頭発表）-"矮新星降着円盤の「第一原理」熱平衡曲線",
廣瀬重信, 連星・変光星・低温度星 研究会
2016, 慶應義塾大学, 2016 年 10 月 31 日（口頭発表）(4) その他（特許，プレス発表，著書等）-

川島朋尚, 嶺重慎, 大須賀健‘横顔が輝く宇宙灯台：謎の超高輝度 X 線パルサーの正体をスパコンがあばく’2016 年 9 月 8 日 (プレスリリース)

9