Measuring the Black Hole and Accretion Parameters of Sagittarius A* from EHT Observations using a Semi-Analytic Model

技術要約：半解析的モデルを用いたEHT観測によるいて座A（Sgr A）のブラックホールおよび降着パラメータの測定**

問題提起
イベント・ホライズン・テレスコープ（EHT）コラボレーションは、天の川銀河の中心にある超大質量ブラックホール、いて座A*（Sgr A*）の見かけ上のシャドウの初となる画像を作成した。しかし、Sgr A*は、観測時間（分から時間単位）よりも短い時間スケールで顕著な構造的可変性を示す。この急速な日内変動のため、地球自転による開口合成を用いて$(u, v)$平面のサンプリングを改善し、光源構造を正確に再構成することが困難となっている。その結果、標準的な撮像技術では、時間とともに変化する放射特性と、強重力および降着環境に由来する安定したシグネチャを分離することに苦慮している。単純な幾何学的モデルを用いた従来の解析では、リング径は制約されたものの、輝度非対称性や物理的な降着条件といった他の形態学的パラメータについては制約が不十分なまま残されていた。さらに、一般相対論的磁気流体力学（GRMHD）シミュレーションは解釈に不可欠であるが、その計算負荷の高さから、疎な可視度データに対する直接的なベイズ推論には適していない。

手法
これらの課題に対処するため、著者らは「スナップショット・モデリング」アプローチとベイズ階層フレームワークを併用している。核となる手法は以下の通りである：

1. 半解析的放射モデル： 本研究では、軽量で堅牢な現象論的「デュアルコーン（二重円錐）」放射モデルを用いている。このモデルは、ブラックホールのスピン軸に沿った、互いに反対方向を向いた2つの薄い円錐内に流れる、シンクロトロン放射を行う磁化プラズマを仮定している。モデルは軸対称であり、カー計量（Kerr geometry）において光線追跡（レイトレーシング）が行われる。含まれるパラメータは、ブラックホールスピン（$a_*$）、観測者傾斜角（$\theta_o$）、円錐開口角（$\theta_s$）、投影されたスピン位置角（p.a.）、特性放射半径（$R_{peak}$）、放射プロファイル指数（$p_1, p_2$）、流体速度（$\beta_v$）、および磁場構成（$\iota, \eta, \chi$）である。質量対距離比（$\theta_g$）は、事前の制約に基づき $5.03\,\mu\text{as}$ に固定されている。
2. スナップショット推論： 2017年4月7日のEHTデータを120秒ごとの「スナップショット」に分割した。各スナップショットについて、著者らはComage.jlパッケージを用いて、複雑な可視度に対するベイズ推論を実行した。これにより、物理モデルのパラメータと、計測器誤差を考慮するための複素ステーションゲイン（振幅と位相）を同時に推論している。
3. サンプリング戦略： 高次元性と潜在的な多峰性の可能性を考慮し、著者らはPigeons.jlを介したパラレル・テンパリング・マルコフ連鎖モンテカルロ（PTMCMC）サンプリングを使用している。この手法により、標準的なMCMC法（メトロポリス・ヘイスティングス法など）が失敗しやすい、低確率領域のピーク間をサンプラーが横断することを可能にしている。
4. 階層的スタッキング： 持続的な光源構造を抽出し、時間的変動を定量化するために、著者らはグローバルなベイズ階層モデルを構築した。全161個のスナップショットの事後分布をスタックして、グローバルな事後分布を形成する。このモデルは、平均的な光源構造（平均 $\mu_p$）と変動分布（標準偏差 $\sigma_p$）を推定し、平均的な構造の不確実性と、時間的変動およびモデルの誤設定による影響を分離する。

主な貢献

• 新規な推論パイプライン： 疎なEHT可視度データに対して、物理ベースの半解析的モデルをスナップショットごとに直接適合させるパイプラインを導入した。これにより、推論ループ中に完全なGRMHDシミュレーションを実行する必要がなくなる。
• 階層的変動分離： スナップショットの事後分布をスタックすることで、Sgr A*特有の急速な時間的変動と、時間平均された持続的な光源構造を分離することに成功した。
• 合成データによる検証： 本パイプラインは、既知の真のパラメータを持つGRMHDシミュレーションから生成された合成EHT観測データを用いて検証されている。このテストにより、推論パイプラインがスピン、傾斜角、または位置角の測定に重大なバイアスを導入しないことが確認された（ただし、スピンは依然として制約されていない）。

結果
この手法を2017年のSgr A*データに適用した結果、以下の制約が得られた（平均的な構造の不確実性 $\pm$ と、時間的変動／モデルの誤設定による不確実性 $\pm_v$ を分離）：

• 観測者傾斜角： 解析は、ブラックホールのスピン軸に対してほぼ正面を向いた観測者傾斜角 $\theta_o = 9.2^\circ \pm 3.6^\circ \pm_v 11.6^\circ$ を強く支持している。この結果はスナップショット間で一貫しており、低傾斜角を示唆する先行研究とも一致している。
• 放射半径： 放射のピークは地平線付近、$R_{peak} = 4.9 \pm 0.1 \pm_v 0.5 \, GM/c^2$ にあると制約される。
• スピン位置角： 投影されたスピン位置角は、垂直方向から反時計回りに $p.a. = 7.3^\circ \pm 7.08^\circ \pm_v 43.5^\circ$ と、ほぼ垂直に近いことが制約された。
• 放射幾何学： 主要な放射は赤道面から $43.4^\circ \pm 2.0^\circ \pm_v 5.9^\circ$ 上方に存在しており、顕著な非赤道放射（厚い円盤またはジェットの可能性）を示唆している。
• スピンと磁場： ブラックホールのスピン（$a_*$）および磁場パラメータ（$\eta, \iota$）は、ほとんど制約されていない。スピンの事後分布は事前分布の全範囲にわたってほぼ一様であり、これは2017年のEHT観測単独では、追加の仮定やデータなしにスピンを制約できないことを示している。
• プラズマ速度： 解析は、最内安定円軌道（ISCO）に対して比較的低速なプラズマ流 $\beta_v = 0.33 \pm 0.22 \pm_v 2.1$ を支持している。

意義と限界
本論文は、ブラックホールのスピン軸に対する観測者の傾斜角の独立した測定を提供しており、Sgr A*をほぼ正面から見ているという仮説を支持している。結果は、降着流が地平線付近および赤道面の上方の放射に支配されていることを示唆している。

しかし、著者らはブラックホールのスピンおよび磁場に関する制約について、データの疎性、パラメータ間の縮退、および非軸対称で時間変動する乱流を捉える上での軸対称デュアルコーンモデルの限界に起因するものとして、控えめな表現を用いている。著者らは、Sgr A*の急速な変動性と、複雑な流体力学を完全に捉えることができないモデルの性質が、重大な不確実性を招いていることを明示的に述べている。彼らは、本手法が合成テストにおいてGRMHDシミュレーションのパラメータを再現できる一方で、実際のデータ解析はモデルの誤設定の影響を受ける可能性があると指摘している。バイアスを低減し、スピンおよび磁場構成に関する制約を改善するためには、光学的厚さの効果、多成分放射（例：ジェット）、および偏光を取り入れる必要があることが今後の課題として挙げられている。