The multimessenger view of Pulsar Timing Array black holes with the Horizon-AGN simulation

問題提起
超大質量ブラックホール（MBH）の起源と進化、特に局所宇宙における合体については依然として不透明である。パルサー・タイミング・アレイ（PTA）は、最近、大質量ブラックホール連星（MBHB）の確率論的な重なり合いと一致する重力波背景放射（GWB）信号を検出したが、観測された振幅と、より低い振幅を支持することが多い理論的予測との間には、持続的な乖ense（テンション）が存在する。さらに、観測データの有限な期間や、背景放射の上に立つ個別に分解可能な連続波（CW）源の存在によって、GWBスペクトルの解釈は複雑になっている。理論モデルをPTAデータと比較する現在の手法は、バイアスを導入する可能性のある簡略化された仮定に依存していることが多い。本論文は、Horizon-AGN宇宙論的シミュレーションを用いて、PTA帯域に寄与するMBHBの特性を研究し、GWBおよびCW予測のための現実的な集団を生成し、タイミング残差からGWBスペクトル特性を推論するために用いられる手法を批判的に評価することを通じて、これらの問題に対処することを目的とする。

方法論
本研究では、標準的な$\Lambda$CDM宇宙論を用いた$142 \, \text{cMpc}^3$の体積をモデリングする、流体力学的なHorizon-AGN宇宙論的シミュレーションを利用している。著者らは、シミュレーションを通じてMBH IDを追跡し、数値的な合体を特定することで、MBH合体のカタログを構築した。数値的な合体（分離距離 $\sim 4 \, \text{kpc}$）と実際の合体との間の物理的な遅延を考慮するため、以下のポストプロセッシング・パイプラインを適用した：

1. 力学的摩擦： 等温球近似内での恒星の力学的摩擦に基づいて計算。
2. 連星の硬化： 確立された処方箋（例：Sesana & Khan 2015; Dotti et al. 2015）に従い、恒星による硬化、円周連星円盤内の粘性トルク、および重力波（GW）放射をモデル化。
3. 離心率の進化： ロスコレ散乱、粘性抵抗、およびGW放射を介して追跡し、先行研究の結果に基づき初期離心率を割り当て。

GWBを推定するために、著者らは2つの異なるアプローチを採用した：

• ガウスアンサンブル（GE）アプローチ： 連続的なソース分布を仮定した特性歪みスペクトルの決定論的な計算であり、文献でよく用いられる手法である。
• 集団アプローチ： 合体カタログに基づくポアソン分布から、インスパイラリング（渦巻状接近）中の連星を抽出することにより、2,000個の「宇宙の実現例」を生成する確率論的手法。これは宇宙論的ゆらぎとソースの離散性を考慮している。

著者らは、これらの集団を用いて、2つの方法でPTAタイミング残差をシミュレートした：

1. ガウスアンサンブル： 平均化されたパルサー応答を介して、特性歪みをタイミング残差のパワースペクトル密度（PSD）に変換。
2. 集団アプローチ： すべての連星の地球項（Earth-term）GW歪みの寄与を合計し、有限の観測期間（$T_{\text{obs}}$）と、sinc関数によるスペクトルリーク効果を明示的にモデル化。

これら理論的予測を、現在のEPTA（25個のパルサー、10.33年）および将来のSKA（80個のパルサー、10.33年）の感度を表す「トイ（玩具）」のPTAデータセットと比較した。また、潜在的なCW候補（信号対雑音比（SNR） $> 3$ の個別の連星）を特定し、AGN放射およびホスト銀河の特性をモデル化することで、それらの電磁気学的特性（光学的、近赤外線、電波、X線）を分析した。

主要な貢献と結果

• GWBスペクトル推論とバイアス： 本研究は、理論モデルをPTAデータと比較する際の重大な方法論的注意点を強調している。無限の観測期間を仮定する標準的なガウスアンサンブル・アプローチは、有限の $T_{\text{obs}}$ によるスペクトルリークを考慮できていない。現実的なPTA解析（タイミングモデルの周辺化を含む）を用いて集団アプローチを用いた場合、推論されるGWBスペクトルは、真の基礎となるスペクトルと比較して、傾きが緩やかになり、$f = 1/\text{yr}$ において振幅が高くなるバイアスが生じる。低周波数ソースからのスペクトルリークによって引き起こされるこのバイアスは、理論的予測と観測されたPTA振幅との間のテンションを部分的に説明できる可能性がある。
• 確率性とソース数： 低周波数（$\sim 3 \, \text{nHz}$）において、GWBは数百から数千の連星によって支配されており、確率論的な背景放射近似の妥当性を裏付けている。しかし、高周波数（$\sim 10\text{--}50 \, \text{nHz}$）では、寄与するソースの数が大幅に減少し（数十個程度）、確率論的な性質の喪失と宇宙論的ゆらぎの増大を招く。離心率を持つ集団は、高次高調波での放射により、円軌道の集団よりも高い周波数まで確率論的な性質を維持する。
• ソースの特性（背景 vs. CW）：
  • 背景： GWBは、赤方偏移 $z \in 0.05 - 1$ におけるチャープ質量 $M_c \sim 10^{8.5} - 10^{9.5} \, M_\odot$ の連星によって支配されている。これらは主に、質量 $\sim 10^{13} \, M_\odot$ のハロー内に位置する、大質量銀河（$M_\star > 10^{11} \, M_\odot$）に宿っている。
  • CW候補： 個別に分解可能なソース（CW候補）は、一般的な背景放射と比較して、より大質量で、より低赤方偏移（$z < 0.4$）に位置する傾向がある。これらは主に、大質量グループやクラスター（$M_{\text{halo}} > 10^{13} \, M_\odot$）の中心銀河に宿っている。
• 電磁気学的対応天体：
  • 性質： 大半のCW候補は、クエーサーではなく活動銀河核（AGN）として現れ、低いエディントン比（$f_{\text{Edd}} < 10^{-2}$）を特徴とする。
  • 検出可能性： 光学および近赤外線（NIR）帯域では、ほとんどのCW候補はホスト銀河に隠れてしまうため、光度曲線の変動探索は困難であり、稀で明るいシステムに偏る。
  • 最適な波長： 電波およびX線帯域が、同定のための最良の展望を提供する。$f_{\text{Edd}} > 10^{-3}$ の候補については、AGN放射がホスト銀河を凌駕する。低降着系であっても、楽観的な見積もりによれば、最大で80%（電波）および30%（X線）がホストに対して検出可能である。
• 検出確率： 10年のベースラインを想定した場合、少なくとも1つのCW候補（$S/N > 3$）を検出する確率は、EPTA級の感度で4%、SKA級の感度で20%と推定される。

意義
本論文の主な意義は、最先端の宇宙論的シミュレーションを用いて、PTA信号を理解するための統一されたマルチメッセンジャーの枠組みを提供することにある。決定的なのは、理論的なGWBモデルを観測データと比較する標準的な手法に異議を唱えている点である。有限の観測期間とスペクトルリークが、推論されるGWB振幅を人工的に増大させ、スペクトルの傾きを平坦化させることを示すことで、著者らは、天体物理学的なモデルとPTA観測との間の見かけ上のテンションが、基礎となる天体物理学の失敗というよりも、むしろ解析手法によるアーティファクト（人工的なもの）である可能性を示唆している。さらに、本研究は電磁気学的なフォローアップに関する具体的な指針を提供しており、探索はクエーサーのような系における光学的な変動のみに頼るのではなく、電波およびX線帯域を優先し、大質量銀河グループやクラスターに焦点を当てるべきであることを示唆している。