第三章 薄膜光学基础：单层膜的反射率与透射率计算、特征矩阵法、导纳概念、等效界面法
光学薄膜设计的本质，是求解光在多层介质界面上的电磁场边界问题。无论是设计增透膜、高反膜还是分光膜，其底层逻辑都建立在单层膜的干涉模型之上。本章从麦克斯韦方程组的边界条件出发，推导单层膜的反射率与透射率公式，进而引入光学导纳与特征矩阵这两个核心工具，最终将多层膜问题简化为等效界面问题。掌握这些基础，是后续进行膜系优化设计的必要前提。
3.1 单层膜的反射率与透射率计算
3.1.1 核心原理：多光束干涉的叠加
当一束光从入射介质（折射率 n₀）射向镀有单层薄膜（折射率 n₁，几何厚度 d₁）的基片（折射率 n_s）时，光会在薄膜的两个界面（界面1：n₀/n₁，界面2：n₁/n_s）上发生多次反射和透射。这些多次反射的光束在无穷远处（或透镜焦平面上）发生相干叠加，形成最终的反射光和透射光。单层膜的光学行为完全由三个参数决定：薄膜的折射率 n₁、薄膜的几何厚度 d₁ 以及 入射角 θ₀。其中，薄膜的相位厚度 δ₁ = 2πn₁d₁cosθ₁/λ，它决定了光在薄膜中往返一次所积累的相位差。
利用振幅反射系数和透射系数的递推关系，可以推导出单层膜的总振幅反射系数 r 和总振幅透射系数 t。对于无吸收介质，反射率 R 和透射率 T 满足能量守恒关系 R + T = 1。在正入射条件下，反射率的表达式可以简化为：
R = (η₀ - Y)² / (η₀ + Y)²其中 Y 是薄膜与基片组合的等效导纳（将在 3.3 节详细讨论）。这个公式表明，只要计算出等效导纳 Y，反射率就迎刃而解。
3.1.2 判断方法：λ/4 与 λ/2 厚度的特殊行为
在实际膜系设计中，λ/4（四分之一波长）和 λ/2（半波长）厚度具有特殊意义。当薄膜的光学厚度 n₁d₁ = λ/4 时，相位厚度 δ₁ = π/2，此时等效导纳 Y = n₁²/n_s，反射率 R = [(n₀ - n₁²/n_s)/(n₀ + n₁²/n_s)]²。若 n₁² = n₀·n_s，则 R = 0，实现完全增透。当薄膜的光学厚度为 λ/2 时，δ₁ = π，等效导纳 Y = n_s，反射率与未镀膜时相同，即半波长膜层对反射率没有影响，这一特性常用于保护膜或间隔层设计。
| 光学厚度 | 相位厚度 δ₁ | 等效导纳 Y | 反射率特征 | 典型应用 |
|---|---|---|---|---|
| λ/4 | π/2 | n₁²/n_s | R 随 n₁ 变化，可增透或增反 | 增透膜、高反膜基础层 |
| λ/2 | π | n_s | R 与裸基片相同 | 保护层、缓冲层 |
| λ/8 | π/4 | 复数 | R 对波长敏感 | 宽带减反射设计 |
3.1.3 案例拆解：单层 MgF₂ 增透膜设计
以 K9 玻璃（n_s = 1.52）为基片，在空气中（n₀ = 1.0）镀制单层 MgF₂ 薄膜（n₁ = 1.38）。目标波长 λ₀ = 550 nm。根据增透条件 n₁ = √(n₀·n_s) = √(1.52) ≈ 1.233，而 MgF₂ 的折射率 1.38 与之接近但不完全匹配。计算 λ/4 厚度：d₁ = λ₀/(4n₁) = 550/(4×1.38) ≈ 99.6 nm。代入反射率公式：Y = n₁²/n_s = 1.38²/1.52 ≈ 1.253，R = [(1.0 - 1.253)/(1.0 + 1.253)]² ≈ 1.27%。相比裸基片的反射率 R₀ = [(1.0 - 1.52)/(1.0 + 1.52)]² ≈ 4.26%，反射率显著降低。但单层膜无法做到零反射，因为 n₁ 无法精确等于 1.233。若改用折射率更低的材料（如 MgF₂ 掺杂），或采用多层膜设计，可进一步逼近零反射。
3.2 特征矩阵法：多层膜计算的通用工具
3.2.1 核心原理：矩阵描述光场传播
特征矩阵法将每一层薄膜视为一个 2×2 矩阵，该矩阵将膜层一侧的电场和磁场切向分量与另一侧的切向分量联系起来。对于第 j 层膜，其特征矩阵为：
M_j = [cosδ_j, (i sinδ_j)/η_j]
[iη_j sinδ_j, cosδ_j]其中 δ_j = 2πn_jd_jcosθ_j/λ 是第 j 层的相位厚度，η_j 是第 j 层的光学导纳。整个多层膜的特征矩阵是各单层矩阵的乘积：M = M₁ × M₂ × ... × M_m。基片的贡献通过边界条件引入，最终得到整个膜系的等效导纳 Y = C/B，其中 [B; C] = M × [1; η_s]。
3.2.2 操作流程：从矩阵到反射率的五步计算
- 确定参数：列出每层膜的折射率 n_j、几何厚度 d_j、入射角 θ₀，计算各层相位厚度 δ_j。
- 计算导纳：对于正入射，η_j = n_j；对于斜入射，需区分 s 偏振（η_j = n_jcosθ_j）和 p 偏振（η_j = n_j/cosθ_j）。
- 构建矩阵：按顺序写出每层膜的特征矩阵 M_j。
- 矩阵相乘：计算总矩阵 M = M₁M₂...M_m，得到 B 和 C 的值。
- 计算反射率：Y = C/B，r = (η₀ - Y)/(η₀ + Y)，R = |r|²。
这一方法可以方便地通过编程实现，是几乎所有光学薄膜设计软件的核心算法。对于无吸收膜系，透射率 T = 1 - R；若存在吸收，则需考虑消光系数 k，此时导纳变为复数 η = n - ik。
3.2.3 案例拆解：双层减反射膜的矩阵计算
设计一个双层减反射膜：在 K9 玻璃上依次镀制 λ/4 的 ZrO₂（n₁ = 2.05）和 λ/4 的 MgF₂（n₂ = 1.38），目标波长 550 nm。计算过程如下：
δ₁ = δ₂ = 90°
M₁ = [0, i/2.05; i×2.05, 0]
M₂ = [0, i/1.38; i×1.38, 0]
M = M₁ × M₂ = [-2.05/1.38, 0; 0, -1.38/2.05] = [-1.485, 0; 0, -0.673]
B = -1.485, C = -0.673×1.52 = -1.023
Y = C/B = 0.689
R = [(1 - 0.689)/(1 + 0.689)]² ≈ 3.4%虽然双层膜在中心波长的反射率（3.4%）略高于单层 MgF₂（1.27%），但双层膜在更宽的波长范围内保持低反射，且可以通过调整各层厚度来优化宽带特性。实际设计中，常使用三层或更多层膜来同时满足带宽和残余反射率的要求。
3.3 导纳概念与等效界面法
3.3.1 核心原理：光学导纳的物理意义
光学导纳 η 定义为电磁场切向分量之比：η = H_t / E_t。在自由空间中，导纳为 Y₀ = 1/377 S（西门子），但在薄膜光学中，通常使用归一化导纳，即 η = n（正入射时）。导纳的概念将电磁场问题转化为阻抗匹配问题：当入射介质的导纳与膜系等效导纳相等时，反射为零。这类似于传输线理论中的阻抗匹配。
等效界面法的核心思想是：将多层膜系统看作一个整体，用一个具有等效导纳 Y 的单一界面来替代。这样，整个多层膜的反射率就简化为两个半无限介质（入射介质和等效介质）之间的单一界面反射问题。等效导纳 Y 可以通过特征矩阵法求得，也可以逐层递推计算。
3.3.2 判断方法：导纳图与反射率极值
在膜系设计中，导纳图（Admittance Diagram）是一种直观的辅助工具。将膜系从基片开始逐层累加，在复平面上画出等效导纳的轨迹。当轨迹穿过实轴时，对应反射率极值；当轨迹接近入射介质导纳 η₀ 时，反射率接近零。例如，对于 λ/4 膜层，导纳轨迹沿实轴移动；对于 λ/2 膜层，轨迹绕一圈回到原点，不改变等效导纳。
| 膜层类型 | 导纳轨迹 | 对反射率的影响 | 设计用途 |
|---|---|---|---|
| λ/4 高折射率层 | 沿实轴向右移动 | 增加反射率 | 高反膜 |
| λ/4 低折射率层 | 沿实轴向左移动 | 降低反射率 | 增透膜 |
| λ/2 任意层 | 闭合曲线 | 不改变反射率 | 相位调节 |
| 非 λ/4 层 | 复平面曲线 | 引入相位变化 | 宽带设计 |
3.3.3 案例拆解：用等效界面法分析三层高反膜
设计一个中心波长 1064 nm 的三层高反膜，结构为：基片 / (HLH) / 空气，其中 H 为 ZnS（n_H = 2.35），L 为 MgF₂（n_L = 1.38），各层均为 λ/4。使用等效界面法逐步分析：
- 第一步：从基片开始，第一层 H（λ/4），等效导纳 Y₁ = n_H²/n_s = 2.35²/1.52 ≈ 3.63。
- 第二步：叠加第二层 L（λ/4），等效导纳 Y₂ = n_L²/Y₁ = 1.38²/3.63 ≈ 0.525。
- 第三步：叠加第三层 H（λ/4），等效导纳 Y₃ = n_H²/Y₂ = 2.35²/0.525 ≈ 10.52。
- 计算反射率：R = [(1 - 10.52)/(1 + 10.52)]² ≈ 70.1%。
若将膜层数增加到 9 层（HLHLHLHLH），等效导纳会更高，反射率可超过 99%。这一递推过程清晰地展示了等效界面法的优势：无需处理复杂的多光束干涉，只需逐层更新等效导纳即可。
3.4 常见误区与总结复盘
3.4.1 误区一：混淆光学厚度与几何厚度
光学厚度是折射率与几何厚度的乘积（nd），它决定了相位延迟。许多初学者在设计 λ/4 膜层时，误将几何厚度设为 λ/4，导致实际相位厚度错误。正确做法是：d = λ/(4n)。例如，对于 550 nm 波长，MgF₂（n=1.38）的 λ/4 几何厚度约为 99.6 nm，而非 137.5 nm。
3.4.2 误区二：忽略色散与吸收
实际材料的折射率随波长变化（色散），且可能存在吸收（消光系数 k ≠ 0）。在宽光谱设计中，必须使用材料的色散数据，而非单一折射率值。对于吸收膜层，导纳变为复数，反射率计算需使用复数公式，且能量守恒变为 R + T + A = 1（A 为吸收率）。
3.4.3 误区三：正入射公式直接套用于斜入射
斜入射时，s 和 p 偏振的导纳不同，反射率也不同。若不加以区分，会导致设计结果与实际严重偏离。在斜入射设计中，应分别计算两种偏振的反射率，并根据实际使用条件（如自然光、线偏振光）进行加权平均。
3.4.4 总结复盘：单层膜设计的核心逻辑
本章的核心逻辑可以归纳为：一个概念（导纳）、一个工具（特征矩阵）、一个方法（等效界面）。导纳将光学问题转化为阻抗匹配问题，特征矩阵提供了通用的计算框架，等效界面法简化了多层膜的分析过程。在实际设计中，这三者相辅相成：用特征矩阵法进行精确计算，用等效界面法进行快速估算和物理理解，用导纳图进行直观判断和优化方向指导。