A Template-Based Search for Large-Scale-Structure--Correlated Anisotropy in the Nanohertz Gravitational-Wave Background Using the Public NANOGrav 15-Year Data Set

この論文は、パルサータイミングアレイ（PTA）を用いたナノヘルツ重力波背景放射（SGWB）の解析において、観測された大規模構造（LSS）と相関する異方性成分を検出するための新しいテンプレートベースの手法を提案し、NANOGrav の 15 年データセットを適用した結果を報告するものです。以下に、問題設定、手法、主要な貢献、結果、そして意義について詳細にまとめます。

1. 問題設定 (Problem)

近年、PTA による解析でナノヘルツ帯の共通スペクトル過程（SGWB 候補）の証拠が報告されています。この信号の起源として最も有力視されているのは、連星超大質量ブラックホール（SMBHB）の集団ですが、宇宙初期の物理過程（相転移や宇宙欠陥など）が起源である可能性も残されています。

• 異方性の重要性: SMBHB に由来する SGWB は、銀河やダークマターハローの分布（大規模構造）を追跡するため、統計的に等方性ではなく、特定の方向に偏った異方性を示すことが期待されます。一方、宇宙初期起源の SGWB は高次まで等方性である可能性が高いです。
• 既存の課題: 従来の異方性検索は、球面調和関数展開による一般的な異方性モデルや、モデル非依存のマッピング手法に依存していました。しかし、有限数の SMBHB に起因する「ショットノイズ（ポアソン揺らぎ）」により、特定の銀河分布テンプレートとの相関が埋もれて検出が困難になる可能性があります。
• 本研究の目的: 観測された銀河カタログ（LSS）を直接テンプレートとして用い、SGWB の異方性成分が「現在の宇宙の物質分布」と相関しているかどうかを、ベイズ推論を用いて厳密にテストすることです。

2. 手法 (Methodology)

本研究では、外部の銀河サーベイデータを PTA の相関関数（Overlap-Reduction Function: ORF）に組み込む新しいベイズ枠組みを構築しました。

• LSS テンプレートの構築:

  • 2MASS 測光赤方偏移カタログ（2MPZ）を使用し、近傍宇宙（$z < 0.2$）の銀河過密度分布をマッピングしました。
  • 赤方偏移を 2 つの帯域（$0 < z \le 0.1$ と $0.1 < z \le 0.2$）に分割し、各帯域で銀河過密度場 $\delta_{gal}$ を構築しました。
  • 天球の欠損（銀河面マスク）を考慮し、低次数の球面調和関数係数（$\ell \le 12$）を MAP（最大事後確率）推定法で再構成しました。これにより、低多極子モードにおける LSS 情報を抽出しました。

• 混合 ORF モデル:

  • PTA のパルサー対間のクロススペクトルを、等方性の Hellings-Downs (HD) 曲線と、LSS テンプレートに比例する異方性成分の線形重ね合わせとしてモデル化しました：
    $$\Gamma_{ab} = \Gamma^{HD}_{ab} + \sum_i \epsilon_i \Gamma^{LSS(i)}_{ab}$$
  • ここで、$\Gamma^{LSS(i)}_{ab}$ は $i$ 番目の赤方偏移帯域の銀河過密度マップに基づいて計算された ORF テンプレートであり、$\epsilon_i$ はその異方性成分の振幅（混合係数）を表します。$\epsilon_i = 0$ は純粋な等方性（HD）を意味します。

• ベイズ推論:

  • NANOGrav 15 年データセット（67 個のパルサー）を解析対象としました。
  • 白色ノイズパラメータは固定し、赤色ノイズと SGWB の振幅・スペクトル指数を同時推定しつつ、$\epsilon_i$ の事後分布をサンプリングしました。
  • 単一帯域モデル（$\epsilon_1$ または $\epsilon_2$ のみ）と、二重帯域モデル（$\epsilon_1$ と $\epsilon_2$ を同時に推定）の両方を検討しました。
  • 証拠比（ベイズファクター）を用いて、LSS 相関モデルと標準的な HD モデルの優劣を比較しました。

3. 主要な貢献 (Key Contributions)

