Theoretical analysis and predictions for the double electron capture of $^{124}$Xe

Cet article présente une analyse théorique complète de la double capture électronique à deux neutrinos dans le $^{124}$Xe en améliorant les calculs de structure nucléaire et atomique, ce qui produit des éléments de matrice nucléaire affinés, prédit des fractions de capture spécifiques pour divers canaux de désintégration (notamment 74 % pour le canal KK et 24 % pour les canaux cumulatifs KL$_1$-KO$_1$) et fournit des énergies de relaxation atomique mises à jour pour la modélisation du bruit de fond dans les expériences de xénon liquide.

L'essentiel

Imaginez l'atome comme une petite ville en pleine effervescence. À l'intérieur de cette ville, le noyau est l'hôtel de ville, rempli de protons et de neutrons. Habituellement, ces citoyens sont très stables, mais parfois, ils décident de se réorganiser pour être plus à l'aise.

Ce document porte sur un événement de « réorganisation » très rare et spécifique se produisant dans une ville appelée Xénon-124. Lors de cet événement, l'hôtel de ville décide de capturer deux de ses propres résidents (des électrons) provenant des quartiers extérieurs pour les attirer à l'intérieur du noyau. Quand cela arrive, la ville se transforme en une nouvelle ville appelée Tellure-124, et elle recrache deux messagers minuscules et invisibles appelés neutrinos.

Les scientifiques appellent cela la Capture Électronique Double (plus précisément la version à deux neutrinos, ou $2\nu$ECEC). C'est comme un double plongeon dans une piscine, mais au lieu de l'eau, il s'agit de particules subatomiques.

Voici ce que les chercheurs ont fait, expliqué simplement :

1. Construire un meilleur plan (La Théorie)

Par le passé, les scientifiques essayaient de prédire la fréquence à laquelle ce « double plongeon » se produit, mais leurs plans étaient un peu rudimentaires. Ils manquaient certains détails sur la façon dont les électrons se déplacent et dont le noyau réagit.

Les auteurs de ce document ont décidé de construire un plan beaucoup plus précis.

• L'analogie de la « Développement de Taylor » : Imaginez essayer de décrire le trajet d'une voiture. Une description simple pourrait simplement dire « elle avance ». Une meilleure description ajoute « elle accélère ». La meilleure description ajoute « elle accélère, puis tourne légèrement, puis ralentit ». Les auteurs ont utilisé un outil mathématique appelé « développement de Taylor » pour ajouter ces couches de détails supplémentaires (jusqu'à la quatrième puissance de l'énergie). Cela leur a permis de voir les « virages et ralentissements » du processus de désintégration que les modèles précédents avaient manqués.
• Les « Nouveaux Ratios » : Parce qu'ils ont ajouté ces détails supplémentaires, ils ont découvert de nouvelles façons de comparer différentes parties du processus (appelées ratios $\xi$). Considérez cela comme de nouveaux points de contrôle sur un circuit de course qui peuvent aider les scientifiques à mieux comprendre la course plus tard.

2. Observer le voisinage (La partie Atomique)

Pour calculer la probabilité que cet événement se produise, vous devez savoir exactement où vivent les électrons.

• La métaphore du « Blocage de Pauli » : Imaginez un ascenseur bondé. Si l'ascenseur est plein, vous ne pouvez pas simplement y pousser quelqu'un d'autre ; vous devez attendre que quelqu'un en sorte. Dans le noyau, les places les plus « intérieures » sont comme un ascenseur plein. Les auteurs ont réalisé que les électrons capturés ne peuvent pas aller n'importe où ; ils sont bloqués par les autres électrons déjà présents. Ils ont pris en compte cette règle d'« encombrement », ce qui modifie le calcul.
• Élargir la recherche : Les études précédentes ne regardaient que les deux quartiers les plus proches du noyau (appelés couches K et L1). Les auteurs ont dit : « Regardons tous les quartiers, même ceux qui sont plus éloignés (jusqu'à la couche O). » Ils ont découvert que, bien que les quartiers extérieurs soient moins susceptibles d'être capturés, ils contribuent tout de même à l'événement total.

3. Simuler l'Hôtel de Ville (La partie Nucléaire)

Le noyau est la partie la plus difficile à simuler car c'est une foule chaotique de particules. Les auteurs ont utilisé deux « moteurs de simulation » différents pour prédire comment le noyau se comporte :

• Moteur A (ISM) : C'est comme une simulation détaillée, pièce par pièce. Ils ont fait tourner la simulation avec différentes règles (appelées « Hamiltoniens ») pour voir si les résultats tena-ient la route. Ils ont constaté que lorsqu'ils incluaient toutes les étapes « intermédiaires » possibles que le noyau emprunte, la « force » prédite de l'événement était plus faible que ce que suggéraient les modèles plus anciens et plus simples.
• Moteur B (pn-QRPA) : Il s'agit d'un type différent de simulation. Ils ont ajusté les paramètres de ce moteur jusqu'à ce qu'il corresponde aux données réelles que nous possédons déjà. Ils ont trouvé que leur nouveau calcul, plus soigneux, donnait une valeur de « force » beaucoup plus petite que les tentatives précédentes utilisant ce moteur.

4. Les Résultats : Qu'ont-ils trouvé ?

En combinant leur meilleur plan, leur carte détaillée du voisinage et leurs deux moteurs de simulation, ils ont fait plusieurs prédictions :

• L'Événement Principal (Canal KK) : Ils prédisent qu'environ 74 % du temps, les deux électrons seront capturés depuis le quartier le plus proche (la couche K). C'est légèrement différent des 72,4 % utilisés dans les expériences précédentes, un petit mais important ajustement.
• Les « Prochains Meilleurs » Événements : Ils prédisent qu'environ 19 % du temps, un électron provient du quartier le plus proche et l'autre du quartier suivant (KL1).
• La Prédiction « Cumulative » : Si l'on additionne tous les événements légèrement moins fréquents (de KL1 jusqu'à KO1), ils représentent environ 24 % de l'événement total. Cela représente environ un tiers de l'événement principal.
• L'Énergie de « Relaxation » : Lorsque les électrons sont capturés, la nouvelle ville (le Tellure) est excitée et doit se calmer. Elle le fait en libérant de l'énergie (comme des rayons X). Les auteurs ont calculé exactement quelle quantité d'énergie est libérée pour chaque type de capture. C'est comme donner aux scientifiques une « empreinte digitale » d'énergie spécifique à rechercher dans leurs détecteurs.

