设计理论基石:单层膜反射率计算,多层膜矩阵法,导纳图解法入门
做高反射膜设计,说白了就是跟光波打交道。你得知道光在每层膜里怎么走、怎么反射、怎么干涉。我刚开始接触镀膜时,觉得理论太枯燥,直接上手试参数。结果呢?镀出来的膜反射率跟设计值差一大截。后来老老实实把这几块基石啃透了,才算是真正入了门。
单层膜反射率计算
先看最简单的单层膜。一束光从空气射到薄膜上,一部分在空气-膜界面反射,一部分进入膜层,在膜-基底界面反射。这两束反射光会干涉,合成最终的反射率。
公式其实不复杂:
R = (η₀ - Y)² / (η₀ + Y)²
其中 η₀ 是入射介质的导纳(空气≈1),Y 是膜层和基底组合的等效导纳。对于单层膜:
Y = (η₂ cos δ + i η₁ sin δ) / (cos δ + i (η₂/η₁) sin δ)
这里 η₁ 是膜层导纳,η₂ 是基底导纳,δ = 2π n d / λ 是膜层相位厚度。
关键点:当膜层光学厚度为 λ/4 时,δ = π/2,公式简化为 Y = η₁² / η₂。此时反射率只取决于 η₁、η₂ 和 η₀。
我记得刚入行时,师傅让我算一个简单的 λ/4 单层增透膜。我套公式算出来反射率 1.5%,镀出来实测 1.8%。排查半天,发现是忽略了膜层色散——折射率随波长变化。从那以后,我算反射率一定带上色散数据。
多层膜矩阵法
单层膜好算,但高反射膜通常几十层甚至上百层。一层层用手算?不现实。这时候矩阵法就派上用场了。
每层膜用一个 2×2 特征矩阵表示:
Mⱼ = [cos δⱼ, (i sin δⱼ)/ηⱼ]
[i ηⱼ sin δⱼ, cos δⱼ]
整个膜系的总矩阵就是各层矩阵的乘积:
M = M₁ × M₂ × M₃ × ... × Mₖ
然后从总矩阵元素算出等效导纳 Y = C/B,再代入反射率公式。
我的习惯:写代码时把矩阵乘法封装成函数,层数再多也不怕。我常用 Python 的 numpy 做矩阵运算,几毫秒就能算完一个 50 层膜系。
我曾经帮一个同事排查设计问题。他设计的 30 层高反膜,反射率死活达不到 99.5%。我用矩阵法逐层检查,发现第 12 层和第 13 层的相位厚度算反了——一个该是 λ/4 的用了 λ/2。改过来后,反射率直接跳到 99.8%。
导纳图解法入门
矩阵法精确,但不够直观。导纳图解法就不一样了——你可以在复平面上画导纳轨迹,一眼看出膜系的行为。
横轴是导纳实部,纵轴是虚部。每镀一层,导纳点就移动一段弧线。起点是基底的导纳,终点是膜系等效导纳。弧线的形状和方向告诉你:这层膜是增反还是增透。
注意:导纳图上的弧线不是随便画的。每层膜的弧线对应一个圆心在实轴上的圆。圆心位置和半径由膜层导纳和相位厚度决定。画错了,分析就全偏了。
我习惯在设计高反射膜时先用导纳图做快速评估。比如要设计 1064nm 的高反膜,基底是 K9 玻璃(导纳 1.52),目标反射率 99.9%。在导纳图上,你需要让终点导纳尽可能接近 0 或无穷大(取决于入射介质)。每镀一对高低折射率层,导纳轨迹就绕一圈。圈数越多,终点越靠近目标点。
下面这张图展示了导纳图解法的核心逻辑:
三种方法对比
| 方法 | 精度 | 直观性 | 适用场景 |
|---|---|---|---|
| 单层膜公式 | 高(单层) | 中 | 快速估算单层膜反射率 |
| 多层膜矩阵法 | 极高 | 低 | 精确计算任意膜系,编程实现 |
| 导纳图解法 | 中 | 高 | 设计初期快速评估,理解膜系行为 |
我的建议:新手先掌握单层膜公式和导纳图解法,建立直观感觉。等需要精确设计时,再用矩阵法编程计算。三种方法配合使用,效率最高。