The NANOGrav 15 yr Data Set: Customized Chromatic Noise Models

技術要約：NANOGrav 15年データセット：カスタマイズされた色収差ノイズモデル

問題提起
パルサー・タイミング・アレイ（PTA）は、超高精度なパルス到着時刻（TOA）の測定を通じて、ナノヘルツ帯の重力波（GW）を探索している。重力波背景放射（GWB）のロバストな推論には、非GWノイズ源の包括的なモデリングが必要である。その重要なクラスが、分散量（DM）の変化、散乱、および太陽風（SW）の遅延といった、星間および惑星間伝播効果に起因する色収差ノイズである。これらの効果は周波数依存的な遅延（$\Delta t \propto \nu^{-\chi}$）をもたらし、これらがモデル化されていない、あるいは誤って特徴付けられた場合、無色（アクロマティック）な赤色ノイズのパラメータへとリークしてしまう。このリークは、GWパラメータの推定にバイアスを与え、GWB信号を模倣したり、あるいは覆い隠したりする可能性がある。従来の解析では、一様な処方箋（例：区分的に定数となるDMXモデル）に依存することが多かったが、これらは複雑でパルサー固有の伝播効果、特に太陽風の変動性と星間物質（ISM）の乱流との相互作用を捉えるには不十分な場合がある。

手法
著者らは、NANOGrav 15年（NG15）データセットに含まれる全67個のパルサーに対し、色収差ノイズモデルのパルサーごとのカスタマイズ化を提示している。その手法は以下の通りである：

1. 拡張されたガウス過程（GP）スイート： 本研究では、色収差による遅延を同時に特徴付けるために、拡張されたGPモデルのスイートを利用している。これには以下が含まれる：

   • フーリエ領域モデル： DMおよび「自由色収差（Free Chromatic; FC）」ノイズに対して、べき乗則のパワースペクトル密度（PSD）を使用する。
   • 時間領域（TD）モデル： 様々なカーネル（二乗指数関数、準周期、リッジ）を用いた線形補間基底を用い、時間相関のある変動をモデル化する。
   • TD-RFモデル： DMの変動が無線周波数間でデコレーション（非相関化）することを許容する多次元モデル。
   • 太陽風（SW）モデル： 太陽電子密度（$n_E$）に関する決定論的なグローバルな区分定数モデルと、個々の視線方向における短周期の変動を捉えるための確率的GP摂動（SWGP）を組み合わせたハイブリッド・アプローチ。

2. モデル選択ワークフロー： 各パルサーに対して厳格な多段階ベイズ・モデル選択プロセスが採用されている：

   • 初期スクリーニング： 3つのベースライン・モデルクラス（Fourier, TD, Ridge）をテストし、カーネルの複雑さと、Savage-Dickeyベイズ因子を用いたFCおよびSW成分の有意性を評価する。
   • 基底サイズの最適化： ベイズ証拠と尤度計算時間の比（$Z/t_L$）を最大化することにより、最適な基底要素の数（周波数カウント $N_f$ または補間間隔 $dt$）を選択し、適合の質と計算コストのバランスを取る。
   • 精緻化： 年周期の色収差変動、過渡的なイベント（例：指数関数的なディップ）、およびバンド依存のDM変動のテストを反復的に行う。
   • グローバル再適合： 個別のパルサーモデルを選択した後、すべての67個のパルサーにわたる因子分解された事後分布アプローチを用いて、グローバルな決定論的SWパラメータを再適合させ、一貫性を確保する。

3. 比較フレームワーク： カスタマイズされたノイズモデル（CNM）を、標準的なNG15ベースライン（DMX）および汎用的な「カスタマイズされていない」フーリエベースのモデルと比較する。指標には、無色赤色ノイズ（RN）のハイパーパラメータ、ホワイトノイズパラメータ（EFAC, EQUAD, ECORR）の変化、自由スペクトル解析、および感度曲線が含まれる。

主な貢献と結果

• カスタマイズされたノイズモデル： 本研究は、全67個のパルサーに対して、最適化された固有の色収差ノイズモデルの導出に成功した。約半数のパルサーが時間領域（TD）モデル（多くの場合RidgeまたはSEカーネルを伴う）を支持し、残りの半分はフーリエ領域モデルを支持している。この分割は、複雑な太陽風や周波数依存のDM変動の存在と相関しており、これらはTD基底によってより良く捉えられる。
• 太陽風の特性評価： ハイブリッドSWモデルは11年周期の太陽サイクルを追跡することに成功しており、推定された電子密度の変動はF10.7 cm太陽フラックスのデータと一致している。微細なSWGP成分の導入により、ISMのDMスペクトル指数が改善され、期待されるコルモゴロフ値（$\gamma_{DM} \approx 8/3$）へとシフトし、モデル化されていない太陽風の変動によるバイアスが減少した。
• 自由色収差（FC）ノイズの検出： 67個中21個のパルサーにおいて、非分散的な色収差ノイズ（FC）の証拠が見つかった。そのうち6個については、PTAの残差における初の報告となる検出である。推定された色収差指数（$\chi$）は幅広い範囲に及んでいる。一部は薄いスクリーンによる散乱の期待値（$\chi \approx 4$）と一致しているが、他のものは大きな値（$\chi > 7$）や小さな値（$\chi < 3$）を示しており、複雑な物理現象、あるいは潜在的なモデルの誤設定を示唆している。
• 無色赤色ノイズへの影響： CNMを適用することで、67個中19個のパルサーにおいて、無色赤色ノイズの推定が大きく変化した。5つのケースでは、以前は有意であった固有の赤色ノイズ（IRN）が有意ではなくなっており、これはその信号が恐らく色収差に起因するものであること（例：散乱や太陽風の影響）を示している。その他のケースでは、赤色ノイズの振幅が減少するか、あるいはスペクトル指数がシフトしている。特に、PSR J1713+0747の赤色ノイズの振幅は、色収差イベントモデルの導入によって大幅に減少し、そのスペクトル指数は増大した。
• ホワイトノイズとECORR： CNMは、特に強い散乱や太陽風感受性を持つパルサーにおいて、大きなECORR（エポック相関ノイズ）パラメータの必要性を一般に減少させる。しかし、一部のパルサーではECORRの有意性が増加しており、これは現在のSWモデルが、特に太陽接近時における短周期の異常値を捉えるための柔軟性がまだ不足していることを示唆している。
• 感度の向上： 感度曲線は、CNMが「ノイズフロア」を低減し、伝達関数を改善することにより、中・高周波数帯におけるPTAの感度を向上させることを示している。ホワイトニングされた残差は、CNMとDMXの両方において概ねガウス分布と一致しているが、特定のパルサーではCNMによって適合の質が向上している。

意義
本論文は、洗練されたパルサー固有の色収差ノイズモデリングが、低周波GW探索の感度と精度を高めるために不可欠であることを示している。色収差のパワーが無色赤色ノイズへとリークすることを軽減することで、CNMはGWBの推論におけるバイアスを減少させ、確率的背景の特性評価を改善する。これらの結果は、「ワンサイズ・フィッツ・オール（一律適用）」の処方箋が、現在のPTA領域で求められる精度には不十分であることを浮き彫りにしている。本研究は、伝播に起因する構造と真の天体物理学的信号を分離するための堅牢な枠組みを提供し、ナノヘルツ帯の重力波背景放射および連続波源の検出・特性評価に向けたNANOGravコラボレーションの取り組みを直接的に支援するものである。さらに、ノイズ成分の分解は、乱流ISM、惑星間物質、およびパルサー磁気圏の変動を研究するための新たな道筋を提供する。