How I stop worrying about non-universality and $b_\phi$: Constraining local $f_{\rm NL}$ with $b_\phi$ priors from HOD posteriors

技術要約：HOD事後分布からの$b_\phi$事前分布を用いた局所的$f_{\text{NL}}$の制約

問題提起
局所的な原始非ガウス性（local PNG）は、振幅パラメータ$f_{\text{NL}}$によって特徴付けられ、銀河のクラスタリングにおけるスケール依存的なバイアスを、$b_\phi f_{\text{NL}}$の積に比例する形で誘起する。ここで、$b_\phi$は長波長の原始ポテンシャルの摂動に対する銀河形成の応答を表す。観測される銀河サーベイから$f_{\text{NL}}$に対してロバストな制約を得る上での主要な障害は、$b_\phi$の不確実性である。$b_\phi$に関する正確な知識がなければ、クラスタリングのデータは退化的な積である$b_\phi f_{\text{NL}}$しか制約できず、$f_{\text{NL}}$自体の中心値や不確実性を信頼性高く推論することが不可能となる。

これらを解決するための一般的なアプローチには、「普遍性関係（universality relation）」（ハロー質量関数を介して$b_\phi$を線形バイアス$b_1$に結びつけるもの）を仮定することや、宇宙論的シミュレーションを用いて$b_\phi$を較正することがある。しかし、普遍性関係は、色やフラックスで選択された観測済みの銀河に対しては近似として不十分であり、シミュレーションに基づく較正も、ターゲットとなるサンプル固有のクラスタリング統計ではなく、一般的な銀河形成の観測量に依存していることが多い。あるいは、トレーサーの存在数（$dn/dz$）の赤方偏移進化から$b_\phi$を直接推論するには、複雑なサーベイ選択効果のモデリングが必要となる。

手法
著者らは、HOD（Halo Occupation Distribution）モデルを介して、ターゲットとなる銀河サンプルの小スケール・クラスタリング統計と宇宙論的シミュレーションを組み合わせる「第三のルート」を提案している。コアとなる手法は以下のステップを含む：

1. HOD事後分布サンプリング： 著者らは、DESI Early Data Release (EDR) One-Percent Survey (SV3) のLuminous Red Galaxies (LRGs) に関するHOD制約（赤方偏移ビン：LRG1: $0.4<z<0.6$, LRG2: $0.6<z<0.8$, LRG3: $0.95<z<1.1$）を利用する。彼らはHODパラメータの事後分布をサンプリングし、DESIの全データセットにおける潜在的な変動を考慮して、不確実性を$3\sigma$（またはLRG3については$1\sigma$）に拡大させている。
2. セパレート・ユニバース・モック生成： サンプリングされたHODパラメータを用いて、AbacusSummit シミュレーション・スイートから一対のモック銀河カタログを生成する。一方は標準的な宇宙論（$c000$, $\sigma_8 = 0.807952$）、もう一方は低$\sigma_8$の宇宙論（$c004$, $\sigma_8 = 0.749999$）に対応しており、他のパラメータは固定されている。
3. $b_\phi$の測定： 各HOD実現値について、「セパレート・ユニバース（Separate Universe）」手法を用いて$b_\phi$を測定する。これは$\sigma_8$の差を長波長モードへの応答のプロキシとして扱い、$b_\phi \approx 2 \frac{\ln N - \ln N_{\text{fid}}}{\ln \sigma_8 - \ln \sigma_{8,\text{fid}}}$ を計算する。
4. 事前分布の構築： 100個のサンプリングされたHOD実現値にわたる測定された$b_\phi$の分布をガウス分布にフィッティングし、データ駆動型の事前分布 $b_\phi \sim \mathcal{N}(\mu, \sigma^2)$ を作成する。
5. 検証と推論： これらの事前分布を、明示的な局所的PNG信号（$f_{\text{NL}} = 0, \pm 30$）を含む AbacusPNG シミュレーションに適用する。著者らは、赤方偏移空間パワースペクトルのモノポールおよびクアドラポールのベイズ推論を行い、以下の3つのモードを比較する：
   • $b_{\text{fnl}}$ モード（事前分布なし）： 平坦な事前分布を用いて、積 $B_i = b_{\phi,i} f_{\text{NL}}$ を制約する。
   • $b_{\text{fnl}}$ モード（事前分布あり）： 共通の$f_{\text{NL}}$を周辺化しつつ、ガウス型の$b_\phi$事前分布を課す。
   • $b$-$p$ モード（普遍性）： 普遍性関係 $b_\phi = 2\delta_c(b_1 - 1)$ を仮定する。
