Forecasting Strongly Lensed Binary Black Hole Rates from the Stochastic Gravitational-Wave Background with TianQin

技術要約：TianQinによる確率的重力波背景放射からの強重力レンズ重力波連星ブラックホール率の予測

問題提起
ミリヘルツ帯の確率的重力波背景放射（SGWB）は、恒星質量連星ブラックホール（BBH）の集団レベルのプローブおよびその宇宙論的な合体履歴を提供します。しかし、SGWBスペクトルはソースの質量や赤方偏移を積分してしまうため、基礎となる集団パラメータ間の縮退が生じやすくなります。TianQinのような計画中の宇宙配置型検出器は、この背景放射を測定できる可能性がありますが、重大な課題が残っています。それは、SGWBの推論に内在する不確mathcalty（不確かさ）を、どのように地上観測量、具体的には強重力レンズ効果を受けたBBH合体へと伝播させるかという点です。先行研究では、未分解の背景放射とレンズ効果を受けた合体率を結びつけていますが、宇宙配置型SGWBの事後分布（検出器ノイズやフォアグラウンドを含む）から、地上検出器のための特定の画像レベルおよびシステムレベルのレンズ観測量へと繋ぐ厳密なフレームワークは欠けていました。

手法
著者は、単一のTianQinコンステレーションを用いて、SGWBから得られたBBH集団の制約を、強重力レンズ合体の予測へと伝播させる定量的なフレームワークを開発しました。この手法は、以下の連結された3つのステップで構成されます。

1. TianQinによるSGWB推論:

   • 検出器モデル: 本研究では、時間遅延干渉法（TDI）を用いて直交する$A, E, T$チャネルを構築し、地球中心のコンステレーションとしてのTianQinをモデル化しています。$A$および$E$チャネルは重力波に対して感度を保持しますが、$T$チャネルは低周波数帯において感度が低くなります。
   • 尤度構成: $A/E$チャネルの$2\times2$複素スペクトル共分散に基づく、位相折り畳み（phase-folded）複素ウィシャート尤度を採用しています。このモデルには以下が含まれます：
     • 計器ノイズ（レーザー周波数ノイズおよび変位ノイズ）。
     • 反応投影された空の異方性テンプレートでモデル化された、未分解の二重白色矮星（DWD）による銀河系フォアグラウンド。
     • 振幅とスペクトルの傾きによってパラメータ化されたBBH SGWB信号。
   • 模擬データ: 1年間の模擬観測が生成されており、特定のデューティサイクル（2つの90日間ウィンドウ）と、周波数帯域 $0.10175$mHz から$0.05$ Hz を備えています。解析には、銀河系フォアグラウンドの周期的な変調を捉えるために32個の軌道位相ブロックが使用されます。

2. 集団の伝播:

   • BBH集団は、GWTC-5.0のPixelPopモデルを用い、局所合体率 $R_0$ にスケーリングして特徴付けられます。
   • 赤方偏移の進化は、ターンオーバーを持つ現象論的なモデルで記述され、高赤方偏移パラメータ（$\alpha_z, z_p$）を固定した上で、低赤方偏移の成長パラメータ（$\lambda_z$）をサンプリングします。
   • 極めて重要な点として、スペクトル振幅 $A_{BBH}$ は独立してサンプリングされるのではなく、集団積分から導出されます。これにより、スペクトルの形状と物理的な合体率との間の整合性が保証されます。

3. レンズおよび選択の予測:

   • レンズモデル: 銀河スケールの単一等温球（SIS）レンズ集団が、速度分散関数によって定義されて使用されます。
   • 選択ルール: 本フレームワークは、2種類の観測量を区別します：
     • イベント率 ($R^{det}_{event}$): 少なくとも1つの像が検出されたレンズ天体の発生率。
     • システム率 ($R^{det}_{sys}$): 両方の像が検出され、かつ（時間窓 $T_{assoc}$ 内で）正常に関連付けられたレンズ系の発生率。
   • 検出基準: 画像の検出可能性は、最適信号対雑音比（SNR）のスケーリングとMarcum-Q関数を用いて計算され、等方的な空の方向および方位について平均化されます。システム増幅率の閾値（$\mu_{sys} \geq 2$）が適用されます。

主な貢献

• エンドツーエンドのフレームワーク: 本論文は、完全な検出器レベルのモデル（ノイズとフォアグラウンドを含む）を用いたミリヘルツ帯SGWBの推論と、地上観測装置によって検出可能な強重力レンズBBHの予測を繋ぐ、最初の定量的なリンクを確立しました。
• 観測量の区別: 「レンズイベント」（単一画像の検出）と「レンズシステム」（二重画像の関連付け）を明示的に区別し、それぞれに対して個別の発生率予測を提供しています。
• 事後分布の伝播: 固定されたSGWB振幅を仮定するのではなく、集団パラメータ（赤方偏移の進化や質量モデルを含む）の完全な事後分布をレンズ計算へと伝播させます。これにより、SGWB推論における不確かさがレンズ効果の発生率にどのように変換されるかを捉えています。
• 感度分析: 周波数サポート、位相情報、クロススペクトルの使用、およびレンズモデルのパラメータに対する予測の感度を厳密にテストしています。

結果
1年間の模擬観測とGWTC-5.0 PixelPop集団モデルを用いた結果、システム増幅率の閾値 $\mu_{sys} \geq 2$ に対して、以下の中央値予測率が報告されています：

• レンズイベント率: $R^{det}_{event} = 5.50 \times 10^{-4} \, \text{yr}^{-1}$
• レンズシステム率: $R^{det}_{sys} = 2.54 \times 10^{-4} \, \text{yr}^{-1}$
• 検出可能システム分率: $F^{det}_{sys} = 3.35 \times 10^{-4}$

本研究は、SGWB測定による情報利得を定量化しています。事前分布のみの伝播（背景振幅が事前分布の平均値に固定されている場合）と比較して、結合された16〜84パーセンタイル幅は以下のように減少します：

• イベント率に対して 6.5%
• システム率に対して 6.4%
• システム分率に対して 4.0%

解析により、SGWB測定が、パラメータ間の相関により赤方偏移の履歴は部分的に特定されるものの、集団に起因するピボット振幅とスペクトルの傾きを制約することが確認されました。

意義と主張
本論文は、宇宙配置型のSGWB測定と地上のレンズ探索を結びつける将来の解析のための「定量的なベースライン」を提供することを目的としています。著者は、報告された結果が以下の特定のモデルに条件付けられていることを強調しています：

• GWTC-5.0 PixelPop質量投影。
• 赤方偏移進化における固定されたターンオーバーパラメータ。
• SISレンズモデルおよび速度分散関数。
• 画像探索および関連付けのルール（例：$T_{assoc} = 365.25$ 日）。
• DWD空間テンプレートおよびフォアグラウンド較正（主に0.5–5 mHz帯で有効）。

著者は、報告された不確かさの減少は、モデルに依存しない測定精度を表すものではないことを明示しています。むしろ、本研究は方法論的な経路を示しているのです。すなわち、位相分解された尤度は、画像、イベント、およびシステムのカウントの区別を維持しながら、物理的に定義された検出可能なレンズイベントへと接続できるということです。このフレームワークにより、データが入手可能になるにつれて、将来の研究が固定された構成要素を、経験的に較正された集団および関連付けモデルへと置き換えることが可能になります。