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

Resumen Técnico: La Herramienta de Inversión IRIS2+

Planteamiento del Problema
Inferir el estado físico del plasma cromosférico a partir de observaciones sigue siendo un desafío significativo en la física solar. Debido a las bajas tasas de colisión en la cromosfera, se requiere una descripción de equilibrio termodinámico no local (no-LTE) para explicar las intensidades emergentes, lo que impide una traducción directa de las intensidades observadas a los parámetros físicos. Si bien los códigos de inversión no-LTE modernos (por ejemplo, STiC) han mejorado la fidelidad al incorporar múltiples líneas espectrales para aumentar la resolución de profundidad, son computacionalmente intensivos. Además, calcular las pérdidas radiativas integradas (IRL)—cruciales para comprender el balance energético y los mecanismos de calentamiento en la cromosfera—requiere un complejo post-procesamiento de los modelos atmosféricos invertidos. Existe la necesidad de una herramienta rápida y fiable capaz de recuperar los parámetros termodinámicos y las pérdidas radiativas para grandes conjuntos de datos observados por el Interface Region Imaging Spectrograph (IRIS) sin el costo computacional prohibitivo de las inversiones iterativas completas para cada píxel.

Metodología
Los autores presentan IRIS2+, una herramienta de inversión mejorada basada en un algoritmo de k-vecinos más cercanos (k-nearest neighbor o k-nn). El núcleo de este método es la base de datos IRIS2+, una extensión de la base de datos IRIS2 anterior, que ahora contiene 135,472 Perfiles Representativos (RPs) sintéticos.

• Construcción de la Base de Datos: La base de datos se construyó seleccionando 550 conjuntos de datos de IRIS que cubren diversas características solares (umbra, tipo poro, penumbra, sol tranquilo, plage). Estos fueron pre-clasificados en cinco clases. Para cada conjunto de datos, se calcularon 320 RPs multilínea utilizando el agrupamiento k-means no supervisado. Cada RP consta de 6 líneas cromosféricas (C II 1334 & 1335 Å, Mg II h&k, triplete Mg II UV) y 6 líneas fotosféricas (Ti II, Fe I, C I, Ni II). Estos RPs fueron invertidos utilizando el código de vanguardia STiC para generar Atmósferas de Modelo Representativas (RMAs) que contienen temperatura ($T$), velocidad en la línea de visión ($v_{los}$), velocidad turbulenta ($v_{turb}$) y densidad electrónica ($n_e$) a través de 39 profundidades ópticas.
• Cálculo de la Pérdida Radiativa: La base de datos incluye valores de Pérdida Radiativa Integrada Representativa (RIRL). Estos se calcularon utilizando STiC (empleando el código RH) integrando la divergencia del flujo radiativo sobre la cromosfera (definida como $4 \le T \le 18$ kK). Para gestionar la complejidad computacional durante la creación de la base de datos, las pérdidas radiativas se computaron tratando las especies (H, Mg II, Ca II) de forma independiente (una en no-LTE, las otras en LTE) en lugar de simultáneamente, un método validado que sobreestima las pérdidas totales en solo $\approx 15\%$ en comparación con los cálculos de no-LTE completo.
• Proceso de Inversión: Para una nueva observación, IRIS2+ transforma el muestreo espectral observado y los pesos para que coincidan con la base de datos. Luego emplea un algoritmo k-nn (con $k=1$) para encontrar el RP sintético más cercano en la base de datos al perfil observado. La RMA y la RIRL correspondientes se asignan a la ubicación observada. La herramienta ofrece opciones de ponderación flexibles ("media", "por defecto", "matriz", "manual") para manejar diferentes combinaciones de líneas y configuraciones de observación.

Contribuciones Clave

