A Multimessenger Search for the Supermassive Black Hole Binary in 3C 66B with the Parkes Pulsar Timing Array

パースク Pulsar Timing Array の第 3 データリリースを用いた 3C 66B における超大質量ブラックホール連星の連続重力波探索において、電磁波観測に基づくパラメータを前提とした結果、検出は確認されなかったものの上限値が設定され、電磁波データとの結合解析手法や標準サイレンとしての応用可能性が示唆された。

要約

1. 舞台設定：宇宙の「双子の巨人」と「鼓手」

まず、物語の登場人物を整理しましょう。

• 3C 66B（さん・せー・ろく・ろく・ビー）:
  地球から約 3 億光年離れた場所にある銀河です。ここには、**「超巨大ブラックホールの双子（連星）」**がいると、電波の観測から疑われています。
  • 比喩: 宇宙の広大なダンスホールで、二人の巨大な巨人が手を取り合い、ゆっくりと回転しながら踊っている状態です。
• パルクス・パルサー・タイミング・アレイ（PPTA）:
  オーストラリアにある「パルクス電波望遠鏡」を使って、32 個の**「パルサー（中性子星）」**を監視するプロジェクトです。
  • 比喩: パルサーは宇宙に点在する**「超正確な鼓手」**です。彼らは一定のリズムで「ドンドン、ドンドン」と光（電波）を放っています。このリズムは原子時計よりも正確です。
• 重力波（CGW）:
  二人の巨人が回転する際、時空（宇宙の布）に波紋が広がります。これが重力波です。
  • 比喩: 巨人が踊るたびに、宇宙の布に「しわ」が寄ります。このしわが鼓手（パルサー）の鼓動のリズムをわずかに乱します。

2. 研究の目的：「さざ波」を探す

これまでの電波観測（光で見える情報）では、3C 66B の中心で二人の巨人が踊っている可能性が高いと推測されていました。
しかし、**「本当に重力波という『しわ』が生まれているのか？」**を、パルサーの鼓動の変化から直接証明しようというのが、この研究の目的です。

もし見つかったら、それは「最終的な 1 パーセク（約 3 光年）の問題」と呼ばれる、ブラックホールがどうやってここまで近づいたのかという謎を解く大発見になります。

3. 調査方法：「耳を澄ます」

研究者たちは、パルクス望遠鏡で集めた 18 年分のデータを分析しました。

• やり方:
  「もし 3C 66B の巨人たちが踊っているなら、パルサーの鼓動はこう乱れるはずだ」という**「予測モデル」**を作り、実際のデータと照らし合わせました。
  • 比喩: 静かな部屋で、特定の場所から「特定の音」が聞こえるかどうかを、ノイズの多い部屋で探しているようなものです。「もし巨人が踊れば、このリズムが少しズレるはずだ」という仮説を立てて、データの中にそのズレがないか探しました。

4. 結果：「見つかった？それとも見つからなかった？」

結論は少し複雑ですが、一言で言えば**「決定的な証拠は見つからなかったが、可能性は完全に消し去れなかった」**というものです。

• 発見はなかった:
  「巨人が踊っている！」という明確な証拠（重力波の信号）は、データの中には見つけられませんでした。
• しかし、否定もできない:
  「巨人は絶対にいない」とも言えません。なぜなら、もし巨人がいたとしても、その重さや距離の組み合わせによっては、今の観測機器では検出できないほど「小さなさざ波」だった可能性があるからです。
  • 比喩: 静かな部屋で「誰かが足音を立てているか？」と探しましたが、足音は聞こえませんでした。しかし、「もし足音がとても小さかったり、壁に遮られていたりしたら、聞こえなくてもおかしくない」ということです。

5. 重要な発見：「巨人の重さ」に制限をかけた

信号は見つかりませんでしたが、研究者たちは**「もし巨人がいたとしたら、どれくらい重くてはいけないか」**という制限（上限値）を導き出しました。

• 結果:
  電波観測（光のデータ）から推定されていた「巨人の重さ」の範囲のうち、**約 60% は「ありえない（重すぎる）」**と排除されました。
  • 比喩: 「もし巨人がいたら、体重は 100kg 以下でなければならない」という制限をかけました。これまでの推定では「100kg〜200kg」の可能性があると言われていましたが、今回の調査で「100kg 超えはあり得ない」と分かったのです。

6. 今後の展望：「宇宙の距離計」としての可能性

この研究のもう一つの大きな意義は、**「もし将来、この巨人の存在が確定したら、宇宙の膨張速度を測るのに使えるかもしれない」**という可能性を示したことです。

• 標準サイレン（Standard Siren）:
  重力波の強さから「距離」が、光の観測から「遠ざかる速度」が分かれば、宇宙の年齢や膨張率（ハッブル定数）を計算できます。
  • 比喩: 遠くで花火が上がったとき、その明るさ（重力波）と色の変化（光）を組み合わせることで、花火までの正確な距離と、宇宙がどれくらい速く広がっているかを測れるようになるかもしれません。

まとめ

この論文は、**「3C 66B という銀河にあるブラックホールの双子が、重力波を放っているか？」**という問いに対して、以下の結論を出しました。

1. 直接の証拠は見つからなかった（まだ「発見」とは言えない）。
2. しかし、**「もしいたとしても、電波観測で言われていたほど重くはない」**という制限を設けることに成功した。
3. 将来的には、このように「場所が分かっている」ブラックホールを、宇宙の距離を測るための**「ものさし（標準サイレン）」**として使えるようになるかもしれない。

