基于X射线反射率（XRR）的薄膜厚度测定方法比较：精度、计算成本与实用性的系统评估

《Journal of Nanotechnology》：Comparison of Thin Film Thickness Determination Methods Based on X-Ray Reflectivity

为解决XRR测定薄膜厚度时不同分析方法的准确性与适用性不清的问题，研究者系统比较了峰位分离法、FFT法、修正Bragg方程法和曲线拟合法。四种方法所得厚度值在1.34%内一致，为不同场景下的方法选择提供了实用指导。

　　

在纳米技术与材料科学的微观世界里，薄膜（厚度通常在纳米到几微米之间）扮演着举足轻重的角色。从生物医学涂层、光学干涉镜，到半导体异质结构、高密度集成电路乃至航空航天与生物材料用的金属陶瓷薄膜，它们的性能往往极度依赖于那层肉眼不可见的纳米级厚度。然而，如何在非破坏的前提下，精准地测量这层薄膜的厚度、密度以及界面的粗糙度，一直是实验物理学家和材料科学家面临的现实难题。

X射线反射率（X-ray reflectivity, XRR）技术应运而生，成为这一领域的非破坏性检测利器。当X射线以极小的掠入射角照射在覆盖有薄膜的基板上时，不仅在薄膜上表面发生反射，还会在薄膜与基板的交界面反射，这些反射波相互干涉，便在反射率-角度曲线上形成了周期性的振荡，即著名的Kiessig fringes。理论上，通过分析这些条纹的周期可推算厚度，通过观察振荡的衰减可评估粗糙度。但在实际操作层面，科研人员面临着多种数据分析路径的选择：是直接测量峰位间距？是对信号进行快速傅里叶变换（FFT）？是利用修正的Bragg方程进行线性拟合？还是通过复杂的曲线拟合（Curve fitting）来模拟整个物理模型？

这些方法各有拥趸，但也各有局限。简单的峰位法受低角度折射影响大；FFT法需要处理非均匀采样和频谱泄露；修正Bragg方程面临干涉级次分配的难题；而曲线拟合虽然信息全面，但计算耗时且对初值敏感。此前缺乏系统的横向对比，研究人员往往不清楚在特定条件下该信赖哪种方法，也难以权衡精度与效率。为此，来自都柏林圣三一学院（Trinity College Dublin）CRANN纳米结构与器件研究中心的团队，以Al₂O₃基板上的Fe/Fe₂O₃薄膜系统为样本，对上述四种主流XRR数据分析方法进行了一次彻头彻尾的“大考”，相关成果发表在《Journal of Nanotechnology》。

为开展此项研究，作者主要采用了以下关键技术方法：使用配备Cu Kα1光源（λ = 0.1540562 nm）的Bruker D8 Discover高精度X射线衍射仪，在0.2001°至6.0001°的2θ范围内以0.02°步长采集XRR原始数据；样本为在Al₂O₃基板上通过分子束外延（MBE）生长的Fe薄膜（表层暴露于空气形成天然Fe₂O₃氧化层）；数据处理与四种厚度测定算法的实现均基于MATLAB平台；其中曲线拟合方法具体采用了基于差分进化（Differential evolution）算法的GenX软件进行多层模型优化。

3.1. 峰位分离法（Peak Separation Method）

该方法基于简化的运动学近似，利用相邻干涉条纹间的角距离（Δ2θ）来计算厚度。研究发现，由于X射线折射率略小于1，在近临界角（θ_c）的低角度区域，折射效应会扭曲峰位，导致早期条纹间距不均。因此，研究者选择较大角度（高2θ）的峰值进行计算。通过对原始数据（0.9°至2.24°区间）进行Savitzky-Golay平滑、归一化、去直流（DC）分量后，提取峰值位置。结果显示，若包含低角度峰值，计算出的厚度误差较大（如88.27 nm或73.56 nm），而仅使用高角度最后四个峰值（对应0.12°至0.14°的Δ2θ）时，平均厚度值为67.9 nm。其不确定度分析得出结果为67.9 ± 3.0 nm。这表明，在忽略复杂折射修正的情况下，峰位分离法若合理选取高角度振荡区，可快速提供可靠的等效厚度估计。

3.2. 傅里叶变换法（Fourier Transform Method, FFT）

