Learning the Universe: Constrained simulations of the Coma galaxy cluster -- I. Radial X-ray and Compton-y signatures

技術要約：宇宙を学習する：Coma銀河団の制約付きシミュレーション — I. 放射状X線およびコンプトン-$y$シグネチャ

問題提起
大規模な体積を持つ宇宙論的シミュレーション（例：IllustrisTNG、MillenniumTNG）は、銀河集団の統計的性質や銀河団のスケーリング関係を再現することに成功しているが、これらは宇宙の限られた体積しかサンプリングしていない。その結果、Coma銀河団のような稀で巨大な系を極めて少数しか含んでおらず、特定の局所銀河団の高精度な観測結果と詳細な、オブジェクトごとの比較を行うには不向きである。一方で、従来の局所銀河団を対象とした「ズームイン」シミュレーションの多くは、単一の実数（realisation）に依存しているか、あるいはIllustrisTNGのような最先端のモデルとは異なる銀河形成物理学を採用していた。したがって、局所の大規模構造（環境）に一致させるための制約付き初期条件と、現代的な銀河形成物理学を用いた高解像度の流体計算を組み合わせ、特定の銀河団の多波長観測を解釈するための、統計的に堅牢な枠組みが必要とされている。

手法
著者らは、Coma銀河団の類似体に対する50個の高精度なズームイン・シミュレーション・スイートを提示している。その手法は、以下の3つの異なる要素を統合している。

1. 制約付き初期条件 (ICs): シミュレーションには、2M++銀河カタログを用いて局所宇宙の広域構造を再構成した、BORG/Manticoreによるベイズ推論（McAlpine et al. 2025）からの事後分布サンプルを利用している。これにより、Comaの類似体の大規模な環境、フィラメント構造、および質量が観測的制約と一致することが保証される。
2. アンサンブル・アプローチ: 単一の「最良」の実数を選択するのではなく、著者らは50個の独立した事後サンプルを進化させている。このアンサンブルは、宇宙論的分散（cosmic variance）と局所的な初期条件の不確実性の両方を捉えており、Coma質量の銀河団が取り得る形成史と力学的状態の統計的分布を提供する。
3. 流体進化: 50個の実数のそれぞれは、AreposコードとIllustrisTNG (TNG) 銀河形成モデルを用いて進化させられている。シミュレーションにはマルチレベルの精細化スキームが採用されており、高解像度領域におけるダークマター粒子の質量は $6.0 \times 10^7 M_\odot$、ガスセルの質量は $1.2 \times 10^7 M_\odot$ に達する。物理プロセスには、磁気流体力学 (MHD)、金属線冷却、恒星フィードバック、および二モードAGNフィードバックが含まれる。

シミュレーションはComaの観測されたスカイフレームに合わせて配置されており、観測量の直接的な投影を可能にしている。著者らは、熱的サンヤエフ・ゼルドビッチ (tSZ) 効果（コンプトン-$y$）とX線表面輝度（0.4–2 keV）について、投影マップおよび方位角平均の放射状プロファイルを算出し、eROSITAおよびPlanck衛星のデータと直接比較している。

主な貢献

• 初のTNGベースの制約付きアンサンブル: 本研究は、BORG/Manticoreの制約付き初期条件から導出された、完全なIllustrisTNG物理スイートを用いて進化した、初のComa類似体のアンサンブルを提供する。
• 局所銀河団のための統計的枠組み: 50個の実数を実行することで、大規模な環境が固定されている場合でも、形成史や合体活動の変化から生じる銀河団の特性（質量、ガス分率、熱力学プロファイル）の固有の散らつきを定量化する。
• 直接的な多波長比較: 論文は、コンプトン-$y$およびX線表面輝度のシミュレーションによる放射状プロファイルを、特定の観測データセット（PlanckおよびeROSITA）と比較するための厳密なパイプラインを確立し、「最良の一致」を示す実数のサブセットを特定している。

結果

