第二章 薄膜干涉基础：光程差与相位差、单层膜反射率计算、多层膜矩阵法简介
各位同学，欢迎来到薄膜干涉的世界。说实话，这一章是整个镀膜设计的根基。你想想看，我们做带通滤光片，说白了就是在玩光波之间的「合纵连横」——让某些波长的光互相加强，让另一些互相抵消。这背后的物理原理，就是今天要讲的薄膜干涉。
2.1 光程差与相位差：干涉的「语言」
先问大家一个问题：为什么薄膜会呈现颜色？
答案很简单——因为光在薄膜上下表面反射后，走了不同的路程。这两束光再相遇时，就有了「光程差」。光程差决定了它们是「握手言和」（相长干涉）还是「互相拆台」（相消干涉）。
我个人习惯把光程差记作 Δ = 2nd cosθ。其中：
- n 是薄膜的折射率
- d 是薄膜的物理厚度
- θ 是光在薄膜内的折射角
为什么是 2 倍？因为光下去又上来，走了个来回。这个公式我建议你刻在脑子里，因为后面所有设计都绕不开它。
光程差和相位差的关系也很直接：
δ = (2π/λ) × Δ
这里的 λ 是真空中的波长。相位差每变化 2π，干涉效果就重复一次。嗯，这里要注意：当光从光疏介质射向光密介质时，反射光会有一个 π 的额外相位突变——也就是「半波损失」。我在项目中遇到过好几次，因为忘了考虑这个，结果设计出来的滤光片中心波长偏了十几纳米。后来我养成了习惯，每次算相位差都先问自己一句：「有没有半波损失？」
核心要点：
- 光程差 Δ = 2nd cosθ
- 相位差 δ = (2π/λ) × Δ
- 注意半波损失（π 相位突变）
2.2 单层膜反射率计算：从公式到直觉
单层膜是最简单的镀膜结构，但别小看它。我记得刚入行时，师傅让我算一个 MgF₂ 增透膜的反射率，我拿着公式算了半天，结果和实测差了 0.3%。后来发现是忽略了基底折射率的色散。从那以后，我对「简单」二字再也不敢掉以轻心。
单层膜反射率的计算公式长这样：
R = (r₁² + r₂² + 2r₁r₂cosδ) / (1 + r₁²r₂² + 2r₁r₂cosδ)其中：
- r₁ = (n₀ - n₁) / (n₀ + n₁) —— 空气与膜层界面的反射系数
- r₂ = (n₁ - nₛ) / (n₁ + nₛ) —— 膜层与基底界面的反射系数
- δ = (4π/λ) × n₁d₁ cosθ —— 膜层引入的相位差
你想想看，当 δ = π（即膜层光学厚度为 λ/4）时，cosδ = -1，公式会简化很多。这时候如果 r₁ 和 r₂ 大小相等、符号相反，反射率可以降到零——这就是经典的 λ/4 增透膜原理。
我的小技巧： 做单层膜设计时，我习惯先算「等效折射率」n_eff = √(n₀ × nₛ)。如果膜层折射率接近这个值，增透效果最好。这个经验公式帮我快速筛选材料，省了不少试错时间。
2.3 多层膜矩阵法简介：把复杂问题「线性化」
单层膜好算，但实际带通滤光片动辄几十层、上百层膜。这时候再用反射率公式硬算，手算到天亮也算不完。怎么办？
答案是——矩阵法。
矩阵法的核心思想很简单：每一层膜用一个 2×2 的矩阵表示，整个膜系就是这些矩阵的乘积。最后从乘积矩阵里提取反射率和透射率。
每一层膜的矩阵长这样：
M = [cosδ, (i sinδ)/η
iη sinδ, cosδ ]其中：
- δ = (2π/λ) × n d cosθ —— 相位厚度
- η 是修正导纳（对 s 偏振和 p 偏振不同）
整个膜系的特征矩阵就是：
M_total = M₁ × M₂ × M₃ × ... × M_k然后反射系数 r 和透射系数 t 可以从 M_total 的元素中直接读出。具体公式我就不展开了，但你可以记住：矩阵法把复杂的多光束干涉问题，变成了简单的矩阵乘法。这就是为什么计算机能快速优化几十层膜的设计。
曾经踩过的坑： 有一次我用矩阵法算一个 60 层的带通滤光片，结果反射率曲线怎么都对不上。查了两天，发现是矩阵相乘的顺序搞反了——光从入射介质到基底，矩阵应该按入射顺序相乘，而不是反过来。从那以后，我每次写代码都会在注释里标清楚「入射侧→基底侧」的顺序。
2.4 本章知识体系
为了帮你理清思路，我画了一张图，把这一章的核心逻辑串起来：
这张图从左到右展示了本章的脉络：从干涉的基本概念出发，分两条路走——单层膜计算和多层膜矩阵法。两条路最终汇合到同一个目标：理解并设计出符合要求的光学薄膜。
好了，这一章的内容就到这里。下一章我们会深入讨论带通滤光片的设计原理，到时候这些基础都会派上用场。记住，基础不牢，地动山摇——尤其是矩阵法，我建议你多写几遍代码练练手。