Dark siren cross-correlations and the sensitivity of $H_0$ to methodological choices

技術的サマリー：ダークサイレンのクロス相関と方法論的選択に対する$H_0$の感度

問題提起
重力波（GW）源は「標準サイレン」として機能し、絶対距離測定を提供することで、ハッブル定数（$H_0$）を制約し、宇宙論パラメータにおける現在の緊張関係を解決するための有望な独立した手段を提供する。「明るいサイレン」（電磁波対応天体を持つもの）と「ダークサイレン」（GW 局所化を銀河カタログと統計的に関連付けるもの）は著しく進展してきたが、ダークサイレンの方法論は、不完全な銀河カタログ、系統的バイアス、大規模構造（LSS）のモデリングに伴う計算の複雑さという課題に直面している。

本論文は、GW イベントと銀河カタログを共有する大規模構造を通じて相関させる角クロス相関法の堅牢性を、様々な方法論的選択に対して検証する。具体的には、共分散の扱い、バイアスのパラメータ化、距離・赤方偏移のビン分割、およびカタログの不完全性が推定される$H_0$値にどのように影響するかを調査する。目的は、適切なモデリングを通じて系統的バイアスを効果的に軽減できるかどうかを判断し、将来の精密宇宙論に対するクロス相関アプローチの有効性を検証することである。

方法論
著者らは、赤方偏移と光度距離のビンにわたる GW イベントと銀河のクラスター化を解析する、トモグラフィックな角クロス相関フレームワークを採用している。

• 形式論： この手法は、GW トレーサーと銀河トレーサー間の角クロスパワースペクトル（$C_{\ell}$）を利用する。計算コストと次元を削減するため、著者らは多重極$\ell$にわたるクロスパワースペクトルの重み付き和として定義される圧縮統計量$S_{ij}$を導入する。この統計量は、2 つのトレーサー間の視線方向のクロス相関を効果的に測定する。
• シミュレーション： 解析は、平坦な$\Lambda$CDM 宇宙論（$H_0 = 70$ km s$^{-1}$ Mpc$^{-1}$、$\Omega_m = 0.25$）の下で生成された 600 個の対数正規モックカタログに基づいている。これらのモックは赤方偏移範囲$0 \le z \le 0.5$をカバーし、等級カットや天空局所化の不確かさなど、現実的な選択関数を組み込んでいる。
• パイプライン： 理論的予測は、Limber 近似の下で CAMB を用いて計算され、選択効果（マルキストバイアスとエディントンバイアス）のフォワードモデリングを投影カーネルに直接組み込んでいる。パラメータ推論は、$H_0$およびバイアスパラメータ（$b_{gal}, b_{GW}$）をフィットするマルコフ連鎖モンテカルロ（MCMC）によって行われる。
• 共分散推定： 本研究は、ガウス仮説から導出された解析的共分散行列と、600 個のモック実現から推定された経験的共分散行列を比較し、有限サンプルサイズを考慮してハートラップおよびパーシバル補正を適用している。

主要な貢献と結果

1. 圧縮統計量と非圧縮統計量：
   著者らは、完全な角パワースペクトルを要約統計量$S_{ij}$に圧縮することが、宇宙論情報の大部分を保持することを示している。$S_{ij}$を使用すると制約力がわずかに低下し（200 個のイベントの場合、完全な$C_{\ell}$と比較して$H_0$の不確かさが約 1.4 倍、1000 個のイベントの場合約 1.5 倍に増加する）が、ベッセル関数の直接積分を回避することで、計算上の利点と数値的安定性を大幅に提供している。

2. 共分散モデリング：
   重要な発見として、パラメータ制約が共分散行列の選択に敏感であることが示された。

   • 大規模サンプル（$\sim 10^3$イベント）の場合、非対角項を含めれば、解析的共分散とモック共分散は一致した結果を与える。
   • 小規模サンプル（200 個のイベント）の場合、解析的共分散はモックから導出された共分散と比較して、不確かさを約 1.5 倍過小評価する。さらに、小規模サンプルの解析的共分散に非対角項を含めると、MCMC 反転において数値的不安定性を引き起こし、収束を妨げる。
   • 本研究は、過度に楽観的な制約を避けるために、特にイベント数が少ない場合、モックに基づく共分散推定が不可欠であると結論づけている。

3. 系統的不確かさ：

   • マルキストおよびエディントンバイアス： 論文は、光度距離の不確かさがエディントンバイアスおよびマルキストバイアスを誘発し、モデル化されていない場合、$H_0$を系統的に高くバイアスさせることを確認している。しかし、これらの効果は、理論的カーネル内で検出確率と選択閾値をフォワードモデリングすることで効果的に軽減される。
   • バイアスパラメータ化： 本研究は、銀河および GW バイアスに対して定数、線形、および 2 次パラメータ化をテストしている。定数パラメータ化は赤方偏移進化を捉えることができないためバイアスを導入するが、線形および 2 次形式はバイアスのない$H_0$推定値を回復する。これは、パラメータ化が十分に柔軟であれば、近未来の観測においてバイアス進化の正確な知識が厳密には必要ではないことを示唆している。
   • カタログの不完全性： 特定された主要な利点は、不完全なカタログに対する手法の堅牢性である。欠落集団のモデリングを必要とする標準的なダークサイレンアプローチとは異なり、クロス相関法は、選択関数を理論的予測で一貫してフォワードモデル化することのみを要求する。これにより、観測された過密度場を直接使用することが可能となり、半径方向の選択効果（例えば、等級制限）があってもバイアスのない$H_0$制約を得ることができる。ただし、天空全体で不完全性が異なるカタログを組み合わせると、グローバルなバイアスパラメータを制約する力が低下し、制約が劣化する可能性がある。

4. ビン分割の影響：
   赤方偏移および光度距離のビン分割の選択は、信号分解能とショットノイズの間のトレードオフを伴う。

   • より細かい赤方偏移ビン分割は、銀河のショットノイズが高分解能でも低く保たれるため、$H_0$制約を大幅に改善する。
   • より細かい光度距離ビン分割は、GW イベントの希薄さにより、限界効用が現れる。これらのビンを分割するとショットノイズが増加し、高分解能の利点を相殺する可能性がある。
   • 最適な戦略は、典型的な GW 距離不確かさと同等かそれ以上のビン幅を採用し、過度の信号希薄化なしに測定誤差を吸収することである。

意義と主張
本論文は、方法論的選択を慎重に管理すれば、角クロス相関法がバイアスのない$H_0$制約を達成できる強力かつ補完的なツールであると主張している。著者らは以下を強調している。

• 距離の不確かさ、バイアス進化、およびカタログの不完全性に起因する系統的バイアスは、欠落集団の明示的な再構成ではなく、フォワードモデリングを通じて効果的に軽減できる。
• 解析的近似が失敗する小規模イベントサンプルの場合、特に信頼性の高い不確かさ推定のために、モックから導出された共分散行列の使用が不可欠である。
• 圧縮統計量$S_{ij}$は、情報の最小限の損失で完全なパワースペクトル解析に代わる計算効率の高い代替手段を提供する。

最終的に、この研究は、適切なモデリングと十分な数の精密な GW イベントのサンプルがあれば、クロス相関法がハッブル緊張に堅牢に対処できることを確立している。これは、選択効果を理論的枠組みに自然に組み込むことで、標準的な銀河カタログアプローチに対する明確な利点を提供する。