Probing globular clusters parameters through gravitational wave lensing with stellar-mass black hole binaries

技術要約：恒星質量ブラックホール連星による重力波レンズを用いた球状星団パラメータの探索

問題提起
球状星団（GC）は、恒星力学および初期銀河形成を研究するための実験室として機能する高密度な恒星系である。それらの内部力学、形成、および進化を理解するための重要なパラメータは、通常、中心部の恒星速度分散（$\sigma_v$）から推測される質量である。現在の質量決定手法は、速度分散の測定や表面輝度プロファイルといった電磁波（EM）観測に依存している。しかし、これらの手法には、視線方向の投影効果、高密度な中心領域における恒星の混雑、明るい巨星によるショットノイズ、および真の星団中心の特定における不確実性といった、重大な限界が存在する。これらの課題は、コア崩壊を起こした星団や、中間質量ブラックホール（IMBH）を宿している可能性のある系においてさらに顕著になる。さらに、次世代の広視野サーベイ（例：LSST、Euclid）は、正確な位置と距離を持つGCカタログを拡大させるであろうが、速度分散のような内部運動学の測定は依然として困難であり、通常は専用の分光観測によるフォローアップを必要とする。

著者らは、補完的なプローブとして、GCによってレンズ効果を受けるコンパクト連星合体（CBC）からの重力波（GW）を利用することを提案している。強重力レンズが複数の像を生み出すのとは異なり、GCによるレンズ効果は波動光学領域（wave-optics regime）で作動する。この領域では、レンズは複数の信号を作成するのではなく、GW信号に対して周波数依存の増幅を刻み込む。本論文では、これらの波動光学的なシグネチャを用いて、レンズとなるGCの固有の特性、特にその中心速度分散を回収できるかどうかを調査している。

手法
本研究では、レンズ効果を受けたGW信号からGCのパラメータの回収可能性を検証するために、シミュレーションおよびベイズパラメータ推定フレームワークを用いている。

1. レンズモデリング： GCは単一等温球（SIS）としてモデル化されている。質量分布は、$\sigma_v$が一定であると仮定した場合、密度プロファイル $\rho(r) = \sigma_v^2 / (2\pi G r^2)$ で近似される。有効なレンズ質量（$M_{SIS}$）は、$\sigma_v$とレンズ距離（$D_l$）の関数として導出される。
2. 信号シミュレーション： 著者らは、レンズ効果を受けたGW150914（恒星質量ブラックホール連星の合体）に類似したGW信号をシミュレートしている。これらの信号は、フレネル・キルヒホッフ回折積分から導出されたSIS増幅率を用いてレンズ効果が付与されている。増幅率は、無次元周波数（$w$）およびインパクトパラメータ（$y$）に依存する。
3. インジェクション研究： レンズ効果を受けた信号は、LIGO Hanford、LIGO Livingston、Virgoの3検出器ネットワークのO4パワースペクトル密度（PSD）から生成されたガウスノイズの中に注入される。
4. パラメータ推定： 注入された信号は、ソースパラメータ（IMRPhenomXPHM近似を用いた15個のバイナリパラメータ）とレンズパラメータ（SISレンズ質量 $M_{SIS}$ およびインパクトパラメータ $y$）の両方を含む波形テンプレートを用いて、ベイズ推論（Bilbyパッケージおよびdynestyネスト・サンプリングを使用）を介して回収される。
5. 候補の特定： 回収された事後分布から、90%信頼区間の空の局在化領域が構築される。この領域は、既知の天の川銀河GCのカタログ（[37]より）と照合され、潜在的なレンズを特定するために使用される。
6. 速度分散の回収： 候補となるGCが特定された場合、カタログから得られる既知の距離（$D_l$）を用いて、$M_{SIS} \propto \sigma_v^4 D_l$ という関係式に基づき、$M_{SIS}$ の事後分布を反転させて$\sigma_v$を推定する。

本研究では、異なる中心速度分散を持つ、天の川銀河カタログから選ばれた3つの特定のGC（Liller-1: $\sigma_v \approx 19.6$ km s$^{-1}$、NGC 6093: $\sigma_v \approx 11.1$ km s$^{-1}$、NGC 2419: $\sigma_v \approx 5.1$ km s$^{-1}$）が選択された。各星団に対して、様々なインパクトパラメータ（$y \in [0.001, 5]$）を用いて20回のインジェクションが行われた。

主な結果

• レンズ質量と速度分散の回収： 有利なソースとレンズのアライメント（小さなインパクトパラメータ）において、注入された速度分散の値は68%信頼区間内で正常に回収された。Liller-1については、中央値として回収された値は19.7 km s$^{-1}$であり、注入された19.6 km s$^{-1}$と密接に一致した。
• インパクトパラメータの影響： 回収の精度は、インパクトパラメータ $y$ に強く依存する。$y$ が増加するにつれて増幅は減少し、信号へのレンズ効果の刻み込みが弱まる。これは、レンズ質量のサンプリングを悪化させ、推論される速度分散の不確実性を増大させる。$y$ が大きい場合、ノイズの実現態によって回収された中央値が真の値からずれることがある。
• 空の局在化とカタログ照合： 注入された星団は、テストされたすべてのケースにおいて、90%信頼区間の空の局在化領域内に正しく特定された。本手法は、ソースが特定のターゲット星団の背後に位置している限り、カタログ内の他の星団から対象を識別することに成功した。
• パラメータの縮退： 本研究では、$y \sim 1$ において、回収されたレンズ質量（$M_{SIS}$）とインパクトパラメータ（$y$）の間に反相関が観察された。これは、点質量レンズに関する先行研究の結果と一致している。

意義と主張
本論文は、レンズ効果を受けたGWが、GCのダイナミクス（特に速度分散および、ひいては質量の推定）のための補完的なプローブとして機能できることを示している。これは、高密度またはコア崩壊領域におけるEMベースの質量決定に固有の限界を克服するための新しい手法を提供する。

著者らは、天の川銀河のGCがGWイベントのレンズとして機能する確率は低い（$\tau \sim 10^{-8}$）と述べており、次世代の検出器がこのようなイベントを発見できる可能性を強調している。本手法は以下のためのツールとして提示されている：

1. 天の川銀河のGCがレンズである場合、GWイベントをレンズ効果を受けた信号として自信を持って特定すること。
2. 星団の速度分散および質量を制約すること。
3. 星団の中心におけるIMBHの存在の可能性を示すこと。

本研究の範囲は、SISモデルを使用し、距離がよく知られている天の川銀河のGCに焦点を当てているため、控えめなものである。著者らは、外部銀河のGC（カタログが不完全である可能性がある場合）においては、レンズ距離と速度分散の間に縮退が生じ、より広範なパラメータ推定が必要になることを認めている。また、本研究ではSISモデルを使用しているが、この手法は将来の研究において、より現実的なモデル（例：コアを持つ等温球、キングモデル）へと拡張可能であることも述べている。結論として、この手法が効果的であるために必要な空の局在化精度を高めるには、同等の感度を持つ少なくとも3つの検出器からなるグローバルなネットワークが不可欠であるとしている。