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

Sintesi Tecnica: Lo Strumento di Inversione IRIS2+

Definizione del Problema
Inferire lo stato fisico del plasma cromosferico dalle osservazioni rimane una sfida significativa nella fisica solare. A causa dei bassi tassi di collisione nella cromosfera, è necessaria una descrizione di equilibrio termodinamico non locale (non-LTE) per spiegare le intensità emergenti, impedendo una traduzione diretta dalle intensità osservate ai parametri fisici. Sebbene i moderni codici di inversione non-LTE (ad es., STiC) abbiano migliorato la fedeltà incorporando molteplici linee spettrali per aumentare la risoluzione in profondità, essi sono computazionalmente intensivi. Inoltre, il calcolo delle perdite radiative integrate (IRL) — cruciali per comprendere il bilancio energetico e i meccanismi di riscaldamento nella cromosfera — richiede una complessa post-elaborazione dei modelli atmosferici invertiti. Vi è la necessità di uno strumento veloce e affidabile capace di recuperare i parametri termodinamici e le perdite radiative per grandi set di dati osservati dall'Interface Region Imaging Spectrograph (IRIS) senza il costo computazionale proibitivo delle inversioni iterative complete per ogni pixel.

Metodologia
Gli autori introducono IRIS2+, uno strumento di inversione migliorato basato su un algoritmo dei k-vicini più prossimi (k-nearest neighbor o k-nn). Il cuore di questo metodo è il database IRIS2+, un'estensione del precedente database IRIS2, che ora contiene 135.472 Profili Rappresentativi (RP) sintetici.

• Costruzione del Database: Il database è stato costruito selezionando 550 dataset IRIS che coprono varie caratteristiche solari (umbra, pore-like, penumbra, quiet sun, plage). Questi sono stati pre-ordinati in cinque classi. Per ogni dataset, sono stati calcolati 320 RP multi-linea utilizzando il clustering k-means non supervisionato. Ogni RP consiste di 6 linee cromosferiche (C II 1334 & 1335 Å, Mg II h&k, Mg II UV triplet) e 6 linee fotosferiche (Ti II, Fe I, C I, Ni II). Questi RP sono stati invertiti utilizzando il codice all'avanguardia STiC per generare Atmosfere Modello Rappresentative (RMA) contenenti temperatura ($T$), velocità lungo la linea di vista ($v_{los}$), velocità turbolenta ($v_{turb}$) e densità elettronica ($n_e$) attraverso 39 profondità ottiche.
• Calcolo della Perdita Radiativa: Il database include valori di Perdita Radiativa Integrata Rappresentativa (RIRL). Questi sono stati calcolati utilizzando STiC (impiegando il codice RH) integrando la divergenza del flusso radiativo sulla cromosfera (definita come $4 \le T \le 18$ kK). Per gestire la complessità computazionale durante la creazione del database, le perdite radiative sono state calcolate trattando le specie (H, Mg II, Ca II) indipendentemente (una in non-LTE, le altre in LTE) anziché simultaneamente, un metodo validato per sovrastimare le perdite totali di solo $\approx 15\%$ rispetto ai calcoli non-LTE completi.
• Processo di Inversione: Per una nuova osservazione, IRIS2+ trasforma il campionamento spettrale osservato e i pesi per farli corrispondere al database. Utilizza quindi un algoritmo k-nn (con $k=1$) per trovare il RP sintetico più vicino nel database rispetto al profilo osservato. La RMA e la RIRL corrispondenti vengono assegnate alla posizione osservata. Lo strumento offre opzioni di pesatura flessibili ("mean", "default", "array", "manual") per gestire diverse combinazioni di linee e configurazioni osservative.

Contributi Chiave

