Testing the inference of kinematics from mock JWST NIRSpec/MSA observations of TNG50 galaxies at $z\sim2-6$

技術要約：TNG50銀河の模擬JWST NIRSpec/MSA観測からの運動学の推論に関する検証

問題提起
銀河の動力学的進化を理解するには、宇宙時間における電離ガスの回転速度（$v$）と速度分散（$\sigma$）を測定する必要がある。IllustrisTNGのような流体シミュレーションは、6次元の空間および運動学的データを提供しているが、特定の観測セットアップに対する明示的な放射伝達解を欠いている。一方、JWSTのNIRSpecマイクロシャッターアレイ（MSA）を用いた観測研究は、2次元スペクトルから運動学を推論するためにフォワードモデリングに依存している。これらのモデルは通常、銀河が軸対称で薄い回転円盤であり、セルシック強度プロファイルとアークタンジェント速度曲線を持つことを仮定している。しかし、高赤方偏移（$z > 2$）の銀河は、しばしば複雑な形態（球状、塊状）や非円運動（流入/流出）を示し、これらの仮定に抵触する場合がある。単純なパラメトリックモデルが、典型的な高赤方偏移銀河の固有の運動学を正確に回収できるかどうかについての統計的な検証、特にIFU（積分視野分光器）と比較した際のMSAの限られた空間分解能やスリットベースの性質を考慮した検証が決定的に不足している。

手法
著者らは、高解像度のTNG50宇宙論的シミュレーションを利用して、赤方偏移 $z = 2, 3, 4, 6$ における4,000以上の星形成銀河（$M_* = 10^8 - 10^{11.5} M_\odot$）に対するエンドツーエンドの模擬観測パイプラインを生成した。

1. 固有の測定値： 著者らは、H$\alpha$放射ガスの有効半径（$r_e$）内における「固有」の運動学的量（$v_{max}$ および $\sigma_e$）を定義した。彼らは、アパーチャサイズ、空間ビンニング（0.5 vs 5 ckpc）、および重み付けスキーム（重みなし vs 星形成率重み付け）を含む測定の選択に対する、固有値の感度を明示的にテストした。
2. 模擬観測：
   • イメージング： 星形成ガスセルの投影、Kennicutt (1998) の関係式によるSFRからH$\alpha$光度への変換、およびNIRCam PSFによる畳み込みにより、模擬NIRCam画像が生成された。ノイズは、深いJADES級の露出を模倣するように加えられた。
   • 分光： 模擬NIRSpec/MSA 2Dスペクトルは、ランダムな軸に沿ってガスセルを投影し、視線方向の速度を割り当て、NIRSpec LSFおよびPSFで畳み込みを行うことで作成された。最も深い既存のJWSTデータと同等のS/N比（S/N $\sim$ 50）を達成するために、ガウスノイズが適用された。
3. 推論パイプライン：
   • 形態： 模擬NIRCam画像は pysersic でフィッティングされ、形態学的事前分布（有効半径 $r_e$、セルシック指数 $n_s$、軸比 $q$、位置角 PA）が導出された。
   • 運動学： 模擬NIRSpecスペクトルは、アークタンジェント速度プロファイルと一定の速度分散を持つ薄い回転円盤モデルを仮定して msafit を用いてフィッティングされた。導出された形態パラメータは強い事前分布として使用される。フィッティングには、塊状の光分布を考慮するための非パラメトリックな行ごとのフラックススケーリングを含む、重み付き $\chi^2$ 最小化が採用されている。

主な貢献と結果
本研究では、模擬スペクトルから導出された最良適合運動学パラメータを、固有のシミュレーション値と比較することで、観測バイアスを定量化している。

• 回収精度： $z=3$ における整列かつ分解されたソースについては、推論された回転速度（$v$）と速度分散（$\sigma$）は、それぞれ固有の測定値に対して概ね因子 $\sim 2$ および $\sim 1.5$ 以内で一致している。$v$ の散らつきは $\sigma$ よりも大きい。
• 系統的バイアス：
  • 形態への依存性： 滑らかで円盤状の系については回収は堅牢である。しかし、本質的に細長い、あるいは球状の系（大きな固有軸比 $C/A$）については、$v$ が過大評価され、$\sigma$ が過小評価される傾向があり、その差は因子2を超える。これは、薄い円盤の仮定が、非回転的な広がりを傾斜角の効果として解釈するようにモデルを強制するためである。
  • 赤方偏移の進化： 高赤方偏移（$z > 3$）では、銀河が薄い円盤の仮定からより強く逸脱するため、回収精度の散らつきが増大する。
  • 個体 vs 集団： 個々の銀河の運動学（$v/\sigma$）は、誤分類される可能性がある（例：分散支持系が回転支持系として現れるなど、因子2〜3の誤差）。しかし、集団の傾向は堅牢に回収される。$\sigma$ および $v/\sigma$ の弱い赤方偏移進化、および全球的なSFRに対する $\sigma$ の依存性は、模擬データにおいて保持されている。
• スケーリング関係： 推論された $\sigma$-SFR関係は、「質量輸送＋フィードバック」の理論モデル（Krumholz et al. 2018）とよく一致しており、個々の不確実性にもかかわらず、運動学的モデリングが銀河集団の物理的駆動要因（乱流の要因としてのフィードバック vs 重力）を信頼性高く制約できることを示唆している。
• 失敗モード： 約5〜10%のソースが薄い円盤モデルのフィッティングに失敗しており、これはしばしば、不規則な形態、急峻な中心速度勾配を持つ高ステラ質量系（アークタンジェントプロファイルでは捉えきれない）、または固有の $\sigma$ がモデルの下限（10 km s$^{-1}$）を下回るソースに起因している。

意義と主張
本論文は、NIRSpec MSA観測を宇宙論的シミュレーションと比較して、電離ガスの運動学の測定におけるバイアスを定量化するための、最初のエンドツーエンドのフレームワークを提供すると主張している。著者らは、単純な薄い円盤モデルは、高赤方偏移銀河の個々の「円盤性」（$v/\sigma$）を決定するには不十分であるが（大きな不確実性を伴うため）、統計的な集団の傾向を回収するには効果的であると結論付けている。

$z \ge 2$ における個々の系の $v/\sigma$ に関する観測的推論は、薄い円盤の幾何学的な仮定や非円運動の存在によってバイアスを受ける可能性があるため、注意深く扱うべきであると本研究は強調している。著者らは、初期銀河の動力学的状態の制約を改善するためには、将来の研究が2次元スリット分光を超えて、3D IFU観測やスリットステッピング戦略へと移行し、高赤方偏移銀河に固有の複雑な3D形態や塊状性を考慮した専用の運動学モデルを開発する必要があると主張している。この研究は、現在の限界が単なる計器的な問題ではなく、単純なパラメトリックモデルと銀河形成物理学の本質的な複雑さとの間のミスマッチに起因していることを強調している。