Modeling the Effect of the Heliospheric Magnetic Field on Cosmic Ray Muon Shadows

Sintesi Tecnica: Modellazione dell'Effetto del Campo Magnetico Eliosferico sulle Ombre dei Muoni dei Raggi Cosmici

Enunciato del Problema
I raggi cosmici (CR), costituiti principalmente da protoni e particelle alfa, vengono deflessi dai campi elettromagnetici mentre attraversano il sistema solare. Sebbene il campo magnetico terrestre (GMF) sia ben mappato, il Campo Magnetico Eliosferico (HMF) è complesso e difficile da modellare. Il Sole e la Luna assorbono i CR incidenti, creando "ombre" nel flusso dei raggi cosmici osservabili dalla Terra. Tuttavia, queste ombre non sono statiche; vengono spostate e distorte dall'HMF e dal GMF. La Collaborazione MINOS+ ha utilizzato dati sui muoni provenienti dal rivelatore lontano MINOS per localizzare queste ombre solari durante diverse fasi del ciclo solare (minimo solare, massimo solare e l'intero arco di 13 anni). Il problema centrale affrontato è se gli spostamenti e le distorsioni osservati dell'ombra solare dei muoni siano coerenti con il modello standard a spirale di Parker dell'HMF combinato con un modello dipolare del GMF, e in caso contrario, quali modifiche ai modelli o ai parametri siano necessarie per conciliare la teoria con l'osservazione.

Metodologia
Lo studio ha utilizzato 13 anni di dati sui muoni (luglio 2003–luglio 2016) provenienti dal rivelatore lontano MINOS, situato a 2070 metri equivalenti d'acqua sottoterra nella miniera di Soudan. L'analisi si è concentrata su eventi di singoli muoni con angoli zenitali incidenti inferiori a 80° e lunghezza di traccia sufficiente per ricostruire le traiettorie.

1. Localizzazione dell'Ombra: Il team ha impiegato un test di probabilità $\Lambda$ (un rapporto di verosimiglianza logaritmica) per localizzare l'ombra. Ciò ha comportato il confronto dei conteggi di particelle osservati ($N_{obs}$) vicino al Sole con un fondo teorico piatto ($N_{back}$) e un modello teorico di ombra ($N_{th}$). Il modello teneva conto delle dimensioni angolari del Sole (0,52°) e degli effetti di scattering nel sovraccarico roccioso.
2. Validazione Statistica: Sono state utilizzate simulazioni Monte Carlo per generare istogrammi di fondo casuali al fine di determinare la significatività delle posizioni osservate delle ombre. È stata calcolata la probabilità che fluttuazioni casuali mimassero l'ombra, ottenendo livelli di confidenza elevati (fino a 6$\sigma$ per l'intero set di dati).
3. Modellazione Computazionale: Per interpretare gli spostamenti delle ombre, gli autori hanno simulato le traiettorie a tempo inverso dei CR dal rivelatore fino al Sole.
   • Modelli di Campo: La simulazione ha sovrapposto il modello a spirale di Parker per l'HMF e un modello dipolare per il GMF. I parametri della spirale di Parker includevano una velocità costante del vento solare ($u_{sowi} = 400$ km/s) e la frequenza di rotazione solare. Il GMF è stato modellato utilizzando un momento di dipolo e convertito in un sistema di coordinate centrato sul Sole.
   • Dinamica delle Particelle: L'equazione della forza di Lorentz è stata integrata utilizzando uno schema Runge-Kutta del 4° ordine. Le particelle sono state trattate come relativistiche, con energie che seguono una distribuzione a legge di potenza ($E^{-2.75}$).
   • Scenari: Le simulazioni hanno testato vari rapporti di polarità del campo magnetico (50/50, 33/67, 67/33) e due tipi di modelli di ombra: un'ombra circolare semplice e un'ombra "distorta" derivata dagli effettivi istogrammi di probabilità $\Lambda$.
   • Variazione dei Parametri: Il taglio energetico primario dei raggi cosmici ($E_0$) è stato trattato come parametro libero per determinare se uno spettro energetico diverso potesse conciliare il modello con le posizioni osservate delle ombre.

Contributi Chiave

