准BIC光栅多层膜COMSOL仿真

关键词：准BIC；光栅多层膜；有限元；传输矩阵；Goos位移

一、文章简要介绍

窄带滤光片、激光器和光学传感器的性能，都卡在谐振峰的 Q 值上：Q 越高，谱线越窄、灵敏度越好。连续域束缚态（BIC）理论上能给出无穷大的 Q，可实际结构因为加工误差或材料损耗，总会退化成准 BIC，Q 值随之打折。这篇 2022 年发表在 Scientific Reports 上的工作（DOI：10.1038/s41598-022-25945-1）给了一个务实的办法：在光栅下面再叠一层多层膜，把谐振 Q 值推到 10 的 7 次方量级，同时拿到大的 Goos-Hänchen 位移。整套分析先用 COMSOL 有限元算本征模定位 BIC，再用传输矩阵法核对反射谱，落地路径清晰。

图1 光栅、多层膜与复合结构的三维模型及对应 Q 值

二、COMSOL 建模与仿真流程

我们按论文搭了一套 COMSOL 模型。光栅层用高折射率 HfO2（n = 1.975），周期 d = 400 nm、高度 h1 = 290 nm；下面紧接一层 160 nm 高的波导层。两侧加 Floquet 周期边界，用模式分析扫投影能带，布里渊区取 2π/d。几何上引入非对称参数 δ = Δd/d0，让光栅在 x-o-y 面保持或打破面内对称，这是后面调 Fano 线型的关键旋钮。

图2 光栅的传输矩阵法分析示意（分 ε1、ε2、ε3 三部分）

多层膜部分换成 GaN（n = 2.45）与 Al2O3（n = 1.77）交替，单层厚 57.25 nm 和 82.47 nm，取两个周期。把光栅和多层膜叠起来就是复合结构。COMSOL 里同时做两件事：一是模式分析求本征频率，直接读出能带和 BIC 的位置；二是给平面波正入射算反射谱，用 Q = λpeak/Δλ 反推 Q 值。能带负责找点，反射谱负责定量，两条线互相印证。

图3 COMSOL 模式分析得到的光栅能带与高 Q BIC 电场分布

三、关键结果与发现

先看单纯光栅。δ = 0 时结构面内对称，反射谱上没有非对称 Fano 峰，这正是 BIC 的特征。把 δ 调到 0.1，光谱里冒出一根极窄的非对称 Fano 线型，反射峰落在 660.9 nm，电场高度局域在波导层里，Q 值随之抬升。这就说明几何非对称是打开准 BIC 辐射通道的开关。

图4 多层膜结构的二维投影能带

再看多层膜。单独算它的能带，同样能在 Γ 点附近找到 BIC。把光栅和多层膜叠成复合结构后，谐振峰移到 674.0 nm，反射接近全反，能带给出的位置也和它对得上，Q 值这次直接超过 10 的 7 次方。

图5 复合结构的 Q 值超过 10 的 7 次方与准 BIC 能带位置

δ 的连续调控值得说一句：δ 从 0.1 往下减小，谐振线宽很快收窄；减到 0，线宽消失，变成一个暗模，Q 值趋于无穷，这才是真正的 BIC。作者用通道干涉解释这套线型：宽通道来自 Fabry-Perot 谐振，窄通道来自导模谐振，两者叠加出 Fano 峰。

图6 不同非对称参数 δ 下的反射谱与 Fano 线型

再看利用准 BIC 拿大 Goos-Hänchen 位移。取 δ = 0.1、λ = 674.0 nm，反射相位角谱在 θ = 2.2 度附近剧烈跳变，位移角谱在同一角度出现尖峰。和传统靠透射型谐振或 Brewster 凹陷增强的 GH 位移不同，这里的最大位移正好落在反射峰上，检测起来方便得多。这种高 Q、强响应、易读出的组合，对高灵敏传感器、波分解复用器和光开关都很有用。

图7 反射相位角谱与 Goos-Hänchen 位移角谱

四、我们提供的服务

我们可提供光子晶体与超表面方向的仿真服务：COMSOL 波动光学建模、Floquet 周期边界与本征模分析、投影能带与 BIC、准 BIC 定位、反射谱与 Fano 线型拟合、Q 因子定量提取、传输矩阵法快速计算、Goos-Hänchen 位移与相位特性评估，以及窄带滤光片、传感器、光开关的方案设计与参数优化。如需相关服务，欢迎私信聊聊。

如有需要，欢迎通过公众号“320科技工作室”联系我们。