A parametric signal plus noise inference framework for short duration non-Gaussian noise transients

技術要約：短時間の非ガウス型ノイズ・トランジェントに対するパラメトリックな信号＋ノイズ推論フレームワーク

問題提起
Advanced LIGO、Virgo、KAGRAといった地上設置型の重力波（GW）検出器は、コンパクト連星合体（CBC）からの信号を日常的に検出している。これらのイベントに関する標準的なベイズ推論は、検出器ノイズが定常かつガウス分布に従うという基本的仮定に基づいている。しかし、地球上のノイズ・アーティファクトである「グリッチ」と呼ばれる現象が、重力波ストレイン・データに頻繁に混入する。これらの非ガウス型トランジェントが、時空間および周波数領域において天体物理学的信号と重なる場合、ガウスノイズの仮定は無効となる。これは、推論される天体特性（成分質量やスピンなど）に重大なバイアスをもたらし、一般相対性理論（GR）からの逸脱という誤った主張を招く可能性がある。BayesWaveやgwsubtractといった既存の緩和技術が存在するものの、それらは簡略化された信号モデルやデータ駆動型のノイズモデルを採用していることが多く、実際のグリッチの現象論を十分に捉えきれなかったり、意図せず天体物理学的な情報を除去してしまったりする場合がある。

手法
著者らは、標準的なベイズ推論ライブラリであるbilbyを拡張し、天体物理学的信号と非ガウス型ノイズ・トランジェントの同時推論を行うフレームワークであるbilby-antiglitchを導入する。

• 結合仮説： 本フレームワークは、「信号のみ」の仮説（$d = h + n$）の下でデータを解析するのではなく、「信号＋ノイズ」の仮説（$d = h + n + g$）を採用する。ここで、$g$はグリッチを表す。
• パラメトリック・グリッチモデル： 本フレームワークは、短時間の非ガウス型トランジェントを記述するために、準物理的かつパラメトリックなアプローチであるAntiGlitchモデルを利用する。データ駆動型モデル（ガウス過程やサイン・ガウス型ウェーブレットなど）とは異なり、AntiGlitchは周波数領域において、インターフェロメータごとに5つのパラメータ（ピーク周波数 $f_0$、逆帯域幅の二乗 $\gamma$、振幅 $A$、位相 $\phi$、中心時刻 $t_0$）によって定義される特定の関数形式を用いる。
• 実装： この手法は、AntiGlitchモデルをbilbyエコシステムに統合しており、最新の波形モデル（例：IMRPhenomXPHM）やネスト・サンプラー（例：dynesty）との連携を可能にしている。このモジュール設計により、将来的な代替グリッチモデルの導入も可能となっている。
• 尤度： グリッチを明示的にモデル化することで、残差データ（$d - M(\theta) - G(\beta)$）は標準的なガウス分布に従うことが期待され、これにより標準的なGW解析で使用されるウィトル（Whittle）尤度の妥当性が回復される。

主な結果
著者らは、理想的なガウスノイズの中に、強力に重なる「ブリップ（blip）」グリッチ（SNR $\rho=70$）が注入されたシミュレーションGW信号（GW150914に類似）を用いて、本フレームワークを検証した。

1. 真の天体特性の回収： 従来の信号のみのbilbyアプローチを用いた場合、グリッチの存在によって深刻なバイアスが生じ、重く非対称な連星であり、かつ極端に近い反並行スピンを持つと推論された。対照的に、bilby-antiglitchは、注入された真の天体特性（質量 $30 M_\odot, 28 M_\odot$ および整列スピン $\chi=0.2$）を高精度に回収することに成功した。
2. 信号およびグリッチの再構成： 本フレームワークは、天体物理学的信号とグリッチの両方を独立して正確に再構成した。推論されたSNRは、注入値（信号 $\rho \approx 30$、グリッチ $\rho \approx 70$）と一致した。
3. 一般相対性理論との整合性： 信号のみの解析では、最終質量とスピンの一致性テスト（$\Delta M_f / \bar{M}_f$ および $\Delta \chi_f / \bar{\chi}_f$）の事後分布が、GRの予測値 (0,0) を90%以上の信頼区間で除外しており、GRからの逸脱を示唆していた。一方、bilby-antiglitchによる解析は、GRとの整合性を回復し、期待される値を90%信頼区間内に回収した。
4. 計算効率： 高次元のパラメータ空間（信号 ＋ グリッチのパラメータ）をサンプリングしているにもかかわらず、bilby-antiglitchは、信号のみのbilby解析（940 CPU時間）と比較して、約半分の計算コスト（500 CPU時間）を要した。著者らは、信号のみのモデルが非ガウス型のグリッチを適合させようとして苦慮し、サンプリングが困難な複雑な尤度曲面を作り出していたためであると考えている。
5. 実データへの適用： 本フレームワークは、実イベントであるGW250114 082203およびGW200129 065458に適用された。GW250114 082203については、LVKの解析結果と一致しており、グリッチの有意な証拠は見られなかった。GW200129 065458については、本フレームワークにより、以前に報告された歳差運動の測定値が短時間の非ガウス型ノイズ・トランジェントに対して堅牢であることが確認された。

意義および主張
本論文は、bilby-antiglitchが、天体物理学的な情報を除去することなく、短時間の非ガウス型ノイズ・トランジェントの影響を緩和する堅牢な手法を提供すると主張している。パラメトリックなグリッチモデルを結合推論フレームワーク内に採用することで、本手法は以下を実現する：

• 強力なグリッチが存在する場合でも、天体パラメータの推論におけるバイアスを防ぐ。
• スプリアスなデータに起因する、一般相対性理論からの誤った逸脱の主張を回避する。
• データ駆動型のノイズモデルと比較して、解釈性が向上し、過学習のリスクが低減される。
• 信号のみの推論が非ガウス型アーティファクトのモデル化を強いられた場合、結合推論の方が計算効率が向上し得ることを示す。

著者らは、このアプローチが、スプリアスなデータにおけるGW信号の特異性を検証し、前処理としての除去手法に代わる実行可能な選択肢を提供し、正確なパラメータ推定のために基礎となる天体物理学的信号が保持されることを保証すると結論付けている。