⚛️ théorie nucléaire

Evaluation of $2\nu\beta\beta$ decay nuclear matrix elements using Lanczos techniques preserving Gamow-Teller strengths from both parent and daughter nuclei

Cet article présente et valide deux approches efficaces basées sur la méthode de Lanczos, DLSF et BLSF, pour le calcul des éléments de matrice de la désintégration double bêta à deux neutrinos qui préservent les intensités de Gamow-Teller dans les noyaux pères et filles, en les validant avec succès sur le $^{48}$ Ca et en les appliquant à d'autres candidats tels que le $^{58}$ Ni, le $^{76}$ Ge et le $^{130}$ Te.

Auteurs originaux : Deepak Patel, Anil Kumar, Noritaka Shimizu, Yutaka Utsuno

Publié 2026-11-03

📖 5 min de lecture🧠 Analyse approfondie

CC BY 4.0

Auteurs originaux : Deepak Patel, Anil Kumar, Noritaka Shimizu, Yutaka Utsuno

Article original sous licence CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ Ceci est une explication générée par l'IA de l'article ci-dessous. Elle n'a pas été rédigée ni approuvée par les auteurs. Pour une précision technique, consultez l'article original. Lire la clause de non-responsabilité complète

Au plus profond du cœur de l'atome, les protons et les neutrons sont liés ensemble dans un équilibre délicat. Parfois, cet équilibre bascule, et un noyau cherche la stabilité en se transformant lui-même. Dans un type de transformation rare et spécifique appelé désintégration double bêta, deux neutrons à l'intérieur d'un noyau se transforment simultanément en deux protons, éjectant deux électrons et deux particules invisibles appelées neutrinos. Ce processus est si lent que pour certains éléments, il ne se produit qu'une fois tous les quelques billions d'années. Bien que les scientifiques aient longtemps observé cette version de la désintégration à deux neutrinos, ils traquent une variante fantomatique et encore plus rare où aucun neutrino n'est émis du tout. Trouver cette version neutrinolessse manquante prouverait que les neutrinos sont leurs propres antiparticules, une découverte qui pourrait expliquer pourquoi l'univers est fait de matière plutôt que d'un vide de matière antimatière. Pour trouver cet événement rare, les physiciens doivent d'abord maîtriser les règles de la version commune à deux neutrinos, en calculant exactement comment le noyau se comporte pendant la transformation. Cela nécessite de cartographier les états d'énergie invisibles de l'atome avec une précision extrême, une tâche qui a longtemps été un goulot d'étranglement computationnel pour les scientifiques tentant de comprendre les blocs fondamentaux de notre monde.

Une équipe de chercheurs a maintenant développé une nouvelle façon hautement efficace d'effectuer ces calculs, contournant les obstacles computationnels massifs qui limitaient auparavant leur précision. Au lieu d'essayer de lister chaque état d'énergie possible par lequel le noyau pourrait passer pendant la désintégration — une tâche qui devient impossible à mesure que l'atome s'alourdit — ils ont utilisé un raccourci mathématique pour capturer le comportement essentiel de l'ensemble du système. En se concentant sur les états « de passage » spécifiques qui initient la transformation, ils ont pu reconstruire l'image complète de la désintégration sans avoir besoin de calculer individuellement des milliers d'étapes intermédiaires. Cette approche leur a permis de préserver la force totale des forces nucléaires impliquées à la fois depuis le point de départ et le point d'arrivée de la désintégration, garantissant qu'aucune information cruciale ne soit perdue dans la simplification. Le résultat est une méthode qui est non seulement plus rapide mais aussi plus fiable pour les atomes lourds où la rare désintégration neutrinolessse pourrait se cacher.

Les chercheurs ont testé leurs nouvelles techniques sur plusieurs atomes candidats, en commençant par le calcium-48, un élément bien étudié souvent utilisé comme référence. En élargissant la « pièce » mathématique dans laquelle ils calculaient le comportement de l'atome pour inclure des mouvements plus complexes de ses particules, ils ont constaté que la force prédite de la désintégration diminuait légèrement, rapprochant la valeur théorique de ce que les expériences mesurent réellement. Cet ajustement suggère que regarder l'atome avec un objectif plus large fournit une description plus précise de sa structure interne. Ils ont ensuite appliqué leurs méthodes au nickel-58, au germanium-76 et au tellure-130, couvrant une gamme de poids atomiques. Pour le nickel-58, ce travail représente la première fois qu'un calcul détaillé du modèle en couches est effectué pour ce type spécifique de désintégration, révélant une valeur théorique bien plus petite que les limites expérimentales actuelles. Cette divergence indique que, bien que la théorie soit robuste, les expériences doivent devenir nettement plus sensibles pour détecter le processus, s'il a lieu.

