Shell model description of the $N=82$ isotonic chain with a new effective interaction

Cet article présente une étude systématique du modèle en couches de la chaîne isotonique $N=82$ ($Z=51$–77) utilisant une interaction effective nouvellement développée via l'analyse en composantes principales, laquelle reproduit avec succès les propriétés nucléaires expérimentales et fournit des prédictions pour les noyaux riches en protons au-delà de la portée expérimentale actuelle.

L'essentiel

Imaginez le noyau atomique non pas comme une boule solide, mais comme un immeuble d'appartements animé à plusieurs étages où vivent de minuscules particules appelées protons et neutrons. Dans cet immeuble, il existe des « étages » ou des niveaux d'énergie spécifiques où ces particules préfèrent se retrouver. Parfois, un étage est complètement rempli, créant un quartier très stable et heureux. En physique nucléaire, nous appelons ces étages pleins des « nombres magiques ».

Cet article porte sur un quartier spécifique où l'étage des neutrons est complètement rempli (le nombre magique 82). Les scientifiques voulaient comprendre comment les protons se comportent dans les étages situés au-dessus de cette base stable, plus précisément dans une gamme d'éléments allant du Tellure à l'Iridium.

Voici une décomposition de ce qu'ils ont fait et découvert, en utilisant des analogies simples :

1. Le Problème : La « Carte » n'était pas parfaite

Les scientifiques ont essayé de dessiner une carte parfaite de la manière dont les protons interagissent entre eux. Les cartes précédentes (appelées « interactions effectives ») étaient correctes, mais elles comportaient des erreurs. Elles étaient comme un GPS qui vous disait parfois de tourner à gauche alors que vous deviez tourner à droite, ou qui prédisait qu'un bâtiment ferait 3 mètres de haut alors qu'il en faisait en réalité 3,6.

Plus précisément, les anciennes cartes peinaient à prédire :

• L'énergie exacte des certains états excités (comme la hauteur d'un rebond de balle).
• Le « spin » ou l'orientation de noyaux lourds et impairs (comme prédire dans quelle direction un toupie va tomber).
• Le comportement de noyaux très lourds et riches en protons, qui sont difficiles à étudier en laboratoire.

2. La Solution : Une Nouvelle Carte plus Intelligente

Les auteurs ont créé une toute nouvelle carte de haute qualité en utilisant une méthode appelée Analyse en Composantes Principales (ACP).

