Revealing mixed modes in compressible hydrodynamical simulations of red giant stars

Resumen Técnico: Revelando los Modos Mixtos en Simulaciones Hidrodinámicas Compresibles de Gigantes Rojas

Planteamiento del Problema
Los modos mixtos en estrellas evolucionadas de baja masa proporcionan diagnósticos críticos para las tasas de rotación del núcleo, las cuales son observadas como significativamente menores de lo predicho por los modelos estándar de evolución estelar. Se ha propuesto que los modos mixtos son un mecanismo para el transporte de momento angular (AM), pero cuantificar su eficiencia requiere un conocimiento preciso de sus amplitudes. Las predicciones empíricas actuales de estas amplitudes dependen de formalismos semiempíricos que están pobremente restringidos por las observaciones para las gigantes rojas evolucionadas. Específicamente, sigue siendo incierto si los modos mixtos pueden transportar suficiente AM para reconciliar las tasas de rotación del núcleo observadas con los modelos teóricos.

Metodología
Este estudio presenta las primeras simulaciones hidrodinámicas compresibles en 2D de modos mixtos en el interior de una estrella de baja masa evolucionada utilizando el código music. Los autores simularon una cuña meridional de una gigante roja de 1.3 $M_\odot$ (basada en el modelo M1 de Belkacem et al. 2015b). Se realizaron dos simulaciones distintas para investigar el impacto del límite exterior en las amplitudes de los modos:

• RGext90: Truncamiento radial exterior en $r_o/r_\star = 0.90$.
• RGext98: Truncación radial exterior en $r_o/r_\star = 0.98$.

Ambos modelos incluyeron una gran región convectiva para excitar estocásticamente los modos y una región sustancial estratificada de forma estable para albergar la componente de modo de gravedad de los modos mixtos. Las simulaciones se ejecutaron hasta alcanzar un estado estacionario dinámico, tras lo cual se capturaron los espectros.

Para validar las simulaciones, los autores compararon las frecuencias y funciones propias recuperadas contra las predicciones de la teoría lineal generadas por dos resolvedores independientes: el código gyre (versión 8.0.1) y un resolvedor de autovalores personalizado escrito en dedalus. El análisis se centró en los espectros de velocidad radial, construyendo diagramas de Échelle y gráficos de frecuencia-radio para grados armónicos $\ell < 5$.

Resultados Clave

1. Acuerdo de Frecuencia y Función Propia:

   • Se encontró un excelente acuerdo entre las simulaciones y la teoría lineal para todos los modos dominados por p.
   • Se observaron discrepancias menores (unos pocos $\mu$Hz) para los modos dominados por g, particularmente en el rango de frecuencia $[60, 240]$ $\mu$Hz.
   • El acuerdo de la función propia fue excelente para todos los modos, excepto para aquellos en el rango de frecuencia mencionado. Los autores atribuyen las discrepancias en este rango a una resolución insuficiente de la parte del modo g en la simulación en comparación con los cálculos de teoría lineal de alto orden.

2. Distribución de Energía Cinética:

   • Contrario a las predicciones empíricas que sitúan la energía cinética máxima alrededor de $\nu_{max} \approx 312.8$ $\mu$Hz, las simulaciones revelaron que las mayores energías cinéticas se localizan en frecuencias bajas ($\nu < 50$ $\mu$Hz).
   • Estos modos de baja frecuencia exhiben energías cinéticas órdenes de magnitud mayores que las predicciones empíricas, con un rápido descenso de la energía a medida que la frecuencia aumenta.
   • La simulación RGext98 (que se extiende más cerca de la superficie) mostró energías cinéticas totales más altas y la aparición de una distribución secundaria de energía, de tipo Gaussiano ancho, centrada alrededor de $\nu = 700$ $\mu$Hz, probablemente debido a la excitación superficial de los modos p.

3. Velocidades Superficiales Extrapoladas:

   • A pesar de sus altas energías cinéticas internas, los modos de baja frecuencia ($\nu < 50$ $\mu$Hz) se extrapolan a velocidades superficiales muy pequeñas (mayormente $< 1$ cm s$^{-1}$), lo que los hace difíciles de observar.
   • La simulación RGext98 produjo velocidades superficiales comparables a las predicciones empíricas para los modos de alta frecuencia ($500-800$ $\mu$Hz), alcanzando un pico en una curva con forma de campana. Sin embargo, el pico está desplazado hacia frecuencias más altas en comparación con la escala empírica.
   • La simulación RGext90 produjo velocidades superficiales insignificantes para los modos de alta frecuencia.

Significancia y Reivindicaciones
El artículo afirma proporcionar las primeras restricciones directas sobre las amplitudes de los modos mixtos en simulaciones hidrodinámicas en 2D de gigantes rojas. La principal significancia radica en el descubrimiento de modos mixtos de baja frecuencia altamente energéticos que no son capturados por las actuales relaciones de escala empíricas.

Los autores postulan que, si bien estos modos de baja frecuencia son difíciles de detectar mediante velocidades superficiales, sus grandes energías cinéticas en lo profundo del interior estelar podrían impactar significativamente el transporte de momento angular. Esto sugiere que los modos mixtos pueden desempeñar un papel más sustancial en el transporte de AM de lo estimado previamente, lo que podría explicar las tasas de rotación del núcleo observadas en las gigantes rojas jóvenes y muy evolucionadas. El estudio demuestra la viabilidad de utilizar simulaciones hidrodinámicas compresibles para capturar modos mixtos, ofreciendo una vía para estimar cuantitativamente la eficiencia del transporte de AM y predecir la observabilidad para futuras misiones como PLATO, siempre que eventualmente se incorporen los efectos 3D, la rotación y los campos magnéticos.