Increasing the sensitivity of full-shape galaxy clustering measurements in configuration-space with three-point statistics

技術要約：3点統計を用いた構成空間におけるフルシェイプ銀河クラスタリング測定の感度向上

問題と動機
大規模構造（LSS）の標準的な宇宙論解析は、2点統計（2PCF）、特にバリオン音響振動（BAO）の特徴や線形・準線形スケールにおけるフルシェイプ・モデリングに大きく依存している。これらの手法は堅牢ではあるが、後期宇宙の重力進化、赤方偏移空間歪み、および銀河バイアスによって生成される非ガウス情報を放棄してしまう。3点相関関数（3PCF）は、これらの高次相関を捉えることができ、非線形構造を探索するための2PCFの自然な拡張となる。しかし、非線形スケールから宇宙論的情報を抽出することは、摂動論の破綻や高次元データベクトルのモデリングに伴う計算の複雑さにより、依然として困難である。さらに、銀河とハローの接続（バイアス）に関する不確実性が、小スケールにおける宇宙論的制約をしばしば不明瞭にする。本研究では、圧縮された視線方向（LOS）依存の3PCF推定器が、将来のRoman宇宙望遠鏡（Roman）やDESI（Dark Energy Spectroscopic Instrument）のような観測において、BAOの特徴よりも小さなスケール（$\lesssim 80 \, h^{-1}\text{Mpc}$）での宇宙論的制約を大幅に強化できるかどうかを調査する。

手法
著者らは、AbacusSummit N体シミュレーションスイートを利用しており、具体的には $z=1.1$ および $z=0.8$ のスナップショットを用い、**ハロー占有分布（HOD）**モデルを用いて銀河を配置している。

• シミュレーション: 基準となる宇宙論（$c000$）はPlanck 2018に基づいている。一連のオフセットシミュレーション（$c100, c102, \dots$）は、宇宙論パラメータ（$h, \omega_{cdm}, n_s, \alpha_s, \sigma_8, N_{ur}, w_0$）を変化させ、追加のシミュレーション（$c901, \dots$）はHODパラメータ（$\log M_{cut}, \sigma, \log M_1, \alpha$）を変化させている。
• 銀河サンプル: 2つの代表的なサンプルをモデリングしている：$z \approx 1.1$ におけるRoman風の放射線放出銀河（ELG）、および $z \approx 0.8$ におけるDESI風の赤色巨星銀河（LRG）。
• 統計量:
  • 2PCF: Corrfunc を用いて計算され、BAOの特徴を避けるためにスケル $2.5 < s < 82.5 \, h^{-1}\text{Mpc}$ に限定されている。モノポール（$\ell=0$）およびクアドラポール（$\ell=2$）の多重極が使用される。
  • 3PCF: 高速な構成空間三角形計数コードである TriCo を用いて測定されている。著者らは、新しいLOS依存の分解を導入している。三角形は、2つの辺の長さ（$s_1, s_2$）と、それらの辺がLOSとなす角度の余弦（$\mu_1, \mu_2$）によってビン分割される。
  • 圧縮: 次元の管理のため、$\mu$ は平行（$\parallel$）、横断（$\perp$）、混合（mixed）の3つのカテゴリにビン分割される。これにより、標準的な2PCF多重極に加えて、3つの異なる3PCF成分（$\bar{\zeta}_{\parallel}, \bar{\zeta}_{\text{mix}}, \bar{\zeta}_{\perp}$）が得られる。
• モデリング: オフセットシミュレーションを用いてデータベクトルのヤコビアンを推定することにより、基準パラメータの周りで線形応答モデルを構築する。HODパラメータは、摂動パラメータとして扱われ、周辺化される。
• 共分散: 1000個の独立した小ボックスシミュレーション（$500 \, h^{-1}\text{Mpc}$）から共分散行列を推定し、全ボックス体積（$2 \, h^{-1}\text{Gpc}$）にリスケールしている。パラメータ依存の尤度問題の複雑さを避けるため、有限のモック数に対するHartlap補正を組み込んだ固定有効共分散を使用している。

主要な貢献

1. 新しい推定器: 赤方偏移空間のクラスタリングに関連する異方性情報を保持しつつ、管理可能なデータベクトルのサイズ（フルセットで440要素）を維持する、圧縮されたLOS依存の3PCF基底の開発。
2. シミュレーションに基づく予測: オフセットシミュレーションを用いてHODの不確実性を明示的に周辺化するシミュレーションベースのフレームワークを適用し、絶対的なバイアス較正に依存することなく、3PCFの相対的な情報利得を頑健に定量化することを可能にした。
3. 検証: 追加のAbacusSummit宇宙論を用いた線形応答モデルの厳密な検証を行い、基準点付近での予測レベルの比較に対するその妥当性を確認した（ただし、最小スケールにおける限界も認めている）。

結果
Roman ELGサンプルおよびDESI LRGサンプルに対して行われたMCMCによる予測は、以下の知見を得た：

• 大幅な制約の改善: 2PCFのモノポールおよびクアドラポールと共にLOS依存の3PCFを組み込むことで、宇宙論的制約が大幅に強化される。Roman ELGサンプルの場合、$\sigma_8$ に対する制約は2PCFのみの場合と比較して約5倍改善する。他の宇宙論パラメータ（$h, \omega_{cdm}, n_s, w_0$）およびHODパラメータについても、同様の、あるいはわずかに小さい改善（2〜4倍）が観察される。
• 累積的な情報利得: この改善は、単一のスケールや特定の三角形の構成によって駆動されるものではない。累積的な $\chi^2$ 診断は、情報が三角形の構成全体に広く分布していることを示しており、この利得が非線形クラスタリング信号の累積的な情報量から生じていることを示唆している。
• LOS依存性の重要性: 3PCFのモノポール（LOS平均）に限定した解析では、利用可能な情報の一部しか捉えられない。完全なLOS依存測定は、モノポール単独よりも約2〜3倍の追加の改善を提供し、三角形の構成の向きが重要な宇宙論的情報を持っていることを実証している。
• バイアスに対する頑健性: HODパラメータを周辺化した後でも利得は持続し、これは3PCFが宇宙論パラメータと銀河バイアスの間の縮退を解くのに役立つことを示唆している。

意義と主張
著者らは、本研究が、LOS依存の3PCF測定を用いることで、フルシェップ銀河クラスタリング解析の感度を大幅に高めることが可能であることを示していると主張している。彼らは、この手法が局所的な特徴への応答ではなく、「真の高次クラスタリング効果」を提供することを強調している。

本論文は、これらの結果が理想的な周期境界ボックス体積に基づく相対的な予測であり、サーベイの幾何学、ファイバースロッティング、その他の系統誤差を考慮する必要がある特定のRomanまたはDESIのサーベイサンプルの正確な制約力を予測するものではないことを明記している。しかし、著者らは、相対的な改善（例：$\sigma_8$ に対する5倍の利得）は、将来のフルシェップ解析に3PCF情報を含める価値を示す頑健な指標であると主張している。本研究は、銀河とハローの接続を周辺化しながら、2点および3点の情報を組み合わせるシミュレーションベースの解析への一歩であり、現在モデリングが困難な非線形スケールを活用するための有望な道筋を提示している。