Reassessing the gallium anomaly using self-consistent electron wave functions

Cet article réévalue l'anomalie du gallium de longue date en calculant les sections efficaces de capture de neutrinos à l'aide de fonctions d'onde électroniques autocohérentes dérivées de l'équation de Dirac-Coulomb, actualisant ainsi la signification globale de la divergence et son interprétation concernant les neutrinos stériles.

L'essentiel

Imaginez que vous essayez de compter combien de gouttes de pluie frappent un seau spécifique. Vous avez une formule mathématique très précise qui prédit exactement combien de gouttes devraient frapper le seau en fonction de la taille de la pluie et de la taille du seau. Cependant, chaque fois que vous comptez réellement les gouttes dans la vraie vie, vous trouvez moins de gouttes que ce que votre formule prédisait. Cette pluie manquante est ce que les physiciens appellent l'« anomalie du Gallium ».

Pendant plus de 30 ans, des expériences utilisant du Gallium-71 (un type de métal) comme « seau » pour capturer des neutrinos (de minuscules particules fantomatiques provenant du soleil ou de sources radioactives) ont systématiquement trouvé moins de neutrinos que prévu. L'écart entre la prédiction et la réalité s'est tellement creusé qu'il est désormais considéré comme un grand mystère de la physique.

Cet article de Cadeddu et ses collègues est comme une équipe de mécaniciens experts décidant de reconstruire entièrement le moteur de cette formule de prédiction pour voir s'ils ont commis une erreur de calcul.

L'Ancien Moteur vs Le Nouveau Moteur

