光学基础：薄膜干涉原理、光学导纳与特征矩阵法在多层膜分析中的应用

滤光片截止带的设计本质上是多层介质薄膜的电磁波响应问题。要精确控制截止带的深度、陡度与波纹，必须从薄膜光学的基本物理图像出发，建立可计算的数学框架。本章围绕三个递进层次展开：首先从单界面的菲涅尔系数引出薄膜干涉的物理图像，然后引入光学导纳这一核心概念将电场与磁场联系起来，最后推导特征矩阵法，使多层膜的光学特性可以通过矩阵连乘快速求解。这三者构成了所有滤光片设计软件与理论分析的基础。

2.1 薄膜干涉的物理图像与相位条件

当一束单色平面波以角度θ₀入射到厚度为d、复折射率为N₁的薄膜上时，光在薄膜的上、下两个界面分别发生反射和透射。薄膜干涉的本质是：上表面反射光与下表面反射光（经薄膜内部传播后再次透射出上表面）之间发生相干叠加。两束反射光之间的相位差由两部分构成：一是光在薄膜中往返传播引入的光程差，二是界面反射时可能引入的附加相位跳变（当光从光疏介质射向光密介质时，反射光发生π相位突变）。

设薄膜的物理厚度为d，入射角为θ₁（薄膜内部），则相邻两束反射光的光程差为：

Δ = 2·N₁·d·cosθ₁

对应的相位差为：

δ = (2π/λ)·Δ = (2π/λ)·2·N₁·d·cosθ₁

当δ = m·π（m为整数）时，若两束反射光相位相同则反射增强，若相位相反则反射相消。对于滤光片设计，我们通常关注的是特定波长范围内反射率尽可能高（截止带）或尽可能低（通带）。例如，在截止带设计中，我们常采用四分之一波长堆栈（QW堆栈），即每层膜的光学厚度N·d = λ₀/4，此时δ = π/2，相邻界面反射光在特定波长λ₀处同相叠加，获得高反射。

2.2 光学导纳：连接电场与磁场的桥梁

光学导纳Y定义为电磁波在介质中传播时，磁场切向分量H与电场切向分量E的比值，即Y = H/E。在光学频段，由于介质的磁导率μ通常等于真空磁导率μ₀，光学导纳可以简化为：

Y = N·Y₀

其中N = n - jk为复折射率（n为折射率，k为消光系数），Y₀为真空导纳（Y₀ ≈ 1/377 S）。对于非磁性介质，光学导纳的实部与折射率成正比，虚部与消光系数相关。引入光学导纳的核心价值在于：它将电磁场的边界条件转化为导纳的递推关系，使得多层膜问题可以像电路中的阻抗传输线一样处理。

在薄膜光学中，我们定义膜层的修正导纳η：

对于s偏振（TE波）：η_s = N·cosθ
对于p偏振（TM波）：η_p = N / cosθ

其中θ为光在介质中的折射角。这一区分在斜入射滤光片设计中至关重要，因为s和p偏振的反射特性不同，会导致偏振分光效应。

利用光学导纳，单层膜与基底组合的等效导纳可以写为：

Y_combo = (Y_sub·cosδ + j·η·sinδ) / (cosδ + j·(Y_sub/η)·sinδ)

这一公式与传输线理论中的输入阻抗公式在数学上完全同构，因此我们可以借用微波工程中的成熟分析工具来理解薄膜光学问题。

2.3 特征矩阵法：多层膜分析的通用工具

特征矩阵法（又称传输矩阵法）是处理任意多层膜最系统、最通用的方法。其核心思想是：将每一层膜用一个2×2矩阵表示，该矩阵描述了膜层两侧电场和磁场切向分量的传递关系。对于第i层膜，其特征矩阵为：

Mᵢ = [ cosδᵢ ,  (j·sinδᵢ)/ηᵢ ]
      [ j·ηᵢ·sinδᵢ ,  cosδᵢ ]

整个多层膜系统的总特征矩阵为各层矩阵按光传播方向的乘积：

M_total = M₁ × M₂ × M₃ × … × Mₘ

设总矩阵为：

