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

Resumo Técnico: Redução da Viscosidade Efetiva por Transporte Anisotrópico e Instabilidades de Plasma nos Núcleos de Aglomerados de Galáxias em Oscilação (Sloshing)

Definição do Problema
O meio intrac de aglomerados de galáxias (ICM) é um plasma quente e fracamente magnetizado onde o caminho livre médio de Coulomb ($\lambda_{mfp}$) é significativamente menor que o tamanho do sistema, mas maior que o raio de Larmor dos íons ($\rho_L$). Embora o ICM seja frequentemente modelado como um fluido, suas propriedades de transporte microfísicas — especificamente a viscosidade — são altamente sensíveis à colisionalidade e à geometria do campo magnético. A viscosidade isotrópica padrão de Spitzer prevê um número de Reynolds ($Re$) que poderia ser baixo o suficiente para suprimir a turbulência, contudo, evidências observacionais de frentes frias (CFs) e dispersões de velocidade sugerem que a viscosidade efetiva é muito menor (1–20% de Spitzer).

Considerações teóricas indicam que, em plasmas de alto-$\beta$, o campo magnético induz tensão viscosa anisotrópica (viscosidade de Braginskii). Além disso, a anisotropia de pressão ($\Delta p = p_\perp - p_\parallel$) impulsionada pela invariância adiabática pode desencadear instabilidades de plasma cinética (firehose, mirror, ciclo-trônico de íons). Essas instabilidades espalham as partículas, efetivamente limitando a anisotropia de pressão a valores marginalmente estáveis e reduzindo a tensão viscosa efetiva. A interação entre a viscosidade de Braginskii, a geometria do campo magnético e os limitadores impulsionados por instabilidades permanece mal restringida no contexto dos movimentos de gás em oscilação (sloshing) observados em aglomerados de núcleo frio.

Metodologia
Os autores realizaram simulações de magnetohidrodinâmica (MHD) paralela de alta resolução usando o código FLASH 4.7.1 para modelar um aglomerado de galáxia de núcleo frio ($M \sim 10^{15} M_\odot$, $T \sim 7$ keV) passando por uma fusão com um subaglomerado. A configuração da simulação inclui:

• Física: As simulações resolvem as equações conservativas de MHD incluindo o tensor de tensão viscosa de Braginskii. A anisotropia de pressão é modelada via uma abordagem de limitador de "parede dura" (hard-wall), onde $\Delta p$ é limitado por limiares de estabilidade marginal derivados da teoria cinética.
• Prescrições de Limite: Três cenários foram testados:
  1. Ilimitado: Sem limites em $\Delta p$ (a viscosidade de Braginskii atua livremente, embora as equações sejam mal postas em escalas cinéticas).
  2. Limitadores v1: $\Delta p$ é restrito pela instabilidade firehose oblíqua ($\delta p \approx -1.4/\beta$) e mirror ($\delta p \approx 1/\beta$).
  3. Limitadores v2: $\Delta p$ é restrito pelos limiares de firehose paralela ($\delta p \approx -2/\beta$) e ciclo-trônico de íons ($\delta p \approx 0.5/\beta^{0.5}$).
• Condições Iniciais: Dois conjuntos de simulações foram executados com diferentes valores iniciais de beta de plasma ($\beta_{ini} = 100$ e $400$) para explorar o efeito de variações na tensão magnética.
• Análise: O estudo analisa a supressão de instabilidades de Kelvin-Helmholtz (KHI) em frentes frias, a distribuição espacial da anisotropia de pressão, fatores de redução de viscosidade efetiva, espectros de amplitude de velocidade (3D e projetados em 2D) e funções de transferência de energia.

Principais Contribuições e Resultados

1. Supressão de Instabilidades de Kelvin-Helmholtz (KHI):
   A viscosidade de Braginskii suprime a KHI nas superfícies das frentes frias, mas o efeito é modesto comparado à viscosidade isotrópica. A supressão depende fortemente da força do campo magnético e da prescrição do limitador. Em simulações com campos magnéticos iniciais mais fortes ($\beta_{ini}=100$), a tensão magnética sozinha fornece supressão significativa de KHI, tornando a adição de viscosidade menos visualmente distinta. Em casos de campos mais fracos ($\beta_{ini}=400$), a viscosidade desempenha um papel maior, embora os limitadores de "parede dura" reduzam a viscosidade efetiva, resultando em menos supressão do que o caso "Ilimitado".

