1. 为什么多层膜解谱比单层膜难那么多
单层膜的情况很简单。你已知膜层材料的折射率 n(λ) 和消光系数 k(λ)，接下来只要做四步：
用菲涅尔方程算出空气/膜层/基底三层结构的理论反射光谱 R(λ, d)
拿这个理论光谱跟实测光谱对比
调整厚度 d，让两者尽量重合
得到膜厚
整个问题只有一个未知数，几百个波长点的数据都在约束它。这种叫超定问题，信息量远大于未知数，解是唯一的。
多层膜就完全是另一回事了。拿 5 层膜举例：每层 1 个厚度未知数，再加上每层材料至少 2-3 个色散参数（如果用 Cauchy 模型），总共 15-25 个自由参数。反射光谱虽然有几百个波长点的数据，但这些数据点高度相关，相邻波长的反射率不是独立的，你增加波长采样点并不能等比增加信息量。
多层膜解谱真正要解决的问题是：自由参数随层数线性增长，但光谱中的独立信息量增长得慢得多。两条线交叉的地方，就是你能解的层数上限。
下面我们从三个环节来拆这件事：怎么算（TMM）、怎么压参数（色散模型）、怎么搜（拟合算法）。
2. 传输矩阵法（TMM）：把多层膜算成一条光谱
2.1 每层膜就是一个 2×2 矩阵
传输矩阵法（Transfer Matrix Method, TMM）是目前算多层膜反射率的标准工具。它的思路很直接：每一层薄膜对光的传播效果，用一个 2×2 的复数矩阵描述。把所有层的矩阵按顺序乘起来，得到整个膜系的总矩阵，从总矩阵里直接提取反射率和透射率。
第 j 层膜（厚度 dⱼ，复折射率 ñⱼ = nⱼ + ikⱼ）的特征矩阵长这样：
    [ i·ηⱼ·sin(δⱼ)   cos(δⱼ) ]
