Tearing Instability in Gyrotropic MHD: Effects of Equilibrium Pressure Anisotropy

Resumen Técnico: Inestabilidad de Desgarro en MHD Girotrópica: Efectos de la Anisotropía de Presión del Equilibrio

Planteamiento del Problema
La inestabilidad de desgarro (tearing instability) es un mecanismo fundamental para la reconexión magnética en láminas de corriente delgadas. La teoría clásica de la Magnetohidrodinámica (MHD) resistiva, establecida por Furth et al. (1963) y extendida por Coppi et al. (1976), asume un tensor de presión de equilibrio isotrópico. En este marco, los escalamientos de orden principal para la tasa de crecimiento y el número de onda del desgarro son independientes de la $\beta$ del plasma (la relación entre la presión térmica y la presión magnética). Sin embargo, muchos plasmas astrofísicos y espaciales son débilmente colisionales, lo que permite mantener una anisotropía de presión significativa ($p_\parallel \neq p_\perp$). Aunque estudios cinéticos y de fluidos han sugerido que la anisotropía influye en la estructura del modo y en las tasas de crecimiento, carecía de un marco analítico sistemático que conectara la anisotropía de equilibrio prescrita con la construcción estándar de la asíntota exterior-interior de la teoría de desgarro. Específicamente, no se habían derivado rigurosamente los efectos separados de la anisotropía en el parámetro de acoplamiento ideal exterior ($\Delta'$) y en el balance de fuerzas de la capa interna resistiva.

Metodología
Los autores desarrollan una teoría lineal de desgarro resistivo dentro de un marco de MHD girotrópica no ideal. El estudio se centra en una lámina de corriente de tipo Harris libre de fuerzas, caracterizada por presiones paralelas y perpendiculares de equilibrio uniformes, definidas por la $\beta$ del plasma ($\beta_0$) y la diferencia de anisotropía ($\Delta\beta_0 = \beta_{\parallel,0} - \beta_{\perp,0}$).

1. Ecuaciones de Gobierno: Los autores utilizan ecuaciones de MHD girotrópicas incompresibles y no ideales, cerradas mediante relaciones politrópicas para las presiones paralelas y perpendiculares ($\gamma_\parallel$ y $\gamma_\perp$). Esto incluye términos de calentamiento óhmico y viscoso para mantener el equilibrio energético, siguiendo la formulación de Hau (2002).
2. Linealización: Se aplican perturbaciones pequeñas a un perfil de equilibrio "congelado". El sistema se linealiza para derivar un conjunto de ecuaciones acopladas para la velocidad, el campo magnético y las perturbaciones de la anisotropía de presión.
3. Análisis Asintótico: Se emplea la construcción estándar de capa exterior-interior:
   • Región Exterior: Una región ideal donde se desprecia la resistividad. Los autores derivan una ecuación diferencial de segundo orden para la perturbación del flujo magnético, identificando cómo la anisotropía modifica la tasa de decaimiento de campo lejano y el índice de estabilidad de desgarro $\Delta'$.
   • Región Interior: Una capa resistiva donde el balance de momento está dominado por la tensión anisotrópica.
   • Acoplamiento (Matching): Las ramas FKR (constante-$\psi$) y Coppi (no constante-$\psi$) se derivan y se acoplan para determinar el modo de crecimiento más rápido.
4. Validación Numérica: Las predicciones analíticas se prueban frente a cálculos numéricos de valores propios utilizando el marco PSECAS (Calculadora de Valores Propios Pseudo-Espectral con un Solucionador Automatizado). Se examinan tres cierres específicos: doble adiabático ($\gamma_\parallel=3, \gamma_\perp=2$), doble politrópico ($\gamma_\parallel=0.5, \gamma_\perp=2$) y doble isotérmico ($\gamma_\parallel=1, \gamma_\perp=1$).

Contribuciones Clave y Resultados

• Dinámica Modificada de la Región Exterior: Los autores derivan que la anisotropía de presión de equilibrio altera la tasa de decaimiento de campo lejano del modo de desgarro. Esto introduce un nuevo corte de estabilidad marginal para modos localizados: $\alpha_c = \sqrt{A/R_0}$, donde $A = 1 - \Delta\beta_0/2$ representa la tensión anisotrópica efectiva y $R_0$ incorpora la respuesta de la presión. Los modos con números de onda $\alpha \geq \alpha_c$ se estabilizan contra el desgarro localizado.
• Dinámica Modificada de la Capa Interior: A diferencia del caso isotrópico, donde la anisotropía afecta principalmente al acoplamiento exterior, la anisotropía de equilibrio prescrita entra en el balance de momento de orden principal en la capa interior a través del factor $A$. Esto modifica directamente el balance de fuerzas que gobierna tanto la rama FKR como la Coppi.
• Leyes de Escalamiento Revisadas: Al acoplar las ramas modificadas de FKR y Coppi, los autores derivan las leyes de escalamiento para el modo de crecimiento más rápido. Aunque se mantienen los exponentes clásicos del número de Lundquist ($S$) ($\gamma_{max} \sim S^{-1/2}$, $k_{max} \sim S^{-1/4}$), los prefactores ya no son universales. Dependen de la anisotropía de equilibrio y del cierre:
  • $\gamma_{max}\tau_A \sim A^{1/2} R_0^{-1/4} S^{-1/2}$
  • $k_{max}a \sim A^{1/4} R_0^{-3/8} S^{-1/4}$
• Confirmación Numérica: Los cálculos de PSECAS confirman el escalamiento de la rama Coppi y la rama FKR cuando se utiliza una aproximación de longitud de onda finita ajustada para $\Delta'$. Los resultados numéricos validan la predicción analítica de que una $\Delta\beta_0$ negativa (dominancia de la presión perpendicular) potencia las tasas de crecimiento y desplaza el modo más rápido hacia números de onda más altos, mientras que una $\Delta\beta_0$ positiva suprime el desgarro y ensancha la capa interna.
• Dependencia del Cierre: El estudio destaca que la magnitud del efecto de la anisotropía depende del cierre politrópico. En el límite doble isotérmico ($\gamma_\parallel=\gamma_\perp=1$), la respuesta de la presión dependiente de la $\beta$ del plasma desaparece, aislando el efecto directo de la anisotropía de equilibrio a través de $A$ y $R_0$.

Significancia y Reivindicaciones
El artículo afirma proporcionar una "extensión controlada" de la teoría clásica de desgarro hacia el régimen de MHD girotrópica. Su principal significancia radica en demostrar que una anisotropía de presión de equilibrio prescrita altera fundamentalmente el inicio lineal del desgarro al modificar tanto las condiciones de acoplamiento ideal exterior como la dinámica de la capa interna resistiva.

Los autores declaran explícitamente que su trabajo es un modelo teórico destinado a clarificar el régimen de MHD girotrópica resistiva lineal, en lugar de proporcionar un modelo cinético directo y completo de la reconexión astrofísica. Señalan que los resultados son condicionales: muestran cómo una anisotropía prescrita modifica el problema de estabilidad, pero no contabilizan la relajación cinética, la física de Hall o la evolución autoconsistente de la anisotropía. Por consiguiente, el artículo argumenta que en entornos de alta $\beta$ o débilmente colisionales, asumir la isotropía de presión puede conducir a errores de orden unidad en la estimación de las tasas de crecimiento de desgarro, las longitudes de onda preferidas y los espesores de la capa interna. El trabajo sirve como una base para futuras extensiones que incorporen compresibilidad, perfiles de presión autoconsistentes y efectos de dos fluidos.