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

GenL 技术摘要：一种用于劳厄振荡和全图谱拟合的可扩展拟合程序

问题陈述
高质量的外延薄膜是关键模型系统，其中 X 射线衍射（XRD）图谱中劳厄振荡的存在是结晶度、厚度均匀性和界面平整度的关键指标。虽然这些振荡的定性存在已得到充分证实，但对其角间距、强度衰减和不对称性的定量分析，能提供关于相干厚度、应变分布和粗糙度的更深层结构洞察。

现有的开源软件（如 SUPREX、GenX、CADEM、InteractiveXRDFit）和商业软件包通常主要侧重于拟合超晶格峰或 X 射线反射率（对电子密度分布敏感），而非劳厄振荡的详细形状。许多现有工具缺乏同时模拟和拟合实验衍射图谱的能力，特别是无法针对原子有序性、应变分布和晶体粗糙度进行专门处理，或者需要大量重新配置才能处理大角度散射数据。因此，亟需一种灵活、开源的程序，既能模拟又能拟合这些特定特征。

方法论
GenL 是一款基于 MATLAB 的程序，旨在模拟和拟合外延薄膜的 X 射线衍射数据。它采用模块化方法，利用两种理论框架：

1. 运动学理论：拟合的主要方法，基于相干散射晶面的散射振幅总和计算衍射强度。其包含：

   • 形状因子：元素特定且依赖于 Q 的形状因子，包括反常色散和吸收修正。
   • 物理修正：偏振因子（考虑单色器位置）、吸收因子（角度相关的衰减）和洛伦兹因子（仪器分辨率和镶嵌性）。
   • 结构参数：晶面间距、相干散射晶面数量以及德拜 - 沃勒因子。
   • 缺陷与界面：层粗糙度模型（厚度的高斯分布）、应变分布模型（从界面呈指数或线性衰减）以及衬底峰模型（通过洛伦兹函数或相干振幅叠加进行建模）。
   • 背景：线性或多项式背景建模，以考虑热漫散射、荧光和探测器噪声。

2. 动力学理论：基于 Holý 和 Fewster 的公式，利用 Parratt 算法进行复数电子密度的矩阵传播的精确计算。与反射率代码中使用的层状模型不同，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 超晶格的衍射图谱，证明其可扩展至复杂多层堆叠。

意义
作者将 GenL 呈现为一种多功能的开源工具（在 GNU 通用公共许可证下发布），填补了外延薄膜分析领域的空白。通过结合模拟原子有序性的能力与针对劳厄振荡的稳健拟合功能，GenL 补充了现有的反射率工具（如 GenX）。该程序通过提取原子晶面间距、相干散射晶面数量、应变分布和晶体粗糙度等参数，使高度结晶薄膜的结构图像更加详细。作者指出，包含并拟合界面和结构参数（如原子台阶和超晶格）的能力将成为研究界的一项宝贵资产，鼓励对衍射图谱进行详细和定量的分析。