An "inside-out" approach to modeling supermassive black hole binary inspiral

技術要約：超巨大ブラックホール連星のインスパイラル（内側から外側へ）のアプローチによるモデリング

問題提起
近年のパルサー・タイミング・アレイ（PTA）の観測は、宇宙における超巨大ブラックホール連星（SMBHB）の集団に由来すると考えられる、ナノヘルツ帯の確率的重力波背景放射（GWB）の強力な証拠を提供している。しかし、これらの信号の解釈は「最終パーセク問題」と、銀河への流入から合体に至るまで、SMBHBが横断しなければならない極めて広大なダイナミックレンジ（およそ10桁）によって困難を極めている。現在のモデルは、NG15Astro（NANOGrav 15年解析）で使用されている二重冪関数（2PL）モデルのように、固定された全タイムスケールを通じて、後期インスパイラル（PTAで観測可能な段階）を初期段階の進化に結びつける「外側から内側へ」のアプローチに依存することが多い。このような結合は、PTAが観測可能な周波数帯におけるバイナリの硬化（hardening）を駆動する特定の天体物理学的プロセスから抽出できる情報を希薄化させてしまう。さらに、純粋な重力波（GW）駆動モデルは、PTAの観測時間が長くなるにつれ、低質量側のバイナリにおいて自己矛盾が生じ始めている。なぜなら、これらのバイナリは、PTA帯の低周波端から出発した場合、GW放出のみではハッブル時間内に合体できないためである。

手法
著者は、PTAデータのパラメータ推論を最適化するために設計された、SMBHBインスパイルのための「内側から外側へ（inside-out）」の解析的フレームワークを提案している。この核心となる前提は、GWBの振幅とスペクトル形状は、バイナリが天体物理学的な硬化からGW支配の硬化へと移行する軌道半長軸 $a_{GW}$ に対して最も敏感であるという点にある。

主な手法の特徴は以下の通りである：

1. $a_{GW}$ のパラメータ化： $a_{GW}$ がグローバルなパラメータの派生的な結果となる従来のモデルとは異なり、$a_{GW}$ を主要な自由パラメータとして扱う。これは重力単位（$R_g$）で定義され、バイナリ質量（$M$）および質量比（$\eta$）を用いて $a_{GW} = a_{GW,9} M_{9}^{\alpha_{GW}} (4\eta)^{\beta_{GW}}$ とスケーリングされる。ここで $a_{GW,9}$ は質量 $10^9 M_\odot$ の等質量バイナリにおける遷移半径である。
2. 内側から外側への構造：
   • 内側領域 ($a < r_{char}$): 硬化率を冪関数（$\dot{a} \propto a^{1-\nu_{in}}$）としてモデル化し、GWレジームと境界 $r_{char}$ を接続する。
   • 外側領域 ($a > r_{char}$): 複雑な天体物理学的進化を明示的にモデル化する代わりに、外側のフェーズを、バイナリ形成から $r_{char}$ に到達するまでの時間を表す一定の遅延時間 $\tau_{out}$ として簡略化する。
   • 境界条件: 物理的な制約、具体的には硬化率 $|\dot{a}|$ が光速 ($c$) を超えてはならないこと、および $r_{char}$ がすべての関連するGW放射を確実に捉えるためにPTA周波数帯の外側に位置しなければならないことを強制する。
3. 実装: 本モデルは、銀河質量関数（GSMF）、SMBH-銀河のスケーリング関係、および銀河合体率を標準的な入力として利用し、集団合成コード holodeck2 に実装されている。

主な貢献

• スケールの分離： 本フレームワークは、後期インスパイル（PTAレジーム）を不確実な初期段階の進化から正常に分離することに成功しており、これによりPTAがGWレジームへの遷移に関する物理を直接制約することを可能にしている。
• クリティカルなパラメータの特定： 本論文は、$a_{GW}$（特に $a_{GW,9}$）が、内側の冪指数（$\nu_{in}$）や外側の遅延時間（$\tau_{out}$）よりもはるかに、GWBのスペクトル形状と振幅を支配する主要なパラメータであることを特定した。
• 物理的一貫性： 本モデルは、純粋なGW駆動モデルの自己矛盾の問題に対処しており、天体物理学的な硬化を明示的に含めることで、関連するすべての質量のバイナリがハッブル時間内に合体できることを保証している。
• 文献との比較： 本モデルを「恒星型」（$\nu_{in} = -1$）および「ガス型」（$\nu_{in} = +2$）の硬化シナリオと比較検証し、異なる天体物理学的駆動力がどのようにGWBスペクトルに現れるかを示している。

