膜厚偏差会怎样改变滤光片:Fabry–Pérot 灵敏度分析

Fabry–Pérot 腔厚扫描:峰位、线宽与制造公差

本篇接续在两面镜子之间选出一种颜色:Fabry–Pérot 窄带滤光片,沿用其中的 550 nm 对称 DBR—腔—DBR 结构。

窄带滤光片的腔层只偏差几纳米,透射峰就可能离开目标波长。这篇教程扫描 Fabry–Pérot 腔厚,量化峰位、半高全宽和 Q 因子的变化,并把光谱容差换算为膜厚公差。

本篇只改变中间 MgF₂ 腔层厚度。两面 DBR、材料折射率、入射角和波长采样全部固定,使峰位变化可以明确归因于腔厚。

Fabry–Pérot 腔内的光在两个反射界面之间多次往返
腔厚改变每次往返积累的相位,因此会移动透射峰Krishnavedala / Wikimedia CommonsCC0 1.0

腔层和镜面共同决定峰位

法向入射时,腔内往返一周的相位条件为

Φ(λ,dc)=4πncdcλ+ϕf(λ)+ϕb(λ)=2πm.\Phi(\lambda,d_c) =\frac{4\pi n_cd_c}{\lambda} +\phi_f(\lambda)+\phi_b(\lambda) =2\pi m.

式中,Φ\Phi 为一周总相位,λ\lambda 为真空波长,ncn_cdcd_c 分别为腔层折射率和物理厚度,ϕf\phi_fϕb\phi_b 为前后 DBR 的反射相位,mm 为正整数谐振级次。相位以弧度表示,λ\lambdadcd_c 使用相同长度单位。

忽略两面 DBR 的反射相位时,一级谐振可简化为

λsimple=2ncdc.\lambda_{\mathrm{simple}}=2n_cd_c.

半波关系适合估算初始腔厚,也能预测“腔越厚,峰越向长波移动”。但 DBR 的反射相位随波长变化,实际灵敏度必须由完整膜系计算。

固定镜面与光学条件

沿用 550 nm Fabry–Pérot 滤光片:

区域层序重复次数
前镜MgF₂ 99.64 nm / TiO₂ 56.12 nm4
腔层MgF₂ 199.28 nm1
后镜TiO₂ 56.12 nm / MgF₂ 99.64 nm4
基底玻璃 1 mm,非相干1
四对 DBR、MgF2 腔和四对镜像 DBR 的完整结构
图 1|灵敏度分析沿用的对称 Fabry–Pérot 结构

前镜和后镜的周期层序互为镜像。分别打开两个 Edit Group,确认 Repeat Count 都为 4。

前镜 Edit Group 对话框中的四对 MgF2 TiO2 周期
图 2|前镜的 MgF₂/TiO₂ 周期单元
后镜 Edit Group 对话框中的四对 TiO2 MgF2 周期
图 3|后镜的 TiO₂/MgF₂ 周期单元
扫描期间不要修改 DBR 的层序、厚度或重复次数。镜面变化也会改变反射相位,使结果不再只表示腔厚灵敏度。

Optics 中设置 450–700 nm、0.25 nm 步长、0° 入射、未偏振光,并启用 ReflectanceTransmittanceAbsorptance

Fabry Perot 灵敏度分析的细波长采样设置
图 4|0.25 nm 步长用于测量窄透射峰

名义腔厚 199.28 nm 时,透射峰位于 550.00 nm,峰值透射率为 95.742%,半高全宽为 3.623 nm。1 nm 步长只能在峰宽内留下三到四个采样点,不足以可靠测量半高交点。

半高全宽(full width at half maximum,FWHM)和品质因子定义为

Δλ=λRλL,Q=λpeakΔλ.\Delta\lambda=\lambda_R-\lambda_L, \qquad Q=\frac{\lambda_{\mathrm{peak}}}{\Delta\lambda}.

式中,λL\lambda_LλR\lambda_R 为透射率下降到峰值一半时的左右交点波长,Δλ\Delta\lambda 为 FWHM,λpeak\lambda_{\mathrm{peak}} 为峰值波长,QQ 为无量纲品质因子。

