Determination of nuclear quadrupole moments for $^{25}$Mg, $^{87}$Sr, and $^{135,137}$Ba via configuration-interaction combined with a coupled-cluster approach

Cet article emploie une approche d'interaction de configuration plus de clusters couplés pour calculer les gradients de champ électrique et les constantes de structure hyperfine de dipôle magnétique pour les états de faible énergie du Mg, du Sr et du Ba, permettant la détermination précise des moments quadrupolaires nucléaires pour $^{25}$Mg, $^{87}$Sr et $^{135,137}$Ba qui révèlent des écarts significatifs avec les valeurs précédemment adoptées pour le strontium et le baryum.

L'essentiel

Imaginez le noyau d'un atome non pas comme une boule parfaitement ronde, mais comme un ballon légèrement écrasé ou étiré. Cette forme n'est pas le fruit du hasard ; c'est une « signature » spécifique de la façon dont les protons et les neutrons sont disposés à l'intérieur. Les scientifiques appellent cette forme le moment quadripolaire nucléaire. Connaître la forme exacte de ce « ballon » est crucial pour comprendre les règles fondamentales de la physique, de la façon dont les atomes s'assemblent au comportement des matériaux.

Cependant, mesurer cet écrasement directement est incroyablement difficile. C'est comme essayer de deviner la forme exacte d'un ballon à l'intérieur d'une boîte scellée et obscure, simplement en écoutant le son qu'il produit lorsqu'on le tapote.

L'expérience : Écouter le « bourdonnement » atomique

Dans cet article, l'auteur, Yong-Bo Tang, agit comme un maître ingénieur du son. Il se concentre sur trois « familles » spécifiques d'atomes : le Magnésium (Mg), le Strontium (Sr) et le Baryum (Ba).

