The IRIS$^{2+}$ inversion tool: recovering the radiative losses and the thermodynamics in the lower solar atmosphere

Resumo Técnico: A Ferramenta de Inversão IRIS2+

Definição do Problema
Inferir o estado físico do plasma cromosférico a partir de observações continua sendo um desafio significativo na física solar. Devido às baixas taxas de colisão na cromosfera, é necessária uma descrição de equilíbrio termodinâmico não local (não-LTE) para explicar as intensidades emergentes, o que impede uma tradução direta das intensidades observadas para parâmetros físicos. Embora os códigos de inversão não-LTE modernos (ex: STiC) tenham melhorado a fidelidade ao incorporar múltiplas linhas espectrais para aumentar a resolução de profundidade, eles são computacionalmente intensivos. Além disso, calcular as perdas radiativas integradas (IRL) — cruciais para entender o balanço de energia e os mecanismos de aquecimento na cromosfera — requer um pós-processamento complexo dos modelos atmosféricos invertidos. Há uma necessidade de uma ferramenta rápida e confiável, capaz de recuperar parâmetros termodinâmicos e perdas radiativas para grandes conjuntos de dados observados pelo Interface Region Imaging Spectrograph (IRIS) sem o custo computacional proibitivo de inversões iterativas completas para cada pixel.

Metodologia
Os autores apresentam o IRIS2+, uma ferramenta de inversão aprimorada baseada em um algoritmo de vizinho mais próximo (k-nearest neighbor ou k-nn). O núcleo deste método é o banco de dados IRIS2+, uma extensão do banco de dados IRIS2 anterior, que agora contém 135.472 Perfis Representativos sintéticos (RPs).

• Construção do Banco de Dados: O banco de dados foi construído selecionando 550 conjuntos de dados do IRIS cobrindo várias características solares (umbra, poros, penumbra, sol quieto, plage). Estes foram pré-classificados em cinco classes. Para cada conjunto de dados, 320 RPs multilinhas foram calculados usando agrupamento k-means não supervisionado. Cada RP consiste em 6 linhas cromosféricas (C II 1334 & 1335 Å, Mg II h&k, triplete Mg II UV) e 6 linhas fotosféricas (Ti II, Fe I, C I, Ni II). Estes RPs foram invertidos usando o código de última geração STiC para gerar Atmosferas de Modelos Representativos (RMAs) contendo temperatura ($T$), velocidade na linha de visada ($v_{los}$), velocidade turbulenta ($v_{turb}$) e densidade eletrônica ($n_e$) através de 39 profundidades ópticas.
• Cálculo da Perda Radiativa: O banco de dados inclui valores de Perda Radiativa Integrada Representativa (RIRL). Estes foram calculados usando o STiC (utilizando o código RH) integrando a divergência do fluxo radiativo sobre a cromosfera (definida como $4 \le T \le 18$ kK). Para gerenciar a complexidade computacional durante a criação do banco de dados, as perdas radiativas foram computadas tratando as espécies (H, Mg II, Ca II) de forma independente (uma em não-LTE, as outras em LTE) em vez de simultaneamente, um método validado por superestimar as perdas totais em apenas $\approx 15\%$ em comparação com cálculos de não-LTE completos.
• Processo de Inversão: Para uma nova observação, o IRIS2+ transforma a amostragem espectral observada e os pesos para corresponder ao banco de dados. Ele então emprega um algoritmo k-nn (com $k=1$) para encontrar o RP sintético mais próximo no banco de dados em relação ao perfil observado. A RMA e a RIRL correspondentes são atribuídas à localização observada. A ferramenta oferece opções flexíveis de ponderação ("média", "padrão", "matriz", "manual") para lidar com diferentes combinações de linhas e configurações observacionais.

Principais Contribuições

