Connecting Diffuse Scattering to Atomic-Site-Resolved Occupancy and Displacement Fields through Fourier Filtering

Résumé technique : Relier la diffusion diffuse à l'occupation des sites atomiques et aux champs de déplacement par filtrage de Fourier

Énoncé du problème
Les corrélations structurelles locales dans les matériaux cristallins, telles que l'ordre chimique et les déplacements atomiques, sont encodées dans les motifs de diffusion diffuse observés en diffraction de rayons X, de neutrons ou d'électrons. Bien que l'analyse cristallographique traditionnelle utilisant la diffraction de Bragg capture la structure moyenne, elle échoue largement à accéder aux déviations locales responsables des propriétés fonctionnelles. Un défi important subsiste dans l'interprétation de la diffusion diffuse : même lorsqu'un modèle structurel atomique (dérivé de raffinements par Monte Carlo inversé, de dynamique moléculaire ou de simulations) reproduit avec succès l'intensité diffuse expérimentale, il reste difficile d'identifier quels motifs atomiques corrélés spécifiques au sein de ce modèle génèrent des caractéristiques diffuses particulières. De plus, distinguer les contributions chevauchantes de l'ordre chimique et des déplacements statiques, ou cartographier ces caractations vers des sites atomiques spécifiques dans des modèles à grande échelle, est non trivial. Les méthodes existantes, telles que la 3D-$\Delta$PDF ou le filtrage de Fourier qualitatif d'images de microscopie, manquent souvent d'un moyen quantitatif pour déterminer les champs de corrélation chimique ou de déplacement avec une résolution de site.

Méthodologie : Le cadre MOSAIC
Les auteurs introduisent MOSAIC (Mapping of Occupancy and Site-resolved Atomic Interactions via Correlation), un cadre computationnel conçu pour lier des caractéristiques de diffusion diffuse spécifiques à des champs d'occupation et de déplacement résolus par site dans une configuration atomique connue. La méthode ne reconstruit pas une structure à partir de données d'intensité seules ; elle interroge un modèle structurel fourni pour calculer l'amplitude de diffusion et isoler des caractéristiques spécifiques.

Le flux de travail procède comme suit :

1. Calcul de l'amplitude de diffusion : À partir d'une configuration atomique, l'amplitude de diffusion complexe totale $A_{tot}(\mathbf{Q})$ est calculée. L'amplitude de diffusion diffuse, $\Delta A(\mathbf{Q})$, est obtenue en soustrayant la contribution de la structure cristalline moyenne.
2. Filtrage de Fourier : Des régions spécifiques de l'espace réciproque (par exemple, pics, bâtonnets ou lobes diffus) sont isolées à l'aide de masques préservant la phase ($W_F(\mathbf{Q})$). L'amplitude masquée subit ensuite une transformée de Fourier inverse pour générer des champs dans l'espace réel.
3. Séparation des modes : Le cadre emploie deux modes distincts pour démêler les contributions chimiques et displacives :
   • Mode Chimique : Une configuration de référence est construite où tous les atomes sont placés sur leurs positions cristallographiques moyennes, mais la distribution des espèces chimiques est conservée. L'amplitude diffuse calculée à partir de cette configuration de « Mode Chimique » isole la diffusion provenant uniquement des variations d'occupation (ordre chimique), supprimant ainsi les effets de déplacement.
   • Mode de Déplacement : Ce mode cible la partie de l'amplitude de diffusion sensible aux déplacements. Comme le champ filtré dans l'espace réel est une convolution des déplacements discrets avec la fonction de étalement ponctuel du masque, les auteurs emploient un « M-décodeur » linéaire. Ce décodeur de données est un opérateur linéaire qui mappe des patchs locaux du champ filtré dans l'espace réel vers des vecteurs de déplacement résolus par site. Le décodeur est entraîné sur des données appariées (champ filtré vs déplacement réel connu) pour la configuration spécifique, puis appliqué uniformément.
