Radiation-hydrodynamics of star-disc collisions: From system parameters to outflows and lightcurves

技術要約：星と円盤の衝突における放射流体力学：システムパラメータからアウトフローおよびライトカーブまで

問題と背景
準周期的な爆発（QPE）は、超巨大ブラックホール（SMBH）に関連する核トランジェントであり、繰り返される軟X線フレアを特徴とする。いくつかのモデルが存在するが、星がタイトで離心率が高く、かつ傾斜した軌道上で降着円盤を繰り返し横切る「星–円盤衝突シナリオ」は、有力な候補の一つである。このモデルは、GSN 069のような天体で観測される再帰時間、交互に現れる「長い・短い」間隔、および「強い・弱い」フレアパターンの性質を自然に説明できる。しかし、これらの衝突の局所的な物理パラメータを、観測可能なフレア特性（光度、持続時間、非対称性）に結びつけるには、膨張する光学的に厚い媒体内における衝撃波、アウトフローの幾何学、および放射輸送の結合した進化を理解する必要がある。先行研究では、1次元、2次元、あるいは非放射的アプローチを用いてこの問題の一部を探求してきたが、システムパラメータとライトカーブを直接結びつける専用の系統的な3次元放射流体力学研究は不足していた。

手法
著者らは、平滑化粒子流体力学（SPH）コードであるPhantomを用い、一連の3次元局所放射流体力学シミュレーションを実施した。シミュレーションは、衝突のタイムスケールが軌道タイムスケールよりもはるかに短いため、全域的な軌道力学、一般相対論、およびガスの自己重力を無視し、円盤中間面を中心とした局所的なデカルト座標系で行われた。

本研究では、$10^6 M_\odot$のSMBH周囲の、光学的に厚く放射圧が支配的な$\alpha$円盤と衝突する主系列星（$R_\star \approx R_\odot$）に基づく標準的なモデルを用いた。星は、速度 $v_\star = 0.1c$ で移動する剛体球としてモデル化された。著者らは、標準モデルに対して以下の5つの主要パラメータを一つずつ変化させながら、系統的に検証を行った：

1. 恒星速度 ($v_\star$): 標準値の0.5倍から1.5倍の範囲。
2. 円盤表面密度 ($\Sigma_d$): 標準的な薄い円盤から高密度の潮汐破壊イベント（TDE）後の円盤まで、標準値の1倍から100倍の範囲。
3. 恒星半径 ($R_\star$): 標準値の0.5倍から4倍の範囲。これは実質的に、円盤スケールハイトと恒星半径の比（$H/R_\star$）を変化させることに相当する。
4. 垂直方向の円盤密度プロファイル: ガウス分布 $\rho(z) \propto e^{-z^2/2\sigma^2}$ としてモデル化し、中間面への質量集中度（$\sigma/R_\odot$）を変化させた。
5. 局所的な衝突角 ($i$): 相対速度ベクトルと円盤中間面とのなす角。$90^\circ$（垂直）から$30^\circ$（斜め）までの範囲。

放射輸送は、電子散乱不透明度を用いたフラックス制限拡散近似によって扱われた。出現するボロメータ光度は、動的に決定されたフォトスフェアを通じて、拡散フラックスと移流フラックスを積分することによって算出された。

主な貢献と結果
本研究は、システムパラメータの関数としてのピーク光度（$L_{\text{peak}}$）およびフレア持続時間（$\Delta t$）に関する経験的なスケーリングを提供している。主な知見は以下の通りである：

