Evidence for Normal Neutrino Ordering with Combined Atmospheric Neutrino Experiments

Cet article étudie l’ordre des masses des neutrinos au moyen d’une analyse combinée de trois expériences consacrées aux neutrinos atmosphériques : Super-Kamiokande, IceCube-DeepCore et KM3NeT/ORCA6. L’objectif principal est de déterminer si le troisième état de masse est plus lourd que les deux autres, conformément à l’ordre normal, ou plus léger, conformément à l’ordre inversé. Les auteurs cherchent également à préciser les paramètres des oscillations atmosphériques — l’angle de mélange θ₂₃ et la différence de masses au carré Δm²₃₁ — ainsi qu’à contraindre la phase de violation de CP, δCP, indépendamment des expériences utilisant des faisceaux d’accélérateurs. L’analyse réunit les données de Super-Kamiokande antérieures à la phase au gadolinium, soit 6 511,3 jours d’observation et une exposition de 484,2 kilotonnes-années, huit années de données de DeepCore et 510 jours de données d’ORCA6, comprenant 5 828 événements sélectionnés. Ces détecteurs présentent des sensibilités complémentaires : Super-Kamiokande observe notamment des événements de basse énergie sensibles à δCP, tandis que DeepCore dispose d’un échantillon statistique important pour mesurer la disparition des neutrinos muoniques. L’ordre des masses est sondé grâce aux effets de la matière terrestre, qui modifient différemment les oscillations des neutrinos et des antineutrinos, particulièrement entre 2 et 10 GeV. L’innovation méthodologique consiste à ajuster simultanément les différents échantillons en considérant comme entièrement corrélées les incertitudes communes sur les flux atmosphériques et les sections efficaces d’interaction, tout en conservant les incertitudes propres à chaque détecteur. Les auteurs utilisent une vraisemblance de Poisson par intervalles, une simulation numérique de la propagation dans la Terre et 143 paramètres de nuisance. Pour Super-Kamiokande, ils développent une simulation Monte-Carlo correspondant à cent fois l’exposition des données et appliquent une correction fournie par la collaboration afin de remédier à un défaut des distributions publiques d’énergie vraie. Les paramètres solaires et θ₁₃ sont fixés à des valeurs externes : l’analyse est donc indépendante des faisceaux d’accélérateurs, mais ne s’affranchit pas de toute information extérieure. Sous l’hypothèse de l’ordre normal, l’ajustement donne Δm²₃₁ = 2,44 × 10⁻³ eV², avec des incertitudes de +0,06 et −0,04 × 10⁻³ eV², ainsi que sin²θ₂₃ = 0,52, avec des incertitudes de +0,03 et −0,04. Les précisions annoncées sont respectivement de 2,2 % et de 6,1 %. Le mélange maximal reste compatible avec les données, et la préférence pour l’octant supérieur est négligeable, avec une différence de χ² de 0,3. La valeur ajustée de δCP est de 1,30π, avec des incertitudes de +0,24π et −0,28π ; les valeurs proches de 0,3π sont défavorisées à plus de 2,7 écarts-types. Le résultat central est une préférence pour l’ordre normal, caractérisée par Δχ² = 11,0. Les deux hypothèses n’étant pas emboîtées, les auteurs évaluent la signification statistique de ce résultat à l’aide de pseudoexpériences. Parmi les 4 397 simulations réalisées sous l’hypothèse de l’ordre inversé, deux atteignent la région extrême correspondant au résultat observé, ce qui donne une valeur p de 4,55 × 10⁻⁴, soit 3,32 écarts-types. Le nombre limité de simulations entraîne une plage d’incertitude de 3,08 à 3,60 écarts-types à un niveau de confiance de 68 %. Le critère alternatif CLs donne, quant à lui, une signification de 2,55 écarts-types. Les auteurs concluent ainsi à un indice substantiel, mais non définitif, en faveur de l’ordre normal. Plusieurs limites invitent toutefois à la prudence : une tension de 2,2 écarts-types entre les expériences, l’absence d’un paramètre d’échelle d’énergie pour ORCA6 et la sensibilité de δCP aux incertitudes systématiques ainsi qu’aux fluctuations des données. La valeur p associée à la qualité globale de l’ajustement passe de 0,011 à environ 0,12 lorsque l’évaluation tient compte de la taille finie des échantillons Monte-Carlo. Enfin, une estimation de la combinaison avec JUNO suggère une signification proche de quatre écarts-types, mais cette perspective devra être confirmée par une analyse conjointe assortie d’un étalonnage statistique adapté.

Neutrinos : trois expériences renforcent la piste de l’ordre normal des masses

Les neutrinos traversent la matière presque sans laisser de trace, mais ils gardent un secret essentiel : dans quel ordre leurs masses sont-elles rangées ? Une analyse combinant trois expériences consacrées aux neutrinos atmosphériques apporte un nouvel indice en faveur de l’ordre dit « normal ». Le résultat est substantiel, sans constituer une réponse définitive.

Dans « Evidence for Normal Neutrino Ordering with Combined Atmospheric Neutrino Experiments », Miaochen Jin, Carlos A. Argüelles, Pablo Fernández-Menéndez et Iván Martínez-Soler réunissent les observations de Super-Kamiokande, IceCube-DeepCore et KM3NeT/ORCA6. Leur travail explore aussi les paramètres des oscillations des neutrinos, indépendamment des expériences utilisant des faisceaux d’accélérateurs.