Pour s'assurer que leurs nouveaux raccourcis n'introduisaient pas d'erreurs, l'équipe a effectué une vérification rigoureuse. Ils ont comparé la somme de toutes les étapes intermédiaires calculées par rapport à une mesure directe de la force totale entre les états de départ et d'arrivée, une valeur qui peut être déterminée indépendamment. Les deux nombres correspondaient parfaitement pour tous les atomes étudiés. Cet accord confirme que leur méthode réussit à récupérer l'image complète de la transformation nucléaire, capturant chaque contribution pertinente sans avoir besoin de calculer explicitement chaque état intermédiaire. Les chercheurs ont constaté que leurs nouvelles approches convergeaient vers la bonne réponse avec beaucoup moins d'étapes de calcul que les méthodes traditionnelles, rendant possible l'étude de noyaux plus lourds et plus complexes qui étaient auparavant hors de portée.

Les implications de ce travail s'étendent au-delà de la simple précision des chiffres pour les désintégrations connues. En fournissant une façon fiable et efficace de calculer les éléments de matrice nucléaire — les chiffres clés qui relient la théorie à l'expérience — cette étude renforce les fondations de la recherche de la désintégration neutrinolessse. Si la version rare et neutrinolessse de ce processus existe, son taux dépend entièrement de ces calculs nucléaires. En démontrant que leurs méthodes peuvent modéliser avec précision la désintégration commune tout en préservant la force totale des forces nucléaires, les chercheurs ont fourni un outil de confiance pour la communauté mondiale. Leur travail suggère que la voie à suivre consiste à étendre les modèles pour inclure des interactions de particules plus complexes, comme on l'a vu dans leurs résultats améliorés pour le calcium-48, et à pousser les détecteurs expérimentaux vers de plus grandes sensibilités, comme le cas du nickel-58 l'exige clairement. La traque des secrets les plus profonds de l'univers continue, mais la carte du paysage atomique est désormais dessinée avec plus de clarté et de confiance.

Résumé technique : Évaluation des éléments de matrice de la désintégration $2\nu\beta\beta$ par les techniques de Lanczos

Énoncé du problème
Le calcul des éléments de matrice nucléaires (NME) pour la désintégration double bêta à deux neutrinos ( $2\nu\beta\beta$ ) et la capture électronique double à deux neutrinos ( $2\nu$ ECEC) est essentiel pour interpréter les données expérimentales et fournir des points de référence pour les modèles de structure nucléaire. Ces NME nécessitent une sommation sur l'ensemble complet des états propres intermédiaires $1^+$ du noyau fils. Cependant, dans les noyaux de masse moyenne et lourde, la dimension de la matrice hamiltonienne croît rapidement, rendant le calcul explicite de tous les états $1^+$ intermédiaire trop coûteux en termes de calcul. Les calculs conventionnels en modèle de couches reposent souvent sur des espaces de modèles tronqués ou un nombre limité d'états intermédiaires, ce qui peut omettre une fraction significative de la force Gamow-Teller (GT) et compromettre la convergence et la fiabilité des NME.

Méthodologie
Les auteurs proposent deux approches efficaces de modèle de couches basées sur les techniques de Lanczos pour évaluer les NME de $2\nu\beta\beta$ sans calculer explicitement l'ensemble complet des états propres $1^+$ intermédiaires :

Fonction de force de Lanczos à double côté (DLSF) : Cette méthode construit des vecteurs de départ en appliquant des opérateurs GT aux états fondamentaux $0^+$ initial et final. Elle effectue des itérations de Lanczos des deux côtés, générant deux ensembles d'états $1^+$ approximatifs. Ces ensembles sont combinés et orthogonalisés à l'aide de l'orthogonalisation de Löwdin, suivis d'une rediagonalisation de l'Hamiltonien dans le sous-espace résultant pour obtenir les vecteurs de Ritz.
Fonction de force de Lanczos par blocs (BLSF) : Cette méthode utilise une procédure de Lanczos par blocs où les vecteurs initiaux (dérivés du noyau parent et du noyau fils) forment un bloc. L'Hamiltonien agit simultanément sur ce bloc, générant un sous-espace de Krylov par blocs. L'Hamiltonien projeté prend une forme tridiagonale symétrique par blocs, qui est ensuite diagonalisée pour produire les vecteurs de Ritz.