• スケーリング関係: 50個のComa類似体のアンサンブルは、TNG-Clusterスイートおよび観測のコンピレーションと一致する標準的な銀河団のスケーリング関係（ガス分率、$Y_{SZ}$–質量、$L_X$–質量、恒星質量–ハロー質量、およびブラックホール・スケーリング関係）を再現している。シェル積分された $Y_{SZ,500}$ およびガス分率 ($f_{gas}$) において、Planckの正規化に対して系統的かつ緩やかな過剰が認められるが、これはTNGのクラスタースケールにおける既知の特性と一致している。
• 放射状プロファイルの散らつき: 50個の実数は、コンプトン-$y$およびX線表面輝度の内側放射状プロファイルにおいて顕著な散らつきを示している。アンサンブルの平均値は観測されたComaのプロファイルを挟み込む形となるが、個々の実数は大きく変動する。
• 最良の一致を示す実数: 10個の実数が、観測された放射状プロファイルに対する最良の一致として特定された。これらの最良一致銀河団は、広い範囲の全ハロー質量 ($M_{200c} \approx 0.87\text{--}1.45 \times 10^{15} M_\odot$) に及び、多様な二次元形態（緩和したものから乱れたものまで）を示しており、1Dの放射状プロファイルへの良好な一致が、2D構造や全質量を一意に決定するものではないことを証明している。
• 中心部の散らつきの起源: 中心部の輝度の散らつきは、主に銀河団内ガス (ICM) の中心エントロピーによって駆動されている。
  • 低エントロピー（クールコア）銀河団（50個中16個）は、中心部のX線およびコンプトン-$y$輝度が高い。
  • 高エントロピー銀河団（50個中34個）は、中心部が暗い。
  • この中心エントロピーの状態は、短時間のAGNバーストや最近の主要な合体によって引き起こされる一時的な特徴ではない（サンプル内の $z < 0.3$ においてこれらは稀である）。むしろ、それは長期的な形成史によって決定されている。
• 形成史との相関: 低エントロピー銀河団は、「早期形成型」であり（中央値形成赤方偏移 $z_f(10\%) \approx 1.59$）、高エントロピー銀河団（$z_f(10\%) \approx 1.15$）と比較して形成時期が早い。これらは主要な合体よりも、「ソフト」な経路（小規模な合体や滑らかな降着）を通じて主に質量を蓄積してきた。この「ソフトに成長した」歴史が、中心ICMが低エントロピー状態を発達させ、維持することを可能にしている。中心ブラックホールの成長もこの歴史を追跡しており、低エントロピーのBCGは $z=1$ においてより大規模なブラックホールを持ち、その後の成長は緩やかである。
• 時間変動性: 大半の銀河団において、放射状プロファイルは過去500 Myrの間、安定している。注目すべき例外は、最も適合度の高い実数の一つであるComa_042であり、これは後期段階のマイナー・マージャー活動によってプロファイルが急速に再形成されており、最良の一致であってもそれが一時的なものである可能性を浮き彫りにしている。

意義と主張
本論文は、制約付き宇宙論的初期条件と最先端の銀河形成物理学を組み合わせることが、特定の銀河団の標的となる観測駆動型の類似体を生成するための効果的な戦略であることを主張している。得られたデータセットは、以下のための貴重なリソースとなる：

1. ICM物理のテスト: フィードバック、合体、および大規模環境が、観測可能なクラスター特性をどのように形作るかを研究するための制御された環境を提供する。
2. スケーリング関係の較正: 銀河団の観測量における散らつきを解釈するための統計的基準を提供する。
3. 観測の解釈: 局所宇宙の事後分布に基づくフォワードモデル・エミュレータを提供することにより、今後のX線およびサンヤエフ・ゼルドビッチ観測（例：XRISMによるもの）の解釈を可能にする。

著者らは、観察されたアンサンブルの散らつきは、較正の不確実性ではなく、局所宇宙のデータに残された自由度を反映した物理的な事後分布の散らつきであることを強調している。彼らは、将来の研究において、制約付き推論の尤度をクラスター規模の熱力学的観測量（コンプトン-$y$およびX線プロファイル）に拡張することで、Coma銀河団の形成史に対する制約をさらに強めることができる可能性を示唆している。