Lorsque ces atomes sont excités (comme une corde de guitare que l'on pince), ils émettent un « bourdonnement » ou une vibration très spécifique appelée structure hyperfine. Ce bourdonnement est causé par l'interaction entre les électrons en rotation à l'extérieur et le noyau écrasé à l'intérieur.

• La partie mesurée : Les scientifiques ont déjà mesuré la hauteur de ce bourdonnement de manière très précise en laboratoire.
• Le maillon manquant : Pour déduire la forme du noyau (le moment quadripolaire) à partir de cette hauteur, il faut savoir exactement comment les électrons sont disposés autour du noyau. Cet arrangement crée un « gradient de champ électrique » (pensez à la pente d'une colline sur laquelle les électrons dévalent).

Le problème : La colline est trop escarpée pour être calculée

Calculer la forme de cette « colline » (le gradient de champ électrique) est un cauchemar pour les ordinateurs. Les électrons ne restent pas immobiles ; ils dansent les uns autour des autres, se poussant et se tirant de manières complexes, ce qu'on appelle la corrélation électronique.

• Si vous ignorez ces danses, votre calcul de la colline est faux.
• Si vous essayez de calculer chaque danse, votre ordinateur plante.

Les tentatives précédentes pour calculer cela étaient comme essayer de cartographier une chaîne de montagnes à l'aide d'une photo satellite floue. Les résultats étaient incohérents. Pour le Strontium et le Baryum, différentes études donnaient des réponses différentes, avec des écarts allant jusqu'à 10 %.

La solution : Une approche hybride de type « Couteau Suisse »

Pour résoudre ce problème, Tang a développé une nouvelle méthode de calcul qui combine deux techniques puissantes :

1. L'Interaction de Configuration (CI) : C'est comme examiner chaque façon possible dont les électrons pourraient se disposer, une par une. C'est approfondi, mais lent.
2. Le Coupled-Cluster (CC) : C'est comme utiliser un raccourci sophistiqué pour prédire comment les électrons s'influencent mutuellement en groupes. C'est rapide, mais cela manque parfois les détails infimes.

La méthode de Tang, CI+CC, est le meilleur des deux mondes. Elle utilise le « raccourci » pour gérer les interactions massives entre les électrons du cœur, puis utilise la méthode « approfondie » pour affiner les détails des électrons périphériques. C'est comme utiliser un drone pour cartographier la forme générale d'une forêt, puis envoyer une équipe de randonneurs pour mesurer la hauteur exacte de chaque arbre spécifique.

Les résultats : Clarifier la confusion

En utilisant ce « couteau suisse » de haute précision, Tang a calculé les gradients de champ électrique pour plusieurs états de basse énergie du Mg, du Sr et du Ba. Il a ensuite combiné ses calculs avec les « bourdonnements » expérimentaux connus pour déterminer les formes nucléaires.

Voici ce qu'il a trouvé :

• Magnésium (25Mg) : Le résultat correspond parfaitement aux expériences précédentes. C'est comme accorder une radio et trouver la station de façon cristalline. La forme calculée concorde avec ce qui a été trouvé via des expériences de « rayons X muoniques » (une autre méthode de haute technologie pour mesurer).
• Strontium (87Sr) : Ici, l'intrigue se corse. Le résultat de Tang suggère que le noyau est environ 10 % plus écrasé que la valeur actuellement acceptée dans les manuels. L'ancienne valeur provenait de l'observation d'un ion de Strontium (un atome ayant perdu un électron), tandis que Tang a observé l'atome neutre. Cette différence suggère que l'ancien calcul a pu manquer certaines danses subtiles des électrons.
• Baryum (135,137Ba) : Similairement au Strontium, les résultats de Tang pour le Baryum diffèrent d'environ 4 % des valeurs actuellement acceptées dérivées des ions de Baryum.

Ce qu'il faut retenir

L'article conclut que, bien que la méthode fonctionne magnifiquement pour le Magnésium, il existe une divergence significative pour le Strontium et le Baryum par rapport aux valeurs considérées comme le « standard d'or » actuellement utilisées par les scientifiques.

Tang suggère que la différence pourrait provenir du fait que les calculs actuels du « standard d'or » ont omis un type spécifique d'interaction électronique appelé excitation triple (où trois électrons interagissent simultanément). Tout comme une chorale sonne différemment si trois chanteurs harmonisent d'une manière que personne n'avait prévue, ces interactions triples pourraient modifier la « hauteur » de la forme de l'atome.

En résumé : L'auteur a construit un meilleur modèle informatique pour mesurer la forme des noyaux atomiques. Pour le Magnésium, le modèle a confirmé ce que nous savions déjà. Pour le Strontium et le Baryum, le modèle suggère que les mesures « officielles » actuelles pourraient être légèrement erronées, laissant entendre que nous devons examiner de plus près comment trois électrons interagissent pour obtenir la véritable forme de ces noyaux atomiques.

Résumé technique

Résumé technique : Détermination des moments quadripolaires nucléaires pour $^{25}\text{Mg}$, $^{87}\text{Sr}$ et $^{135,137}\text{Ba}$

Énoncé du problème
Le moment quadripolaire électrique nucléaire ($Q$) est un paramètre fondamental décrivant la symétrie de la structure nucléaire, avec des implications significatives en physique atomique, nucléaire et de la matière condensée. Bien que les constantes de structure hyperfine expérimentales ($B$) puissent être mesurées avec une grande précision, l'extraction de valeurs précises de $Q$ nécessite des calculs théoriques des gradients de champ électrique ($q$) tout aussi précis. La littérature existante révèle des écarts notables dans les valeurs recommandées pour les moments quadripolaires nucléaires de $^{87}\text{Sr}$ et $^{137}\text{Ba}$, particulièrement lors de la comparaison entre les valeurs dérivées d'atomes neutres et d'ions. Par exemple, les valeurs dérivées de l'atome de Sr neutre sont environ 8 % plus grandes que celles recommandées à partir des données d'ions, et des incohérences similaires existent pour le Ba. De plus, les approches théoriques standards peinent souvent à rendre compte des effets complexes de corrélation électronique nécessaires pour des calculs de gradient de haute précision dans les systèmes alcalino-terreux lourds.

Méthodologie
Les auteurs utilisent une approche hybride relativiste combinant la méthode de l'Interaction de Configuration (CI) avec la méthode de Coupled-Cluster (CC), appelée méthode RCI+all-order. Ce cadre est conçu pour rendre compte de manière exhaustive des corrélations électroniques cœur-cœur, cœur-valence et valence-valence, qui sont décisives pour des calculs précis de gradient de champ électrique.

Composants méthodologiques clés :

• Hamiltonien et interaction : L'interaction hyperfine est modélisée à l'aide d'opérateurs de tenseurs sphériques pour les interactions dipolaires magnétiques ($k=1$) et quadripolaires électriques ($k=2$). Le noyau est traité comme une sphère uniformément magnétisée et chargée.
• Schéma de calcul hybride : La méthode utilise un ensemble de base fini de fonctions de type gaussiennes à température uniforme pour étendre les fonctions de Dirac. Les potentiels de corrélation sont dérivés via cinq stratégies de calcul distinctes pour évaluer l'incertitude :
  1. CI + Théorie de la perturbation de corps multiples du second ordre (MBPT(2)).
  2. CI + Linearized Coupled-Cluster Singles and Doubles (LCCSD).
  3. CI + Full Coupled-Cluster Singles and Doubles (CCSD).
  4. CI + LCCSD avec paramètres de mise à l'échelle (CI+LCCSDs).
  5. CI + CCSD avec paramètres de mise à l'échelle (CI+CCSDs).
• Corrections : Les éléments de matrice de transition incluent l'approximation de la phase aléatoire (RPA), les effets de Brueckner du cœur, la radiation structurelle et les corrections de normalisation (collectivement appelées corrections d'ordre élevé ou HO), ainsi que les interactions à deux particules (TP) jusqu'au second ordre.
• Extraction de Q : Le moment quadripolaire nucléaire est extrait en utilisant la relation $Q = B / (234.9648867 \cdot q)$, où $B$ est la constante de structure hyperfine mesurée expérimentalement et $q$ est le gradient de champ électrique calculé théoriquement.

Contributions clés et résultats
L'étude calcule systématiquement les gradients de champ électrique et les constantes de structure hyperfine dipolaire magnétique ($A$) pour les états de faible énergie du $^{25}\text{Mg}$, $^{87}\text{Sr}$ et $^{135,137}\text{Ba}$ neutres. Les états spécifiques analysés sont :

• Mg : $3s3p \ ^3P_{1,2}$
• Sr : $5s4d \ ^1D_2$ et $5s5p \ ^3P_{1,2}$
• Ba : $6s5d \ ^3D_{1,2,3}$, $6s5d \ ^1D_2$, et $6s6p \ ^1P_1$

Énergie et constantes hyperfines :
Les énergies calculées et les constantes dipolaires magnétiques ($A$) montrent un excellent accord avec les données expérimentales (NIST) et d'autres résultats théoriques. La méthode CI+LCCSD est identifiée comme la plus fiable, les valeurs finales recommandées pour $A$ correspondant aux mesures expérimentales dans leurs incertitudes (typiquement à moins de 2 % pour la plupart des états). Cela valide l'approche hybride CI+CC pour les systèmes alcalino-terreux.

Moments quadripolaires nucléaires ($Q$) :
En combinant les gradients calculés avec les valeurs expérimentales de $B$, les auteurs déterminent les moments quadripolaires nucléaires suivants :

• $^{25}\text{Mg}$ : $0,203(2)$ b. Ce résultat est en parfait accord avec l'expérience de rayons X mésoniques ($0,201(3)$ b) et les déterminations antérieures basées sur l'hyperfine.
• $^{87}\text{Sr}$ : $0,336(4)$ b. Cette valeur est environ 10 % plus grande que la valeur actuellement adoptée de $0,305(2)$ b (dérivée de l'état $4d_{5/2}$ de Sr$^+$), mais s'aligne bien avec les valeurs dérivées de l'état $5p_{3/2}$ de Sr$^+$ et des états de Sr neutre ($5s5p \ ^3P_{1,2}$).
• $^{137}\text{Ba}$ : $0,246(4)$ b. C'est environ 4 % de plus que la valeur recommandée de $0,236(3)$ b (dérivée des états $5d_{3/2,5/2}$ de Ba$^+$). Cependant, cela concorde excellemment avec les valeurs dérivées de l'état $6p_{3/2}$ de Ba$^+$.
• $^{135}\text{Ba}$ : $0,161(2)$ b.

Signification et revendications
L'article affirme que l'approche d'interaction de configuration plus coupled-cluster fournit un cadre robuste et précis pour déterminer les moments quadripolaires nucléaires, particulièrement pour les noyaux instables ou lourds où la mesure directe est difficile.

Les auteurs soulignent un écart significatif entre les valeurs dérivées des atomes neutres et celles des ions pour le Sr et le Ba. Plus précisément, les résultats pour $^{87}\text{Sr}$ et $^{137}\text{Ba}$ dérivés des atomes neutres dans ce travail diffèrent d'environ 10 % et 4 %, respectivement, des valeurs actuellement adoptées dans la littérature (qui reposent sur des données d'ions). Les auteurs suggèrent que ces écarts peuvent provenir de l'omission des contributions de triple excitation dans les calculs standard de coupled-cluster, notant que des travaux récents sur les ions indiquent que les triples excitations contribuent à hauteur de 1 à 2 % aux constantes hyperfines.

L'étude conclut que, bien que leurs résultats pour $^{25}\text{Mg}$ soient cohérents avec les données établies, leurs conclusions pour $^{87}\text{Sr}$ et $^{137}\text{Ba}$ remettent en question les valeurs actuellement adoptées dérivées des états d'ions. Les auteurs affirment que leurs extractions basées sur les atomes neutres sont fiables, soutenues par la haute précision de leurs calculs de constantes dipolaires magnétiques, mais reconnaissent que des travaux théoriques supplémentaires incluant les triples excitations et les interactions à trois corps sont nécessaires pour résoudre pleinement les différences entre les déterminations basées sur les atomes neutres et les ions.