Testing Lepton Wave Function Factorization in $^{71}\text{Ge}$ Electron Capture

Résumé Technique : Test de la factorisation de la fonction d'onde leptonique dans la capture électronique du $^{71}$Ge

Énoncé du problème
Les détecteurs au germanium de haute pureté sont essentiels pour les recherches d'événements rares (neutrinos et matière noire), mais ils souffrent de l'activation de fond produisant l'isotope $^{71}$Ge. Cet isotope subit une capture électronique (CE) vers l'état fondamental du $^{71}$Ga, un processus qui sert à la fois de source de bruit de fond et d'outil de calibration. Une description théorique robuste de ce processus de CE est essentielle pour interpréter les données expérimentales. Alors que la théorie standard de la CE orbitale factorise généralement les fonctions d'onde leptoniques des éléments de matrice nucléaire (en les développant à l'ordre dominant au sein du volume nucléaire), des discussions récentes suggèrent que cette factorisation peut être insuffisante pour interpréter certaines observations expérimentales, particulièrement dans le contexte de « l'anomalie du gallium ». Cet article étudie si les effets de non-factorisation — spécifiquement l'interaction entre la densité de transition nucléaire et les fonctions d'onde leptoniques exactes — modifient de manière significative les taux de CE et les rapports prédits pour le $^{71}$Ge.

Méthodologie
Les auteurs développent un cadre général pour prédire les taux de CE du $^{71}$Ge en utilisant des fonctions d'onde radiales d'électrons exactes de type Dirac-Hartree-Fock-Slater (DHFS) obtenues via le package RADIAL, plutôt que de s'appuyer sur des approximations simplifiées.

• Formalisme : Le taux de capture est calculé comme une somme incohérente de fonctions d'onde partielles de neutrinos. L'amplitude de transition est traitée comme le chevauchement spatial exact des états initiaux et finaux sur l'opérateur d'interaction faible, intégrant la densité de transition nucléaire ($\rho^{TD}_{1,1}$) avec les fonctions d'onde leptoniques (électron et neutrino).
• Traitement non-factorisé : Contrairement à l'approche factorisée standard où la structure nucléaire s'annule dans les rapports, ce travail conserve la pleine dépendance radiale de la densité de transition. Les auteurs analysent le comportement radial des noyaux leptoniques (produits des fonctions d'onde électron et neutrino) pour diverses couches atomiques (K, L, M).
• Contraintes phénoménologiques : Pour traiter la complexité du calcul de la densité de transition nucléaire à partir de principes premiers pour les noyaux de taille moyenne, les auteurs adoptent des densités de transition contraintes phénoménologiquement (spécifiquement des paramétrisations Double-Gaussienne et Double-Gaussienne modifiée) précédemment proposées dans la Réf. [7] pour résoudre l'anomalie du gallium.
• Corrections atomiques : Le calcul inclut les corrections d'échange et de recouvrement ($B_x$) ainsi que les effets de shake-up/shake-off ($S_x$), évalués à l'aide des fonctions d'onde électroniques spécifiques dérivées dans ce travail.

Contributions clés

1. Amplitudes radiales explicites : L'article dérive des amplitudes radiales explicites pour toutes les couches atomiques pertinentes (K, $L_{1-3}$, $M_{1-5}$) en évaluant les contributions dominantes des canaux de neutrinos. Il identifie que les orbitales $s_{1/2}$ et $p_{1/2}$ sont dominées par la fonction de Bessel d'ordre zéro ($j_0$), tandis que les couches supérieures impliquent des fonctions de Bessel d'ordres supérieurs ($j_1, j_2$, etc.) qui sont significativement supprimées près de l'origine.
2. Analyse de la validité de la factorisation : Les auteurs démontrent que les couches dominantes (K, $L_1$, $L_2$, $M_1$, $M_2$) possèdent des amplitudes leptoniques effectivement plates à l'intérieur du rayon nucléaire. Par conséquent, ces couches sondent des moments radiaux similaires de la densité de transition nucléaire. En revanche, les couches subdominantes comme $L_3$ présentent une dépendance radiale différente ($r^2$), rendant leur intensité effective sensible à la forme spécifique de la densité de transition.
3. Réévaluation des rapports : L'étude réévalue les rapports de capture L/K, M/K et M/L en utilisant à la fois l'approximation standard factorisée à l'ordre dominant (LO) et le schéma complet non-factorisé avec les densités de transition phénoménologiques.

Résultats

• Impact de la non-factorisation : L'introduction du schéma complet non-factorisé, même avec des densités de transition variant considérablement en forme (par exemple, s'étendant au-delà du rayon nucléaire), affecte à peine les prédictions pour les rapports dominants L/K et M/K. Les différences entre les schémas factorisés et non-factorisés restent dans les incertitudes théoriques pour ces rapports dominants.
• Couches subdominantes : L'approche non-factorisée permet de prédire les contributions individuelles des couches subdominantes (par exemple, $L_3$, $M_3$). Bien que le profil spécifique de la densité de transition puisse modifier ces prédictions subdominantes d'un facteur allant jusqu'à deux, ces contributions sont hautement supprimées (ordres de grandeur plus petites que les couches dominantes) et restent impossibles à mesurer avec la sensibilité expérimentale actuelle ou future proche.
• Comparaison avec les données :
  • Les prédictions théoriques pour le rapport L/K dans le schéma factorisé montrent une divergence de $3,6\sigma$ avec la moyenne mondiale expérimentale actuelle.
  • Les divergences pour les rapports M/K et M/L sont d'environ $2\sigma$.
  • Cependant, les prédictions théoriques sont en excellent accord avec la dernière mesure expérimentale de CONUS+ [29], qui rapporte un rapport M/K précis de $0,0214(4)$ et un rapport L/K de $0,119(1)$.
  • La tension entre la prédiction théorique et l'ancienne moyenne mondiale expérimentale persiste même en tenant compte des effets de non-factorisation (réduisant la divergence de seulement $3,6\sigma$ à $3,3\sigma$).

Signification et affirmations
L'article conclut que l'approximation standard factorisée à l'ordre dominant (LO) est robuste pour interpréter les mesures actuelles des rapports de capture électronique du $^{71}$Ge. L'analyse confirme que les effets de structure nucléaire, médiés par la densité de transition, n'altèrent pas de manière significative les rapports de capture dominants (L/K, M/K) car les couches atomiques dominantes partagent une dépendance radiale similaire qui annule efficacement les variations de la structure nucléaire dans le rapport.

Les auteurs affirment que la tension de $3,6\sigma$ observée entre leurs prédictions théoriques et les anciennes moyennes mondiales expérimentales est peu probable d'être résolue par des effets de structure nucléaire ou des mécanismes de non-factorisation. Au lieu de cela, l'accord avec les données récentes de CONUS+ suggère que les anciens jeux de données pourraient contenir des problèmes systématiques non pris en compte. L'article souligne que de nouvelles mesures de haute précision sont essentielles pour clarifier le tableau expérimental, car les données actuelles restent complexes et contradictoires. Ce travail fournit une évaluation actualisée de l'état théorique et expérimental de la CE du $^{71}$Ge, validant l'utilisation des modèles factorisés standards pour les couches dominantes tout en offrant un cadre rigoureux pour de futures investigations sur les effets subdominants.