Non-linear Dynamical Stability of Magnetic Polytropes

Resumen Técnico: Estabilidad Dinámica No Lineal de Politropos Magnéticos

Planteamiento del Problema
Si bien la estabilidad dinámica lineal de los equilibrios de gas ideal está bien establecida (estable para $\gamma > 4/3$), los límites de estabilidad no lineal para sistemas sujetos a perturbaciones de amplitud finita siguen siendo menos comprendidos. Esto es particularmente relevante para las estrellas masivas evolucionadas, que a menudo operan cerca de sus límites de estabilidad lineal. Los modelos existentes para los campos magnéticos estelares internos suelen sacrificar el rigor en favor de la tratabilidad o dependen de supuestos aximétricos que pueden no capturar la complejidad total de los campos turbulentos y desordenados. Además, los modelos simplificados frecuentemente fallan al aplicar de manera autoconsistente la ecuación de movimiento radial promediada esféricamente, lo que conduce a posibles inexactitudes en el tratamiento de la fuerza de Lorentz. Este trabajo investiga la estabilidad dinámica no lineal de los politropos de gas ideal bajo flujo homólogo, incorporando un modelo de campo medio autoconsistente para los campos magnéticos internos para determinar las condiciones bajo las cuales perturbaciones radiales pequeñas pero finitas pueden desvincular una estrella.

Metodología
El estudio emplea un modelo de campo medio de orden reducido, unidimensional, derivado de las ecuaciones de magnetohidrodinámica (MHD) ideal para un gas autogravitante y adiabático. La derivación se basa en dos simplificaciones primarias:

1. La aproximación de Cowling: Negligir la perturbación al potencial gravitatorio.
2. Flujo radial dominante: Asumir que el campo de velocidad es primordialmente radial y promediado esféricamente, ignorando efectivamente las fluctuaciones de velocidad no radiales.

El núcleo de la metodología consiste en derivar la ecuación de movimiento radial promediada esféricamente. A diferencia de los supuestos estándar donde la fuerza de Lorentz radial se trata simplemente como un gradiente de presión ($-\partial P_B/\partial r$), este trabajo deriva la forma exacta de la densidad de la fuerza de Lorentz radial promediada esféricamente. Esto conduce a un modelo de campo medio donde el campo magnético se caracteriza por la relación entre la energía magnética radial promedio y la energía magnética total, $\chi(r,t)$. El modelo asume un flujo homólogo ($r = r_0 a(t)$), lo que permite la integración de las ecuaciones de continuidad, inducción y Poisson para obtener las ecuaciones de evolución temporal para el factor de escala $a(t)$ y la distribución de la energía magnética.

Contribuciones Clave
El artículo introduce varios resultados teóricos novedosos respecto a la interacción de los campos magnéticos y la dinámica estelar:

• Formulación Corregida de la Fuerza de Lorentz: El trabajo demuestra que para un campo magnético isotrópico, la fuerza de Lorentz radial promedio es $-\frac{1}{3} \frac{dP_B}{dr}$, en lugar del típicamente asumido $-\frac{dP_B}{dr}$.
• Tensión Magnética e Isotropía: Establece que un pico central en el campo magnético requiere isotropía ($\chi = 1/3$) en el origen; todas las demás configuraciones deben desvanecerse en el origen debido a la tensión magnética.
• Soluciones de Entalpía Armónica: Para un $\gamma$ general, el análisis identifica una clase de soluciones con un perfil de entalpía espacialmente armónico. Estas soluciones difieren de los modelos estándar de Lane-Emden y permiten pulsaciones radiales no lineales.
• Umbrales de Estabilidad No Lineal: El artículo deriva condiciones analíticas para la desvinculación de un politropo. Encuentra que la presión de radiación puede desestabilizar un politropo linealmente estable ($\gamma > 4/3$) si la sobrepresión $\delta = \Delta P_0 / P_{eq}$ excede un umbral específico dependiente de $\gamma$, la relación de presión de radiación a material $\alpha_0$, y la distribución de la energía magnética.

Resultados
El análisis arroja dos clases amplias de soluciones homólogas:

1. Soluciones con $\gamma = 4/3$: Estas se reducen a perfiles espaciales similares a los de Lane-Emden, pero exhiben un comportamiento temporal restringido ya sea al colapso (para perturbaciones hacia adentro) o a la fuga al infinito (para perturbaciones hacia afuera). No existen pulsaciones estables para este índice.
2. Soluciones de Entalpía Armónica ($\gamma \neq 4/3$): Estas soluciones permiten pulsaciones radiales no lineales persistentes. El paisaje de estabilidad se mapea en el plano $(\gamma, \delta)$.
   • Para $\gamma > 4/3$, el sistema es linealmente estable. La fuga no lineal requiere una sobrepresión $\delta > \delta_{esc}$, donde $\delta_{esc} = \frac{3\gamma - 4}{1 + 3(\gamma - 1)\alpha_0}$.
   • A medida que la presión de radiación ($\alpha_0$) aumenta, el umbral de fuga $\delta_{esc}$ disminuye, haciendo que la estrella sea más susceptible a la desvinculación por perturbaciones finitas.
   • Para un politropo $n=3$ totalmente ionizado con $2/3$ de su energía magnética en el componente radial, el umbral de fuga se aproxima como $\delta \gtrsim 0.15 \mu^{-1} (300 M_\odot / M)^{1/2}$.

Adicionalmente, el artículo deriva relaciones masa-radio y campos magnéticos máximos para materia degenerada (enanas blancas) bajo este formalismo, encontrando masas máximas super-Chandrasekhar y campos pico que escalan como $\sim 10^{12} (M/M_\odot)^2$ G.

Significancia y Reivindicaciones
El artículo postula que el umbral de inestabilidad dinámica no lineal derivado puede tener relevancia para la pérdida de masa y la dispersión de estrellas masivas evolucionadas. Específicamente, sugiere que las perturbaciones radiales del orden del 1–10% podrían desvincular estrellas dominadas por la presión de radiación. Los autores señalan que, si bien el modelo depende de supuestos idealizados (flujo homólogo, perfil de entalpía armónica), las estimaciones derivadas de las amplitudes de inestabilidad pueden aplicarse de manera más general. El trabajo también especula que tales perturbaciones podrían ser "natales" (derivadas de la formación estelar) y persistir debido a los largos tiempos de difusión de la radiación en estrellas de baja densidad y gran tamaño.

Los autores enfatizan que la utilidad principal de este trabajo es servir como un modelo idealizado para cuantificar los límites de estabilidad no lineal. Reconocen que la aplicación directa a estrellas reales está limitada por los supuestos idealizados, pero argumentan que el rigor del modelo magnético —específicamente el tratamiento autoconsistente de la fuerza de Lorentz promediada esféricamente— ofrece una mejora sobre los enfoques analíticos existentes y proporciona un marco para guiar futuras búsquedas numéricas de equilibrios magnéticos estables.