扫描腔层厚度

Sweep 中选择 MgF2 Cavity → Thickness,从 190 nm 扫描到 210 nm,步长 5 nm。

运行前先预测:腔厚增加时透射峰应向长波移动;若 DBR 反射相位不可忽略,实际斜率会偏离简单估算值 2nc2n_c
Sweep 页面中 190 到 210 nm 的 MgF2 腔厚扫描
图 5|MgF₂ 腔厚扫描

运行扫描后,在 Transmittance 中比较五条曲线,并在 500–600 nm 内跟踪同一个腔模。

五种 MgF2 腔厚对应的 Fabry Perot 透射峰
图 6|腔厚变化引起的透射峰漂移
腔厚半波公式预测仿真峰位峰值 TTFWHMQ
190 nm524.39 nm538.50 nm95.514%3.602 nm149.5
195 nm538.19 nm544.75 nm95.764%3.596 nm151.5
200 nm551.99 nm551.00 nm95.355%3.643 nm151.2
205 nm565.79 nm557.00 nm95.764%3.686 nm151.1
210 nm579.59 nm563.25 nm95.642%3.769 nm149.4

腔厚增加 20 nm 后,透射峰向长波移动了 24.75 nm;峰值仍约为 95.4%–95.8%,FWHM 保持在 3.60–3.77 nm。这个范围内,腔厚主要移动峰位,没有显著破坏透射峰。

MgF2 腔厚与透射峰波长的仿真结果和半波估算
图 7|半波公式给出正确方向,但高估了实际峰位灵敏度

全部扫描结果的最大能量闭合误差小于 1.0×10141.0\times10^{-14}

从峰位漂移得到膜厚公差

对五个仿真点做线性拟合,得到当前设计附近的局部灵敏度

S=dλpeakddc1.235 nmnm.S=\frac{\mathrm{d}\lambda_{\mathrm{peak}}}{\mathrm{d}d_c} \approx1.235\ \frac{\mathrm{nm}}{\mathrm{nm}}.

式中,SS 为峰位对腔厚的局部灵敏度,dλpeak\mathrm{d}\lambda_{\mathrm{peak}}ddc\mathrm{d}d_c 分别为峰位与腔厚的微小变化。腔厚变化 1 nm,峰位在本设计附近约移动 1.24 nm。

半波公式给出的 Ssimple=2nc=2.760 nm/nmS_{\mathrm{simple}}=2n_c=2.760\ \mathrm{nm}/\mathrm{nm} 高出一倍以上,原因是它忽略了 DBR 的反射相位与等效穿透深度。

若峰位允许偏差为 Δλallow\Delta\lambda_{\mathrm{allow}},可估算腔厚偏差:

ΔdcΔλallowS.|\Delta d_c|\leq\frac{\Delta\lambda_{\mathrm{allow}}}{S}.

式中,Δdc\Delta d_c 为腔厚偏差,Δλallow\Delta\lambda_{\mathrm{allow}} 为允许的峰位偏差,SS 为上述局部灵敏度。

允许峰位偏差腔厚偏差上限
±1 nm约 ±0.81 nm
±2 nm约 ±1.62 nm
±5 nm约 ±4.05 nm

这是一项确定性的单参数灵敏度分析,不是随机误差或成品率计算。改变 DBR 周期数、材料、入射角或工作波长后,需要重新扫描并拟合 SS

设计判断

设计动作主要影响重新验收
调整腔厚移动透射峰峰位和膜厚公差
增加 DBR 周期数通常使峰更窄峰值、FWHM、Q 与吸收损耗
改变入射角峰位发生角度漂移使用角度与偏振
改用色散材料光学厚度随波长变化全波段峰位与线宽

名义设计命中 550 nm 只是第一步。可制造的滤光片还要把光谱规格转换成厚度规格,并确认沉积过程能够满足它。

变化题

把腔厚扫描改为 195–203 nm、步长 1 nm,其他设置不变。重新拟合名义厚度附近的 SS,比较细扫结果与本篇五点拟合的 1.235 nm/nm 是否一致。


← 返回教程目录 · 下一篇:层数越多越好吗

Copyright © 2026 Dreapex