• 衝突ダイナミクス: 星はボウショック（弓状衝撃波）を駆動し、円盤の両側に2つのアウトフロー（前方および後方）を生成する。前方アウトフローは一般に、より質量が大きく、よりエネルギーが高い。ボウショックおよびアウトフローの定性的な形態は、$v_\star$や$\Sigma_d$の変化に対して概ね不変であるが、熱力学的な正規化定数は変化する。
• 恒星速度 ($v_\star$): 高速の星は、より高い衝撃加熱率により、より明るいフレア（$L_{\text{peak}} \propto v_\star^2$）を生み出す。また、より多くのエネルギーが注入されるため、より長いフレア（$\Delta t \propto v_\star^{1.4-2.2}$）を生み出すが、強い断熱損失により光度はより速く減衰する。
• 円盤表面密度 ($\Sigma_d$): 密度の高い円盤はアウトフローの質量と運動量を増加させるが、ピーク光度への影響は弱い（前方では $L_{\text{peak}} \propto \Sigma_d^{-0.1}$、後方では $L_{\text{peak}} \propto \Sigma_d^{0.2}$）。これは、衝撃加熱の増加と、光学的に厚い放出物における強い断熱損失との間の競合によるものである。しかし、密度の高い円盤は、光子の拡散時間が長くなるため、フレアの持続時間を大幅に延長させる（$\Delta t \propto \Sigma_d^{0.5-0.6}$）。
• 恒星半径 ($R_\star$): 大きな星はより多くのガスと相互作用し、より明るい（前方では $L_{\text{peak}} \propto R_\star^{1.5}$、後方では $L_{\text{peak}} \propto R_\star^{0.7}$）かつ、より長い（$\Delta t \propto R_\star^{1.5}$）フレアを生み出す。$R_\star$の増加は、前方と後方のアウトフロー間の運動量の非対称性を強化する。
• 垂直密度プロファイル ($\sigma$): より中心に集中した円盤（より小さい $\sigma$）は、より明るく短いフレアをもたらす。高密度の中心層はより多くの衝撃エネルギーを吸収する一方で、低密度の外層は断熱損失を抑えつつ、放射をより効率的に脱出させる。
• 衝突角 ($i$): より斜めの衝突（より小さい $i$）は、2つのアウトフロー間の運動量および光度の非対称性を減少させる。後方アウトフローは著しく明るくなり（$L_{\text{peak}} \propto (\sin i)^{-0.8}$）、相互作用時間が長くなることと、衝撃を受けたガスの再分配がより対称的になることにより、フレアはより長くなる（$\Delta t \propto (\sin i)^{-0.9 \text{ to } -1}$）。

GSN 069への適用
著者らは、これらの経験的なスケーリングを、よく特性化されている天体であるGSN 069に適用した。観測されたピーク光度（$\sim 3 \times 10^{42}$ erg s$^{-1}$）、持続時間（$\sim 1$ 時間）、および強・弱の振幅比に一致するようにパラメータ空間をサンプリングすることで、2つの実行可能な領域を特定した：

1. 推奨される解: 密度が高く、垂直方向に集中したポストTDE円盤（$\Sigma_d \sim 45 \times$ 標準値、$\sigma/R_\odot \sim 0.3$）と衝突する逆行軌道。このモデルは観測されたフレア特性を良好に再現している。
2. 代替解: 標準的な薄い円盤密度を持つが、極めて浅い衝突幾何学（$i \sim 6^\circ$）を持つ逆行軌道。この解は精度が低く、極端なパラメータを必要とする。

著者らは、高密度で垂直方向に集中した円盤を好むという結果が、GSN 069におけるTDE的な活動を示す観測証拠と一致していると述べている。

意義と限界
本論文は、局所的な衝突パラメータとQPEのライトカーブを結びつける最初の系統的な3次元放射流体力学研究であると主張している。導出された経験的なスケーリングは、局所的なシミュレーションと全域的なQPEモデルとの間の架け橋として機能し、システムパラメータに基づいたフレア特性の予測を可能にする。著者らは、特定のパラメータの組み合わせが、QPEの特徴的な振幅、持続時間、およびデューティサイクルを再現できることを示唆している。

本研究は、モデルが分離可能な単一パラメータのスケーリングを仮定していること、円盤の回転と剪断を無視していること、星を剛体球として扱っていること、および視線方向の効果やスペクトル進化（例：コンプトン化）を考慮していないことを認め、その限界を謙虚に述べている。したがって、結果は一意で決定的なベストフィットモデルとしてではなく、物理的に動機付けられたパラメータ空間の領域として提示されている。著者らは、将来の研究において、現実的な恒星構造、円盤ダイナミクス、および視線方向依存の光度を取り入れることで、多様な観測対象となるQPEソースを再現するモデルの能力を完全に定量化できることを示唆している。