結果

• $a_{GW}$ への感度: GWBは $a_{GW,9}$ に対して非常に敏感である。観測可能なGWBは、$a_{GW,9}$ が $\sim 20$ から $6000 R_g$ の範囲にある場合にのみ生成される。この狭い範囲（約30〜100倍の制約内）において、GWBの振幅とスペクトル形状は劇的に変化する。
  • 最適範囲: 最大のGWB振幅と、理想的な $f^{-2/3}$ 冪則に近いスペクトルが得られるのは、$a_{GW,9} \sim 500-2000 R_g$ のときである。このレジームでは、天体物理学がバイナリをPTA帯へと効率的に運び、その後GWが最終的なインスパイルを駆動する。
  • 低い $a_{GW,9}$: 強力な天体物理学的硬化をもたらし、GWBの低周波側での減衰（ターンオーバー）を引き起こす。
  • 高い $a_{GW,9}$: バイナリをPTAレジームの外側で停滞させ、特に高周波側でのGWB振幅の急激な低下を招く。
• 他のパラメータへの非感度:
  • GWBは質量比のスケーリング（$\beta_{GW}$）にはほとんど影響を受けず、質量スケーリング（$\alpha_{GW}$）にも弱い影響しか受けない。
  • 境界 $r_{char}$ は、PTA帯の外側にある限り、GWBの形状に最小限の影響しか与えない。
  • 外側の遅延時間 $\tau_{out}$ は、$\tau_{out} \gtrsim 3$ Gyr でない限り、主に（$z=0$ までに合体するバイナリの数を通じて）全体の振幅に影響を与えるが、スペクトル形状には影響しない。
• 天体物理学的駆動要因:
  • 恒星型モデル ($\nu_{in} = -1$): $a_{GW}$ が停滞を避けるのに十分大きくあれば、低周波のターンオーバーを最小限に抑えた高振幅のGWBを生み出すことができる。
  • ガス型モデル ($\nu_{in} = +2$): GW放出が非常に小さな分離（$\sim 100 R_g$）においてのみ支配的となるため、より低振幅で平坦なスペクトルを生じる傾向がある。
• 質量スケーリング: GWBは $\alpha_{GW}$ や $\beta_{GW}$ に対しては非感受であるが、全合体率およびLISAイベント率はこれらのスケーリングに対して非常に敏感である。

意義と主張
本論文は、$a_{GW}$ を自由パラメータとして扱うことが、PTAデータから物理的な制約を抽出するために不可欠であると主張している。「内側から外側へ」のアプローチは、電磁波による制約が乏しい後期インスパイルの段階において、最も重要なパラメータを孤立させた、計算効率が高く解析可能なフレームワークを提供する。

著者は、本モデルが以下の点において優れていると主張している：

1. 複雑な天体物理学的モデリングの必要性と、GWB解析における単純なパラメータ推論の必要性との間の緊張関係を解消している。
2. 現在のPTAデータが、全SMBHB集団に対して純粋なGW駆動モデルが自己矛盾を起こす領域をすでに超えていることを示しており、天体物理学的な硬化を含める必要があることを実証している。
3. 将来のPTAデータセットが、$a_{GW}$ と $\nu_{in}$ を測定することで、支配的な天体物理学的硬化メカニズム（恒星 vs ガス）を制約できる可能性を示唆している。
4. LISAイベント率およびマルチメッセンジャー源の特性を予測するための基礎を提供しており、GWBは特定の質量スケーリングパラメータには非感受である一方で、それらのパラメータこそが個々の連続波源の検出率や合体率を予測する上で極めて重要であることを指摘している。

本論文は、モデルがバイナリの硬化における複雑な天体物理学を簡略化しているものの、出現しつつあるGWB信号を解釈するために必要な本質的な物理を捉えており、将来の観測戦略を導くものであるとして、控えめに結論づけている。