2. Redução da Viscosidade Efetiva:
   O efeito combinado da anisotropia do campo magnético e dos limitadores de instabilidade de plasma reduz drasticamente a viscosidade efetiva.

• Supressão Geométrica: A geometria emaranhada do campo magnético sozinha reduz o fluxo viscoso em relação ao caso isotrópico (fator $f_{geom}$).
• Supressão por Instabilidade: Os limitadores reduzem ainda mais o fluxo (fator $f_{limit}$).
• Resultado Líquido: No núcleo do aglomerado, a viscosidade efetiva é reduzida para $\sim 20\%$ do valor de Spitzer em 57–75% do volume, e para $\sim 50\%$ em 90–96% do volume. Essa redução é mais pronunciada em regiões onde o campo magnético é mais fraco e ainda não foi amplificado pelo cisalhamento (shear).

3. Anisotropia de Pressão e Instabilidades:
   Os movimentos de oscilação (sloshing) geram anisotropias de pressão significativas ($|\delta p| \sim 0.03-0.07$ para $\beta \sim 100-1000$). Em simulações "Ilimitadas", uma fração do volume (até 30% para os limiares mirror/firehose oblíqua) excede os limites de instabilidade. Quando os limitadores são aplicados, a anisotropia de pressão é contida pela estabilidade marginal, com uma fração não negligenciável do volume situando-se exatamente nos limiares mirror ou firehose. Isso tem implicações para o transporte de raios cósmicos via espalhamento de ângulo de pitch.

4. Espectros de Velocidade e Turbulência:

• Inclinação Espectral: Em simulações sem viscosidade (inviscid), o espectro de amplitude de velocidade 3D segue uma inclinação do tipo Kolmogorov ($k^{-1/3}$). A inclusão da viscosidade de Brag']] sugere um endurecimento do espectro para $\sim k^{-2/3}$ a $k^{-1}$.
• Interpretação: Os autores sugerem que este endurecimento pode resultar do amortecimento de movimentos turbulentos pela viscosidade, ou pode refletir os movimentos de oscilação de grande escala subjacentes que ainda não decaíram totalmente em turbulência. O endurecimento é menos pronunciado quando os limitadores são aplicados, pois a viscosidade efetiva é menor.
• Observáveis Projetados: Mapas projetados de velocidade bulk e dispersão de velocidade mostram que características de pequena escala são suavizadas em execuções viscosas em comparação com as sem viscosidade, embora as distribuições de grande escala permaneçam semelhantes.

5. Magneto-Imutabilidade:
   As simulações fornecem evidências para a dinâmica "magneto-imutável". As funções de distribuição de probabilidade (PDFs) da taxa de deformação (strain) alinhada ao campo ($S$) são significativamente mais estreitas em simulações de Braginskii do que em execuções de MHD sem viscosidade. Isso indica um processo de auto-organização onde fluxos turbulentos rearranjam para minimizar a deformação alinhada ao campo magnético, suprimindo ainda mais a dissipação viscosa além do que a geometria sozinha preveria. No entanto, os autores observam que isso pode ser uma realização "parcial", já que as funções de transferência de energia mostram dissipação não negligenciável.

Significância e Alegações
O artigo afirma que, se o transporte de momento viscoso no ICM for controlado pela anisotropia de pressão e limitado por instabilidades de plasma, suas assinaturas observacionais serão sutis.

• Frentes Frias: É difícil restringir a viscosidade efetiva apenas pela suavidade das frentes frias, pois a tensão magnética e a viscosidade anisotrópica produzem supressão de KHI semelhante.
• Campos de Velocidade: A assinatura observacional mais promissora reside nas propriedades do campo de velocidade acessíveis via instrumentos de microcalorímetros (ex: XRISM). O endurecimento do espectro de amplitude de velocidade é um indicador potencial de viscosidade anisotrópica, embora a origem física da inclinação específica ($k^{-2/3}$) exija investigação adicional.
• Auto-organização: Os resultados apoiam a existência de dinâmicas magneto-imutáveis nos núcleos de aglomerados, o que aumenta a supressão da viscosidade e pode estender as cascata de turbulência para escalas menores do que anteriormente pensado.

Os autores concluem que, embora suas simulações capturem a física essencial do transporte anisotrópico e do limite de instabilidade, trabalhos futuros com maior resolução, turbulência impulsionada em caixas estratificadas e tratamentos cinéticos mais avançados (ex: CGL-MHD) são necessários para resolver totalmente a origem das inclinações espectrais e a extensão da magneto-imutabilidade.