M_total = [ m₁₁ , m₁₂ ]
          [ m₂₁ , m₂₂ ]

则多层膜系统的反射系数和透射系数可以通过以下公式计算：

r = (η₀·m₁₁ + η₀·η_sub·m₁₂ - m₂₁ - η_sub·m₂₂) / (η₀·m₁₁ + η₀·η_sub·m₁₂ + m₂₁ + η_sub·m₂₂)
t = 2·η₀ / (η₀·m₁₁ + η₀·η_sub·m₁₂ + m₂₁ + η_sub·m₂₂)

其中η₀为入射介质的修正导纳，η_sub为基底（出射介质）的修正导纳。反射率R = |r|²，透射率T = (η_sub/η₀)·|t|²。

2.4 案例拆解：单层增透膜与双层截止带的矩阵分析

场景设定：假设我们需要在玻璃基底（n_sub = 1.52）上设计一个中心波长λ₀ = 550nm的增透膜，目标是在可见光中心波段将反射率从4.2%降至0.5%以下。同时，我们还需要评估一个简单的双层膜结构（高/低折射率组合）在截止带设计中的初步效果。

初始参数：

• 基底：N_sub = 1.52（无吸收）
• 入射介质：空气，N₀ = 1.0
• 可选膜料：MgF₂（n = 1.38）、TiO₂（n = 2.35）、SiO₂（n = 1.46）
• 设计波长：λ₀ = 550nm

第一步：单层增透膜计算

对于单层膜，最优增透条件为膜层光学厚度N₁·d₁ = λ₀/4，且膜层折射率满足n₁ = √(n₀·n_sub) = √(1.0×1.52) ≈ 1.23。但实际材料中折射率接近1.23的稳定薄膜很少，因此我们选用MgF₂（n = 1.38）作为近似。利用特征矩阵法计算：

δ = 2π·1.38·d/550nm，当d = 550/(4×1.38) ≈ 99.6nm时，δ = π/2
η₀ = 1.0，η₁ = 1.38，η_sub = 1.52
M₁ = [ cos(π/2) , j·sin(π/2)/1.38 ] = [ 0 , j/1.38 ]
     [ j·1.38·sin(π/2) , cos(π/2) ]   [ j·1.38 , 0 ]

代入反射系数公式：
r = (1.0×0 + 1.0×1.52×j/1.38 - j·1.38 - 1.52×0) / (1.0×0 + 1.0×1.52×j/1.38 + j·1.38 + 1.52×0)
  = (j·1.10 - j·1.38) / (j·1.10 + j·1.38)
  = (-j·0.28) / (j·2.48)
  = -0.113

R = |r|² = 0.0128 ≈ 1.28%

计算结果显示，单层MgF₂增透膜在550nm处可将反射率从4.2%降至约1.28%，但未达到0.5%的目标。这是因为MgF₂的折射率（1.38）与理想值（1.23）存在偏差。

第二步：双层膜结构分析

为了进一步降低反射率，我们采用双层膜结构：在基底上先镀一层λ₀/4的TiO₂（高折射率），再镀一层λ₀/4的MgF₂（低折射率）。利用特征矩阵法：

M₁ (TiO₂) = [ 0 , j/2.35 ]，M₂ (MgF₂) = [ 0 , j/1.38 ]
           [ j·2.35 , 0 ]          [ j·1.38 , 0 ]

M_total = M₁ × M₂ = [ -2.35/1.38 , 0 ] = [ -1.70 , 0 ]
                    [ 0 , -1.38/2.35 ]   [ 0 , -0.587 ]

代入反射系数公式：
r = (1.0×(-1.70) + 1.0×1.52×0 - 0 - 1.52×(-0.587)) / (1.0×(-1.70) + 1.0×1.52×0 + 0 + 1.52×(-0.587))
  = (-1.70 + 0.892) / (-1.70 - 0.892)
  = (-0.808) / (-2.592)
  = 0.312

R = |r|² = 0.097 ≈ 9.7%

