The NANOGrav 15 yr Data Set: Impacts of Customized Chromatic Noise Models on Gravitational Wave Analyses

技術要約：NANOGrav 15年データセット：カスタマイズされた色収差ノイズモデルが重力波解析に与える影響

問題提起
パルサー・タイミング・アレイ（PTA）は、ナノヘルツ重力波（GW）を検出するために、ミリ秒パルサーの回転安定性に依存している。しかし、PTAのデータセットは、固有のスピンノイズ、分散量（DM）の変動による星間物質（ISM）の変化、太陽風の影響、およびタイミングモデルの誤差を含む、多くのノイズ源によって汚染されている。これらの色収差的（周波数依存的な）ノイズ源を誤ってモデリングすることは、共通スペクトル過程の推論にバイアスを与え、確率的重力波背景放射（GWB）のヘリングス・ダウンズ（HD）相関シグネチャを増幅または抑制させる可能性がある。従来のNANOGrav 15年（NG15）データセットの解析では、標準的なノイズモデル（DM変動に対するDMXなど）が使用されていたが、近年の研究では、カスタマイズされた色収差ノイズモデル（CNM）がこれらの系統誤差をより適切に軽減できることが示唆されている。本論文では、CNMを用いることがGWBの有意性、スペクトル特性、ならびに連続波（CW）および代替重力波偏極の探索にどのように影響するかを調査する。

手法
著者らは、Larsen, Baier et al. (2026) によって開発されたCNMを用いてNG15データセットを再解析した。主な手法の構成要素は以下の通りである：

• ノイズ・モデリング： 分散量に関する標準的な「DMX」モデルを、ガウス過程（DMGP）モデルに置き換えた。さらに、CNMフレームワークは、各パルサー（67個）に対して、ベイズ証拠に基づき、カスタマイズされた色収差信号（自由な色収差変動、確率的太陽風、過渡事象など）をテストし、適用する。
• 尤度実装： パルサーごとの色収差モデルに伴う計算コストに対処するため、著者らは「キャッシュされた色収差尤度（Cached-Chromatic Likelihood）」を実装した。この手法は、2段階の周辺化プロセスを$N+2$段階のプロセスへと拡張するものであり、色収差パラメータ（単一パルサー解析における最大事後確率値に固定）を、メインのGW解析の前に周辺化することを可能にする。これにより、従来の構成と比較して、尤度評価時間は、共通の無相関赤ノイズ（CURN）に対しては約2倍、HD相関に対しては約4倍短縮された。
• 解析範囲： 本研究は、Agazie et al. (2023a, 2023d, 2024) および Agarwal et al. (2025) による主要なGW解析を再現した。これには以下が含まれる：
  • GWBのスペクトル特性（べき乗則、破断べき乗則、および自由スペクトル）。
  • HD相関の有意性（ベイズ因子および最適統計量）。
  • スカラー横波（ST）偏極モードの探索。
  • 全天およびターゲットを絞った連続波（CW）探索。
  • GWBの天体物理学的および宇宙論的解釈。

主な貢献と結果

1. 強化されたHD相関の有意性：

   • 14個のフーリエモードを用いたCNMを使用した場合、CURN過程に対するHD相関のベイズ因子は、約185（Agazie et al. 2023aにおける値）から $1571 \pm 14$ へと、約8倍増加した。
   • 最適統計量を用いた頻度論的解析では、信号対雑音比（S/N）が5.0から5.3へと向上した。
   • ベイズ $p$ 値は、$7.9 \times 10^{-4}$ ($3.16\sigma$) から $4.5 \times 10^{-4}$ ($3.32\sigma$) へと減少した。
   • 著者らは、この増加は新しいデータの追加ではなく、特に10–30 nHzの周波数領域におけるノイズ軽減によるものであるとしている。

2. スペクトル特性：

   • 推定されたGWB振幅は、$A_{\text{GWB}} = 2.1^{+0.6}_{-0.5} \times 10^{-15}$（$\gamma = 13/3$ で固定）へと減少し、以前の推定値から約12%減少した。
   • 変動-$\gamma$解析において、スペクトル指数は $\gamma_{\text{GWB}} = 3.5^{+0.7}_{-0.6}$ へと上昇し、超大質量ブラックホール連星（SMBHB）の理論的期待値である $13/3$ に近づいた。
   • Agazie et al. (2025a) で報告された、16 nHzにおけるパワースペクトルの顕著な偏りは、PSR J1713+0747に関連するノイズモデル依存のアーティファクトであることが判明し、CNMの下ではもはや有意ではない。

3. 代替偏極および連続波：

   • 非アインシュタイン的なスカラー横波（ST）モードの証拠は減少した。HDに対するSTのベイズ因子は、約2.3から約158へと増加し、ST振幅の95%上限は減少した（約26%減少）。
   • 全天CW探索において、低周波の候補（$\sim4.2$ nHz）の有意性は、HD相関を持つ背景（$B_{\text{HD+CW}}^{\text{HD}} = 1.69 \pm 0.06$）に対して統計的に有意ではないままとなった。
   • 特定のAGN候補（SDSS J0729+4008 および SDSS J1536+0441）を対象としたターゲットCW探索では、有意性が低下し、ベイズ因子は一般的な候補集団と一致するレベルまで減少した。
   • 中・高周波数帯におけるホワイトノイズフロアの減少により、NANOGravアレイのCW探索における検出体積は3.2倍増加した。

4. 赤ノイズ・バジェットと検出器感度：

   • 共通過程とともに有意な固有赤ノイズ（IRN）を示すパルサーの数は、12個（Agazie et al. 2023b）から6個へと減少した。これは、以前のIRN検出が、おそらく誤ったモデリングによる色収差ノイズであったことを示唆している。
   • 感度曲線は、CNMが中・高周波数帯において大幅な感度向上を提供することを示しており、最も感度の高い周波数帯を、約5 nHz（DMX）から約6.5 nHzへとシフトさせている。

5. 天体物理学的および宇宙論的解釈：

   • 天体物理学的（SMBHB）： 更新されたスペクトル特性は、SMBHB起源と矛盾しない。現象論的な集団モデルへの適合では、わずかな変化しか見られず、銀河の特性質量および内部硬化スロープのわずかな低下傾向が示された。
   • 宇宙論的： 宇宙インフレーション、一次相転移、および安定した宇宙ひもに関する再解釈では、パラメータの後ろ分布にわずかな定量的変化しか見られなかった。宇宙ひもについては、最良適合テンション $G\mu$ が約35%減少し、これはひもを生成する相転移のエネルギースケールがより低いことを意味する。

意義と主張
本論文は、CNMの主な影響は色収差ノイズの改善された軽減であり、それがHD相関シグネチャのより堅牢で統計的に有意な検出につながっていると主張している。著者らは、GWBのスペクトル特性が現在「テンションが少ない（less tension）」状態であり、振幅は以前よりも静かで急峻になっていると述べている。

決定的な点として、著者らは、GWBの統計的証拠は強まったものの、天体物理学的および宇宙論的な解釈は概ね変わっていないことは、それらが現在「事前分布支配的（prior dominated）」であるためであると強調している。本論文は、詳細な色収差のモデリングが、将来のデータセットにおいて検出感度を最適化し、ナノヘルツ重力波背景放射を正確に特徴付けるために不可欠であることを結論付けているが、GWBシグナルと理論モデルとの間のすべてのテンションを解決したと主張するものではない。本研究は、将来のNANOGrav解析のためのCNMフレームワークの検証として機能するものである。