3、薄膜光学基础：单层膜反射与透射、菲涅尔公式、导纳法、特征矩阵法、等效界面概念
光学薄膜设计的本质，是在基底表面构建一层或多层具有特定光学厚度的介质或金属薄膜，通过界面处的多次反射与透射干涉，实现对入射光振幅、相位和偏振态的控制。要理解这一过程，必须从单层膜的电磁场行为入手。本章将系统梳理单层膜分析的核心工具：从菲涅尔公式出发，引入导纳法和特征矩阵法，并建立等效界面的物理图像。这些概念不仅是多层膜设计的数学基础，更是工程中判断膜系灵敏度、公差分配和工艺容差的理论依据。
上图展示了单层膜分析的核心逻辑链：入射介质、薄膜层、基底三层结构，以及从菲涅尔公式到特征矩阵、再到等效界面的递进关系。下面我们将逐一展开这些概念，并说明它们如何在实际膜系设计中发挥作用。
3.1 菲涅尔公式与单界面的反射透射
菲涅尔公式描述的是平面电磁波在两种均匀介质界面上的反射和透射行为。对于光学薄膜设计，我们通常关注振幅反射系数 r 和振幅透射系数 t，它们与界面两侧的光学导纳 η 直接相关。光学导纳定义为切向磁场分量与切向电场分量之比，对于非磁性材料，η = N·Y₀，其中 N = n - jk 为复折射率，Y₀ 为自由空间导纳（约 1/377 S）。在工程计算中，常将 Y₀ 归一化为 1，因此 η 在数值上等于复折射率 N。
对于光从介质 0 垂直入射到介质 1 的单一界面，振幅反射系数为：
r = (η₀ - η₁) / (η₀ + η₁)能量反射率 R = |r|²，透射率 T = 1 - R（无吸收时）。当入射角 θ₀ 不为零时，需要区分 s 偏振（TE）和 p 偏振（TM）分量。s 偏振的导纳为 η_s = N·cosθ，p 偏振的导纳为 η_p = N / cosθ。这一区分在斜入射膜系设计（如偏振分光镜、消偏振膜系）中至关重要。
3.2 单层膜的导纳分析法
导纳法是一种直观的图解与计算工具。对于单层膜，我们可以将薄膜和基底视为一个整体，用一个“等效导纳” Y 来表征。这个等效导纳是薄膜厚度、折射率以及基底导纳的函数。推导过程基于薄膜上下两个界面的电磁场边界条件，最终得到：
Y = (η₂·cosδ + j·η₁·sinδ) / (cosδ + j·(η₂/η₁)·sinδ)其中 δ = 2π·N₁·d₁·cosθ₁ / λ 是薄膜的相位厚度，θ₁ 是薄膜内的折射角。得到等效导纳 Y 后，整个单层膜系统的反射系数即为 r = (η₀ - Y) / (η₀ + Y)。
导纳法的优势在于其物理图像清晰：它把复杂的多次反射干涉问题，简化为一个“虚拟界面”的反射问题。在膜系设计初期，工程师常利用导纳图（以实部为横轴、虚部为纵轴）来观察膜层参数变化时等效导纳的轨迹，从而快速判断膜系对折射率或厚度误差的敏感度。
| 分析方法 | 核心物理量 | 适用场景 | 优点 | 局限 |
|---|---|---|---|---|
| 菲涅尔公式 | r, t（单界面） | 单界面、粗略估算 | 简单、物理意义明确 | 无法处理多层干涉 |
| 导纳法 | 等效导纳 Y | 单层膜、双层膜快速评估 | 图解直观、便于灵敏度分析 | 层数多时计算繁琐 |
| 特征矩阵法 | 2×2 矩阵 M | 任意层数膜系 | 适合编程、易于级联 | 物理图像相对抽象 |
| 等效界面 | 组合导纳 Yeq | 膜系简化、概念理解 | 将多层简化为单层 | 仅对特定波长/角度严格成立 |
3.3 特征矩阵法：多层膜的统一数学框架
特征矩阵法（又称传输矩阵法）是处理任意多层膜的标准方法。每一层薄膜对应一个 2×2 的特征矩阵，该矩阵将薄膜入射界面的电场和磁场与出射界面的电场和磁场联系起来。对于第 j 层膜，其特征矩阵为：
M_j = [ cosδ_j, (j·sinδ_j)/η_j ]
[ j·η_j·sinδ_j, cosδ_j ]整个膜系的总特征矩阵是各层矩阵按顺序相乘：M_total = M₁ × M₂ × ... × M_m。然后，膜系与基底组合后的等效导纳为 Y = C/B，其中 [B; C] = M_total × [1; η_substrate]。最终反射率为 R = |(η₀ - Y)/(η₀ + Y)|²。
特征矩阵法的强大之处在于其通用性和计算效率。无论是介质膜、金属膜、还是含有吸收或色散的材料，只需将对应的复折射率代入 δ 和 η 即可。在工程实践中，所有主流光学薄膜设计软件（如 TFCalc、Essential Macleod、OptiLayer）的核心算法都是基于特征矩阵法。
3.4 等效界面概念与膜系简化
等效界面概念是连接“单层膜分析”与“多层膜设计直觉”的桥梁。它的核心思想是：任意一个多层膜系，在特定波长和入射角下，都可以被看作一个具有某个等效复导纳 Yeq 的单一界面。这个等效导纳由整个膜系的特征矩阵决定。因此，我们可以把复杂的多层膜“压缩”成一个虚拟的基底，从而用单层膜的公式来理解其反射特性。
这一概念在以下场景中尤为实用：
- 减反射膜设计：对于单层 MgF₂ 增透膜（n≈1.38，λ/4 厚度），其等效导纳 Y = n₁²/n_sub。当 n₁² = n₀·n_sub 时，反射率为零。这就是经典单层增透膜的“零反射条件”。
- 高反射膜设计：对于 λ/4 堆栈（交替高低折射率），等效导纳随层数增加而趋向于 (n_H/n_L)^(2p)·n_sub，从而获得极高的反射率。
- 膜系公差分析：当我们需要评估某一层厚度偏差对整体光谱的影响时，可以将该层视为“扰动”，利用等效导纳的导数来快速估算灵敏度，而不必重新计算整个膜系。
3.5 案例拆解：单层 MgF₂ 增透膜的导纳与矩阵分析
我们以一个最经典的工程案例——在玻璃基底（n=1.52）上镀制单层 MgF₂（n=1.38）增透膜——来对比上述三种方法。
背景与目标：设计波长 λ₀ = 550nm，要求垂直入射时反射率尽可能低。MgF₂ 膜厚设计为 λ/4，即光学厚度 n₁d₁ = λ₀/4 = 137.5nm，物理厚度 d₁ = 137.5/1.38 ≈ 99.6nm。
方法一：菲涅尔公式（仅考虑两个界面，忽略多次反射）
界面1（空气→MgF₂）：r₁ = (1.0 - 1.38)/(1.0 + 1.38) = -0.1597
界面2（MgF₂→玻璃）：r₂ = (1.38 - 1.52)/(1.38 + 1.52) = -0.0483
若忽略干涉，总反射率 ≈ |r₁|² + |r₂|² ≈ 0.0255 + 0.0023 = 0.0278（即 2.78%）这个结果高估了实际反射率，因为它没有考虑两束反射光之间的相位关系。
方法二：导纳法
相位厚度 δ = 2π·n₁·d₁/λ₀ = π/2（因为 n₁d₁ = λ₀/4）
等效导纳 Y = (η₂·cosδ + j·η₁·sinδ) / (cosδ + j·(η₂/η₁)·sinδ)
代入 cosδ=0, sinδ=1：
Y = (j·η₁) / (j·(η₂/η₁)) = η₁² / η₂ = 1.38² / 1.52 = 1.2526
反射率 R = |(1.0 - 1.2526)/(1.0 + 1.2526)|² = |(-0.2526)/2.2526|² = 0.0126（即 1.26%）方法三：特征矩阵法
M₁ = [ cos(π/2), j·sin(π/2)/1.38 ] = [ 0, j/1.38 ]
[ j·1.38·sin(π/2), cos(π/2) ] [ j·1.38, 0 ]
基底向量 [1; η_sub] = [1; 1.52]
[B; C] = M₁ × [1; 1.52] = [ j·1.52/1.38 ; j·1.38 ]
等效导纳 Y = C/B = (j·1.38) / (j·1.52/1.38) = 1.38²/1.52 = 1.2526
R = |(1 - 1.2526)/(1 + 1.2526)|² = 1.26%结果对比与复盘：导纳法和特征矩阵法得到完全一致的结果（1.26%），而忽略干涉的菲涅尔近似给出 2.78%，误差超过一倍。这说明在 λ/4 膜厚下，多次反射干涉是决定反射率的主导因素。实际镀制中，由于 MgF₂ 膜层可能存在的柱状结构导致折射率略低（如 1.37），实测反射率通常在 1.3%-1.5% 之间，与理论值吻合良好。若想进一步降低反射率，则需要采用双层或多层膜系，例如 MgF₂ + ZrO₂ 双层增透膜，可将残余反射率降至 0.5% 以下。
3.6 常见误区与设计复盘
在单层膜分析中，工程师常犯的错误包括：
- 混淆光学厚度与物理厚度：相位厚度 δ 使用的是光学厚度（n·d），而非物理厚度。在镀膜机监控中，晶振监控的是物理厚度，而光学监控（如光电极值法）监控的是光学厚度。两者换算必须考虑膜层折射率的色散。
- 忽略折射率色散：MgF₂ 在 400nm 处折射率约为 1.39，在 700nm 处约为 1.37。若用 550nm 的折射率计算全波段光谱，会导致短波和长波处的反射率计算偏差。
- 将等效导纳视为常数：如前所述，等效导纳是波长和角度的函数。在变角度光谱拟合中，必须逐点计算。
- 对吸收膜层使用实数导纳：金属膜（如 Ag、Al）在可见光波段 k 值很大，必须使用复折射率。否则计算出的反射率和透射率会严重偏离实测值。