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

Technische Zusammenfassung von GenL: Ein erweiterbares Anpassungsprogramm für Laue-Oszillationen und Ganzmuster-Anpassung

Problemstellung
Hochwertige epitaktische Dünnschichten sind kritische Modellsysteme, bei denen das Vorhandensein von Laue-Oszillationen in Röntgenbeugungsmustern (XRD) als Schlüsselindikator für Kristallinität, Dickenhomogenität und Grenzflächenglätte dient. Während das qualitative Vorhandensein dieser Oszillationen gut etabliert ist, bietet die quantitative Analyse ihres Winkelabstands, ihrer Intensitätsabnahme und ihrer Asymmetrie tiefere strukturelle Einblicke hinsichtlich kohärenter Dicke, Spannungsprofile und Rauheit.

Bestehende Open-Source-Software (z. B. SUPREX, GenX, CADEM, InteractiveXRDFit) und kommerzielle Pakete konzentrieren sich häufig primär auf die Anpassung von Supergitterpeaks oder Röntgenreflektivität (empfindlich gegenüber Elektronendichteprofilen) und weniger auf die detaillierte Form von Laue-Oszillationen. Viele bestehende Werkzeuge verfügen nicht über die Fähigkeit, experimentelle Beugungsmuster unter spezifischer Berücksichtigung atomarer Ordnung, Spannungsprofile und Kristallrauheit gleichzeitig zu simulieren und anzupassen, oder sie erfordern erhebliche Neukonfigurationen, um Hochwinkel-Streudaten zu verarbeiten. Es besteht ein identifizierter Bedarf an einem flexiblen, Open-Source-Programm, das sowohl die Simulation als auch die Anpassung dieser spezifischen Merkmale ermöglicht.

Methodik
GenL ist ein MATLAB-basiertes Programm, das zur Simulation und Anpassung von Röntgenbeugungsdaten aus epitaktischen Dünnschichten entwickelt wurde. Es verfolgt einen modularen Ansatz unter Verwendung zweier theoretischer Rahmenwerke:

1. Kinematische Theorie: Der primäre Ansatz für die Anpassung, der die diffraktierte Intensität basierend auf der Summation der Streuamplituden von kohärent streuenden Gitterebenen berechnet. Sie umfasst:

   • Formfaktoren: Element-spezifische, Q-abhängige Formfaktoren einschließlich anomaler Dispersion und Absorptionskorrekturen.
   • Physikalische Korrekturen: Polarisationsfaktoren (unter Berücksichtigung der Monochromatorposition), Absorptionsfaktoren (winkelabhängige Dämpfung) und Lorentzfaktoren (instrumentelle Auflösung und Mosaizität).
   • Strukturparameter: Gitterebenenabstände, Anzahl der kohärent streuenden Ebenen und Debye-Waller-Faktoren.
   • Defekte und Grenzflächen: Modelle für Schichtrauheit (Gaußsche Verteilung der Dicke), Spannungsprofile (exponentieller oder linearer Abfall von Grenzflächen) und Substratpeaks (modelliert als Lorentz-Funktionen oder durch kohärente Amplitudenaddition).
   • Untergrund: Lineare oder polynomiale Untergrundmodellierung zur Berücksichtigung von thermischer diffuser Streuung, Fluoreszenz und Detektorrauschen.

2. Dynamische Theorie: Eine exakte Berechnung basierend auf der Formulierung von Holý und Fewster unter Verwendung des Parratt-Algorithmus für die Matrixpropagation komplexer Elektronendichte. Im Gegensatz zu Schichtmodellen, die in Reflektivitätscodes verwendet werden, unterteilt GenL Einheitszellen in zahlreiche Schichten (typischerweise 100), um die Dichte im Realraum zu berechnen, was die Konzeptionalisierung von Spannung und Rauheit ermöglicht. Dieser Ansatz berücksichtigt Brechung, Darwin-Verbreiterung und Totalreflexion, Effekte, die in der kinematischen Näherung fehlen.

