eROSITA cosmology with galaxy groups: hot gas budget out to the virial radius

技術要約：eROSITAによる銀河群を用いた宇宙論：ビリアル半径に至る高温ガスの総量

問題提起
非重力的なプロセス、具体的には活動銀河核（AGN）フィードバックは、高温ガスを銀河群や銀河団のビリアル領域の外側へと排出させることが知られており、メガパーセクスケールにおける物質分布を再形成する。銀河群の$R_{500}$（500倍の臨界密度を含む半径）内における高温ガス分率（$f_{gas}$）が、宇宙の平均値よりも低いことは確立されているが、$R_{200}$といったより大きな半径における測定は、高温ガス（ICM）の表面輝度の低さと背景放射レベルの上昇により、依然として乏しく不確実である。この不確実性は、フィードバックの兆候（加熱および排出されたガス）が最も顕著に現れるはずの領域である外縁部において極めて重要である。さらに、観測的な制約と、流体解析シミュレーション（FLAMINGO、BAHAMAS、SIMBAなど）との間には、バリオンの再分配におけるフィードバックモデルの有効性に関する相違が存在する。これらの相違は、次世代の大規模構造サーベイにおける主要な系統誤差となる、物質パワースペクトル（MPS）の抑制に関する予測に直接的な影響を与える。

手法
著者らは、eROSITA-DEの第1回公開データ（eRASS1）から選定され、2MASS赤外線銀河サーベイ（2MRS）の光学銀河カタログとクロスマッチされた、25個の銀河群からなる完全なサンプルを利用している。サンプルの空間分解能を確保し、未分解の銀河放射による汚染を最小限に抑えるため、低赤方偏移（$0.01 < z < 0.025$）に限定している。

主な手法ステップは以下の通りである：

1. データ処理： 著者らは、インストゥルメンタルな輝線や「ライトリーク」のアーティファクトを避けるため、[0.6–2.3] keVエネルギーバンドにフィルタリングした上で、eSASSソフトウェアシステムを用いてeRASS1イベントリストを処理する。光源検出と点源のマスキングには、マルチスケール・ウェーブレット分解を採用している。
2. プロファイル抽出： 背景放射を差し引いた画像から、$R_{200}$までの表面輝度（$SB_x$）プロファイルを抽出する。外縁部での信号対雑音比（S/N比）を確保するため、25個の銀河群のプロファイルをスタッキング（重ね合わせ）してフィッティングを行う。
3. 質量導出：
   • 高温ガス質量（$M_{gas}$）： $L_X-T$関係に基づく一定温度の仮定を用いて、$SB_x$プロファイルを放射量（EM）プロファイルに変換することで導出する。マルチスケール・フォワードフィッティング法を用い、球対称を仮定して2次元のEMプロファイルを3次元の密度プロファイルへとデプロジェクション（投影逆変換）する。
   • 全質量（$M_{tot}$）： 速度分散－質量関係に対して検証された、2MRSの力学的質量推定値を採用する。
   • ガス分率（$f_{gas}$）： 特定の開口部（$R_{500}$および$R_{200}$）内における$M_{gas}$と$M_{tot}$の比として計算される。
4. 統計分析： 著者らは、スケーリング関係（$M_{gas}-M_{tot}$、$L_X-M_{tot}$、$L_X-M_{gas}$）をフィッティングするためにベイズ階層分析を行う。極めて重要な点として、積分された物理量が積分半径に依存することによって生じる人工的な関係の平坦化に対処するため、**開口部共分散（aperture covariance）**の補正を実装している。
5. 宇宙論的推論： 測定されたバリオン分率とSP(k)モデル（ANTILLESシミュレーション・スイートに基づく）を用い、ダークマターのみの宇宙に対する物質パワースペクトルの抑制を推定する。

主な結果

