基于Comsol的光学仿真：SPR、光子晶体光纤与模式分析

comsol光学仿真 SPR comsol光纤仿真 光子晶体光纤 折射率传感器 温度传感器 rsoft光纤仿真 表面等离子体共振光纤 spr-pcf SPR 论文复现图是仿的一个spr传感器和一个三芯分束器图左原文图右仿的结果 基于SPR的光纤传感器 光子晶体光纤偏振分束器 光子晶体光纤仿真 模式分析 计算等效折射率限制损耗模式色散有效模面积深夜两点盯着屏幕里的模式场分布我第N次确认这鬼畜的倏逝波耦合终于对上了文献里的趋势。搞光纤仿真的兄弟肯定懂SPR这玩意儿调起来有多虐——金属层厚度差个5纳米灵敏度直接扑街。不过话说回来玩转表面等离子体共振SPR和光子晶体光纤PCF的组合确实真香去年复现那篇SPR-PCF温度传感器论文时发现把空气孔排列改成梅花阵型TM模的损耗峰竟然出现了双共振谷这骚操作直接让温度灵敏度飙到-3.5nm/℃。在COMSOL里建模三芯分束器时边界条件设置绝对是门玄学。有次偷懒直接用了完美电导体边界结果模式耦合强度比原文低了40%后来换成散射边界条件配合自定义端口才算救回来。贴段模式分析的脚本核心部分% 导入COMSOL模式数据 mode_data mphgetquad(model,ewfd.neff,dataset,dset1); neff real(mode_data); % 等效折射率实部 alpha 8.686*2*pi*1e6*imag(mode_data)/log(10); % 限制损耗计算 plot(wavelength*1e9, alpha,LineWidth,2);这里有个坑——COMSOL导出的neff其实是归一化频率乘以真空波数得用mphinterp函数做后处理。算限制损耗时别漏了那个2π因子去年有个学弟在这卡了两星期差点退学别问我是怎么知道的。comsol光学仿真 SPR comsol光纤仿真 光子晶体光纤 折射率传感器 温度传感器 rsoft光纤仿真 表面等离子体共振光纤 spr-pcf SPR 论文复现图是仿的一个spr传感器和一个三芯分束器图左原文图右仿的结果 基于SPR的光纤传感器 光子晶体光纤偏振分束器 光子晶体光纤仿真 模式分析 计算等效折射率限制损耗模式色散有效模面积说到RSoft它的BPM算法跑传统光纤是真快。上次用BeamPROP模块仿偏振分束器30秒就出结果但遇到SPR这种需要精确计算倏逝波的情况还是得靠COMSOL的有限元法。特别是做参数扫描时在RSoft里写个循环脚本自动改金属层厚度比COMSOL的APP开发省事多了for (thickness 40; thickness 60; thickness 2) { set_material(Au, thickness, thicknessnm); run_simulation(); export_data(SPR_curve_thickness.txt); }最近在复现那篇三芯光纤分束器的论文时发现有效模面积的计算必须用矢量模式积分。COMSOL的内置变量ewfd.E有时候会漏掉横向分量正确的姿势应该是sqrt(comp1.integrate(emw.Ex^2 emw.Ey^2 emw.Ez^2)) / max(emw.normE)调通这个公式的瞬间看着仿真结果和论文里的模式轮廓完美重合突然觉得被虐的三个月都值了——直到发现作者把工作波长写反了小数点850nm写成85nm这种事真的存在啊喂。玩光子晶体光纤仿真最重要的就是吃透模式奇偶特性。有次把PCF的六个空气孔环改成渐变直径TM模的等效折射率突然出现跳变后来发现是出现了轨道角动量模式。这种意外收获比中彩票还刺激虽然大概率会被导师骂不务正业...最后给新人提个醒仿真和实验对不上时先检查材料库的介电常数曲线。去年用默认的Drude模型仿金膜SPR结果共振波长比实测偏了130nm换成Johnson Christy的实验数据立马对齐。这行当的真理永远是——代码可以糙材料数据必须骚。