Pourquoi est-ce important ?

Ce document ne prétend pas guérir des maladies ou alimenter des villes. Au lieu de cela, il agit comme une carte raffinée pour les explorateurs.

De grandes expériences utilisant du xénon liquide (comme celles qui recherchent la Matière Noire) surveillent constamment ce type d'événement de « double plongeon ». Cependant, cet événement peut ressembler à un « bruit de fond » qui confond les données. En fournissant une carte plus précise de la fréquence à laquelle cela se produit, de l'énergie qu'il libère et des « quartiers » dont proviennent les électrons, les auteurs aident les expérimentateurs à distinguer un signal réel du bruit de fond.

En bref, ils ont pris une photo floue et à basse résolution d'un événement atomique rare et l'ont transformée en un modèle 3D haute définition, aidant les scientifiques à savoir exactement quoi chercher dans leurs détecteurs.

Résumé technique

Résumé technique : Analyse théorique et prédictions pour la double capture électronique du $^{124}$Xe

Énoncé du problème
Le processus de double capture électronique (ECEC), en particulier le mode deux neutrinos ($2\nu$ECEC), est une transition nucléaire rare d'un intérêt majeur pour sonder la structure nucléaire et les propriétés des neutrinos. Bien que la désintégration $2\nu\beta^-\beta^-$ ait été largement étudiée, les modes ECEC ont reçu moins d'attention en raison de probabilités de transition plus faibles. Cependant, le $^{124}$Xe est un cas unique car il possède la plus grande valeur Q parmi les isotopes capables de subir des transitions ECEC, ce qui en fait une cible primaire pour les expériences à grande échelle utilisant le xénon liquide (par exemple, celles recherchant la matière noire et la désintégration double bêta sans neutrino).

Les calculs théoriques précédents pour la demi-vie du $2\nu$ECEC du $^{124}$Xe souffraient de limitations tant dans la modélisation nucléaire qu'atomique. Plus précisément, les travaux antérieurs reposaient souvent sur l'approximation de fermeture pour les éléments de matrice nucléaire (NME), ne considéraient que les captures des couches K et L$_1$, et utilisaient des traitements moins rigoureux de l'écrantage atomique et du blocage de Pauli. Avec les mesures expérimentales récentes de la demi-vie du $^{124}$Xe devenant disponibles, une description théorique plus rigoureuse et complète est nécessaire pour interpréter ces données et prédire les canaux de désintégration non observés.

Méthodologie
Les auteurs fournissent un cadre théorique complet améliorant à la fois le formalisme et les calculs spécifiques pour le $^{124}$Xe :

1. Formalisme amélioré : Le taux de désintégration est dérivé en utilisant une méthode de développement de Taylor, conservant des termes jusqu'à la quatrième puissance des énergies leptoniques. Cette approche introduit des corrections d'ordre supérieur et définit de nouveaux ratios de NME, $\xi_{31}^{2\nu\text{ECEC}}$ et $\xi_{51}^{2\nu\text{ECEC}}$, qui dépendent des dénominateurs d'énergie des états intermédiaires $1^+$.
2. Calculs de structure nucléaire :
   • Modèle de couches interagissant (ISM) : Les calculs ont été effectués dans l'espace de modèles $jj55$ (orbitales $0g_{7/2}, 1d_{5/2}, 1d_{3/2}, 2s_{1/2}, 0h_{11/2}$). Contrairement aux études précédentes qui supposaient une dominance d'état unique (SSD), ce travail effectue une somme complète sur les états intermédiaires $1^+$ en utilisant une approche de fonction de force. Deux hamiltoniens effectifs (GCN5082 et SVD) ont été employés avec des facteurs de trempe ($q_H$) calibrés pour reproduire les données $2\nu\beta^-\beta^-$ du $^{136}$Xe. Des troncatures permettant jusqu'à quatre excitations de nucléons à partir d'orbitales inférieures ont été incluses.
   • Approximation de l'état fondamental aléatoire de quasi-particules proton-neutrons (pn-QRPA) : Les calculs ont utilisé la restauration de l'isospin. La force d'interaction particule-particule ($g_{pp}^{T=0}$) a été ajustée pour reproduire la demi-vie expérimentale du $^{124}$Xe. L'étude a utilisé des potentiels nucléon-nucléon réalistes (Argonne V18 et CD-Bonn) et a pris en compte le recouvrement du vide BCS.
