Constraining Orbital Eccentricity of a Supermassive Black Hole Binary Candidate PKS 2131-0211

技術要約：超巨大ブラックホール連星候補 PKS 2131-021 の軌道離心率の制約

問題の所在
ブレーザー PKS 2131-021 は、数十年にわたる電波フラックス密度の正弦波状の変動を示しており、これは以前、円軌道を持つ超巨大ブラックホール（SMBH）連星系による相対論的ドップラー・ブースティングとして解釈されてきた（O'Neill et al. 2022）。円軌道モデルは観測された変動に対して運動学的な説明を提供するが、この連星系に残留的な軌道離心率が存在するかどうかは未解決の問題である。このような離心率を検出することは、銀河核における SMBH 連星の力学的進化や、ナノヘルツ重力波の源としての可能性を理解する上で極めて重要である。しかし、小さな離心率の信号を、クエーサーの光度曲線に特有の固有の確率的変動（レッドノイズ）から区別することは、分析上の大きな課題となる。

手法
著者らは、軌道離心率の影響を組み込むことで、既存の PKS 2131-021 に対する運動学的軌道モデルを拡張している。

• 解析的枠組み： パルサー連星の極めて小さな離心率に対して通常用いられる ELL1 タイミングモデル（Lange et al. 2001）を応用し、著者らは電波フラックス密度の時間変化に関する完全な解析的表現を導出した。著者らは、軌道を記述するために Laplace-Lagrange パラメータ（$\epsilon_1 \equiv e \sin \omega'$ および $\epsilon_2 \equiv e \cos \omega'$）を利用しており、これにより小離心率（$e \ll 1$）における数値的安定性を確保している。導出されたフラクショナル・フラックス変化 $\delta \ln S$ の式には、純粋な正弦波からの偏差を検出できるよう、軌道位相の第2高調波までの項が含まれている。
• データの活用： 解析には、Haystack 天文台、ミシガン大学 ラジオ天文台（UMRAO）、および Owens Valley Radio Observatory（OVRO）のアーカイブデータを含む、1975 年から 2021 年にわたる編纂された電波光度曲線を用いている。データセットは3つのエポックに分割されており、周期的な変動を示す期間に焦点を当てるため、エポック2（静穏期）は除外されている。
• 統計的手法： 著者らは、Nautilus Python パッケージを用いて、ベイズパラメータ推定およびモデル選択を行っている。彼らは以下の2つの主要なモデルを比較している：
  1. 円軌道モデル。
  2. $\epsilon_1$ および $\epsilon_2$ を持つ離心軌道モデル（ELL1）。
     決定的なことに、解析は、クエーサーの固有のレッドノイズをモデル化するための減衰ランダムウォーク（DRW）プロセスを含める場合と含まない場合の両方で行われている。DRW モデルは、固有の分散、減衰タイムスケール、および標準的な DRW 行動からの偏差を考慮するための指数に関する追加パラメータを導入している。

主な貢献

1. 解析的拡張： 本論文は、PKS 2131-021 に関連する「極小離心率」領域に特化した、SMBH 連星のドップラー・ブースティング・モデルに軌道離心率を組み込むための、計算効率の高い完全な解析的処方を提供している。
2. 厳密なノイズモデリング： 本研究は、レッドノイズの影響を調査することで、離心率の制約を系統的に検証している。DRW プロセスを用いた場合と用いない場合を比較することにより、著者らは、確率的変動がいかにして小さな離心率の信号を模倣、あるいは隠蔽し得るかを実証している。
3. 包括的なデータセット解析： 本研究は、PKS 2131-021 に関して利用可能な最も広範なデータセット（45年以上）を利用しており、結果の堅牢性をテストするために様々な観測装置データの組み合わせを解析している。

結果

• レッドノイズモデリングなし： レッドノイズを考慮せずにデータを解析した場合、ベイズ証拠は円軌道モデルよりも離心軌道 ELL1 モデルを支持した。全結合データセット（OVRO + Haystack + UMRAO）において、解析の結果、残留軌道離心率は $e = 0.053 \pm 0.015$ となった。自然対数によるベイズ因子（$\ln \mathcal{B}$）は $+3.15$ であり、この特定の文脈において離心軌道モデルを強く支持している。
• レッドノイズモデリングあり（DRW）： 確率的変動を考慮するために DRW プロセスを組み込んだ場合、傾向が変化した。すべてのデータセットの組み合わせにおいて、円軌道モデルが離心軌道モデルよりも統計的に好ましいものとなった。離心＋DRW モデルと円軌道＋DRX モデルとのベイズ因子は負の値（$-0.13$から$-0.20$の範囲）であり、離心モデルへの支持はないことを示している。この領域において、著者らは離心率の上限値を設定することしかできず、最も保守的な 95$e < 0.15$ である。
• 周期の安定性： DRW プロセスの導入とそれに伴うパラメータの不確実性の増大（特に軌道周期について）にもかかわらず、Circular+DRW モデルは一貫してコヒーレントな周期信号を回復している。軌道周期は、すべてのデータセットにおいて明確に定義されており（観測者系で約 1731 日）、信号の周期的な性質を強化している。

意義と主張
本論文は、レッドノイズを無視した場合、データは当初小さな非ゼロの離心率（$e \approx 0.05$）を示唆していたものの、より厳密なデータの確率的特性の取り扱いにより、円軌道モデルが好まれることを主張している。著者らは、PKS 2131-021 の好ましい記述は、強い確率的レッドノイズの中に埋め込まれた円軌道の SMBH 連星であると結論づけている。

本研究の意義は、観測された正弦波状の変調が現実的なノイズモデルの導入に対しても堅牢であることを実証した点にあり、これにより PKS 2131-021 がナノヘルツ重力波の信頼できる源であるという説をさらに補強している。著者らは、現在の離心率の証拠は検出基準（$\ln \mathcal{B} > 8$）を満たしておらず、また、より高次の項として現れるはずの近日点移動の効果も、まだ測定可能ではないと述べている。彼らは、離心率や近日点移動率を確実に測定するためには、OJ 287 に対して行われたようなジェットの位置角の測定を含む、将来的な取り組みが必要であると考えている。本研究は、ブラーザーの光度曲線における周期信号を解釈する際に、離心率に関する偽陽性を避けるためにレッドノイズをモデル化することの必要性を強調している。