A model-independent measurement of the Hubble constant from gravitational-wave standard sirens and electromagnetic observations

技術要約：重力波標準サイレンと電磁波観測に基づく、モデルに依存しないハッブル定数の測定

問題提起
本論文は、「ハッブル・テンション（Hubble tension）」に対処している。これは、$\Lambda$CDMモデルの下での宇宙マイクロ波背景放射（CMB）から導出された初期宇宙のハッブル定数（$H_0 = 67.4 \pm 0.5$ km s$^{-1}$ Mpc$^{-1}$）と、Ia型超新星（SNe Ia）および局所的な距離梯子を用いた後期宇宙の直接測定（$H_0 = 73.04 \pm 1.04$ km s$^{-1}$ Mpc$^{-1}$）との間の、約$5\sigma - 6\sigma$の有意な不一致である。従来の測定手法は、特定の宇宙論モデルや音響地平線のスケールに関する仮定に依存することが多い。著者らは、これら理論的なバイアスを回避するため、電磁波（EM）と重力波（GW）の両方のデータセットを利用し、純粋に観測データのみに基づいた、完全にモデルに依存しない$H_0$の決定を提供することを目指している。

手法
本研究では、ガウス過程（GP）回帰、すなわち非パラメトリックなベイズ法を用いて、あらかじめ定義された関数形式（例：$\Lambda$CDM）を仮定することなく、観測データから光度距離関数 $D_L(z)$ およびその一次導関数を直接再構成する。

1. フレームワーク: 著者らは Gapp Pythonコードを利用している。$D_L(z)$ を、平均関数と共分散カーネルによって定義されるガウス過程としてモデル化する。Matérnカーネルよりも安定性と導関数の振る舞いが良好であることから、二乗指数関数カーネルが選択されている。
2. $H_0$ の導出: ハッブル定数は、再構成された関数から以下の関係式を用いて解析的に導出される：
   $$H_0 = c \left[ \frac{d}{dz} \left( \frac{D_L(z)}{1+z} \right) \bigg|_{z=0} \right]^{-1}$$
   これには、$D_L(z)$ とその $z=0$ における導関数の両方を正確に再構成することが求められる。
3. データの統合: 解析には以下を組み合わせている：
   • 重力波: GWTC-3カタログ（$p_{astro} > 0.5$ の90個の候補）から得られたデータ。これらは、宇宙の距離梯子とは独立した絶対光度距離を提供する。
   • 電磁波プローブ:
     • コズミック・クロノメーター（CC）：直接的な $H(z)$ 測定値を $D_L(z)$ に変換したもの。
     • DESI BAO DR2：バリオン音響振動データ。モデルへの依存性を排除するため、音響地平線 $r_d$ は、モデルに依存しない性質を維持するために、摂動パラメータとして扱い、周辺化（marginalize）している。
     • Ia型超新星：3つのコンピレーション（Pantheon+SH0ES、Union 3.0、および DES-SN5YR）。
     • ガンマ線バースト（GRB）：高赤方偏移の指標として、アマティ関係（Amati relation）を介して使用。
4. 誤差処理: 低赤方偏移（$z < 0.05$）における特異速度誤差を考慮し、これらを観測誤差に二乗和の平方根として加算している。極めて重要な点として、GPフレームワーク内で、再構成された関数 $D_L(z)$ とその導関数 $D'_L(z)$ の間の共分散を明示的に計算・伝播させることで、統計的に一貫した不確実性推定を実現している。

主な貢献

• ハイブリッドなモデル非依存アプローチ: 本研究は、音響地平線のスケールや特定のダークエネルギーモデルに依存することなく、$H_0$ を制約するために、多様なEMデータセットと共に GWTC-3 標準サイレン・データを GP フレームワークに統合した最初期の試みの一つである。
• GWの影響の定量化: 本研究は、$H_0$ の精度に対する GW データの寄与を明示的に分離している。GWデータを含む場合と含まない場合の GP 実行結果を比較することで、GWTC-3 データを追加することにより、$H_0$ の $1\sigma$ 不確かさが（特定の組み合わせにおいて $1.02$から$0.58$ km s$^{-1}$ Mpc$^{-1}$ へと）約15–20%減少することを実証している。
• 共分散の取り扱い: GP フレームワーク内での解析的な微分の実装により、距離関数と導関数の間の相関を厳密に扱うことが可能となり、これが $z=0$ におけるバイアスを最小化するための決定的な要因となっている。

結果
モデルに依存しない再構成の結果、得られた $H_0$ の値は、Planck の結果と SH0ES の結果の間に位置するが、SH0ES の値（$H_0 = 73.04 \pm 1.04$ km s$^{-1}$ Mpc$^{-1}$）を強く支持している。

• ほとんどのデータ組み合わせにおいて、導出された $H_0$ と SH0ES の結果との間のテンションは $2\sigma$ 未満であり、統計的に有意な不一致は認められない。
• GW データの導入は、EM のみの解析と比較して、制約を大幅に強めている。
• 特定の組み合わせ（GWTC3+CC+DESI+DESY5）では、Planck の測定値をわずかに支持する傾向が見られるが、全体的な傾向としては局所的な距離梯子の結果と一致している。

意義と主張
著者らは、本研究がハッブル・テンションを決定的に解決するものではないことを控えめに主張しているが、観測入力のみを用いて SH0ES の決定を支持する、堅牢で一貫した $H_0$ 推定値を提供している。本研究の意義は以下の点にある：

1. GW データの検証: 現時点での限界（低赤方偏移のカバー範囲や統計的不確かさ）があるものの、GW 標準サイレンが、EM データと組み合わせた際に $H_0$ 測定の精度を意味のある形で向上させる絶対的な較正手段となることを示した。
2. 整合性の確認: モデルに依存しないフレームワークにおいて、GW と EM の測定値が互いに矛盾しないことを確認した。
3. 将来の有用性: 次世代の検出器によるより大規模な GW カタログや、より精密な EM データセットへと拡張可能な手法を確立しており、理論的バイアスなしにハッブル・テンションの起源を調査するための経路を提示している。

著者らは、電磁波の対応天体を持たない連星ブラックホール（BBH）イベントにおける赤方偏移推定のための集団レベルの仮定への依存といった現在の限界を認めており、これらの系統誤差の厳密な取り扱いが今後の課題であることを述べている。