Reduced Effective Viscosity from Anisotropic Transport and Plasma Instabilities in the Sloshing Cores of Galaxy Clusters

Resumen Técnico: Reducción de la Viscosidad Efectiva por Transporte Anisotrópico e Inestabilidades de Plasma en los Núcleos de Remoción (Sloshing) de Cúmulos de Galaxias

Planteamiento del Problema
El medio intracúmulo (ICM) de los cúmulos de galaxias es un plasma caliente y débilmente magnetizado donde el camino libre medio de Coulomb ($\lambda_{mfp}$) es significativamente menor que el tamaño del sistema pero mayor que el radio de Larmor de los iones ($\rho_L$). Aunque el ICM se modela a menudo como un fluido, sus propiedades de transporte microfísico —específicamente la viscosidad— son altamente sensibles a la colisionalidad y a la geometría del campo magnético. La viscosidad isotrópica estándar de Spitzer predice un número de Reynolds ($Re$) que podría ser lo suficientemente bajo como para suprimir la turbulencia, sin embargo, la evidencia observacional de frentes fríos (CFs) y dispersiones de velocidad sugiere que la viscosidad efectiva es mucho menor (1–20% de la de Spitzer).

Las consideraciones teóricas indican que, en plasmas de alto $\beta$, el campo magnético induce un esfuerzo viscoso anisotrópico (viscosidad de Braginskii). Además, la anisotropía de presión ($\Delta p = p_\perp - p_\parallel$) impulsada por la invariancia adiabática puede desencadenar inestabilidades de plasma cinéticas (firehose, mirror, ion-cyclotron). Estas inestabilidades dispersan las partículas, limitando efectivamente la anisotropía de presión a valores de estabilidad marginal y reduciendo el esfuerzo viscoso efectivo. La interacción entre la viscosidad de Braginskii, la geometría del campo magnético y los limitadores impulsados por estas inestabilidades permanece poco restringida en el contexto de los movimientos de gas por remoción (sloshing) observados en los núcleos fríos de los cúmulos.

Metodología
Los autores realizaron simulaciones de magnetohidrodinámica (MHD) paralela de alta resolución utilizando el código FLASH 4.7.1 para modelar un cúmulo de galaxias de núcleo frío ($M \sim 10^{15} M_\odot$, $T \sim 7$ keV) que experimenta una fusión con un subclúster. La configuración de la simulación incluye:

• Física: Las simulaciones resuelven las ecuaciones conservativas de MHD incluyendo el tensor de esfuerzo viscoso de Braginskii. La anisotropía de presión se modela mediante un enfoque de limitador de "pared dura" (hard-wall), donde $\Delta p$ se restringe a los umbrales de estabilidad marginal derivados de la teoría cinética.
• Prescripciones de Limitación: Se probaron tres escenarios:
  1. Sin límites (Unlimited): Sin límites en $\Delta p$ (la viscosidad de Braginskii actúa libremente, aunque las ecuaciones están mal planteadas en escalas cinéticas).
  2. Limitadores v1: $\Delta p$ es restringido por los umbrales de inestabilidad firehose oblicua ($\delta p \approx -1.4/\beta$) y mirror ($\delta p \approx 1/\beta$).
  3. Limitadores v2: $\Delta p$ es restringido por los umbrales de inestabilidad firehose paralela ($\delta p \approx -2/\beta$) e ion-cyclotron ($\delta p \approx 0.5/\beta^{0.5}$).
• Condiciones Iniciales: Se ejecutaron dos conjuntos de simulaciones con diferentes valores iniciales de beta de plasma ($\beta_{ini} = 100$ y $400$) para explorar el efecto de variar la tensión magnética.
• Análisis: El estudio analiza la supresión de las inestabilidades de Kelvin-Helmholtz (KHI) en los frentes fríos, la distribución espacial de la anisotropía de presión, los factores de reducción de la viscosidad efectiva, los espectros de amplitud de velocidad (3D y proyectados en 2D) y las funciones de transferencia de energía.

Contribuciones Clave y Resultados

1. Supresión de las Inestabilidades de Kelvin-Helmholtz (KHI):
   La viscosidad de Braginskii suprime la KHI en las superficies de los frentes fríos, pero el efecto es modesto comparado con la viscosidad isotrópica. La supresión depende fuertmente de la fuerza del campo magnético y de la prescripción de limitación. En simulaciones con campos magnéticos iniciales más fuertes ($\beta_{ini}=100$), la tensión magnética por sí sola proporciona una supresión significativa de la KHI, haciendo que la adición de viscosidad sea menos distintiva visualmente. En los casos de campos más débiles ($\beta_{ini}=400$), la viscosidad juega un papel mayor, aunque los limitadores de "pared dura" reducen la viscosidad efectiva, resultando en una menor supresión que en el caso "Sin límites".