Les deux méthodes sont conçues pour préserver simultanément les forces GT totales associées aux noyaux initial et final. Elles sont implémentées dans le code KSHELL. Pour valider les méthodes, les auteurs comparent les sommes des NME sans dénominateur d'énergie contre les éléments de matrice double-Gamow-Teller (DGT) exacts calculés directement entre les états fondamentaux initial et final en utilisant les relations de clôture.

Résultats clés
Les méthodes ont été testées et appliquées à plusieurs noyaux dans différentes régions de masse :

$^{48}$ Ca (couche fp) : En utilisant l'interaction GXPF1A, la DLSF et la BLSF ont toutes deux convergé vers un NME de 0,0511 en 30 itérations, correspondant au résultat de la diagonalisation exacte utilisant 2 500 états intermédiaires. Les méthodes ont préservé avec succès les forces GT pour les transitions $^{48}$ Ca $\to$ $^{48}$ Sc et $^{48}$ Ti $\to$ $^{48}$ Sc, une caractéristique qui n'était pas atteinte simultanément par la précédente méthode de Lanczos à un seul côté (SLSF).
$^{48}$ Ca (couche sdpf) : L'extension de l'espace de modèle pour inclure les excitations inter-couches jusqu'à $2\hbar\omega$ en utilisant l'interaction SDPF-MU a produit un NME convergé de 0,0463. Cette valeur est plus proche de l'estimation expérimentale (dérivée des données de demi-vie) que le résultat de la couche fp, suggérant que l'élargissement de l'espace de modèle améliore la description de la désintégration.
$^{58}$ Ni (couche fp) : Les auteurs présentent le premier calcul de modèle de couches du NME $2\nu\beta^+EC/2\nu$ ECEC pour $^{58}$ Ni. Le NME convergé est d'environ 0,1250. Comparé à la limite inférieure expérimentale actuelle sur la demi-vie, le NME théorique est nettement plus petit, indiquant que les expériences actuelles manquent de sensibilité pour contraindre efficacement le NME.
$^{76}$ Ge et $^{130}$ Te : Pour $^{76}$ Ge (interaction JUN45) et $^{130}$ Te (interaction SN100PN), les méthodes ont atteint une convergence complète (NME de 0,1656 et 0,0795, respectivement) avec relativement peu d'itérations (75 et 60). Dans le cas de $^{130}$ Te, la diagonalisation exacte jusqu'à 1 000 états a montré une convergence lente, soulignant la nécessité des approches de Lanczos pour les espaces de modèles plus larges.
Validation : Pour tous les noyaux testés, les sommes obtenues via DLSF et BLSF sans dénominateur d'énergie concordaient précisément avec les éléments de matrice DGT exacts calculés directement. Cela confirme que les méthodes récupèrent correctement l'ensemble des contributions des états intermédiaires.

Signification et revendications
L'article affirme que les méthodes DLSF et BLSF fournissent un cadre robuste et efficace sur le plan computationnel pour calculer les NME de $2\nu\beta\beta$ dans de grands espaces de modèles de couches. En préservant simultanément les forces GT des noyaux parent et fils, ces méthodes offrent un traitement plus symétrique et cohérent des contributions des états intermédiaires par rapport aux approches précédentes à un seul côté. La reproduction réussie des éléments de matrice DGT exacts sert de validation indépendante que les méthodes capturent toute la force GT accessible dans les espaces de modèles adoptés. Les résultats pour $^{48}$ Ca suggèrent que l'élargissement de l'espace de modèle pour inclure les excitations inter-couches produit des NME plus proches des valeurs expérimentales, soutenant ainsi l'utilisation de tels espaces élargis pour prédire les NME de double bêta sans neutrino ( $0\nu\beta\beta$ ). L'étude souligne également la nécessité de mesures expérimentales plus sensibles de la demi-vie $2\nu\beta^+EC$ de $^{58}$ Ni pour fournir des contraintes significatives sur les modèles théoriques.

Articles similaires