Exploring non-Poisson satellite occupation in HOD models and its impact on 2- and 3-point galaxy clustering

技術的サマリー：HOD モデルにおける非ポアソン型衛星銀河の占拠の探求と、2 点および 3 点銀河クラスターリングへの影響

問題提起
銀河とダークマターハロの間の関係は現代宇宙論の基盤であり、通常、ハロ占拠分布（HOD）フレームワークを用いてモデル化される。HOD モデルにおける標準的な仮定は、固定質量のハロ内の衛星銀河の数がポアソン分布に従うというものである。しかし、数値シミュレーションおよび観測解析は、ハロの組立バイアス、潮汐ストリッピング、相関した落下履歴などの効果により、衛星集団が厳密なポアソン統計からの逸脱を示す可能性を示唆している。これらの逸脱は、衛星数のサブポアソン型またはスーパーポアソン型の散乱として現れ、占拠数の第 2 階乗モーメントを変化させる。このモーメントは直接、クラスターリング統計における 1 ハロ項の振幅を制御するため、非ポアソン挙動を無視することは、DESI や Euclid などのステージ IV 調査から導き出される小規模クラスターリング解析および宇宙論的制約に系統的誤差を導入する可能性がある。

手法
これに対処するため、著者らは標準的なポアソンモデルの最小かつ柔軟な拡張として、コンウェイ・マクスウェル・ポアソン（CMP）分布を導入する。CMP 分布は、確率質量関数 $P(N_{sat}=k | \lambda, \nu) = \frac{\lambda^k}{(k!)^\nu Z(\lambda, \nu)}$ によって定義され、ここで $\nu$ は分散パラメータである：$\nu=1$ はポアソン分布を回復し、$\nu < 1$ はスーパーポアソン型挙動（分散 > 平均）をもたらし、$\nu > 1$ はサブポアソン型挙動（分散 < 平均）をもたらす。

手法は 3 つの主要な技術的構成要素を含む：

1. 解析的近似と数値スキーム： 正規化定数 $Z(\lambda, \nu)$ に閉じた形がないため、著者らは 2 つの領域における平均と分散の解析的近似を導出した：
   • 小 $\lambda$（$\lambda \to 0$）： 平均と分散を近似するために、$O(\lambda^4)$ までの級数展開が導出される。
   • 大 $\lambda$（$\lambda \to \infty$）： 鞍点法に基づく漸近展開が利用される。
   • グローバルスケーリング： これらの領域を接続する滑らかな補間が誤差関数を用いて構築され、ターゲット平均占拠 $\langle N_{sat}(M) \rangle$ と CMP パラメータ $\lambda$ の間の連続的なマッピングが確保され、$0.5 < \nu < 2$ に対して約 5% の精度が達成される。
2. モックカタログの生成： HODDIES パッケージと AbacusSummit N 体シミュレーション（Planck 2018 $\Lambda$CDM 宇宙論）を用いて、著者らは輝赤色銀河（LRG）サンプルを表す 16 のモック銀河カタログを生成する。標準的な HOD パラメータは固定され、CMP 分散パラメータ $\nu$ は 0.5 から 2.0 まで変化させられる。
3. クラスターリング統計の解析： $\nu$ を変化させる影響は以下の指標を用いて定量化される：
   • 2 点統計： 投影相関関数（$w_p$）、および赤方偏移空間相関関数のモノポール（$\xi_0$）と四重極（$\xi_2$）。
   • 3 点統計： $k$ 範囲 $[0, 0.3] \, h/\text{Mpc}$ において Sugiyama 推定量を用いて計算された銀河ビスペクトル（モノポール $B_{000}$ と四重極 $B_{202}$）。
   • 高次モーメント： 局所的な銀河環境と分布の高次モーメントを探るための円柱内カウント（CiC）統計。

主要な貢献と結果

• HOD における CMP の実装： 本論文は、HOD フレームワーク内で CMP 分布を実装するための堅牢な数値スキームを提供し、物理的な平均 $\langle N_{sat} \rangle$ を維持しつつ、パラメータ $\nu$ を介して平均衛星占拠と分散を成功裡に分離する。
• 占拠統計への影響： $\nu$ を変化させることは、衛星数の分布を著しく変える。$\nu < 1$ の場合、分布は広がり（スーパーポアソン型）、高占拠ハロの確率が増加する。$\nu > 1$ の場合、分布は狭まる（サブポアソン型）。衛星の平均数はポアソンの場合に対して最大で $\pm 5\%$ 変動し、この逸脱は低 $\lambda$ における CMP 分布の非対称性に起因すると考えられる。
• 2 点クラスターリングへの影響： 非ポアソン統計は主に 1 ハロ項（小スケール、$r_p \lesssim 1 \, h^{-1}\text{Mpc}$）に影響を与える。
  • 投影相関関数 $w_p$ は、ポアソンの場合に対して最大 10% の逸脱を示す。
  • モノポールと四重極は最大 5% の逸脱を示す。
  • これらの変化は分散と相関しており、より大きな分散（より低い $\nu$）はクラスターリングパワーを増強し、より小さな分散はそれを抑制する。
• 3 点クラスターリング（ビスペクトル）への影響： 樹木レベルの銀河ビスペクトル（Sugiyama 基底において）は、$k_{\text{max}} = 0.3$ の範囲で非ポアソン仮定に対して頑健であることが判明した。
  • ビスペクトルモノポールの変動は 2% 未満である。
  • 著者らは、現在および近未来の大規模解析（$k < 0.3$）において、非ポアソン型衛星統計がビスペクトルから導き出される宇宙論的制約にバイアスをかける可能性は低いと結論づける。なぜなら、その信号は予想されるデータ誤差よりも小さいからである。
• 円柱内カウント（CiC）の感度： CiC 統計は衛星占拠の分散に対して高い感度を示す。
  • ポアソンの場合からの逸脱は最大 30% に達する。
  • これは、CiC が HOD パラメータの縮退を破り、非ポアソン挙動を検出するための強力なプローブであることを確認し、この点において標準的な 2 点統計を上回ることを示している。

意義と主張
本論文は、CMP 分布が HOD フレームワーク内でポアソン仮説を検証するための「最小かつ制御されたメカニズム」を提供すると主張する。著者らは、非ポアソン統計が小規模クラスターリングの正確なモデル化（1 ハロ項に著しく影響する）にとって重要である一方、大規模（$k < 0.3$）におけるパワースペクトルおよびビスペクトル測定から導き出される宇宙論的制約への影響は限定的であると述べている。

この研究は、観測精度が高まるにつれて、銀河 - ハロの関係を理解する上で、系統的不確実性を制御するためにポアソン仮定を超えて進化させる必要があることを強調している。著者らは、ビスペクトルが頑健である可能性がある一方で、分散パラメータ $\nu$ を制約するためには CiC や銀河多重体などの高次統計が不可欠であると示唆している。論文は、将来の宇宙論的測定のための堅牢な系統的誤差予算を導出するために、そのような拡張を取り入れることが必要であると結論づけているが、ビスペクトルを用いた宇宙論的パラメータ推論への影響の詳細な解析は今後の課題として残されていると控えめに注記している。