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

GenL の技術的概要：ラウエ振動および全パターンフィッティング用の拡張可能フィッティングプログラム

問題提起
高品質なエピタキシャル薄膜は重要なモデル系であり、X 線回折（XRD）パターンにおけるラウエ振動の存在は、結晶性、厚さの均一性、および界面の滑らかさの重要な指標となる。これらの振動の定性的な存在はよく確立されているが、その角間隔、強度減衰、および非対称性の定量的分析は、コヒーレント厚さ、ひずみプロファイル、および粗さに関するより深い構造的洞察を提供する。

既存のオープンソースソフトウェア（SUPREX、GenX、CADEM、InteractiveXRDFit など）および商用パッケージは、主に超格子ピークまたは X 線反射率（電子密度プロファイルに敏感）のフィッティングに焦点を当てており、ラウエ振動の詳細な形状には重点が置かれていないことが多い。多くの既存ツールは、原子配列、ひずみプロファイル、および結晶粗さに特段の注意を払って、実験的回折パターンを同時にシミュレーションおよびフィッティングする能力を欠いており、または高角散乱データを処理するために大幅な再構成を必要とする。これらの特定の機能をシミュレーションおよびフィッティングできる柔軟なオープンソースプログラムへの必要性が確認されている。

手法
GenL は、エピタキシャル薄膜からの X 線回折データをシミュレーションおよびフィッティングするために設計された MATLAB ベースのプログラムである。これは、2 つの理論的枠組みを利用するモジュール式のアプローチを採用している。

1. 運動学理論：フィッティングの主要なアプローチであり、コヒーレントに散乱する格子面からの散乱振幅の総和に基づいて回折強度を計算する。これには以下が含まれる：

   • 原子散乱因子：異常分散および吸収補正を含む、元素固有の Q 依存性原子散乱因子。
   • 物理的補正：偏光因子（モノクロメータの位置を考慮）、吸収因子（角度依存の減衰）、およびローレンツ因子（装置分解能およびモザイク性）。
   • 構造パラメータ：格子面間隔、コヒーレントに散乱する面の数、およびデバイ - ワッラー因子。
   • 欠陥および界面：層粗さ（厚さのガウス分布）、ひずみプロファイル（界面からの指数関数的または線形減衰）、および基板ピーク（ローレンツ関数によるモデル化またはコヒーレント振幅の加算によるモデル化）のモデル。
   • 背景：熱拡散散乱、蛍光、および検出器ノイズを説明するための線形または多項式背景モデル。

2. 動力学理論：Parratt のアルゴリズムを用いた複素電子密度の行列伝播に基づく Holý および Fewster の定式化に基づく厳密な計算である。反射率コードで使用されるスラブモデルとは異なり、GenL は単位格子を多数のスラブ（通常 100）に分割して実空間で密度を計算し、ひずみおよび粗さの概念化を可能にする。このアプローチは、運動学近似には存在しない屈折、ダーウィン幅広がり、および全反射を考慮する。

フィッティングアルゴリズム
本プログラムは、実験データとシミュレーション間の適合度指標を最小化するために、遺伝的アルゴリズムの枠組み内で微分進化アルゴリズムを利用する。フィッティングプロセスは、格子面間隔、コヒーレント厚さ、ひずみパラメータ、粗さ、および背景係数などのパラメータを調整する。アルゴリズムは、停止基準（最大反復回数）が満たされるまで、パラメータベクトルの集団を反復する。

主な貢献と結果

• デュアルモード動作：GenL は、運動学近似を用いて単一層をフィッティングするためのグラフィカルユーザーインターフェース（GUI）と、コマンドライン版の両方を提供する。コマンドライン版は、運動学的または動力的アプローチのいずれかを使用して、任意の構造（多層および超格子を含む）を処理可能である。
• ひずみおよび粗さのモデル化：本プログラムは、面外ひずみプロファイル（指数関数的または線形）および層粗さを明示的にモデル化し、ひずみ勾配または界面欠陥にしばしば関連する非対称なラウエ振動の分析を可能にする。
• 検証と比較：
  • GaAs バルク結晶：GaAs に対する動力学計算は、反射率領域において GenX（Parratt 形式）と優れた一致を示し、ダーウィン幅を正しく再現することで、動力学実装を検証した。
  • MgO 上のバナジウム：105 Å の V 薄膜のフィッティングは、GenL が引張面外ひずみおよび粗さを考慮する能力を実証した。得られたコヒーレント厚さ（86.3 Å）は、以前の CADEM シミュレーション（84 Å）と同等であり、GenL はより詳細なひずみプロファイルを提供した。
  • Al₂O₃上のタングステン：W(110) および W(220) ピークのフィッティングは、本プログラムが重なり合う基板ピークを処理し、界面粗さ（4 Å）を推定する能力を実証した。これは GenX の反射率フィッティングと一致する。
  • MgAl₂O₄上の鉄：Fe 薄膜の分析により、単純なひずみプロファイルでは振動の非対称性を完全に捉えられないことが明らかになった。二層モデル（bct Fe/bcc Fe）は非対称性を成功裡に再現し、薄膜限界における成長モードの調査における本プログラムの有用性を浮き彫りにした。
  • 超格子：超格子データのフィッティングは GUI ではまだ実装されていないが、コマンドライン版は Fe/V 超格子の回折パターンを成功裡にシミュレーションし、複雑な多層スタックへの拡張可能性を実証した。

意義
著者らは、エピタキシャル薄膜の分析におけるギャップを埋める汎用性の高いオープンソースツール（GNU 一般公衆ライセンスの下でリリース）として GenL を提示する。原子配列のシミュレーション能力とラウエ振動に対する堅牢なフィッティング能力を組み合わせることで、GenL は GenX などの既存の反射率ツールを補完する。本プログラムは、原子格子面間隔、コヒーレント散乱面数、ひずみプロファイル、および結晶粗さなどのパラメータを抽出することで、高結晶性薄膜のより詳細な構造的像を可能にする。著者らは、界面および構造パラメータ（原子ステップや超格子など）を含めてフィッティングする能力が研究コミュニティにとって有用な資産となり、回折パターンの詳細かつ定量的な分析を促進すると示唆している。