General Relativistic Shock Wave Solutions with Black Hole Formation: The Singular Isothermal Sphere Case

技術要約：ブラックホール形成を伴う一般相対論的衝撃波解

問題提起
JWSTによって明らかになった高赤方偏移（$z \gtrsim 6$）における超大質量ブラックホール（SMBH）の急速な出現は、効率的な形成経路を必要としている。特異等温球（SIS）がブラックホールへと滑らかに連続的に崩壊する過程は、Cai & Shu (2005, 以下CS05) によって完全な一般相対論（GR）の枠組みで解析的に解かれているが、そのような崩壊に伴って必然的に発生する衝撃波については、これまで解析的に導出されてこなかった。恒星質量程度の種（シード）からのエディントン限界内での降着モデルでは、観測された大質量ブラックホールの存在量を、最初の500 Myr以内に説明することに苦慮している。さらに、これらの種（シード）の形成に伴うエネルギー放出——潜在的に、誕生したばかりのブラックホールの周囲に熱放射や高密度な繭（コクーン）を駆動するもの——は、GR領域において解析的に未探索のままである。

手法
著者らは、不連続な流れを含むためにCS05の自己相似フレームワークを拡張している。本研究では以下の手法を用いている：

1. 自己相似フレームワーク： 著者らは、変数 $\zeta \equiv r/t$ に依存する計量を用いた自己相似シュバルツシルト（SSS）定式化を利用している。流体は等温状態の方程式 $P = \gamma\rho$ （ここで $\gamma \equiv c_s^2$ は音速の二乗）に従うと仮定される。
2. 衝撃波ジャンプ条件： 本論文では、不連続な面を横切る等温流体の一般相対論的ジャンプ条件を導出している。衝撃波静止系（SRF）におけるバリオン数フラックスおよび運動量フラックスの保存則を適用し、それをSSSテトラド・フレームへ変換することにより、著者らは相対論的なランク・ウゴニオ条件の一般化を導出している。これにより、許容される衝撃波解を決定する前段（プレショック）速度に関する四次方程式が得られる。
3. 座標整合： 外部エンベロープ（外層）と内部崩壊の間の隔たりを埋めるため、著者らは「ゼロ速度面」（SSSにおいて流体速度 $v=0$ であり、かつ共変量としての計量の微分 $\Omega=0$ となる場所）を用いた座標整合技術を導入している。この一意な軌跡により、共変量による記述の積分定数をSSS解によって固定することが可能となり、エンベロープから特異点に至るまでの時空の完全な記述が可能となる。
4. 数値積分： 著者らは、臨界曲線（滑らかな通過が行われる場所）から外側に向かって、導出されたジャンプ条件を介して、様々な外部境界条件（静水圧、膨張、または崩壊するエンベロープ）と一致するように流体方程式を積分している。

主要な貢献

• GRジャンプ条件の導出： 本論文は、完全な一般相対論における等温流体の衝撃波ジャンプ条件に関する初の解析的導出を提供しており、ニュートン極限（Tsai & Hsu 1995 および Shu et al. 2002 の結果を回収）との整合性を検証している。
• 衝撃波解： 著者らは、内部崩壊と外部エンベロープを接続する、豊かな自己相似衝撃波解の族を特定した。これには、静水圧エンベロープを持つ解のほか、膨張（EE）、崩壊（EC）、あるいは「ブリーズ（微風）」挙動を示す動的なエンベロープが含まれる。

3. 座標の架け橋： ゼロ速度面を用いることで、シュバルツシルト記述と共変量自己相似記述を一意に接続する新しい手法が提示されており、成長するブラックホールの時空の特性付けを完結させている。
4. エネルギー収支分析： 著者らは、衝撃波面におけるエネルギー放出量と、滑らかな崩壊シナリオと比較した中心部への降着率の抑制を定量化している。

結果

• 衝撃波速度とマッハ数： 衝撃波解は光速の最大 $\sim 40\%$ の速度で伝播する。マッハ数は、音速（$\gamma$）が増加するにつれて単調に減少する。明確な相対論的限界が見出された。静水圧エンベロープの場合、$\gamma \gtrsim 0.154$（$c_s \approx 0.39c$ に相当）では、上流の流れが超音速を維持できなくなるため、衝撃波解は存在しない。
• 外部エンベロープの型： 利用可能な外部解の種類は、$\gamma$ の増加とともに減少する。$\gamma \lesssim 0.13$ では、3つの動的エンベロープ型（EC、EE、ブリーズ）のすべてが可能である。$\gamma \gtrsim 0.14$ では、エンベロープ崩壊（EC）解のみが残る。
• 降着率の抑制： 中心部の降着率は、内部の崩壊のみによって決定される。衝撃波の存在は、CS05の滑らかな膨張波解（EWS）と比較して、ブラックホールへの質量供給を $\sim 5\text{--}7$ 倍抑制する。これは物理的には、シェルが衝撃波によって外側に押し出され、減速した後、後に再び落下してくるという「噴水（ファウンテン）」効果として解釈される。
• エネルギー抽出効率： 衝撃波面で放射されるエネルギーは、特定のパラメータ範囲（例：$\gamma \approx 0.14$ のECエンベロープ）において、包含される静止質量の $\sim 10\%$ に達する（$E_{shock} \approx 0.1 M_{enc} c^2$）。この効率は、シュバルツシルト降着の放射効率（ISCOにおける5.7%）のほぼ2倍である。
• 音速への依存性： エネルギー抽出の効率は $\gamma$ とともに増加し、低 $\gamma$ での $\sim 0.5\%$ から、音速カットオフ付近での $\sim 10\%$ まで変化する。

意義と主張
本論文は、衝撃波を伴う直接崩壊ブラックホール形成に関する初の解析的なエネルギー収支を提供すると主張している。結果は以下を示唆している：

1. SMBHシードの組み立て： 衝撃波メカニズムは、「ヘビーシード」形成への経路を提供する。ここでは、重力結合エネルギーの大部分が中心のブラックホールではなく、周囲のガスに注入される。これにより、破滅的な過剰供給を防ぎつつ、後の成長のための高密度なガス貯留層を維持し、降着率を制御できる可能性がある。
2. 高赤方偏移の観測： これらの解によって予測される、加熱され衝撃波に支えられたエンベロープは、JWSTによって観測されている「リトル・レッド・ドット（LRD）」の周囲に推定される高密度な繭（コクーン）を説明できる可能性があり、それらの熱放射を駆動し、周囲のガスの分散を遅らせる可能性がある。
3. 相対論的トランジェント： 星質量崩壊シナリオ（例：コラプサー）において、衝撃波によるエネルギー放出（$M \sim 10 M_\odot$ の場合 $\sim 10^{54}$ erg）は、磁気駆動またはニュートリノ駆動のメカニズムを補完する、ガンマ線バーストのエネルギー論に寄与する可能性がある。
4. 解析的ベンチマーク： これらの解は、強重力場における一般相対論的流体力学（GRHD）コードのベンチマークを行うための重要な解析的テストベッドとして機能する。数値シミュレーションでなければアクセスできない衝撃波速度やエネルギー放出の閉形式の表現を提供している。

著者らは、これらの結果が球対称性と等温状態の方程式という理想化された仮定の下で導出されたものであることを強調し、実際の天体物理学的崩壊には回転、磁場、非平衡熱力学が関与することを認めている。しかし、彼らはSISの枠組みが、ブラックホール形成のエネルギー論を理解するための、必要不可欠な第一原理に基づく基準を提供すると主張している。