二、光学基础回顾：光的干涉、薄膜光学的基本方程（特征矩阵法）、反射率与透射率的计算基础

窄带滤光膜设计的本质，是在微观尺度上精确操控光的相位与振幅。无论是设计一个中心波长 532nm、半宽度 0.3nm 的激光线滤光片，还是一个用于荧光显微成像的多带通滤光片，其物理基础都离不开光的干涉现象与薄膜光学中的特征矩阵计算。本章将系统梳理干涉的物理图像、薄膜特征矩阵的推导逻辑，以及反射率/透射率的工程计算方法，为后续的膜系设计、优化与制备工艺奠定数学与物理基础。

2.1 光的干涉：薄膜光学的基础物理图像

核心原理：多光束干涉的叠加逻辑

当一束单色平面波以角度 θ₀ 入射到一层厚度为 d、折射率为 n₁ 的介质膜上时，光会在薄膜的上表面和下表面分别发生反射与透射。透射光在薄膜内部传播后，又会在下表面反射，再次到达上表面并透射出去。这一过程会无限重复，形成一系列振幅递减、相位依次延迟的反射光束和透射光束。这些光束在远场发生相干叠加，其合成振幅取决于各光束之间的相位差。

相邻两束反射光的光程差为 Δ = 2·n₁·d·cosθ₁，其中 θ₁ 是光在薄膜内的折射角。对应的相位差为 δ = (2π/λ)·Δ = (2π/λ)·2·n₁·d·cosθ₁。当 δ = 2mπ（m 为整数）时，反射光干涉增强；当 δ = (2m+1)π 时，反射光干涉相消。对于窄带滤光膜，正是利用这种干涉增强/相消效应，在特定波长处实现高反射或高透射。

判断方法：干涉增强/相消的快速判断

• 半波损失判断：当光从折射率较小的介质（光疏）射向折射率较大的介质（光密）时，反射光发生 π 相位突变，等效于额外增加 λ/2 的光程。反之，从光密到光疏反射时无相位突变。
• 光学厚度匹配：对于 1/4 波长堆栈（QWOT），每层的光学厚度 n·d = λ₀/4，此时相邻界面反射光相位差为 π，叠加后反射率最高。
• 偏振态影响：斜入射时，TE 波（s 偏振）和 TM 波（p 偏振）的有效导纳不同，导致干涉条件出现差异。窄带滤光片在斜入射时会出现中心波长漂移和偏振分离，这是设计时必须考虑的。

案例拆解：单层增透膜的干涉分析

以 K9 玻璃基底（n_sub = 1.52）上镀制单层 MgF₂ 增透膜（n₁ = 1.38）为例。目标是在中心波长 550nm 处实现零反射。根据干涉相消条件，膜层的光学厚度应满足 n₁·d = λ₀/4 = 137.5nm，即物理厚度 d = 137.5/1.38 ≈ 99.6nm。此时，上表面反射光（n₀=1 → n₁=1.38，有半波损失）与下表面反射光（n₁=1.38 → n_sub=1.52，有半波损失）的相位差为：δ = 2π·(n₁·d·2)/λ₀ = π，加上两次半波损失共 2π，等效相位差为 π，干涉相消。实际计算得到剩余反射率约为 1.3%，而未镀膜时反射率约为 4.2%。

常见误区与复盘

2.2 薄膜光学的基本方程：特征矩阵法

核心原理：从 Maxwell 方程组到矩阵连乘

特征矩阵法的核心思想是：将每一层薄膜视为一个二端口网络，其输入输出关系由该层的特征矩阵 Mⱼ 描述。对于第 j 层膜，其折射率为 nⱼ，物理厚度为 dⱼ，光在其中的折射角为 θⱼ，则相位厚度 δⱼ = (2π/λ)·nⱼ·dⱼ·cosθⱼ。该层的特征矩阵为：

Mⱼ = [cosδⱼ          i·sinδⱼ/ηⱼ]
      [i·ηⱼ·sinδⱼ       cosδⱼ    ]

其中 ηⱼ 为第 j 层的光学导纳。对于 TE 波（s 偏振），ηⱼ = nⱼ·cosθⱼ；对于 TM 波（p 偏振），ηⱼ = nⱼ/cosθⱼ。整个膜系的特征矩阵为各层矩阵的连乘：M_total = M₁ × M₂ × ... × M_N。基底和入射介质的导纳分别为 η_sub 和 η₀。

膜系与基底组合后的等效导纳定义为 Y = C/B，其中 [B; C] = M_total × [1; η_sub]。反射系数 r = (η₀ - Y)/(η₀ + Y)，反射率 R = |r|²。透射率 T = (4·η₀·η_sub)/(|η₀·B + C|²)。这一方法适用于任意层数、任意折射率分布的膜系，是膜系设计软件（如 TFCalc、Essential Macleod）的核心算法。

操作流程：特征矩阵法的计算步骤