1. Banco de Dados Expandido: O banco de dados IRIS2+ é mais de três vezes maior que o seu predecessor (135.472 vs. 40.320 RPs), cobrindo uma gama mais ampla de condições solares, da fotosfera à região de transição.
2. Perdas Radiativas Integradas: Ao contrário de ferramentas anteriores, o IRIS2+ fornece diretamente as perdas radiativas integradas para a cromosfera juntamente com os parâmetros termodinâmicos, derivados do banco de dados RMA.
3. Eficiência Computacional: A ferramenta é implementada em Python e utiliza bibliotecas padrão de aprendizado de máquina (Scikit-learn). Ela reduz o tempo de execução em ordens de magnitude ($O(10^5-10^6)$) em comparação ao STiC, permitindo que varreduras de raster típicas sejam invertidas em minutos.
4. Ferramenta de Inspeção Interativa: Uma ferramenta complementar permite que os usuários visualizem perfis observados vs. sintéticos, variáveis termodinâmicas ao longo da profundidade óptica, mapas de IRL e qualidade de ajuste ($\chi^2$).

Resultados
Os autores compararam os resultados do IRIS2+ com inversões do STiC para 10 conjuntos de dados do IRIS:

• Perdas Radiativas (IRL): Os resultados do IRIS2+ e do STiC mostram uma forte correlação linear ($r_p \approx 0.9$) após a filtragem de outliers. Os mapas espaciais são visualmente equivalentes, com discrepâncias atribuídas em grande parte a outliers ou eventos episódicos que exigem esquemas de inversão mais sofisticados. A ferramenta tem melhor desempenho quando as linhas de C II e Mg II são invertidas simultaneamente.
• Termodinâmica:
  • Temperatura ($T$): As correlações são de moderadas a fortes ($r_p \gtrsim 0.5$) entre $\log(\tau) = -7$ e $-1$. O IRIS2+ tende a produzir valores ligeiramente mais frios em altas alturas cromosféricas ($\log(\tau) = -6$), mas valores ligeiramente mais quentes na cromosfera média e na fotosfera ($\log(\tau) = -4, -1$), provavelmente devido a distribuições bimodais no ajuste do banco de dados.
  • Velocidades ($v_{los}, v_{turb}$): Existem correlações fortes na cromosfera alta e média ($r_p \approx 0.7$). No entanto, na cromosfera baixa e na fotosfera, as correlações enfraquecem, e o IRIS2+ geralmente produz valores menores que o STiC. Os autores atribuem isso ao limite artificial de $15$ km s$^{-1}$ imposto ao $v_{turb}$ no banco de dados e à sensibilidade limitada das linhas selecionadas às mudanças termodinâmicas na fotosfera.
• Sensibilidade de Ponderação: O estudo confirma que o uso de todas as linhas disponíveis (cromosféricas e fotosféricas) produz o melhor acordo com o STiC. Depender apenas de linhas fotosféricas ou apenas de linhas de C II resulta em correlações pobres para parâmetros cromosféricos.

Significância e Alegações
O artigo afirma que o IRIS2+ é uma ferramenta de inversão rápida e confiável, capaz de fornecer valores aproximados do estado termodinâmico e das perdas radiativas para uma ampla variedade de cenas solares observadas pelo IRIS. Os autores enfatizam que, embora o IRIS2+ não seja um substituto para o rigoroso código iterativo ST easily em todos os cenários, ele fornece resultados que são "comparáveis na maioria dos casos" ao STiC quando todas as linhas são invertidas simultaneamente.

A significância reside na capacidade de processar grandes volumes de dados do IRIS de forma eficiente, permitindo estudos estatísticos de aquecimento cromosférico e perdas radiativas que eram anteriormente proibitivos computacionalmente. Os autores observam modestamente que a precisão da ferramenta é limitada pela sensibilidade das linhas observadas e pela diversidade do banco de dados; eles sugerem que expandir o banco de dados para 300.000–500.000 elementos poderia melhorar ainda mais a cobertura e reduzir a degeneração. A ferramenta é apresentada como uma solução prática para a comunidade de física solar acessar dados termodinâmicos e de perda radiativa sem a pesada carga computacional das inversões de não-LTE completas.