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

Résumé technique de GenL : Un programme d'ajustement extensible pour les oscillations de Laue et l'ajustement de motifs complets

Énoncé du problème
Les films minces épitaxiés de haute qualité constituent des systèmes modèles essentiels, où la présence d'oscillations de Laue dans les diagrammes de diffraction des rayons X (DRX) sert d'indicateur clé de la cristallinité, de l'uniformité de l'épaisseur et de la régularité des interfaces. Bien que la présence qualitative de ces oscillations soit bien établie, l'analyse quantitative de leur espacement angulaire, de leur décroissance d'intensité et de leur asymétrie offre des insights structurels plus profonds concernant l'épaisseur cohérente, les profils de contrainte et la rugosité.

Les logiciels open source existants (par exemple, SUPREX, GenX, CADEM, InteractiveXRDFit) et les packages commerciaux se concentrent souvent principalement sur l'ajustement des pics de super-réseau ou de la réflexion des rayons X (sensible aux profils de densité électronique) plutôt que sur la forme détaillée des oscillations de Laue. De nombreux outils existants manquent de la capacité de simuler et d'ajuster simultanément les diagrammes de diffraction expérimentaux avec une attention particulière à l'ordre atomique, aux profils de contrainte et à la rugosité cristalline, ou ils nécessitent une reconfiguration significative pour traiter les données de diffusion à grand angle. Il existe un besoin identifié d'un programme flexible et open source capable de simuler et d'ajuster ces caractéristiques spécifiques.

Méthodologie
GenL est un programme basé sur MATLAB conçu pour simuler et ajuster les données de diffraction des rayons X provenant de films minces épitaxiés. Il adopte une approche modulaire utilisant deux cadres théoriques :

1. Théorie cinématique : L'approche principale pour l'ajustement, calculant l'intensité diffractée basée sur la sommation des amplitudes de diffusion provenant des plans réticulaires diffusant de manière cohérente. Elle intègre :

   • Facteurs de forme : Facteurs de forme dépendants de Q et spécifiques à l'élément, incluant les corrections de dispersion anomale et d'absorption.
   • Corrections physiques : Facteurs de polarisation (tenant compte de la position du monochromateur), facteurs d'absorption (atténuation dépendante de l'angle) et facteurs de Lorentz (résolution instrumentale et mosaïcité).
   • Paramètres structuraux : Espacements interréticulaires, nombre de plans diffusant de manière cohérente et facteurs de Debye-Waller.
   • Défauts et interfaces : Modèles de rugosité de couche (distribution gaussienne de l'épaisseur), profils de contrainte (décroissance exponentielle ou linéaire depuis les interfaces) et pics de substrat (modélisés par des fonctions de Lorentz ou via l'addition d'amplitudes cohérentes).
   • Fond : Modélisation du fond linéaire ou polynomial pour tenir compte de la diffusion thermique diffuse, de la fluorescence et du bruit du détecteur.

2. Théorie dynamique : Un calcul exact basé sur la formulation de Holý et Fewster utilisant l'algorithme de Parratt pour la propagation matricielle de la densité électronique complexe. Contrairement aux modèles de plaques utilisés dans les codes de réflexion, GenL divise les mailles élémentaires en de nombreuses plaques (généralement 100) pour calculer la densité dans l'espace réel, permettant la conceptualisation de la contrainte et de la rugosité. Cette approche prend en compte la réfraction, l'élargissement de Darwin et la réflexion totale, des effets absents dans l'approximation cinématique.

Algorithme d'ajustement
Le programme utilise un algorithme d'évolution différentielle dans le cadre d'un algorithme génétique pour minimiser la figure de mérite entre les données expérimentales et la simulation. Le processus d'ajustement ajuste des paramètres tels que l'espacement interréticulaire, l'épaisseur cohérente, les paramètres de contrainte, la rugosité et les coefficients de fond. L'algorithme itère à travers une population de vecteurs de paramètres jusqu'à ce qu'un critère d'arrêt (nombre maximum d'itérations) soit atteint.

Contributions et résultats clés

• Fonctionnement en double mode : GenL offre à la fois une interface graphique (GUI) pour l'ajustement de couches simples utilisant l'approximation cinématique et une version en ligne de commande capable de gérer des structures arbitraires (y compris les multicouches et les super-réseaux) utilisant soit l'approche cinématique, soit l'approche dynamique.
• Modélisation de la contrainte et de la rugosité : Le programme modélise explicitement les profils de contrainte hors plan (exponentiels ou linéaires) et la rugosité de couche, permettant l'analyse des oscillations de Laue asymétriques souvent associées à des gradients de contrainte ou à des défauts d'interface.
• Validation et comparaison :
  • Cristal massif de GaAs : Les calculs dynamiques pour le GaAs ont montré un excellent accord avec GenX (formalisme de Parratt) dans la région de réflexion et ont correctement reproduit les largeurs de Darwin, validant ainsi l'implémentation dynamique.
  • Vanadium sur MgO : Un ajustement d'un film de V de 105 Å a démontré la capacité de GenL à prendre en compte la contrainte hors plan en traction et la rugosité. L'épaisseur cohérente résultante (86,3 Å) était comparable aux simulations CADEM précédentes (84 Å), GenL fournissant un profil de contrainte plus détaillé.
  • Tungstène sur Al₂O₃ : Les ajustements des pics W (110) et W (220) ont démontré la capacité du programme à gérer les pics de substrat se chevauchant et à estimer la rugosité d'interface (4 Å), cohérent avec les ajustements de réflexion de GenX.
  • Fer sur MgAl₂O₄ : L'analyse des films de Fe a révélé qu'un profil de contrainte simple ne pouvait pas capturer entièrement l'asymétrie des oscillations. Un modèle bicouche (Fe bct/Fe bcc) a reproduit avec succès l'asymétrie, mettant en évidence l'utilité du programme pour investiguer les modes de croissance à la limite du film mince.
  • Super-réseaux : Bien que l'ajustement des données de super-réseaux ne soit pas encore implémenté dans l'interface graphique, la version en ligne de commande a simulé avec succès les diagrammes de diffraction pour les super-réseaux Fe/V, démontrant l'extensibilité à des empilements multicouches complexes.

Signification
Les auteurs présentent GenL comme un outil open source polyvalent (publié sous la Licence publique générale GNU) qui comble une lacune dans l'analyse des films minces épitaxiés. En combinant la capacité de simuler l'ordre atomique avec des capacités d'ajustement robustes pour les oscillations de Laue, GenL complète les outils de réflexion existants comme GenX. Le programme permet une image structurelle plus détaillée des films hautement cristallins en extrayant des paramètres tels que les espacements interréticulaires atomiques, les nombres de plans de diffusion cohérente, les profils de contrainte et la rugosité cristalline. Les auteurs suggèrent que la capacité d'inclure et d'ajuster les paramètres d'interface et structuraux (tels que les marches atomiques et les super-réseaux) constituera un atout utile pour la communauté de recherche, encourageant une analyse détaillée et quantitative des diagrammes de diffraction.