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

Technische Samenvatting: Het koppelen van diffuse verstrooiing aan bezettings- en verplaatsingsvelden op atomaire locatie via Fourier-filtering

Probleemstelling
Lokale structurele correlaties in kristallijne materialen, zoals chemische ordening en atomaire verplaatsingen, zijn gecodeerd in diffuse verstrooiingspatronen die worden waargenomen bij röntgen-, neutronen- of elektronen diffractie. Hoewel traditionele kristallografische analyse met behulp van Bragg-diffractie de gemiddelde structuur vastlegt, slaagt deze grotendeels niet in het ontsluiten van de lokale afwijkingen die verantwoordelijk zijn voor functionele eigenschappen. Er bestaat een aanzienlijke uitdaging bij het interpreteren van diffuse verstrooiing: zelfs wanneer een atomistisch structureel model (afgeleid van Reverse Monte Carlo-verfijningen, moleculaire dynamica of simulaties) de experimentele diffuse intensiteit succesvol reproduceert, blijft het moeilijk om te identificeren welke specifieke gecorreleerde atomische motieven binnen dat model bepaalde diffuse kenmerken genereren. Bovendien is het onderscheiden tussen overlappende bijdragen van chemische ordening en statische verplaatsingen, of het mappen van deze kenmerken naar specifieke atomaire locaties in grootschalige modellen, niet triviaal. Bestaande methoden, zoals 3D-$\Delta$PDF of kwalitatieve Fourier-filtering van microscopiebeelden, missen vaak een kwantitatieve methode om de onderliggende chemische of verplaatsingscorrelatievelden met resolutie op de locatie van de atomaire site te bepalen.

Methodologie: Het MOSAIC-raamwerk
De auteurs introduceren MOSAIC (Mapping of Occupancy and Site-resolved Atomic Interactions via Correlation), een computationeel raamwerk dat is ontworpen om specifieke diffuse-verstrooiingskenmerken te koppelen aan bezettings- en verplaatsingsvelden op de locatie van de atomaire site in een bekende atomaire configuratie. De methode reconstrueert niet een structuur uit intensiteitsgegevens alleen; in plaats daarvan onderzoekt het een aangeleverd structureel model om de verstrooiingsamplitude te berekenen en specifieke kenmerken te isoleren.

De workflow verloopt als volgt:

1. Berekening van de verstrooiingsamplitude: Vertrekkend vanuit een atomaire configuratie wordt de totale complexe verstrooiingsamplitude $A_{tot}(\mathbf{Q})$ berekend. De diffuse verstrooiingsamplitude, $\Delta A(\mathbf{Q})$, wordt verkregen door de bijdrage van de gemiddelde kristallografische structuur af te trekken.
2. Fourier-filtering: Specifieke regio's in de reciproke ruimte (bijv. diffuse pieken, staven of lobben) worden geïsoleerd met behulp van fase-behoudende maskers ($W_F(\mathbf{Q})$). De gemaskeerde amplitude wordt vervolgens invers Fourier getransformeerd om real-space velden te genereren.
3. Modus-scheiding: Het raamwerk maakt gebruik van twee afzonderlijke modi om chemische en displacieve bijdragen te ontwarren:
   • Chemische modus: Een referentieconfiguratie wordt geconstrueerd waarbij alle atomen op hun gemiddelde kristallografische posities worden geplaatst, maar de verdeling van de chemische soorten behouden blijft. De berekende diffuse amplitude uit deze "Chemische Modus"-configuratie isoleert de verstrooiing die uitsluitend voortkomt uit bezettingsvariaties (chemische ordening), waardoor effecten van verplaatsing worden onderdrukt.
   • Verplaatsingsmodus: Deze modus richt zich op het verplaatsingsgevoelige deel van de diffuse amplitude. Omdat het inverse getransformeerde veld een convolutie is van discrete verplaatsingen met de puntspreidingsfunctie van het masker, gebruiken de auteurs een lineaire "M-decoder". Deze data-gestuurde lineaire operator brengt lokale segmenten van het gefilterde real-space veld in kaart naar site-geresolveerde verplaatsingsvectoren. De decoder wordt getraind op gepaarde data (gefilterd veld versus bekende ware verplaatsing) voor de specifieke configuratie en wordt vervolgens uniform toegepast.
