A HINSA view of cosmic-ray ionization in IC 348 and NGC 1333: evidence for a strong low-energy cosmic-ray disparity

Sintesi Tecnica: Una visione HINSA dell'ionizzazione da raggi cosmici in IC 348 e NGC 1333

Problematica
Il tasso di ionizzazione da raggi cosmici (CRIR, $\zeta_{ion}$) è un parametro fondamentale che governa l'evoluzione chimica e dinamica delle nubi molecolari, in particolare regolando l'accoppiamento tra campi magnetici e gas. Sebbene osservazioni recenti abbiano rivelato valori di CRIR elevati nelle regioni di formazione stellare massiccia e vicino ai protostelle, le origini fisiche e i meccanismi di accelerazione dei raggi cosmici a bassa energia (LECR, $E < 1$ GeV) rimangono incerti. Inoltre, esiste una discrepanza nella vicina complessa di nubi molecolari di Perseo: mentre NGC 1333 mostra un deficit nei residui di raggi $\gamma$ sub-GeV diffusi, la sua vicina IC 348 presenta un residuo positivo e uno spettro energetico dei raggi cosmici più duro, nonostante NGC 1333 presenti un'attività di formazione stellare più elevata. L'entità con cui le sorgenti locali di raggi cosmici sovrastano il fondo galattico in queste regioni, e la natura della disparità tra esse, richiede chiarimenti.

Metodologia
Questo studio presenta nuove stime del CRIR verso IC 348 e NGC 1333 utilizzando la tecnica HINSA (H i narrow self-absorption). La metodologia prevede:

1. Osservazioni: Utilizzo di dati H i provenienti dal telescopio FAST (NGC 1333) e dal survey GALFA-H i (IC 348), combinati con dati d'archivio $^{13}$CO ($J=1-0$ e $3-2$) della survey COMPLETE e JCMT, e dati del continuo di polvere della Herschel Gould Belt Survey.
2. Analisi HINSA: Derivazione della colonna densità di idrogeno atomico $N(\text{HINSA})$ tramite l'adattamento di derivate del secondo ordine delle caratteristiche di assorbimento H i rispetto al campo di velocità $^{13}$CO. La frazione HINSA ($f_{\text{HINSA}} = N(\text{HINSA})/N(\text{H}_2)$) è calcolata per regioni con $N(\text{H}_2) > 3 \times 10^{21}$ cm$^{-2}$ per minimizzare gli effetti di fotodissociazione.
3. Vincoli di Densità: Stima della densità volumetrica del gas $n(\text{H}_2)$ mediante modellazione di trasferimento radiativo non-LTE (RADEX) vincolata dalle intensità delle linee $^{13}$CO e campionamento MCMC (Markov Chain Monte Carlo).
4. Calcolo del CRIR: Derivazione di $\zeta_{ion}$ dal bilancio in stato stazionario tra la formazione di H$_2$ sui grani di polvere e la distruzione da parte dei raggi cosmici (inclusi processi di ionizzazione, dissociazione e UV indotti dai raggi cosmici).
5. Modellazione Spettrale: Modellazione dello spettro LECR all'interno di un framework di attenuazione a lastra finita. Il modello incorpora vincoli dagli spettri di raggi cosmici ad alta energia inferiti dalle osservazioni di raggi $\gamma$ di Fermi per adattarsi ai profili di $\zeta_{ion}$ osservati, variando la pendenza spettrale a bassa energia ($\alpha$) e il break di momento ($p_0$).

Risultati Chiave

• Disparità del CRIR: Lo studio trova una pronunciata differenza di CRIR tra le due nubi. In IC 348, il CRIR medio ponderato varia da $(5,35 \pm 1,49) \times 10^{-17}$ s$^{-1}$ a $(2,33 \pm 0,55) \times 10^{-16}$ s$^{-1}$. Al contrario, NGC 1333 presenta valori circa un ordine di grandezza inferiori, che vanno da $(3,94 \pm 2,46) \times 10^{-18}$ s$^{-1}$ a $(1,68 \pm 0,25) \times 10^{-17}$ s$^{-1}$.
• Trend della Colonna Densità: In entrambe le nubi, $\zeta_{ion}$ diminuisce con l'aumentare della colonna densità di H$_2$, in linea con i modelli di attenuazione privi di sorgenti di accelerazione interne.
• Modellazione Spettrale: I migliori spettri di raggi cosmici a bassa energia riproducono i profili osservati. Il modello indica che il fattore di $\sim$10 di differenza nei tassi di ionizzazione è guidato da una corrispondente differenza nel flusso di protoni a bassa energia ($\sim$10 MeV). Lo spettro di IC 348 è significativamente più duro e intenso rispetto a quello di NGC 1333, che mostra un deficit rispetto al fondo galattico alle energie GeV.
• Analisi di Densità e Struttura: Gli autori investigano potenziali bias derivanti dalla stima della densità del gas e dalla struttura della nube (tramite simulazioni 3D-PDR). Sebbene le incertezze in $n(\text{H}_2)$ e le strutture frattali delle nubi possano influenzare i valori assoluti di $\zeta_{ion}$ (potenzialmente sottostimando a bassi $N(\text{H}_2)$ e sovrastimando ad alti $N(\text{H}_2)$), l'analisi conclude che questi fattori non possono spiegare la disparità di un ordine di grandezza tra le due regioni.

Significatività e Rivendicazioni
L'articolo sostiene che la disparità osservata nel CRIR sia intrinseca alle proprietà fisiche delle nubi piuttosto che un artefatto dei tempi evolutivi o dei bias di misurazione.

• Sorgenti di Accelerazione Esterne: L'alto CRIR in IC 348, unito al suo spettro $\gamma$ più duro, suggerisce la presenza di sorgenti locali di accelerazione LECR esterne ai nuclei densi delle nubi, probabilmente distinte dall'attività protostellare che domina NGC 1333. Gli autori propongono che gli shock di terminazione del vento stellare nel cluster di IC 348, più evoluto e meno immerso (contenente più stelle di tipo B), siano una fonte plausibile.
• Contributo Protostellare: I risultati suggeriscono che, sebbene i getti protostellari possano aumentare il CRIR su piccola scala (come visto nella mappatura di NGC 1333), essi non guidano il CRIR su larga scala in questi ambienti, specialmente dato che NGC 1333 ha un'attività di formazione stellare più elevata pur avendo un CRIR globale inferiore.
• Implicazioni per la Formazione Stellare: Il CRIR globalmente più alto in IC 348 implica una frazione di ionizzazione più elevata, portando a un più forte accoppiamento tra gas e campo magnetico. Gli autori notano che questo può potenziare il braking magnetico durante il collasso del nucleo, spiegando potenzialmente le scoperte osservative secondo cui i dischi protostellari in IC 348 sono mediamente più piccoli rispetto a quelli in NGC 1333.
• Validazione Metodologica: Lo studio valida la tecnica HINSA come sonda per il CRIR nelle regioni di formazione stellare a bassa massa, sottolineando la necessità di future modellazioni 3D per raffinare i valori assoluti pur confermando la robustezza delle disparità relative.