3. 薄膜光学基础

薄膜光学设计的核心任务是回答两个问题：给定膜系结构，其光学特性如何计算？给定目标光学特性，膜系结构如何确定？本章从物理图像出发，建立单层膜与多层膜干涉的统一数学框架——特征矩阵法，并引入导纳图解法与等效界面概念，为后续膜系设计打下坚实的理论基础。理解这些工具的内在联系，比机械记忆公式更为重要。

3.1 单层膜的干涉原理与特征矩阵

当一束光从入射介质（折射率 n₀）垂直或斜入射到厚度为 d₁、折射率为 n₁ 的单层薄膜上，薄膜下方是折射率为 n₂ 的基底时，光在薄膜的两个界面（入射介质/薄膜界面、薄膜/基底界面）上分别发生反射和透射。由于薄膜厚度通常在波长量级，上下界面反射的光束之间会产生固定的相位差，从而发生多光束干涉。这种干涉的强度取决于薄膜的光学厚度（n₁d₁）与波长的比值，以及各界面的反射系数。

为了定量描述这一过程，我们引入光学导纳 η。对于非磁性介质，导纳定义为：

η = N × y₀（垂直入射时）
η = N × y₀ / cosθ（TM偏振，即P光）
η = N × y₀ × cosθ（TE偏振，即S光）

其中 N = n - ik 为复折射率，y₀ 为自由空间导纳（约 1/377 S）。在薄膜光学中，通常将导纳归一化，使自由空间导纳为 1，此时 η = N（垂直入射）。单层膜的特征矩阵将薄膜上下界面的电磁场联系起来：

[B]   [ cosδ₁      (i sinδ₁)/η₁ ] [1]
[C] = [ iη₁ sinδ₁      cosδ₁     ] [η₂]

其中 δ₁ = 2π n₁ d₁ cosθ₁ / λ 为薄膜的相位厚度，θ₁ 为光在薄膜内的折射角。矩阵中的 2×2 矩阵即为单层膜的特征矩阵。通过该矩阵，我们可以求得整个膜系的等效导纳 Y = C/B，进而得到反射系数 r = (η₀ - Y)/(η₀ + Y) 和反射率 R = |r|²。

3.2 多层膜的特征矩阵法：递推与连乘

对于由 m 层薄膜组成的多层膜系，只需将每一层的特征矩阵按光传播的先后顺序连乘，即可得到整个膜系的特征矩阵：

[B]   [ m                    ] [1]
[C] = [ Π Mⱼ                 ] [η_sub]
      [ j=1                 ]

其中 Mⱼ 为第 j 层膜的特征矩阵，η_sub 为基底导纳。矩阵连乘的顺序不能颠倒——先入射到的层在左边。得到总矩阵后，膜系的等效导纳 Y = C/B，反射系数和透射系数分别为：

r = (η₀ - Y)/(η₀ + Y)
t = 2η₀ / (η₀B + C)

透射率 T = Re(η_sub) × |t|² / η₀（垂直入射时简化为 T = 4η₀η_sub / |η₀B + C|²）。对于吸收膜层，折射率为复数，矩阵元素中的 η 和 δ 均需使用复数形式，此时 R + T + A = 1，吸收率 A 可由此求出。

在实际工程中，膜层数可能多达数十层甚至上百层，手算矩阵连乘不现实，通常使用专业软件（如 TFCalc、Essential Macleod、OptiLayer）进行数值计算。但理解矩阵法的物理意义，有助于工程师判断软件输出结果的合理性，以及在膜系优化时选择正确的初始结构。

3.3 导纳图解法：直观判断与设计

导纳图解法是一种在复平面上绘制膜系等效导纳轨迹的图形化方法。对于无吸收的介质膜，等效导纳 Y 始终位于复平面的实轴上（因为无吸收时导纳为实数），但对于吸收膜或斜入射时，Y 可能为复数。导纳图的核心用途是帮助设计者直观地理解膜层对导纳的“搬运”作用。

以单层膜为例，从基底导纳 η₂ 出发，经过相位厚度 δ₁ 的薄膜后，输入端的导纳 Y 满足：

Y = η₁ × (η₂ + iη₁ tanδ₁) / (η₁ + iη₂ tanδ₁)

当 δ₁ = π/2（即光学厚度为 λ/4）时，tanδ₁ → ∞，上式简化为 Y = η₁²/η₂。这就是著名的 λ/4 层导纳变换公式。若 η₁² = η₀η₂，则 Y = η₀，实现零反射——这正是单层减反射膜的设计原理。当 δ₁ = π（即光学厚度为 λ/2）时，Y = η₂，薄膜如同不存在，称为“虚设层”或“半波层”。

