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

Résumé technique : Tomographie de l'intérieur de la Terre par les neutrinos de haute énergie avec IceCube

Problème et motivation
La structure interne de la Terre, plus précisément sa distribution de densité radiale, est fondamentale pour comprendre l'évolution géologique, la géodynamo et la formation planétaire. Historiquement, ces propriétés ont été déduites principalement par la sismologie (analyse de la propagation des ondes sismiques) et la gravimétrie, qui sondent la réponse macroscopique de la matière aux forces élastiques et gravitationnelles. Bien que le modèle de référence préliminaire de la Terre (PREM) constitue un standard robuste basé sur ces méthodes, il repose sur des inférences indirectes. Une approche physique distincte consiste à sonder la densité de la Terre en utilisant l'interaction faible. Les neutrinos de haute énergie, qui n'interagissent que via la force faible, traversent la Terre mais subissent une atténuation dépendant de la densité de colonne de matière qu'ils rencontrent. À mesure que l'énergie des neutrinos augmente, la section efficace d'interaction croît, rendant la Terre partiellement opaque aux neutrinos au-dessus de $\sim$10 TeV. Cette atténuation est sensible au profil de densité le long de la trajectoire du neutrino, offrant ainsi une sonde nouvelle et indépendante des intérieurs planétaires.

Méthodologie
Cette étude utilise 10,7 années de données (2011–2022) de l'observatoire de neutrinos IceCube au pôle Sud. L'analyse se concentre sur 368 071 événements de neutrinos muoniques ascendants dans la gamme d'énergie de 500 GeV à 100 TeV. Ces événements comprennent à la fois des neutrinos atmosphériques (produits par les interactions des rayons cosmiques dans l'atmosphère) et un flux astrophysique diffus subdominant.

La méthodologie procède comme suit :

1. Sélection des événements : Le jeu de données comprend des traces de muons « de départ » (starting) et « traversantes » (throughgoing). Les traces de départ (interactions à l'intérieur du volume instrumenté) offrent une reconstruction d'énergie supérieure, tandis que les traces traversantes (interactions à l'extérieur du détecteur) offrent des statistiques élevées.
2. Modélisation : Le profil de densité de la Terre est paramétré comme un ensemble de couches concentriques à densité uniforme. L'étude utilise principalement une configuration à cinq couches (noyau interne, noyau externe et couches manteau/croûte) et valide les résultats avec une extension à huit couches.
3. Simulation et vraisemblance : Des simulations de Monte Carlo modélisent l'acceptation du détecteur IceCube, les incertitudes systématiques (y compris les propriétés de la glace glaciaire) et le flux de neutrinos attendu. L'analyse tient compte de la section efficace neutrino-nucléon du Modèle Standard, qui est bien contrainte par les fonctions de distribution de partons (PDF) et les données de diffusion inélastique profonde.
4. Inférence : Un ajustement conjoint est effectué à l'aide d'une vraisemblance de Poisson groupée pour contraindre les densités des couches, la normalisation du flux, les incertitudes de la section efficace et les paramètres de nuisance du détecteur. Les distributions postérieures sont échantillonnées à l'aide d'un algorithme de Monte Carlo par chaîne de Markov (MCMC).
5. Analyse en aveugle : Pour éviter les biais, l'analyse a été menée en deux étapes : un ajustement en aveugle excluant les événements proches de l'horizon ($-1,0 < \cos\theta < -0,1$), où les profondeurs de colonne sont les plus faibles et où la modélisation erronée systématique est la plus probable, suivi d'un ajustement utilisant toutes les données disponibles ($-1,0 < \cos\theta < 0,0$).

Résultats clés
L'analyse produit les résultats primaires suivants :

• Profil de densité : Le profil de densité radiale inféré, dérivé du modèle à cinq couches, est cohérent avec le modèle PREM dans les intervalles de crédibilité de la plus haute densité postérieure (HPD) de 68 %. Les résultats ne montrent aucune déviation significative par rapport au profil PREM, validant la sonde par interaction faible.
  • Couche 1 (Noyau interne) : Mode entre 0,60 et 0,73 fois la densité PREM.
  • Couches 2 et 3 (Noyau externe) : Les modes varient de 0,72 à 1,1 fois la densité PREM.
  • Couches 4 et 5 (Manteau/Croûte) : Les modes varient de 1,3 à 1,6 fois la densité PREM.
  • L'inclusion des événements proches de l'horizon déplace légèrement le mode postérieur de la couche la plus externe (Couche 5) vers l'unité, suggérant que ces événements contribuent de manière constructive aux contraintes.
• Quantités globales dérivées :
  • Masse de la Terre : Le mode postérieur pour la masse de la Terre est $M = 7,25 \times 10^{24}$ kg, avec un intervalle HPD de 68 % de $[6,31, 8,31] \times 10^{24}$ kg. Bien que le mode soit supérieur à la valeur de référence gravitationnelle ($5,97 \times 10^{24}$ kg), la valeur de référence se situe dans l'intervalle HPD de 95 %. L'hypothèse de masse nulle est exclue à plus de $5\sigma$.
  • Moment d'inertie polaire : Le mode postérieur est $I = 1,05 \times 10^{38}$ kg m$^2$ (HPD 68 % : $[9,04, 12,4] \times 10^{37}$ kg m$^2$). La valeur de référence ($8,01 \times 10^{37}$ kg m$^2$) est contenue dans l'intervalle HPD de 95 %. Le moment d'inertie normalisé $I/(MR_\oplus^2)$ est calculé à environ 0,357, comparable à la référence PREM de 0,331.

Signification et affirmations
L'article affirme présenter la détermination la plus précise de l'intérieur de la Terre basée sur l'interaction faible à ce jour. En utilisant un ordre de grandeur de données supplémentaires par rapport aux tentatives précédentes et en mettant en œuvre un traitement des incertitudes systématiques nettement mis à jour, l'étude démontre que l'atténuation des neutrinos peut servir de sonde viable et complémentaire à la sismologie et à la gravimétrie.

Les résultats confirment que le profil de densité de la Terre, tel que mesuré via l'atténuation des neutrinos, est cohérent avec les modèles établis dérivés de la gravité et des ondes sismiques. Les auteurs soulignent que ce travail établit les neutrinos comme un outil physique distinct pour la tomographie planétaire. Bien que les incertitudes actuelles soient dominées par la modélisation du flux, la connaissance de la section efficace et les systématiques du détecteur/glace, l'article postule que les améliorations futures de la calibration des détecteurs (par exemple, l'IceCube Upgrade) et les prochains grands télescopes de nouvelle génération transformeront cette technique, passant d'une vérification complémentaire à une sonde de précision de la structure planétaire.