1. Base de Datos Expandida: La base de datos IRIS2+ es más de tres veces más grande que su predecesora (135,472 frente a 40,320 RPs), cubriendo un rango más amplio de condiciones solares desde la fotosfera hasta la región de transición.
2. Pérdidas Radiativas Integradas: A diferencia de herramientas anteriores, IRIS2+ proporciona directamente las pérdidas radiativas integradas para la cromosfera junto con los parámetros termodinámicos, derivados de la base de datos RMA.
3. Eficiencia Computacional: La herramienta está implementada en Python y utiliza librerías estándar de aprendizaje automático (Scikit-learn). Reduce el tiempo de ejecución en órdenes de magnitud ($O(10^5-10^6)$) en comparación con STiC, permitiendo que los escaneos de rastro típicos se inviertan en minutos.
4. Herramienta de Inspección Interactiva: Una herramienta complementaria permite a los usuarios visualizar perfiles observados vs. sintéticos, variables termodinámicas a lo largo de la profundidad óptica, mapas de IRL y la calidad del ajuste ($\chi^2$).

Resultados
Los autores compararon los resultados de IRIS2+ contra las inversiones de STiC para 10 conjuntos de datos de IRIS:

• Pérdidas Radiativas (IRL): Los resultados de IRIS2+ y STiC muestran una fuerte correlación lineal ($r_p \approx 0.9$) tras filtrar los valores atípicos. Los mapas espaciales son visualmente equivalentes, con discrepancias atribuidas en gran medida a valores atípicos o eventos episódicos que requieren esquemas de inversión más sofisticados. La herramienta funciona mejor cuando las líneas de C II y Mg II se invierten simultáneamente.
• Termodinámica:
  • Temperatura ($T$): Las correlaciones son de moderadas a fuertes ($r_p \gtrsim 0.5$) entre $\log(\tau) = -7$ y $-1$. IRIS2+ tiende a producir valores ligeramente más fríos en alturas cromosféricas elevadas ($\log(\tau) = -6$) pero valores ligeramente más calientes en la cromosfera media y la fotosfera ($\log(\tau) = -4, -1$), probablemente debido a distribuciones bimodales en el ajuste de la base de datos.
  • Velocidades ($v_{los}, v_{turb}$): Existen correlaciones fuertes en la cromosfera alta y media ($r_p \approx 0.7$). Sin embargo, en la cromosfera baja y la fotosfera, las correlaciones se debilitan y el IRIS2+ generalmente arroja valores más bajos que STiC. Los autores atribuyen esto al límite artificial de $15$ km s$^{-1}$ impuesto a $v_{turb}$ en la base de datos y a la sensibilidad limitada de las líneas seleccionadas a los cambios termodinámicos en la fotosfera.
• Sensibilidad de Ponderación: El estudio confirma que el uso de todas las líneas disponibles (cromosféricas y fotosféricas) produce la mejor concordancia con STiC. Depender únicamente de las líneas fotosféricas o solo de las líneas de C II resulta en pobres correlaciones para los parámetros cromosféricos.

Significancia y Reivindicaciones
El artículo afirma que IRIS2+ es una herramienta de inversión rápida y fiable capaz de proporcionar valores aproximados del estado termodinámico y las pérdidas radiativas para una amplia variedad de escenas solares observadas por IRIS. Los autores enfatizan que, si bien IRIS2+ no es un reemplazo para el riguroso código iterativo STiC en todos los escenarios, proporciona resultados que son "comparables en la mayoría de los casos" a STiC cuando todas las líneas se invierten simultáneamente.

La significancia radica en la capacidad de procesar grandes volúmenes de datos de IRIS de manera eficiente, permitiendo estudios estadísticos del calentamiento cromosférico y las pérdidas radiativas que anteriormente eran computacionalmente prohibitivos. Los autores señalan modestamente que la precisión de la herramienta está limitada por la sensibilidad de las líneas observadas y la diversidad de la base de datos; sugieren que expandir la base de datos a 300,000–500,000 elementos podría mejorar aún más la cobertura y reducir la degeneración. La herramienta se presenta como una solución práctica para que la comunidad de física solar acceda a datos termodinámicos y de pérdida radiativa sin la pesada carga computacional de las inversiones de no-LTE completas.