这个结果出乎意料——双层λ₀/4膜堆在中心波长处反射率反而升高到9.7%。这是因为高/低折射率交替的λ₀/4堆栈在中心波长处形成的是高反射结构（类似于截止带的基本单元），而非增透结构。这提醒我们：膜层结构的选择必须与设计目标匹配，λ₀/4堆栈是截止带的基本构建单元，而增透膜需要采用折射率渐变或非λ₀/4厚度组合。

第三步：截止带初步评估

将上述双层膜结构视为一个最简单的截止带单元（H/L堆栈），我们计算其反射率随波长的变化。在λ₀ = 550nm处R = 9.7%，但在偏离中心波长的位置，反射率会下降。要获得宽截止带，需要多个H/L交替层。例如，8层（4对）H/L堆栈在中心波长处的反射率可计算为：

对于4对H/L堆栈，总矩阵为 (M_H × M_L)⁴
在中心波长处，每对矩阵为 [ -n_H/n_L , 0 ] = [ -1.70 , 0 ]
                          [ 0 , -n_L/n_H ]   [ 0 , -0.587 ]

(M_H × M_L)⁴ = [ (-1.70)⁴ , 0 ] = [ 8.35 , 0 ]
               [ 0 , (-0.587)⁴ ]  [ 0 , 0.119 ]

r = (1.0×8.35 + 1.52×0.119 - 0 - 1.52×0) / (1.0×8.35 + 1.52×0.119 + 0 + 1.52×0)
  = (8.35 + 0.181) / (8.35 + 0.181)
  = 1.0

R ≈ 100%

计算表明，仅4对H/L堆栈就可在中心波长处获得接近100%的反射率，这正是截止带设计的基础。实际设计中，我们通常需要10-20对堆栈来获得足够宽的截止带和足够的截止深度。

2.5 特征矩阵法的数值实现与精度控制

在实际工程中，特征矩阵法通常通过计算机程序实现。以下是一个简化的Python实现框架，展示了核心计算逻辑：

import numpy as np

def layer_matrix(n, d, theta, lam, polarization):
    """计算单层膜的特征矩阵"""
    # 计算折射角
    theta = np.arcsin(np.sin(theta) / n.real)
    # 修正导纳
    if polarization == 's':
        eta = n * np.cos(theta)
    else:  # p偏振
        eta = n / np.cos(theta)
    # 相位厚度
    delta = 2 * np.pi * n * d * np.cos(theta) / lam
    # 特征矩阵
    M = np.array([[np.cos(delta), 1j * np.sin(delta) / eta],
                  [1j * eta * np.sin(delta), np.cos(delta)]])
    return M

def multilayer_reflectance(layers, n0, n_sub, lam, theta=0, polarization='s'):
    """计算多层膜反射率"""
    M_total = np.eye(2)
    for n, d in layers:
        M_total = np.dot(M_total, layer_matrix(n, d, theta, lam, polarization))
    
    eta0 = n0 * np.cos(theta) if polarization == 's' else n0 / np.cos(theta)
    eta_sub = n_sub * np.cos(theta) if polarization == 's' else n_sub / np.cos(theta)
    
    r = ((eta0 * M_total[0,0] + eta0 * eta_sub * M_total[0,1] - M_total[1,0] - eta_sub * M_total[1,1]) /
         (eta0 * M_total[0,0] + eta0 * eta_sub * M_total[0,1] + M_total[1,0] + eta_sub * M_total[1,1]))
    return np.abs(r)**2

在数值实现中，需要注意以下几点以保证计算精度：

1. 波长网格密度：对于窄带滤光片，波长步长应小于带宽的1/10；对于宽带截止滤光片，建议步长不超过1nm。
2. 复数精度：当膜层数超过50层时，矩阵元素可能跨越多个数量级，建议使用双精度浮点数，并考虑采用散射矩阵法（S矩阵）替代传输矩阵法以提高数值稳定性。
3. 色散模型：实际膜料的折射率需通过柯西公式、塞尔迈耶公式或实测数据拟合获得，不能简单使用常数。
4. 入射角影响：在斜入射时，s和p偏振的修正导纳不同，需分别计算；对于非偏振光，反射率取两者的平均值。