[論文レビュー] Rapid inference of gravitational-wave signals in the time domain using a heterodyned likelihood

Neha Sharma, A. Vijaykumar+1 more|arXiv (Cornell University)|Jan 16, 2026

Pulsars and Gravitational Waves Research被引用数 0

ひとこと要約

時域での異周波数性尤度を導入し、波形モードをダウンサンプリングして事前計算済みの概要データを使用することで、長い重力波信号のベイズ推定を劇的に高速化。後方分布は劣化せず、多量の処理量を削減する。

ABSTRACT

Parameter estimation of gravitational wave signals is computationally intensive and typically requires millions of likelihood evaluations to construct posterior probability distributions. This computational cost increases significantly in the time domain, which requires non-diagonal covariance matrices to compute the likelihood. Consequently, parameter estimation of long-duration gravitational wave signals, such as binary neutron star mergers, becomes computationally infeasible in time domain. In this work, we detail a framework for the heterodyned likelihood that enables rapid inference in the time domain. Our method is applicable to signals with arbitrary mode content, and leverages the smoothness of the ratio of complex-valued waveform modes, approximating the ratio as a linear function within appropriately chosen time bins. This allows downsampling of the waveform modes and a reformulation of the likelihood, such that it depends only on the bin edges. We demonstrate that this likelihood recovers posteriors that are indistinguishable from those obtained using the standard likelihood in the time domain. We also observe dramatic improvement in speed - for a 128 seconds-long gravitational wave signal, our method is at least $\sim 400$ times faster than the standard time-domain analysis, reducing the wall clock time to just a few hours. We also demonstrate the reliability and unbiasedness of the likelihood using percentile-percentile tests for binary black hole and binary neutron star injections. We use the Gohberg-Semencul representation of the inverse of Toeplitz covariance matrix to accelerate matrix-vector products, which has potential applications even in non heterodyned time-domain inference.

研究の動機と目的

• 時間領域の GW 分析における高い計算コストに対処するための迅速なパラメータ推定の必要性を動機づける。
• 時間ビン内の波形モード比の線形近似を用いる異周波数性尤度フレームワークを開発する。
• 尤度評価を加速しつつ精度を保持する概要データとビニング戦略を定式化する。
• BBHおよびBNSシナリオにおけるインジェクションとP-Pテストを通じて精度と速度向上を実証する。

提案手法

• 尤度を複素波形モンペル h^{3lm} と基準波形 h^{3lm}_o の比の形として表現し、時ビン内でこれらの比を線形近似（異周波数化）する。
• ビン端でのみ波形モードをダウンサンプリングし、必要時に補間して完全な波形を再構成する。
• 採取中の大規模な行列積を避けるため、基準モード、検出機器ノイズ、およびデータに依存する概要データを事前計算する。
• 対数尤度を比 r^{3lm} および概要データを用いて表現し、反復ごとの行列逆行列計算と積を回避する。
• 吸い込み段階での適応的な時ビニング方式を導入し、内挿増減で近接部の精度と速度のバランスをとる。

実験結果

リサーチクエスチョン

• RQ1 異周波数時系列尤度は、全時系列尤度と同一の後方分布を再現できるか。
• RQ2 時間領域の異周波数化を用いた場合、信号長が異なる GW 信号でどの程度の速度アップが得られるか。
• RQ3 特に連結時の merger 近傍で、計算を最小化しつつ精度を維持するにはビンをどう選ぶべきか。
• RQ4 サブ支配モードおよび事前回転（プリセッション）波形への拡張性はあるか。
• RQ5 BBH および BNS インジェクションでのパーセンタイル・パーセンタイル検定は、方法が有偏でないことをどう示すか。

主な発見

Duration of GW signal [s]	Number of bins	Summary data computation time [s]	Heterodyned Likelihood [ms]	Full Likelihood [ms]	Speedup factors	Number of bins (Frequency Domain)	Heterodyned Likelihood [ms]	Full Likelihood [ms]	(Frequency domain speedup)
2	191	8.6989	2.290	42.9296	120	119	1.445	4.596
4	193	53.9353	2.418	89.97518	121	121	1.724	8.662
8	302	116.1305	3.860	216.558	121	121	1.812	14.14
16	382	561.3549	6.378	791.557	124.10740	121	1.801	27.35
32	476	967.96419	8.580	884.42	103.07940	121	1.848	51.18
64	483	4103.1448	9.140	4163.54	455.5295	121	1.817	105.3
128	486	20012.4119	9.191	3129.884	340.5379	121	1.849	209.6

• 異周波数時域尤度からの後方分布は、試験的なインジェクションで全尤度の後方分布と区別できない。
• 2秒信号では、完全な時域法と比較して尤度評価が約19倍高速。
• 16秒および128秒の信号では、単一尤度評価が約6.4 msおよび約9.2 msとなり、それぞれ約124倍および約340倍の速度アップ。
• 信号長に応じて速度アップが拡大し、週から月単位の解析を数時間以内に完了可能。
• 適応的で吸い込みを優先し merger付近で絞るビニング戦略により、精度を保ちながら波形評価回数を劇的に削減。
• BBHおよびBNSインジェクションに対するP-P検査は後方分布の有偏性がないことを示し、信頼性を検証。