光学薄膜基础：光的干涉与薄膜光学原理、单层膜与多层膜、特征矩阵法、反射率与透射率计算基础
光学薄膜设计的本质，是借助光在多层介质界面上反射与透射时的相位叠加，实现对光谱能量分布的精确调控。无论是增透膜、高反膜、分光膜还是滤光片，其底层逻辑都建立在光的干涉理论之上。本章将从物理图像出发，建立单层膜与多层膜的数学模型，并重点掌握特征矩阵法这一通用计算工具，为后续软件操作中的膜系设计、优化与公差分析打下坚实的理论基础。
核心原理：光的干涉与薄膜光学基础
光在两种折射率不同的介质界面上发生反射和折射时，反射光与透射光的振幅比例由菲涅尔系数决定。当介质层厚度与光波长可比拟时，从薄膜上下表面反射的光会发生相位差，进而产生干涉。设薄膜物理厚度为 d，折射率为 n，入射角为 θ，则相邻两束反射光的光程差为 2nd·cosθ。当光程差为波长的整数倍时，反射光加强；为半波长的奇数倍时，反射光减弱。
薄膜光学中，通常引入光学厚度（n·d）和相位厚度 δ = (2π/λ)·n·d·cosθ 来描述膜层对光波的影响。对于垂直入射（θ=0），相位厚度简化为 δ = 2πnd/λ。当 nd = λ/4 时，δ = π/2，这是薄膜设计中最常用的四分之一波长堆叠条件。理解这一物理图像是后续使用软件进行膜系设计的基础，因为绝大多数光学薄膜的初始结构都基于 QWOT（Quarter-Wave Optical Thickness）堆叠。
- 光学厚度（OT）：n·d，决定相位延迟。
- 相位厚度（δ）：2π·n·d·cosθ/λ，决定干涉状态。
- 导纳（η）：对 s 偏振 η = n·cosθ，对 p 偏振 η = n/cosθ。
- 等效导纳（Y）：多层膜整体对入射光的导纳响应。
单层膜与多层膜：从物理图像到数学模型
单层膜是最简单的薄膜系统。设入射介质导纳为 η₀，膜层导纳为 η₁，基底导纳为 η₂，则单层膜的特征矩阵可表示为：
M₁ = [cosδ₁, (i·sinδ₁)/η₁]
[i·η₁·sinδ₁, cosδ₁]整个膜系的特征矩阵为 M = M₁，则膜系的等效导纳 Y = C/B，其中 [B; C] = M · [1; η₂]。反射率 R = |(η₀ - Y)/(η₀ + Y)|²。对于四分之一波长膜层（δ=π/2），矩阵简化为：
M₁ = [0, i/η₁]
[i·η₁, 0]此时等效导纳 Y = η₁²/η₂。若 η₁ < η₂，则 Y < η₂，反射率降低，表现为增透效果；若 η₁ > η₂，则 Y > η₂，反射率升高，表现为高反效果。这就是单层增透膜和高反膜的基本原理。
多层膜是单层膜的推广。对于 q 层膜系，总特征矩阵为各层矩阵的连乘：
[B] [cosδ₁ i·sinδ₁/η₁] [cosδ₂ i·sinδ₂/η₂] [1]
[C] = [i·η₁·sinδ₁ cosδ₁] [i·η₂·sinδ₂ cosδ₂] ··· [η_s]矩阵连乘的顺序必须与光传播方向一致，即从入射侧到出射侧依次相乘。这一方法将复杂的多次反射问题转化为简单的矩阵乘法，是几乎所有光学薄膜设计软件（如 TFCalc、Essential Macleod、OptiLayer）的核心算法基础。
| 膜系类型 | 典型结构 | 特征矩阵特点 | 应用场景 |
|---|---|---|---|
| 单层增透膜 | 基底 / λ/4 低折射率层 | 单矩阵，Y = η₁²/η₂ | 眼镜片、显示盖板 |
| 双层增透膜 | λ/4 高折射率 + λ/4 低折射率 | 双矩阵连乘，可调带宽 | 精密光学镜头 |
| 高反膜 | (HL)^p H 或 (HL)^p | 周期结构，反射带平坦 | 激光腔镜、反射镜 |
| 带通滤光片 | 腔体结构 (HL)^p H (LH)^p | 法布里-珀罗腔效应 | 光谱仪、通信滤光片 |
特征矩阵法：计算流程与操作步骤
在软件中，特征矩阵法的计算流程是自动完成的，但理解其步骤有助于排查设计问题和理解优化结果。标准流程如下：
- 定义膜系结构：确定每层膜的折射率 n、物理厚度 d（或光学厚度）、入射角 θ。
- 计算各层相位厚度：δⱼ = (2π/λ)·nⱼ·dⱼ·cosθⱼ，注意色散时 n 随 λ 变化。
- 构建各层特征矩阵：按上述公式生成 2×2 复数矩阵。
- 矩阵连乘：按入射顺序依次相乘，得到总矩阵 [B; C]。
- 计算等效导纳：Y = C/B。
- 计算反射率与透射率：R = |(η₀ - Y)/(η₀ + Y)|²，T = 1 - R（无吸收时）。
在软件操作中，用户只需输入膜层参数和波长范围，软件自动完成上述计算。但需要注意：矩阵连乘的顺序错误是初学者在手动验证时最常见的错误。例如，设计一个 5 层膜系，从空气入射到玻璃基底，矩阵顺序应为 M₁·M₂·M₃·M₄·M₅，其中 M₁ 是最靠近空气的膜层。
案例拆解：单层增透膜的设计与验证
以常见的 532nm 激光增透膜为例，基底为 BK7 玻璃（n_s = 1.52），目标是在空气（n₀ = 1.0）中实现反射率低于 0.5%。选择 MgF₂ 作为镀膜材料（n₁ = 1.38）。
初始设计：采用单层 λ/4 膜，光学厚度 nd = λ/4 = 133nm，物理厚度 d = 133/1.38 ≈ 96.4nm。理论反射率计算：
Y = n₁²/n_s = 1.38²/1.52 = 1.253
R = |(1.0 - 1.253)/(1.0 + 1.253)|² = |(-0.253)/2.253|² ≈ 0.0126 → 1.26%该结果虽优于未镀膜的 4.2% 反射率，但未达到 0.5% 的目标。此时需要引入多层膜设计。在软件中，我们可以使用双层增透膜结构：在基底上先镀一层 λ/4 的 Al₂O₃（n=1.65），再镀一层 λ/4 的 MgF₂（n=1.38）。
软件操作步骤：
- 新建膜系，设置入射介质为空气（n=1.0），基底为 BK7（n=1.52）。
- 添加第一层：Al₂O₃，光学厚度 133nm（λ/4 @ 532nm）。
- 添加第二层：MgF₂，光学厚度 133nm。
- 设置波长范围 400-700nm，步长 1nm。
- 运行计算，查看 532nm 处的反射率。
结果对比：双层膜在 532nm 处的反射率降至约 0.3%，满足设计要求。但带宽较窄，在 450nm 和 650nm 处反射率上升至 1% 以上。若需要更宽的低反射带，可进一步优化各层厚度，或增加第三层膜。
常见误区与总结复盘
在薄膜光学基础学习中，以下几个误区需要特别警惕：
- 混淆光学厚度与物理厚度：软件中默认输入的是物理厚度（nm），但设计时通常以光学厚度（QWOT）为单位。在 TFCalc 中，可以通过设置“Optical Thickness”模式来直接输入 QWOT 值。
- 忽略色散：折射率随波长变化，尤其在可见光到近红外波段。设计时必须使用材料的色散数据，否则计算出的光谱曲线会与实际镀膜结果偏差较大。
- 入射角的影响：斜入射时，s 和 p 偏振的导纳不同，反射率也不同。如果系统工作在非垂直入射状态，必须在软件中设置正确的入射角，并分别查看 s、p 偏振的平均反射率。
- 吸收与散射：特征矩阵法默认膜层无吸收。对于金属膜或含吸收材料的膜系，需要使用复折射率（n - ik），此时透射率 T ≠ 1 - R，还需考虑吸收率 A。