Towards direct imaging and orbital parameter estimation of supermassive black hole binaries with spaceborne VLBI

技術要約：宇宙VLBIによる超巨大ブラックホール連星の直接撮像および軌道パラメータ推定に向けて

問題提起
サブパーセク規模の超巨大ブラックホール連星（SMBHB）の直接的な電磁波観測は、銀河の進化と重力波背景放射の起源を理解する上で、依然として極めて重要なミッシングリンクとなっている。重力波の検出や理論モデルはこれらの系の存在を示唆しているが、決定的な電磁波による証拠はまだ得られていない。地上ベースの超長基線電波干渉法（VLBI）は、地球の直径と大気の不透明度によって根本的に制限されており、角分解能は約15 µasに制限されている。この分解能は、サブパーセク連星の特性的な分離（しばしくらいは1 µas未満）を解像したり、連星であることを確認するために必要な軌道運動を検出したりするには、しばしば不十分である。したがって、これらの系を検出するための宇宙VLビの展望を評価し、その軌道パラメータを推定するための手法を開発する必要がある。

手法
著者らは、Black Hole Explorer (BHEX) ミッションを主要なケーススタディとして、SMBHBを検出し、その軌道をフィッティングするための包括的なフレームワークを提案している。この手法は、主に以下の3つのコンポーネントで構成される。

1. 連星検出基準： 本研究では、D'Orazio & Loeb (2018) に由来する条件を採用しており、連星の各成分が独立して検出可能（共にミリ波で明るい）であり、かつ観測される軌道周期（$P_{obs}$）が、運動を検出する可能性を高めるために10年未満であることを要求している。検出は、基線上の相関フラックス密度の不一致に基づいて、連星を単一のガウス光源から区別することと定義される。
2. 軌道フィッティング・パイプライン： 動的ネスト・サンプリング（dynesty）を用いたベイズ推論アプローチが開発された。このパイプラインは、連星の運動を伝播させるためのポスト・ニュートン（PN）軌道モデル（3.5PN次まで）を利用する。モデルは相対的な光行差による時間遅延を考慮し、以下のパラメータをフィッティングする：全質量（$M$）、質量比（$q$）、軌道長半径（$a$）、離心率（$e$）、軌道傾斜角（$i$）、昇交点幅角（$\Omega$）、近日点引数（$\omega$）、および近日点通過の時期（$\tau$）。
3. 合成データシミュレーション： その有効性を評価するため、著者らはBHEXの観測（中地球軌道における3.4 mアンテナを用い、約320 GHzおよび86 GHzで動作）を、地上アレイ（次世代イベント・ホライズン・テレスコープ、ngEHTを含む）と共にシミュレートしている。可視度データを生成するために eht-imaging と ngehtsim を使用し、熱雑音、大気の影響、およびモデルの誤設定誤差（加法的および分数的なノイズ）を組み込んでいる。パイプラインは、曲線的な軌道運動を直線的な運動（帰無仮説）から区別する能力、および軌道要素を制約する能力をテストする。

主な貢献と結果

• BHEXの検出閾値： 本研究は、BHEXが全フラックス密度 $\ge$ 40 mJy かつ最小分離 $\sim$2 µas（その公称分解能6 µasを下回る）の連星系を検出できることを確立した。全フラックスが高い（0.5 Jy）場合、フラックス密度の比（$f$）が0.045という低さであっても検出は可能だが、全フラックスが低い場合には感度が低下する。
• 軌道パラメータの制約： 2年間のミッション中に3回の年次観測を行う合成データを用いて、著者らは、特定のノイズ仮定（10%の分数雑音および10 mJyの加算雑音）が満たされる場合、BHEXが観測周期 $P_{obs} \le$ 10年の連星の軌道長半径（$a$）および離心率（$e$）を、真の値から0.06 dex以内に制約できることを示している。
• アストロメトリ精度： BHEXは、地上のみのアレイ（ngEHT）と比較して、相対的なアストロメトリの不確かさを大幅に改善する。分離が10 µas付近の場合、BHEXはノイズ条件に応じて、アストロメトリの精度を3〜9倍向上させる。この精度により、分離が5から15 µas、周期が7年から28年の連星において、非線形（曲線的）な運動を自信を持って検出することが可能となる。
• 縮退と限界： 解析の結果、$a$ と $e$ はよく制約される一方で、$M$ と $i$ は縮退のために事後分布が広くなることが明らかになった。質量比（$q$）は、追加の事前分布がない限り、ほとんど制約されない。本研究は、短いミッション期間（2年）がサンプリングされる軌道の弧を制限し、長周期軌道（$>30$年）の検出を困難にし、しばしば非線形運動に対する高い偽陽性率（FAR）をもたらすことを強調している。
• 将来のシステム要件： 数十の系を超える統計的に有意なSMBHBのサーベイを実施するためには、将来の宇宙VLBIシステムには、角分解能の桁違いの向上（THEZAで提案されている1 µasを目標とする）と感度の向上が必要であると、本論文は主張している。そのようなシステムは、地上アレイから独立した密な $(u, v)$ 被覆を実現するために、複数の宇宙素子を必要とする可能性が高い。

意義
本論文は、BHEXが主にM87やSgr Aのフォトロン・リングを撮像するために設計されているものの、その能力は理論的に、サブパーセク規模のSMBHBの最初の直接的な電磁波検出を提供する可能性まで及ぶと主張している。これは、そのような系の存在を示す最初の決定的な証拠となる。しかし、著者らはその成果の範囲について控えめな表現を用いている。単一の連星を検出することは大きな節目となるが、真の科学的価値は、進化モデルを制約し「最終パーセク問題」を解決するために、統計的に有意な集団を特徴付けることにある。

本研究は、BHEXは変革的なステップではあるものの、SMBHBの宇宙論的進化を完全に理解するために必要なサーベイを実行するには、より専用で、より高性能な宇宙VLBIミッションが必要であると結論付けている。提示された軌道フィッティングの手法は、これらの候補を対象とする将来の解析パイプラインのための基礎的なツールとして機能する。