对于多层膜，导纳图解法通常采用“史密斯圆图”或“导纳轨迹图”的形式。每一层膜对应一段圆弧（无吸收时为实轴上的线段），设计者可以观察导纳轨迹是否接近目标点（通常是 η₀）。这种方法在早期没有计算机时非常流行，现在虽然计算已由软件完成，但导纳图仍然是理解膜系行为、诊断膜系问题的有力工具。例如，当反射率曲线出现异常峰时，通过导纳图可以快速定位是哪一层膜的相位厚度偏离了设计值。

3.4 等效界面概念：化繁为简的思维

等效界面概念是薄膜光学中一种重要的简化思维。它的核心思想是：任意一个多层膜系，从外部看，都可以等效为一个具有特定导纳和相位特性的“虚拟界面”。这个等效界面的反射特性与原多层膜完全相同。具体来说，将基底和其上所有膜层视为一个整体，用等效导纳 Y 来描述，那么入射介质与这个等效界面之间的反射系数就是 r = (η₀ - Y)/(η₀ + Y)。

等效界面概念的价值在于：

• 简化分析：将复杂的多层膜问题转化为单界面问题，便于快速估算反射率。
• 模块化设计：可以将一个多层膜系视为一个“超级层”，与其他膜系组合时只需考虑其等效导纳。
• 理解膜系作用：例如，高反膜堆可以等效为一个具有高反射率的界面，减反射膜可以等效为一个导纳匹配层。

值得注意的是，等效导纳 Y 是波长的函数，因此等效界面只在特定波长或波段内成立。在宽光谱设计中，需要将不同波长的等效导纳分别考虑，这正是导纳图解法中“色散轨迹”的由来。

3.5 案例拆解：单层 MgF₂ 减反射膜的设计与验证

以最常见的可见光减反射膜为例：在折射率 n₂ = 1.52 的 K9 玻璃基底上，镀制单层 MgF₂ 薄膜（n₁ = 1.38），入射介质为空气（n₀ = 1.0）。设计目标是使 550nm 处的反射率最低。

设计计算：根据 λ/4 条件，最佳光学厚度为 n₁d₁ = λ/4 = 137.5nm，因此物理厚度 d₁ = 137.5 / 1.38 ≈ 99.6nm。此时等效导纳 Y = n₁²/n₂ = 1.38²/1.52 ≈ 1.253。反射率 R = |(1 - 1.253)/(1 + 1.253)|² ≈ 1.26%。相比未镀膜时 R₀ = |(1 - 1.52)/(1 + 1.52)|² ≈ 4.26%，反射率显著降低。

实际镀膜与测试：采用电子束蒸发镀膜，沉积速率控制在 0.3 nm/s，基底温度 250°C，真空度优于 2×10⁻³ Pa。镀膜后用分光光度计测量 400-700nm 波段的反射率。实测结果：550nm 处反射率约 1.5%，与理论值接近；但 450nm 处反射率约 2.8%，偏离理论值。分析原因：MgF₂ 薄膜在短波段的折射率色散（n₁ 随波长减小而增大）导致 λ/4 条件不再严格满足。

优化方案：将设计波长改为 520nm，使 450nm 和 650nm 处的反射率更均衡。重新计算厚度 d₁ = 520/(4×1.38) ≈ 94.2nm。镀膜后实测 450nm 反射率降至 2.1%，650nm 反射率约 1.8%，全波段平均反射率优于 2%。

复盘经验：单层减反射膜只能在一个波长附近实现低反射，且受材料折射率限制（n₁ 需接近 √n₀n₂ ≈ 1.23，但常见材料中 MgF₂ 的 1.38 已是最低）。若要进一步降低反射率或拓宽波段，必须采用多层膜设计。此外，镀膜工艺中的厚度控制误差（±1nm）会导致反射率偏移约 0.1%-0.3%，在精密光学设计中需引入监控系统。

3.6 常见误区与总结

在薄膜光学基础学习中，以下误区需要特别警惕：

• 混淆光学厚度与物理厚度：干涉条件取决于 n×d，而非 d 本身。设计时务必明确是光学厚度还是物理厚度。
• 忽略色散：折射率随波长变化，λ/4 条件只在单一波长严格成立。宽光谱设计必须考虑色散。
• 斜入射时忘记偏振：S 光和 P 光的导纳不同，反射率也不同。在 45° 或更大角度入射时，必须分别计算。
• 矩阵连乘顺序错误：从入射介质一侧开始，先入射的层在左。顺序颠倒会导致结果完全错误。
• 将等效导纳视为常数：等效导纳是波长的函数，不能用于宽光谱的简单估算。