• 初の直接埋め込みアプローチ: 観測された銀河分布（LSS）を、PTA の尤度関数における ORF 成分として直接埋め込んだ初めてのベイズ解析を行いました。これにより、理論的な予測に依存せず、実際の宇宙の構造と SGWB の相関をテストする道筋を開きました。
• 低赤方偏移スライスによる解析: 2MPZ データを用いて、$z < 0.2$ の領域を 2 つの赤方偏移スライスに分割し、各スライスごとの LSS 追跡成分の制約を定量化しました。
• 再現可能なフレームワークの提示: 銀河カタログから LSS テンプレートを構築し、PTA 解析に組み込むための完全なパイプライン（マスク処理、球面調和関数再構成、ORF 計算、ベイズ推論）を提示しました。

4. 結果 (Results)

NANOGrav 15 年データセットに対する解析結果は以下の通りです。

• 統計的有意性の欠如:

  • 単一帯域解析および二重帯域解析の両方において、LSS 相関成分の振幅 $\epsilon_i$ は統計的にゼロと矛盾しません。
  • 推定値（中央値と 68% 信頼区間）は以下の通りでした：
    • 単一帯域 ($0 < z \le 0.1$): $\epsilon_1 = 0.20^{+1.68}_{-1.66}$
    • 単一帯域 ($0.1 < z \le 0.2$): $\epsilon_2 = -0.11^{+2.04}_{-1.83}$
    • 二重帯域モデル: $\epsilon_1 = 0.37^{+1.80}_{-1.91}$, $\epsilon_2 = -0.20^{+2.26}_{-2.14}$
  • 誤差範囲が広大であり、$\epsilon_i = 0$ を排除する証拠は見つかりませんでした。

• ベイズファクター:

  • 異方性モデル（HD + LSS）と等方性モデル（HD のみ）を比較したベイズファクターはすべて 1 未満でした（例：$B_{HD+LSS1, HD} = 0.40$）。
  • これは、現在のデータでは LSS 相関成分を追加する必要性（モデルの複雑さに見合う証拠）がないことを示しており、単純な等方性 HD モデルが支持されています。

• パラメータ間の相関:

  • 二重帯域モデルでは、$\epsilon_1$ と $\epsilon_2$ の間に強い負の相関（$\rho \approx -0.61$）が観測されました。これは、現在の PTA のパルサー配置（天球上のサンプリングが粗く不均一）では、2 つの低赤方偏移テンプレートが直交せず、線形結合としてしか区別できないためです。

• SGWB パラメータへの影響:

  • LSS 成分を許容しても、SGWB の振幅やスペクトル指数の推定値には大きな変化は見られませんでした。

5. 意義と結論 (Significance and Conclusion)

• 現状の限界と将来展望:
  • 本研究は、現在の PTA データでは、SMBHB に由来する LSS 追跡異方性を検出するには感度が不足していることを示しました。検出には、より多くのパルサー、より均一な天球カバレッジ、および高 SNR のデータが必要です。
  • 得られた結果は、SGWB の異方性成分が 2MPZ で追跡される低赤方偏移の銀河分布と相関している可能性に対する「上限値」を設定したことになります。
• 手法論的進展:
  • 観測された LSS 情報を PTA 解析に統合する確立された手法を提供しました。これは将来の PTA データ（より多くのパルサー、より長い観測期間）や、より深さのある銀河サーベイ（例：LSST, Euclid）と組み合わせることで、SGWB の起源（天体物理的か宇宙論的か）を決定づける強力なツールとなり得ます。
• 結論:
  • NANOGrav 15 年データを用いた解析では、LSS と相関する SGWB 異方性の統計的有意な証拠は見つかりませんでした。しかし、この研究で確立されたフレームワークは、次世代の PTA 解析において、SGWB が宇宙の物質分布（コズミックウェブ）を追跡しているかどうかを検証するための重要な基盤となります。

総じて、この論文は「SGWB の異方性が銀河分布と一致するか」という具体的な仮説を検証するための堅牢な方法論を提案し、現状のデータではその証拠は見つからなかったが、将来の検出に向けた道筋を示した重要な研究です。