• Caratterizzazione dell'Ombra: Il documento ha localizzato con successo le ombre solari dei muoni con alta significatività statistica in tre distinti intervalli temporali: l'intero set di dati di 13 anni, il minimo solare (2007–2009) e il massimo solare (2012–2014).
• Sviluppo di Modelli: Gli autori hanno introdotto un metodo per creare modelli di campionamento di ombre "distorte" basati sulle distribuzioni di probabilità $\Lambda$ osservate, andando oltre l'assunzione di un'ombra perfettamente circolare.
• Test del Modello: Lo studio ha testato rigorosamente il modello a spirale di Parker contro i dati osservativi, quantificando l'accordo tra gli spostamenti delle ombre simulati e osservati in Ascensione Retta ($\Delta RA$) e Declinazione ($\Delta DEC$).
• Analisi dello Spettro Energetico: Il documento ha indagato la sensibilità delle posizioni delle ombre allo spettro energetico primario dei CR, testando specificamente l'impatto del taglio a bassa energia ($E_0$).

Risultati

• Spostamenti dell'Ombra: Le ombre osservate hanno mostrato spostamenti distinti a seconda della fase del ciclo solare. L'ombra del minimo solare è stata spostata diversamente rispetto all'ombra del massimo solare, e l'intero set di dati rappresenta una combinazione di questi spostamenti.
• Coerenza del Modello:
  • Il modello a spirale di Parker ha mostrato il migliore accordo con i dati del minimo solare.
  • Il modello ha mostrato il minore accordo con i dati del massimo solare, in particolare per quanto riguarda lo spostamento $\Delta DEC$.
  • I modelli di ombra distorta hanno funzionato meglio per i periodi di massimo solare e per l'intero set di dati, mentre i modelli rotondi semplici sono stati più efficaci per i dati del minimo solare.
• Effetti di Polarità: Le simulazioni hanno confermato che il rapporto di polarità del campo magnetico influisce direttamente sulla posizione $\Delta DEC$ dell'ombra. I dati sembravano favorire un HMF dominato da polarità positiva, specialmente durante il massimo solare.
• Discrepanza del Taglio Energetico: Per allineare le medie delle ombre simulate con il centro teorico (0,0) utilizzando il modello a spirale di Parker, il taglio energetico primario dei CR ($E_0$) ha dovuto essere aggiustato a 22,8 $\pm$ 1,2 TeV. Questo valore è significativamente più alto (2–3 volte) rispetto al valore sperimentalmente determinato di $\approx$10 TeV.
• Correzione Zenitale: È stata applicata una correzione per la dipendenza dell'angolo zenitale del taglio energetico. Sebbene ciò abbia aggiustato lo spettro per corrispondere meglio alla distribuzione energetica osservata dei muoni, è stato insufficiente a spiegare completamente lo spostamento $\Delta RA$ osservato.

Significato e Conclusioni
Il documento conclude che le ombre solari dei muoni dei raggi cosmici sono osservabili e possono essere utilizzate per sondare la struttura dell'HMF. I risultati indicano che l'ombra del minimo solare è la più coerente con il modello standard a spirale di Parker combinato con un GMF dipolare. Tuttavia, il modello non riesce a riprodurre completamente l'ombra del massimo solare e l'intero set di dati senza invocare un taglio energetico irrealisticamente alto.

Gli autori affermano che la discrepanza suggerisce che il modello a spirale di Parker è una semplificazione eccessiva per le condizioni presenti durante il massimo solare. Propongono che fattori non modellati, come l'inclinazione del foglio di corrente neutro (il "Modello Ballerina"), i fronti d'urto e le complesse variazioni di polarità, contribuiscano probabilmente agli spostamenti osservati. Lo studio stabilisce che, sebbene la direzione generale della deflessione sia coerente con la spirale di Parker, la struttura dettagliata dell'HMF durante il massimo solare richiede una modellazione più complessa di quella fornita dal semplice modello a spirale. Il lavoro dimostra l'utilità delle osservazioni delle ombre dei muoni per vincolare i modelli HMF, ma evidenzia la necessità di simulazioni più sofisticate per tenere conto delle variazioni del ciclo solare.