1. 从两个参数到 16 个元素：穆勒矩阵补上了哪部分信息
椭偏测量的基本思路是看偏振态的变化。一束偏振光打到样品上，反射或透射之后偏振态被改变，变化的程度由膜层厚度、折射率与结构决定。基础椭偏测量把这个改变归结为两个量：Ψ 描述 p 偏振与 s 偏振的振幅比，Δ 描述两者的相位差。
这两个量之所以够用，是因为默认了一个前提：样品不改变光的偏振度。光进去是偏振的，出来仍是同样程度的偏振，只是振动方向和相位变了。在这个前提下，样品的偏振变换可以用一个 2×2 的琼斯矩阵（Jones matrix）完整描述，Ψ 和 Δ 就是琼斯矩阵归一化之后的两个角度量。
实际样品不总是守这个前提。粗糙界面会把偏振光散射成部分偏振光，厚衬底的背反射会把两束偏振态不一致的光叠在一起，周期性纳米结构会让不同位置的偏振响应互相混合。这些过程都会降低出射光的偏振度，也就是产生退偏振。一旦出现退偏振，2×2 的琼斯矩阵就描述不了，Ψ 和 Δ 也不再是完备的描述。
穆勒矩阵（Mueller matrix）解决的就是这件事。它是一个 4×4 的实数矩阵，作用在斯托克斯矢量（Stokes vector）上：
Sout = M · Sin
斯托克斯矢量的四个分量分别对应总光强、水平与垂直线偏振的强度差、±45° 线偏振的强度差、左右旋圆偏振的强度差。穆勒矩阵把入射的四个分量映射到出射的四个分量，因此能描述部分偏振与退偏振，这是它比琼斯矩阵多出来的能力。
| 对比项 | Ψ、Δ 基础椭偏 | 穆勒矩阵椭偏 |
|---|---|---|
| 输出参数 | 2 个角度量 | 4×4 共 16 个元素 |
| 数学形式 | 2×2 复数琼斯矩阵的归一化 | 4×4 实数矩阵 |
| 能否描述退偏振 | 不能 | 能 |
| 能否描述各向异性 | 仅能处理简单情形 | 可完整表征单轴、双轴各向异性 |
| 建模复杂度 | 低，参数少、收敛快 | 高，需处理元素之间的耦合 |
| 典型样品 | 平整的各向同性介电薄膜 | 液晶、取向聚合物、纳米结构、粗糙或多层散射体 |
可以这样理解两者的关系：Ψ 和 Δ 是穆勒矩阵在无退偏振、各向同性这一子集里的简化形式。样品越简单，简化描述越够用；样品一旦复杂，简化假设失效，就需要回到完整的 16 个元素。技术底座可参考椭偏仪原理深度解析：偏振态测量与薄膜光学常数反演。
2. 16 个自由度的分工：1 个强度、3 个二向色性、3 个相位延迟、9 个退偏振
16 个元素看着多，其实分工是清楚的。按照偏振光学领域的通行划分，一个穆勒矩阵的 16 个自由度可以拆成四组：
| 自由度 | 个数 | 对应的物理量 | 常见来源 |
|---|---|---|---|
| 强度 | 1 | 总透过率或总反射率，与偏振状态无关的那部分 | 材料吸收、膜层透过损失 |
| 二向色性（diattenuation） | 3 | 两个正交偏振态被不等量衰减的程度 | 取向吸收材料、金属界面、偏振器件 |
| 相位延迟（retardance） | 3 | 两个正交偏振态之间的相位差，也就是双折射 | 各向异性膜、液晶、应力双折射 |
| 退偏振（depolarization） | 9 | 偏振度被降低的程度 | 粗糙界面、背反射、散射结构、多层叠加 |
这组划分的价值在于，它回答了选型时经常被问到的一个问题：多出来的 16 个元素里，真正带来新信息的是哪部分。答案是有 12 个：3 个二向色性加 3 个相位延迟，再加上 9 个退偏振。只测 Ψ 和 Δ 的方案，能拿到的主要是相位延迟里的线性部分，另外两类信息基本拿不到。
这里有一件事要提前说清楚，避免误读：这个自由度划分不等于元素与参数的逐一对应。样品同时存在多种偏振效应时，各个元素会互相纠缠，单看某一个元素的数值很难直接下结论。比如线性双折射和线性二向色性会同时出现在右下角的 3×3 子矩阵里，粗糙界面的散射也会叠加进来。要分开，必须做矩阵分解，这一点在第 5 节展开。
3. 按区块读矩阵：二向色性、偏振度、相位延迟分别落在哪
把 4×4 矩阵按位置分块，可以看出不同区块承担的角色。下面的对应关系适用于归一化之后（每个元素除以左上角的 M11，消除总强度的影响）、退偏振不强的样品；带明显退偏振的样品要先分解再判读。
| 元素位置 | 主要对应的偏振信息 | 判读说明 |
|---|---|---|
| M11 | 非偏振光照射时的总透过或总反射 | 归一化基准，反映整体透过损失 |
| 首行 M12、M13、M14 | 二向色性：样品对不同偏振入射光的衰减差异 | M12、M13 对应线性二向色性，M14 对应圆二向色性 |
| 首列 M21、M31、M41 | 偏振度（polarizance）：非偏振光入射时样品产生偏振的能力 | 数值越大，样品对入射光的偏振化作用越明显 |
| 右下 3×3 子矩阵 | 相位延迟与退偏振，两者混在一起 | 其中线性双折射主要落在 M22、M23、M32、M33 |
| M24、M34、M42、M43 | 圆双折射，也就是旋光 | 旋光角可由延迟矩阵的 M24 与 M43 求出 |
| M44 | 圆退偏振相关项 | 与 M22、M33 一起反映退偏振的整体水平 |
常用指标也可以从矩阵直接算出来。以归一化矩阵为基础，线性二向色性对应首行中间两个元素的平方和开方，圆二向色性对应 M14 的绝对值，总二向色性再把 M14 一并计入。偏振度则由首列三个元素算出。旋光角取延迟矩阵 M24 与 M43 之比的反正切。
在实测数据上，这些分块会表现为不同的变化规律。以研究级双旋转补偿器椭偏仪 RC2 的测量为例，旋转测量时右下角九个元素会出现各向异性样品的相位延迟信号，不同区域的对比度分别与衰减和偏振相关，未旋转的相位延迟则集中在另一片区域。这类图案本身就是判断样品各向异性是否被正确采集的直观依据。
这里要提醒一个常见误判：不能拿对角线元素的数值直接当厚度或折射率用。对角线元素在空气直射测量中接近 1，说明仪器定标正常；在样品测量中出现偏离，可能是双折射，也可能是退偏振或二向色性，来源不止一种。真正要拿到膜厚与光学常数，仍然要走建模拟合这条路。
各向异性材料与液晶盒的穆勒矩阵数据示意
4. 四类必须用穆勒矩阵的样品
4.1 各向异性材料：单轴与双轴
取向聚合物、有机半导体膜、部分晶体在不同方向上折射率不同，p 偏振与 s 偏振之外还引入了新的方向差异。这类样品的光学响应无法用单一的 Ψ 和 Δ 描述，需要完整的各向异性层模型，而模型的输入数据就是穆勒矩阵或广义椭偏数据。RC2 的测量类型中，通用椭偏模式对应各向异性样品的完整 2×2 琼斯矩阵，穆勒矩阵模式则给出全部 16 个元素。
4.2 液晶器件与相位延迟器件
液晶盒是穆勒矩阵测量很典型的场景。扭转向列液晶的光轴沿厚度方向连续变化，各向异性和退偏振同时存在，盒厚通常有几微米，层内的偏振混合过程复杂。原厂资料给出的示例层结构是玻璃基底、氧化铟锡 19.6 nm、聚酰亚胺 127 nm、扭转液晶 4.8 μm，再对称叠回去。这种结构里，完整穆勒矩阵可以用来表征光轴的扭转与预倾斜，以及液晶的各向异性折射率。
4.3 带退偏振的样品
退偏振的来源不止一种，常见的包括：
粗糙界面：界面起伏把相干光散射成非相干光，偏振度下降
厚衬底背反射：衬底背面反射的光与表面反射的光叠加，两者偏振态不一致
周期性纳米结构：纳米光栅、超表面这类结构在不同位置产生不同的偏振响应，整体表现为退偏振
多层散射叠层：有机叠层、多孔材料在厚度方向上的不均匀性
这类样品如果只用 Ψ 和 Δ，拟合残差会长期偏大，厚度与折射率出现系统偏差，而且偏差不易被发现。穆勒矩阵测量能把这些退偏振信息记录下来，配合合适的模型分辨它是来自界面、衬底还是结构本身。
4.4 具有手性或旋光的材料
手性分子、螺旋结构、部分生物材料会让线偏振面发生旋转，对应的偏振效应是圆双折射与圆二向色性，主要落在矩阵的 M24、M34、M42、M43 与 M14、M41 这些位置。只测 Ψ 和 Δ 的方案基本看不出这部分信息。
4.5 反过来看：什么情况不必用穆勒矩阵
如果样品是平整的各向同性介电薄膜，单层或常规多层，没有退偏振也没有明显的界面散射，那么 Ψ 和 Δ 足够，建模更快、参数更少、产线通量也更高。穆勒矩阵不是所有场景的默认选项，它是当简化假设不成立时的必要选项。这条边界建议在开测之前跟应用工程师确认一次，避免为了用满功能而把项目复杂化。
5. 拿到 16 个数字之后：极分解与拟合判读
实验测出的穆勒矩阵，各个元素是混在一起的，直接把元素值往材料参数上对应，很容易得出错误结论。偏振光学里通行的做法是做极分解（polar decomposition，工程上也称 Lu-Chipman 分解），把一个复杂的穆勒矩阵拆成三个基本矩阵的乘积：
M = MΔ · MR · MD
三个因子分别对应纯退偏振、纯相位延迟、纯二向色性三种理想元件。拆开之后，三类偏振效应各自的量值就可以单独读出来了：
| 指标 | 含义 | 取值范围与判读 |
|---|---|---|
| 退偏系数 Δ | 由退偏振矩阵的特征值计算，反映整体退偏能力 | 0 表示完全不退偏，1 表示完全退偏 |
| 总相位延迟 R | 由延迟矩阵的迹计算的双折射总量 | 再拆分为线性延迟与圆延迟（旋光） |
| 二向色性 D | 由矩阵首行算出的偏振相关衰减 | 可拆分为线性与圆两部分 |
| 偏振度 P | 由首列算出的偏振化能力 | 与二向色性不一定相等，要分开看 |
分解之后，判读的逻辑就清楚了：相位延迟对应双折射，与材料的取向和各向异性有关；二向色性对应偏振相关的吸收，常出现在有取向的吸收材料或金属界面；退偏系数则提示样品内部存在散射、粗糙或多重反射。
落到拟合环节，还有一条经验值得记住：退偏振如果不建进模型，误差会转移到厚度和折射率上。粗糙界面和背反射导致的退偏振，若在模型里被当作单纯的膜层处理，拟合算法会把无法解释的残差分配给它能调的参数，也就是厚度和 n、k。表现是残差曲线在某些波段持续偏离，厚度值则对测量位置和测量条件特别敏感。把退偏振来源建模，或引入相应的退偏振参数之后，厚度与折射率的重复性通常会明显改善。
RC2 配套的 CompleteEASE 软件提供穆勒矩阵数据的可视化与分解工具，可以直接观察各元素随波长与角度的变化，定位各向异性与退偏振的来源，再把对应结构写进层模型。
CompleteEASE 穆勒矩阵图形工具界面
6. 采全 16 个元素需要什么条件
穆勒矩阵测量不是软件功能，它对硬件有实际要求。要在一次测量中获得足够的偏振信息，需要在光路里引入多组已知的偏振调制状态。
6.1 双旋转补偿器架构
RC2 采用的是样品前后各一个补偿器、两者同步旋转的方案。补偿器持续改变入射光的偏振态，出射侧再分析偏振态，通过多组调制状态组合，一次采集即可得到 4×4 矩阵的全部 16 个元素，不需要更换硬件或分次测量。这一架构也是研究级椭偏仪处理各向异性与退偏振样品时的常见选择。
6.2 用一个指标判断矩阵测量准不准
判断穆勒矩阵测量是否可靠，有一个很直接的办法：测空气直射。此时样品对偏振没有任何影响，理想结果是归一化矩阵等于单位矩阵。RC2 的实测数据是对角线元素为 1 ± 0.002，非对角线元素为 0 ± 0.002。这组数字反映的是矩阵本身的测量准确性，比只看单点厚度重复性更能说明仪器做穆勒测量的能力。同一台设备上，薄氧化膜的厚度重复性可优于 0.005 nm，说明高准确度的矩阵数据能直接转化为厚度结果的稳定。
6.3 光谱与角度配置
RC2 支持多种光谱配置，覆盖范围与波长点数不同，按被测材料的吸收特性与目标膜厚选择：
| 配置 | 光谱范围 | 同步波长点数 | 适用场景 |
|---|---|---|---|
| U / X | 210-1000 nm | 790 | 可见至近红外的常规薄膜表征 |
| D | 193-1000 nm | 800 | 需要覆盖深紫外，适合高 k 介质与光刻相关膜层 |
| UI / XI | 210-1690 nm | 1065 | 含近红外波段，适合厚膜与近红外吸收材料 |
| DI | 193-1690 nm | 1075 | 深紫外到近红外全段覆盖，兼顾超薄膜与厚膜 |
| XI+ | 210-2500 nm | 1040 | 扩展至 2500 nm 近红外，适合窄带隙半导体与厚膜 |
角度方面，RC2 提供固定 65° 入射角、水平自动变角 45°-90° 与垂直自动变角 20°-90° 三种方式。变角的意义在于，各向异性材料的响应随入射角和方位角变化，多角度数据能给层模型提供更多约束；对简单薄膜，固定角测量即可满足。
6.4 采集速度与测量范围
RC2 在全光谱范围内同步采集 1000 个以上波长点，全光谱测量约需三分之一秒。配合聚焦组件，测量光斑可小至 25 × 60 μm；配置自动载物台后可在 300 × 300 mm 范围内做 Mapping，用于评估各向异性或退偏振特征在片内的分布。
7. 总结与延伸
把前面几节收拢成一条判断路径：
先看样品是不是各向同性、无退偏振。是，用 Ψ 和 Δ 的方案，建模更快、通量更高。
再看有没有退偏振的来源。粗糙界面、厚衬底背反射、周期纳米结构都算，有的话需要能采穆勒矩阵的仪器。
确认要测哪些偏振信息。要各向异性轴和相位延迟，看重右下角 3×3 子矩阵；要退偏振，看退偏系数；要旋光，看 M24 与 M43 这类圆偏振相关项。
数据到手后先分解，再建模。极分解把二向色性、相位延迟、退偏振分开，避免把不同来源的偏差都算到厚度上。
用空气直射的归一化矩阵检查数据质量。对角线偏离 1 有多大，决定了后面建模的可信度。
穆勒矩阵的价值不在于元素多，而在于它把退偏振和各向异性这两类信息纳入了可测量的范围。对于遇到复杂样品、拟合长期不稳定、厚度结果对测量条件格外敏感的项目，这套数据往往是打开局面的关键。选型时可以先明确测量对象与关注的信息类型，再确定需要的测量模式和光谱配置。
优尼康科技代理 J.A. Woollam 椭偏仪，可提供 RC2、M-2000、theta-SE 三款的选型咨询、安装培训与售后技术支持，也支持寄送样品做适用性实测。RC2 穆勒矩阵椭偏仪完整采集 4×4 矩阵的 16 个元素，椭偏仪产品线页面列出了全部在售机型与光谱配置。
