[論文レビュー] A fully precessing higher-mode surrogate model of effective-one-body waveforms

B. U. Gadre, M. Pürrer+3 more|arXiv (Cornell University)|Mar 1, 2022

Pulsars and Gravitational Waves Research被引用数 5

ひとこと要約

本稿では、完全にプレセッシングされた高モードを含むSEOBNRv4PHM効果的1体模型重力波形の、高速かつ高精度な代替モデルを提示する。この代替モデルは、パラメータ空間の大部分で1%未満の不一致誤差を達成し、重力波天文学におけるバイアス推定のための高速かつ実用的な大規模パラメータ推定を可能にする。

ABSTRACT

We present a surrogate model of exttt{SEOBNRv4PHM}, a fully precessing time-domain effective-one-body waveform model including subdominant modes. We follow an approach similar to that used to build recent numerical relativity surrogate models. Our surrogate is 5000M in duration, covers mass-ratios up to 1:20 and dimensionless spin magnitudes up to 0.8. Validating the surrogate against an independent test set we find that the bulk of the surrogate errors is less than $\sim 1\%$ in mismatch, which is similar to the modelling error of exttt{SEOBNRv4PHM} itself. At high total mass a few percent of configurations can exceed this threshold if they are highly precessing and they exceed a mass-ratio of 1:4. This surrogate is nearly two orders of magnitude faster than the underlying time-domain exttt{SEOBNRv4PHM} model and can be evaluated in $\sim 50$ ms. Bayesian inference analyses with exttt{SEOBNRv4PHM} are typically very computationally demanding and can take from weeks to months to complete. The two order of magnitude speedup attained by our surrogate model enables practical parameter estimation analyses with this waveform family. This is \emph{crucial} because Bayesian inference allows us to recover the masses and spins of binary black hole mergers given a model of the emitted gravitational waveform along with a description of the noise.

研究の動機と目的

• 完全にプレセッシングされたSEOBNRv4PHM波形モデルに、従属モードを含めた計算的に効率的な代替モデルを開発すること。
• バイアス推定における使用を制限する、元のSEOBNRv4PHMモデルの高い計算コストを克服すること。
• 質量比が1:20まで、スピンの大きさが0.8までの広範なパラメータ空間において、高い精度を維持すること。
• 評価時間を数週間から数分に短縮することで、大規模なパラメータ推定解析を実用的に行えるようにすること。
• 独立したテストセットとの比較を通じて代替モデルの正確性を検証し、現実のデータ解析シナリオにおける性能を評価すること。

提案手法

• 5000Mの時間領域におけるSEOBNRv4PHM波形のトレーニングセットを用いて代替モデルを構築し、質量比は1:20まで、スピンの大きさは0.8までをカバーする。
• パラメータ空間を質量比とスピンに基づいて部分領域に分割するドメイン分解アプローチを採用し、精度と収束性を向上させる。
• 波形のズレを補正するための動的時間ワープ（DTW）技術を用いた径路基底関数補間を実施し、フィッティング誤差を低減する。
• 数値相対論にインspiredされた技術と効果的1体模型物理学を組み合わせてトレーニングすることで、物理的一致性を保持する。
• 1000波形の独立したテストセットを用いて代替モデルを検証し、不一致誤差と位相誤差を計算して正確性を定量化する。
• 元のSEOBNRv4PHMモデルおよび最新のモデル（例：IMRPhenomXPHM）と比較して、代替モデルの性能をベンチマークする。

Figure 1: This figure shows the 14 subdomains in mass-ratio $(q)$ and effective spin parameter $(\chi_{\rm{eff}})$ of the SEOBNRv4PHM surrogate. Different colors correspond to the division along $q$ while shading types defines boundaries along $\chi_{\rm{eff}}$ .

実験結果

リサーチクエスチョン

• RQ1 完全にプレセッシングされた高モードを含むSEOBNRv4PHM波形に対して、高精度を維持しながら、計算速度を1000倍以上に向上させる代替モデルを構築可能か？
• RQ2 特に高プレセッシングおよび高質量比領域において、代替モデルの不一致誤差と位相誤差はパラメータ空間全体でどのように振る舞うか？
• RQ3 トレーニングドメインを越えて外挿する場合、精度の著しい低下を伴わずにどの程度まで外挿可能か？
• RQ4 IMRPhenomXPHMなどの既存の波形モデルと比較して、代替モデルの速度と正確性はどの程度か？
• RQ5 GW150914のような現実の重力波データ解析において、代替モデルはバイアスのない効率的なベイズパラメータ推定を可能にするか？

主な発見

• トレーニングドメイン全体で中央値の不一致誤差が1%未満であり、99%の構成で不一致誤差が1%未満である。
• 質量比が1:4を超える高プレセッシング系や高全質量系では、不一致誤差が1%を超える構成が最大3%程度にのぼるが、これは元のSEOBNRv4PHMモデルのモデル誤差範囲内に収まる。
• 代替モデルは元のSEOBNRv4PHMモデルと比較して約50倍高速であり、1波形あたりの評価時間は約50 msである。
• GW150914のパラメータ推定において、代替モデルはバイアスのない推定を可能とし、過去のLVC解析と一致するがわずかに狭い事後分布を得た。
• 80,000Mの長時間代替モデルも構築されたが、重複領域では類似した正確性を示したが、高時間分解能の要件のため、スケール上で計算が非現実的であることが判明した。
• トレーニングドメインを越えて外挿する（例：スピン>0.8）と、a=0.9では不一致誤差が最大5%、a=0.99では20%にまで上昇し、特に高質量比領域で顕著になる。これは外挿の限界を明確に示している。

Figure 2: Histograms of errors due to the errors in respective surrogate data pieces on the training set. Left: RMS errors due to dynamics quantities: spin vectors (solid violet and dash-dotted brown), unit quaternions (red dashed), and orbital phase (yellow long dash-dotted); Right: $L^{2}$ -type e