第二节 光学基础：光的干涉、反射、透射、吸收原理，薄膜光学的基本方程与特征矩阵
光学镀膜的本质，是在基片表面构建一层或多层具有特定厚度与折射率的薄膜，利用光在这些界面上的行为来实现增透、增反、分光或滤光。要真正掌握镀膜工艺，首先必须透彻理解光与薄膜相互作用的物理基础。本节将系统梳理光的干涉、反射、透射与吸收四大基本现象，并引入薄膜光学中描述多层膜行为的两大核心工具——特征矩阵与基本方程。这些内容不仅是膜系设计的理论起点，更是分析镀膜故障（如膜层脱落、光谱漂移、吸收异常）时的底层判断依据。
上图概括了本节的核心逻辑：入射光进入薄膜系统后，在每一个界面同时发生反射与透射，膜层内部还存在吸收损耗；而特征矩阵则将这些复杂的光学行为浓缩为简洁的数学表达，使多层膜的计算成为可能。下面我们将逐一展开这些概念。
2.1 光的干涉、反射、透射与吸收——镀膜世界的四大基石
2.1.1 核心原理：光与物质的三种基本交互
当一束光照射到薄膜表面时，能量会按比例分配到三个通道：反射（R）、透射（T）和吸收（A），且满足能量守恒关系：R + T + A = 1。在镀膜领域，我们最关心的是如何通过膜层设计来精确控制R和T的分配比例，同时尽量降低A。
- 反射：光在折射率突变的界面上发生方向改变。反射率由菲涅尔公式决定，与界面两侧的折射率差直接相关。折射率差越大，反射越强。例如，空气（n≈1.0）与玻璃（n≈1.52）界面的自然反射率约为4%。
- 透射：光穿过界面进入另一种介质。透射率与反射率互补，且受膜层厚度和折射率周期性影响。
- 吸收：光能量被膜层材料转化为热能或其他形式能量。吸收与材料的消光系数k直接相关，k值越大，吸收越严重。金属膜（如Al、Ag）的k值较大，而介质膜（如SiO₂、TiO₂）的k值通常极小（<0.001）。
- 干涉：当两束或多束相干光在空间叠加时，产生明暗相间的强度分布。在薄膜中，界面1的反射光与界面2的反射光（经膜层内部传播后）发生干涉，从而改变整体反射率。
2.1.2 判断方法：如何快速评估膜层的光学行为
在实际镀膜工艺中，我们通常通过以下步骤判断膜层的光学特性：
- 确认材料折射率（n）和消光系数（k）：通过椭偏仪或分光光度计测量。n决定相位变化，k决定吸收损耗。
- 计算光学厚度（nd）：光学厚度 = 物理厚度 × 折射率。λ/4膜层（nd = λ/4）是最基本的膜层单位，其相位延迟为90°。
- 判断干涉类型：若两束反射光的相位差为π的奇数倍，则相消干涉（反射率降低）；若相位差为2π的整数倍，则相长干涉（反射率升高）。
- 评估吸收影响：当k > 0.01时，吸收不可忽略，需在膜系设计中预留补偿。
| 现象 | 物理条件 | 镀膜应用 | 典型材料 |
|---|---|---|---|
| 相消干涉 | 2nd·cosθ = (m+1/2)λ | 增透膜、减反膜 | MgF₂, SiO₂ |
| 相长干涉 | 2nd·cosθ = mλ | 高反膜、反射镜 | TiO₂, Ta₂O₅ |
| 吸收损耗 | k > 0.01 | 中性密度滤光片 | Ni, Cr, Ge |
| 偏振分离 | θ > 45° | 偏振分光棱镜 | 多层介质膜 |
2.1.3 案例拆解：单层增透膜的干涉机制
以最常见的MgF₂单层增透膜为例（n=1.38，镀于n=1.52的BK7玻璃上），设计波长λ=550nm：
- 初始现象：未镀膜时，玻璃表面反射率约4.2%。
- 膜层设计：选择光学厚度nd = λ/4 = 137.5nm，物理厚度d = 137.5/1.38 ≈ 99.6nm。
- 干涉分析：界面1（空气/MgF₂）反射光与界面2（MgF₂/玻璃）反射光的光程差为2nd = λ/2，相位差为π，两束反射光相消干涉。
- 结果：中心波长550nm处反射率降至约1.3%，透射率提升至98.7%。
- 代价：单层增透膜仅在中心波长附近有效，偏离中心波长时反射率迅速上升。
2.2 薄膜光学的基本方程与特征矩阵——多层膜计算的数学引擎
2.2.1 核心原理：从单层到多层的数学抽象
单层膜的光学行为可以用一个2×2矩阵完整描述，这就是特征矩阵。对于多层膜，只需将各层矩阵按顺序相乘，即可得到整个膜系的等效矩阵。这一方法由Abeles于1940年代系统化，至今仍是光学薄膜设计的标准工具。
单层膜的特征矩阵为：
M = [cosδ, (i·sinδ)/η]
[i·η·sinδ, cosδ]其中，δ = 2π·n·d·cosθ/λ 为相位厚度，η为修正导纳（对TE偏振η = n/cosθ，对TM偏振η = n·cosθ）。
对于由p层膜组成的多层膜系，总特征矩阵为：
M_total = M₁ × M₂ × ... × M_p = [A, B]
[C, D]膜系与基片组合后的导纳为Y = C/B，反射系数r = (η₀ - Y)/(η₀ + Y)，反射率R = |r|²。透射率T = 4η₀·Re(Y)/|η₀ + Y|²（假设基片无吸收）。
2.2.2 操作流程：用特征矩阵计算膜系光谱
- 输入参数：各层膜的折射率n、物理厚度d、入射角θ、波长λ、基片折射率n_s。
- 计算各层相位厚度：δ_j = 2π·n_j·d_j·cosθ_j/λ，其中θ_j由斯涅尔定律确定。
- 构建各层特征矩阵：按上述公式写出每层的M_j。
- 矩阵相乘：按膜层顺序依次相乘，得到总矩阵元素A、B、C、D。
- 计算导纳与反射率：Y = C/B，r = (η₀ - Y)/(η₀ + Y)，R = |r|²。
- 扫描波长：在目标波段内重复步骤2-5，即可得到光谱曲线。
2.2.3 案例拆解：双层增透膜的矩阵计算
设计一个双层增透膜：MgF₂（n₁=1.38，d₁=99.6nm）/ ZrO₂（n₂=2.05，d₂=67.1nm），镀于BK7玻璃（n_s=1.52）上，设计波长λ=550nm，正入射。
- 初始目标：将550nm处反射率降至0.5%以下。
- 计算过程：
- δ₁ = 2π×1.38×99.6/550 = 90°（λ/4层）
- δ₂ = 2π×2.05×67.1/550 = 90°（λ/4层）
- M₁ = [0, i/1.38; i×1.38, 0]
- M₂ = [0, i/2.05; i×2.05, 0]
- M_total = M₁ × M₂ = [-1.38/2.05, 0; 0, -2.05/1.38]
- Y = C/B = (-2.05/1.38)/(-1.38/2.05) = (2.05/1.38)² ≈ 2.21
- r = (1 - 2.21)/(1 + 2.21) = -0.377，R = 14.2%
- 发现问题：计算结果显示反射率高达14.2%，远高于目标。原因在于两层λ/4膜的组合并不自动产生增透效果，需要调整厚度。
- 优化方案：将ZrO₂层厚度调整为λ/2（d₂=134.2nm），此时M₂ = [-1, 0; 0, -1]，总矩阵M_total = M₁ × (-I) = -M₁，Y = 1.38²/1.52 ≈ 1.25，r = (1-1.25)/(1+1.25) = -0.111，R = 1.23%。
- 结果对比：调整后反射率从14.2%降至1.23%，接近目标。
- 复盘经验：λ/4-λ/4组合并非万能，必须通过矩阵计算验证。实际设计中常采用λ/4-λ/2或非规整厚度组合来获得更宽带的增透效果。
2.3 常见误区与总结复盘
| 误区 | 错误认知 | 正确理解 |
|---|---|---|
| 误区一 | 认为膜层越厚，反射率越高 | 反射率由干涉条件决定，与厚度呈周期性关系，并非单调递增 |
| 误区二 | 忽略吸收的影响 | 在红外或紫外波段，即使k=0.001的介质膜，多层叠加后吸收也可能显著 |
| 误区三 | 正入射结果直接套用于斜入射 | 斜入射时S/P偏振分离，需分别计算，且膜层有效厚度随角度变化 |
| 误区四 | 认为特征矩阵只能用于规整膜系 | 特征矩阵适用于任意厚度组合，是通用计算方法 |
总结复盘：本节内容为后续所有膜系设计奠定了理论基础。核心要点可归纳为：干涉是手段，矩阵是工具，吸收是约束。在实际镀膜工艺中，工程师应养成“先计算、后镀膜、再验证”的习惯，利用特征矩阵法在设计阶段预测光谱性能，避免盲目试镀。同时，要始终关注材料的光学常数（n、k）随工艺条件的变化，这是理论计算与工程实际之间最重要的桥梁。