L'Ancien Moteur (L'Approximation) :
Auparavant, les scientifiques calculaient comment les neutrinos interagissent avec les atomes de Gallium en utilisant une version « brouillon » de la mathématique. Ils traitaens les électrons (les minuscules particules orbitant autour du noyau de l'atome) comme s'ils étaient des ondes simples et lisses qui ne changeaient pas beaucoup à l'intérieur de l'atome. C'était comme estimer la forme d'une route accidentée en regardant simplement une carte plate. Ils supposaient que l'onde de l'électron était la même partout à l'intérieur du minuscule noyau.

Le Nouveau Moteur (La Solution Exacte) :
Dans cette nouvelle étude, les auteurs ont décidé d'arrêter d'utiliser une carte plate. Au lieu de cela, ils ont utilisé un GPS haute définition pour résoudre les équations exactes (appelées l'équation de Dirac-Coulomb) qui décrivent comment les électrons se comportent réellement.

• L'Analogie : Imaginez que le noyau est une piste de danse bondée. L'ancienne méthode supposait que tout le monde sur la piste de danse se tenait immobile en un cercle parfait. La nouvelle méthode compte réellement chaque danseur, en tenant compte de la façon dont ils se cognent les uns aux autres et se déplacent dans l'espace bondé. Ils ont résolu les mathématiques pour les électrons coincés à l'atome et pour ceux qui s'envolent, en utilisant un programme informatique spécialisé pour obtenir la forme exacte de l'« onde » de l'électron.

Le Truc de l'« Moyenne »

Un autre changement clé dans cet article est la manière dont ils gèrent la taille du noyau.

• L'Ancienne Méthode : Ils choisissaient un point unique au centre du noyau (comme mesurer la température d'une pièce en enfonçant un thermomètre exactement au milieu) et supposaient que cela représentait l'ensemble.
• La Nouvelle Méthode : Ils ont réalisé que le noyau a une taille, ils ont donc « moyenné » le comportement de l'électron sur l'ensemble du volume du noyau. C'est comme prendre la température à chaque endroit de la pièce et trouver la véritable moyenne, plutôt que de simplement deviner à partir du centre.

Qu'ont-ils trouvé ?

Lorsqu'ils ont lancé leurs nouveaux calculs, plus précis :

1. La Prédiction a Changé : Leur nouvelle formule, plus précise, prédisait que moins de neutrinos devraient être capturés que l'ancienne formule.
2. L'Écart s'est Creusé : Parce que leur nouvelle prédiction est plus basse, la différence entre ce que nous attendons de voir et ce que les expériences ont réellement vu est devenue encore plus grande.
3. Le Résultat : Les « neutrinos manquants » sont désormais un problème de 5,5 sigma. Dans le monde de la science, le « sigma » est une mesure de confiance. Un résultat de 5 sigma est l'étalon-or d'une découverte, signifiant qu'il y a moins d'une chance sur 3,5 millions que cette divergence soit simplement un coup de chance aléatoire.

L'Hypothèse du « Neutrino Stérile »

Les physiciens ont une théorie favorite pour expliquer cette pluie manquante : les Neutrinos Stériles.

• La Métaphore : Imaginez que les neutrinos sont comme un vol d'oiseaux volant vers le seau. La théorie suggère que certains de ces oiseaux sont « invisibles » (stériles). Ils n'interagissent pas du tout avec le seau ; ils passent simplement à travers. Si ces oiseaux invisibles existent, ils expliqueraient pourquoi le seau est plus vide que prévu.

Les auteurs ont mis à jour les mathématiques pour voir si cette théorie de l'« oiseau invisible » tient toujours la route.

• La Bonne Nouvelle : Les mathématiques permettent toujours l'existence de ces oiseaux invisibles. Les données du Gallium pointent toujours fortement vers leur présence.
• La Mauvaise Nouvelle : D'autres expériences (observant les neutrinos de réacteurs, les neutrinos solaires et les mesures de masse) ont installé des « clôtures » qui disent que ces oiseaux invisibles ne devraient pas pouvoir voler de la manière dont les données du Gallium le suggèrent. Les données du Gallium veulent que les oiseaux soient très actifs, mais les autres clôtures disent qu'ils doivent être très timides.

La Conclusion

Les auteurs n'ont pas résolu le mystère ; ils l'ont rendu plus mystérieux. En utilisant une meilleure mathématique et des modèles d'électrons plus précis, ils ont confirmé que les « neutrinos manquants » sont un problème réel et robuste, et non une erreur de calcul.

Ils concluent que, bien que l'idée du « Neutrino Stérile » soit toujours le principal suspect, elle est actuellement dans un bras de fer avec d'autres preuves expérimentales. Le mystère reste entier, et les auteurs suggèrent qu'une nouvelle expérience avec un type de détecteur différent pourrait être nécessaire pour enfin attraper le coupable.

Résumé technique

Résumé technique : Réévaluation de l'anomalie du gallium à l'aide de fonctions d'onde électroniques auto-cohérentes

Énoncé du problème
L'article traite de l'« anomalie du gallium », un écart persistant observé dans les expériences de capture de neutrinos utilisant le $^{71}\text{Ga}$ comme cible. Les expériences employant des sources radioactives de $^{51}\text{Cr}$ et de $^{37}\text{Ar}$ (GALLEX, SAGE, et récemment BEST) ont systématiquement mesuré des taux de désintégration bêta inverse (IBD) induite par neutrinos plus faibles que les prédictions théoriques. Ce déficit, qui dépasse désormais une signification de $4\sigma$, remet en question la robustesse des calculs de section efficace théoriques et a alimenté des hypothèses concernant une physique au-delà du Modèle Standard, telles que les neutrinos stériles. Les auteurs identent que les évaluations théoriques précédentes reposaient sur des approximations d'ordre de grandeur (LO) pour les fonctions d'onde électroniques et sur des traitements incohérents des éléments de matrice nucléaires et leptoniques.

Méthodologie
Les auteurs entreprennent une révision complète de la section efficace de capture de neutrinos sur le $^{71}\text{Ga}$ en abandonnant les approximations conventionnelles de type LO au profit d'une approche numérique auto-cohérente. Les composantes méthodologiques clés incluent :

1. Solutions exactes des fonctions d'onde : Au lieu d'utiliser des fonctions d'onde approximatives, les auteurs résolvent numériquement l'équation de Dirac-Coulomb pour les états électroniques liés (capture électronique, EC) et de continuum (IBD). Ceci est réalisé à l'aide d'un progiciel personnalisé basé sur la bibliothèque RADIAL, employant les équations de Dirac-Hartree-Fock-Slater (DHFS).
2. Potentiels auto-cohérents : Les calculs utilisent un potentiel DHFS qui inclut les potentiels nucléaire, électronique et d'échange. Le potentiel nucléaire est modélisé à l'aide d'une distribution de Fermi à deux paramètres, incorporant la première mesure directe du rayon de charge nucléaire du $^{71}\text{Ge}$ ($R_{ch} = 4,032 \pm 0,002$ fm). Le paramètre du potentiel d'échange ($C_{ex}$) est moyenné à $1,25$ pour tenir compte des variations de modèle entre les approximations de Slater ($3/2$) et de Kohn-Sham ($1$).
3. Procédure de moyennage : Dépassant la convention consistant à évaluer les fonctions d'onde en un point unique (par exemple, $r_0 = 0$ ou $r_0 = R_{box}$), les auteurs proposent de moyenner les fonctions d'onde électroniques et la fonction de Fermi sur la densité de charge nucléaire ($\rho_{ch}(r)$). Cette méthode est appliquée tant au calcul de la section efficace d'IBD qu'au taux de capture électronique (EC) utilisé pour extraire l'élément de matrice nucléaire via le principe du bilan détaillé.
4. Données d'entrée actualisées : L'analyse incorpore des données expérimentales récentes, notamment :
   • Une moyenne pondérée de quatre mesures indépendantes de la demi-vie de capture électronique du $^{71}\text{Ge}$ ($t_{1/2} = 11,465 \pm 0,003$ j).
   • Une valeur Q actualisée pour le processus EC ($232,47 \pm 0,09$ keV).
   • Des rapports de probabilité de capture électronique affinés ($P_L/P_K$ et $P_M/P_K$) ainsi que les énergies de liaison.
   • Les forces Gamow-Teller pour les états excités dérivées de données de diffusion $(p,n)$ et $(^3\text{He}, ^3\text{H})$.

Contributions clés et résultats

• Sections efficaces révisées : Le nouveau calcul produit une section efficace de l'état fondamental environ 2,8 % plus faible que les estimations précédentes en raison de la procédure de moyennage. En incluant les contributions des états excités, les sections efficaces totales pour les sources de $^{51}\text{Cr}$ et de $^{37}\text{Ar}$ sont environ 2 % plus petites que celles rapportées dans la littérature la plus récente (Réf. [27]).
• Signification de l'anomalie : En comparant les prédictions théoriques actualisées avec les données expérimentales de GALLEX, SAGE et BEST, les auteurs réévaluent la signification de l'anomalie :
  • En utilisant les données de diffusion $(p,n)$ pour les contributions des états excités, le écart atteint $5,5\sigma$ ($R = 0,835^{+0,030}_{-0,048}$).
  • En utilisant les données de diffusion $(^3\text{He}, ^3\text{H})$, le écart est de $4,7\sigma$ ($R = 0,811^{+0,037}_{-0,035}$).
  • Une estimation conservatrice utilisant uniquement la section efficace de l'état fondamental produit une anomalie de $3,8\sigma$.
• Auto-cohérence : Les auteurs soulignent qu'en utilisant le même logiciel et les mêmes hypothèses sous-jacentes pour la fonction de Fermi (IBD) et la densité électronique au noyau (EC), les décalages systématiques s'annulent dans le rapport, rendant le résultat plus robuste.

Signification et interprétation
L'article affirme que l'anomalie du gallium est désormais plus significative qu'indiqué précédemment, atteignant jusqu'à $5,5\sigma$. Les auteurs réexaminent l'interprétation de cette anomalie en termes d'oscillations de neutrinos actifs-stériles à courte base (modèle 3+1). Leurs résultats confirment que les données du gallium, incluant les mesures précises de BEST, nécessitent un grand angle de mélange actif-stérile ($0,16 \lesssim \sin^2 2\vartheta_{ee} \lesssim 0,54$ à $2\sigma$ pour $\Delta m^2_{41} \gtrsim 0,98 \text{ eV}^2$).

Cependant, les auteurs notent que ce grand mélange requis est en tension avec les contraintes des données de neutrinos antineutrinos de réacteurs, les limites des neutrinos solaires et la mesure de la masse du neutrino par KATRIN. Par conséquent, l'article conclut que bien que la signification statistique de l'anomalie ait augmenté, l'explication par les neutrinos stériles reste en tension avec les autres données expérimentales. Les auteurs déclarent que la résolution de l'anomalie du gallium demeure un problème ouvert, nécessitant potentiellement une nouvelle expérience de source avec un détecteur différent pour résoudre la discrepancy.