Theoretical ab initio Evolution of Satellite Intensity near Threshold for Cu K-shell transitions

Résumé Technique : Évolution Théorique Ab Initio de l'Intensité des Satellites près du Seuil pour les Transitions de la Coquille K de Cu

Énoncé du Problème
Les calculs de structure atomique théorique pour les systèmes complexes reposent souvent sur l'approximation de la soudaineté (SA - sudden approximation), qui suppose que l'ionisation ou la désintégration nucléaire se produit sur un intervalle de temps beaucoup plus court que la relaxation électronique. Bien que la SA fournisse une excellente description des processus à haute énergie, sa validité est discutable près des seuils d'ionisation. Dans ce régime de basse énergie, l'électron émis possède une énergie cinétique faible et reste à proximité du système atomique pendant un temps comparativement long, violant ainsi les hypothèses fondamentales de la SA. Par conséquent, le système approche d'une limite adiabatique où la structure électronique évolue de manière continue. Les observations expérimentales antérieures de cet effet, particulièrement en spectroscopie d'absorption des rayons X (XAS), ont été limitées par des contraintes de résolution. De plus, bien que des modèles théoriques tels que ceux de Thomas et Roy existent pour prédire l'évolution de l'intensité de secousse (shake intensity ou raies satellites) du régime adiabatique au régime soudain, il est nécessaire de valider ces modèles par rapport à des données d'émission à haute résolution et de déterminer si les méthodes ab initio standards peuvent simuler avec précision cette évolution sans paramètres d'ajustement empiriques pour l'intensité.

Méthodologie
Les auteurs ont employé une approche ab initio complète utilisant des méthodes de configuration multi-diélectronique (MCDF) de pointe, utilisant spécifiquement le code MCDFGME. La méthodologie comprenait plusieurs étapes clés :

1. Calculs de Structure Atomique : Toutes les transitions de la coquille K du Cu, y compris les taux radiatifs et non radiatifs, ont été calculées. Les fonctions d'onde atomiques ont été optimisées par la méthode du niveau optimisé avec relaxation complète pour les taux radiatifs, tandis que les taux non radiatifs utilisaient des fonctions d'onde d'électrons du continuum sans relaxation.
2. Détermination de la Probabilité de Secousse (Shake Probability) : La probabilité des processus de montée (shake-up) et de décroissance (shake-off) a été calculée théoriquement. La probabilité totale de secousse a été dérivée du produit scalaire interne des fonctions d'onde radiales des états initial et final. Les auteurs ont distingué les probabilités de secousse à un seul électron (en utilisant une approche de distribution binomiale) des probabilités de secousse totales, en tenant compte des excitations monopolaires.
3. Génération de Spectre Synthétique : Un spectre synthétique pour les transitions K$\alpha_{1,2}$ du Cu a été généré en pondérant tous les chemins de désintégration possibles par leurs taux et poids statistiques calculés.
4. Modélisation de l'Évolution Énergétique : Pour simuler l'évolution de l'intensité des raies près du seuil, les auteurs ont développé une fonction de modulation. Cela impliquait la convolution de la largeur lorentzienne naturelle des niveaux ionisés avec un profil gaussien représentant la distribution d'énergie du faisceau d'excitation. Un « Profil Lorentzien Effectif » a été utilisé pour rendre compte de la nature asymétrique de la section efficace d'ionisation près du seuil.
5. Incorporation des Excitations Résonnantes : Reconnaissant que les données d'ionisation seules ne pouvaient expliquer les résultats expérimentaux dans la région proche du seuil, les auteurs ont incorporé des excitations résonnantes provenant des phases d'oxydes de Cu(I) et Cu(II) (spécifiquement les transitions $1s \to 3d$ et $1s \to 4p$) identifiées via des données de structure fine d'absorption des rayons X (XANES).
6. Ajustement et Analyse : Les résultats théoriques ont été comparés aux données de Galambosi et al. [8]. Une analyse bayésienne utilisant l'échantillonnage imbriqué (nested sampling) a été employée pour déterminer les décalages d'énergie et les multiplicateurs de largeur les plus probables pour différents processus de désintégration (diagramme, shake-off, shake-up et excitation) afin d'aligner la simulation avec les observations expérimentales.