δⱼ 是相位厚度，ηⱼ 是修正光学导纳：
ηⱼ = ñⱼ · cos(θⱼ)  （s偏振）
ηⱼ = ñⱼ / cos(θⱼ)  （p偏振）
N 层膜的总传输矩阵等于各层矩阵的乘积：
从总矩阵的四个元素 m₁₁, m₁₂, m₂₁, m₂₂ 可以直接算出反射率：
2.2 为什么大家都用 TMM
TMM 有几个明显优势。第一是数值稳定，不像某些递推方法（比如 Rouard 法）在厚膜或高吸收膜里会发散。第二是层数不受限，理论上多少层都能算，实际瓶颈只在于拟合能不能收敛。第三是角度和偏振天然支持，改一下 δⱼ 和 ηⱼ 的表达式就行。第四是梯度可以解析推导，这对拟合算法的效率很关键。
韩国标准科学研究院（KRISS）在 2025 年的一篇综述里提到，基于 TMM 的多层膜反射模型加优化算法，已经能在 10 层以内的透明介质膜系上做到高精度解谱，前提是各层材料的光学常数事先标定好，层厚在 20 nm 到 5 μm 之间。
3. 色散模型：3个参数替掉100个未知数
3.1 参数是怎么炸掉的
如果没有色散模型，每一层膜在每一个波长处的 n 和 k 都是独立未知数。380-1050 nm 范围、1 nm 步长，每层有 671×2 = 1342 个未知数，5 层膜就是 6710 个。从一条反射光谱里解出几千个参数是不可能的。
色散模型做的事就是用 3-6 个物理参数来描述整个波长范围内的 n(λ) 和 k(λ)，把未知数从几千压缩到个位数。
3.2 三种模型，各管一摊
| 模型 | 参数数量 | 适用材料 | n(λ) 表达式 | 局限 |
|---|---|---|---|---|
| Cauchy | 3（A, B, C） | 透明介质（k≈0）：SiO₂、Al₂O₃、MgF₂、光刻胶 | n(λ) = A + B/λ² + C/λ⁴ | 不能描述吸收（k=0）；带边附近不准 |
| Sellmeier | 3-5（B₁,B₂,B₃, C₁,C₂,C₃） | 宽透明波段材料：SiO₂、CaF₂、GaN、ITO | n²(λ) = 1 + Σ Bᵢλ²/(λ²−Cᵢ²) | 不含 k；参数对拟合初值敏感 |
| Tauc-Lorentz | 4-6（A, E₀, C, Eg, ε∞...） | 半导体、吸收膜：a-Si、TiO₂、Ta₂O₅、SiNx | ε₂(E) 用 Lorentz 振子 + Tauc 带边，ε₁ 由 K-K 积分 | 参数多、拟合难度大 |
3.3 色散模型的坑
色散模型减少了参数，但也带来了风险：如果模型的函数形式跟材料的真实色散对不上，拟合出的厚度会有系统性偏差。比如用 Cauchy 模型拟合 TiO₂ 薄膜在 400 nm 附近（接近带边，有微弱吸收）的光谱——Cauchy 模型假定 k=0，实际 k 不为零，拟合出来的厚度可能偏小 1-3%。
一个更稳妥的做法是：先用椭偏仪在单层校准片上建立每种材料的光学常数数据库（n/k vs λ），把这些数据作为固定输入代入多层膜 TMM 模型，只让厚度参与拟合。业内管这叫光学常数预标定，是多层膜解谱里最有效的降维手段。
4. 拟合策略：从局部最优跳到全局
4.1 最小二乘法能干活，但不够
拟合的数学目标是最小化均方误差：
p₁...pm 是待拟合的参数——各层厚度加上色散参数。
薄膜拟合最常用的是 Levenberg-Marquardt（LM）算法。它在远离最优点时像梯度下降，稳；靠近最优点时像 Gauss-Newton，快。但 LM 有一个硬伤：它只能找到离你给的初始值最近的那个局部最优，不保证是全局最优。
4.2 多层膜的 MSE 地形有多复杂
多层膜的 MSE 在参数空间里有大量局部极小值。举个例子，两层膜如果折射率相近（像 SiO₂/SiNx），你把第一层加厚一点、第二层减薄一点，总光学厚度可能互相抵消，反射光谱看起来差不多，这就是参数耦合。
参数耦合让 MSE 空间里形成很多浅而长的沟壑。沿着沟壑方向 MSE 变化很小，LM 算法在里面走不动，精度大打折扣。
4.3 需要全局搜索的时候
| 方法 | 原理 | 优势 | 局限 | 典型耗时（10层膜） |
|---|---|---|---|---|
| LM（局部） | 梯度引导的迭代优化 | 收敛快、精度高 | 依赖初始值，容易卡在局部最优 | ＜1 秒 |
| 差分进化（DE） | 种群随机搜索 + 变异/交叉/选择 | 不依赖梯度，全局搜索强 | 收敛慢，精度不如 LM | 10-60 秒 |
| 模拟退火（SA） | 概率性接受更差解，逐步降温 | 理论上能逃出任何局部最优 | 降温策略调参敏感 | 30-120 秒 |
| 遗传算法（GA） | 选择/交叉/变异 | 适合离散+连续混合参数 | 计算量大，收敛不稳定 | 60-300 秒 |
| DE + LM 混合 | DE 粗搜 + LM 精修 | 兼顾全局搜索和局部精度 | 需要合理设 DE 终止条件 | 5-30 秒 |
实际推荐 DE + LM 两步走。先用差分进化在宽范围里跑一轮全局搜索，拿到一个接近全局最优的粗解；然后把这个粗解丢给 LM 做精修。在 5-8 层膜的典型场景下，这套组合拳可以把收敛成功率从 30-50% 拉到 90% 以上。Filmetrics 软件里的自动拟合模式，逻辑跟这个类似。
4.4 光看 MSE 不够
MSE 低不代表解对了。还需要检查几件事：
残差图：实测减拟合，看残差随波长的分布。如果残差是系统性的波动而不是随机噪声，模型可能漏掉了什么——比如没建模的界面层或者色散模型不对。
参数敏感性：把某个参数改 ±5%，光谱变化大不大？如果几乎不变，这个参数没被光谱约束，拟合出来的值不可靠。
物理合理性：拟合出的折射率在那个材料的正常范围里吗？厚度是正的吗？如果跑出 d = −3 nm，MSE 再小也是错的。
参数相关性：两个参数的相关性系数超过 0.95，说明它们高度耦合，调一个就能补偿另一个。这时候应该固定其中一个（比如预先标定好光学常数）。
5. 层数陷阱：为什么30层比10层难100倍
5.1 不是参数多，是信息不够
很多人以为多层膜解谱难是因为"层数多了算不动"。算力不是瓶颈，瓶颈在信息量。反射光谱里真正独立的特征，干涉峰、谷、拐点。数量是有限的，并不随层数增加而增加。
380-1050 nm 波段内，一条典型反射光谱大约有 5-15 个明显的干涉振荡周期，每个周期大约能约束 1-2 个独立参数。所以光谱能可靠提取的参数上限大概在 10-30 个。一旦你的自由参数超过这个数，光谱就开始"不够用了"。
| 膜层数 | 未知数（厚度+Cauchy参数） | 光谱独立信息量 | 可解性 | 建议策略 |
|---|---|---|---|---|
| 1-2 层 | 2-6 | 10-20 | 高度可解 | 直接拟合，结果可信 |
| 3-5 层 | 6-20 | 10-20 | 可解 | 光学常数预标定，减少自由参数 |
| 6-10 层 | 12-40 | 10-20 | 条件可解 | 必须固定光学常数，需要合理初值 |
| 11-20 层 | 22-80 | 10-20 | 困难 | 仅当膜系有周期性（λ/4 堆叠）时可解 |
| 20+ 层 | 40+ | 10-20 | 通常不可解 | 需椭偏仪或多角度反射补充信息 |
5.2 规整膜系是例外
有一个重要的例外，规整膜系，比如 λ/4 高反膜堆叠。拿 (HL)¹⁵H 的 31 层高反膜来说，所有 H 层（Ta₂O₅）厚度理论上一样，所有 L 层（SiO₂）厚度也一样。31 层膜实际上只有 2 个厚度参数 + 2 组色散参数 ≈ 8 个自由参数，远小于 31×3 = 93。利用这个周期性约束，反射光谱完全可以解。
这也解释了为什么镀膜工艺监控中，发现成品光谱不对的时候，工程师一般不会去逐层解厚度，而是用等效折射率或整体厚度缩放因子来诊断。
6. 实操中管用的几条经验
以下几条是实际工作中反复验证过的，按优先级排。
1. 先标定光学常数，再解厚度
镀多层膜之前，用同样工艺先镀一片单层校准片。拿椭偏仪或反射法 + Cauchy 模型把材料在你设备上的实际 n(λ) 测出来。然后把这个 n(λ) 作为固定参数代入多层膜拟合，只让厚度参与拟合。这一步减少的参数最多（每层减 2-3 个），对拟合稳定性的提升也最明显。
2. 逐层镀、逐层测
如果工艺允许，每镀一层就测一次反射光谱。这样每次面对的都是单层膜 + 已知基底的问题，不用同时拟合多层。退一步，也可以在镀完 2-3 层后测一次，只拟合最近沉积的层，之前的层作为固定参数。
3. 给拟合算法一个好起点
多层膜拟合对初始值非常敏感。初值来源有几档：
QCM 读数，±5% 的偏差已经够 LM 收敛到正确解
镀膜时间 × 标定速率，更粗糙但能用
SEM 断面，破坏性但绝对可靠，适合工艺开发阶段建立基准
4. 扩展波长范围
反射光谱的独立信息量跟波长范围成正比。标准 380-1050 nm 不够用的时候，扩到 190-1700 nm（紫外+近红外）可以把独立信息量提升 50-100%。紫外区的吸收特征和近红外区的长周期干涉能提供额外的约束。
5. 看残差形状，别看残差值
拟合完了一定要画残差图（实测 − 拟合 vs 波长）。随机分布的残差说明模型没问题；周期性波动的残差暗示有没建模的膜层，比如界面层或表面污染；特定波段残差突然变大，可能是色散模型在那个波段不适用。
6. LM 不收敛就上 DE 开路
LM 反复不收敛或每次结果不一样，说明初值离全局最优太远。换差分进化跑 100-500 轮全局搜索，把 DE 找到的解作为 LM 的初值。
7. 承认解不了，换方法
面对 15+ 层、非规整、各层材料光学常数相近的多层膜，而且没有预标定的光学常数，垂直入射反射法大概率解不出来。别强行拟合，那只会得到一个数学上完美但物理上错的答案。考虑换椭偏仪（双通道信息）、做多角度反射、或者直接上截面 SEM/TEM。
7. 椭偏仪还是反射法？看层数选
这是被问到最多的问题之一。直接给结论。
光谱反射法（SR）
多层膜场景中的表现：
速度快（＜1秒/点），适合 Mapping
光路简单，对振动和样品位置不敏感
1-5 层膜且光学常数已知，精度够
6-10 层膜需要预标定光学常数
10+ 层非规整膜系信息量不够
超薄膜（＜10 nm）信号弱
光谱椭偏仪（SE）
多层膜场景中的表现：
Psi + Delta 双通道，信息量翻倍
同时获得 n/k 和厚度
超薄膜灵敏（＜1 nm 可测）
10-30 层膜系仍有解析能力
速度较慢（2-30秒/点）
光路复杂，需要精确对焦和角度校准
按层数选：
≤ 5 层，光学常数已知 → 反射法，快且够用
≤ 5 层，光学常数未知 → 椭偏仪，同时解 n/k 和厚度
6-15 层 → 椭偏仪优先；反射法只在光学常数已预标定且有合理初值时可用
15-30 层（规整膜系）→ 椭偏仪，利用周期性约束减参数
15-30 层（非规整）→ 椭偏仪 + 多角度测量，反射法不推荐
30+ 层 → 单一方法都不够，多方法联合（截面分析 + 光学方法交叉验证）
2026 年 7 月仪器信息网上有一篇实测对比值得参考：Si 衬底上的 SiO₂ 和 Al₂O₃ 薄膜（10 nm - 2000 nm），椭偏仪和反射仪的结果在所有样品上一致，偏差都在各自的不确定度范围内。但椭偏仪在薄层（＜500 nm）上不确定度更低，反射仪在 1000 nm 以上的厚膜上反而占优。
8. 总结
多层膜反射光谱解谱本质上是一个信息量管理问题，不是调参技巧问题。自由参数随层数线性增长，但光谱中独立信息量增长缓慢，两条线的交叉点就是你能解的层数上限。
解决思路分三个方向：
减参数：用色散模型把光学常数压缩到 3-6 个参数；做好光学常数预标定，干脆不让色散参数参与拟合；利用规整膜系的周期性把 N 层膜的独立厚度参数压缩到 2-3 个。
换算法：DE + LM 两步走，全局搜索打头阵，局部精修收尾；合理设参数边界（厚度 > 0、折射率在物理合理范围内）；能用解析梯度就不用数值差分。
加信息：扩展波长范围（UV → NIR）；增加测量角度；上椭偏仪拿双通道数据；逐层沉积逐层测量。
比会解更重要的是知道什么时候解不了。15 层非规整膜系 + 未知光学常数 + 一条垂直入射反射光谱，信息根本不够。这时候强行拟合出来的 MSE=0.001，只能说明过拟合了。
日常工作中最划算的投入是：花半天时间，给每种镀膜材料建一个工艺-光学常数数据库。镀一片单层校准片，做一次椭偏仪测量，这个数据以后每次多层膜拟合都能用。它把"几乎不可解"变成了"可靠可解"。
参考来源：Kim et al., "A Review of Thin-film Thickness Measurements using Optical Methods", Int. J. Precis. Eng. Manuf. (2024)；KRISS 多层膜测量不确定度评估方法研究 (2025)；复享光学 MetronFilm 33层复合多层膜实测案例 (2025)；Horiba Cauchy & Sellmeier Dispersion Technical Notes；J.A. Woollam Ellipsometry Tutorial — Optical Constants；仪器信息网《纳米级测厚：光谱反射与椭圆偏振》(2026)；Flexfilm 全光谱椭偏仪技术综述 (2025)。
