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

技术摘要：关于 SASI II 物理机制的解析洞察：旋转恒星核中前向螺旋模的螺旋不稳定性

问题陈述
在旋转大质量恒星的核塌缩过程中，驻波吸积冲击波不稳定性（SASI）驱动非轴对称运动，并在中微子和引力波信号中留下特定的特征印记。虽然先前的研究已经确定旋转有利于前向螺旋模（prograde modes）的增长，但其背后的物理机制仍未得到解释，特别是在旋转速率慢于 SASI 频率的机制中。此外，现有的摄动分析通常依赖于圆柱或球赤道近似，忽略了流体的极向结构，这限制了其在解释旋转原恒星数值模拟方面的准确性。此外，旋转对主导 SASI 模振荡频率的影响缺乏稳健的解析框架。

研究方法
本研究扩展了在球对称几何下对停滞吸积冲击波的摄动分析，以包含旋转和微分旋转的影响。该方法包括：

1. 定常流建模： 构建了一个简化的赤道面模型，假设为理想气体，绝热指数 $\gamma = 4/3$，由相对论电子简并压力主导。流体包含依赖于密度和压力的冷却函数，引力势在牛顿近似下进行处理。
2. 线性摄动分析： 线性摄动通过具有极向数 $\ell$ 和方位数 $m$ 的球面谐函数 $Y^m_\ell$ 进行表征。分析利用了一组关于质量通量、熵、比角动量和涡度的耦合微分方程。
3. 守恒特性： 为了克服在旋转原恒星中求解极向结构的困难，本研究利用了径向涡度的守恒特性。这使得系统得以简化，并可以使用射击法（shooting method）进行数值求解。
4. 绝热近似： 开发了一个绝热模型以隔离绝热过程对失稳性的贡献。这降低了微分系统的阶数，并促进了积分方程形式的构建，将后冲击声学区域视为一个自强迫振荡器（self-forced oscillator）。
5. 解析近似： 通过在不同冲击半径和旋转速率范围内拟合数值结果，推导出了主导 SASI 模（$\ell=m=1$ 和 $\ell=m=2$）特征频率的显式解析表达式。

核心贡献与结果

• 解析特征频率公式： 本文提供了主导螺旋模（$\ell=m=1$ 和 $\ell=m=2$）特征频率实部（$\omega_r$）的解析近似，当冲击半径 $r_{sh}$ 大于约 $1.5$倍的最大减速半径$r_\nabla$时，其精度约为 10$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 频率方面的重要性，强调了以往圆柱或赤道近似的局限性。