Analytic insight into the physics of SASI II. Spiral instability of the prograde mode in a rotating stellar core

技術要約：SASI II の物理学に関する解析的知見：回転する恒星核における順行モードの螺旋不安定性

問題提起
回転する大質量星のコア崩壊において、定常降着衝撃波不安定性（SASI）は非軸対称運動を駆動し、ニュートリノおよび重力波信号に特定のシグネチャを刻み込む。これまでの研究により、回転が順行螺旋モードの成長を促進することが確立されているが、その不安定化の背後にある物理的メカニズムは、特に回転速度がSASI周波数よりも遅い領域において、未解明のままであった。さらに、既存の摂動解析の多くは、円筒座標系または球状赤道近似に依存しており、流れの極方向（ポロイダル）構造を無視していたため、回転する前駆星の数値シミュレーションを解釈する上での精度に限界があった。加えて、回転が支配的なSASIモードの振動周波数に与える影響についても、強固な解析的枠組みが欠如していた。

手法
本研究は、停滞した降着衝撃波の摂動解析を、回転および差動回転の効果を含める形で球状幾何学へと拡張する。そのアプローチは以下の通りである：

1. 定常流モデリング： 理想気体を仮定し、相対論的電子縮退圧が支配的な断熱指数 $\gamma = 4/3$ を用いて、赤道面の簡略化されたモデルを構築する。流れには密度と圧力に依存する冷却関数を含み、重力ポテンシャルはニュートン近似で扱う。
2. 線形摂動解析： 線形摂動を、極方向数 $\ell$ と方位角数 $m$ を持つ球面調和関数 $Y^m_\ell$ によって特徴付ける。解析では、質量フラックス、エントロピー、比角運動量、および渦度に関する一連の結合微分方程式を用いる。
3. 保存特性： 回転する前駆星における極方向構造の解法の困難さを克服するため、本研究では径方向渦度の保存特性を利用する。これにより、系を簡約化し、シューティング法を用いて数値的に解くことが可能となる。
4. 断熱近似： 不安定化に対する断熱過程の寄与を分離するために、断熱モデルを開発する。これにより、微分方程式系の次数を下げ、ポストショックのアコースティック領域を自己強制振動子として扱う積分方程式の定式化を容易にする。
5. 解析的近似： 数値結果を様々な衝撃波半径および回転率に対してフィッティングすることで、支配的なSAS Iモード（$\ell=m=1$ および $\ell=m=2$）の固有周波数の明示的な解析的近似式を導出する。

主な貢献と結果

• 解析的固有周波数公式： 本論文は、衝撃波半径 $r_{sh}$ が最大減速半径 $r_\nabla$ の約1.5倍を超える場合に、約10%の精度で、支配的な螺旋モード（$\ell=m=1$ および $\ell=m=2$）の固有周波数の実部（$\omega_r$）に対する解析的近似を提供している。周波数は、衝撃波半径 $r_{sh}$ と比角運動量 $L$ の関数として次のように表される：
  $$\omega_{fit}^r(\ell, m) \sim \frac{|v_{sh}|}{r_{sh}} F_{\ell,m}\left(\frac{r_{sh}}{r_\nabla}\right) + \frac{mL}{r_{sh}r_\nabla} G_{\ell,m}\left(\frac{r_{sh}}{r_\nabla}\right)$$
  ここで、$F$ および $G$ は数値的フィッティングから導かれた無次元関数である。
• 不安定化のメカニズム： 差動回転による順行モードの不安定化は、移流強制とアコースティック構造との間の**位相混合（phase mixing）**によって引き起こされることを説明する。差動回転は、順行モードに関連する渦度の径方向波長を増大させる。これにより、アコースティック構造の大きな径方向スケールとの位相整合が改善され、結合効率が向上する。
• 回転率の役割：
  • 中程度の回転領域（共回転半径が存在しない領域）では、回転周波数の急激な内向きの増加（$\propto 1/r^2$）により、位相混合への悪影響が小半径において顕著となる。
  • 共回転半径を含む強い回転の領域では、定常位相近似により、主要な移流・アコースティック結合が共回転帯と衝撃波の間で発生することが特定される。
• モードの優位性： 解析により、順行モードの成長率は比角運動量に対してほぼ線形に増加することが確認された。回転が増加するにつれ、基本モード $\ell=m=2$ が $\ell=m=1$ に対して支配的になる可能性がある。両者の周波数の比は回転率に依存し、非回転の場合と比較して最大2.3倍まで増加する可能性がある。
• 冷却の詳細への独立性： 支配的なモードの振動周波数は、衝撃波半径が十分に大きい場合（$r_{sh} \gtrsim 1.5 r_\nabla$）、ニュートリノ放出（冷却関数のパラメータ）の詳細な定式化には弱く依存することが判明した。特性は、プロト中性子星の半径ではなく、$r_\nabla$ の単位で衝撃波半径を測定することで最も適切に捉えられる。

意義
本論文は、回転の効果を明示的に含む、SASI固有周波数の初の解析的近似を提供すると主張している。 「自己強制振動子」の枠組みを回転流へと拡張することで、差動回転による順行モードの成長率の線形的な増幅という、これまで数値計算では観察されていたものの解析的には理解されていなかった現象に対し、物理的な説明を与えている。導出された公式は、観測されたニュートリノおよび重力波の周波数から、恒星の特性（衝撃波半径や比角運動量など）を復元できる可能性を示唆している。また、本研究は、回転する前駆星におけるSASI周波数を正確に解釈する上で、完全な極方向構造の重要性を実証することにより、従来の円筒座標系や赤道近似の限界を浮き彫りにしている。