The Impact of Galaxy Formation on Galaxy Biasing, and Implications for Primordial non-Gaussianity Constraints

技術要約：銀河形成が銀河バイアスに与える影響と、原始非ガウス性制約への示唆

問題提起
パラメータ $f_{\text{NL}}$ は、宇宙の初期エネルギー揺らぎにおける局所的な原始非ガウス性（PNG）を定量化するものである。非ゼロの検出が得られれば、単一場または多場インフレーションモデルを区別することができる。次世代の銀河サーベイは、統計誤差 $\sigma_{f_{\text{NL}}} \sim 1$ という精度で $f_{\text{NL}}$ を制約することを目指しているが、現在、系統誤差がこれらの制約を支配している。具体的には、$f_{\text{NL}}$ の測定は、銀河バイアスパラメータである $b_1$（質量密度への応答）および $b_\phi$（原始ポテンシャル摂動への応答）と縮退している。

現在の制約は、ダークマターハローから導出された「普遍性関係」（例：$b_\phi = 2\delta_c(b_1 - 1)$）を仮定する「普遍性関係」に依存することが多い。しかし、近年の研究（例：Barreira et al. 2020）は、この関係が銀河については成立しないことを示している。なぜなら、$b_\phi$ は銀河の選択や銀河形成の具体的な物理学に依存するためである。銀河形成のモデリングの複雑さは、$b_\phi$ に重大な不確実性を導入し、将来の $f_{\text{NL}}$ 制約を脅かす要因となっている。

手法
銀河形成モデリングに由来する $b_\phi$ の不確実性に対処するため、著者らは、セパレート・ユニバース（separate-universe）N体シミュレーション内に組み込まれた Santa Cruz Semi-Analytic Model (SC-SAM) を利用している。

• シミュレーション: 3つの AREPO N体シミュレーションが、ボックスサイズ $L = 205 h^{-1}$ cMpc、粒子数 $N = 1280^3$ で実行された。これらのシミュレーションは同一のガウス初期条件を共有しているが、セパレート・ユニバース手法を用いて、長波長摂動に対する銀河形成の応答を効率的に測定するために、原始パワースペクトルの振幅（$\sigma_8$）が異なっている。
• モデルのバリエーション: 著者らは56個のユニークな SC-SAM インスタンスを実行した。これらのバリエーションは、以下の3つのフィードバックパラメータを変更することを含む：
  1. $\epsilon_{\text{SN}}$ （恒星駆動による質量放出振幅、因子 $A_{\text{SN1}}$ により変化）。
  2. $\alpha_{\text{rh}}$ （質量放出関係の指数、因子 $A_{\text{SN2}}$ により変化）。
  3. $\kappa_{\text{radio}}$ （AGN ラジオモード・フィードバック振幅、因子 $A_{\text{AGN}}$ により変化）。
• 測定:
  • $b_\phi$ は、高 $\sigma_8$、標準、低 $\sigma_8$ のシミュレーション間の銀河数カウントの有限差分によって計算される。
  • $b_1$ は、銀河・物質クロスパワースペクトルから導出される。
  • 銀河は、赤方偏移 $z=0, 1, 2$ において、恒星質量 ($M_*$)、星形成率 (SFR)、および比星形成率 (sSFR) に基づいて選択された。
  • 各モデルには、観測された $z=0$ の恒星質量関数（Bernardi et al. 2018）との一致度に基づいた「リアリズム」指標（$\chi^2$）が割り当てられた。
  • 本研究では、銀河を同様の質量を持つが形成履歴が異なるハローへとシャッフルすることによる、銀河のアセンブリ・バイアスの影響も評価した。

主要な貢献