FFT法的核心在于将XRR振荡信号转换到频率域，从而直接读取对应厚度的频率分量。文章详细梳理了该方法的物理推导：在远场近似下（θ >> θ_c），反射率表达式可简化为含余弦项的谐波函数，其“准波矢量”x = 2√(sin²θ - sin²θ_c)/λ 作为变换轴（类比时间轴），变换后的峰值位置则直接对应厚度（nm）。研究中特别强调，此前文献中“准波矢量”定义常遗漏因子2，会导致厚度计算出入两倍，本文明确了应包含该因子。数据预处理包括：提取振荡区（0.45°< θ <1.12°）、归一化、去DC分量、通过线性插值将非均匀x轴转换为均匀采样（以满足FFT要求）、加汉宁（Hamming）窗以减少频谱泄露，以及补零（Zero padding）至2的幂次（如67点补至128点）。最终FFT的功率谱密度（PSD）显示出尖锐峰值，对应厚度为67.52 nm（675.2 ?），不确定度为67.52 ± 6.0 nm。该方法虽预处理繁琐，但能有效分离各层对应的频率成分。

3.3. 基于全反射条件修正的Bragg方程法（Bragg Equation Modified by the Total Reflection Condition Method）

此方法将干涉极值（极大值和极小值）位置与厚度通过修正的Bragg方程联系起来：θ_i²= (n_i+ Δn)λ / (2t)²+ θ_c²。通过绘制θ_i²对 (n_i+ Δn) 的线性关系进行最小二乘拟合，斜率用于求厚度，截距用于求临界角θ_c。关键挑战在于干涉级次n_i的起始赋值，需迭代验证使截距符合材料电子密度理论值。本研究同时利用了11个极大值和11个极小值的位点。拟合得到斜率为1.34398×10^-6，截距为3.5662×10^-5，算出厚度t = λ / (2√斜率) = 66.45 nm，临界角θ_c= 0.342°。误差传播分析基于斜率的标准误，得出结果为66.45 ± 0.09 nm。该方法通过线性拟合有效修正了折射影响，能同时给出厚度与临界角。

3.4. 曲线拟合法（Curve Fitting Method）

作为最严密但也最耗时的手段，研究使用了GenX软件，基于Parratt形式化理论和差分进化算法进行全局优化。与前述视系统为单一等效层不同，拟合模型允许设置多层结构。针对Fe膜自然氧化问题，模型设为两层：底层Fe（密度可调）和表层Fe₂O₃氧化层（密度可调），各层均设厚度与粗糙度参数。拟合结果精准匹配了0.2°至2.24°（2θ）的全范围数据，揭示出：Fe层厚度为50 nm，Fe₂O₃层厚度为17 nm，总计67 nm；同时给出了各层密度与粗糙度。单个层的不确定度通过χ²景观分析得出（Fe: 50_-0.4^+0.5nm；Fe₂O₃: 17_-0.1^+0.2nm），总厚度误差按高斯传播计算得67 ± 0.43 nm。此方法虽计算重，但能解耦亚层结构，是唯一能明确分辨出17 nm氧化层的方法。

结论与讨论

综合全文四种方法的定量结果——峰位分离法（67.9 ± 3.0 nm）、FFT法（67.52 ± 6.0 nm）、修正Bragg方程法（66.45 ± 0.09 nm）以及曲线拟合法（67 ± 0.43 nm）——所有值均在1.34%的偏差范围内高度一致，验证了XRR技术在薄膜厚度测定上的可靠性。

研究进一步明确了各方法的定位：峰位分离法最为快捷，仅需有限个高角干涉峰即可在预筛选或快速检测中提供有效估计；FFT法虽需繁复的预处理（平滑、归一化、插值、加窗、补零），但能解析出各层对应的频率分量，且结果（67.52 nm）极接近曲线拟合的真值，适合需分层信息但又想避开复杂建模的场景；修正Bragg方程法通过线性拟合极值位置，不仅能以约0.8%的偏差获知总厚度，还可一并确定临界角θ_c，是平衡速度与精度的有力工具；曲线拟合法（如GenX）虽最耗时，却是处理复杂多层膜、获取厚度/密度/粗糙度甚至探测表面氧化层（如本例中的17 nm Fe₂O₃）不可或缺的金标准。此外，将前三种方法得出的厚度作为初值代入曲线拟合，能显著降低优化算法的复杂度。

总之，这项系统比较不仅为XRR薄膜表征领域提供了一份清晰的“方法选型指南”，也再次印证了无论采用何种简化解析或复杂模拟，只要处理得当，XRR都能交出一份高精度的纳米级厚度答卷。对于纳米制造与薄膜工业应用而言，理解这些方法的trade-off（取舍），意味着能在研发效率与数据深度之间做出更明智的把控。