Voyez cela comme l'accordage d'un instrument de musique massif et complexe possédant 165 cordes différentes (les paramètres d'interaction). Au lieu d'essayer d'accorder chaque corde parfaitement en devinant, ils ont utilisé un algorithme intelligent pour trouver les 30 cordes les plus importantes qui changent réellement le son de la musique. Ils ont ensuite « accordé » ces 30 cordes en écoutant 204 notes réelles (des données provenant de noyaux réels) et en ajustant la carte jusqu'à ce que la musique corresponde parfaitement.

Le résultat ? Une carte incroyablement précise. La différence entre leurs prédictions et les mesures réelles est infime — environ la largeur d'un seul noyau atomique (102 keV).

3. Ce qu'ils ont Découvert

Avec cette nouvelle carte précise, ils ont pu décrire le « quartier » avec un grand détail :

• La « Sub-fermeture » à Z=64 : Ils ont confirmé qu'à un nombre de protons spécifique (64, qui est le Gadolinium), il existe un « sous-étage » spécial qui agit comme un mini-mur. Cela rend le noyau extra stable et plus difficile à exciter, un peu comme un bâtiment avec un plancher en béton armé au milieu. Leur carte a parfaitement décrit cela.
• Prédire l'Invisible : Parce que leur carte est si fiable, ils l'ont utilisée pour prédire les propriétés de noyaux si lourds et instables que les scientifiques n'ont pas encore pu mesurer. Ils ont fait des prédictions spécifiques pour des noyaux tels que le Tantale-155, le Tungstène-156, le Rhénium-157, l'Osmium-158 et l'Iridium-159. Ils ont prédit des choses comme si ces noyaux resteraient ensemble ou s'effondreraient (émission de protons).
• Résoudre des Mystères : Ils ont résolu un puzzle de longue date concernant l'« état fondamental » (la position de repos) de certains noyaux lourds. Les anciennes cartes se trompaient de direction pour le spin de certains de ces noyaux ; la nouvelle carte réussit à chaque fois.

4. Ce qu'il faut retenir

Cet article consiste essentiellement à construire un « livre de règles » meilleur et plus fiable pour la façon dont les protons se comportent dans une région spécifique du monde atomique. En utilisant une approche mathématique plus intelligente pour ajuster les données, ils ont créé un outil qui non seulement explique ce que nous savons déjà, mais prédit aussi avec confiance ce que nous n'avons pas encore vu.

Ils n'ont pas seulement corrigé les chiffres ; ils ont fourni une image claire de la structure sous-jacente de ces atomes, montrant exactement sur quels « étages » les protons vivent et comment ils interagissent avec leurs voisins. Ce nouveau livre de règles est maintenant disponible pour que d'autres scientifiques puissent l'utiliser lors de futures études sur ces éléments lourds.

Résumé technique

Résumé technique : Description par le modèle en couches de la chaîne isotonique N = 82 avec une nouvelle interaction effective

Énoncé du problème
La structure nucléaire de la chaîne isotonique $N=82$, caractérisée par des protons de valence dans la couche $g_{7/2}d_{5/2}s_{1/2}d_{3/2}h_{11/2}$ au-dessus du cœur $^{132}\text{Sn}$, a fait l'objet d'études approfondies. Bien que les interactions effectives précédentes (par exemple, JJ56PNA, KH5082, CW5082) aient fourni des informations, elles présentent des divergences systématiques par rapport aux données expérimentales. Notamment, ces interactions ont tendance à surestimer les énergies d'excitation des états $3^-$ et ne parviennent pas à prédire correctement les spins et les parités des états fondamentaux pour les isotones $A$ impairs plus lourds ($Z > 64$). De plus, des prédictions théoriques fiables sont nécessaires pour les noyaux riches en protons proches et au-delà de la ligne de drip de protons (par exemple, $^{155}\text{Ta}$, $^{156}\text{W}$), où les données expérimentales sont rares ou inexistantes.

Méthodologie
Les auteurs emploient le modèle en couches à interaction complète (full-configuration-interaction) via le code BIGSTICK dans un espace de modèle comprenant tous les orbitales protoniques entre $Z=50$ et $Z=82$. L'innovation majeure réside dans le développement d'une nouvelle interaction effective de haute qualité, optimisée via l'Analyse en Composantes Principales (ACP).

1. Espace des paramètres : Le Hamiltonien comprend 5 énergies de particules uniques (SPE) et 160 éléments de matrice à deux corps (TBME) proton-proton, soit un total de 165 paramètres d'entrée.
2. Stratégie d'optimisation : Au lieu de l'ajustement direct de tous les paramètres, les auteurs utilisent l'ACP pour tronquer l'espace des paramètres. Ils effectuent un ajustement par moindres carrés sur 204 points de données expérimentaux, qui incluent :
   • 21 énergies de liaison.
   • 101 énergies d'excitation dans les noyaux pair-pair (se concentrant sur les états yrast et des états non-yrast spécifiques comme $2^+_2, 3^-$, etc.).
   • 82 énergies d'excitation dans les noyaux $A$ impairs.
3. Implémentation de l'ACP : La matrice d'erreur de l'ajustement est diagonalisée. L'espace des paramètres est tronqué aux 30 composantes principales (combinaisons linéaires des paramètres originaux) associées aux plus petites valeurs propres (plus petites incertitudes). Cette approche limite le surapprentissage tout en conservant les paramètres les plus sensibles.
4. Contraintes : Des facteurs d'échelle de masse standard sont appliqués aux TBME. Les effets Coulombiens sont incorporés implicitement via la procédure d'ajustement. Des charges effectives ($e_\pi = 1.5e$) et des facteurs de trempe ($g_s = 5.586 \times 0.7$) sont utilisés pour les propriétés électromagnétiques.

Contributions clés et résultats

• Optimisation de l'interaction : La nouvelle interaction atteint un écart quadratique moyen (RMSD) de 102 keV sur les 204 points de données. Cela représente une amélioration significative par rapport aux interactions précédentes, qui présentaient souvent des écarts de l'ordre du MeV pour les énergies de liaison.
• Énergies de liaison et de séparation : Le modèle reproduit les énergies de liaison de $^{134}\text{Te}$ à $^{151}\text{Tm}$ avec un RMSD de 59 keV. Il prédit avec succès les excès de masse atomique et les énergies de séparation pour les noyaux au-delà de la portée expérimentale actuelle, incluant $^{156}\text{W}$ jusqu'à $^{159}\text{Ir}$. Les calculs suggèrent que $^{153}\text{Lu}$ est au-delà de la ligne de drip de protons, tandis que $^{154}\text{Hf}$ reste lié. Pour $^{155}\text{Ta}$, le modèle prédit une émission de proton unique, et pour $^{156}\text{W}$, il prédit une énergie de séparation de un proton positive mais une énergie de séparation de deux protons négative, ce qui est cohérent avec la suppression observée de la désintégration à deux protons.
• Spectres de basse énergie :
  • Noyaux pair-pair : Le modèle reproduit la tendance systématique des énergies d'excitation $2^+_1$, identifiant correctement le pic à $^{146}\text{Gd}$ ($Z=64$) et le maximum local à $^{140}\text{Ce}$, confirmant la fermeture de la sous-couche $Z=64$. Il résout les divergences précédentes en prédisant avec précision les énergies d'excitation $3^-$, particulièrement le minimum observé dans $^{146}\text{Gd}$.
  • Noyaux $A$ impairs : L'interaction reproduit correctement les spins et parités des états fondamentaux sur toute la chaîne, y compris les états difficiles de spin $1/2^+$ pour $^{147}\text{Tb}$ et de $3/2^+$ pour $^{155}\text{Ta}$ et $^{157}\text{Re}$, qui étaient auparavant difficiles à décrire.
• Propriétés électromagnétiques :
  • Valeurs B(E2) : Le modèle concorde généralement bien avec les valeurs expérimentales de $B(E2)$. Cependant, des écarts sont notés pour certaines transitions dans $^{138}\text{Ba}$ et $^{139}\text{La}$, suggérant un possible renversement dans l'ordonnancement des états $9/2^+$ calculés dans $^{139}\text{La}$.
  • Moments magnétiques : Les moments magnétiques dipolaires ($\mu$) calculés concordent bien avec l'expérience (écarts $< 0.5 \mu_N$).
  • Forces B(E3) : Le modèle sous-estime les forces $B(E3; 3^- \to 0^+)$ d'un facteur d'environ deux. Les auteurs attribuent cela à des contributions manquantes des excitations de neutrons à travers la lacune de la couche $N=82$, notant qu'augmenter la charge effective à 2.0 permet d'obtenir concordance avec les résultats.
• Analyse structurelle :
  • Écarts de sous-couches : Les énergies de particules uniques de champ moyen (MSPE) révèlent un écart de sous-couche de 1,788 MeV à $Z=64$ (entre $1d_{5/2}$ et $0h_{11/2}$) et de 1,062 MeV à $Z=58$ (entre $0g_{7/2}$ et $1d_{5/2}$).
  • Seniorité : L'étude analyse la structure de seniorité des états yrast. Les états $6^+_1$ dans les isotones légers sont dominés par des configurations $(0g_{7/2})^2$ ou $0g_{7/2}1d_{5/2}$, tandis que ceux des isotones plus lourds ($Z > 64$) correspondent à $(0h_{11/2})^2$.
  • Isomérie : Le modèle prédit plusieurs états excités de très basse énergie dans les noyaux $A$ impairs (par exemple, un état $1/2^+$ dans $^{155}\text{Ta}$ à 4 keV), suggérant une potentielle isomérie.

Signification et revendications
L'article affirme que la nouvelle interaction optimisée fournit une description systématique et fiable de la chaîne isotonique $N=82$, résolvant les divergences de longue date concernant les énergies d'excitation et les spins fondamentaux qui ont entravé les modèles précédents. En utilisant l'ACP pour optimiser l'interaction contre un large ensemble de données, les auteurs démontrent qu'une description à une seule couche majeure avec une interaction optimisée est suffisante pour capturer des caractéristiques structurelles complexes, telles que la fermeture de la sous-couche $Z=64$ et l'évolution des motifs de seniorité, sans avoir besoin de recourir à de larges excitations à travers les couches majeures.

Les auteurs soulignent que cette interaction sert de fondation solide pour les études futures des noyaux proches de $N=82$. Ils citent explicitement son application (et son utilisation continue) dans :

1. Les études de l'évolution de la forme nucléaire et des transitions de phase quantique dans les noyaux des terres rares.
2. Les investigations systématiques des isotones $N=80$ pair-pair.
3. La première mesure spectroscopique de $^{150}\text{Yb}$.

L'étude conclut en fournissant l'interaction optimisée comme matériel supplémentaire pour faciliter la recherche ultérieure dans cette région de masse.