On the Nonlinear Dependence of Underground Muon Rate on Atmospheric Temperature Observed at Daya Bay

Sintesi Tecnica: Sulla dipendenza non lineare del tasso di muoni sotterranei dalla temperatura atmosferica osservata a Daya Bay

Definizione del problema
La modulazione del flusso di muoni cosmici sotterranei da parte della temperatura atmosferica è un fenomeno ben consolidato, tradizionalmente spiegato dalle teorie di Barrett, Gaisser e altri. Queste teorie postulano una correlazione positiva tra il tasso di muoni e la temperatura atmosferica, guidata dalla competizione tra il decadimento (pioni e kaoni) e l'interazione dei mesoni. Tuttavia, i dati dell'esperimento di neutrini di Daya Bay, in particolare nelle sale sperimentali meno profonde EH1 ed EH2, mostrano un effetto di modulazione non lineare: la correlazione appare più forte durante i mesi caldi rispetto a quelli freddi. Gli attuali quadri teorici, che spesso si basano su assunzioni di atmosfera isotermica o su proiezioni asintotiche di perturbazioni di temperatura su pesi di profondità monodimensionali, non riescono a spiegare questa non linearità quando applicati a profili di temperatura realistici e non isotermi.

Metodologia
Gli autori propongono un quadro teorico rivisto e una strategia di verifica numerica:

1. Soluzione analitica migliorata: Il documento deriva una soluzione migliorata delle equazioni di evoluzione del cascata atmosferica per gli spettri energetici di mesoni e muoni. A differenza delle precedenti approssimazioni che collegano i limiti di alta e bassa energia o assumono condizioni isoterme, questa soluzione mantiene la piena dipendenza non locale dal profilo di temperatura atmosferica. Nello specifico, tratta la lunghezza di decadimento come un funzionale esplicito della relazione densità-temperatura, preservando l'integrale della storia della temperatura lungo la traiettoria del mesone dal momento della produzione al decadimento.
2. Nuove definizioni: Sulla base di questa soluzione migliorata, gli autori ridefiniscono la temperatura effettiva ($T_{eff}$), la sua variazione ($\Delta T_{eff}$) e il coefficiente di temperatura ($\alpha_T$). Invece di utilizzare quantità isoterme proiettate, queste sono definite tramite la derivata funzionale del tasso di muoni rispetto al profilo di temperatura reale e variabile nel tempo $T(Y)$. Questo approccio decompone la risposta in un termine locale (probabilità di decadimento allo strato di produzione) e un termine di propagazione (storia di sopravvivenza del mesone).
3. Simulazione numerica (MCEq): Per validare la teoria, gli autori utilizzano lo strumento Matrix Cascade Equations (MCEq). Inseriscono dati atmosferici reali dal database ERA5 (coprendo il periodo 2011–2021) e dati specifici del terreno per le sale sperimentali di Daya Bay (EH1, EH2, EH3).
   • La simulazione tiene conto del flusso di raggi cosmici primari (modello Hillas Gaisser 2012), delle interazioni adroniche (SIBYLL2.3C) e del sovraccarico di roccia (overburden).
   • Lo studio riproduce le variazioni del tasso di muoni nell'arco di 11 anni.
   • Due percorsi di calcolo vengono confrontati: uno utilizzando il convenzionale metodo di pesatura di Grashorn (che proietta le perturbazioni su un peso di profondità 1D) e uno utilizzando il nuovo metodo della derivata funzionale (calcolato numericamente tramite differenze finite in MCEq).

Risultati chiave

• Riproduzione della non linearità: Utilizzando le soluzioni di cascata convenzionali e il metodo standard di pesatura di Grashorn, gli autori riescono a riprodurre la correlazione non lineare tra la variazione relativa del tasso di muoni sotterranei e la variazione relativa della temperatura effettiva osservata a Daya Bay. Si osserva che il coefficiente di temperatura $\alpha_T$ è maggiore alle temperature più elevate.
• Ripristino della linearità: Quando la variazione del tasso di muoni viene calcolata utilizzando i nuovi pesi basati sulla derivata funzionale (che tengono conto della completa storia non locale della temperatura), la relazione tra il tasso di muoni e la temperatura effettiva diventa lineare.
• Validazione: Il calcolo numerico della derivata funzionale tramite MCEq conferma le previsioni teoriche derivate nella Sezione 3. La dipendenza lineare è recuperata attraverso diversi modelli di interazione (SIBYLL2.3C, 2.1, 2.3E) e soglie energetiche.

Significatività e affermazioni
Il documento afferma che la non linearità osservata a Daya Bay non è un breakdown fondamentale della fisica della correlazione temperatura-muoni, ma piuttosto un artefatto delle approssimazioni utilizzate nelle definizioni convenzionali di temperatura effettiva. Trattando il profilo di temperatura come un input funzionale anziché proiettarlo su un peso semplificato, gli autori dimostrano che la relazione fisica sottostante è lineare.

Il lavoro fornisce:

1. Una soluzione più generale delle equazioni di cascata che descrive pienamente come l'intero profilo di temperatura influenzi il tasso finale di muoni.
2. Definizioni raffinate per la temperatura effettiva e il coefficiente di temperatura che sono valide per atmosfere non isoterme.
3. Un quadro di calcolo che può spiegare gli effetti non lineari riscontrati a Daya Bay e offre un metodo più preciso per calcolare i coefficienti di temperatura per altri esperimenti sotterranei.

Gli autori concludono che il loro metodo recupera la prevista dipendenza lineare, suggerendo che le precedenti interpretazioni della modulazione non lineare erano limitate dalle assunzioni di proiezione isotermica inerenti ai trattamenti analitici standard.