设计理论基石:单层膜反射率计算,多层膜矩阵法,导纳图解法入门

做高反射膜设计,说白了就是跟光波打交道。你得知道光在每层膜里怎么走、怎么反射、怎么干涉。我刚开始接触镀膜时,觉得理论太枯燥,直接上手试参数。结果呢?镀出来的膜反射率跟设计值差一大截。后来老老实实把这几块基石啃透了,才算是真正入了门。

单层膜反射率计算

先看最简单的单层膜。一束光从空气射到薄膜上,一部分在空气-膜界面反射,一部分进入膜层,在膜-基底界面反射。这两束反射光会干涉,合成最终的反射率。

公式其实不复杂:

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 或无穷大(取决于入射介质)。每镀一对高低折射率层,导纳轨迹就绕一圈。圈数越多,终点越靠近目标点。

下面这张图展示了导纳图解法的核心逻辑:

导纳图解法核心流程 基底导纳 η₂ 第1层弧线 第2层弧线 终点 Y 弧线规律: • 每层膜对应一段圆弧,圆心在实轴上 • 弧线长度由相位厚度 δ 决定,λ/4 对应半圆 • 弧线方向:高折射率层顺时针,低折射率层逆时针 示例:4层高反膜导纳轨迹(λ/4 堆栈) 起点 (1.52, 0) → 第1层 H → 第2层 L → 第3层 H → 第4层 L 轨迹:实轴上来回摆动,逐渐靠近原点 终点越靠近原点,反射率越高 4层后反射率约 95%,继续加层可逼近 99.9%+

三种方法对比

方法 精度 直观性 适用场景
单层膜公式 高(单层) 中 快速估算单层膜反射率
多层膜矩阵法 极高 低 精确计算任意膜系,编程实现
导纳图解法 中 高 设计初期快速评估,理解膜系行为

我的建议:新手先掌握单层膜公式和导纳图解法,建立直观感觉。等需要精确设计时,再用矩阵法编程计算。三种方法配合使用,效率最高。