Ab initio calculations of nuclear charge radii across and beyond ${}^{132}$Sn: Putting chiral EFT nuclear interactions to the test

Cette étude utilise des calculs *ab initio* de type Bogoliubov coupled cluster pour tester les interactions nucléaires de la théorie effective de champ chirale sur les rayons de charge de la chaîne de l'étain, révélant que les modèles actuels échouent à reproduire simultanément les tendances observées et soulignant la nécessité de mesures expérimentales supplémentaires et d'améliorations théoriques pour mieux contraindre ces interactions au-delà de ${}^{132}$Sn.

L'essentiel

🧪 Le Grand Défi : Mesurer le "Gros" des Atomes

Imaginez que vous essayez de comprendre comment sont construits les Lego de l'univers, les noyaux atomiques. Plus précisément, les chercheurs s'intéressent à une famille d'atomes appelée l'étain (Sn).

Ces atomes d'étain existent sous différentes formes (des isotopes), un peu comme des frères et sœurs qui ont le même visage (le même nombre de protons) mais des poids différents (un nombre variable de neutrons).

Le but de cette étude ? Mesurer la taille de ces noyaux (leur "rayon de charge") et voir comment cette taille évolue quand on ajoute des neutrons. C'est un peu comme essayer de deviner si une famille grandit de manière régulière, ou si elle fait des bonds soudains.

🛠️ L'Outil : La "Recette" de l'Univers

Pour prédire la taille de ces noyaux sans les mesurer en laboratoire, les scientifiques utilisent des supercalculateurs et une "recette" théorique appelée Théorie du Champ Effectif Chiral (χEFT).

Imaginez cette théorie comme une recette de cuisine pour faire des gâteaux (les noyaux).

• Les ingrédients sont les forces qui lient les protons et les neutrons entre eux.
• Les chercheurs ont testé trois versions différentes de cette recette (trois Hamiltoniens) pour voir laquelle donne le meilleur gâteau.

1. La recette classique (1.8/2.0 EM) : Elle est bonne pour le goût (l'énergie), mais le gâteau est trop petit (elle sous-estime la taille des noyaux).
2. La recette améliorée (ΔNNLOGO) : Elle a été ajustée pour mieux coller à la réalité, mais elle a encore des défauts.
3. La recette "magique" (1.8/2.0 EM7.5) : C'est la star du moment. Elle a été finement réglée pour donner la bonne taille aux gâteaux les plus simples (comme l'oxygène).

📉 Le Test : La Courbe en Parole et le "Coup de Couteau"

Les chercheurs ont fait tourner leurs calculs sur toute la chaîne de l'étain, de l'atome le plus léger (100Sn) au plus lourd (150Sn).

Ils cherchaient deux choses spécifiques :

1. La courbe en parabole : Entre deux points clés (les nombres magiques 50 et 82), la taille des noyaux devrait suivre une courbe douce, comme une colline.
2. Le "Coup de Couteau" (Kink) : À un moment précis, avec l'atome 132Sn, la taille devrait faire un saut brusque, comme si on donnait un coup de couteau dans la courbe. C'est un phénomène très important qui indique un changement dans la structure interne du noyau.

🎭 Le Résultat : Qui a gagné ?

C'est ici que ça devient intéressant, car aucune recette n'est parfaite :

• La recette classique a raté le coup : elle a prédit des noyaux trop petits et n'a pas vu le "coup de couteau" à 132Sn.
• La recette "magique" (EM7.5) a réussi à reproduire le "coup de couteau" à 132Sn ! C'est une victoire. MAIS, il y a un piège.
  • Pour obtenir ce résultat, cette recette a inventé une structure interne bizarre pour les atomes plus lourds (au-delà de 132Sn). Elle prédit un autre "coup de couteau" inversé à 142Sn qui semble faux.
  • L'analogie : C'est comme si un architecte avait réussi à construire un pont qui ne s'effondre pas, mais en utilisant des matériaux qui feraient s'effondrer le bâtiment d'à côté. Il a eu la bonne réponse pour la mauvaise raison.

🔍 Pourquoi c'est important ?

Cette étude nous apprend deux choses cruciales :

