\texttt{calypso}: a Parameter-Conditioned Stochastic Surrogate Model for Circumbinary Accretion Time-Series

技術的概要：calypso – 連星周回降着の時間系列のためのパラメータ条件付き確率論的代理モデル

問題提起
Vera C. ルビン天文台の「宇宙と時間の遺産調査（LSST）」のような遷移天体調査における超大質量ブラックホール連星（MBHBs）の同定は、活動銀河核（AGN）の光度曲線における周期的な変光を検出することに依存している。現在の探索手法は、しばしば変光が準正弦波であると仮定している。しかし、流体力学シミュレーションは、連星周回降着率が、特に離心率を持つ連星において、極めて非正弦波的で、バースト的、かつ鋸歯状の形態を示すことを実証している。さらに、降着の変動は、短期的な確率的変動と、円盤歳差運動に起因する長期的な変調の両方の影響を受ける。既存のモデルは、連星パラメータ（離心率 $e_b$ と質量比 $q_b$）に対するこの複雑で非正弦波的な依存関係を捉えられておらず、検出の見落としや誤分類につながる可能性がある。集団研究やパラメータ推論を促進するために、$(e_b, q_b)$ パラメータ空間全体にわたって現実的な降着時間系列を生成できる、高速な確率論的代理モデルが必要である。

手法
著者らは、2 次元流体力学シミュレーションの suite 上で訓練されたパラメータ条件付き確率論的代理モデル calypso（"Circumbinary Accretion Lightcurves Yielded via Probabilistic Spectral Operators"）を提示する。

1. 訓練データ: このモデルは、$e_b \in [0, 0.8]$ および $q_b \in [0.1, 1.0]$ をカバーする 100 件の流体力学シミュレーション上で訓練されている。各シミュレーションは、5,000 連星軌道にわたる連星成分への降着率（$\dot{M}_1, \dot{M}_2$）および総降着率（$\dot{M}_b$）を追跡している。
2. データ前処理: 長期的な円盤歳差運動のドリフトから短期的な変動を分離するために、時系列はそれぞれ 10 軌道分の 500 のウィンドウに分割される。各ウィンドウ内では、データはパワースペクトルの支配的な周波数に基づいてトレンド除去され、位相合わせが行われる。$\dot{M}_b, \dot{M}_1,$ および $\dot{M}_2$ の時系列は特徴ベクトルとして連結される。
3. PCA 次元削減: 連結された訓練データに対して、グローバルな主成分分析（PCA）が実行される。著者らは、3,000 時間ステップから管理可能な潜在空間へ次元を削減しつつ、分散の 90% 以上を捉えるために $k=142$ の主成分を保持すると結論づけている。
4. 確率論的モデリング: 訓練セット内の各 $(e_b, q_b)$ 対について、500 のウィンドウにわたる PCA 係数は、多変量ガウス分布からの抽出としてモデル化される。これは、降着流の内在的な確率論的性質と歳差運動に駆動される変調による アレイトリック不確実性 を捉える。
5. 補間: 訓練ドメイン内の任意の $(e_b, q_b)$ 対に対する降着率を予測するために、共分散行列の平均ベクトルとコレスキー因子が線形補間される。
6. エピステミック不確実性: 著者らは、平均係数の空間的変動に基づく分散項を追加することで、限られた訓練データと補間に起因する エピステミック不確実性 をモデル化しようとした。しかし、検証により、この項が予測分散を実際のテストデータの散らばりを超えて過大評価することが示された。その結果、最終モデルはエピステミック項を除外し、アレイトリック共分散のみに依存している。
7. 推論能力: 付録 A に示される閉形式のガウス尤度が導出され、観測された降着率時間系列から $(e_b, q_b)$ の直接ベイズ推論を可能にする。

主要な結果と検証

• 性能: このモデルは、パラメータ空間を横断する 13 件の保持されたテストシミュレーションで検証された。z スコア化された平均二乗誤差（zMSE）指標を用いて、エミュレータは流体力学シミュレーションの変動パターン、周期性、および振幅を成功裡に再現した。
• 成分効率: PCA 基底は、一次元（$\dot{M}_1$）の降着率を、二次元（$\dot{M}_2$）または総降着率（$\dot{M}_b$）よりも効率的に捉えている。これは、一次元の振幅変動が小さく、パターンがより規則的であるためと考えられる。
• パラメータ空間のカバレッジ: エミュレータは、$(e_b, q_b)$ 空間の大部分で良好に機能する。特定の領域（例：$e_b=0.4, q_b=0.25$）ではわずかに高いバイアスが観測され、他の領域（例：$e_b=0.25, q_b=0.85$）ではより高い分散が見られたが、系統的な失敗領域は特定されなかった。
• ガウス性仮定: ウォッサーシュタイン距離を用いた PCA 係数分布の分析（付録 B）は、多変量ガウス性の仮定が訓練データにとって「十分良い」近似であることを確認し、偏差は一般的に標準偏差の範囲内であった。
• 合成光度曲線: 著者らは、エミュレータ出力を合成 LSST $g$ 帯見かけ等級に変換するエンドツーエンドのパイプラインを実証した。変動振幅と検出可能性は $(e_b, q_b)$ に強く依存し、離心率が高く質量が不均等な系は鋭いピークを示し、近点通過時に調査の検出限界を下回る可能性があることを示した。

意義と主張
本論文は、calypso が、離心率を持ち質量が不均等な連星パラメータ空間全体における降着変動の最初の確率論的エミュレータであると主張している。その主な意義は以下の点にある：

1. 現実性: 正弦波近似を超え、連星周回降着のバースト的で非正弦波的な性質を捉えること。
2. 確率論的性質: LSST およびパルサータイミングアレイ（PTA）検出の時代における集団研究やパラメータ推論に不可欠な、降着流に内在する不確実性を捉える原理的な方法を提供すること。
3. 有用性: 調査計画のための合成時系列の生成や、観測データからの連星パラメータの直接推論を行うための、計算的に効率的なツール（Python パッケージ経由）を提供すること。

限界と留保事項
著者らは明示的にいくつかの限界を指摘している：

• シミュレーション依存性: このモデルは、固定されたパラメータ（局所等温状態方程式、固定された円盤アスペクト比 $H/R=0.1$、固定された粘性 $\alpha=0.1$）を持つ特定の 2 次元流体力学シミュレーション上で訓練されている。実際の円盤は熱力学的に異なる可能性がある。
• 補間対外挿: このモデルは訓練グリッド内を補間するが、その外へは外挿できない。
• 時系列の長さ: エミュレータは固定された 10 軌道ウィンドウを生成する。より長い時系列は独立したウィンドウから継ぎ接ぎされる必要があるため、$P_b \gg$ の時間スケールでの長期的なドリフトは捉えられない。
• 放射パイプライン: 等級への変換は、再処理、風、相対論的効果なしの理想化された薄円盤放射モデルを使用する。これらの合成等級は、教条的な予測というよりも、転送関数の例示として意図されている。
• 時間的共分散: 主要な性能指標（z スコア化された MSE）は、各時間ビンにおける周辺分布を評価するが、時系列の時間的共分散構造における誤差には感応しない。

本論文は、calypso が重要な前進である一方で、将来的な研究は 3 次元/MHD シミュレーション、より広範な円盤パラメータ、および長期的な歳差運動変調の明示的なモデル化を取り入れることを目指すと結論づけている。