[論文レビュー] Non-parametric inference of the population of compact binaries from gravitational wave observations using binned Gaussian processes

Anarya Ray, Ignacio Magaña Hernandez+3 more|arXiv (Cornell University)|Apr 17, 2023

Pulsars and Gravitational Waves Research参考文献 71被引用数 4

ひとこと要約

本論文は、重力波データからのコンpakト二重星集団の非パラメトリック推論のためのボーリングされたガウス過程フレームワークを導入し、同時に成分質量分布と赤方偏移分布をモデル化しながら、それらの間の相関を捉える。この手法により、特定の関数形を仮定しない最小限の前提で、質量ピークや赤方偏移の進化といった特徴を強固に検出可能であり、シミュレートされたデータおよび実際の GWTC-3 データで検証された結果、ボーリング解像度に関わらず一貫性があり正確な再構成が得られた。

ABSTRACT

The observation of gravitational waves from multiple compact binary coalescences by the LIGO-Virgo-KAGRA detector networks has enabled us to infer the underlying distribution of compact binaries across a wide range of masses, spins, and redshifts. In light of the new features found in the mass spectrum of binary black holes and the uncertainty regarding binary formation models, non-parametric population inference has become increasingly popular. In this work, we develop a data-driven clustering framework that can identify features in the component mass distribution of compact binaries simultaneously with those in the corresponding redshift distribution, from gravitational wave data in the presence of significant measurement uncertainties, while making very few assumptions on the functional form of these distributions. Our generalized model is capable of inferring correlations among various population properties such as the redshift evolution of the shape of the mass distribution itself, in contrast to most existing non-parametric inference schemes. We test our model on simulated data and demonstrate the accuracy with which it can re-construct the underlying distributions of component masses and redshifts. We also re-analyze public LIGO-Virgo-KAGRA data from events in GWTC-3 using our model and compare our results with those from some alternative parametric and non-parametric population inference approaches. Finally, we investigate the potential presence of correlations between mass and redshift in the population of binary black holes in GWTC-3 (those observed by the LIGO-Virgo-KAGRA detector network in their first 3 observing runs), without making any assumptions about the specific nature of these correlations.

研究の動機と目的

• 重力波観測からコンパクト二重星の質量と赤方偏移の集団分布を推論する非パラメトリックでデータ駆動型の手法を開発すること。
• 特定の関数形を仮定せずに、質量分布と赤方偏移分布を同時にモデル化し、それらの間の相関を捉えること。
• 潜在的な集団分布に存在する未モデル化された特徴を逃がす可能性のあるパラメトリック仮定への依存を低減すること。
• 顕著な測定不確実性とモンテカルロサンプルノイズが存在する状況でも正確な推論を可能にすること。
• シミュレートされたデータ上で手法を検証し、公開済みの GWTC-3 データを再分析することで、従来のパラメトリックおよび非パラメトリック手法と比較すること。

提案手法

• 成分質量と赤方偏移の連合分布を、対数質量および赤方偏移空間におけるボックスに分けたガウス過程モデルを用いて表現する。
• 重力波データにおける測定不確実性と検出感度を反映するために、非ガウス型の尤度関数を用いる。
• ガウス過程のハイパーパramータ（滑らかさを制御する長さスケールを含む）を推論するために階層ベイズフレームワークを適用する。
• 計算的に効率的な方法でモンテカルロ不確実性を周辺化し、尤度のカットを避け、耐障害性を確保する。
• 宇宙論的体積と赤方偏移の進化に整合する人口レート密度に対する事前分布を用いる。
• 既知の潜在的分布を持つシミュレートされたデータ上で手法をテストし、代替の推論手法と比較することで検証を行う。

実験結果

リサーチクエスチョン

• RQ1 特定の関数形を仮定せずに、非パラメトリック手法がコンパクト二重星の成分質量分布に特徴を検出できるか？
• RQ2 この手法は、合体率の赤方偏移進化とその質量分布との相関をどの程度正確に再構成できるか？
• RQ3 ボーリング解像度を増加させた場合、モデルが安定し、推論が収束しているか？
• RQ4 質量と赤方偏移の相関が集団推論に与える影響は何か？また、事前の仮定なしにそれらを検出できるか？
• RQ5 実際の GWTC-3 データにおいて、この手法は従来のパラメトリックおよび非パラメトリック集団推論手法と比較して、どの程度正確かつ一貫性があるか？

主な発見

• ボーリングされたガウス過程モデルは、顕著な測定不確実性が存在するシミュレートされたデータにおいても、潜在的質量分布および赤方偏移分布を高い精度で再構成できた。
• ボーリング解像度を2倍にした場合でも、GPの長さスケールの事後分布が安定しており、収束性と耐障害性を示している。
• 異なるボーリング選択に対して一貫した集団制約が得られたことから、特徴検出を越える解像度が推論を変えることはないことが確認された。
• GWTC-3 データの再分析では、ブラックホール二重星集団において質量と赤方偏移の強い相関の証拠は得られず、モデルの感度範囲内では有意に検出されなかった。
• 計算的ペナルティなしにモンテカルロ不確実性を正しく周辺化できており、一部の既存手法とは異なり、耐障害性を有していた。
• このフレームワークは、パラメトリック仮定なしに質量スペクトルおよび赤方偏移進化の部分構造を検出可能であり、天体物理学的特徴のデータ駆動型発見を支援する。