6. アセンブリ・バイアス・テスト： ハローのアセンブリ・バイアス（固定質量におけるハロー環境への銀河占有の依存性）に対して、デコレーテッドHODモデルを用いて手法の堅牢性をテストする。

主な結果

• $f_{\text{NL}}$の回収： 已知の$f_{\text{NL}}$を持つAbacusPNGモックに適用した際、本手法は真の入力値を正常に回収した。HOD由来の$b_\phi$事前分布を含めることで、$f_{\text{NL}}$に対する$68\%$信頼区間は、事前分布なしの場合と比較して**9%から55%**タイトになった。最大の改善は、線形バイアスが高い高赤方偏移（LRG3）で見られた。
• 普遍性との比較： HOD事前分布を用いた制約は、厳密な普遍性関係（$p=1$）を仮定した場合よりもわずかに広くなった。しかし、著者らは、観測された銀河は普遍性関係に従わないことが理論的に予想されるため、HOD事前分布はより柔軟で物理的に動機付けられた代替案であり、検証されていない仮定を回避できると述べている。
• 結合制約： 3つのLRG赤方偏移ビン（LRG1–LRG3）のデータを組み合わせることで、$f_{\text{NL}}$の制約はLRG3のみを使用した場合と比較して約**65%**改善した。これは、有効体積の増加と宇宙論的分散の減少によるものである。
• アセンブリ・バイアスへの堅牢性： データベクトルにアセンブリ・バイアス効果が含まれている場合でも、推論された$f_{\text{NL}}$は安定している。解析の結果、アセンブリ・バイアスは主にガウス線形バイアスパラメータ$b_1$を（$2\text{--}3\sigma$）シフトさせるものの、スケール依存的なPNG項（$b_\phi f_{\text{NL}}$）には影響を与えないことが示されており、この手法がこの特定の系統誤差に対して堅牢であることを示している。

意義と主張
本論文は、シミュレーションの較正と直接的な観測的推論の間のギャップを埋める、局所的PNGを制約するための「第三のルート」を提供すると主張している。小スケール・クラスタリング（HOD経由）を大スケール・バイアス（$b_\phi$）へとマッピングすることで、普遍性関係や複雑な$dn/dz$モデリングに依存することなく、サンプル固有の事前分布を生成する。

著者らは、このアプローチが普遍性の仮定に対する柔軟な代替案を提供することを強調している。普遍性関係はよりタイトな制約をもたらすが、その制約は観測された銀河が従わないことが予想される関係に基づいている。HOD事前分布は、モデルの仮定に対する堅牢性と引き換えに、わずかな精度を妥協するものである。

本研究は、以下の制限を持つパイロット研究として提示されている：

• 現在はLRGに限定されており、HODの理解が進んでいないELG（Emission Line Galaxies）やクエーサーへの拡張は容易ではない。
• 制約は、仮説的な有効体積を持つ立方体モックから導出されており、実際のDESIデータではない（ただし、DESI DR1/DR2への適合が自然な次のステップとして特定されている）。
• 事前分布の構築において、HODパラメータを独立したガウス分布として扱っており、共分散を無視している。これは著者らも近似の源泉として認めている。

最終的に、本論文は、小スケール・クラスタリングの測定や銀河形成の理解が進むにつれ、これらの$b_\phi$事前分布が鋭くなり、Stage-Vサーベイ（MUST、Spec-S5、WSTなど）が局所的PNGを通じてインフレーション物理学を探索する展望を高めるであろうと結論付けている。