COMSOL模拟等离子共振结构超材料:多波段高吸收率文献复现案例
2026/6/10 2:16:07 网站建设 项目流程

COMSOL等离子共振结构超材料完美吸收体的多波段高吸收率,该案例为文献复现。

打开COMSOL时总有种在玩精密乐高的错觉,尤其是搞等离子体超材料这种需要微操的模型。今天要复现的这个多波段完美吸收体案例,核心在于金属-介质叠层纳米结构产生的局域表面等离子共振——说人话就是让特定频率的光波在微观结构里"卡住"。

先甩个建模关键代码片段:

model.component("comp1").geom("geom1").create("sq1", "Square"); model.component("comp1").geom("geom1").feature("sq1").set("size", ["Lx", "Ly"]); model.component("comp1").geom("geom1").feature("sq1").set("pos", ["-Lx/2", "-Ly/2"]);

这个方块几何用来定义金属谐振单元,Lx/Ly参数化方便后续调参。重点在于材料设置时要用到Drude色散模型:

epsilon_inf = 3.7; omega_p = 1.38e16; % 等离子体频率 gamma = 1.08e14; % 碰撞频率 epsilon = epsilon_inf - (omega_p^2)/(omega^2 + 1i*gamma*omega);

注意这里的虚部项决定了损耗机制,直接关系到吸收带宽。有次手滑把gamma设小了三个量级,结果吸收曲线像被刀切过一样尖锐,差点以为发现新物理现象...

边界条件设置是灵魂操作,上下用完美匹配层(PML)模拟开放空间,横向用Floquet周期边界。重点来了——多波段实现的关键在结构堆叠:

Layer1: 金纳米盘 (厚度30nm) Layer2: SiO2间隔层 (60nm) Layer3: 金纳米十字架 (45nm)

这种三明治结构会在可见光和近红外区激发不同模式的耦合共振。跑频域扫描时建议用参数化扫描代替连续扫频,毕竟在3.7THz和157THz附近各有个吸收峰,分段计算能省30%内存。

后处理阶段有个骚操作:在电磁场分布图中按住Ctrl键拖拽切面,实时观察不同频率下的热点位置。1.2μm波长时电场主要聚集在纳米盘边缘,切换到850nm则转移到十字架交叉点,这肉眼可见的模式切换比看论文示意图带感多了。

最后上张暴力调参记录表:

纳米盘直径(nm)吸收峰1位置(nm)吸收率
12082096.3%
14089098.1%
16095092.7%

看出规律没?尺寸每增加20nm,峰位红移约70nm,但超过临界值会引发模式杂化导致吸收率下降。这种非线性关系用COMSOL的优化模块自动调参比手动香,特别是结合全局搜索算法时。

模型验证时和原文献数据对不上?先检查网格尺寸是否小于最小特征长度的1/5,再确认是否忘了勾选"包含空间谐波"。有次因为漏选这个选项,计算结果在斜入射时完全跑偏,差点把键盘给砸了...

需要专业的网站建设服务?

联系我们获取免费的网站建设咨询和方案报价,让我们帮助您实现业务目标

立即咨询