Ca ii 854.2 nm in an enhanced network region simulated with MURaM-ChE

Título: Ca II λ854.2 nm em uma região de rede aprimorada simulada com MURaM-ChE

Autores: P.A. Ondratschek, D. Przybylski, H.N. Smitha, R.H. Cameron, S.K. Solanki
Instituição: Instituto Max Planck de Pesquisa do Sistema Solar, Göttingen, Alemanha.

1. Problema e Contexto

A linha espectral do Cálcio II em 854,2 nm (Ca II λ854.2 nm) é uma ferramenta diagnóstica fundamental para estudar a cromosfera solar. Observações do Sol quieto (QS) mostram que a média espacial desta linha apresenta uma assimetria vermelha e um desvio para o vermelho (redshift) no centro da linha.

• O Desafio: Modelos numéricos anteriores falharam em reproduzir simultaneamente a largura da linha e a assimetria observada. A literatura sugere que o efeito de divisão isotópica (isotopic splitting) é crucial para explicar a assimetria, mas modelos anteriores ainda não conseguiam combinar isso com a dinâmica atmosférica correta para obter a largura da linha observada.
• Objetivo: Investigar quão bem uma simulação computada com a extensão cromosférica do código MURaM (MURaM-ChE) reproduz o perfil médio espacial da linha Ca II λ854.2 nm, determinando as contribuições relativas da divisão isotópica versus a dinâmica atmosférica.

2. Metodologia

Os autores realizaram uma modelagem avançada combinando simulações magnetohidrodinâmicas (MHD) com cálculos de transferência radiativa (RT).

• Simulação Atmosférica:

  • Utilizou-se o código MURaM-ChE (extensão cromosférica do MURaM), que trata o hidrogênio fora do equilíbrio termodinâmico (NE) e inclui perdas radiativas e aquecimento por UV extremo.
  • O domínio simula uma região de rede aprimorada (enhanced network) com um campo magnético bipolar de grande escala.
  • A simulação abrange 24 Mm x 24 Mm x 24 Mm, com resolução horizontal de ~23 km e vertical de 20 km.
  • Foram analisados 10 minutos de série temporal após o sistema atingir um estado estacionário.

• Modelagem de Transferência Radiativa (Forward Modeling):

  • Utilizou-se o código RH1.5D para sintetizar a linha espectral em 10 colunas verticais independentes (abordagem 1.5D).
  • Três abordagens de modelagem atômica foram comparadas:
    1. Modelo de Isótopos (IM): Considera explicitamente os 6 isótopos estáveis do cálcio (40, 42, 43, 44, 46, 48Ca) como átomos separados em NLTE.
    2. Modelo de Átomo Composto (CM): Uma aproximação onde um único átomo "composto" mimetiza a presença de todos os isótopos, somando perfis Voigt ponderados pela abundância.
    3. Modelo de Isótopo Único (SIM): Considera apenas o isótopo mais abundante (40Ca), ignorando a divisão isotópica.

• Comparação: Os perfis sintetizados foram comparados com o atlas do espectrógrafo de Fourier de Hamburgo (Hamburg FTS atlas) e analisados através de bisetores de linha (line bisectors).

3. Resultados Principais

• Compatibilidade com Observações:

  • O modelo MURaM-ChE produziu um perfil médio espacial e temporal que coincide bem com as observações do Sol quieto em termos de profundidade, largura e forma da linha.
  • A largura da linha observada só foi reproduzida porque a atmosfera simulada possui velocidades dinâmicas (rms) suficientemente altas. O desvio padrão das velocidades Doppler nas linhas individuais foi de 4,56 km/s, cerca de 4 vezes maior do que em simulações anteriores (Leenaarts et al. 2009), o que explica a largura da linha média.

• Papel da Divisão Isotópica:

  • Assimetria: A assimetria vermelha característica (formato de "C invertido" no bisetor) só foi reproduzida quando a divisão isotópica foi incluída (Modelos IM e CM). O modelo de isótopo único (SIM) resultou em uma linha quase simétrica, falhando em capturar a assimetria observada.
  • Modelo Composto vs. Completo: O modelo de átomo composto (CM) foi uma boa aproximação do modelo completo de isótopos (IM) para a largura geral e assimetria, embora mostrasse pequenas diferenças no núcleo da linha e na intensidade exata, indicando que a suposição de densidades populacionais constantes entre isótopos não é perfeita em toda a cromosfera.

• Desvio para o Vermelho (Redshift):

  • O centro da linha média apresentou um leve desvio para o vermelho. Isso foi atribuído a um desequilíbrio de fluxos verticais na altura de formação do centro da linha: há uma cobertura de área maior de downflows (57%) em comparação com upflows (43%), resultando em um fluxo líquido descendente na altura de formação, embora isso não implique um fluxo de massa médio descendente global.

• Dependência das Condições Atmosféricas:

  • A análise por regiões (campos magnéticos fortes vs. fracos, upflows vs. downflows) mostrou que a forma do perfil e a assimetria dependem fortemente do local de formação.
  • Em regiões de campo magnético forte, mesmo o modelo SIM (sem isótopos) mostrou uma leve assimetria, sugerindo que a estrutura do fluxo acima das características magnéticas pode mimetizar parcialmente o efeito isotópico. No entanto, em regiões de campo fraco e downflows, a divisão isotópica é essencial para a assimetria correta.

4. Contribuições e Significância

1. Validação do MURaM-ChE: O estudo demonstra que o código MURaM-ChE é capaz de gerar modelos atmosféricos dinâmicos que reproduzem com alta fidelidade as propriedades das linhas cromosféricas observadas, superando limitações de modelos anteriores que subestimavam a largura da linha.
2. Confirmação da Necessidade de Isótopos: Reafirma-se que a inclusão da divisão isotópica do cálcio é obrigatória para modelar corretamente a assimetria da linha Ca II λ854.2 nm. Ignorar esse efeito leva a erros na interpretação física.
3. Implicações para Inversões: O trabalho alerta que, se os efeitos de múltiplos isótopos não forem considerados durante a inversão de perfis de linha (para inferir campos magnéticos e velocidades), os códigos de inversão tentarão compensar a falta de assimetria ajustando erroneamente as estruturas de velocidade e campo magnético, introduzindo viés nos resultados.
4. Dinâmica Cromosférica: Os resultados sugerem que a cromosfera do Sol quieto possui movimentos dinâmicos mais intensos do que o previsto em modelos anteriores, o que é crucial para explicar o alargamento das linhas espectrais.

Em resumo, o artigo estabelece que a combinação de uma atmosfera suficientemente dinâmica e a inclusão rigorosa do efeito de divisão isotópica são os dois pilares necessários para modelar com sucesso a linha Ca II λ854.2 nm e reproduzir as observações do Sol quieto.