2、光学基础回顾：光的干涉、反射、透射、吸收基本定律，以及薄膜光学中的特征矩阵方法

膜系构建的核心，本质上是对光与多层介质界面相互作用的精确控制。无论是设计增透膜、高反膜还是分光膜，都必须回归到光学的最基本规律：光在界面处的反射与透射、光在介质中的传播与吸收，以及多光束叠加产生的干涉效应。本章将系统梳理这些基础定律，并重点引入薄膜光学中最为强大的数学工具——特征矩阵方法。该方法将复杂的多层膜问题转化为矩阵连乘运算，是后续进行膜系优化设计与光谱性能分析的理论基石。

上图展示了特征矩阵方法的核心逻辑：将每一层薄膜视为一个2×2矩阵，整个膜系的总矩阵由各层矩阵按顺序连乘得到。这一方法将物理光学问题转化为线性代数问题，极大地简化了多层膜的计算与分析。下面我们将从基础定律出发，逐步推导并理解这一强大工具。

2.1 光与介质相互作用的基本定律

光在穿越不同介质界面时，会发生反射、透射和吸收。理解这些基本过程是分析任何光学薄膜的起点。对于薄膜设计而言，我们最关心的是能量如何分配，以及相位如何变化。

2.1.1 菲涅尔公式与界面反射率

当一束平面波入射到两种各向同性介质的界面时，反射光和透射光的振幅与入射角、偏振态以及两侧介质的折射率有关。这一关系由菲涅尔公式精确描述。对于垂直入射（θ=0°）这一在薄膜计算中最常见的情形，反射系数可简化为：

r = (n₀ - n₁) / (n₀ + n₁)
R = |r|² = ((n₀ - n₁) / (n₀ + n₁))²

其中 n₀ 为入射介质折射率，n₁ 为薄膜折射率。例如，在空气中（n₀=1.0）镀制一层 MgF₂（n₁=1.38）增透膜，单界面的反射率约为 2.6%。但实际增透膜的反射率远低于此值，这正是利用了多光束干涉效应，而非单纯依赖界面反射。

对于斜入射情形，需要区分 s 偏振（TE波）和 p 偏振（TM波）。两种偏振态的有效导纳不同，导致反射率也不同。在膜系设计中，斜入射时偏振分离效应是必须考虑的因素，例如在激光分光镜设计中，通常需要控制偏振相关损耗（PDL）。

2.1.2 吸收与消光系数

实际材料在透光波段之外往往存在吸收。描述吸收的物理量是消光系数 k，它与吸收系数 α 的关系为 α = 4πk/λ。在薄膜光学中，我们通常使用复折射率 N = n - jk 来统一描述介质的折射和吸收特性。当 k 不为零时，光在薄膜内的传播会产生衰减，特征矩阵中的相位厚度 δ 变为复数，导致反射率和透射率的计算更为复杂。

2.2 多光束干涉与薄膜的特征矩阵

单层薄膜的反射率并非简单地由两个界面反射率相加。光在薄膜内部会经历多次反射和透射，这些光束在无穷远处叠加，形成干涉。特征矩阵方法正是处理这种多光束干涉问题的优雅数学工具。

2.2.1 单层膜的特征矩阵推导

考虑一层厚度为 d₁、折射率为 n₁ 的薄膜，位于折射率为 n₀ 的入射介质和 n_s 的基底之间。光在薄膜上表面（界面1）和下表面（界面2）发生多次反射。利用电磁场边界条件，可以将薄膜上下表面的电场和磁场联系起来：

[E₁]   [cos δ₁      j sin δ₁ / η₁] [E₂]
[H₁] = [j η₁ sin δ₁    cos δ₁    ] [H₂]

其中 δ₁ = 2πn₁d₁cosθ₁/λ 为相位厚度，η₁ 为薄膜的有效导纳（对 s 偏振 η = n·cosθ，对 p 偏振 η = n/cosθ）。这个 2×2 矩阵就是单层膜的特征矩阵 M₁。它的物理意义在于：它将薄膜下表面的电磁场状态“传递”到上表面，包含了光在薄膜中传播的相位积累和阻抗变换信息。

2.2.2 多层膜的总特征矩阵

对于由 q 层薄膜组成的膜系，总特征矩阵为各层矩阵的连乘：

M_total = M₁ · M₂ · M₃ · ... · M_q

矩阵乘法不满足交换律，因此各层顺序必须严格对应光传播的先后次序。得到总矩阵后，膜系的等效导纳 Y = C/B，其中 B 和 C 由下式给出：

[B] = M_total · [1]
[C]              [η_s]

最终，膜系的反射系数和透射系数为：

r = (η₀B - C) / (η₀B + C)
t = 2η₀ / (η₀B + C)

反射率 R = |r|²，透射率 T = (η_s/η₀)·|t|²。对于无吸收膜系，R + T = 1 恒成立，这可以作为计算正确性的一个校验条件。

2.3 特征矩阵方法的工程应用与案例拆解

特征矩阵方法不仅是理论分析的工具，更是膜系设计软件（如 TFCalc、Essential Macleod）的核心算法。理解其物理意义，有助于工程师在设计出现问题时快速定位原因。

2.3.1 案例：单层增透膜的设计与验证

背景与目标：在 n_s = 1.52 的 K9 玻璃基底上设计单层增透膜，目标是在 550nm 中心波长处反射率尽可能低。

关键输入参数：入射介质 n₀ = 1.0（空气），薄膜材料选择 MgF₂（n₁ = 1.38），设计波长 λ₀ = 550nm。

分析过程：根据增透条件，薄膜的光学厚度应满足 n₁d₁ = λ₀/4，即相位厚度 δ₁ = π/2。此时 cos δ₁ = 0，sin δ₁ = 1，特征矩阵简化为：

M₁ = [0      j/η₁]
     [jη₁    0   ]

代入公式计算可得，当 n₁ = √(n₀·n_s) ≈ 1.23 时，反射率理论上可降为零。但 MgF₂ 的折射率 1.38 高于理想值，因此无法完全消除反射。实际计算得到 R ≈ 1.3%，相比未镀膜时的 4.2% 已显著降低。

结果与复盘：使用特征矩阵方法计算得到的反射率曲线在 550nm 处呈 V 形，偏离中心波长时反射率迅速上升。这一现象说明单层增透膜的带宽有限，仅适用于窄带场景。若需宽带增透，必须采用多层膜系设计。

2.3.2 常见误区与注意事项

• 矩阵顺序错误：总特征矩阵的连乘顺序必须从入射侧开始，依次到基底侧。颠倒顺序会导致计算结果完全错误。
• 忽略偏振态：在斜入射时，s 和 p 偏振的有效导纳不同，必须分别计算。若设计规格未明确偏振要求，通常取两偏振的平均值或按最差情况考虑。
• 吸收介质的处理：当薄膜材料有吸收时，折射率变为复数，特征矩阵元素变为复数，计算量增大。此时应使用复数的矩阵运算，并注意 R + T < 1，差值即为吸收率。
• 相位厚度与物理厚度的混淆：相位厚度 δ = 2πnd/λ 中的 n 是折射率，d 是物理厚度。设计时通常先确定光学厚度 nd，再根据实际材料折射率换算物理厚度。

2.4 知识结构总结与对比

为了帮助系统化理解，下表对比了不同分析方法的适用场景与特点：