[論文レビュー] Accelerated parameter estimation in Bilby with relative binning

K. Murali Krishna, A. Vijaykumar+6 more|arXiv (Cornell University)|Dec 10, 2023

Pulsars and Gravitational Waves Research被引用数 7

ひとこと要約

この論文では、オープンソースのベイズ推論ソフトウェアBilbyに相対的ビニングの実装を提示し、重力波パラメータ推定の高速化を図っている。周波数帯域ごとに波形比を区分的線形関数で近似することで、尤度評価のコストを低減し、1つのCPUコアでGW170817のパラメータ回復を14時間で実現しながらも、精度を維持し、統計的検定に合格している。

ABSTRACT

We describe an implementation of the relative binning technique to speed up parameter estimation of gravitational-wave signals. We first give a pedagogical overview of relative binning, discussing also the expressions for the likelihood marginalized over phase and distance. Then, we describe the details of the code in exttt{Bilby}, an open-source software package commonly used for parameter estimation of gravitational-wave sources. Our code is able to reproduce the parameters of GW170817 in 14 hours on a single-core CPU, performs well on simulated signals, and passes the percentile-percentile (p-p) tests. We also illustrate that relative binning is an ideal technique to estimate the parameters of signals in next-generation gravitational wave detectors.

研究の動機と目的

• 長期間の重力波信号、特に二進 neutron 星系からのベイズ的パラメータ推定の高い計算コストに対処すること。
• 精度を損なわず、尤度評価に要する時間を短縮すること。
• 次世代の重力波観測所における長期間の帯域内信号と高感度を想定した、より迅速な推定を可能にすること。
• Bilbyソフトウェアパッケージに相対的ビニングを実装し、効率的でスケーラブルかつ再現可能なパラメータ推定を可能にすること。
• 後続の統計的検証手法（例：百分位数-百分位数（p-p）検定）との互換性を確保し、信頼性を確認すること。

提案手法

• 周波数帯域ごとに、パrameter空間における重力波形比を区分的線形関数で近似することで、必要な波形評価回数を削減する。
• 帯域端での波形比の基本係数を事前に計算し、サンプリング中に波形を再評価せずに高速な補間を可能にする。
• 補間された波形比と事前に計算された要約統計量にのみ依存する修正された尤度定式化を採用し、リアルタイム計算を最小限に抑える。
• Bilbyに統合し、既存のベイズ推論とサンプリングのインfraを活用する。
• 時間マージナライゼーションの2通りの代替手法を実装：1つは波形比の補間、もう1つは修正された要約統計量。ただし、後者は計算的に非効率であることが示された。
• 実データおよびシミュレートされた信号にこの手法を適用し、p-p検定を用いて結果を検証し、標準的な尤度手法と比較した。

実験結果

リサーチクエスチョン

• RQ1 相対的ビニングは、精度を損なわず、重力波パラメータ推定の計算コストを著しく低減できるか？
• RQ2 GW170817のような実際の重力波イベントにおいて、相対的ビニングの実行時間とパラメータ回復性能はどのように評価できるか？
• RQ3 相対的ビニングを用いた場合、時間マージナライゼーションは実現可能で有益か、それとも計算上のオーバーヘッドを引き起こすか？
• RQ4 相対的ビニングはBilbyソフトウェアパッケージに効果的に統合可能で、広くコミュニティが利用できるか？
• RQ5 特に次世代の観測所（例：Cosmic Explorer や Einstein Telescope）における信号長の延長と感度向上を想定した場合、この手法はどのようにスケーリングするか？

主な発見

• Bilbyにおける相対的ビニングの実装により、1つのCPUコアでGW170817のパラメータをわずか14時間で回復可能となり、標準的手法に比べて顕著な高速化が達成された。
• この手法は百分位数-百分位数（p-p）検定に合格し、事後分布が統計的に妥当で適切にキャリブレーションされていることが確認された。
• シミュレートされた信号に対しても精度が維持されており、多様な信号パラメータとノイズ実現の下でも頑健であることが示された。
• 波形比の補間による時間マージナライゼーションは実現可能ではあるが、標準的な相対的ビニングよりも計算コストが高いため、全体としての効率性は劣る。
• 信号長の延長に伴い、この手法は良好にスケーリングされ、信号長が延びる帯域内信号を観測する次世代の観測所（例：Cosmic Explorer や Einstein Telescope）に適している。
• この実装はBilbyに統合されており、最小限の実行時オーバーヘッドで再現可能でオープンソースのパラメータ推定を可能としている。