4. Schaalbaarheid: Om grote atomaire ensembles (miljoenen atomen) en uitgebreide regio's in de reciproke ruimte aan te kunnen, maakt de implementatie gebruik van een Map–Reduce-architectuur en Non-Uniform Fast Fourier Transforms (NUFFTs). Dit voorkomt de geheugenbottlenecks van globale grid-gebaseerde FFT's en maakt efficiënte verwerking van willekeurige Q-punt meshes mogelijk.

Belangrijkste resultaten en validatie
Het artikel valideert het MOSAIC-raamwerk door middel van simulaties en toepassingen op experimentele modellen:

• Scheiding van chemische ordening: In een gesimuleerd $\gamma$-LiFeO$_2$-systeem dat coëxisterende kort-bereik orde (SRO) en nanoscale L1$_1$-ordening bevat, slaagde MOSAIC erin deze twee effecten te scheiden. Door maskers toe te passen op specifieke regio's in de reciproke ruimte (de L-punt reflecties versus het continue diffuse oppervlak), reconstrueerde de methode onderscheidende real-space kaarten: één die de L1$_1$-nanoclusters accentueerde en de andere die de matrix SRO onthulde, ondanks het feit dat beide bijdragen aan hetzelfde totale verstrooiingspatroon.
• Ontwarren van distorties: In een gesimuleerd ReO$_3$-type perovskiet met coëxisterende octaëdrische rotaties en ademhalingsdistorties (breathing distortions), scheidde het raamwerk de overlappende diffuse kenmerken. Cilindrische maskers isoleerden de rotationele wanorde (diffuse staven), terwijl sferische maskers de ademhalingsmodus (superlattice pieken) isoleerden. De lineaire additiviteit van de methode werd bevestigd, aangezien de som van de gescheiden verplaatsingsvelden de oorspronkelijke totale distortieveld reconstrueerde met een verwaarloosbare residuele fout.
• Experimentele toepassing: Het raamwerk werd toegepast op RMC-verfijnde modellen van relaxor-ferro-elektrische materialen PbMg$_{1/3}$Nb$_{2/3}$O$_3$ (PMN) en PMN–PbTiO$_3$ (PMN-PT).
  • Chemisch: De Chemische Modus onthulde continue variaties in de mate van rochtsalt-type ordening en identificeerde anti-fase grenzen, wat een site-geresolveerd beeld gaf van chemische omgevingen die in de ruwe modellen verborgen waren.
  • Verplaatsing: De Verplaatsingsmodus filterde ongecorreleerde ruis uit de Pb-verplaatsingsvelden, waardoor nanoscale polaire texturen en draaiende polarisatiepatronen werden onthuld die consistent zijn met eerdere observaties.

Betekenis en claims
De auteurs beweren dat MOSAIC een robuuste, kwantitatieve methode biedt voor het toeschrijven van specifieke diffuse-verstrooiingskenmerken aan onderliggende atomaire-locatie-geresolveerde bezettings- en verplaatsingsvelden. Belangrijke bijdragen zijn:

• Kwantitatieve ontkoppeling: Het vermogen om chemische en displacieve correlaties te scheiden die overlappende diffuse verstrooiingskenmerken produceren, een taak die met traditionele analyse moeilijk is.
• Lineariteit en additiviteit: Het raamwerk garandeert strikte lineariteit, wat betekent dat bijdragen van disjuncte regio's in de reciproke ruimte exact optellen in de real-space. Dit maakt de rigoureuze decompositie van complexe verstrooiingspatronen in constituerende structurele motieven mogelijk.
• Schaalbaarheid: Het gebruik van NUFFTs en Map–Reduce-verwerking maakt de analyse van grootschalige structurele modellen (bijv. uit RMC of MD) mogelijk die voorheen computationeel onhaalbaar waren voor dergelijke gedetailleerde filtering van kenmerken.
• Veelzijdigheid: De methode is toepasbaar op zowel 3D atomaire configuraties afgeleid van verstrooiingsgegevens of simulaties, als op 2D structurele projecties van atomaire-resolutie elektronenmicroscopie.

Het artikel benadrukt dat MOSAIC een hulpmiddel is voor het analyseren van bekende configuraties om de oorsprong van specifieke verstrooiingskenmerken te begrijpen, in plaats van een directe inversiemethode om structuren op te lossen uit enkel intensiteitsgegevens. Het biedt een pad om lokale structurele motieven kwantitatief te koppelen aan macroscopische functionele responsen in complexe materialen.