1. Database Espanso: Il database IRIS2+ è più di tre volte più grande del suo predecessore (135.472 contro 40.320 RP), coprendo una gamma più ampia di condizioni solari dalla fotosfera alla regione di transizione.
2. Perdite Radiative Integrate: A differenza degli strumenti precedenti, IRIS2+ fornisce direttamente le perdite radiative integrate per la cromosfera insieme ai parametri termodinamici, derivati dal database RMA.
3. Efficienza Computazionale: Lo strumento è implementato in Python e utilizza librerie standard di machine learning (Scikit-learn). Riduce il tempo di esecuzione di ordini di grandezza ($O(10^5-10^6)$) rispetto a STiC, permettendo che i tipici raster scan siano invertiti in pochi minuti.
4. Strumento di Ispezione Interattivo: Uno strumento complementare permette agli utenti di visualizzare i profili osservati vs sintetici, le variabili termodinamiche lungo la profondità ottica, le mappe IRL e la qualità del fit ($\chi^2$).

Risultati
Gli autori hanno confrontato i risultati di IRIS2+ con le inversioni di STiC per 10 dataset IRIS:

• Perdite Radiative (IRL): I risultati di IRIS2+ e STiC mostrano una forte correlazione lineare ($r_p \approx 0.9$) dopo il filtraggio degli outlier. Le mappe spaziali sono visivamente equivalenti, con discrepanze attribuite in gran parte a outlier o eventi episodici che richiedono schemi di inversione più sofisticati. Lo strumento performa meglio quando le linee C II e Mg II sono invertite simultaneamente.
• Termodinamica:
  • Temperatura ($T$): Le correlazioni sono da moderate a forti ($r_p \gtrsim 0.5$) tra $\log(\tau) = -7$ e $-1$. IRIS2+ tende a produrre valori leggermente più freddi ad alte quote cromosferiche ($\log(\tau) = -6$) ma valori leggermente più caldi nella media cromosfera e nella fotosfera ($\log(\tau) = -4, -1$), probabilmente a causa di distribuzioni bimodali nel fitting del database.
  • Velocità ($v_{los}, v_{turb}$): Esistono forti correlazioni nella parte alta e media della cromosfera ($r_p \approx 0.7$). Tuttavia, nella bassa cromosfera e nella fotosfera, le correlazioni si indeboliscono e IRIS2+ generalmente produce valori inferiori rispetto a STC. Gli autori attribuiscono ciò al limite artificiale di $15$ km s$^{-1}$ imposto su $v_{turb}$ nel database e alla limitata sensibilità delle linee selezionate ai cambiamenti termodinamici nella fotosfera.
• Sensibilità della Pesatura: Lo studio conferma che l'uso di tutte le linee disponibili (cromosferiche e fotosferiche) produce l'accordo migliore con STiC. Affidarsi esclusivamente alle linee fotosferiche o solo alle linee C II risulta in scarse correlazioni per i parametri cromosferici.

Significato e Rivendicazioni
L'articolo afferma che IRIS2+ è uno strumento di inversione veloce e affidabile capace di fornire valori approssimativi dello stato termodinamico e delle perdite radiative per una vasta gamma di scene solari osservate da IRIS. Gli autori sottolineano che, sebbene IRIS2+ non sia un sostituto del rigoroso codice iterativo STiC in tutti gli scenari, fornisce risultati che sono "comparabili nella maggior parte dei casi" a STiC quando tutte le linee sono invertite simultaneamente.

Il significato risiede nella capacità di elaborare grandi volumi di dati IRIS in modo efficiente, consentendo studi statistici sul riscaldamento cromosferico e sulle perdite radiative che prima erano computazionalmente proibitivi. Gli autori notano con modestia che l'accuratezza dello strumento è limitata dalla sensibilità delle linee osservate e dalla diversità del database; suggeriscono che espandere il database a 300.000–500.000 elementi potrebbe migliorare ulteriormente la copertura e ridurre la degenerazione. Lo strumento è presentato come una soluzione pratica per la comunità della fisica solare per accedere ai dati termodinamici e di perdita radiativa senza l'elevato carico computazionale delle inversioni non-LTE complete.