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

Resumen Técnico: Conexión de la Dispersión Difusa con Campos de Ocupación y Desplazamiento Resueltos por Sitio Atómico mediante Filtrado de Fourier

Planteamiento del Problema
Las correlaciones estructurales locales en materiales cristalinos, tales como el ordenamiento químico y los desplazamientos atómicos, están codificadas en los patrones de dispersión difusa observados en la difracción de rayos X, neutrones o electrones. Mientras que el análisis cristalográfico tradicional utilizando la difracción de Bragg captura la estructura promedio, este falla en gran medida al intentar acceder a las desviaciones locales responsables de las propiedades funcionales. Existe un desafío significativo en la interpretación de la dispersión difusa: incluso cuando un modelo estructural atomístico (derivado de refinamientos de Monte Carlo Inverso, dinámica molecular o simulaciones) reproduce con éxito la intensidad difusa experimental, sigue siendo difícil identificar qué motivos atómicos correlacionados específicos dentro de ese modelo generan características difusas particulares. Además, distinguir entre las contribuciones superpuestas del ordenamiento químico y los desplazamientos estáticos, o mapear estas características a sitios atómicos específicos en modelos de gran escala, no es trivial. Los métodos existentes, como el 3D-$\Delta$PDF o el filtrado de Fourier cualitativo de imágenes de microscopía, a menudo carecen de un medio cuantitativo para determinar los campos de correlación química o de desplazamiento con resolución de sitio.

Metodología: El Marco MOSAIC
Los autores presentan MOSAIC (Mapping of Occupancy and Site-resolved Atomic Interactions via Correlation), un marco computacional diseñado para vincular características específicas de la dispersión difusa con campos de ocupación y de desplazamiento resueltos por sitio en una configuración atómica conocida. El método no reconstruye una estructura a partir de datos de intensidad únicamente; más bien, interroga una configuración estructural suministrada para calcular la amplitud de dispersión e aislar características específicas.

El flujo de trabajo procede de la siguiente manera:

1. Cálculo de la Amplitud de Dispersión: Partiendo de una configuración atómica, se calcula la amplitud de dispersión compleja total $A_{tot}(\mathbf{Q})$. La amplitud de dispersión difusa, $\Delta A(\mathbf{Q})$, se obtiene restando la contribución de la estructura cristalográfica promedio.
2. Filtrado de Fourier: Se aíslan regiones específicas del espacio recíproco (por ejemplo, picos difusos, varillas o lóbulos) utilizando máscaras que preservan la fase ($W_F(\mathbf{Q})$). La amplitud enmascarada se somete entonces a una transformada de Fourier inversa para generar campos en el espacio real.
3. Separación de Modos: El marco emplea dos modos distintos para desenredar las contribuciones químicas y displacivas:
   • Modo Químico: Se construye una configuración de referencia donde todos los átomos se colocan en sus posiciones cristalográficas promedio, pero se retiene la distribución de las especies químicas. La amplitud difusa calculada a partir de esta configuración de "Modo Químico" aísla la dispersión que surge únicamente de las variaciones de ocupación (ordenamiento químico), suprimiendo los efectos de desplazamiento.
   • Modo de Desplazamiento: Este modo se dirige a la parte de la amplitud de dispersión sensible al desplazamiento. Debido a que la transformada inversa del campo es una convolución de desplazamientos discretos con la función de dispersión de punto de la máscara, los autores emplean un "decodificador M" lineal. Este operador lineal basado en datos mapea parches locales del campo en el espacio real filtrado hacia vectores de desplazamiento resueltos por sitio. El decodificador es entrenado con datos emparejados (campo filtrado vs. desplazamiento real conocido) para la configuración específica y luego se aplica de manera uniforme.
4. Escalabilidad: Para manejar grandes conjuntos atómicos (millones de átomos) y regiones extensas del espacio recíproco, la implementación utiliza una arquitectura Map–Reduce y Transformadas Rápidas de Fourier No Uniformes (NUFFT). Esto evita los cuellos de botella de memoria de las FFT globales basadas en rejillas y permite el procesamiento eficiente de mallas de puntos Q arbitrarias.

Resultos Clave y Validación
El artículo valida el marco MOSAIC mediante simulaciones y aplicaciones a modelos experimentales:

• Separación del Ordenamiento Químico: En un sistema simulado de $\gamma$-LiFeO$_2$ que contiene coexistencia de orden de corto alcance (SRO) y ordenamiento L1$_1$ a nanoescala, MOSAIC separó con éxito ambos efectos. Al aplicar máscaras a regiones específicas del espacio recíproco (las reflexiones del punto L frente a la superficie difusa continua), el método reconstruyó mapas distintos en el espacio real: uno que resaltaba los nanocúmulos L1$_1$ y otro que revelaba el SRO de la matriz, a pesar de que ambos contribuyen al mismo patrón de dispersión general.
• Desenredando Distorsiones: En una perovskita tipo ReO$_3$ simulada con coexistencia de rotaciones octaédricas y distorsiones de respiración (breathing), el marco separó las características difusas superpuestas. Máscaras cilíndricas aislaron el desorden rotacional (varillas difusas), mientras que máscaras esféricas aislaron el modo de respiración (picos de superred). La linealidad aditiva del método fue confirmada, ya que la suma de los campos de desplazamiento separados reconstruyó el campo de distorsión total original con un error residual insignificante.
• Aplicación Experimental: El marco se aplicó a modelos refinados por RMC de ferroeléctricos relaxores PbMg$_{1/3}$Nb$_{2/3}$O$_3$ (PMN) y PMN–PbTiO$_3$ (PMN–PT).
  • Químico: El Modo Químico reveló variaciones continuas en el grado de ordenamiento tipo sal de roca e identificó límites de fase antiparalelos, proporcionando una visión resuelta por sitio de los entornos químicos que estaba oscurecida en los modelos originales.
  • Desplazamiento: El Modo de Desplazamiento filtró el ruido no correlacionado de los campos de desplazamiento de Pb, revelando texturas polares a nanoescala y patrones de polarización giratorios consistentes con observaciones previas.

Significación y Reivindicaciones
Los autores afirman que MOSAIC proporciona un método robusto y cuantitativo para atribuir características específicas de la dispersión difusa a campos de ocupación y de desplazamiento resueltos por sitio atómico subyacente. Las contribuciones clave incluyen:

• Desacoplamiento Cuantitativo: La capacidad de separar correlaciones químicas y displacivas que producen características de dispersión superpuestas, una tarea difícil con el análisis tradicional.
• Linealidad y Aditividad: El marco garantiza una linealidad estricta, lo que significa que las contribuciones de regiones disjuntas del espacio recíproco se suman exactamente en el espacio real. Esto permite la descomposición rigurosa de patrones de dispersión complejos en motivos estructurales constituyentes.
• Escalabilidad: El uso de NUFFTs y el procesamiento Map–Reduce permite el análisis de modelos estructurales a gran escala (por ejemplo, de RMC o MD) que anteriormente eran computacionalmente prohibitivos para este tipo de filtrado detallado de características.
• Versatilidad: El método es aplicable tanto a configuraciones atómicas 3D derivadas de datos de dispersión o simulaciones, como a proyecciones estructurales 2D de microscopía electrónica de resolución atómica.

El artículo enfatiza que MOSAIC es una herramienta para analizar configuraciones conocidas para comprender los orígenes de características de dispersión específicas, en lugar de ser un método de inversión directa para resolver estructuras a partir de datos de intensidad únicamente. Ofrece una vía para vincular cuantitativamente los motivos estructurales locales con las respuestas funcionales macroscópicas en materiales complejos.