上个月帮学生复现一篇基于Fano共振的折射率传感器论文学生在COMSOL里搭了一整天模型跑出来的透射谱始终是光滑的单谷完全没有论文里那种不对称的尖锐线型。他自己也说不清问题出在哪儿是几何不对、材料参数不对还是网格不够密。这类问题我在做表面等离激元结构仿真时遇到过太多次了。这篇内容就围绕“COMSOL仿真复现Fano共振折射率传感器”展开从物理原理、建模流程、结果判读到排错思路把整个链路完整梳理一遍。适合正在复现相关文献、或者打算用COMSOL做等离激元传感仿真的朋友作为一份带坑位地图的实操参考。1. Fano共振的物理基础与传感指标先弄懂原理再动手1.1 什么是Fano共振两条路径才能干涉出不对称线型Fano共振不是一种独立的“新型共振”而是一类由离散态与连续态干涉产生的非对称线型。物理图像可以这样理解你站在路边观察车流所有车都在快车道上正常行驶——这是连续态。突然有一辆车驶入一条很窄的岔路又绕回来继续归入主路车流——这是离散态与连续态的耦合。当岔路里的车回到主路时与原来的车流发生相位差的叠加某些速度区间相消、某些区间相长整个车流分布就出现了不对称。光子学中的Fano共振也是同一个道理。在MIM波导加侧耦合stub腔的结构里主波导中传播的连续模式相当于“连续态”stub腔中支持的分立局域模式相当于“离散态”。入射光沿波导传播时一条路径是直接透过另一条路径是耦合进stub腔再耦合回来。两条路径在输出端发生干涉透射谱就呈现出典型的非对称Fano线型。如果只有一条路径主导谱线就是对称的Lorentz线型失去Fano结构。数学上Fano线型常用这个式子描述T(ε) T0 × (q ε)² / (1 ε²)其中ε (E - E0)/(Γ/2)是归一化能量E0是共振能量Γ是线宽q是Fano参数。当|q|很大时线型接近对称当q接近0时线型反转为谷状而中间值则形成不对称的“峰-谷”结构。1.2 为什么非对称线型更适合做折射率传感折射率传感器的基本逻辑很简单周围介质折射率变化导致结构的光学共振波长移动通过检测波长移动量来反推折射率变化。这里有两个关键指标灵敏度S和品质因数FOMS Δλ / ΔnFOM S / FWHM灵敏度S代表折射率每变化1个RIU共振波长移动多少纳米。FOM则进一步把线宽FWHM考虑进来因为线宽越窄波长移动越容易分辨检测极限越低。Fano共振在这里的优势非常明显。普通Lorentz共振线型对称两侧斜率较小检测峰位置时定位精度有限。Fano线型一侧非常陡峭波长略有变化透射强度就会发生剧烈变化。这就意味着两种检测模式都很好用一是直接追踪峰谷波长算灵敏度二是固定探测波长看强度变化后者在微纳传感领域更灵敏。所以文献里大量基于Fano共振的传感器设计并不是偶然而是确有物理优势。1.3 复现文献前要先定目标参数我见过不少同学打开COMSOL就开始建几何完全没想清楚自己要复现什么指标。结果仿真跑完不知道结果对不对甚至不知道自己在找什么。复现一篇Fano传感文献动手前至少要把这些定下来目标工作波长范围、待复现的灵敏度数值、谱线性质共振谷的线宽和深度。这些数值是后续验证模型是否正确的标尺。比如你选的文献说“灵敏度560 nm/RIU共振波长1550 nm附近Fano线宽约8 nm”。那你就要把模型调到这里附近跑到最后对比才有意义。如果文献中还给出了某个特定折射率下的透射谱和电场分布图那就更好了这些都能作为逐项比对的“标准答案”。复现不是把图做出来就完事而是能用物理逻辑解释每一步结果。2. 仿真前的准备工作物理场、材料模型与结构选型2.1 用哪个物理场接口波动光学还是射频在COMSOL里做这类结构通常选择波动光学模块中的“电磁波频域”接口。它和射频模块本质都在解Maxwell方程组但波动光学的默认设置更适合光学频段材料属性的定义方式、边界条件的封装如端口、散射边界条件、后处理的S参数提取都更顺手。对于二维MIM波导结构模型选择“二维”即可。因为几何在z方向可以视为无限长模式在x-y平面内传播这样二维近似就能抓住核心物理。如果你要研究z方向有限高度的真实结构可以用三维但计算量会上升一到两个数量级初期完全不建议一上来就做三维。先把二维模型跑通再考虑扩展。这里要提醒一个概念二维电磁波仿真中默认求解的是面内电场和面外磁场这正是MIM波导中gap SPP模式的主要分量。换言之二维模型天然适配这个物理场景不需要做额外简化。2.2 材料模型为什么推荐银和Drude模型金属的介电常数模型是这类仿真里最容易出问题的地方。表面等离激元结构对金属折射率的虚部非常敏感虚部决定损耗损耗直接影响共振深度和线宽。常用方案是在COMSOL中直接设置Drude模型。对银我惯用的参数是ε_Ag(ω) ε∞ - ωp² / (ω² iγω)取ε∞ 3.7光子能量形式的等效公式更直观用变量Eb hbar_const*omega/e_const表示光子能量单位取eV那么银的相对介电常数可以写成eps_Ag 3.7 - 9.1^2 / (Eb^2 i0.018Eb)这里的9.1 eV对应银的等离子体频率0.018 eV对应散射率。需要注意这套参数在近红外波段表现良好但在可见光短波区域银的介电常数实部会有偏差如果工作波长在600 nm以下建议改用Johnson Christy实验数据插值。用金做材料时Drude模型在500到800 nm范围内受带间跃迁影响很大直接套用会失真。所以除非你的工作波段早已落在近红外否则别偷懒用简易Drude。很多复现不上文献结果问题就出在材料模型太粗糙。2.3 结构选型MIM波导加侧耦合stub腔MIM波导结构是文献中复现Fano共振最简单的载体一层金属一层纳米级介质间隙再一层金属光以gap SPP模式在介质间隙中传播。然后在其中一侧金属壁上开一个垂直的stub腔形成T字形通道就能产生离散模式与连续背景的干涉。具体几何参数建议这样起步工作波段1300 nm到1800 nm中心尽量靠近1550 nm通信波段波导间隙宽度w_gap50 nmstub腔宽度w_stub50 nmstub腔深度L_stub200到260 nm这是调节共振波长的主参数金属层厚度h_metal300 nm以上确保电磁场在到达外边界前已经衰减干净这样的参数范围在二维仿真中网格量适中物理图像清晰参数扫描也不会太慢。实际复现时结构可能有双stub或者加环形腔但核心的光学机制依然是连续态与离散态的干涉建议先用最简结构跑通再逐步加复杂部件。2.4 参数命名和几何草图的细节COMSOL中所有几何参数、材料参数都可以用全局变量管理。建议一开始在“全局定义-参数”里把上述尺寸全部定义好单位写清楚例如w_gap 50[nm]。如果不写单位COMSOL默认采用国际单位制50就会被当成50米然后网格划分直接报错这是新手常见问题。草图设计时把仿真域理解为一个大的金属矩形中间挖出T形通道。这样建模型的好处是材料分配非常简单整个域默认是银T形通道区域单独设置为传感介质。后续做折射率传感参数扫描时只需要修改传感介质的折射率属性不需要动几何结构。3. COMSOL建模步骤详解从几何到频域扫描3.1 几何创建与布尔运算在COMSOL“几何”节点下用矩形工具创建三个区域波导间隙矩形横向长条、stub腔矩形纵向短条、整体金属矩形。然后做布尔差集从金属矩形中减去两个通道矩形得到挖空后的金属基底。操作顺序建议这样创建Rect_metal大小覆盖整个仿真域比如宽度3000 nm、高度1500 nm创建Rect_gap宽度2000 nm、高度w_gapy方向位置居中创建Rect_stub宽度w_stub、高度L_stub从波导间隙的顶部边界向上延伸使用“差集”操作用Rect_metal减去Rect_gap和Rect_stub得到金属区域T形通道本身不需要额外建模它自然形成介质域。布尔运算完成后检查一下几何是否有细小斑块或未闭合边界。有时差集之后会在直角连接处产生零厚度边对网格划分影响虽小但会干扰模式分析建议用“删除实体”清理干净。stub腔与主波导交界处的直角在仿真中可以保留物理上它带来局部场增强有利于Fano干涉。如果追求更贴近真实工艺也可以给直角加5 nm圆角但不要一开始就加会显著增加网格量。3.2 材料赋值与变量定义材料设置方面需要新建两个材料一个是银另一个是传感介质。银的相对介电常数用函数表达式定义传感介质直接用折射率定义。在二维模型中COMSOL允许你直接用相对介电常数εr作为材料属性也可以输入折射率n和消光系数k软件会自动转换。银的材料设置方法新建一个“空材料”在“电磁属性”中把相对介电常数设置为表达式。为了保证表达式里的omega和波长同步更新建议在全局参数里定义lambda0 1550[nm]omega 2*pi*c_const/lambda0Eb hbar_const*omega/e_const然后银的εr表达式写成3.7 - 9.1^2/(Eb^2 i*0.018*Eb)。传感介质的折射率同样用全局参数n_s控制初始值可以取1.33水溶液环境后续做传感扫描时就把n_s当作扫描变量。这里特别想提醒COMSOL的虚数单位是i不要用j。表达式写错一个小细节结果就会完全不对。3.3 端口边界与外部边界设置在MIM波导仿真中左端和右端需要设置端口边界条件用于激励入射光和提取透射/反射。在COMSOL中给左侧边界添加“端口”边界条件设置为激励端口模式类型选择波导模式。右侧边界同样设置端口但不勾选激励。如果软件版本支持建议使用端口的“边界模式分析”功能让COMSOL自动计算波导截面上的模式分布并确认该模式的主要能量集中在波导间隙里。跑模式分析时如果看到的模式场分布不在间隙中而是跑到金属表面甚至计算域外侧说明端口模式没选对这会影响后续所有透射谱结果。外部边界方面由于金属层厚度达到300 nm远大于银在近红外波段约22 nm的趋肤深度电磁场在穿过金属层到达外部边界前基本衰减完毕因此外部边界用默认的连续边界即可。如果你希望更严格可以在整个计算域外圈加一层完美匹配层PML但会增加网格量和求解时间。对于本文这个结构的参数范围我认为不加PML不影响主要结论。3.4 网格划分等离激元结构的关键命门网格划分是这种仿真中最影响结果精度的环节。表面等离激元场的特征是场在金属-介质界面处急剧衰减波动主要被约束在纳米级间隙内如果网格太粗数值色散会严重扭曲模式。我的推荐配置是波导间隙和stub腔区域内使用映射网格或自由三角形网格最大单元尺寸不超过5 nm金属-介质界面添加边界层网格第一层厚度1 nm增长因子1.2层数5到6金属远端区域网格可以放松到20到30 nm减少总网格量整体最大网格单元尺寸不要超过λ/10在1550 nm波段就是155 nm实际从紧不从松。曾经有一次我把间隙中网格从5 nm放宽到10 nm共振波长直接蓝移了差不多15 nmFano峰形也变钝了。这提醒我们网格密度对等离激元结构的影响远大于普通介质波导做任何结论之前都要先确认网格收敛性。3.5 频域扫描参数与求解器在“研究”节点中选择“频域”研究在“频率”设置里选择“波长”扫描。例如先做粗扫描波长范围range(1300[nm], 5[nm], 1800[nm])跑完后查看透射谱找到共振谷的大致位置再把扫描范围缩小步长改到0.2 nm甚至0.1 nm精细重建Fano线型。比较重要的是求解器选择。二维模型的网格数如果控制在几十万量级用COMSOL默认的直接求解器即可。如果模型较大内存吃紧可以尝试PARDISO求解器。有一点容易被忽略频域扫描时每个波长点的求解相对独立COMSOL会逐个计算。建议先用一个波长点试算确认模型能够正常求解、端口模式正确、内存占用可接受再提交全谱扫描。否则一个粗心导致整个扫描白跑几小时非常浪费时间。4. 仿真结果的判读与传感性能评估4.1 透射谱的提取与Fano线型判断频域扫描完成后在“派生值”里选择“全局计算”表达式写abs(ewfd.S21)^2得到的就是透射率T随波长的变化。将结果绘制成一维曲线正常情况下应该在某个波长位置出现一个不对称的“峰-谷”结构。怎么判断这确实是Fano线型而不是数值噪音或普通共振谷两个方法用Fano公式或不对称Lorentz函数对曲线拟合拟合出的q参数是否有明显的非零值观察谷两侧斜率一侧明显比另一侧陡形成不对称的过渡带。如果你跑出来的谱线是完全对称的单一谷说明离散态和连续态之间的干涉没有形成有效Fano结构。这时先别急着加密网格而是检查stub腔的共振模式是否被有效激发以及端口模式是否真的是gap SPP模式。4.2 电场分布验证模式的物理图像透射谱只能告诉你哪里发生了共振电场分布能告诉你为什么共振。在共振谷波长处绘制“电场模”图应该看到stub腔内出现强烈的场增强场被局域在stub末端或者stub与主波导交界区域。这是离散模式被激发的直接证据。在非共振波长处画一张电场图场应该主要沿波导间隙传播stub腔内几乎没有能量。通过对比两张图就能直观理解“离散态与连续态耦合”这个抽象概念在结构中的具体表现。很多文献的审稿人也会要求提供这样的对比图作为Fano机理的证据。如果你发现共振波长处场最强的地方在金属外表面而不是stub腔内说明端口激励的不是gap SPP模式而是某种表面辐射模式。这种情况回到端口设置重新分析模式不要继续往下做。4.3 折射率参数扫描灵敏度S的计算明确了共振结构之后就可以开始做折射率传感评估。把传感介质的折射率n_s作为扫描参数在典型生物检测范围内取1.33到1.38步长0.01。每组折射率跑一次频域精细扫描记录共振谷或共振峰的波长位置。把记录到的共振波长对折射率作图应该接近一条直线直线斜率就是灵敏度S。为了减少偶然误差建议至少取5个折射率点并用最小二乘拟合。如果结果在折射率范围内不是完全线性可以报告局部灵敏度dλ/dn但大多数文献中S写成线性的斜率。常见情况下这类结构的S可以做到500到2000 nm/RIU具体数值取决于共振模式和线宽。如果你的灵敏度远低于文献多半是stub腔的模式面积太大、对介质折射率响应减弱可以尝试缩小stub腔尺寸或改变波导宽度。4.4 品质因数FOM与强度检测模式灵敏度S高不等于检测能力强还要看线宽。窄线宽的意义在于共振峰移动1 nm时容易分辨出这个移动如果峰很宽1 nm的移动几乎看不出变化。因此FOM S/FWHM是更综合的指标。在计算FOM时FWHM需要从透射谷的谱线中量取。注意Fano线型不对称谷的两侧半高度位置可能不对称通常取谷底到两侧半极大值点之间的波长差即总半高全宽。用拟合出的Lorentz/Fano线宽也可以但要保证和S的统计口径一致。除了波长追踪法还有强度检测模式。