• 表面輝度プロファイル： スタッキングされたプロファイルは、$R_{500}$内では一様に平坦な$\beta$パラメータ（$\beta = 0.38 \pm 0.04$）を示すが、$R_{500}$を超えて$R_{200}$までは、より急峻な$\beta = 0.76 \pm 0.19$へと変化する。この急峻化は、バリオン放出の弱まりを示唆しており、外縁部のプロファイルをより平坦にすると予測する強いフィードバックモデルとは対照的である。
• 高温ガス分率：
  • 中央質量（$M_{500} = 2.54 \times 10^{13} M_\odot$）において、$f_{gas,500} = 4.32 \pm 0.42\%$である。
  • $R_{200}$（中央質量 $M_{200} = 3.69 \times 10^{13} M_\odot$）において、この分率は$f_{gas,200} = 5.78 \pm 0.69\%$へと増加する。
  • 両方の値とも、宇宙のバリオン分率（約15.7%）を大幅に下回っており、大規模なバリオンの変位を裏付けている。
• スケーリング関係： 開口部共分散を補正した関係は以下の通りである：
  • $M_{gas}-M_{tot}$：$R_{500}$における傾きは$1.59^{+0.34}_{-0.28}$であり、自己相似的な期待値よりもわずかに急である。
  • $L_X-M_{tot}$：観測された光度－質量関係は、標準的なFLAMINGO（L1_m9）およびBAHAMASシミュレーションによって最もよく再現される。強いフィードバックを持つFLAMINGOのバリアント（$f_{gas}-4\sigma$, $f_{gas}-8\sigma$）は、光度を$6.5\sigma$から$8.6\sigma$過小評価している。
• シミュレーションとの比較：
  • 標準モデル： 標準的なFLAMINGOおよびBAHAMASシミュレーションは、$f_{gas}-M_{tot}$および$L_X-M_{tot}$関係に最も近い一致を示す。
  • 強いフィードバックモデル： FLAMINGOの強いフィードバック・バリアントおよびSIMBAシミュレーションは、$L_X-M_{tot}$関係および$R_{500}$以降の密度プロファイルの急峻化に関して、データと顕著な不一致（テンション）を示す。
  • その他のシミュレーション： IllustrisTNG100、EAGLE、およびFABLEは、中央質量におけるガス分率を過大評価する傾向がある。
• 物質パワースペクトル： SP(k)モデルを用い、ダークマターのみの宇宙に対する$k = 5 h \text{ Mpc}^{-1}$での物質パワースペクトルの減少を、$10\%–15\%$（$0.85 \le P(k)/P_{DMO}(k) \le 0.90$）と推定した。この結果は標準的なFLAMINGOおよびBAHAMASと一致するが、小さなスケールにおいては強いフィードバック・バリアントとの間で不一致が大きくなることを明らかにしている。

意義および主張
本論文は、光学的に選別された完全な銀河群サンプルについて、$R_{200}$における$M_{gas}-M_{tot}$関係に関する、最初の個別のeROSITAによる制約を提示したと主張している。最近のスタッキング研究の2倍の統計的精度を達成することで、本研究は、宇宙論的流体解析シミュレーションにおけるAGNフィードバックの実装に対する強力なテストを提供する。

著者らは、以下の点を強調している：

1. フィードバックモデルの制約： $R_{500}$以降の密度プロファイルの急峻化と特定のガス分率は、コアの密度を抑制する一方で外縁部を標準的な値に留めるような、強いフィードバックモデルに課題を突きつけている。
2. シミュレーションの検証： 標準的なFLAMINGOおよびBAHAMASとの一致は、これらの特定のフィードバック較正を支持するものである。一方で、強いフィードバック・バリアントとの不一致は、これらのモデルがガスを排出しすぎているか、あるいは高温ガスの熱力学的状態を正確に捉えられていない可能性を示唆している。
3. 宇宙論への影響： 導出されたバリオン分率は、「バリオニフィケーション（baryonification）」補正スキームへの不可欠な入力となり、Euclidやナンシー・グレース・ローマン宇宙望遠鏡のような次世代のサーベイにおける物質パワースペクトルの抑制予測に直接的な情報を提供する。結果は、フィードバックがパワーを抑制する一方で、その抑制の大きさはシミュレーションに実装された具体的なフィードバック物理学によって厳密に制約されることを示唆している。