Anpassungsalgorithmus
Das Programm verwendet einen Differential-Evolutions-Algorithmus innerhalb eines genetischen Algorithmus-Rahmens, um die Gütefunktion zwischen experimentellen Daten und Simulation zu minimieren. Der Anpassungsprozess justiert Parameter wie Gitterebenenabstand, kohärente Dicke, Spannungsparameter, Rauheit und Untergrundkoeffizienten. Der Algorithmus iteriert durch eine Population von Parametervektoren, bis ein Abbruchkriterium (maximale Iterationen) erreicht ist.

Hauptbeiträge und Ergebnisse

• Dual-Mode-Betrieb: GenL bietet sowohl eine grafische Benutzeroberfläche (GUI) zur Anpassung einzelner Schichten unter Verwendung der kinematischen Näherung als auch eine Kommandozeilenversion, die beliebige Strukturen (einschließlich Mehrschichten und Supergitter) unter Verwendung kinematischer oder dynamischer Ansätze verarbeiten kann.
• Spannungs- und Rauheitsmodellierung: Das Programm modelliert explizit Spannungsprofile senkrecht zur Ebene (exponentiell oder linear) und Schichtrauheit, was die Analyse asymmetrischer Laue-Oszillationen ermöglicht, die häufig mit Spannungsgradienten oder Grenzflächendefekten verbunden sind.
• Validierung und Vergleich:
  • GaAs-Blockkristall: Dynamische Berechnungen für GaAs zeigten eine hervorragende Übereinstimmung mit GenX (Parratt-Formalismus) im Reflektivitätsbereich und reproduzierten korrekt die Darwin-Breiten, was die dynamische Implementierung validiert.
  • Vanadium auf MgO: Eine Anpassung eines 105 Å dicken V-Films demonstrierte die Fähigkeit von GenL, Zugspannung senkrecht zur Ebene und Rauheit zu berücksichtigen. Die resultierende kohärente Dicke (86,3 Å) war mit früheren CADEM-Simulationen (84 Å) vergleichbar, wobei GenL ein detaillierteres Spannungsprofil lieferte.
  • Wolfram auf Al₂O₃: Anpassungen der W(110)- und W(220)-Peaks demonstrierten die Fähigkeit des Programms, überlappende Substratpeaks zu handhaben und die Grenzflächengrauheit (4 Å) abzuschätzen, was mit GenX-Reflektivitätsanpassungen konsistent ist.
  • Eisen auf MgAl₂O₄: Die Analyse von Fe-Filmen ergab, dass ein einfaches Spannungsprofil die Asymmetrie der Oszillationen nicht vollständig erfassen konnte. Ein Zweischichtenmodell (bct Fe/bcc Fe) reproduzierte die Asymmetrie erfolgreich und unterstrich die Nützlichkeit des Programms bei der Untersuchung von Wachstumsmoden im Dünnschichtlimit.
  • Supergitter: Obwohl die Anpassung von Supergitterdaten in der GUI noch nicht implementiert ist, simulierte die Kommandozeilenversion erfolgreich Beugungsmuster für Fe/V-Supergitter und demonstrierte die Erweiterbarkeit auf komplexe Mehrschichtstapel.

Bedeutung
Die Autoren präsentieren GenL als vielseitiges, Open-Source-Werkzeug (unter der GNU General Public License veröffentlicht), das eine Lücke in der Analyse epitaktischer Dünnschichten schließt. Durch die Kombination der Fähigkeit, atomare Ordnung zu simulieren, mit robusten Anpassungsfähigkeiten für Laue-Oszillationen ergänzt GenL bestehende Reflektivitätswerkzeuge wie GenX. Das Programm ermöglicht ein detaillierteres strukturelles Bild hochkristalliner Filme durch die Extraktion von Parametern wie atomaren Gitterebenenabständen, Anzahlen kohärent streuender Ebenen, Spannungsprofilen und Kristallrauheit. Die Autoren schlagen vor, dass die Fähigkeit, Grenzflächen- und Strukturparameter (wie atomare Stufen und Supergitter) einzubeziehen und anzupassen, als nützliches Asset für die Forschungscommunity dienen wird und eine detaillierte sowie quantitative Analyse von Beugungsmustern fördert.