1. フィードバックの系統的探索: 単一の流体シミュレーションや固定された SAM を用いることが多い従来の研究とは異なり、本研究では恒星および AGN フィードバックパラメータを系統的に変化させ、その結果生じる $b_\phi$ と $b_1$ のランドスケープをマッピングしている。
2. セパレート・ユニバース SAM の適用: 著者らはセパレート・ユニバース手法を半解析的モデル（SAM）に適用しており、これにより、フル流体シミュレーションでは計算量的に困難な広範なパラメータ空間にわたって $b_\phi$ を効率的に測定することを可能にした。
3. 選択依存性の分析: 本論文は、異なる銀河選択基準（$M_*$, SFR, sSFR）が銀河形成の物理の変化にどのように応答するかを明示的に比較し、モデリングの不確実性に対してどの選択が最も堅牢であるかを特定している。
4. 他モデルとの検証: 結果は、IllustrisTNG 流体シミュレーション、Munich SAM、および Galacticus SAM と比較され、SC-SAM のパラメータ空間が他の銀河形成モデルの挙動を包含しているかどうかが検討されている。

結果

• 普遍性からの逸脱: SC-SAM モデルは、普遍性関係（$b_\phi = 2\delta_c(b_1 - 1)$）が銀河については正確ではないことを裏付けている。$b_\phi(b_1)$ の関係は、銀河の選択および特定のフィードバックパラメータに依存して大きく変化する。
• 恒星質量選択: $b_\phi$ は一般に恒星質量とともに増加する。しかし、$b_\phi(b_1)$ の関係は IllustrisTNG (Barreira et al. 2020) よりも急峻であり、特に低質量選択において顕著である。この関係は恒星フィードバックパラメータ（$A_{\text{SN1}}, A_{\text{SN2}}$）に敏感である。
• SFR 選択: SFR で選択された銀河の $b_\phi$ 値は混沌としており、しばしばゼロに近い。この関係はフィードバックパラメータに対して非常に敏感である。なぜなら、異なるパラメータが SFR–ハロー質量関係を複雑な形で変化させるためである。
• sSFR 選択（主要な知見）: 比星形成率 (sSFR) に基づく選択は、銀河形成モデリングの変動に対して極めて堅牢である。56 個すべてのモデルにおいて、$b_\phi$ 対 sSFR および $b_\phi(b_1)$ の関係は、タイトで一貫した経路を辿る。これは、SC-SAM の自己調節的な性質により、フィードバックパラメータの変動に関わらず、一貫した sSFR–ハロー質量関係が維持されるためである。
• アセンブリ・バイアス: 本研究では、銀河のアセンブリ・バイアスが $b_1$ の測定に影響を与え（恒星質量選択ではクラスタリングを減少させ、SFR/sSFR 選択では形成履歴をシャッフルした際にクラスタリングを増加させる）、しかし、これらのアセンブリ・バイアスの信号が測定された $b_\phi(b_1)$ 関係には最小限の影響しか与えないことを見出した。
• 他モデルとの比較: SC-SAM のバリエーションは、Munich SAM および Galacticus で測定された $b_\phi$ の挙動を包含している。IllustrisTNG の結果とは（おそらく SAM に欠落しているバリオン効果やマニーツリー構築の違いによる）不一致があるものの、SC-SAM のバリエーションは、他のモデルと整合する広範な $b_\phi$ の値を提供している。

意義と主張
本論文は、銀河形成のモデリングにおける不確実性が、将来の $f_{\text{NL}}$ 制約における決定的な系統誤差であることを主張している。$b_\phi$ は普遍的ではなく、選択とモデリングに強く依存することを実証することで、著者らは、固定された普遍性関係を仮定することは、偏った制約をもたらすリスクがあると論じている。

本研究の主な意義は、sSFR を $f_{\text{NL}}$ 研究のための堅牢な銀河選択として特定したことにある。sSFR 選択は、幅広い銀河形成物理学にわたって一貫した $b_\phi$ 値をもたらすため、将来のサーベイにおける系統誤差を軽減する経路を提供する。著者らは、どの単一のモデルも宇宙を完全に記述することはできないが、sSFR によって銀河を選択することは、原始非ガウス性の制約に対する銀河形成の詳細な物理への感度を最小限に抑えることができると結論付けている。さらに、本研究は、観測データに基づいた $b_\phi$ の情報的事前分布を構築するための枠組みを提供し、物理的に動機付けられた原始非ガウス性の制約へと向かう道筋を示している。