Pourquoi l’ordre des masses compte

Les neutrinos existent sous trois saveurs : électronique, muonique et tauique. Lorsqu’ils se déplacent, ils peuvent passer d’une saveur à l’autre : c’est le phénomène d’oscillation. Ces transformations dépendent notamment de trois états de masse, qui ne correspondent pas directement aux trois saveurs.

Deux arrangements restent possibles. Dans l’ordre normal, le troisième état de masse est plus lourd que les deux autres. Dans l’ordre inversé, il est plus léger. Distinguer ces configurations constitue une étape importante pour comprendre les propriétés des neutrinos.

Les neutrinos atmosphériques, produits lorsque les rayons cosmiques interagissent avec l’atmosphère, offrent un laboratoire naturel. Certains traversent la Terre avant d’atteindre un détecteur. La matière terrestre modifie alors leurs oscillations différemment selon qu’il s’agit de neutrinos ou d’antineutrinos. Particulièrement utiles entre 2 et 10 GeV, ces effets permettent de sonder l’ordre des masses.

Trois détecteurs aux atouts complémentaires

L’étude rassemble des échantillons importants, mais différents par leur couverture énergétique et leur sensibilité :

• Super-Kamiokande apporte 6 511,3 jours d’observation avant sa phase au gadolinium, soit une exposition de 484,2 kilotonnes-années.
• IceCube-DeepCore fournit huit années de données et un grand échantillon pour étudier la disparition des neutrinos muoniques.
• KM3NeT/ORCA6 contribue avec 510 jours de données comprenant 5 828 événements sélectionnés.

Super-Kamiokande observe notamment des événements de basse énergie sensibles à la phase δCP. Ce paramètre est lié à une éventuelle violation de la symétrie CP dans les oscillations, c’est-à-dire à des différences entre le comportement des neutrinos et celui des antineutrinos. La complémentarité des instruments permet ainsi d’examiner plusieurs propriétés à la fois.

Une combinaison qui partage les incertitudes communes

L’innovation méthodologique consiste à ajuster simultanément les échantillons en traitant comme entièrement corrélées les incertitudes communes sur les flux atmosphériques et les sections efficaces d’interaction. Les incertitudes propres à chaque détecteur restent distinctes. Cette approche évite de traiter les mêmes sources d’incertitude comme indépendantes d’une expérience à l’autre.

L’analyse utilise une vraisemblance de Poisson par intervalles, une simulation de la propagation dans la Terre et 143 paramètres de nuisance. Pour Super-Kamiokande, les auteurs développent une simulation Monte-Carlo équivalant à cent fois l’exposition réelle, puis appliquent une correction fournie par la collaboration aux distributions publiques d’énergie vraie.

Les paramètres solaires et l’angle θ₁₃ sont fixés à des valeurs externes. L’étude est donc indépendante des faisceaux d’accélérateurs, mais pas de toute information extérieure.

Un indice significatif, pas une découverte

Le résultat central favorise l’ordre normal avec une différence de χ² de 11,0. Toutefois, les deux hypothèses ne sont pas emboîtées : convertir cette différence en signification statistique exige un étalonnage adapté. Les auteurs utilisent donc des pseudoexpériences pour déterminer à quel point leur résultat serait inhabituel si l’ordre inversé était correct.

Sur 4 397 simulations réalisées sous cette hypothèse, deux atteignent la région extrême correspondant au résultat observé. La valeur p obtenue est de 4,55 × 10⁻⁴, soit 3,32 écarts-types. Ce nombre ne représente pas la probabilité que l’ordre inversé soit vrai.

La quantité limitée de simulations entraîne une incertitude : à 68 % de confiance, la signification s’étend de 3,08 à 3,60 écarts-types. Le critère alternatif CLs donne une signification plus modérée, de 2,55 écarts-types. L’ensemble constitue donc un indice substantiel, mais non définitif.

Des paramètres mieux cernés, avec des réserves

Sous l’hypothèse de l’ordre normal, l’ajustement donne Δm²₃₁ = 2,44 × 10⁻³ eV² et sin²θ₂₃ = 0,52, avec des précisions annoncées de 2,2 % et 6,1 %. Le mélange maximal reste compatible avec les observations, et la préférence pour l’octant supérieur de θ₂₃ est négligeable.

La phase δCP est ajustée à 1,30π, avec des incertitudes de +0,24π et −0,28π. Les valeurs proches de 0,3π sont défavorisées à plus de 2,7 écarts-types, mais cette mesure reste sensible aux incertitudes systématiques et aux fluctuations des données.

D’autres réserves comptent : une tension de 2,2 écarts-types entre expériences et l’absence de paramètre d’échelle d’énergie pour ORCA6. La valeur p de qualité globale de l’ajustement passe toutefois de 0,011 à environ 0,12 lorsque la taille finie des simulations est prise en compte.

Enfin, une estimation associant ces résultats à JUNO suggère une sensibilité proche de quatre écarts-types. Cette perspective devra être confirmée par une véritable analyse conjointe et un étalonnage statistique adapté : la piste se renforce, mais l’enquête reste ouverte.