Prospects of resolving and localising individual supermassive black hole binaries with pulsar timing arrays: the host ranking challenge

技術要約：パルサー・タイミング・アレイによる個別の超大質量ブラックホール連星の解明と局在化の見通し

問題提起
パルサー・タイミング・アレイ（PTA）は、近年、超大質量ブラックホール（SMBH）の集団に由来すると考えられる確率的重力波（GW）背景放射を検出した。次の大きな節目は、この背景放射の中から個別の連星を解明することであり、それが実現すればマルチメッセンジャー天文学が可能となる。しかし、主要な課題は、PTAの空の局在化能力が限られていることである。最初に個別に解明される連星の90%信頼区間局在領域（$A_{90}$）は、数百平方度に及ぶことが予想され、そこには数十万個の候補銀河が含まれる可能性がある。さらに、既存の電磁波（EM）銀河カタログは、恒星や塵の混入が多い領域（例：銀河面）において不完全であり、「欠落した」ホスト候補が生じている。膨大な数の候補をGWの後験分布およびEM特性に基づいて効率的にランク付けし、追跡観測の優先順位を決定するための標準化された手法は、現在存在しない。

手法
著者らは、holodeckソフトウェアを用いて、NANOGrav 15年データセット（NG15）で測定されたGW背景放射と整合する、現実的なSMBH連星集団をシミュレートしている。彼らは1,000個の連星集団の実現例を生成し、将来のPTA構成では検出可能であるがNG15では検出不可能である、各集団の中で最も明るい連星を選択した。

将来の感度を代表するものとして、3つのPTA構成をシミュレートしている：

• IPTA_20: 116個のパルサーを用いた20年のベースライン（次期IPTA第3データリリースに類似）。
• IPTA_25 (基準モデル): 158個のパルサーを用いた25年のベースライン（IPTA_20に42個の新しいパルサーを追加）。
• IPTA_30: 200個のパルサーを用いた30年のベースライン。

著者らは、選択された連星を、確率的GW背景を含むシミュレーションされたタイミングデータに注入し、QuickCWパイプラインを使用して、連星パラメータ（空の位置、チャープ質量、振幅、周波数）を制約した。彼らは、「局在化可能な」波源を、位置が十分に制約されるために通常必要とされる閾値である信号対雑音比（S/N）$\ge 8$を持つものと定義した。

ホストを特定するために、著者らは得られた$A_{90}$局在マップを、2つの包括的な全天EMカタログと照合した：

1. ETGカタログ: 2MPZおよびWISExSCOSから導出された光度赤方偏移と恒星質量を持つ、約1,940万個の静穏な早期型銀河（ETG）を含む。
2. AGNカタログ: W1バンドのW1帯域見かけ等級を持つ、約450万個の活動銀河核（AGN）候補（AssefらによるR90カタログ）を含む。

ホスト・ランキング・システムが開発され、各候補銀河にスコア（$H$）を割り当てる。このスコアは、天球位置の確率（$p_{sky}$）と、銀河のEM特性がGWの後験分布と一致する確率（$p_{cat}$）の積である。

• ETGの場合、$p_{cat}$は銀河の赤方偏移（光度距離に変換）および、恒星質量からスケーリング関係を通じて導出された推定SMBH質量から導かれる。
• AGNの場合、$p_{cat}$は、連星の全質量と光度距離の関数である、見かけのW1等級から導かれる。
  著者らは、EMカタログのカバー範囲外となっている$A_{90}$領域の割合を推定することにより、「欠落した」ホストの数を明示的に計算している。

主な結果

• 検出可能性と局在化:

  • 次の0〜5年以内（IPTA_20）において、連星を解明できる確率（$S/N \ge 4$）は約21%であるが、十分に制約された局在化を達成できる確率（$S/N \ge 8$）は約0.3%である。
  • 次の5〜10年以内（IPTA_25）において、連星を解明できる確率は約38%に上昇し、局在化できる確率は約3.8%となる。
  • 次の10年以内（IPTA_30）において、解明の確率は約51%に達し、局在化できる確率は約14.1%となる。
  • 局在化可能な波源は、比較的近傍（中央値 $d_{lum} \approx 400$ Mpc）にあり、かつ非常に大質量（$M_{tot} \gtrsim 10^{9.8} M_{\odot}$）であることが予想される。

• 局在領域とホスト数:

  • 基準構成（IPTA_25）において、潜在的に局在化可能な38個の連星のうち30個は、$A_{90} < 2,000$ deg$^2$である。中央値の$A_{90}$は408 deg$^2$である。
  • これらの領域には、中央値で約19万個のETGと約4万個のAGNが含まれている。
  • IPTA_30では、中央値の$A_{90}$は212 deg$^2$に縮小し、中央値のホスト数は約40%減少する。

• 不完全なカバー範囲:

  • 不完全なEMカバー範囲は、重大な制限要因である。基準構成において、欠落しているETGホストの中央値は約2万個、欠落しているAGNホストの中央値は約6,000個である。
  • 一部の波源については、$A_{90}$がカタログのフットプリント（例：銀河面との重複）の外側に部分的に、あるいは完全に位置しており、新しいサーベイなしでは真のホストを特定できない。

• ランキングの効率性:

  • ランキング手法は、特性がGW後験分布と矛盾する銀河を効果的に除外することに成功した。
  • ETG: 赤方偏移またはSMBH質量の不一致により、約47%（IPTA_25）から61%（IPCTA_30）のETGが除外される。ランキング効率（$E_{95-50}$、95パーセンタイル・スコアと中央値の比）は、IPTA_25で約12.4、IPTA_30で約17.2である。
  • AGN: 約0.3%（IPTA_25）から3.8%（IPTA_30）のAGNのみが除外される。ランキング効率は著しく低く（IPTA_25で約3.7）、これは見かけの等級が見かけ上の質量と距離の縮退した組み合わせであること、およびチャープ質量がPTAによって十分に制約されていないことに起因する。

意義と主張
本論文は、最初のPTA解明SMBH連星に関する、候補ホストの数と不完全な空のカバー範囲の影響を定量化した最初の現実的な研究であると主張している。著者らは、ランキングシステム（GW後験分布とEM特性に基づく）を用いることで、特に赤方偏移と恒星質量が利用可能なETGにおいて、探索空間を効果的に縮小できることを示している。

本研究は以下の点を強調している：

1. マルチメッセンジャー発見は可能だが困難である： 10年以内に波源を局在化できる確率は無視できない（約14%）が、局在領域は依然として大きい。
2. カタログの完全性が極めて重要である： サーベイのギャップ（例：銀河面）による数千の欠落したホストの存在は、ホスト特定における偽陰性のリスクをもたらす。
3. データの質がランキングの成功を左右する： ランキングの効率は、独立した赤方偏移および質量推定の利用可能性に強く依存している。現在のAGNのランキングが効率的に行えない理由は、AGNカタログにおける赤方偏移や質量のデータの欠如、およびPTAによるチャープ質量の制約の弱さに帰せられる。
4. 将来の改善： 著者らは、パルサーの距離測定の向上（GW信号のパルサー項を利用するため）および、赤方偏移と質量推定を備えたAGNカタログの開発が、局在化およびランキング能力を大幅に向上させると示唆している。

本論文は、最初のマルチメッセンジャー検出は射程圏内にあるものの、その成功はPTAの感度だけでなく、EMカタログの完全性と一様性、および追跡観測の優先順位を決定するための効率的なランキングパイプラインの開発にかかっていると結論付けている。