椭圆偏振光谱仪测量薄膜厚度与光学常数的原理与方法
椭圆偏振光谱仪（Spectroscopic Ellipsometry, SE）通过测量偏振光在薄膜表面反射后偏振态的变化，采用模型拟合逆向求解薄膜的厚度 d、折射率 n 和消光系数 k。作为一种非接触、非破坏性、亚纳米级精度的光学测量技术，椭偏仪广泛应用于半导体、光伏、光学镀膜等领域的薄膜表征。本文面向薄膜工艺工程师和材料表征人员，涵盖从基础物理原理到实际建模操作的完整技术链条，所有讨论均以光谱椭偏仪的实际使用场景为前提。
1 基本原理
1.1 椭圆偏振光的形成
入射光束可分解为两个正交分量：
- p 分量：电场振动方向平行于入射面
- s 分量：电场振动方向垂直于入射面（源自德文 senkrecht）
当线偏振光以一定角度照射到薄膜表面时，p 分量和 s 分量在界面 1（空气/薄膜）和界面 2（薄膜/衬底）发生多次反射和折射。各反射光之间产生多光束干涉，其干涉结果反映薄膜的光学特性。
菲涅尔反射系数 rp 和 rs 分别描述 p 分量和 s 分量在界面处的振幅和相位变化。对于光从介质 i 入射到介质 j 的单一界面，菲涅尔系数为：
rp(ij) = (Nj cos φi − Ni cos φj) / (Nj cos φi + Ni cos φj)
rs(ij) = (Ni cos φi − Nj cos φj) / (Ni cos φi + Nj cos φj)
其中 Ni 和 Nj 为两种介质的复折射率（N = n + ik），φi 和 φj 分别为入射角和折射角，满足 Snell 定律 Ni sin φi = Nj sin φj。
对于单层薄膜结构，光在薄膜内部经历无限多次反射，总反射系数由各次反射的叠加构成：
R = (r01 + r12 e^(−i2β)) / (1 + r01 r12 e^(−i2β))
其中 r01 为空气/薄膜界面的菲涅尔系数，r12 为薄膜/衬底界面的菲涅尔系数，β = 2π(d/λ) N1 cos φ1 为薄膜内部的相位厚度。β 中的 d 和 N1 正是待求的膜厚和复折射率。
由于 p 分量和 s 分量的菲涅尔反射系数不同——包括振幅和相位变化的差异——反射光中两个分量的合成轨迹呈椭圆状，故称椭圆偏振光。
1.2 椭偏参数 Ψ 和 Δ
椭偏仪直接测量的量是菲涅尔反射系数比 ρ：
ρ = Rp / Rs = tanΨ · e^(iΔ)
| 参数 | 定义 | 物理意义 |
|---|---|---|
| Ψ（Psi） | tanΨ = |Rp| / |Rs| | 反射后 p 分量与 s 分量的振幅比 |
| Δ（Delta） | Δ = δp − δs | 反射后 p 分量与 s 分量的相位差 |
Ψ 的取值范围为 0° ~ 90°，Δ 的取值范围为 0° ~ 360°。当 Δ = 0° 或 Δ = 180° 时，反射光为线偏振光；当 Δ = 90° 或 Δ = 270° 且 Ψ = 45° 时，反射光为圆偏振光——这是两种极限情况。
Ψ 和 Δ 与薄膜厚度 d、折射率 n 和消光系数 k 之间存在隐函数关系：
(Ψ, Δ) = f(N1, N2, d, λ, φ₁)
其中 N1 和 N2 分别为薄膜和衬底的复折射率，λ 为入射光波长，φ₁ 为入射角。由此可知，测量得到 Ψ 和 Δ 后，便可逆向求解薄膜的光学常数和厚度。
Ψ 和 Δ 对薄膜参数变化的灵敏度并非均匀——厚度变化引起的 Δ 变化远大于 Ψ 的变化，而折射率的变化则在 Ψ 上体现更为显著。这一特性意味着对于超薄膜（< 5 nm），厚度和折射率之间存在较强的参数耦合，必须通过多角度或多波长测量实现解耦。
1.3 为何不能直接求解
椭偏方程是非线性超越方程，无法直接解析求出 n、k、d。必须采用"假设模型 + 正向计算 + 最小二乘拟合"的迭代方法：
- 建立薄膜物理结构模型（层数、各层材料类型）
- 选择色散模型（Cauchy、Sellmeier、Lorentz 等）
- 给定初始试探解（n₀, k₀, d₀）
- 利用菲涅尔公式和特征矩阵法正向计算理论 (Ψ_calc, Δ_calc)
- 与实测 (Ψ_meas, Δ_meas) 对比，最小化均方误差 MSE
- 迭代优化至收敛，输出 n、k、d 最优解
模型选择是拟合成功的关键——错误的模型会导致拟合结果与真实值严重偏离。
对于多层膜系统，正向计算需借助特征矩阵法。第 j 层薄膜的特征矩阵为：
Mj = [[cos βj, (i/qj) sin βj], [i qj sin βj, cos βj]]
其中 βj = 2π(dj/λ) Nj cos φj，qj 对 p 偏振为 Nj / cos φj、对 s 偏振为 Nj cos φj。多层膜的总特征矩阵为各层矩阵的乘积：
M_total = M1 · M2 · … · Mm
由此计算总反射系数 Rp 和 Rs，进而获得理论椭偏参数用于拟合。特征矩阵法的优势在于可以处理任意层数的薄膜结构，且每层的厚度和光学常数作为独立参数参与拟合。但对于超过 5 层的复杂膜系，自由参数过多将显著增加拟合不确定度，需通过固定已知层参数（如衬底光学常数由裸片测量获得）来约束求解。
2 仪器构造
椭偏仪由以下核心组件构成：
| 组件 | 功能 |
|---|---|
| 光源 | Xe 灯、Hg-Xe 灯或激光（HeNe），提供稳定宽谱/单色光束 |
| 起偏器 | 将自然光转换为线偏振光 |
| 补偿器（1/4 波片） | 将线偏振光转换为椭圆偏振光，控制入射偏振态 |
| 样品台 | 固定待测样品，控制入射角（45° ~ 90°） |
| 检偏器 | 分析反射后光的偏振态变化 |
| 探测器 | 光电倍增管、硅光电池或 InGaAs，记录光强 |
| 数据处理系统 | 光谱拟合与参数求解 |
光谱椭偏仪在宽波长范围（250 ~ 1700 nm）内扫描，获得 Ψ(λ) 和 Δ(λ) 光谱，比单波长椭偏仪具有更高的信息量和准确性。穆勒矩阵椭偏仪可以分析各向异性材料的完整光学响应。
按测量原理，椭偏仪可分为三种主要构型：
消光式椭偏仪：手动旋转起偏器和检偏器至探测器接收光强最小时的方位角，直接读取 Ψ 和 Δ。原理简单、无需探测器校准，但逐点测量速度慢，不适用于光谱测量。历史上是最早的椭偏仪构型。
旋转检偏器式椭偏仪：起偏器固定，检偏器以恒定角速度旋转，探测器记录光强随检偏器角度的调制波形，通过傅里叶分析提取 Ψ 和 Δ。光谱测量速度快，但 Δ 在 0° 和 180° 附近灵敏度显著下降。
旋转补偿器式椭偏仪：在起偏器和检偏器之间插入以恒定角速度旋转的补偿器，探测器信号同时包含 Ψ 和 Δ 的完整信息。可测量全范围 Δ（0° ~ 360°），且 Δ 在 0°/180° 附近无灵敏度零点，特别适合透明衬底上的超薄膜测量——后者在旋转检偏器构型中恰处于灵敏度盲区。
穆勒矩阵椭偏仪使用双旋转补偿器设计，可测量样品完整 4×4 穆勒矩阵的 16 个元素，适用于各向异性材料（单晶硅不同晶向、应力诱导双折射）、手性材料及具有退偏效应的粗糙表面。
3 常用色散模型
不同材料需要选择匹配的色散模型：
| 模型 | 适用材料 | 公式形式 |
|---|---|---|
| Cauchy | 透明材料（SiO₂、Al₂O₃、SiNₓ、MgF₂、氧化铟锡（ITO）） | n(λ) = A + B/λ² + C/λ⁴, k = 0 |
| Sellmeier | 透明 + 红外吸收材料（Ge、Si） | n²(λ) = 1 + Σ(Aᵢλ²)/(λ² − Bᵢ²) |
| Lorentz | 吸收材料（a-Si、金属） | ε = ε∞ + Σ(AᵢE₀ᵢ²)/(E₀ᵢ² − E² − iΓᵢE) |
| Tauc-Lorentz | 非晶半导体 | 组合 Tauc 带隙和 Lorentz 振子 |
| EMA | 混合/多孔材料 | 有效介质近似，三种模型可选 |
| Drude | 金属薄膜 | 自由电子气模型 |
对于光伏电池中常见的 SiNₓ:H 减反膜，通常采用 Cauchy 模型（透明波段）+ Tauc-Lorentz 模型（紫外吸收区）的组合。
有效介质近似（EMA）模型用于描述由两种或多种材料物理混合而成的薄膜——如多孔 SiO₂（SiO₂ + 空气）、粗糙表面层（薄膜材料 + 空气）或共溅射复合膜。EMA 将混合膜的宏观介电函数 ε 与各组分介电函数 εi 和体积分数 fi 关联。三种标准 EMA 公式的适用场景不同：
Maxwell-Garnett 模型：ε = εh [1 + 2f(εi − εh)/(εi + 2εh)] / [1 − f(εi − εh)/(εi + 2εh)]，适用于主相材料（介电函数 εh）中分散少量第二相（介电函数 εi）的体系，如金属纳米颗粒嵌入介质基体。
Bruggeman 模型：Σ fi (εi − ε)/(εi + 2ε) = 0，假设各组分以统计对称的方式混合，不存在明确的主相/分散相区分。适用于多孔材料、表面粗糙层的椭偏建模——实际工作中粗糙层通常设为 50% 薄膜材料 + 50% 空穴（void）的 Bruggeman EMA 层。
Linear 模型：ε = f εi + (1 − f) εj，各组分介电函数按体积分数线性叠加，适用于折射率对比度较低的混合体系。
EMA 模型中体积分数 fi 作为额外的拟合参数引入，需谨慎使用——增加一个自由参数将增大拟合不确定度，应在有物理依据时引入（如截面 SEM/TEM 确认存在孔隙或粗糙层）。
模型选择对拟合结果的影响可通过 MSE 的定量比较来评估。两个嵌套模型的 MSE 差异若在统计上显著（F-test），则接受参数更多的模型；否则选择参数更少的模型以避免过拟合。拟合过程中还需关注参数间的相关系数矩阵——两参数相关系数绝对值 > 0.9 表明二者强耦合，拟合结果可能不稳定，需通过固定其中一参数或增加多角度测量数据来实现解耦。
4 技术特点
| 特点 | 说明 |
|---|---|
| 纳米级精度 | 可探测单原子层厚度变化（< 0.1 nm），重复性 0.02 nm |
| 非接触无损 | 纯光学测量，不接触样品，无损伤 |
| 多参数同时获取 | 一次测量获得 d、n、k，以及粗糙度、组分、应力等 |
| 测量范围广 | 1 nm ~ 数十 μm，最小测量区域 ~0.3 mm |
| 快速 | 20 ms/点，适用于在线生产检测 |
| 无需参考样品 | 测量 p/s 分量比值，不受光源波动影响 |
| 可测量多层膜 | 特征矩阵法处理任意膜层数 |
椭偏测量对入射角的精度要求极高。入射角误差 0.1° 即可引入可观的膜厚拟合偏差——透明薄膜的影响尤为显著，因为透明膜厚度的信息主要来源于膜内干涉引起的 Δ 振荡周期，该周期对入射角高度敏感。在线检测系统需定期使用标准样品（热氧化 SiO₂/Si，厚度经 TEM 标定）校准入射角。环境温度波动引起的样品台机械热膨胀同样会导致入射角漂移，稳定温控是获得 0.02 nm 重复性精度的重要前提。
对于超薄膜（< 5 nm），Ψ 和 Δ 对厚度的灵敏度下降，参数耦合加剧，测量不确定度增大。在此区间，多入射角测量（如 65°、70°、75° 三角度）可显著改善厚度与折射率的解耦效果。使用穆勒矩阵椭偏仪获取退偏信息可区分真实超薄膜厚度与等效表面粗糙度。
5 在光伏领域的应用
5.1 SiNₓ 减反膜检测
PECVD 沉积的 SiNₓ:H 减反膜（75 ~ 85 nm）的厚度和折射率是 PERC/TOPCon 电池的关键工艺控制指标。椭偏仪可同时测量膜厚 d 和折射率 n，反馈调节 SiH₄/NH₃ 流量比和沉积时间。SE 400adv PV 等激光椭偏仪专为 SiNₓ 单层膜检测设计。SiNₓ:H 折射率与硅氮比和氢含量直接关联——n@633 nm 从 1.9 增至 2.3 对应于富硅程度增加，该参数影响减反膜的吸光特性和钝化效果。在线椭偏系统通过多点测量（晶圆中心和四角）监控膜厚横向均匀性，均匀性偏移超过 ±3% 时触发 PECVD 腔体维护。
5.2 超薄 SiOₓ 隧穿氧化层
TOPCon 电池的隧穿氧化层厚度仅 1 ~ 2 nm。椭偏仪是该厚度范围唯一可非接触精确表征的手段。测量时需注意衬底光学模型的准确性——1 nm SiO₂ 的光学响应极易与表面粗糙度混叠。
5.3 非晶硅/多晶硅薄膜
通过 Tauc-Lorentz 模型分析 a-Si:H 和 poly-Si 薄膜的厚度、折射率和消光系数。椭偏结果可反映非晶硅的氢含量（通过折射率关联）和多晶硅的结晶度。非晶硅折射率随氢含量增加而降低——n@633 nm 从无氢 a-Si 的 ~4.3 降至高氢含量（~15 at%）的 ~3.5，该关联可用于非破坏性监测 PECVD 沉积的氢含量趋势。
5.4 i-a-SiOₓ:H 钝化层
在 HJT 和 TOPCon 电池中，超薄 i-a-SiOₓ:H 层的厚度和光学常数对钝化性能至关重要。椭偏仪测量需区分薄膜和界面层的贡献，采用多层模型（Si substrate / SiOₓ interface / SiOₓ:H layer / roughness）拟合。
5.5 透明导电氧化物（TCO）表征
HJT 电池正面的 ITO 膜（通常 80 ~ 100 nm）和钙钛矿-硅叠层电池中的透明电极需要精确控制光学常数。ITO 的折射率随退火温度升高而增大（结晶度提高），消光系数在近红外区受自由载流子吸收影响显著。椭偏仪通过 Drude + Lorentz 组合模型同时获取 ITO 的膜厚、可见光折射率和近红外消光系数，为光学设计和电学性能优化提供量化依据。
5.6 钙钛矿薄膜光学表征
钙钛矿吸光层（通常 400 ~ 800 nm）的折射率和消光系数是器件光学模拟的必要输入参数。椭偏测量在钙钛矿材料中面临额外挑战：钙钛矿组分复杂（如 FA₁₋ₓCsₓPbI₃₋ᵧBrᵧ 四元体系）、薄膜表面粗糙度高（RMS 可达 20 ~ 50 nm）、且部分材料在紫外光下发生相分离。采用 Bruggeman EMA 处理表面粗糙层、多振子 Lorentz/Tauc-Lorentz 模型拟合吸收特征，并控制测量光剂量以避免光致退化。
6 与其他膜厚测量方法的对比
| 方法 | 原理 | 精度 | 优势 | 局限 |
|---|---|---|---|---|
| 椭偏仪 | 偏振态变化 | 亚纳米 | 多参数、非接触、多层膜 | 需模型拟合、不透明膜受限 |
| 光谱反射仪 | 干涉振荡周期 | 纳米级 | 快速、简单 | 单参数（膜厚）、需已知 n |
| 四探针 | 方块电阻 + 电阻率 | 间接 | 电路兼容 | 需已知电阻率、接触式 |
| 台阶仪 | 机械探针扫描 | 亚纳米 | 直接测量 | 接触式、需台阶 |
| SEM/TEM | 截面成像 | 原子级 | 直接观测 | 破坏性、耗时长 |
椭偏仪与反射仪的互补使用是工程中的常见策略：反射仪对透明薄膜膜厚的测量速度快、模型简单，适合大批量在线检测；椭偏仪则在需要同时获取 n 和 k 的场合（如工艺开发阶段确定 SiNₓ:H 折射率目标值）无可替代。
7 总结
椭圆偏振光谱仪通过测量偏振光与薄膜相互作用后 Ψ 和 Δ 的变化，结合物理模型和数学拟合，实现薄膜厚度和光学常数的非接触高精度表征。核心技术挑战在于光学模型的准确建立——模型选择直接影响测量结果的可信度。从菲涅尔系数和特征矩阵法出发的物理建模，到 Cauchy/Lorentz/Tauc-Lorentz/EMA 等色散模型的选择，再到多角度测量和参数解耦策略，椭偏测量是一个完整的方法论链条而非简单的仪器操作。在光伏产业中，椭偏仪是 SiNₓ 减反膜、SiOₓ 隧穿层、a-Si/poly-Si 钝化层、TCO 电极及钙钛矿吸光层等关键薄膜的标准化检测工具，贯穿从研发到量产质量控制的完整流程。