平面等离子体手性纳米材料结构-comsol模型打开COMSOL时先别急着点鼠标——搞手性结构最怕的就是手滑。这玩意儿讲究的是几何不对称带来的光学不对称咱们今天要折腾的是经典的风车结构四个金属臂绕着中心转45度。别小看这旋转角度搞不好能让左旋光和右旋光打得头破血流。先在几何模块用矩形工具画个十字架参数化尺寸是关键。把臂长设为lp宽度wp厚度tp用全局参数控制lp 100e-9; % 纳米级尺寸得记得换算单位 wp 30e-9; tp 20e-9; rotation_angle [0,90,180,270]45; % 这才是手性的灵魂接下来用旋转命令处理每个臂注意别搞成中心对称。有个坑是旋转中心点设置一定要锁定在原点坐标0,0。见过有人在这里栽跟头结果做出来的结构比军训队列还对称。平面等离子体手性纳米材料结构-comsol模型材料库选金的Drude模型时别直接用默认参数。真正的老手会把碰撞频率调高20%epsilon_inf 3.7; omega_p 2.3e15; % 等离子体频率 gamma 1.2*0.07e15; % 手动加阻尼实测更接近实验数据边界条件设置要够骚气。背景场用散射边界但记得把入射光的极化方向改成椭圆极化。用参数化扫描把波长从400nm扫到800nm同时左右旋光分开算——这才是CD谱的玄机所在。网格划分别犯强迫症等离子体区域需要边界层网格。试过三次的人都知道在尖锐拐角处加三层边界层电场增强才能算准mesh1 createMesh(physics,emw); mesh1.geom(geom1).feature().create(bl1,BoundaryLayer); mesh1.geom(geom1).feature(bl1).set(thickness,0.1e-9);求解器设置记得勾选波形分解不然白算半天拿不到偏振信息。后处理时用交叉偏振分量算椭圆率角这才是手性响应的核心指标theta atan2(real(Ey),real(Ex)); % 电场矢量旋转角度 CD (abs(E_LCP).^2 - abs(E_RCP).^2)./(abs(E_LCP).^2 abs(E_RCP).^2);最后提醒模拟结果别直接跟文献对比先把结构导出STL文件用电子显微镜插件验证下几何形状。见过有人对着错误的结构调了三个月参数结果发现是旋转角度设反了——这种悲剧比咖啡洒在键盘还惨。