The Effects of Cosmic Ray Protons on Galactic Nonthermal Filaments

Sintesi Tecnica: Gli Effetti dei Protoni di Raggi Cosmici sui Filamenti Nontermici Galattici

Problematica
Il Centro Galattico (GC) ospita una popolazione di filamenti nontermici (NTF) che si estendono per decine di parsec di lunghezza, illuminati dalla radiazione di sincrotrone da parte dei raggi cosmici (CR). L'origine fisica di queste strutture rimane irrisolta. Sono proposte due ipotesi principali per il meccanismo di iniezione: (1) gli NTF sono alimentati da getti provenienti da nebulose di vento di pulsar, con una popolazione dominata da leptoni; oppure (2) gli NTF sono alimentati da particelle accelerate da shock interstellari, con una popolazione dominata da protoni (con un rapporto protone-leptone di circa 100:1). Sebbene la popolazione di CR galattici esibisca una distribuzione di energia a legge di potenza coerente con le osservazioni degli NTF, il collegamento causale specifico tra la sorgente di iniezione e la morfologia del filamento rimane poco chiaro. Inoltre, la topologia del campo magnetico poloidale all'interno della Central Molecular Zone (CMZ) è complessa, con molti filamenti che appaiono non lineari, suggerendo il potenziale ruolo della turbolenza interstellare. Questo studio mira a distinguere tra gli scenari dominati dai leptoni e quelli dominati dai protoni simulando la propagazione dei CR, il riscaldamento e il flusso di plasma sotto diverse condizioni del Centro Galattico.

Metodologia
Gli autori hanno impiegato il codice di magnetoidrodinamica (MHD) Athena++, modificato per includere termini di perdita radiativi e collisionali, per simulare la propagazione dei CR in un tubo di flusso 1D che rappresenta un filamento. La simulazione modella un mezzo di fondo uniforme perturbato dall'iniezione di energia dei CR nel punto medio del filamento.

• Processi Fisici: Il modello incorpora la diffusione (scattering su onde Alfvén) e lo streaming (convezione alla velocità di Alfvén). I meccanismi di perdita per i leptoni includono la radiazione di sincrotrone, lo scattering Compton inverso, le perdite Coulomb e il bremsstrahlung. La produzione secondaria è modellata dove i protoni primari interagiscono con i protoni ambientali per generare elettroni e positroni secondari.
• Parametri: Lo studio varia parametri chiave tra cui l'intensità del campo magnetico ($B$), la densità del plasma ($n_p$), il coefficiente di diffusione dei CR ($\kappa$) e l'energia totale iniettata ($E_{\pm}$).
• Set di Simulazioni: Ogni set di parametri è stato simulato due volte: una volta con soli leptoni e una volta con una miscela protone-leptone (scalando il tasso di iniezione totale di 100 per tenere conto del rapporto protone-leptone). I parametri fiduciali adottati sono $B = 200\,\mu\text{G}$, $n_p = 100\,\text{cm}^{-3}$ e $\kappa = 1 \times 10^{26}\,\text{cm}^2\,\text{s}^{-1}$.
• Analisi: Gli autori hanno analizzato i profili spaziali della densità di energia, i flussi di gas indotti, i tassi di riscaldamento e la risultante luminosità di sincrotrone e gli indici spettrali nel tempo (fino a 40 kyr).

Risultati Chiave

1. Indistinguibilità dei Meccanismi di Iniezione: Le simulazioni hanno rivelato scarse differenze osservabili tra i casi dominati dai protoni e quelli dominati dai leptoni. Sebbene l'iniezione di protoni guidi flussi di gas più forti e fornisca maggiore riscaldamento a causa della maggiore densità di energia totale, questi effetti rimangono al di sotto delle soglie rilevabili. Anche alle densità più basse considerate, le velocità indotte ($\sim 2.5\,\text{km s}^{-1}$) e i tassi di riscaldamento sono trascurabili rispetto al raffreddamento radiativo e alle tipiche larghezze di riga termica.
2. Vincoli di Propagazione: L'osservata luminosità superficiale costante e uniforme degli NTF suggerisce un coefficiente di diffusione significativamente inferiore al valore galattico tipico ($\kappa \sim 10^{28}\,\text{cm}^2\,\text{s}^{-1}$). Gli autori riscontrano che un coefficiente di diffusione più basso ($\kappa \sim 10^{26}\,\text{cm}^2\,\text{s}^{-1}$) è necessario per riprodurre i profili di luminosità osservati.
3. Lunghezza e Età del Filamento: Le simulazioni indicano che, senza un meccanismo per troncare la propagazione, i filamenti continuerebbero ad espandersi indefinitamente. La lunghezza di un NTF dipende quindi dalla sua età. Le perdite radiative e collisionali, anche se combinate con la produzione secondaria, non creano una lunghezza in stato stazionario perché i leptoni secondari nascono più lontano dal sito di iniezione e si diffondono in modo simile alle particelle primarie.
4. Variabilità dell'Indice Spettrale: Lo studio ha esplorato i meccanismi per la variazione osservata degli indici spettrali di sincrotrone (che variano da -2 a +0.5).
   • Diffusione dipendente dall'energia: Se la diffusività scala con l'energia ($\kappa \propto \gamma^q$), l'indice spettrale varia lungo il filamento, generalmente indurendosi allontanandosi dal sito di iniezione.
   • Invecchiamento da sincrotrone: Le particelle ad alta energia perdono energia più velocemente, portando a un irrigidimento spettrale nel tempo.
   • Produzione secondaria: Indurisce lo spettro nei tempi tardivi.
   • Nonostante questi meccanismi, gli autori notano che nessuna combinazione di questi effetti spiega completamente le variazioni di indice spettrale irregolari e asimmetriche osservate nei reali NTF.

Significato e Conclusioni
Il contributo principale di questo lavoro è la dimostazione del fatto che la presenza di una vasta popolazione di protoni cosmici non produce firme distinte e osservabili nella morfologia, nel flusso o nel riscaldamento degli NTF che permettano agli astronomi di distinguere tra origini guidate da pulsar (leptoni) o da shock (protoni) basandosi sulle attuali simulazioni MHD 1D.

Di conseguenza, gli autori sostengono che l'incapacità di distinguere tra questi due modelli di iniezione standard, unita alle imperfezioni strutturali osservate negli NTF (pieghe, divisioni e allineamenti non lineari), motiva la considerazione di un terzo meccanismo di formazione. Il documento propone che gli NTF possano derivare da strutture intermittenti della turbolenza del Centro Galattico. In questo scenario, la turbolenza MHD disomogenea potrebbe amplificare i campi magnetici e stirare gli eddy in sottili filamenti, creando le strutture osservate senza richiedere una specifica sorgente di iniezione di CR localizzata. Gli autori concludono che il lavoro futuro dovrebbe concentrarsi sulla simulazione di questi effetti turbolenti per comprendere meglio l'origine degli NTF.