3. Structure atomique et facteurs d'espace de phase (PSF) :
   • Le cadre auto-cohérent Dirac-Hartree-Fock-Slater (DHFS) a été utilisé pour générer les fonctions d'onde électroniques.
   • Les améliorations incluent un traitement rigoureux de l'écrantage atomique, des corrections de surface nucléaire diffuse, une densité de charge nucléaire réaliste et des effets d'échange-corrélation électronique.
   • Crucialement, le calcul prend en compte le blocage de Pauli des états nucléaires les plus internes et étend la somme de capture à tous les électrons d'onde s disponibles (jusqu'à la sous-couche O$_1$), plutôt que de se limiter aux couches K et L$_1$ comme dans les études antérieures.
   • Les énergies de relaxation atomique pour l'atome final de $^{124}$Te ont été calculées pour aider à la modélisation du bruit de fond.

Résultats clés

• Éléments de matrice nucléaire (NME) :
  • Les calculs ISM produisent des NME qui sont robustes face aux changements de troncature et d'hamiltonien effectif. La somme complète sur les états intermédiaires (sans SSD) a réduit le NME pour l'hamiltonien GCN5082 d'environ 25 % par rapport aux approches tronquées.
  • Les résultats pn-QRPA, lorsqu'ils sont contraints par la demi-vie expérimentale, produisent des valeurs de NME significativement plus basses que les calculs pn-QRPA précédents (qui reposaient souvent sur différentes stratégies d'ajustement pour $g_{pp}$). Les nouveaux NME pour le $^{124}$Xe sont trouvés comparables à ceux de la désintégration $2\nu\beta^-\beta^-$ du $^{128,130}$Te, reflétant des caractéristiques d'appariement similaires.
• Prédictions de demi-vie :
  • En utilisant l'ISM, la demi-vie totale calculée pour le $^{124}$Xe concorde avec la valeur expérimentale de $(1,1 \pm 0,2_{\text{stat}} \pm 0,1_{\text{sys}}) \times 10^{22}$ ans dans un facteur deux.
  • L'inclusion des termes de développement de Taylor et des PSF raffinés a contribué à rapprocher les prédictions théoriques des données expérimentales.
• Fractions de capture et canaux :
  • La fraction de capture prédite pour le canal KK dominant est de 74,1 %, différant des 72,4 % utilisés dans les analyses expérimentales précédentes. Cette différence provient de la considération exclusive des électrons d'onde s dans le modèle théorique par rapport à l'inclusion de tous les électrons dans le modèle de signal expérimental.
  • Le canal suivant le plus probable est KL$_1$ avec une fraction d'environ 19 %.
  • La fraction de capture cumulative pour les canaux KL$_1$ à KO$_1$ est prédite à environ 24 %.
• Énergies de relaxation atomique :
  • Des énergies de relaxation précises ont été calculées pour diverses configurations de trous (par exemple, KK à 64,62 keV, KL$_1$ à 37,05 keV). Ces valeurs sont destinées à être utilisées dans la modélisation du bruit de fond pour les expériences de xénon liquide.

Signification et affirmations
L'article affirme fournir la première description théorique complète du $2\nu$ECEC dans le $^{124}$Xe qui améliore simultanément le traitement de la structure nucléaire (via la somme complète et la restauration de l'isospin) et la physique atomique (via le DHFS, le blocage de Pauli et l'extension de la capture orbitale).

Les auteurs affirment que leurs résultats ISM sont cohérents avec les données expérimentales, validant la fiabilité des NME calculés et des nouvelles fractions de capture prédites. Une contribution primaire est la prédiction que les canaux cumulés KL$_1$-KO$_1$ (représentant ~24 % des désintégrations) devraient être expérimentalement observables dans l'intervalle d'énergie 31,93–37,05 keV. De plus, ce travail suggère que le ratio des demi-vies partielles (par exemple, $T_{KK}^{1/2}/T_{KL_1}^{1/2}$) pourrait servir de contrainte expérimentale sur le paramètre $\xi_{31}^{2\nu\text{ECEC}}$, offrant une nouvelle voie pour sonder l'interaction entre les états intermédiaires de faible énergie et de plus haute énergie. L'étude souligne que les estimations théoriques précédentes de demi-vie étaient souvent incompatibles avec les PSF les plus précis en raison de l'utilisation d'approximations de fermeture et de sommations d'orbitales incomplètes.