Stochastic modelling of cosmic-ray sources for Galactic diffuse emissions

Questo studio utilizza una simulazione Monte Carlo stocastica per dimostrare che la natura discreta dei resti di supernova come sorgenti di raggi cosmici introduce incertezze significative nelle emissioni diffuse galattiche, specialmente a energie superiori a decine di TeV nel caso di diffusione dipendente dal tempo, offrendo così un potenziale strumento per conciliare i modelli teorici con le osservazioni di LHAASO e per localizzare le sorgenti cosmiche.

L'essenziale

🌌 La Galassia come un Gigantesco Concerto di Fuochi d'Artificio

Immagina la nostra Galassia, la Via Lattea, non come una nebbia uniforme e tranquilla, ma come un enorme concerto di fuochi d'artificio che esplode nel cielo notturno.

Ogni "fuoco d'artificio" è una Supernova (l'esplosione di una stella morente). Queste esplosioni sono le "fabbriche" che lanciano nello spazio particelle incredibilmente veloci chiamate Raggi Cosmici. Quando queste particelle viaggiano e colpiscono il gas che riempie lo spazio tra le stelle, producono una luce speciale: i raggi gamma (e anche i neutrini). Questa luce diffusa è ciò che gli astronomi chiamano "emissione diffusa galattica".

Finora, gli scienziati hanno spesso trattato questi fuochi d'artificio come se fossero una nebbia luminosa continua, come se la luce fosse distribuita in modo perfetto e uniforme in tutto il cielo. È un modello comodo, come se la galassia fosse un unico grande lampadario acceso.

Ma la realtà è diversa.
In questo studio, gli autori (Anton e Philipp) ci dicono: "Aspetta un attimo! I fuochi d'artificio non sono una nebbia. Sono singoli, distinti e sparsi in modo casuale. E questo cambia tutto."

🔍 Il Problema: La "Macchia" Casuale

Immagina di guardare un muro bianco illuminato da una sola lampadina. Se ci sono molte lampadine vicine, il muro sembra illuminato in modo uniforme. Ma se ci sono solo poche lampadine e sono lontane tra loro, vedrai delle zone molto luminose (dove c'è la lampadina) e delle zone buie (dove non c'è).

Gli scienziati hanno usato un potente simulatore al computer (un "mondo virtuale") per lanciare migliaia di "lampadine" (supernove) in posizioni casuali nella galassia e vedere cosa succede alla luce che ne risulta. Hanno scoperto tre cose fondamentali:

1. La Galassia non è "liscia", è "a macchie"

Quando guardi la galassia con occhi molto sensibili (come i telescopi moderni LHAASO o IceCube), non vedi una luce uniforme. Vedi delle macchie e delle strisce più luminose.

• L'analogia: È come guardare la pioggia. Se ne vedi una pioggia fitta e continua, sembra una cortina d'acqua. Ma se guardi una goccia alla volta, vedi che sono oggetti distinti che cadono in punti specifici.
• La scoperta: In certi modelli, queste macchie possono essere così intense da far sembrare che la luce in quel punto sia doppia o tripla rispetto a quanto ci si aspetterebbe se la galassia fosse "liscia".

2. Tre modi diversi di "lanciare" le particelle

Gli scienziati hanno testato tre scenari su come queste "lampadine" lanciano le particelle:

• Scenario A (Il "Burst"): Tutte le particelle escono tutte insieme, subito dopo l'esplosione. È come un singolo scatto di flash. Qui, le macchie di luce sono visibili, ma non sono troppo strane.
• Scenario B (La "Fuga Energetica"): Le particelle più veloci scappano prima, quelle più lente dopo. È come se il fuoco d'artificio avesse colori che escono in momenti diversi. Questo rende le macchie di luce ancora più difficili da prevedere e meno intense alle energie più alte.
• Scenario C (Il "Labirinto"): Questo è il più interessante! Immagina che intorno a ogni esplosione ci sia una zona di nebbia fitta che intrappola le particelle per un po' di tempo prima di lasciarle andare. In questo caso, le particelle rimangono ammassate vicino alla fonte.
  • Il risultato: Le macchie di luce diventano enormi e mostruose. In alcuni punti del cielo, la luce può essere migliaia di volte più intensa di quanto previsto dai modelli "lisci". È come se il fuoco d'artificio non si fosse spento, ma avesse acceso un riflettore potente proprio lì.

3. Perché questo è importante?

Perché gli astronomi stanno guardando il cielo con telescopi sempre più potenti (come LHAASO in Cina o IceCube in Antartide).

• Se il modello "liscio" è sbagliato, le previsioni degli scienziati non combaciano con ciò che vedono.
• Questo studio dice: "Non preoccupatevi se i vostri dati non corrispondono perfettamente alle previsioni vecchie. Potrebbe essere solo perché una supernova vicina e giovane sta creando una 'macchia' di luce che il modello liscio non vedeva!"

🧩 Il Messaggio Finale: Cosa dobbiamo imparare?

1. La casualità conta: Non possiamo più trattare la galassia come un fluido uniforme. Dobbiamo considerare che le fonti sono discrete e casuali, come i granelli di sabbia su una spiaggia.
2. L'incertezza è reale: C'è una nuova fonte di "errore" o "incertezza" nelle nostre previsioni. Non è un errore di calcolo, ma una caratteristica naturale dell'universo: a volte siamo fortunati e abbiamo una supernova vicina che illumina tutto; altre volte siamo nel buio.
3. Il futuro è nelle "macchie": Più i telescopi saranno capaci di vedere dettagli piccoli (alta risoluzione), più vedranno queste macchie. Questo ci permetterà di trovare le supernove che non abbiamo ancora visto e capire come funzionano queste "fabbriche" di particelle.

In sintesi:
Questo studio ci insegna che la nostra galassia è un luogo più "disordinato" e "spettacolare" di quanto pensassimo. Invece di una luce uniforme, abbiamo un mosaico di lampi brillanti causati da esplosioni stellari vicine. Capire questo "disordine" è la chiave per decifrare i messaggi che la galassia ci invia con i raggi cosmici.

Sommario tecnico

Ecco un riassunto tecnico dettagliato del paper "Stochastic modelling of cosmic-ray sources for Galactic diffuse emissions", presentato in italiano.

Titolo: Modellazione stocastica delle sorgenti di raggi cosmici per le emissioni diffuse galattiche

Autori: Anton Stalla e Philipp Mertsch (RWTH Aachen University)
Contesto: Studio delle emissioni diffuse galattiche (GDE) in raggi gamma e neutrini, con focus sull'impatto della natura discreta delle sorgenti di raggi cosmici (CR).

---

1. Il Problema

Le emissioni diffuse galattiche (GDE) in raggi gamma e neutrini derivano dall'interazione dei raggi cosmici (CR) con il mezzo interstellare (ISM). Queste emissioni sono fondamentali per sondare l'intensità dei raggi cosmici al di fuori del Sistema Solare. Tuttavia, le previsioni dei modelli teorici sono influenzate dalle proprietà delle sorgenti dei raggi cosmici.

• Limitazione attuale: La maggior parte dei modelli tratta le sorgenti (principalmente i Resti di Supernova, SNR) come una distribuzione continua e liscia.
• La realtà fisica: Le sorgenti sono entità discrete, con posizioni e tempi di iniezione specifici. Poiché i tempi di fuga dei CR aumentano con l'energia, alle alte energie (decine di TeV e oltre) solo poche sorgenti vicine contribuiscono significativamente al flusso osservato.
• Domanda di ricerca: Quanto la natura discreta e stocastica delle sorgenti influisce sulle previsioni delle GDE? È possibile che questa "stocasticità delle sorgenti" spieghi discrepanze tra modelli e dati osservativi (es. LHAASO, IceCube)?

2. Metodologia

Gli autori hanno sviluppato un modello Monte Carlo stocastico per calcolare le GDE hadroniche (prodotte da protoni) considerando la natura puntuale delle sorgenti.

• Modello di Trasporto:

  • Si assume un volume di trasporto galattico con un alone di altezza $H = 4$ kpc.
  • Le sorgenti sono distribuite secondo un modello a spirale con barra galattica, con età estratte da una distribuzione uniforme.
  • L'equazione di trasporto dei CR è risolta analiticamente utilizzando la funzione di Green, assumendo condizioni al contorno di fuga libera.
  • È stata implementata una condizione di "cono di luce" per evitare violazioni di causalità e un'estensione finita delle sorgenti (10 pc) per regolarizzare le divergenze.

• Scenari di Iniezione e Trasporto Vicino alla Sorgente:
  Sono stati studiati tre scenari principali per testare diverse ipotesi fisiche:

  1. Scenario "Burst-like" (Impulsivo): Tutte le rigidità fuggono dalla sorgente contemporaneamente; il coefficiente di diffusione è costante nel tempo.
  2. Fuga Dipendente dall'Energia (Energy-dependent escape): Le particelle ad alta rigidità fuggono prima di quelle a bassa rigidità (modello CREDIT), ma la diffusione rimane costante nel tempo.
  3. Diffusione Dipendente dal Tempo (Time-dependent Diffusion - TDD): Tutte le particelle fuggono insieme, ma subiscono inizialmente una diffusione soppressa (coefficiente costante e basso) per un periodo di 10 kyr a causa della turbolenza auto-generata, prima di passare alla diffusione standard dipendente dalla rigidità.

• Calcolo delle Emissioni:

  • Le intensità dei CR sono calcolate su una griglia spaziale (HEALPix) per 1.000 realizzazioni diverse della distribuzione delle sorgenti.
  • Le emissioni gamma sono ottenute integrando l'intensità dei CR lungo la linea di vista, ponderata per la densità del gas (HI e H2) e le sezioni d'urto di interazione.
  • È stato incluso l'effetto di assorbimento dei fotoni ad alte energie (>10 TeV) dovuto alla produzione di coppie $\gamma\gamma \to e^+e^-$.
  • Efficienza Computazionale: L'uso di GPU e librerie come jax ha permesso di accelerare i calcoli di un ordine di grandezza, rendendo fattibile l'analisi di migliaia di realizzazioni.

3. Contributi Chiave e Risultati

A. Distribuzioni Statistiche delle Intensità

• Le distribuzioni delle intensità delle GDE non seguono una semplice legge normale (Gaussiana).
• Legge Stabile: Le distribuzioni sono meglio descritte da leggi stabili (stable laws) con code a potenza (divergenza della varianza), specialmente lungo le linee di vista nel piano galattico.
  • Per intensità vicine alla media, la distribuzione è meglio descritta da una legge stabile con indice $\alpha = 4/3$ (distribuzione 2D delle sorgenti).
  • Per le code ad alta intensità (eventi estremi), $\alpha = 5/3$ (distribuzione 3D) è più appropriato.
• Modelli Misti: Per le linee di vista fuori dal piano galattico, la distribuzione è una miscela di una legge stabile (contributo di sorgenti vicine) e una distribuzione Gaussiana (contributo di sorgenti lontane).

B. Deviazioni Massime e Morfologia

• Scenario Burst-like e Fuga Dipendente dall'Energia: Le deviazioni massime rispetto al modello liscio sono tipicamente dell'ordine di decine di percentuali (es. 20-50% a 100 TeV). Le strutture sono visibili ma meno estreme.
• Scenario Diffusione Dipendente dal Tempo (TDD): Le deviazioni sono estreme, raggiungendo ordini di grandezza (fino a diverse migliaia di percentuali a 100 TeV). La soppressione della diffusione vicino alle sorgenti crea "macchie" di emissione molto intense e localizzate, sfocando il confine tra emissione diffusa e sorgenti puntuali.

C. Confronto con i Dati (LHAASO)

• Finestre Celesti (Sky Windows): Quando le emissioni sono mediate su grandi regioni del cielo (come le finestre di LHAASO), l'incertezza dovuta alla stocasticità delle sorgenti è subdominante negli scenari Burst-like e Fuga Dipendente dall'Energia (sotto il 10% a 100 TeV). In questi casi, la stocasticità non può spiegare il disaccordo tra i modelli lisci e i dati di LHAASO.
• Scenario TDD: In questo scenario, l'incertezza diventa significativa sopra alcune decine di TeV. Le bande di incertezza si allargano notevolmente, permettendo al modello di conciliare le previsioni con le misure di LHAASO senza bisogno di parametri aggiuntivi.

D. Correlazioni Energetiche

• Nei modelli lisci, le mappe a 10 GeV e 100 TeV sono fortemente correlate.
• Nei modelli discreti, questa correlazione si rompe per le deviazioni ($\Delta J$) rispetto alla media, poiché le sorgenti dominanti cambiano con l'energia.
• Tuttavia, nello scenario TDD, le deviazioni a diverse energie rimangono correlate perché il trasporto iniziale è indipendente dall'energia.
• La correlazione delle intensità totali ($J$) rimane alta negli scenari moderati, ma si rompe nello scenario TDD a causa dell'impatto estremo delle sorgenti vicine.

4. Significato e Conclusioni

• Nuova Fonte di Incertezza: La stocasticità delle sorgenti rappresenta una fonte di incertezza modellistica non trascurabile, specialmente alle alte energie (>10-100 TeV) e in scenari di trasporto complessi (come la diffusione dipendente dal tempo).
• Strumento Diagnostico: Le misure delle GDE con alta risoluzione spaziale (prossime generazioni di osservatori come SWGO o analisi avanzate di LHAASO) possono essere utilizzate non solo per mappare il gas, ma per vincolare i modelli di iniezione e trasporto delle sorgenti.
• Rilevanza per i PeVatroni: La mancata osservazione di "macchie" estreme di emissione diffusa potrebbe indicare che il tasso di sorgenti PeVatroni è inferiore a quanto ipotizzato, oppure che lo scenario di diffusione dipendente dal tempo non è corretto.
• Neutrini: Poiché le emissioni gamma e neutrini sono correlate (rapporto $\pi^0/\pi^\pm$), queste conclusioni si applicano anche alle previsioni per i neutrini galattici (IceCube).

In sintesi, il lavoro dimostra che ignorare la natura discreta delle sorgenti può portare a sottostimare le incertezze nelle previsioni delle emissioni diffuse ad alta energia, e che l'analisi delle fluttuazioni spaziali offre una via promettente per comprendere la fisica dell'accelerazione e del confinamento dei raggi cosmici nella nostra Galassia.