4. Scalabilité : Pour gérer de grands ensembles atomiques (des millions d'atomes) et des régions étendues de l'espace réciproque, l'implémentation utilise une architecture Map–Reduce et des transformées de Fourier rapides non uniformes (NUFFT). Cela évite les goulots d'étranglement de la mémoire des FFT globales basées sur une grille et permet un traitement efficace de maillages de points Q arbitraires.

Résultats clés et validation
L'article valide le cadre MOSAIC à travers des simulations et des applications à des modèles expérimentaux :

• Séparation de l'ordre chimique : Dans un système simulé de $\gamma$-LiFeO$_2$ contenant la coexistence d'un ordre à courte portée (SRO) et d'un ordre L1$_1$ à l'échelle nanométrique, MOSAIC a réussi à séparer ces deux effets. En appliquant des masques à des régions spécifiques de l'espace réciproque (les réflexions au point L vs la surface diffuse continue), la méthode a reconstruit des cartes distinctes dans l'espace réel : l'une mettant en évidence les nanocristaux L1$_1$ et l'autre révélant le SRO de la matrice, bien que les deux contribuent au même motif de diffusion globale.
• Démêlage des distorsions : Dans un type de pérovskite de type ReO$_3$ présentant la coexistence de rotations octaédriques et de distorsions de type « breathing », le cadre a séparé les caractéristiques diffuses chevauchantes. Des masques cylindriques ont isolé le désordre rotationnel (bâtonnets diffus), tandis que des masches sphériques ont isolé le mode de respiration (pics de surstructure). La linéarité additive de la méthode a été confirmée, car la somme des champs de déplacement séparés a reconstruit le champ de distorsion total d'origine avec un résidu d'erreur négligeable.
• Application expérimentale : Le cadre a été appliqué à des modèles affinés par RMC de ferroélectriques relaxeurs PbMg$_{1/3}$Nb$_{2/3}$O$_3$ (PMN) et PMN–PbTiO$_3$ (PMN-PT).
  • Chimique : Le Mode Chimique a révélé des variations continues du degré d'ordre de type rocsalt et a identifié des joints de grains anti-phase, fournissant une vue résolue par site des environnements chimiques qui étaient obscurcis dans les modèles bruts.
  • Déplacement : Le Mode de Déplacement a filtré le bruit non corrélé des champs de déplacement de Pb, révélant des textures polaires à l'échelle nanométrique et des motifs de polarisation tourbillonnants cohérents avec les observations précédentes.

Signification et affirmations
Les auteurs affirment que MOSAIC fournit une méthode robuste et quantitative pour attribuer des caractéristiques de diffusion diffuse spécifiques à des champs d'occupation et de déplacement résolus par site atomique. Les principales contributions sont :

• Découplage quantitatif : La capacité de séparer les corrélations chimiques et displacives qui produisent des caractéristiques de diffusion chevauchantes, une tâche difficile avec l'analyse traditionnelle.
• Linéarité et additivité : Le cadre garantit une linéarité stricte, ce qui signifie que les contributions de régions disjointes de l'espace réciproque s'ajoutent exactement dans l'espace réel. Cela permet la décomposition rigoureuse de motifs structurels complexes en constituants.
• Scalabilité : L'utilisation des NUFFT et du traitement Map–Reduce permet l'analyse de modèles structurels à grande échelle (provenant par exemple de RMC ou de MD) qui étaient auparavant informatiquement prohibitifs pour un filtrage de caractéristiques aussi détaillé.
• Polyvalence : La méthode est applicable à la fois à des configurations atomiques 3D dérivées de données de diffusion ou de simulations, ainsi qu'à des projections structurelles 2D de microscopie électronique à résolution atomique.

L'article souligne que MOSAIC est un outil d'analyse de configurations connues pour comprendre les origines de caractéristiques de diffusion spécifiques, plutôt qu'une méthode d'inversion directe pour résoudre des structures à partir de données d'intensité seules. Il offre une voie pour lier quantitativement les motifs structurels locaux aux réponses fonctionnelles macroscopiques dans les matériaux complexes.