Comsol光学仿真模型：纳米球/柱Mie散射多级分解

Comsol光学仿真模型：包括纳米球/柱 Mie散射多级分解

在COMSOL里折腾纳米颗粒的光学响应总让我想起小时候拆收音机的经历——表面看起来是个简单的金属疙瘩，内部却藏着复杂的电磁场舞蹈。这次咱们重点聊怎么用多级分解的手法，把纳米球和纳米柱的Mie散射掰开了揉碎了看。先来个直径100nm的金球镇楼，在可见光波段蹦跶的等离激元共振可不是吃素的。

建模第一步得学会玩转参数化扫描，毕竟尺寸效应在纳米尺度下比女朋友的脸色还敏感。在COMSOL的JavaAPI里可以这样浪：

double[] diameters = {80, 100, 120}; for (double d : diameters) { model.geom("geom1").feature().create("sph"+d, "Sphere"); model.geom("geom1").feature("sph"+d).set("r", d/2 + "nm"); }

这串代码批量生成不同直径的纳米球，比手动点鼠标快三倍咖啡时间。注意单位换算这个老六，COMSOL默认是m为单位，手滑少个nano前缀能直接送你个黑洞级别的仿真结果。

散射场的多级展开才是重头戏。在频域研究的后处理里插入这个表达式：

sqrt(abs(sum(emw.SrLmodes^2 + emw.SrMmodes^2, mode, 1, 10)))

这里mode参数控制展开的项数，1到10对应从电偶极到24极的展开过程。前阵子做银纳米柱时就发现，当波长接近400nm时，四级项的贡献突然暴走，直接把散射截面顶高了两个数量级，这货的电磁模式切换比川剧变脸还刺激。

Comsol光学仿真模型：包括纳米球/柱 Mie散射多级分解

网格划分的讲究在于表面等离子体的捕捉，这里有个邪道技巧——在材料界面处塞个边界层：

model.mesh("mesh1").feature().create("boundLyr", "BoundaryLayers"); model.mesh("mesh1").feature("boundLyr").set("boundary", "nanosphere_surf"); model.mesh("mesh1").feature("boundLyr").set("thickness", "lambda0/20");

这种操作能让电场增强区域的解析度提升五倍，代价是计算时间多喝两杯拿铁的功夫。记得把lambda0替换成工作波长，别傻乎乎直接输字符串。

当把纳米球拉长成纳米柱，散射模式立马开启精分模式。用旋转对称节省计算量时，记得在方位角模态数设置里玩俄罗斯套娃：

model.physics("emw").feature("azimMode").set("m", range(0,5));

这种多模式叠加的骚操作，能同时捕捉到偶极震荡和环形电流模式之间的量子纠缠——哦不，是经典电磁耦合。上次模拟金纳米柱阵列时，第三阶模式突然和相邻颗粒的基模共振，那散射图谱扭曲得毕加索都直呼内行。

最后提个醒：多级分解的结果别直接用默认配色方案，在MATLAB后处理里搞个彩虹渐变映射，不同级次用不同透明度叠加，这样在组会上展示时，审稿人想找茬都得先被美学暴击三连。毕竟搞仿真的不仅要会算数，还得是个合格的颜色控，你说是不？