[論文レビュー] The Need For Speed: Rapid Refitting Techniques for Bayesian Spectral Characterization of the Gravitational Wave Background Using PTAs

William G. Lamb, Stephen R. Taylor+1 more|arXiv (Cornell University)|Mar 27, 2023

Pulsars and Gravitational Waves Research被引用数 4

ひとこと要約

本稿では、事前に計算されたタイミングパワー スペクトルの後部推定値を用いて、パulsar 時刻配列（PTA）における重力波背景のスペクトル特徴付けを高速化するベイジアン再適合技術を導入する。この手法は、フルパイプラインの結果と10%未満のヘリンジャー距離を達成しながらも、100〜10,000倍高速であるため、大規模PTAデータセットにおける効率的なモデル選択とパラメータ推定を可能にする。

ABSTRACT

Pulsar timing arrays (PTAs) have recently found evidence for a nanohertz-frequency stochastic gravitational-wave background (SGWB). Constraining its spectral characteristics will reveal its origins. To achieve this, we must understand how data and modeling conditions in each pulsar influence the precision and accuracy of SGWB spectral recovery, typically requiring many Bayesian analyses on real data sets and large-scale simulations that are slow and computationally taxing. To combat this, we have developed several new rapid approaches that operate on intermediate SGWB analysis products. These techniques refit SGWB spectral models to previously computed Bayesian posterior estimates of the timing power spectra. We test our new techniques on simulated PTA data sets and the NANOGrav 12.5-year data set, where in the latter our refit posterior achieves a Hellinger distance -- bounded between 0 for identical distributions and 1 for zero overlap -- from the current full production-level pipeline that is < 0.1. Our techniques are ~ $10^2$--$10^4$ times faster than the production-level likelihood and scale much more favorably (sub-linearly) as a PTA is expanded with new pulsars or observations. Our techniques also allow us to demonstrate conclusively that SGWB spectral characterization in PTA data sets is driven by the longest-timed pulsars and the best-measured power spectral densities. Indeed, the common-process spectral properties found in the NANOGrav 12.5-year data set are given by analyzing only the ~14 longest-timed pulsars out of the full 45 pulsar array, and we find that the 'shallowing' of the common-process power-law model occurs when gravitational-wave frequencies higher than ~50 nanohertz are included. The implementation of our techniques is openly available as a software suite to allow fast and flexible PTA SGWB spectral characterization and model selection.

研究の動機と目的

• 大規模データセット上で繰り返しフル尤度評価を要する、パulsar 時刻配列（PTA）におけるベイジアンスペクトル特徴付けの計算的ボトル neck を解消すること。
• スペクトルモデルフィッティングの精度を保ちつつ、フルベイジアン推論の代替として高速でスケーラブルな手法を開発すること。
• SGWBスペクトル回復の精度と正確性に最も影響を与えるパルサーおよびデータ特性を特定すること。
• PTAデータからの天体物理学的・宇宙論的推論のための迅速で柔軟なモデル選択を可能にすること。
• 効率的な重力波背景のスペクトル解析のためのオープンソースソフトウェアスイートを提供すること。

提案手法

• 本手法は、タイミングパワー スペクトルの途中段階のベイジアン後部推定値を用い、それらをスペクトルモデルの再適合のための十分統計量として扱う。
• 最適化されたバンド幅（シータージョーンズ法を用いて）を用いたカーネル密度推定（KDE）を用いて、パワー スペクトルの後部分布を表現する。
• KDEから得られる後部確率を用いて、スペクトルパラメータ推定のための高速で近似尤度を構築する。
• パラメータ推定にはMCMC（PTMCMC）、モデル選択にはネストドサンプリング（UltraNest）をそれぞれ用いる。
• シミュレーテッドPTAデータおよびNANOGrav 12.5年分のデータセットを用いて手法を検証し、プロダクションレベルのパイプラインと比較する。
• ソフトウェアスイートceffylは、エンタープライズおよびPTArcadeと統合可能であり、ChainConsumerを介して効率的な後部分布比較を可能にする。

実験結果

リサーチクエスチョン

• RQ1 タイミングパワー スペクトルの事前計算済み後部推定値へのスペクトルモデルの再適合における、計算コストと精度のトレードオフは何か？
• RQ2 どのパルサーおよびデータ特性がSGWBスペクトル回復の精度と正確性に最も顕著に影響を与えるか？
• RQ3 迅速な再適合手法は、フルベイジアンパイプラインと同等の精度を達成しつつ、計算時間を桁違いに短縮できるか？
• RQ4 高周波数の重力波成分を含めると、共通プロセス信号の推定スペクトル形状にどのような影響が生じるか？
• RQ5 支配的である天体物理学的信号の下に、副次的な宇宙論的重力波背景を検出・特徴付けるための条件は何か？

主な発見

• NANOGrav 12.5年分のデータセットにおいて、迅速な再適合手法はフルプロダクションレベルパイプラインとヘリンジャー距離≤0.1を達成しており、後部分布がほぼ同一であることを示している。
• 本手法はフル尤度パイプライン比で10²〜10⁴倍高速であり、PTAサイズの増加に伴い非線形的にスケーリングするため、より大きなアレイへの効率的スケーリングが可能である。
• スペクトル特徴付けは、多数のパルサー間の相互相関ではなく、最も長時間測定されたパルサーおよび最も正確に測定されたパワー スペクトル密度に支配されている。
• NANOGrav 12.5年分のデータセットでは、45個のパルサーのうち約14個（主に最も長時間測定されたもの）の情報のみで、共通プロセスのスペクトル形状が完全に再現可能である。
• 重力波周波数が約50 nHzを超える成分を含めると、共通プロセスのパワーロー・モデルの「なだらかさ」が生じ、高周波数寄与によるスペクトルの曲がりが示唆される。
• オープンソースソフトウェアスイートceffylは、将来のPTA解析における迅速で柔軟かつ再現可能なスペクトル特徴付けとモデル選択を可能にする。