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

Resumen Técnico: Perspectiva Analítica de la Física del SASI II. Inestabilidad Espiral del Modo Progrado en un Núcleo Estelar Rotante

Planteamiento del Problema
Durante el colapso del núcleo de una estrella masiva en rotación, la Inestabilidad de Choque de Acreción Estacionaria (SASI, por sus siglas en inglés) impulsa movimientos no axisimétricos que imprimen firmas específicas en las señales de neutrinos y de ondas gravitacionales. Aunque estudios previos establecieron que la rotación favorece el crecimiento de los modos espirales progrados, el mecanismo físico detrás de esta desestabilización permanecía sin explicación, particularmente en regímenes donde las tasas de rotación son más lentas que la frecuencia del SASI. Además, los análisis perturbativos existentes a menudo dependían de aproximaciones cilíndricas o ecuatoriales esféricas que omitían la estructura poloidal del flujo, limitando su precisión para interpretar simulaciones numéricas de progenitores rotantes. Asimismo, el impacto de la rotación en las frecuencias de oscilación de los modos dominantes del SASI carecía de un marco analítico robusto.

Metodología
Este estudio extiende el análisis perturbativo de choques de acreción estancados en geometría esférica para incluir los efectos de la rotación y la rotación diferencial. El enfoque consiste en:

1. Modelado de Flujo Estacionario: Se construye un modelo simplificado del plano ecuatorial, asumiendo un gas perfecto con un índice adiabático $\gamma = 4/3$, dominado por la presión de degeneración electrónica relativista. El flujo incluye una función de enfriamiento dependiente de la densidad y la presión, y el potencial gravitatorio se trata en la aproximación newtoniana.
2. Análisis de Perturbación Lineal: Las perturbaciones lineales se caracterizan mediante armónicos esféricos $Y^m_\ell$ con número poloidal $\ell$ y número azimutal $m$. El análisis utiliza un conjunto de ecuaciones diferenciales acopladas para el flujo de masa, la entropía, el momento angular específico y la vorticidad.
3. Propiedades de Conservación: Para superar la dificultad de resolver la estructura poloidal en un progenitor rotante, el estudio utiliza una propiedad de conservación de la vorticidad radial. Esto permite reducir el sistema y resolverlo numéricamente mediante un método de disparo (shooting method).
4. Aproximación Adiabática: Se desarrolla un modelo adiabático para aislar la contribución de los procesos adiabáticos a la desestabilización. Esto reduce el orden del sistema diferencial y facilita una formulación de ecuación integral, tratando la región acústica post-choque como un oscilador auto-forzado.
5. Aproximación Analítica: El estudio deriva expresiones analíticas explícitas para las frecuencias propias de los modos dominantes del SASI ($\ell=m=1$ y $\ell=m=2$) mediante el ajuste de resultados numéricos a través de un rango de radios de choque y tasas de rotación.

Contribuciones Clave y Resultados

• Fórmulas de Frecuencia Propia Analíticas: El artículo proporciona aproximaciones analíticas para la parte real de las frecuencias propias ($\omega_r$) de los modos espirales dominantes ($\ell=m=1$ y $\ell=m=2$) con una precisión aproximada del 10% cuando el radio del choque $r_{sh}$ excede $\sim 1.5$ veces el radio de máxima deceleración ($r_\nabla$). La frecuencia se expresa como una función del radio de choque y el momento angular específico $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)$$
  donde $F$ y $G$ son funciones adimensionales derivadas de ajustes numéricos.
• Mecanismo de Desestabilización: El estudio explica que la desestabilización de los modos progrados por la rotación diferencial es impulsada por la mezcla de fases (phase mixing) entre el forzamiento advectivo y la estructura acústica. La rotación diferencial aumenta la longitud de onda radial de las perturbaciones de vorticidad asociadas al modo progrado. Esto conduce a un mejor ajuste de fase con la gran escala radial de la estructura acústica, aumentando así la eficiencia del acoplamiento.
• Papel de la Tasa de Rotación:
  • En regímenes de rotación modesta (donde no existe un radio de corotación), el efecto adverso sobre la mezcla de fases es significativo a radios pequeños debido al aumento abrupto hacia el interior de la frecuencia de rotación ($\propto 1/r^2$).
  • En regímenes de rotación más fuerte que involucran un radio de corotación, la aproximación de fase estacionaria identifica que el acoplamiento advectivo-acústico dominante ocurre entre la zona de corotación y el choque.
• Dominancia de Modos: El análisis confirma que la tasa de crecimiento del modo progrado aumenta aproximadamente de forma lineal con el momento angular específico. A medida que la rotación aumenta, el modo fundamental $\ell=m=2$ puede volverse dominante sobre el $\ell=m=1$. La relación de sus frecuencias depende de la tasa de rotación, pudiendo aumentar hasta un factor de 2.3 con la rotación en comparación con el caso no rotante.
• Independencia de los Detalles de Enfriamiento: Se encuentra que la frecuencia de oscilación del modo dominante depende débilmente de la formulación detallada de la emisión de neutrinos (parámetros de la función de enfriamiento), siempre que el radio del choque sea suficientemente grande ($r_{sh} \gtrsim 1.5 r_\nabla$). Las propiedades se capturan mejor cuando se mide el radio del choque en unidades de $r_\nabla$ en lugar del radio de la proto-estrella de neutrones.

Significancia
El artículo afirma proporcionar las primeras aproximaciones analíticas para las frecuencias propias del SASI que incluyen explícitamente el efecto de la rotación. Al extender el marco del "oscilador auto-forzado" a flujos rotantes, ofrece una explicación física para el aumento lineal de las tasas de crecimiento de los modos progrados por la rotación diferencial, un fenómeno observado previamente de forma numérica pero no comprendido analíticamente. Las fórmulas derivadas permiten la recuperación potencial de las propiedades estelares (como el radio del choque y el momento angular específico) a partir de las frecuencias observadas de neutrinos y ondas gravitacionales. El estudio también destaca las limitaciones de las aproximaciones cilíndricas o ecuatoriales previas, demostrando la importancia de la estructura poloidal completa para interpretar con precisión las frecuencias del SASI en progenitores rotantes.