Contributions Clés

• Validation des Méthodes Ab Initio : L'étude démontre que les méthodes ab initio standards de pointe, sans paramètres d'ajustement d'intensité réglables pour l'ionisation, peuvent simuler avec succès l'évolution de l'intensité des satellites près des seuils d'ionisation.
• Quantification des Probabilités de Secousse : Les auteurs ont calculé une probabilité totale de secousse pour l'ionisation de la coquille K du Cu d'environ 25,97 %, avec une décomposition de $\approx$13,08 % pour le shake-off et $\approx$12,89 % pour le shake-up.
• Identification des Origines de Seuil : Le travail identifie que l'intensité de satellite « sous le seuil » observée dans les expériences ne provient pas uniquement des processus de secousse, mais est significativement influencée par les excitations résonnantes dans les phases d'oxydes de Cu(I) ($1s \to 4p$) et Cu(II) ($1s \to 3d$).
• Raffinement des Modèles Théoriques : Le papier constitue un test rigoureux des modèles de Thomas et Roy. Il confirme la forme fonctionnelle du modèle de Thomas pour la transition adiabatique-soudaine et souligne les divergences des paramètres du modèle de Roy lorsqu'ils sont ajustés aux données expérimentales, suggérant que les données expérimentales contiennent des processus physiques (tels que des effets de l'état solide) non pleinement capturés par les modèles atomiques simples.

Résultats

• Accord avec l'Expérience : Les rapports d'intensité simulés des raies satellites par rapport à l'intensité totale de K$\alpha_{1,2}$ ont montré un bon accord avec les données expérimentales pour les énergies d'excitation supérieures à 9050 eV, où tous les canaux de désintégration sont ouverts. Le rapport simulé à 10 keV a été trouvé à $0,262 \pm 0,00004$, ce qui est cohérent avec la probabilité totale de secousse calculée dans l'approximation de la soudaineté.
• Comportement de Seuil : En dessous de 9050 eV, la simulation basée uniquement sur l'ionisation n'a pas réussi à reproduire la montée de l'intensité des satellites observée expérimentalement. Cependant, en incluant les excitations résonnantes des oxydes de Cu, la simulation a réussi à reproduire les données expérimentales dans la plage 8975–9060 eV.
• Ajustement de Modèle : L'ajustement du modèle de Thomas à la simulation a produit des paramètres ($\Delta E$, $r$, $\mu_\infty$) qui sont en bon accord avec l'ajustement expérimental et les calculs théoriques. L'ajustement du modèle de Roy a montré une certaine proximité avec les valeurs théoriques mais des écarts significatifs dans les paramètres d'énergie de liaison, indiquant les limites du modèle pour capturer la pleine complexité du spectre expérimental.
• Caractéristiques Spectrales : L'inclusion de transitions de shake-up d'ordre supérieur (par exemple, shake-off de $3s$, $3p$, $2p$ et shake-up de $3d$) a considérablement modifié la forme de la raie satellite théorique par rapport aux calculs Dirac-Fock précédents, soulignant la nécessité de calculs de secousse détaillés pour la métrologie à haute résolution.

Signification
L'article affirme que sa principale importance réside dans la démonstration que les calculs théoriques standards peuvent prédire l'évolution moyenne de l'intensité spectrale près des seuils d'ionisation sans dépendre de l'ajustement empirique de l'intensité. En séparant avec succès l'intensité de satellite (provenant des processus de secousse) de l'intensité directe du diagramme et en identifiant la contribution des excitations résonnantes dans les phases d'oxydes, les auteurs fournissent une image physique plus complète de la région proche du seuil. Ce travail souligne l'importance des calculs détaillés de shake-up et de shake-off pour la spectroscopie de rayons X de pointe et suggère que les divergences entre théorie et expérience dans ce régime peuvent souvent être attribuées à des effets de l'état solide (tels que la présence d'oxydes) et aux formes de raie spécifiques utilisées dans l'analyse des données, plutôt qu'à un échec des modèles de physique atomique sous-jacents.