Radiation-hydrodynamics of star-disc collisions for quasi-periodic eruptions

技術要約：準周期爆発（QPE）に関する放射流体力学的な星・円盤衝突の研究

問題提起
準周期爆発（QPE）は、超巨大ブラックホール（SMBH）付近で発生する高輝度なX線トランジェントであり、数時間から数日の間隔で明るいフレアが発生することを特徴とする。有力な理論モデルの一つは、傾いた軌道を持つ恒星が降着円盤と繰り返し衝突しているというものである。これまでの研究では星と円盤の相互作用が探求されてきたが、詳細な放射輸送を欠いていたり、恒星が著しく摂動を受ける領域や、ガスが恒星の周囲を流れるには衝突時間が短すぎる領域に焦点を当てていたりすることが多かった。さらに、一部のQPE天体で見られる「強・弱」の交互のフレアパターンの物理的な起源は、衝突ダイナミクスによって完全には説明されていない。本研究は、恒星がほとんど摂動を受けず、衝撃波を受けたガスが恒星の周囲を流れるのに十分な時間がある領域において、星と円盤の衝突に対する3次元（3D）放射流体力学的な取り扱いを行う必要性に応えるものである。

手法
著者らは、平滑化粒子法（SPH）コードであるPhantomを用いて、3次元放射流体力学シミュレーションを行った。シミュレーションは、局所的な衝突ダイナミクスを分離するため、降着円盤と共に移動する参照系で行われた。

• 設定: 硬い球状の恒星（$R_\star = R_\odot$）が、一様な、光学的に厚く、放射圧が支配的な降着円盤（$H = 3R_\star$, $\rho_d = 3.7 \times 10^{-8}$ g cm$^{-3}$）を、速度 $v_\star = 0.1c$ で垂直に通過する。システムパラメータは、SMBH質量が $10^6 M_\odot$、QPE周期が4時間であることを対応させている。
• 物理: ガスは断熱状態方程式（$\Gamma = 5/3$）に従う。放射輸送は、局所熱平衡（LTE）を仮定したフラックス制限拡散近似を用いて扱われる。恒星は、ガス粒子と弾性的に相互作用する剛体としてモデル化されている。衝突のタイムスケールは軌道タイムスケールよりもはるかに短いため、SMBHからの重力および自己重力は無視される。
• 範囲: シミュレーションは、グローバルなシミュレーションの計算コストを抑えつつ、衝撃波構造とアウトフローの形成を解像するために、円盤の局所的なセクション（$12R_\star \times 2H$）に焦点を当てている。

主要な貢献
本研究は、多くのQPEに関連するパラメータ空間（円盤が十分に厚く、ガスが恒星の周囲を流れることができる領域）を特に対象とした、星と円盤の衝突に関する初の3次元放射流体力学シミュレーションを提示する。本研究は、以下の詳細な分析を提供する：

1. ボウショック（弓状衝撃波）および二次衝撃波の形成と進化。
2. 恒星によって注入された運動量とエネルギーの非対称な再分配。
3. 前方および後方アウトフローの明確な特性。
4. 結果として得られるボロメトリック光度曲線および実効温度。

結果

• 衝撃波ダイナミクス: 衝突により、恒星と共に移動し、横方向に拡大する、ほぼ放物線状のボウショックが発生する。恒星の航跡（ウェイク）では、ガスの流れが収束することによって二次衝撃波が形成される。衝撃波面の横方向への拡大は、べき乗則のスケーリング（$r_\perp \propto t^{0.4}$）に従う。これは定常状態の予測である $t^{0.5}$ よりもわずかに緩やかであり、システムが依然として定常状態に向かって進化していることを示している。
• 非対称なアウトフロー: 衝突は、円盤の反対側に向かって2つの異なるアウトフローを生み出す：
  • 前方アウトフロー（運動方向）: この流れは、恒星が円盤を脱出する際の準球状の衝撃波突破（ショックブレイクアウト）によって駆動される。これはより速く拡大し、後方フローよりも大幅に多くの質量（約10倍）を含み、より多くのエネルギーと運動量を運び、後方アウトフローよりも約2倍明るい。
  • 後方アウトフロー: この流れには、準球状の衝撃波突破成分が含まれない。これは、恒星の周囲を流れるガスが二次衝撃波によって方向を変えられたもの、あるいは円盤の上端から脱出したガスによって形成される。これは質量および光度がより低い。
• 光度と光度曲線: 前方アウトフローは、約 $10^{41}$ erg s$^{-1}$ のピークボロメトリック光度に達し、これは後方アウトフローの約2倍である。前方光度曲線は急速な上昇を示した後、緩やかな減少へと変化する。これは、初期の衝撃波突破に伴う放出物からの放射から、より長い遅延を経てキャビティから脱出するガスからの放射への移行として解釈される。実効温度は、前方フローで約 29 eV、後方フローではそれよりわずかに低い値にピークに達する。
• 運動量注入: 恒星は、その進行方向に優先的に運動量を注入する。衝撃を受けたガスの全垂直運動量は $P_z \approx -0.45 M_d v_\star$ に漸近し、これはエネルギー注入の異方性と、ガス反射の斜めの角度を反映している。

意義と主張
著者らは、衝突ダイナミクスの本質的な非対称性が、いくつかのQPE天体（例：GSN 069）で観測される「強・弱」のフレアパターンを自然に説明できると主張している。なぜなら、前方アウトフローは後方アウトフローよりも著しく明るいため、特定の角度からシステムを観測している観測者は、恒星の連続的な円盤通過に対応する、交互の明るいおよび暗いフレアを見ることになるからである。

本研究は、恒星が摂動を受けない領域において、星と円盤の衝突モデルが、観測されたピーク光度（$10^{41}-10^{43}$ erg s$^{-1}$）、タイムスケール（数時間）、およびQPEの交互の光度パターンを再現できると結論付けている。著者らは、制限事項として、剛体星の使用（質量剥離や変形を無視していること）、LTE仮定（光子が不足する領域では崩壊する可能性があること）、および全域的な円盤回転とSMBH重力の無視を挙げている。彼らは、非LTE物理学や現実的な恒星構造を取り入れることが、モデルを完全に制約し、スペクトルの硬さに関する不一致（シミュレーションの $\sim 30$ eV 対 観測の $\sim 100-200$ eV）に対処するために必要であると示唆している。