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

Résumé technique : Sur la dépendance non linéaire du taux de muons souterrains par rapport à la température atmosphérique observée à Daya Bay

Énoncé du problème
La modulation du flux de muons de rayons cosmiques souterrains par la température atmosphérique est un phénomène bien établi, traditionnellement expliqué par les théories de Barrett, Gaisser et autres. Ces théories postulent une corrélation positive entre le taux de muons et la température atmosphérique, pilotée par la compétition entre la désintégration des mésons (pions et kaons) et leur interaction. Cependant, les données de l'expérience Daya Bay, particulièrement dans les salles expérimentales moins profondes 1 et 2, présentent un effet de modulation non linéaire : la corrélation semble plus forte durant les mois chauds que durant les mois froids. Les cadres théoriques existants, qui reposent souvent sur des hypothèses d'atmosphère isotherme ou sur des projections asymptotiques de perturbations de température sur des poids de profondeur unidimensionnels, ne parviennent pas à rendre compte de cette non linéarité lorsqu'ils sont appliqués à des profils de température réalistes et non isothermes.

Méthodologie
Les auteurs proposent un cadre théorique révisé et une stratégie de vérification numérique :

1. Solution analytique améliorée : L'article dérive une solution améliorée des équations d'évolution de la cascade atmosphérique pour les spectres d'énergie des mésons et des muons. Contraãirement aux approximations précédentes qui lient les limites de haute et basse énergie ou supposent des conditions isothermes, cette solution conserve la pleine dépendance non locale du profil de température atmosphérique. Plus précisément, elle traite la longueur de désintégration comme une fonction explicite de la relation densité-température, préservant l'intégrale de l'historique de température le long de la trajectoire du méson, de la production à la désintégration.
2. Nouvelles définitions : Sur la base de cette solution améliorée, les auteurs redéfinissent la température effective ($T_{eff}$), sa variation ($\Delta T_{eff}$) et le coefficient de température ($\alpha_T$). Au lieu d'utiliser des quantités projetées isothermes, celles-ci sont définies via la dérivée fonctionnelle du taux de muons par rapport au profil de température réel et variable dans le temps $T(Y)$. Cette approche décompose la réponse en un terme local (probabilité de désintégration au niveau de la couche de production) et un terme de propagation (historique de survie du méson).
3. Simulation numérique (MCEq) : Pour valider la théorie, les auteurs utilisent l'outil Matrix Cascade Equations (MCEq). Ils injectent des données atmosphériques réelles provenant de la base de données ERA5 (couvrant la période 2011–2021) ainsi que des données de terrain spécifiques pour les salles expérimentales de Daya Bay (EH1, EH2, EH3).
   • La simulation prend en compte le flux de rayons cosmiques primaires (modèle Hillas Gaisser 2012), les interactions hadroniques (SIBYLL2.3C) et la couverture rocheuse.
   • L'étude reproduit les variations du taux de muons sur 11 ans.
   • Deux voies de calcul sont comparées : l'une utilisant la pondération conventionnelle de type Grashorn (qui projette les perturbations sur un poids de profondeur 1D) et l'autre utilisant la nouvelle méthode de dérivée fonctionnelle (calculée numériquement par différences finies dans MCEq).

Résultats clés

• Reproduction de la non linéarité : En utilisant les solutions de cascade conventionnelles et la méthode de pondération standard de Grashorn, les auteurs parviennent à reproduire la corrélation non linéaire entre la variation relative du taux de muons souterrains et la variation relative de la température effective observée à Daya Bay. Il est démontré que le coefficient de température $\alpha_T$ est plus élevé à des températures plus hautes.
• Restauration de la linéarité : Lorsque la variation du taux de muons est calculée en utilisant les nouveaux poids basés sur la dérivée fonctionnelle (qui tiennent compte de l'historique complet de la température non locale), la relation entre le taux de muons et la température effective devient linéaire.
• Validation : Le calcul numérique de la dérivée fonctionnelle à l'aide de MCEq confirme les prédictions théoriques dérivées de la Section 3. La dépendance linéaire est récupérée à travers différents modèles d'interaction (SIBYLL2.3C, 2.1, 2.3E) et seuils d'énergie.

Signification et revendications
L'article affirme que la non linéarité observée à Daya Bay n'est pas une rupture fondamentale de la physique de la corrélation température-muon, mais plutôt un artefact des approximations utilisées dans les définitions conventionnelles de la température effective. En traitant le profil de température comme une entrée fonctionnelle plutôt qu'en le projetant sur un poids simplifié, les auteurs démontrent que la relation physique sous-jacente est linéaire.

Ce travail fournit :

1. Une solution plus générale des équations de cascade qui décrit pleinement comment l'ensemble du profil de température influence le taux de muons final.
2. Des définitions raffinées de la température effective et du coefficient de température qui sont valables pour des atmosphères non isothermes.
3. Un cadre de calcul capable d'expliquer les effets non linéaires trouvés à Daya Bay et offrant une méthode plus précise pour calculer les coefficients de température pour d'autres expériences souterraines.

Les auteurs concluent que leur méthode restaure la dépendance linéaire attendue, suggérant que les interprétations précédentes de la modulation non linéaire étaient limitées par les hypothèses de projection isotherme inhérentes aux traitements analytiques standards.