2、光学薄膜基础：光的干涉、薄膜反射与透射的基本概念、特征矩阵法简介

光学薄膜的本质是在基底表面沉积一层或多层具有特定厚度与折射率的介质薄层，利用光在薄膜界面上的反射与透射行为，通过多光束干涉实现对光谱特性的调控。F-P型滤光膜正是基于这一原理，通过两个高反射膜层之间的谐振腔结构，形成窄带透射峰。理解光的干涉、薄膜反射与透射的基本概念，以及掌握特征矩阵法这一核心分析工具，是后续推导F-P型滤光膜带宽计算公式的前提。本节将从波动光学的基本模型出发，逐步建立薄膜光学分析的理论框架。

2.1 光的干涉：薄膜光学的基础物理机制

光在薄膜中传播时，会在上下两个界面分别发生反射和折射。当薄膜厚度与入射光波长处于同一量级时，从薄膜上下表面反射（或透射）的多束光波在空间叠加，由于光程差的存在，产生相长干涉或相消干涉，从而改变整个膜系的光谱响应。这一现象是薄膜光学器件的物理基础。

设一束波长为 λ 的单色光以角度 θ₀ 从折射率为 n₀ 的介质入射到厚度为 d、折射率为 n₁ 的薄膜上，薄膜下方为折射率为 n₂ 的基底。光在薄膜内传播时，相邻两束反射光之间的光程差为：

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

其中 θ₁ 为光在薄膜内的折射角，满足斯涅尔定律 n₀·sinθ₀ = n₁·sinθ₁。对应的相位差为：

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

当 δ = 2mπ（m 为整数）时，反射光发生相长干涉；当 δ = (2m+1)π 时，反射光发生相消干涉。对于正入射情况（θ₀ = 0），上述公式简化为 δ = 4π·n₁·d/λ。

2.2 薄膜反射与透射的基本概念

当光从一种介质传播到另一种介质时，在界面上会发生反射和透射。对于单层薄膜，其反射系数和透射系数不仅取决于上下界面的菲涅尔系数，还取决于薄膜的厚度和折射率。设 r₁、t₁ 为入射介质与薄膜界面的菲涅尔反射和透射系数，r₂、t₂ 为薄膜与基底界面的菲涅尔系数，则整个单层膜系的等效反射系数为：

r = (r₁ + r₂·e^(-i2δ)) / (1 + r₁·r₂·e^(-i2δ))

其中 δ 为薄膜的相位厚度。对应的反射率 R = |r|²，透射率 T = |t|²。对于无吸收介质，能量守恒要求 R + T = 1。

在实际薄膜设计中，我们更关心多层膜系的行为。多层膜的总反射和透射特性不能简单地用单层公式叠加，而需要借助特征矩阵法进行系统计算。下表对比了单层膜与多层膜分析方法的差异：

对比维度	单层膜	多层膜
分析方法	菲涅尔系数叠加	特征矩阵连乘
计算复杂度	解析公式直接求解	矩阵乘法，适合编程
适用场景	减反射膜、单层保护膜	F-P滤光膜、高反膜、分光膜
物理图像	双光束干涉	多光束干涉，等效导纳

2.3 特征矩阵法简介

特征矩阵法（Characteristic Matrix Method）是处理多层光学薄膜最常用的数学工具。其核心思想是将每一层薄膜等效为一个 2×2 矩阵，整个膜系的特征矩阵为各层矩阵的乘积，从而将复杂的多光束干涉问题转化为矩阵运算。

对于第 j 层薄膜（折射率 nⱼ，厚度 dⱼ，入射角 θⱼ），其特征矩阵为：

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

其中 δⱼ = 2π·nⱼ·dⱼ·cosθⱼ/λ 为该层的相位厚度，ηⱼ 为导纳（对 s 偏振 ηⱼ = nⱼ·cosθⱼ，对 p 偏振 ηⱼ = nⱼ/cosθⱼ）。

整个膜系的特征矩阵为：

