Radio sirens: inferring $H_0$ with binary black holes and neutral hydrogen in the era of the Einstein Telescope and the SKA Observatory

技術的サマリー：電波サイレン：連星ブラックホールと中性水素を用いた $H_0$ の推定

問題提起
重力波（GW）宇宙論は、宇宙の膨張率を測定するために「標準サイレン」に依存している。連星ブラックホール（BBH）合体からの GW シグナルは光度距離の直接測定を提供するが、宇宙論パラメータの抽出には独立した赤方偏移情報が必要である。「明るいサイレン」（電磁波対応天体を持つ事象）は極めて稀であり、これまで一度しか観測されていない。したがって、この分野は統計的に赤方偏移を推定する「暗いサイレン」に依存している。現在の暗いサイレン手法は、通常、視線方向（LoS）の赤方偏移事前分布を構築するために銀河カタログを利用するか、相互相関のような要約統計量を採用している。しかし、これらのアプローチは、カタログの完全性、フラックス限界、特に高赤方偏移における銀河光度関数の複雑なモデリングの必要性という課題に直面している。本論文は、離散的で不完全な可能性のある銀河カタログに依存することなく、宇宙の大規模構造（LSS）を活用する堅牢な分光学的赤方偏移事前分布の必要性に対処する。

手法
著者は、アインシュタイン望遠鏡（ET）と平方キロメートルアレイ観測所（SKAO）からの模擬観測を相乗させる「電波サイレン」と呼ばれる新たな枠組みを提案する。この手法は、重力波データと中性水素（HI）強度マッピング（IM）データを組み合わせる階層ベイズ推論（HBI）枠組みを採用している。

1. シミュレーション:

   • ダークマター（DM）および HI マップ: 高解像度のダークエネルギーおよび巨大ニュートリノ宇宙（DEMNUni）シミュレーションを用いて、架空の LSS マップを構築した。これらは、DM 場に対する線形バイアス関係（$b_{HI}=1$）を仮定して HI 密度揺らぎマップ（$\delta_{HI}$）に変換された。これらのマップは、SKA-Mid の単一パラボラモードの能力に対応し、赤方偏移 $z \sim 3$ までの $20,000 \text{ deg}^2$ の天域をカバーしている。
   • GW 事象: ET 向けに $10^4$ 件の BBH 合体事象をシミュレーションした。合体率は、重力波が基礎的な物質分布を追跡すると仮定（$b_{GW}=1, \alpha_{GW}=0$）し、DM 密度揺らぎの関数としてモデル化された。信号対雑音比（SNR）が $> 150$ の事象のみを含むようにフィルタリングされ、約 3,100 件の検出可能事象が得られた。2 つのカタログが生成された。「クラスター化 GW（Clustered GWs）」（DM 場に従って分布）と「等方性 GW（Isotropic GWs）」（一様に分布）である。

2. 推論枠組み:

   • 分析には、階層ベイズ分析を実行する icarogw ソフトウェアが使用された。
   • 中核的な革新は、3 次元 HI 密度場を使用して GW 事象の LoS 赤方偏移事前分布を構築することである。離散的な銀河に確率を割り当てるのではなく、この手法は連続的な HI 密度揺らぎを用いて、特定の視線方向に沿った特定の赤方偏移に GW 源が存在する確率を重み付けする。
   • 分析では、HI 事前分布の制約力を孤立させ、特定の天体物理的質量モデルへの依存を避けるため、源固有の質量情報（スペクトルサイレン手法）を明示的に除外した。

主な貢献

• 電波サイレンの概念: 本論文は、暗いサイレンに対する分光学的赤方偏移事前分布として HI 強度マッピングを利用する「電波サイレン」アプローチを導入し、検証した。これは、フラックス制限された銀河カタログに固有の完全性の問題を回避する、LSS のトモグラフィックな視点を提供する。
• 概念実証の検証: 本研究は、ET による BBH 検出と SKAO による HI マップを組み合わせることで、宇宙論パラメータに対する制約を大幅に強化できることを実証した。
• 系統誤差リスクの特定: 著者は、仮定されたクラスター化モデルに対する手法の感度を厳密にテストした。クラスター化された GW 分布と等方性 GW 分布を比較することで、手法の有効性は GW 集団が HI 場を正確に追跡しているかどうかに依存することを浮き彫りにした。

結果

• $H_0$ への制約: 「クラスター化 GW」カタログを HI 事前分布とともに分析した際、ハッブル定数は 68% 信頼度で $H_0 = 65.8^{+4.9}_{-5.2} \text{ km s}^{-1} \text{ Mpc}^{-1}$ と制約された。これは約 8% の精度を表す。
• 改善率: この結果は、HI 情報なしで同じクラスター化 GW カタログを分析した場合（$H_0$ に対する事後分布はほとんど情報を持たないことが判明）と比較して、約 90% の改善を構成する。
• バイアス感度: 「等方性 GW」シナリオ（GW が HI 場を追跡しない場合）において、HI 事前分布を適用すると、$H_0$ を過大評価する顕著なバイアスが生じた。これは、GW と LSS の間の仮定されたクラスター化関係が正しくない場合、手法はバイアスされた結果をもたらすことを確認するものである。
• 膨張履歴: HI データの組み込みにより、$z \sim 3$ までの後期宇宙の膨張履歴 $H(z)$ の再構築が可能となった。制約は HI 密度揺らぎが最も高い低赤方偏移（$z \lesssim 0.5$）で最も厳密であるが、この手法は主に膨張履歴の詳細な形状（例：$\Omega_{m,0}$ の制約は依然として弱い）ではなく、全体の正規化（$H_0$）を制約する。
• 自己整合性: 枠組みは、バイアスパラメータ（$b_{GW}, \alpha_{GW}$）を変動させることを許可された場合、注入されたパラメータを正常に回復し、宇宙論およびクラスター化の同時推論に対して手法が自己整合的であることを示した。

意義と主張
本論文は、第 3 世代 GW 検出器および次世代電波望遠鏡の時代において、高精度宇宙論への強力かつ独立した道筋として「電波サイレン」概念の基本的な実現可能性を確立すると主張している。著者は、このアプローチが銀河カタログに必要な完全性補正なしに、宇宙の 3 次元構造を利用する物理的に動機付けられた赤方偏移事前分布を提供することを強調している。

ただし、著者は現在の限界に関しては控えめな立場を維持している。本研究は、以下の楽観的な仮定に基づく概念実証であることを認めている。

• HI による DM の完全な追跡（$b_{HI}=1$）および確率論的要素の欠如。
• GW パラメータ推定に対するフィッシャー行列近似の使用。
• HI マップにおけるビーム効果、前景除去、較正問題などの観測的系统誤差の無視。
• BBH 質量スペクトル情報の除外。

論文は結論として、電波サイレン手法が約 3,000 件の高 SNR 事象を用いて $H_0$ を約 8% の精度で制約する可能性を示しているが、今後の研究では、BBH 環境に対する現実的な天体物理モデル、IM における観測的系统誤差、および質量スペクトルとクラスター化の同時モデリングに取り組む必要があると述べている。