Cross-Correlating the Universe: The Gravitational Wave Background and Large-Scale Structure

技術要約：宇宙の相互相関：重力波背景放射と大規模構造

問題提起
パルサー・タイミング・アレイ（PTA）によって最近検出された、ナノヘルツ帯の重力波背景放射（GWB）の性質は、宇宙論における極めて重要な未解決問題である。主要な起源は、銀河合体に伴う超大質量ブラックホール連星（SMBHB）からの信号の重ね合わせであると仮定されているが、その起源には宇宙論的なソースが含まれる可能性もある。鍵となる課題は、GWBが大規模構造（LSS）を追跡するもの（SMBHBによるもの）である場合と、一様で等方的なソースの分布によるものである場合を、どのように区別するかである。さらに、現在のPTAの感度は、信号を支配してしまう「うるさい（loud）」個別のソースによって制限されており、それが宇宙網（cosmic web）によって誘起される微細な異方性を覆い隠してしまう可能性がある。本論文では、GWBの異方性と銀河分布を相互相関させることで、LSSの刻印を分離し、異なるSMBHB個体群モデルを区別できるかどうかを検証する。

手法
著者らは、純粋な理論的予測ではなく、宇宙論的分散（cosmic variance）とソースの実現における確率的性質を明示的に考慮した、厳密なシミュレーションベースのフレームワークを採用している。

1. シミュレーション・パイプライン:

   • LSSトレーサー: 著者らは、N体シミュレーションのAbacusSummitスイートを用いて、ハロー・ライトコーンを生成している。銀河は、観測された銀河のクラスタリングに一致するように較正されたAbacusHOD（ハロー占有分布）モデルを用いて配置されている。恒星質量はUniverseMachineフレームワークを用いて割り当てられている。
   • SMBHB個体群: 銀河には、NANOGrav PTAの測定値に較正された半解析的な主要合体モデルを用いて、SMBHBが配置されている。これにより、固有の分散を定量化するために1,000個の独立したSMBHB個体群の実現（realization）を生成している。
   • マップ構築: 銀河の数密度とGWBの特性歪み（$h_c$）の全天マップを作成する。GWBマップは、分解能閾値（$h_{res}$）を超えるソースを除外した上で、すべてのソースの周波数ビン内での歪みの寄与を合計することによって構築される。

2. 統計的枠組み:

   • 角パワースペクトル: 著者らは、銀河の過密場とGWBの過密場を球面調和係数（$a_{\ell m}$）を用いて解析し、角パワースペクトル（$C_\ell$）および相互スペクトル（$C^{gal, GW}_\ell$）を算出する。
   • 識別性指標（$D$）: 特定の実現（例：LSSに起因する異方性を持つもの）が、帰無仮説（一様な分布）と区別可能かどうかを判断するために、識別性統計量$D$を定義する。これは、モデルの実現が持つ固有の分散によって正規化された上で、ある実現の$\chi^2$を基準モデルと比較するものである。
   • ソース除去: 「うるさい」ソース（歪み $h > h_{res}$ を持つソース）が、異方性信号を支配するポアソン的なノイズとして機能するのを防ぐため、これらを除去することが重要なステップとなる。著者らは、残りのGWB場の統計量がウィシャート分布（ガウス型ランダム場と一致する分布）に従うようになる閾値（$h_{res}$）を決定するために、コルモゴロフ・スミルノフ検定を実施している。

主な貢献

• 「うるさい」ソースの影響の定量化: 本論文は、現在のおよび近い将来のPTAの感度において、分解されていない「うるさい」ソース（歪み $h \gtrsim 10^{-15}$）がGWBの異方性信号を支配することを実証している。この領域では、角パワースペクトルの分布はウィシャート分布から大きく逸脱し、標準的な$C_\ell$統計では場の特徴を捉えるには不十分となる。
• 相互相関の必要性: 本研究は、GWBの自己相関だけでは、LSSに基づいたSMBHBの分布と一様な分布を区別するには不十分であることを確立している。LSSの刻印を検出するための統計的な力は、GWBと銀河トレーサーとの間の相互相関に完全に依存している。
• 宇宙論的分散の扱い: 従来の理論的予測とは異なり、本研究は個々の実現を明示的にシミュレートすることで、宇宙論的分散によって生じる「散らばり（scatter）」を定量化し、単一の宇宙における検出可能性に関するより現実的な推定値を提供している。

結果

• 解像度の要件: LSSの刻印を識別する能力は、PTA実験の最大角解像度（$\ell_{max}$）に強く依存する。
  • $\ell_{max} \ge 20$ で、識別性は $1\sigma$ に達する。
  • $\ell_{max} \ge 30$ で、識別性は $2\sigma$ に達する。
  • $\ell_{max} \ge 42$ で、識別性は $3\sigma$ に達する。
  • $\ell_{max} \ge 72$ で、識別性は（楽観的な設定において）$5\sigma$ に達する。
• ソース除去閾値: これらの結果を得るためには、歪み $h \gtrsim 10^{-16}$ のソースを分解し、除去しなければならない。この閾値を下回ると、残りの未分解の個体群はLSSの公正なトレーサーとして機能し、場の統計量はガウス型の期待値と一致するようになる。
• バイアス依存性: 本研究は、識別性係数 $D$ がGWBバイアス（$b_{GW}$）に対して単調にスケールすることを確認している。クラスタリング強度がより高いほど、相互相関信号は強くなる。

意義および主張
著者らは、彼らのアプローチが、GWBを大規模構造のトレーサーとして用いるための新たな道を切り開くものであると主張している。主な意義は、相互相関がGWBの天体物理学的起源（SMBHB）を示す「決定的な証拠（smoking gun）」となり、個体群モデルを検証できることを示した点にある。

本論文は、これらの結果が、理想的なマップ再構成と完全なソース除去を仮定した「理論的な情報量」を表していることを控えめに述べている。実際のPTA実験における感度は、観測機器のノイズ、タイミングの不規則性、およびマップ作成アルゴリズムの精度に依存することになるが、これらは今後の課題である。しかし、本フレームワークは、必要な$\ell_{max}$に到達するために必要なパルサー密度と等方的な空の被覆を提供するであろうSKAOのような次世代の実験を用いて、GWBの性質をテストし、SMBHB個体群モデルを制約するための堅牢な手法を提供している。