M = M₁ · M₂ · ... · Mₙ = [ m₁₁  m₁₂ ]
                          [ m₂₁  m₂₂ ]

设入射介质导纳为 η₀，基底导纳为 ηₛ，则膜系的反射系数和透射系数分别为：

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

反射率 R = |r|²，透射率 T = (ηₛ/η₀)·|t|²。对于无吸收膜系，R + T = 1。

2.4 案例拆解：单层 SiO₂ 薄膜的反射率计算

为加深对特征矩阵法的理解，我们以一个实际案例进行演示。设基底为 K9 玻璃（nₛ = 1.52），在其上沉积一层 SiO₂ 薄膜（n₁ = 1.46），膜厚 d = 200 nm。计算在可见光范围（400~700 nm）正入射时的反射率变化。

初始条件：入射介质为空气（n₀ = 1.0），正入射（θ₀ = 0），因此 θ₁ = 0，cosθ₁ = 1。SiO₂ 薄膜的相位厚度为 δ = 2π·n₁·d/λ。

计算步骤：

1. 在波长 λ = 550 nm 处，δ = 2π × 1.46 × 200 / 550 ≈ 3.34 rad
2. 计算特征矩阵元素：m₁₁ = cos(3.34) ≈ -0.983，m₁₂ = i·sin(3.34)/1.46 ≈ i·0.184/1.46 ≈ i·0.126
3. 代入反射系数公式，得到 r ≈ 0.12，R ≈ 1.4%
4. 对比未镀膜玻璃的反射率 R₀ = ((1.52-1)/(1.52+1))² ≈ 4.2%

结果分析：在 λ = 550 nm 处，200 nm 的 SiO₂ 薄膜将反射率从 4.2% 降低到约 1.4%，起到了减反射作用。但这一效果随波长变化明显——在 λ = 584 nm（δ = π）处反射率达到最大值，在 λ = 292 nm（δ = 2π）处反射率回到未镀膜水平。这说明单层减反射膜的工作带宽有限，实际应用中常需多层膜系来拓宽减反射波段。

2.5 特征矩阵法在 F-P 滤光膜分析中的延伸

F-P 型滤光膜的基本结构为：反射镜1 / 腔层 / 反射镜2。其中反射镜通常由多层高/低折射率交替膜系（如 TiO₂/SiO₂）构成，腔层为具有一定光学厚度的介质层。利用特征矩阵法，我们可以将两个反射镜和腔层分别用矩阵表示，然后连乘得到整个 F-P 滤光膜的特征矩阵：

M_FP = M_mirror1 · M_cavity · M_mirror2

通过分析 M_FP 的矩阵元素，可以提取 F-P 滤光膜的关键性能参数：

• 中心波长 λ₀：由腔层光学厚度决定，满足 2·n_c·d_c·cosθ = m·λ₀（m 为干涉级次）
• 峰值透射率 T_max：取决于反射镜的反射率 R₁、R₂ 和膜系吸收损耗
• 半高全宽（FWHM）：与反射镜反射率的关系为 FWHM ≈ λ₀²/(2π·n_c·d_c·√(R₁·R₂))

下表总结了特征矩阵法在 F-P 滤光膜设计中的关键应用：

设计参数	矩阵法中的作用	对滤光膜性能的影响
反射镜层数	增加矩阵连乘次数	提高反射率，收窄带宽
腔层厚度	改变腔层矩阵的相位厚度	决定中心波长位置
入射角	影响各层导纳和相位厚度	中心波长随入射角漂移
材料折射率比	决定反射镜矩阵的反射率	折射率比越大，带宽越窄

在实际工程应用中，特征矩阵法不仅用于理论计算，还常与优化算法（如遗传算法、模拟退火）结合，自动搜索满足目标光谱响应的膜系结构。对于 F-P 型滤光膜，设计者通常需要权衡带宽、峰值透射率、截止深度和工艺可实现性等多个指标，特征矩阵法为这一多目标优化提供了高效的计算基础。