Hydrodynamical simulation of wind production from hot accretion flows in tidal disruption events

技術的サマリー：潮汐破壊現象における高温降着流からの風生成の流体力学シミュレーション

問題提起
潮汐破壊現象（TDE）における超エディントン降着相からの風はよく研究されているが、その後の長期間持続する準エディントン相からの風の特性は未だ十分に探求されていない。この後期段階では、降着流は標準的な薄円盤から高温降着流（例えば、移流優勢降着流：ADAF）へと遷移する。アクティブ銀河核（AGN）や X 線連星（XRB）における高温降着流からの風は広範にモデル化されているが、TDE は明確に異なる物理シナリオを示す：恒星の破片が直接ブラックホールの直近に注入され、大半径からの定常供給ではなく、一時的でコンパクトな流を形成するからである。本論文は、この特定の準エディントン TDE 領域における風の形成と特性を調査し、ブラックホール質量、粘性係数（$\alpha$）、および恒星破片の温度が流出ダイナミクスにどのように影響するかに焦点を当てる。

手法
著者らは、球座標系における ZEUS-MP コードを用いて、2 次元軸対称流体力学シミュレーションを実施した。シミュレーションでは、太陽型恒星（$M_* = M_\odot, R_* = R_\odot$）が $10^6 M_\odot$ および $10^7 M_\odot$ の超大質量ブラックホール（SMBH）によって破壊された後の高温降着流の進化をモデル化した。

• 初期条件: 本研究は、降着率が $\sim 0.01 \dot{M}_{\text{Edd}}$ まで低下した遅延段階（破壊後数十年から数百年）に焦点を当てた。恒星破片は円盤化半径（$R_c = 2R_t$）で注入され、そのフォールバック率は標準的な $t^{-5/3}$ 減衰に従う。
• パラメータ空間: 3 つの主要パラメータを変化させて 8 つのモデルを実行した。
  1. ブラックホール質量（$M_{\text{BH}}$）: $10^6 M_\odot$ および $10^7 M_\odot$。
  2. 粘性パラメータ（$\alpha$）: $0.01$および$0.1$。
  3. 破片注入温度（$T_{\text{inj}}$）: $1 \times 10^4$ K（断熱進化）および $2 \times 10^7$ K（等温/静水圧平衡）。
• 解析: 著者らは流入率、流出率、正味の降着率を定義し、「真の」非束縛風（ベルヌーイパラメータ $B_e > 0$）と束縛対流流出（$B_e < 0$）を区別した。速度、質量フラックス、運動量フラックス、および運動エネルギーの時間平均された半径方向および角度方向のプロファイルを解析した。

主要な結果

1. 準定常状態: 降着流は数個の粘性時間スケール以内に準定常状態へ急速に進化し、準エディントン TDE のモデリングにおける準定常近似の妥当性を裏付けた。
2. ブラックホール質量の影響: より質量の大きなブラックホール（$10^7 M_\odot$）は、降着した質量の割合を高め、より高速な風を発生させる。これは、破片がシュワルツシルト半径でスケーリングされたより小さな半径（$R_c/R_S$）で注入されるため、ガスがより強く束縛されるからである。したがって、脱出するガスはより深いポテンシャル井戸を克服するために、より高い速度を達成しなければならない。
3. 粘性（$\alpha$）の影響:
   • $\alpha = 0.1$: 効率的な角運動量移動により、**非束縛で軽度の相対論的風（$\sim 0.1c$）**が放出される。これらの風は主に注入領域外の赤道面から放出される。これらの風の運動エネルギーは約 $\sim 10^{-4} L_{\text{Edd}}$ である。
   • $\alpha = 0.01$: 非効率的な角運動量移動により、対流的に束縛された流出が生じる。ガスは渦構造を示し、負のベルヌーイパラメータを持つため、無限遠へ脱出するのではなく、流れ内で循環する。これは、同様の $\alpha$ における典型的な AGN/XRB 高温流（しばしば非束縛風を示す ADIOS シナリオ）とは対照的である。著者らは、この相違を TDE 固有の初期条件、すなわち、低い初期ベルヌーイパラメータを持つ冷たい破片がポテンシャル井戸の深部へ注入されることに起因すると結論づけた。
4. 温度の影響: 恒星破片の初期温度（$10^4$ K 対 $10^7$ K）は、降着流の構造にほとんど影響を与えない。これは、降着中に解放される重力エネルギーがガスの初期内部エネルギーを支配するためである。
5. 風の幾何学: $\alpha = 0.1$ の場合、非束縛風は広い角度範囲（$10^6 M_\odot$ の場合、極から $40^\circ$ を超える）にわたって放出され、極方向の角度と比較して、エッジオン（赤道）角度においてより多くの質量、運動量、およびエネルギーを運ぶ。

意義と応用
本論文は、これらの知見が TDE の特定の観測的特徴を説明するメカニズムを提供し、物理パラメータを制約する新たな手法を提供すると主張している。

• 遅延ラジオ輝光増大: $\alpha \approx 0.1$ における非束縛風の運動エネルギー（$\sim 10^{-4} L_{\text{Edd}}$）は、ASASSN-14ae、AT2018fyk、および AT2019teq などの TDE で観測される、破壊後 $\sim 10^3$ 日における遅延ラジオ輝光増大を駆動するのに十分である。これは、これらの特定の TDE における粘性パラメータはおそらく $0.01 < \alpha < 0.1$ の範囲にあることを示唆している。
• 中間質量ブラックホール（IMBH）の検出: シミュレーションは、IMBH（$10^4 - 10^5 M_\odot$）を伴う TDE は極めて低い降着率を生み出し、容易に風を形成することを示唆している。著者らは、球状星団や矮小銀河のラジオおよび X 線サーベイが、これらの持続的な低レベル活動を検出できる可能性を提案しており、IMBH を発見する有望な手段を提供している。
• 粘性の制約: 低 $\alpha$ における TDE と AGN/XRB の流出構造の違い（束縛対非束縛）は、後期段階の風観測を通じて TDE における粘性パラメータを制約するための新たな診断ツールを提供する。

著者らは結論として、粘性の磁場起源は未解決の課題であるが、TDE 風の流体力学的挙動は、恒星破片の初期注入条件に起因して定常状態の AGN 降着と根本的に異なる、これらの一過性現象における降着物理の明確なシグネチャを提供すると述べている。