2. Reducción de la Viscosidad Efectiva:
   El efecto combinado de la anisotropía del campo magnético y los limitadores de inestabilidad de plasma reduce drásticamente la viscosidad efectiva.

   • Supresión Geométrica: La geometría enredada del campo magnético por sí sola reduce el flujo viscoso respecto al caso isotrópico (factor $f_{geom}$).
   • Supresión por Inestabilidad: Los limitadores reducen aún más el flujo (factor $f_{limit}$).
   • Resultado Neto: En el núcleo del cúmulo, la viscosidad efectiva se reduce a $\sim 20\%$ del valor de Spitzer en el 57–75% del volumen, y a $\sim 50\%$ en el 90–96% del volumen. Esta reducción es más pronunciada en regiones donde el campo magnético es más débil y aún no ha sido amplificado por la cizalla (shear).

3. Anisotropía de Presión e Inestabilidades:
   Los movimientos de remoción (sloshing) generan anisotropías de presión significativas ($|\delta p| \sim 0.03-0.07$ para $\beta \sim 100-1000$). En las simulaciones "Sin límites", una fracción del volumen (hasta un 30% para los umbrales de inestabilidad mirror/firehose oblicua) excede los límites de inestabilidad. Cuando se aplican los limitadores, la anisotropía de presión se restringe a la estabilidad marginal, con una fracción no despreciable del volumen situándose exactamente en los umbrales de mirror o firehose. Esto tiene implicaciones para el transporte de rayos cósmicos mediante la dispersión de ángulo de pitch.

4. Espectros de Velocidad y Turbulencia:

   • Pendiente Espectral: En simulaciones sin viscosidad, el espectro de amplitud de velocidad 3D sigue una pendiente tipo Kolmogorov ($k^{-1/3}$). La inclusión de la viscosidad de Braginskii inclina el espectro hacia $k^{-2/3}$ a $k^{-1}$.
   • Interpretación: Los autores sugieren que este endurecimiento puede resultar del amortiguamiento de los movimientos turbulentos por la viscosidad, o puede reflejar los movimientos de remoción a gran escala subyacentes que no se han decaído completamente en turbulencia. El endurecimiento es menos pronunciado cuando se aplican los limitadores, ya que la viscosidad efectiva es menor.
   • Observables Proyectados: Los mapas de velocidad bulk proyectada y de dispersión de velocidad muestran que las características de pequeña escala se suavizan en las ejecuciones viscosas en comparación con las no viscosas, aunque las distribuciones a gran escala permanecen similares.

5. Magneto-Inmutabilidad:
   Las simulaciones proporcionan evidencia de la dinámica "magneto-inmutable". Las funciones de distribución de probabilidad (PDF) de la tasa de deformación ($S$) alineada con el campo son significativamente más estrechas en las simulaciones de Braginskii que en las ejecuciones de MHD sin viscosidad. Esto indica un proceso de autoorganización donde los flujos turbulentos reordenan el gas para minimizar la deformación alineada con el campo magnético, suprimiendo aún más la disipación viscosa de lo que la geometría sola predeciría. Sin embargo, los autores señalan que esto puede ser una realización "parcial", ya que las funciones de transferencia de energía muestran una disipación no despreciable.

Significancia y Reivindicaciones
El artículo sostiene que si el transporte de momento viscoso en el ICM está controlado por la anisotropía de presión y limitado por inestabilidades de plasma, sus firmas observacionales son sutiles.

• Frentes Fríos: Es difícil constreñir la viscosidad efectiva basándose únicamente en la suavidad de los frentes fríos, ya que la tensión magnética y la viscosidad anisotrópica producen una supresión de la KHI similar.
• Campos de Velocidad: La propiedad más prometedora de los campos de velocidad, accesible mediante instrumentos de microcalorimetría (ej. XRISM), es el endurecimiento del espectro de amplitud de velocidad, lo cual es un indicador potencial de la viscosidad anisotrópica, aunque el origen físico de la pendiente específica ($k^{-2/3}$) requiere mayor investigación.
• Autoorganización: Los resultados apoyan la existencia de dinámicas magneto-inmutables en los núcleos de los cúmulos, lo que potencia la supresión de la viscosidad y puede extender las cascadas turbulentas a escalas más pequeñas de lo que se pensaba anteriormente.

Los autores concluyen que, si bien sus simulaciones capturan la física esencial del transporte anisotrópico y la limitación por inestabilidad, es necesario un trabajo futuro con mayor resolución, turbulencia forzada en cajas estratificadas y tratamientos cinéticos más avanzados (ej. CGL-MHD) para resolver plenamente el origen de las pendientes espectrales y el alcance de la magneto-inmutabilidad.