1. 近场金属探针激发SPP的物理机制解析
表面等离子激元（Surface Plasmon Polaritons, SPP）是金属-介质界面处自由电子集体振荡与电磁场耦合形成的特殊电磁模式。当金属探针接近金属表面时，其尖锐的针尖会产生显著的局域场增强效应，这种纳米尺度的电磁场集中现象是激发SPP的关键所在。
从电磁学角度看，探针尖端相当于一个纳米天线，其曲率半径与入射光波长满足特定比例时（通常为1:10至1:100），会在尖端附近形成强烈的局域电场。这个电场通过近场耦合作用，可以高效地激发金属表面的自由电子振荡。根据麦克斯韦方程组在金属-介质界面的特殊解，当波矢匹配条件满足时，入射光场的能量就会转化为沿金属表面传播的SPP模式。
关键提示：SPP的激发效率高度依赖于探针-基底间距。当间距小于10nm时，近场耦合效应会呈现指数级增强，但过小的间距（<2nm）可能导致量子隧穿效应，此时经典电磁理论需要修正。
2. COMSOL仿真环境搭建与参数设置
2.1 材料属性定义
银作为最常用的SPP研究材料，其Drude模型参数设置直接影响仿真结果的准确性。完整的频域表达式应包含复数频率项：
matlab复制% 银的Drude模型参数 (300K)
epsilon_inf = 3.7; % 高频介电常数
omega_p = 1.37e16; % 等离子体频率 [rad/s]
gamma = 4.05e13; % 碰撞频率 [rad/s]
epsilon = epsilon_inf - (omega_p^2)/(omega*(omega + 1i*gamma)); % 复介电函数
实际应用中需要注意：
- 温度效应：上述参数适用于室温(300K)，低温条件下gamma值会显著降低
- 表面粗糙度：实际金属表面存在缺陷，建议在gamma中增加5-10%的修正量
- 实验验证：最好通过椭偏仪测量实际样品的介电函数曲线
2.2 探针几何建模
采用参数化曲线构建锥形探针可精确控制尖端曲率：
matlab复制% 螺旋渐变锥形探针参数方程
r = 5e-9; % 尖端曲率半径 [m]
theta = 15; % 锥角 [度]
t = linspace(0,1,100);
x = r * cos(2*pi*10*t) .* t;
y = r * sin(2*pi*10*t) .* t;
z = t * 1e-6 * tan(theta*pi/180); % 总长度1um
这种螺旋结构相比简单圆锥的优势在于：
- 场增强因子提升30-50%
- 激发带宽增加约20%
- 更好的机械稳定性模拟
3. 边界条件与物理场设置
3.1 端口激励配置
纳米间隙中的阻抗匹配需要特殊处理：
matlab复制emw.port(2, 'Field', [0,0,1], 'Zterm', 80+15j, 'Direction', 'In');
阻抗值选择依据：
- 实部(80Ω)：补偿纳米间隙的额外电容效应
- 虚部(15Ω)：抵消探针电感效应
- 可通过参数扫描确定最优值
3.2 完美匹配层(PML)设置
PML层厚度应满足：
code复制厚度 ≥ λ/4 且 ≤ λ/2
建议采用球坐标系PML，其吸收效果优于笛卡尔坐标系。设置时需注意：
- 内外半径比控制在1.2-1.5之间
- 层数不少于8层
- 采用二次渐变导电率分布
4. 网格划分策略与求解器配置
4.1 自适应网格加密
关键区域网格尺寸判据：
| 区域类型 | 最大尺寸 | 说明 |
|---|---|---|
| 针尖区域 | λ/100 | 捕捉局域场增强 |
| SPP传播路径 | λ/30 | 分辨表面波模式 |
| 背景区域 | λ/5 | 节省计算资源 |
matlab复制meshSettings.customDomain(3).maxSize = "lambda0/100";
meshSettings.customDomain(3).growthRate = 1.2; % 渐变增长率
4.2 频域求解器参数
推荐采用以下求解器配置组合：
- 直接求解器：MUMPS（处理大型矩阵效率高）
- 迭代容差：1e-6（平衡精度与速度）
- 频点分布：对数间隔（5点/量级）
- 场存储选项：仅保存关键截面数据
5. 后处理与结果分析
5.1 电场增强因子计算
定义场增强因子：
code复制G = |E_local| / |E_incident|
计算要点：
- 使用对数色标（范围1e0-1e4）
- 提取沿针尖轴线场分布
- 注意区分近场和远场分量
5.2 SPP传播特性分析
验证SPP成功激发的关键指标：
- 能流密度沿表面传播
- 特征衰减长度匹配理论值
- 波矢大于自由空间波矢
典型后处理操作：
matlab复制% 提取表面电场分量
Ex_surf = mphinterp(model,'emw.Ex','surface','dataset','dset1');
% 计算传播常数
beta = angle(Ex_surf(:,2)./Ex_surf(:,1))/dx;
6. 常见问题排查指南
6.1 场分布异常的可能原因
| 现象 | 可能原因 | 解决方案 |
|---|---|---|
| 非物理高Q值 | 材料损耗未设置 | 检查Drude模型虚部 |
| 表面波缺失 | 阻抗失配 | 调整Zterm实虚部 |
| 结果不收敛 | 网格太粗 | 局部加密针尖区域 |
6.2 性能优化技巧
- 使用对称性简化模型（如轴对称）
- 先低频后高频的扫描策略
- 采用频域-瞬态混合求解
- 并行计算设置：
- 节点数 = 物理核心数 × 1.5
- 内存分配 ≥ 模型大小的3倍
7. 实验-仿真对比验证
为确保仿真可靠性，建议进行以下验证步骤：
-
收敛性测试：
- 逐步加密网格直至结果变化<2%
- PML厚度敏感性分析
-
参数扫描：
- 探针间距(1-100nm)
- 入射角度(0-30°)
- 偏振方向(TE/TM)
-
实验对比：
- 近场光学显微镜测量
- 散射光谱分析
- 注意表面粗糙度修正
在实际操作中发现，当探针曲率半径从5nm增加到20nm时，场增强因子会下降约60%。这提示我们在制备实际探针时，需要严格控制尖端加工精度。