1. Nos théories sont encore imparfaites. Même les meilleures "recettes" actuelles ne peuvent pas tout prédire parfaitement en même temps (l'énergie, la taille, et la stabilité).
2. L'avenir est excitant. La région au-delà de l'atome 132Sn est un "terrain de jeu" idéal. C'est là que les différences entre les théories sont les plus grandes. Si on mesure la taille des atomes d'étain encore plus lourds (comme le 134Sn ou au-delà) et plus légers (vers 100Sn) en laboratoire, on pourra dire quelle recette est la vraie.

🚀 Conclusion Simple

Les physiciens ont utilisé des supercalculateurs pour simuler la taille des atomes d'étain. Ils ont découvert que leurs meilleures théories actuelles sont comme des lunettes qui corrigent la vue, mais qui déforment encore un peu les couleurs.

Pour avoir une vision parfaite de l'univers, il faut :

• De nouvelles mesures en laboratoire (pour avoir la vérité).
• Des calculs encore plus précis (en ajoutant des détails complexes que l'on ignorait).
• Et surtout, comprendre pourquoi certaines théories réussissent par hasard sur un point mais échouent sur le suivant.

C'est un peu comme essayer de deviner la recette exacte d'un plat secret en goûtant différents échantillons : on est proche, mais il manque encore un ingrédient clé pour que tout soit parfait !

Résumé technique

1. Problématique et Contexte

La description des noyaux atomiques à partir des principes premiers (ab initio) repose sur la théorie du champ effectif chiral ($\chi$EFT) pour les interactions nucléon-nucléon. Bien que ces méthodes aient réussi à décrire avec succès de nombreux observables dans les noyaux légers et de masse moyenne (énergies de liaison, spectres), la reproduction des rayons de charge reste un défi majeur.

• Le problème spécifique : Les interactions $\chi$EFT standard ont tendance à sous-estimer systématiquement les rayons de charge absolus. De plus, elles peinent à reproduire des comportements complexes tels que la courbure parabolique inversée entre $^{40}$Ca et $^{48}$Ca, ou la forte augmentation observée entre $^{48}$Ca et $^{52}$Ca.
• L'objectif : Tester la capacité des interactions $\chi$EFT de pointe à reproduire l'évolution des rayons de charge sur une longue chaîne isotopique lourde, en particulier à travers et au-delà du noyau doublement magique $^{132}$Sn (Z=50, N=82), où des mesures récentes ont révélé un "coudé" (kink) prononcé.

2. Méthodologie

Les auteurs ont employé une approche ab initio rigoureuse combinant plusieurs éléments théoriques et numériques :

• Formalisme : Utilisation de la théorie du Coupled Cluster de Bogoliubov (BCC). Cette méthode est adaptée aux noyaux à couches ouvertes (open-shell) où les corrélations d'appariement sont cruciales.
  • L'état fondamental est paramétré comme $|\Psi_0\rangle = e^T |\Phi\rangle$, où $|\Phi\rangle$ est un vide de Bogoliubov (solution des équations HFB) et $T$ est l'opérateur de cluster.
  • Les calculs sont effectués au niveau BCCSD (Singles and Doubles), incluant les excitations à 1 et 2 quasi-particules.
• Hamiltoniens : Trois interactions $\chi$EFT finement ajustées ("fine-tuned") ont été comparées :
  1. 1.8/2.0 (EM) : Précise pour les énergies de liaison mais sous-estime les rayons de charge (~5 %).
  2. $\Delta$NNLOGO : Inclut des degrés de liberté $\Delta$ et est ajusté sur les propriétés de saturation de la matière nucléaire pour mieux reproduire les rayons.
  3. 1.8/2.0 (EM7.5) : Une version récente où la constante de basse énergie $c_D$ a été refittée pour reproduire l'énergie et le rayon de $^{16}$O, améliorant ainsi les rayons dans les noyaux magiques.
• Opérateurs et Base :
  • Calculs dans une base d'oscillateur harmonique sphérique (13 couches majeures, $e_{max}=12$).
  • Le rayon quadratique moyen de charge est calculé en incluant les corrections de centre de masse, les tailles finies des protons et neutrons, les corrections spin-orbite et le terme relativiste Darwin-Foldy.
  • Les incertitudes de l'espace de modèle sont évaluées en variant la fréquence de l'oscillateur ($\hbar\omega$) et la taille de la base ($e_{max}$).

3. Résultats Clés

A. Rayons de charge absolus et écarts isotopiques

• Accord global : Les résultats BCCSD sont très cohérents avec des calculs antérieurs utilisant l'IMSRG (VS-IMSRG).
• Performance des Hamiltoniens :
  • L'interaction 1.8/2.0 (EM) sous-estime fortement les rayons expérimentaux (erreur RMS de 5,6 %).
  • Les interactions $\Delta$NNLOGO et 1.8/2.0 (EM7.5) corrigent cette déficience, réduisant l'erreur RMS à 1,2 % et 0,7 % respectivement.
• Rôle des corrélations : Bien que les corrélations BCCSD ne contribuent qu'à environ 1 % du rayon total, elles sont systématiquement nécessaires pour améliorer l'accord avec l'expérience. Elles augmentent les rayons lorsque le champ moyen est trop petit (cas 1.8/2.0 EM) et les diminuent lorsqu'il est trop grand (cas 1.8/2.0 EM7.5).

B. Le "Kink" à $N=82$ ($^{132}$Sn) et au-delà

• Reproduction du kink : L'augmentation brutale du rayon entre $^{132}$Sn et $^{134}$Sn (mesurée expérimentalement à 0,22 fm$^2$) n'est reproduite correctement que par l'interaction 1.8/2.0 (EM7.5) (0,23 fm$^2$). Les deux autres interactions ne capturent qu'environ la moitié de cette augmentation.
• Structure en couches critique : Cette réussite pour $^{134}$Sn est liée à une structure en couches spécifique prédite par 1.8/2.0 (EM7.5) : les neutrons au-delà de $N=82$ occupent l'orbitale $1h_{9/2}$, qui est spatialement plus étendue, tirant ainsi les protons vers l'extérieur.
• Le paradoxe de $^{142}$Sn : L'interaction 1.8/2.0 (EM7.5) prédit un "kink inversé" à $^{142}$Sn (une baisse du rayon isotopique) corrélée à une augmentation suspecte de l'énergie de séparation de deux neutrons ($S_{2n}$). Cela suggère que la reproduction correcte du kink à $^{132}$Sn par cet Hamiltonien est obtenue pour les mauvaises raisons (une structure de couches incorrecte au-delà de $N=82$).
• Sensibilité accrue : Au-delà de $^{132}$Sn, la sensibilité des résultats à l'Hamiltonien choisi est fortement amplifiée, offrant un test plus strict que dans la région $N=50-82$.

C. Composantes linéaires et paraboliques

• Entre $^{100}$Sn et $^{132}$Sn, les écarts isotopiques peuvent être décomposés en une composante linéaire et une composante quasi-parabolique.
• Les Hamiltoniens $\Delta$NNLOGO et 1.8/2.0 (EM7.5) reproduisent bien la pente linéaire, tandis que 1.8/2.0 (EM) est moins performant.
• La composante parabolique est bien reproduite par 1.8/2.0 (EM) et $\Delta$NNLOGO, mais sous-estimée par 1.8/2.0 (EM7.5). Des irrégularités locales (maxima/minima) apparaissent dans les calculs théoriques (liées aux sous-couches) qui ne sont pas visibles dans les données expérimentales, suggérant l'absence de corrélations collectives dans les calculs actuels.

4. Contributions et Significativité

• Test rigoureux des interactions $\chi$EFT : L'article démontre que la chaîne isotopique de l'étain, en particulier au-delà de $^{132}$Sn, constitue un "terrain de jeu" unique pour contraindre les paramètres des interactions nucléaires. Aucune des trois interactions testées ne parvient à capturer simultanément toutes les caractéristiques clés (énergies de liaison, rayons absolus, kinks, et structure de couches au-delà de $N=82$).
• Limites des approximations actuelles : La nécessité d'inclure des corrections de type "triples" (excitations à 3 quasi-particules) et des corrections de densité de charge à deux corps est soulignée pour atteindre une précision inférieure au pourcent, surtout pour les noyaux à couches ouvertes lourds.
• Appel à l'action expérimentale : Les auteurs recommandent vivement de mesurer les écarts isotopiques vers $^{100}$Sn et au-delà de $^{134}$Sn pour valider ou invalider les prédictions théoriques concernant la structure en couches dans ces régions extrêmes.

Conclusion

Cette étude met en lumière les défis persistants de la physique nucléaire ab initio pour les noyaux lourds. Bien que les interactions $\chi$EFT aient fait des progrès significatifs, la reproduction simultanée des énergies de liaison et des rayons de charge, ainsi que la prédiction correcte de l'évolution des structures de couches au-delà des nombres magiques, reste un objectif non atteint. Les résultats suggèrent que les propriétés critiques des interactions nucléon-nucléon doivent encore être affinées, en particulier dans le secteur des interactions à trois corps et des constantes de basse énergie.