📚 课程目录第3章多层膜理论基础

3. 多层膜理论基础：传输矩阵法(TMM)推导、导纳轨迹图、等效界面法

多层膜设计，说白了就是跟光波打交道。你叠一层膜，光会反射一部分；再叠一层，反射又变一点。怎么算得准？怎么设计得快？这就是本章要聊的核心。

我个人习惯把多层膜理论拆成三块：传输矩阵法（TMM）是数学工具，导纳轨迹图是可视化手段，等效界面法是简化思路。三块合在一起，你就能看懂大部分高反射膜的设计逻辑。

先问个问题：一束光穿过好几层膜，每层折射率不同，厚度也不同，你怎么算总反射率？

最直接的办法——把每层膜当成一个“黑盒子”，用矩阵描述它对光波的影响。然后一层层乘起来，就得到整个膜系的效果。这就是传输矩阵法的基本思想。

对于第 j 层膜，折射率为 n_j，厚度为 d_j，光波入射角为 θ_j。这层膜的特性可以用一个 2×2 矩阵表示：

M_j = [cos(δ_j)          (i sin(δ_j))/η_j]
      [i η_j sin(δ_j)     cos(δ_j)       ]

其中：

我的小习惯：写代码时我会把 δ 和 η 分开算，避免公式写成一长串。调试时也方便检查中间值。

假设膜系有 m 层，从入射介质到基底依次编号 1, 2, ..., m。整个膜系的传输矩阵就是各层矩阵的乘积：

M_total = M_1 × M_2 × ... × M_m = [m11  m12]
                                   [m21  m22]

注意顺序——先入射的膜放左边，后入射的放右边。我刚开始做时搞反过，算出来的反射率完全不对。

有了总矩阵，反射系数和透射系数就能直接算：

r = (η_0 m11 + η_0 η_s m12 - m21 - η_s m22) / (η_0 m11 + η_0 η_s m12 + m21 + η_s m22)
t = (2 η_0) / (η_0 m11 + η_0 η_s m12 + m21 + η_s m22)

其中 η₀ 是入射介质的导纳（通常是空气，η₀ ≈ 1），η_s 是基底的导纳。

反射率 R = |r|²，透射率 T = (η_s/η₀) |t|²。

  关键点：TMM 是精确解法，没有近似。只要你知道每层的 n、d 和入射角，就能算出任意膜系的 R 和 T。这是所有镀膜设计软件的底层算法。

TMM 算得准，但不够直观。你想想看，设计一个 20 层的高反膜，光看矩阵数字能看出什么？

导纳轨迹图就是解决这个问题的。它把每层膜对光波的影响画在复平面上，让你一眼看出膜系在往哪个方向走。

光学导纳 Y 定义为磁场强度与电场强度的比值。对于单层膜，从膜层上表面看进去的等效导纳 Y 会随着膜厚变化，在复平面上画出一个圆。

这个圆的圆心在实轴上，半径由膜层折射率和基底导纳决定。每增加一层膜，导纳就沿着圆走一段弧线。

画导纳轨迹图时，记住几条规则：

我曾经踩过的坑：画轨迹图时忘了区分 s 和 p 偏振的导纳。结果设计了一个在 0° 入射时反射率 99.9% 的膜，45° 入射时直接掉到 80%。后来才意识到，斜入射时 s 和 p 的导纳不同，轨迹图要分开画。

导纳轨迹图最直观的用处：

等效界面法，说白了就是把多层膜“压缩”成一个等效界面。你不需要关心每层膜内部发生了什么，只需要知道从外面看进去的等效导纳是多少。

对于任意多层膜，从最顶层上表面看进去，整个膜系可以等效为一个具有特定导纳 Y_eq 的界面。反射率直接由 Y_eq 和入射介质导纳 η₀ 决定：

R = |(η_0 - Y_eq) / (η_0 + Y_eq)|²

这个 Y_eq 怎么求？用 TMM 算出来的总矩阵元素就能得到：

Y_eq = (m21 + m22 η_s) / (m11 + m12 η_s)

我个人觉得等效界面法最有用的时候，是设计周期膜堆。比如经典的 λ/4 膜堆 (HL)^p，每对 HL 可以等效为一个“大界面”，导纳会逐渐远离 η₀。你算几对就能知道反射率大概到多少。

快速估算技巧：用 TiO₂ (n≈2.35) 和 SiO₂ (n≈1.46) 做 λ/4 膜堆，n_H/n_L ≈ 1.61。每增加一对 HL，反射率大约提升一个数量级。3 对就到 99%，5 对到 99.9%。这个估算在中心波长附近很准。

这两个工具其实是同一件事的两种表达。导纳轨迹图是“过程”，等效界面法是“结果”。

设计时我通常先用等效界面法估算层数，再用导纳轨迹图微调。比如发现轨迹在某个区域绕圈，说明那几层膜对反射率贡献不大，可以优化掉。

注意：等效界面法只对单一波长有效。如果你要设计宽带高反膜，不能只看中心波长的等效导纳，还要考虑色散和相位匹配。我做过一个 400-700nm 的宽带高反膜，中心波长 550nm 算出来 R=99.5%，结果 400nm 端只有 95%。后来加了匹配层才拉平。

嗯，到这里三种方法都聊完了。我个人建议你先把 TMM 的代码写一遍，跑几个例子看看结果。然后拿导纳轨迹图去验证，最后用等效界面法快速估算。三招齐下，多层膜设计就没那么神秘了。