Comsol电磁波模型：金属超表面光栅，TE/TM偏振下斜入射不同衍射级反射光谱计算
打开COMSOL时总会被电磁波模块的选项搞得头大？今天咱们用金属超表面光栅模型实战演练。这种结构在超透镜和全息显示里常见，关键就在于精确计算不同偏振下的衍射效率。
先建个二维模型试试水。画个金属光栅周期阵列，这里有个偷懒诀窍——用矩形堆叠代替真实纳米结构，参数化设置周期d=500nm，槽宽w=300nm，厚度t=50nm。材料库直接调用Gold光学常数，注意要勾选"色散材料"选项，毕竟金属在可见光波段可是有强色散特性的。
// COMSOL Java API示例 model.param().set("d", "500[nm]"); model.geom("geom1").feature().create("r1", "Rectangle"); model.geom("geom1").feature("r1").set("size", new String[]{"w", "t"}); model.geom("geom1").feature("r1").set("pos", new String[]{"0", "-t"});边界条件设置是重头戏。上下边界用散射边界条件（SBC），左右设周期性条件。注意TM波需要调整磁场方向，这里容易搞反。入射端口设置有个坑——斜入射时要勾选相位延迟补偿，否则算出来的相位分布会错乱。
// 端口设置关键参数 model.physics("emw").feature("port1").set("ThetaSpec", "theta_inc"); model.physics("emw").feature("port1").set("PhaseShiftType", "user");求解器配置建议用频域扫掠，从400nm到800nm扫波长。内存不够的话，试试波矢扫描替代频域扫描。后处理里用emw.ref约耳这个隐藏变量直接提取各衍射级次效率，比手动积分省事多了。
Comsol电磁波模型：金属超表面光栅，TE/TM偏振下斜入射不同衍射级反射光谱计算
观察结果时发现个有趣现象：TE波在600nm附近出现异常反射抑制，而TM波在同样位置反射增强。这可能与表面等离激元共振有关，改天可以专门聊聊这个。注意衍射级次是否出现消失模，当入射角过大时高阶衍射会突然消失，这时候需要检查布儒斯特角条件。
最后附上验证代码的小技巧：把金属换成理想导体对比结果，反射率应该趋近100%。如果出现异常衰减，八成是材料设置或边界条件漏能量了。记住保存参数化扫描数据，下次做优化可以直接调取数据集。
折腾一整天终于搞定这个模型，顺手泡的咖啡都凉了三次。下次试试三维螺旋结构，估计又要跟网格剖分大战三百回合了...（查看完整代码请戳评论区置顶链接）