つまり、**「今回は『見つかりませんでした』という結果でしたが、その『見つからなかった』という事実自体が、宇宙の謎を解くための重要な手がかりになった」**という、非常に前向きな研究でした。

技術概要

論文要約：パルサータイミングアレイ（PPTA）DR3 による 3C 66B における超大質量ブラックホール連星の探索

1. 研究の背景と課題 (Problem)

• 背景: パルサータイミングアレイ（PTA）は、ナノヘルツ帯の重力波背景（GWB）の検出に成功しつつあるが、個々の超大質量ブラックホール連星（SMBHB）からの連続重力波（CGW）の検出は次の重要なマイルストーンである。
• 対象: 銀河 3C 66B の中心にある SMBHB 候補は、電磁波観測（VLBI 等）から軌道運動が確認されており、PTA 周波数帯で最も強い重力波源候補の一つとされている。
• 課題: 過去の NANOGrav などの探索では検出に至らず、電磁波観測から推定されたパラメータ（特に chirp mass: 連星の質量）と PTA の感度の整合性が議論されていた。また、電磁波データと重力波データを統合的に解析する手法の確立が求められていた。

2. 手法とデータ (Methodology)

• データセット: パークス・パルサータイミングアレイ（PPTA）の第 3 回データリリース（DR3）を使用。18 年間にわたる 31 個のミリ秒パルサーの観測データ（2004 年〜2022 年）を解析対象とした。
• 信号モデル:
  • 円軌道を描く非スピン SMBHB から発せられる連続重力波（CGW）をモデル化。
  • 「地球項（Earth term）」と「パルサー項（Pulsar term）」の両方を考慮した解析を行い、特にパルサー項を含めることで感度を向上させた。
  • 静止位相近似（Stationary Phase Approximation）を採用し、観測期間中の周波数変化は無視可能と判断した。
• 統計的解析:
  • ベイズ推論: 電磁波観測から得られたパラメータ（位置、周波数、質量など）を事前分布（Prior）として設定し、モデル選択（検出の有無）と上限値の設定を行った。
  • 頻度論的アプローチ: F 統計量（$F_e$）を用いた検定を行い、偽検出の可能性を評価。
  • 電磁波データの再解析: Iguchi et al. (2010) の 3mm フラックス変動データをベイズ的に再解析し、重力波事後分布との直接比較を可能にした。
• ノイズモデル: PPTA DR3 の標準的なノイズモデルに加え、重力波背景（GWB）に起因する共通スペクトルノイズを考慮した。

3. 主要な結果 (Key Results)

• 検出の有無:
  • ベイズ因子（$\ln B$）は $-0.0027(7)$ であり、信号の検出は確認されなかった。
  • 頻度論的検定でも $F_e = 4.90$（偽警報確率 FAP = 0.50, $p \approx 0.04$）となり、ノイズ変動と一致する結果となった。
  • 結論として、3C 66B からの CGW 信号は検出されず、かつ電磁波モデルを完全に否定するものでもない（「確認も排除もできない」状態）。
• 上限値の設定（95% 信頼度）:
  • Chirp Mass ($\mathcal{M}$): $M < 6.90 \times 10^8 M_\odot$（パルサー項を含む場合）。
  • 特性ひずみ振幅 ($h_0$): $\log_{10}(h_0) < -14.44$。
  • これらの上限値は、電磁波観測から推定されたモデル（$\mathcal{M} \approx 7.9^{+3.8}_{-4.5} \times 10^8 M_\odot$）の不確実性領域の約 58% を排除した。
• 周波数依存性: 3C 66B の推定周波数（約 60 nHz）において、PPTA は最も感度が高いことが確認された。

4. 新たな貢献と提言 (Key Contributions & Significance)

• マルチメッセンジャー・ジョイント・尤度法の提案:
  • 電磁波データと重力波データが独立であるという性質を利用し、両者の尤度を掛け合わせた「ジョイント・尤度（Joint Likelihood）」を構築する新しい手法を提案した。
  • この手法により、電磁波観測で推定されたパラメータ（質量など）を重力波解析と直接統合し、より厳密な制約を課すことが可能になる。
• 標準サイレン（Standard Siren）としての応用可能性:
  • 天体の位置が既知である「標的探索（Targeted Search）」を用いれば、SMBHB 候補を標準サイレンとして扱える可能性を示唆した。
  • 重力波信号の検出（または非検出による距離の制約）と既知の赤方偏移を組み合わせることで、ハッブル定数（$H_0$）への制約が可能になることを示した（ただし、今回は検出されなかったため、電磁波モデルが正しい場合の仮説的な制約を示すにとどめた）。
• 宇宙論への示唆:
  • 3C 66B の電磁波モデルが予測するひずみ振幅は、期待される重力波背景（GWB）の振幅よりも大きい「明るい源」である可能性が高い。今回の非検出は、GWB の振幅モデルや SMBHB の数密度に関する議論にも寄与する。

5. 結論 (Conclusion)

本研究は、PPTA DR3 データを用いて 3C 66B における SMBHB 候補を初めて標的探索したものである。結果として重力波信号は検出されなかったが、電磁波モデルの質量パラメータ空間の半分近くを排除し、より厳しい上限値を設定することに成功した。また、電磁波と重力波データを統合する新しい解析手法を提案し、将来の PTA による高精度宇宙論（ハッブル定数の測定）への道筋を示した。今後のより長い観測データや、より高感度な SKA（Square Kilometre Array）時代において、この手法が確実な SMBHB 検出と宇宙論への応用につながることが期待される。