GenL: An extensible fitting program for Laue oscillations and whole pattern fitting

Resumen Técnico de GenL: Un Programa de Ajuste Extensible para Oscilaciones de Laue y Ajuste de Patrones Completos

Planteamiento del Problema
Las películas delgadas epitaxiales de alta calidad son sistemas modelo críticos, donde la presencia de oscilaciones de Laue en los patrones de difracción de rayos X (XRD) sirve como un indicador clave de cristalinidad, uniformidad de espesor y suavidad de la interfaz. Si bien la presencia cualitativa de estas oscilaciones está bien establecida, el análisis cuantitativo de su espaciado angular, decaimiento de intensidad y asimetría ofrece perspectivas estructurales más profundas sobre el espesor coherente, los perfiles de deformación y la rugosidad.

El software de código abierto existente (por ejemplo, SUPREX, GenX, CADEM, InteractiveXRDFit) y los paquetes comerciales a menudo se centran principalmente en el ajuste de picos de superred o de reflectividad de rayos X (sensibles a los perfiles de densidad electrónica) en lugar de la forma detallada de las oscilaciones de Laue. Muchas herramientas existentes carecen de la capacidad de simular y ajustar simultáneamente patrones de difracción experimentales con una atención específica al ordenamiento atómico, los perfiles de deformación y la rugosidad cristalina, o requieren una reconfiguración significativa para manejar datos de dispersión a altos ángulos. Existe una necesidad identificada de un programa flexible y de código abierto capaz de simular y ajustar estas características específicas.

Metodología
GenL es un programa basado en MATLAB diseñado para simular y ajustar datos de difracción de rayos X de películas delgadas epitaxiales. Emplea un enfoque modular que utiliza dos marcos teóricos:

1. Teoría Cinemática: El enfoque principal para el ajuste, calculando la intensidad difractada basada en la suma de las amplitudes de dispersión de los planos de la red que dispersan coherentemente. Incorpora:

   • Factores de Forma: Factores de forma dependientes del elemento y de Q, que incluyen correcciones por dispersión anómala y absorción.
   • Correcciones Físicas: Factores de polarización (que tienen en cuenta la posición del monocromador), factores de absorción (atenuación dependiente del ángulo) y factores de Lorentz (resolución instrumental y mosaico).
   • Parámetros Estructurales: Espaciamientos interplanares, número de planos que dispersan coherentemente y factores de Debye-Waller.
   • Defectos e Interfases: Modelos para la rugosidad de la capa (distribución gaussiana del espesor), perfiles de deformación (decaimiento exponencial o lineal desde las interfases) y picos del sustrato (modelados como funciones de Lorentz o mediante adición de amplitudes coherentes).
   • Fondo: Modelado de fondo lineal o polinomial para tener en cuenta la dispersión térmica difusa, la fluorescencia y el ruido del detector.

2. Teoría Dinámica: Un cálculo exacto basado en la formulación de Holý y Fewster que utiliza el algoritmo de Parratt para la propagación matricial de la densidad electrónica compleja. A diferencia de los modelos de capas utilizadas en los códigos de reflectividad, GenL divide las celdas unitarias en numerosas capas (típicamente 100) para calcular la densidad en el espacio real, permitiendo la conceptualización de la deformación y la rugosidad. Este enfoque tiene en cuenta la refracción, el ensanchamiento de Darwin y la reflexión total, efectos ausentes en la aproximación cinemática.

Algoritmo de Ajuste
El programa utiliza un algoritmo de evolución diferencial dentro de un marco de algoritmo genético para minimizar la figura de mérito entre los datos experimentales y la simulación. El proceso de ajuste optimiza parámetros como el espaciamiento interplanar, el espesor coherente, los parámetros de deformación, la rugosidad y los coeficientes de fondo. El algoritmo itera a través de una población de vectores de parámetros hasta que se cumple un criterio de parada (máximo de iteraciones).

Contribuciones y Resultados Clave

• Operación en Doble Modo: GenL ofrece tanto una Interfaz Gráfica de Usuario (GUI) para ajustar capas individuales utilizando la aproximación cinemática como una versión de línea de comandos capaz de manejar estructuras arbitrarias (incluyendo multicapas y superredes) utilizando enfoques cinemáticos o dinámicos.
• Modelado de Deformación y Rugosidad: El programa modela explícitamente perfiles de deformación fuera del plano (exponenciales o lineales) y rugosidad de la capa, permitiendo el análisis de oscilaciones de Laue asimétricas a menudo asociadas con gradientes de deformación o defectos de interfaz.
• Validación y Comparación:
  • Cristal Masivo de GaAs: Los cálculos dinámicos para GaAs mostraron un excelente acuerdo con GenX (formalismo de Parratt) en la región de reflectividad y reprodujeron correctamente los anchos de Darwin, validando la implementación dinámica.
  • Vanadio sobre MgO: Un ajuste de una película de V de 105 Å demostró la capacidad de GenL para tener en cuenta la deformación de tracción fuera del plano y la rugosidad. El espesor coherente resultante (86.3 Å) fue comparable a las simulaciones anteriores de CADEM (84 Å), proporcionando GenL un perfil de deformación más detallado.
  • Tungsteno sobre Al₂O₃: Los ajustes de los picos W (110) y W (220) demostraron la capacidad del programa para manejar picos de sustrato superpuestos y estimar la rugosidad de la interfaz (4 Å), consistente con los ajustes de reflectividad de GenX.
  • Hierro sobre MgAl₂O₄: El análisis de películas de Fe reveló que un perfil de deformación simple no podía capturar completamente la asimetría de las oscilaciones. Un modelo de bicapa (Fe bct/Fe bcc) reprodujo con éxito la asimetría, destacando la utilidad del programa para investigar modos de crecimiento en el límite de la película delgada.
  • Superredes: Si bien el ajuste de datos de superredes aún no está implementado en la GUI, la versión de línea de comandos simuló con éxito patrones de difracción para superredes de Fe/V, demostrando la extensibilidad a pilas multicapa complejas.

Significado
Los autores presentan GenL como una herramienta versátil y de código abierto (publicada bajo la Licencia Pública General GNU) que llena un vacío en el análisis de películas delgadas epitaxiales. Al combinar la capacidad de simular el ordenamiento atómico con capacidades de ajuste robustas para las oscilaciones de Laue, GenL complementa las herramientas de reflectividad existentes como GenX. El programa permite una imagen estructural más detallada de películas altamente cristalinas al extraer parámetros como espaciamientos interplanares atómicos, conteos de planos de dispersión coherente, perfiles de deformación y rugosidad cristalina. Los autores sugieren que la capacidad de incluir y ajustar parámetros de interfaz y estructurales (como escalones atómicos y superredes) servirá como un activo útil para la comunidad de investigación, fomentando un análisis detallado y cuantitativo de los patrones de difracción.