High-Energy Neutrino Tomography of the Earth's Interior with IceCube

Sintesi Tecnica: Tomografia ad Alta Energia del Neutrino dell'Interno della Terra con IceCube

Problema e Motivazione
La struttura interna della Terra, specificamente la sua distribuzione radiale della densità, è fondamentale per comprendere l'evoluzione geologica, il geodinamismo e la formazione planetaria. Storicamente, queste proprietà sono state inferite principalmente attraverso la sismologia (analisi della propagazione delle onde sismiche) e la gravimetria, che sondano la risposta macroscopica della materia alle forze elastiche e gravitazionali. Sebbene il Modello di Riferimento Preliminare della Terra (PREM) fornisca uno standard robusto basato su questi metodi, esso si basa su inferenze indirette. Un approccio fisico distinto consiste nel sondare la densità della Terra utilizzando l'interazione debole. I neutrini ad alta energia, che interagiscono solo tramite la forza debole, attraversano la Terra ma subiscono un'attenuazione dipendente dalla densità di colonna della materia che incontrano. All'aumentare dell'energia del neutrino, la sezione d'urto aumenta, rendendo la Terra parzialmente opaca ai neutrini sopra i $\sim$10 TeV. Questa attenuazione è sensibile al profilo di densità lungo la traiettoria del neutrino, offrendo una sonda innovativa e indipendente dei profili planetari.

Metodologia
Questo studio utilizza 10,7 anni di dati (2011–2022) dell'Osservatorio Neutrinico IceCube al Polo Sud. L'analisi si concentra su 368.071 eventi di muon-neutrino ascendenti nell'intervallo di energia compreso tra 500 GeV e 100 TeV. Questi eventi includono sia neutrini atmosferici (prodotti dalle interazioni dei raggi cosmici nell'atmosfera) sia un flusso astrofisico diffuso subdominante.

La metodologia procede come segue:

1. Selezione degli Eventi: Il dataset comprende tracce di muoni "partenti" (starting) e "attraversanti" (throughgoing). Le tracce partenti (interazioni all'interno del volume strumentato) forniscono una ricostruzione dell'energia superiore, mentre le tracce attraversanti (interazioni all'esterno del rivelatore) offrono un'alta statistica.
2. Modellazione: Il profilo di densità della Terra è parametrizzato come un insieme di gusci concentrici a densità uniforme. Lo studio impiega principalmente una configurazione a cinque gusci (nucleo interno, nucleo esterno e strati mantello/crosta) e valida i risultati con un'estensione a otto gusci.
3. Simulazione e Verosimiglianza: Simulazioni Monte Carlo modellano l'accettanza del rivelatore IceCube, le incertezze sistematiche (incluse le proprietà del ghiaccio glaciale) e il flusso di neutrini atteso. L'analisi tiene conto della sezione d'urto neutrino-nucleone del Modello Standard, che è ben vincolata dalle funzioni di distribuzione dei partoni (PDF) e dai dati di scattering profondo incompleto (DIS).
4. Inferenza: Viene eseguita un'analisi di fit congiunto utilizzando una verosimiglianza di Poisson a classi (binned Poisson likelihood) per vincolare le densità dei gusci, la normalizzazione del flusso, le incertezze della sezione d'urto e i parametri di disturbo (nuisance parameters) del rivelatore. Le distribuzioni a posteriori vengono campionate utilizzando un algoritmo Markov Chain Monte Carlo (MCMC).
5. Analisi Cieca (Blind Analysis): Per prevenire bias, l'analisi è stata condotta in due fasi: un fit cieco escludendo gli eventi quasi orizzontali ($-1,0 < \cos\theta < -0,1$), dove le profondità di colonna sono minime e la mismodellazione sistematica è più probabile, seguito da un fit utilizzando tutti i dati disponibili ($-1,0 < \cos\theta < 0,0$).

Risultati Chiave
L'analisi produce i seguenti risultati primari:

• Profilo di Densità: Il profilo di densità radiale inferito, derivato dal modello a cinque gusci, è coerente con il modello PREM entro gli intervalli di credibilità dell'HPD (highest-posterior-density) del 68%. I risultati non mostrano deviazioni significative dal profilo PREM, validando la sonda tramite l'interazione debole.
  • Guscio 1 (Nucleo Interno): Moda tra 0,60–0,73 volte la densità PREM.
  • Gusci 2 e 3 (Nucleo Esterno): Le mode variano tra 0,72 e 1,1 volte la densità PREM.
  • Gusci 4 e 5 (Mantello/Crosta): Le mode variano tra 1,3 e 1,6 volte la densità PREM.
  • L'inclusione degli eventi quasi orizzontali sposta leggermente la moda a posteriori dello strato più esterno (Guscio 5) più vicino all'unità, suggerendo che questi eventi contribuiscono costruttivamente ai vincoli.
• Quantità Globali Derivate:
  • Massa della Terra: La moda a posteriori per la massa della Terra è $M = 7,25 \times 10^{24}$ kg, con un intervallo HPD al 68% di $[6,31, 8,31] \times 10^{24}$ kg. Sebbene la moda sia superiore al valore di riferimento gravitazionale ($5,97 \times 10^{24}$ kg), il valore di riferimento rientra nell'intervallo HPD al 95%. L'ipotesi di massa nulla è esclusa a $>5\sigma$.
  • Momento d'Inerzia Polare: La moda a posteriori è $I = 1,05 \times 10^{38}$ kg m$^2$ (68% HPD: $[9,04, 12,4] \times 10^{37}$ kg m$^2$). Il valore di riferimento ($8,01 \times 10^{37}$ kg m$^2$) è contenuto nell'intervallo HPD al 95%. Il momento d'inerzia normalizzato $I/(MR_\oplus^2)$ è calcolato come $\approx 0,357$, comparabile al riferimento PREM di 0,331.

Significatività e Rivendicazioni
Il documento afferma di presentare la determinazione più precisa dell'interno della Terra basata sull'interazione debole ad oggi. Utilizzando un ordine di grandezza di dati in più rispetto ai tentativi precedenti e implementando un trattamento significativamente aggiornato delle incertezze sistematiche, lo studio dimostra che l'attenuazione dei neutrini può servire come sonda valida e complementare alla sismologia e alla gravimetria.

I risultati confermano che il profilo di densità della Terra, misurato tramite l'attenuazione dei neutrini, è coerente con i modelli stabiliti derivati dalla gravità e dalle onde sismiche. Gli autori sottolineano che questo lavoro stabilisce i neutrini come uno strumento fisico distinto per la tomografia planetaria. Sebbene le incertezze attuali siano dominate dalla modellazione del flusso, dalla conoscenza della sezione d'urto e dalle sistematiche del rivelatore/ghiaccio, il documento postula che i futuri miglioramenti nella calibrazione dei rivelatori (ad esempio, l'IceCube Upgrade) e i prossimi grandi telescopi trasformeranno questa tecnica da un controllo incrociato complementare a una sonda di precisione della struttura planetaria.