1. 输入参数：确定入射介质折射率 n₀、基底折射率 n_sub、各层折射率 nⱼ 和物理厚度 dⱼ、入射角 θ₀、工作波长 λ。
2. 计算折射角：根据 Snell 定律 n₀·sinθ₀ = nⱼ·sinθⱼ，逐层求出 θⱼ。
3. 计算导纳：根据偏振态计算每层的 ηⱼ 和基底的 η_sub。
4. 计算相位厚度：δⱼ = (2π/λ)·nⱼ·dⱼ·cosθⱼ。
5. 构建特征矩阵：写出每层的 Mⱼ，然后按顺序连乘得到 M_total。
6. 求解等效导纳：Y = C/B。
7. 计算 R 和 T：代入反射系数和透射系数公式。

案例拆解：三层减反射膜的特征矩阵计算

设计一个三层减反射膜，用于 400-700nm 宽光谱增透。膜系结构为：基底 K9（n=1.52）/ M₂O₃（n=1.65, d=120nm）/ MgF₂（n=1.38, d=90nm）/ 空气。在 λ=550nm 处，垂直入射，TE 波。计算过程如下：

δ₁ = 2π/550 × 1.65 × 120 = 2.26 rad
δ₂ = 2π/550 × 1.38 × 90 = 1.42 rad
η₀ = 1.0, η₁ = 1.65, η₂ = 1.38, η_sub = 1.52

M₁ = [cos(2.26)    i·sin(2.26)/1.65]
     [i·1.65·sin(2.26)   cos(2.26) ]

M₂ = [cos(1.42)    i·sin(1.42)/1.38]
     [i·1.38·sin(1.42)   cos(1.42) ]

M_total = M₁ × M₂
B = M_total[0][0]×1 + M_total[0][1]×1.52
C = M_total[1][0]×1 + M_total[1][1]×1.52
Y = C/B
r = (1.0 - Y)/(1.0 + Y)
R = |r|²

计算得到 R ≈ 0.8%，相比单层膜进一步降低。该案例展示了特征矩阵法如何通过矩阵连乘快速评估多层膜的光学性能。

常见误区与复盘

2.3 反射率与透射率的计算基础

核心原理：能量守恒与边界条件

反射率和透射率是描述膜系光学性能的基本参数。根据能量守恒定律，R + T + A = 1，其中 A 为吸收率。对于无吸收介质膜（消光系数 k=0），A=0，此时 R + T = 1。但在实际镀膜材料中，尤其是金属膜或含有缺陷的介质膜，吸收不可忽略。

透射率的计算需要特别注意：T = (η_sub/η₀)·|t|²，其中 t 为透射系数，t = 2·η₀/(η₀·B + C)。因此，T = (4·η₀·η_sub)/|η₀·B + C|²。这个公式适用于无吸收膜系。对于有吸收的膜系，需要引入复折射率 N = n - ik，此时导纳变为复数，计算更为复杂。

判断方法：窄带滤光片的关键性能指标

• 中心波长（CWL）：透射率最大值对应的波长，通常要求与设计值偏差小于 ±0.1%。
• 峰值透射率（T_max）：中心波长处的透射率，窄带滤光片通常要求 T_max > 90%。
• 半宽度（FWHM）：透射率下降到峰值一半时的波长范围，窄带滤光片 FWHM 通常为 0.1-5nm。
• 截止深度（OD值）：截止带内透射率的最小值，OD = -log₁₀(T_min)，例如 OD6 表示 T_min = 10⁻⁶。

案例拆解：法布里-珀罗窄带滤光片的性能计算

一个典型的法布里-珀罗（F-P）窄带滤光片结构为：基底 / (HL)^p H (LH)^p / 空气，其中 H 为高折射率层（如 Ta₂O₅，n=2.1），L 为低折射率层（如 SiO₂，n=1.46），p 为周期数。以 p=4 为例，设计中心波长 532nm，垂直入射。通过特征矩阵法计算得到：

该案例表明，通过调整周期数 p 和间隔层厚度，可以精确控制 FWHM 和截止深度。p 越大，FWHM 越窄，但峰值透射率会下降，需要权衡。

常见误区与复盘

2.4 本章总结与工程应用要点

特征矩阵法是薄膜光学设计的基石，它将复杂的电磁场问题转化为矩阵运算，使得多层膜的光谱计算变得高效且精确。在实际工程中，设计师需要熟练掌握以下要点：

• 干涉条件的物理直觉：快速判断膜层厚度与波长的关系，预估光谱的基本形状。
• 特征矩阵的编程实现：能够用 Python、MATLAB 或商业软件快速计算任意膜系的光谱。
• 偏振与斜入射的处理：在非垂直入射场景（如 45° 入射的滤光片）中，必须分别计算 TE/TM 波。
• 误差分析与工艺容差：在设计阶段就考虑镀膜误差的影响，确保批量生产的良率。