COMSOL金属SPP降维降损
金属表面等离子体激元（SPP）的模拟总让人又爱又恨——高局域场增强的特性是真香，但三维全波仿
真动不动就内存爆炸也是真头疼。最近在COMSOL里折腾SPP降维模型时发现，只要玩点几何骚操作，计算量
能砍半，损耗还能压得更低。这里分享几个野路子，代码直接贴出来大家随便薅。
### 把三维拍扁成二维
传统做法是在三维空间里建个金属-介质交界面，但光场沿着表面指数衰减的特性其实给了降维空
间。直接上代码片段：
```matlab
% COMSOL弱形式PDE设置
model.physics('w').feature('wfeq').set('weakExpression',...
'-epsilon*test(Ex)*Ex - test(Ey)*Ey + 1i*k0*Zs*Hz*test(Hz)');
```
这里把z方向场量用倏逝波形式隐式处理，用复波数k0和表面阻抗Zs代替完整三维计算。实测在155
0nm波段，内存占用从32G直降到8G，但TM模式的传播长度误差控制在3%以内。
### 边界条件魔改
金属损耗大头来自欧姆损耗，在二维模型中容易算得比实际偏高。解决方法是在边界条件里埋个陷
阱：
```matlab
% 自定义阻抗边界
surfaceImpedance = (1+1i)*sqrt(mu0/(epsilon0*omega*sigma));
model.physics('emw').feature('surfImp1').set('Zs', surfaceImpedance);
```
这里把经典表面阻抗公式改成含频率ω和电导率σ的动态表达式。当金属厚度小于趋肤深度时，这个
修正项能把损耗计算误差从15%压到5%以下。
### 几何参数调参玄学
波导宽度和介质层厚度存在神秘最优解。用COMSOL的全局扫参配合MATLAB脚本自动寻优：
```matlab
for w = linspace(50e-9,200e-9,20)
model.param.set('wg_width', num2str(w));
model.sol('sol1').runAll;
