Wind-mediated Eddington-limited emission in a $10^{4}M_\odot$ Black Hole Tidal Disruption Event

技術要約：$10^4 M_\odot$ ブラックホールにおける潮汐破壊イベントでの風によるエディントン限界放射

問題提起
潮汐破壊イベント（TDE）は、最も明るい熱的トランジェントの一つであるが、初期段階の光度ピークの物理的な起源は依然として不明である。古典的なモデルは、恒星の破片がコンパクトな降着円盤へと効率的に円形化することを予測しているが、観測ではしばしば、より低い黒体温度、より大きな半径、そして拘束された破片のフォールバック率を厳密には追従しない光度が示される。さらに、中間質量ブラックホール（IMBH、$10^4 - 10^6 M_\odot$）の存在とそのデモグラフィクス（人口統計学的特性）は、この質量領域におけるTDEに関する堅牢なモデルの欠如により、依然として議論の対象となっている。従来のシミュレーションは、グローバルな軌道進化と、近点におけるノズル・ショックや放射伝達といった決定的な流体力学的プロセスを同時に解像するために必要な、極端なダイナミックレンジの解決に苦慮してきた。

手法
著者らは、主系列星（$M_\star = 0.5 M_\odot$, $R_\star = 0.47 R_\odot$, polytropic index $n=1.5$）が非回転のIMBH（$M_{\rm BH} = 10^4 M_\odot$）によって破壊される、初のエンド・ツー・エンドの3次元放射流体力学シミュレーションを提示している。シミュレーションにはRICHコードを用い、グレイ放射輸送のためのフラックス制限拡散（FLD）スキームと、水素およびヘリウムのイオン化・再結合を含む現実的な状態方程式を採用している。

主な手法の特徴は以下の通りである：

• 重力： パチンスキー・ウィタ（Paczyński–Wiita）ポテンシャルを用いて相対論的な近日点移動を近似しており、計算効率のために $r < 0.6 r_p$ 内で人工的な平滑化を適用している。
• 解像度： 著者らは適応格子細分化（AMR）スキームを採用し、グローバルな3D TDEシミュレーションにおいて前例のない解像度を実現した。収束性を評価するために、3つの解像度レベル（Low, Middle, High）をテストした。高解像度ランでは最大 $4 \times 10^6$ セルを使用し、最小セルサイズは約 $0.04 R_\odot$、軌道面における最小セル質量は約 $1.5 \times 10^{-15} M_\odot$ に達した。
• 物理： 本シミュレーションは、軌道エネルギーから熱および放射エネルギーへの変換を追跡し、不可逆的な散逸率と光球の進化を明示的に計算している。

主要な結果
シミュレーションの結果、IMBHのTDEは超巨大ブラックホール（SMBH）のイベントと比較して、明確に異なる進化経路を辿ることが明らかになった：

1. 非効率な円形化： 恒星の破片は効率的に円形化できない。近点でのノズル・ショックの形成にもかかわらず、軌道エネルギーの散逸は、シミュレーション期間内（$\approx 2.2 t_{\rm fb}$）にコンパクトな降着円盤を形成するには不十分であった。拘束された物質は高い離心率（$e \approx 0.93$）を保持し、高度に楕円の軌道に留まっている。
2. 風による媒介放射： 円盤の代わりに、近点付近に低密度の放射駆動風が形成される。この流出物は準球状に膨張し、より小さな半径で生成された放射を外側へと運ぶ。この風は、シミュレーションの終了時までに、非拘束状態かつ放射優勢（$E_{\rm rad} / E_{\rm th} \sim 10$）となる。
3. エディントン限界の光度： ボロメトリック光度曲線は、$t \approx 1.5 t_{\rm fb}$ 付近で約 $2 L_{\rm Edd}$（$\sim 8 \times 10^{41}$ erg/s）に達した後、エディントン限界（$L_{\rm Edd} \approx 3 \times 10^{41}$ erg/s）付近に落ち着く。
   • 決定的なことに、総散逸率は放出される光度を1桁以上上回っている。半径 $\approx 10^{13}$ cm（$\approx 10^{14}$ cmまで膨張）に位置する光球がボトルネックとして機能し、散逸されたエネルギーの即時的な放射を妨げている。
   • ピーク光度は、素朴な推定値（$L \sim 0.1 \dot{M}_{\rm fb} c^2 \sim 10^{46}$ erg/s）よりも数桁低くなっており、この領域において衝撃散逸が極めて非効率であることを裏付けている。
4. 収束性： 近点におけるノズル・ショックは未解像である可能性が高いが、グローバルな指標（光球半径、光度曲線の形状、総散逸率）は、3つの解像度レベルにわたって定性的に収束しており、概ね定量的に堅牢である。解像度の主な違いは、高解像度ランにおける飽和と質量放出の開始のタイミングの遅れである。

意義および観測的含意
本論文は、この研究がIMBH TDEの初期放射に対する自己整合的な説明を提供するものであると主張しており、その原因を降着円盤の形成ではなく、衝撃駆動の風による媒介放射に求めている。結果は以下を示唆している：

• IMBHの検出可能性： シミュレートされた光度曲線と温度（$T \sim \text{数} \times 10^4$ K）は、光学・UV観測と一致している。著者らは、次世代のサーベイであるVera Rubin天文台のLSST（gおよびrバンド）がこのようなイベントを $z \approx 0.1$ まで、ULTRASAT（近紫外）が $z \approx 0.06$ まで検出できる可能性があると見積もっている。
• モデルの制約： 本知見は、IMBH TDEにおいて効率的な円形化を仮定するモデルに疑問を投げかけており、円形化の欠如がこれらのイベントの堅牢な特徴であることを示唆している。これは、SMBHのTDEよりも質量注入シミュレーションスキームの適用が容易であることを意味する可能性がある。
• LFBOTとの関連性： 導出されたボロメトリック光度は、Luminous Fast Blue Optical Transients (LFBOTs) のピークよりも著しく低いため、初期段階のIMBH TDEがLFBOTの支配的な説明であるという説には否定的な立場をとっている。ただし、後期段階の円盤モデルの可能性は残されている。

著者らは、本モデルが簡略化された重力およびグレイ放射伝達に依存しているものの、「IMBH TDEは、非効率な円形化を伴う、膨張する放射駆動風を介したエディントン限界の放射を生み出す」という定性的な結論は堅牢であると結論づけている。今後の課題として、ノズル部分の高解像度化、マルチグループ放射伝達、および磁気流体力学（MHD）への対応が挙げられている。