An equal mass ratio supermassive binary black holes in Q J0158-4325 with periodic microlensing signature?

技術要約：Q J0158-4325における周期的なマイクロレンズ・シグネチャーを伴う等質量比超大質量連星ブラックホール

問題提起
超大質量連星ブラックホール（SMBBH）の検出は、階層的な銀河形成モデルによって予測されているものの、直接的に解像することが困難であるため、天体物理学における極めて重要な課題となっている。レンズ効果を受けたクエーサー Q J0158-4325（赤方偏移 $z=1.29$）は、光度曲線において約173日の周期的な振動を示しており、これは Millon et al. (2022) によって、主円盤内の軌道ホットスポットによるマイクロレンズ・シグネチャーであると解釈されていた。本研究の目的は、より物理的に現実的な「トリプル・ディスク（三層円盤）」構造（サーカムバイナリ円盤と2つのミニ円盤）を持つSMBBHシステムを用いて、観測された周期的なマイクロレンズ変動と、光学およびUVデータを含むシステムのスペクトルエネルギー分布（SED）を同時に説明できるかどうかを厳密に検証することである。

手法
著者らは、幾何学的モデリング、放射伝達、およびマイクロレンズ・シミュレーションを組み合わせた統一的な理論的枠組みを構築している。

1. トリプル・ディスク・モデル： システムは、質量比 $q = M_2/M_1$ と分離距離 $a_{BBH}$ を持つ連星ブラックホールとしてモデル化されている。降着構造は、サーカムバイナリ円盤（内縁半径 $r_{in,c}$ で切断される）と、各ブラックホールの周囲にある2つのミニ円盤（ロッシュ・ローブ半径で切断される）で構成される。放射は、標準的な薄い円盤モデル（Shakura & Sunyaev 1973）を用い、内縁切断半径とエディントン比（$f_E$）を考慮した多色黒体放射として計算される。
2. マイクロレンズ・シミュレーション： レンズパラメータ（$\kappa, \gamma$）と恒星質量関数を取り入れたレイ・シューティング法を用いて、レンズ像Bの増幅マップを作成する。回転するトリプル・ディスク・システムの時間依存的な表面輝度は、これらのマップと畳み込まれ、マイクロレンズ光度曲線がシミュレートされる。
3. スペクトル・モデリング： 3つの円盤からのフラックスを積分することで、理論的なSEDを生成する。モデルを観測と比較するため、銀河系内の消散、マクロレンズ増幅、および Fe II 擬似連続体放射（I Zw 1 のテンプレートを使用）に対する補正を含めている。
4. 統計解析： パラメータ（$M_{\bullet}, q, f_E$, 円盤サイズ）を制約するために、ベイズ解析およびマルコフ連鎖モンテカルロ法（emcee を使用）を用いた MCMC フィッティングを採用している。研究では、2つの軌道周期仮説、$P_{orb} = P_{obs}$（低質量比）と $P_{orb} = 2P_{obs}$（高質量比、変動が両方のミニ円盤の結合放射に由来する場合）を比較している。

主な結果

• 光度曲線のフィッティング： モデルは、約173日の周期的なマイクロレンズ変動を再現することに成功した。低解像度のフィッティングでは複数のモデルが同等であることが示唆されるが、高解像度のフィッティングでは、$P_{orb} = 2P_{obs}$（$q > 0.5$ を示唆）のモデルが、特に交互に現れるピーク構造において、光度曲線の形状に対して一般により優れた適合を示すことが明らかになった。ただし、データの欠落により、光度曲線のみからの強い制約には限界がある。
• SED による制約： 光学スペクトルのみのフィッティングでは、低質量比モデルと高質量比モデルの両方が許容される。しかし、HST の UV データを取り入れることで、低質量比モデル（$q < 0.5$）は UV フラックスを過大評価してしまうため、明確に否定される。物理的な事前分布を適用した際のベイズ因子（$BF_s \approx 14$）は、$P_{orb} = 2P_{obs}$ のシナリオ、すなわち高質量比のシナリオを支持する強い証拠を与えている。
• 最適フィット・パラメータ： 最も可能性の高いモデル（例：Model Ma）は、全ブラックホール質量 $M_{\bullet} \approx 2.8 \times 10^9 M_{\odot}$ と、ほぼ等質量な質量比（$q \approx 1$）を示している。連星の分離距離は約 $30 r_g$（重力半径）であり、これは重力波放射による合体までのタイムスケールが数十年であることを意味する。
• X線予測： モデルは、X線マイクロレンズが光学変動と同じ周期と位相を示すが、よりコンパクトな X 線放射領域のために振幅が大きくなることを予測している。しかし、現在の X 線観測は、これらの予測を確認したり追加の制約を課したりするための、時間分解能や精度が不足している。

意義と結論
本論文は、多波長マイクロレンズ・シグネチャーと詳細なスペクトル・モデリングを組み合わせることが、光度曲線解析単独では達成できない SMBBH システムへの強固な制約を提供することを実証している。本研究は、Q J0158-4325 が、おそらく等質量な SMBBH システムであり、トリプル・ディスク構造を持っていると結論付けており、これは合体の最終段階（数十年から数世紀）にあるシステムの稀な観測例である。

著者らは、このような短命な進化段階の検出は統計的に稀であるが、開発された枠組みは他のレンズ・クエーサーにも適用可能であることを強調している。また、将来の時系列サーベイ、特にルービン天文台の LSST（Legacy Survey of Space and Time）が、同様の周期的信号を系統的に探索できる可能性についても言及している。さらに、この特定のシステムの重力波周波数は現在 LISA の帯域を下回っているが、本手法は、将来の低周波重力波源に対する電磁波対応天体を特定するための経路を確立するものである。本研究は、SMBBH の特性決定における残存する不確実性を解消するために、高頻度かつ多波長でのモニタリングが必要であることを強調している。