COMSOL太赫兹超表面BIC与能带折叠

comsol太赫兹超表面BIC与能带折叠。超表面结构里藏着不少反直觉的物理现象特别是当能带折叠遇上BIC连续谱中的束缚态总能在仿真结果里搞出些让人挠头的惊喜。最近用COMSOL折腾太赫兹频段的超表面时发现这两个机制的耦合效应特别适合用来设计高Q值谐振器今天就拎几个实操要点出来聊聊。咱们先理清楚两个概念能带折叠说白了就是周期性结构把原本的高频模式压缩到低频段显示就像把折纸反复对折后图案叠加而BIC则是谐振模式完美避开辐射损耗的玄学状态在参数空间里像幽灵一样时隐时现。这俩货在超表面里碰头时仿真结果经常会出现诡异的双峰劈裂这时候千万别急着怀疑人生——八成是参数扫描步长设大了没抓住关键点。comsol太赫兹超表面BIC与能带折叠。建模时建议先搭个最简结构练手。比如用硅基板刻蚀周期排列的十字形金属结构COMSOL脚本里初始化模型参数可以这么写model ModelUtil.create(THzMetasurface); model.param().set(p, 50e-6, 晶格常数); model.param().set(w, 12e-6, 金属线宽); model.component().create(comp1, true);这段代码里有个坑要注意参数单位必须显式声明别以为默认单位是米COMSOL有时候会抽风用厘米当基准。几何建模阶段务必开启周期端口边界条件设置里这个开关尤其关键physics.create(flprl, PeriodicPort, geom1); physics.feature(flprl).set(sweeptype, all);这里sweeptype设为all才能同时激发TE和TM模式否则可能漏掉某些谐振模式。网格划分建议用用户控制网格配合边界层特别是金属-介质界面附近至少要布三层边界元否则计算场增强时会得到离谱的数值。跑完频域扫描后重点看S参数曲线的凹陷点。当发现某个谐振峰在参数调整时Q值突然飙升到10^4量级基本可以判定触发了BIC机制。这时候用参数化扫描配合特征频率分析能在能带图上清晰看到模式交叉点——就像下面这个能带折叠的典型特征model.study.create(std2); model.study(std2).create(param, Parametric); model.study(std2).create(eig, Eigenfrequency); model.study(std2).feature(param).set(pname, new String[]{p}); model.study(std2).feature(param).set(plist, linspace(45e-6,55e-6,50));这段脚本会扫描晶格常数从45到55微米的变化每个步长对应一个特征频率解。后处理时用矩阵图绘制频率-波矢关系如果看到原本平直的能带突然出现Z字型弯折说明能带折叠开始和BIC产生耦合。最后给个实用技巧当需要精准定位BIC点时可以祭出渐近分析法。在接近临界参数的区域改用对数坐标扫描同时开启场分布监控。当观察到谐振腔内的电磁场突然从偶极振荡转变为四极对称模式时赶紧收窄扫描范围——这时候的Q值曲线会像坐火箭似的往上蹿抓住那个转折点就是设计高性能传感器的黄金参数。