从理论到仿真用COMSOL把波导BIC从“听起来很玄”变成“看得见摸得着”做光子学仿真的人多少都听过BIC连续谱束缚态Bound state in the continuum这个词。它听起来像是量子力学里的概念实际上在光学、声学、水波等领域都成立。简单说普通的束缚态能量在连续谱之外而BIC特殊在它藏在连续谱里面却能保持完全不辐射、不泄漏品质因子理论上无穷大。这个特性在激光器、传感器、非线性光学、滤波器和慢光器件里都有很大的想象空间。我最早接触这个概念是在做波导耦合器的时候当时想找一种“高Q但不依赖微环”的结构来做片上滤波。翻文献发现光栅波导结构里的对称保护型BIC是个很好的切入点但想复现文献里的结果手动调参数根本不现实。后来我用COMSOL Multiphysics的波动光学模块完整做了仿真把BIC从“理论图上的一条线”变成了实实在在的Q值、场图和色散数据。这篇文章就是把我自己在COMSOL里探索波导BIC的全过程整理出来包括物理基础、几何设计、边界条件的细节、参数扫描的思路以及我踩过的坑。适合有一定COMSOL基础、正在做光子晶体或集成光学仿真的人参考完全是实操向的分享。1. 为什么BIC值得做仿真物理背景与选题思路1.1 先搞清楚BIC到底是什么BIC这个概念最早是1929年von Neumann和Wigner在量子力学框架下提出的但在光子学里真正火起来是最近十几年的事。一个直观的理解方式是这样的想象一栋楼里有一条完全封闭的走廊束缚态它不跟外界连通所以人只能在这条走廊里活动。大部分束缚态对应的能量低于楼外的世界连续谱所以“困住”很好理解。但BIC这种束缚态它的能量在楼外世界的能量范围内理论上人可以直接走出去但因为它和外界没有任何耦合通道所以依然一动不动被困住。在波导结构里BIC的物理表现就是模式处于辐射连续谱的能级范围内但辐射通道的耦合系数为零。文献里最常见的做法是用对称保护型BIC即利用结构的几何对称性让某个模式的辐射场在远场完全相消干涉从而实现不辐射。这种BIC不需要精细调参只要对称性没被破坏它就在某个k点波矢点严格存在。从仿真角度看BIC意味着什么意味着在该点辐射Q值趋向无穷大损耗为零。但COMSOL里没法直接仿真“无穷大”所以实际操作时会在BIC附近取参数点观察Q值的增长趋势或者用微小扰动打破对称性得到“准BIC”quasi-BIC——这才是实际器件里最常用的形态。1.2 为什么选COMSOL而不是其他工具做BIC仿真有人用Lumerical的FDTD有人用RCWA也有人用Python自己写本征模式展开。我最终选择COMSOL主要基于几个方面的考虑。FDTD方法做时域仿真天然适合宽带响应和远场辐射的计算但BIC的Q值动辄上万甚至百万时域仿真里要等场衰减到足够小需要极长的仿真时间而且数字色散误差在低群速度区域很让人头疼。RCWA在周期性结构里效率很高但是它本质上是傅里叶展开对高折射率对比度、非矩形截面结构的收敛速度慢处理波导这种衬底辐射更麻烦。COMSOL做BIC使用的是特征频率或特征值求解Eigenfrequency/Eigenvalue study本质上是频域本征问题。它的优势有两个一个是可以直接求解复特征频率实部对应谐振频率虚部对应辐射损耗Q值直接从虚部算出来另一个是网格和边界条件的灵活性高可以逐层做局部加密处理任意几何截面。对于波导BIC这种“需要精细控制辐射通道”的问题这个特性非常关键。1.3 仿真目标怎么定复现一个对称保护型BIC拿到题目之后我没急着建模先把目标拆成了三个子问题能不能在文献中找到可参考的波导BIC结构能不能用COMSOL复现其基本特征色散关系、Q值、场分布能不能通过参数扫描找到准BIC的Q值变化规律。最终我选了一个在文献里非常常见的结构悬浮薄膜光栅波导suspended membrane grating waveguide。它的结构很简单一层高折射率薄膜比如SiN或硅上面刻蚀出周期性光栅两侧是空气。这种结构在中心对称位置存在对称保护型BIC模式在k 0附近完全不向外辐射。选它的原因也很实际几何简单网格好做对称边界条件直接可用而且文献数据充足便于交叉验证。2. 建模前的准备几何、材料、边界条件怎么选2.1 几何参数的设计逻辑波导BIC仿真的几何通常不是“画一个波导看它怎么传光”这么简单。因为我们关注的是周期性结构所以用的是单元模型unit cell。我选的结构参数如下波导厚度膜层220 nm这是硅光子的标准厚度光栅周期p650 nm大致对应1550 nm通信波段附近的带边光栅占空比刻蚀宽度与周期之比约50%后续会做扫描刻蚀深度120 nm即部分刻蚀保留底部平板波导层。这里有个关键点如果是严格的一维光栅可以做一个二维截面模型等效为y方向无限长如果你还想研究周期结构的二维效应比如带结构的k_x - k_y色散那就需要三维模型。我先从二维模型入手因为它的计算量小网格质量容易控制而且对称保护型BIC在二维模型里就能复现。实话说在二维里先搞懂物理再升级三维能少踩很多坑。2.2 材料参数与色散设置材料设置其实是个容易忽视但影响精度的环节。COMSOL的材料库里有常见的Si硅和Si₃N₄氮化硅但要注意光波导仿真通常需要指定折射率n和消光系数k而高压测的半导体材料折射率随波长变化明显直接用材料库默认常数会有误差。我更推荐手动设置解析色散模型。在COMSOL的光学模块里你可以把材料的折射率定义成随波长变化的插值函数数据可以从文献或从Palik手册里导入。对硅来说在1550 nm附近n大概等于3.48k在低损耗区可以设成0对SiNn大概在2.0左右。周围介质是空气n1不做色散。实际操作时可以在“全局定义”里建立解析函数然后材料属性引用这些函数。听起来复杂但现在这些色散数据COMSOL都内置好了直接选择“折射率实数”加上类型为“来自材料库”即可。只要你把仿真波长锁定在1550 nm附近用常数折射率也能得到合理结果如果做宽带扫描才需要完整色散。2.3 边界条件的关键思考PML还是周期性边界BIC仿真最核心的边界条件问题怎么处理辐射通道。对波导结构垂直方向通常是z方向上能量可以辐射到空气或衬底中所以必须用吸收边界。COMSOL里常用的是散射边界条件SBC或完美匹配层PML。SBC只对平面波完美吸收对高Q模式的近场辐射有一定反射PML虽然是“完美”的但设置不当也会有寄生反射。我的经验是做BIC仿真优先用PML并把它设置得尽量远离结构表面。PML的厚度一般取仿真波长的1/2以上比如1550 nm波长PML至少给800 nm到1 μm。更关键的是PML内部的网格需要随层数渐变否则即使吸收了也会有数值反射。水平方向是周期边界。COMSOL里有周期条件可以设置Floquet周期——对应不同的布洛赫波矢k。做BIC色散扫描时这个k就是扫描参数。要注意COMSOL里周期的设置是“边界选择周期类型”两个相对边界需要一致。如果你把模型只做了半个周期还要设置反对称/对称条件这个必须和模式的空间对称性匹配否则会漏掉模式或者多出伪模式。3. 实操过程从几何建模到BIC模式求解3.1 第一步画出几何结构我用的是COMSOL 6.0波动光学模块所有步骤在“模型向导”里选择“2D模型”物理场选“电磁波频域ewfd”研究类型选“特征频率”。几何层面注意单位统一COMSOL默认单位是米先改成nm更直观。第一步画底部波导层一个矩形宽等于周期p高为220 nm。第二步画光栅层两个矩形占空比50%的情况分别放在薄膜上方两端高度为刻蚀深度120 nm。如果占空比是50%周期650 nm那光栅齿宽就是325 nm空气槽也是325 nm。注意这里的坐标要精确建议直接在矩形属性里输入坐标而不是拖动鼠标。画完之后做布尔运算把光栅矩形和波导层联合成一个整体。COMSOL里这一步叫“并集”Union注意保留域边界不要把物理层之间的边界删掉因为后续设置材料时你还需要区分不同域。3.2 第二步材料指派与物理场设置选中底部波导层和光栅层指定为高折射率材料我用Si空气区域保持默认。物理场设置里需要注意以下几个点第一默认的“电磁波频域”求解的是亥姆霍兹方程特征频率求解会自动求解波动方程的本征值。但这个本征值对应的是复频率COMSOL默认显示的是实部单位是Hz。如果你要得到波长需要用c除以频率。第二面的边界条件默认是完美电导体这显然不对。你需要给所有外边界设置默认的“连续性”只有模型边缘设置PML。新建一个PML域在最外层加一个矩形框和模型区域之间保持适当的间隙我通常留300500 nm空气间隔然后在材料里设成“空气”再在物理场里添加一个“完美匹配层”特征选择这层矩形域。注意PML域不能和结构域重叠。第三周期边界。COMSOL里有“周期条件”设置需要选中左右两条边界并设置周期类型为“Floquet”或“自动”。这时候出现的是布洛赫矢的实部和虚部kx, ky在二维模型里只需要设置kx 0对应Γ点即k0然后后续做参数扫描时把这个kx改成扫描变量。这里有个细节COMSOL的Floquet周期条件需要给边界指定“源边界”和“目标边界”源和目标的方向必须匹配。我第一次做的时候搞错了方向结果算出来的模式都是乱掉的——后来看文档才发现COMSOL有个箭头符号标识源和目标需要确保箭头朝向一致。3.3 第三步网格划分的关键决策BIC仿真成败很大程度取决于网格尤其是周期结构里高Q模式模拟任何数值误差都会“唤醒”本应被抑制的辐射通道。我自己用的网格策略是这样的周期性结构内部最大单元尺寸设为波长的1/10到1/12比如1550 nm波长单元不超过150 nm光栅边缘和薄膜上下表面用边界层网格Boundary Layer加密厚度渐变PML区域要单独划分PML内单元沿径向逐渐变粗形成渐变网格这样吸收效果最好不要用全局“超细”网格否则计算时间爆炸。用局部加密把资源花在关键区域。网格划分之后先做一个快速测试在非BIC点比如kx偏离0的某个值求解一个模式看看收敛性和Q值是否合理。如果Q值随网格加密变化剧烈说明网格还没收敛需要继续加密。3.4 第四步特征频率求解与模式判别一切就绪后开始求解。COMSOL默认会给出研究配置里的模式数量可以自己指定我一般设为6个左右因为低阶模式才有物理意义。求解完成后先看结果的实部和虚部。在“一维绘图”里把频率实部列出来或直接在“派生值”里查看本征频率。COMSOL显示的是“频率复值”格式类似complex。需要用实部取倒数换算到波长虚部计算QQ ≈ f_real / (2 * |f_imag|)这里面有个容易混淆的点辐射损耗对应的是虚部为负值Q只跟虚部的绝对值有关。有些模式虚部很大说明损耗大是辐射模式不在关注范围内。真正要追踪的是虚部很小的那些点它们往往是BIC或准BIC。如何判断找对了模式一个简单的方法打开电场分布图2DBIC模式通常在结构内部形成高度局域化的场分布光场并没有漏到PML里去。如果你看到的场图里光在上下空间有大幅传播的样子那基本就是辐射模式。此外BIC对称性很明确——在中心对称的几何里模式要么是偶对称要么是奇对称场图会体现出明显的对称性。还可以把kx作为参数扫描观察模式的频率实部如何随kx变化。如果有一条平带群速度低或者带边处频率极值并且Q在某个kx值附近爆炸式增长基本可以断定这就是BIC点。4. 参数扫描与BIC点验证如何“逼出”准BIC4.1 参数化扫描设计在COMSOL里做参数扫描有两种常用方式一种是把kx设置成全局参数然后在研究里用“参数扫描”步骤扫kx另一种是用“辅助扫描”直接扫描几何参数。我建议先做kx扫描因为BIC的位置主要由波矢决定几何参数先固定。扫描范围从Γ点kx0开始向第一布里渊区边界移动一步然后逐渐加精细步长。一个典型的步长序列可以是0, 0.01(2π/p), 0.02, ..., 0.1然后把最接近BIC的区域细化。注意如果kx离BIC太远模式可能从束缚态过渡到漏模Q值急剧下降这时追踪模式的计算机逻辑可能会跳号——COMSOL的“特征值跟踪”功能可以帮你在扫描中持续追踪同一模式但需要在“特征频率”设置里打开“寻找接近指定值的特征值”。4.2 用Q值曲线定位准BIC理论上BIC处Q值趋近无穷大但仿真中是有限的。你观察到的趋势应该是随着kx趋近BIC点Q值持续上升比如从几百到几千再到几万。它并不会在某一个点突然变成无穷大而是表现为“尖峰型”曲线。我在自己的仿真里得到过比较漂亮的趋势从kx0.05(2π/p)开始这个模式的Q值从几百一路升到几万场图逐步向波导中央收缩表面辐射越来越小。这正是准BIC的典型表现。这里有个判断技巧如果Q值已经到几十万但kx稍微偏离BIC点Q值就掉回几千说明找到的是一个宽度极窄的准BIC——这类模式对结构扰动非常敏感实际加工难度高。如果你的目标不是理论研究而是做器件建议不要死守在严格BIC点而是稍微偏离一点取一个Q几千到几万的准BIC来使用。这就像调琴弦——完全共振最“完美”但实际演奏时总会有一点偏差。4.3 通过对称性微扰工程化Q值准BIC在实际中的重要操作通过破坏结构对称性把BIC变成可控的准BIC。做法很简单把光栅一侧的齿宽稍微改变比如让左右两齿宽度从325 nm变为320 nm和330 nm。这个不对称性打开了辐射通道Q值从无穷大变为有限大。在COMSOL里我把齿宽定义成两个参数w1和w2让它们之差是一个新的扫描参数。扫描不对称程度你会得到一条光滑的Q⁻¹关于不对称度平方的曲线——这个1/Q²规律是准BIC的经典特征文献里经常看到。做这一步的好处是你不只是在“复现”文献你实际上在设计一个可调Q值的器件。以后做传感或者滤波的时候这个调节能力就是核心。5. 常见问题与排查技巧实录5.1 Q值异常偏低怎么排查最常遇到的问题就是文献里明明说某个点Q值很高但自己的仿真结果Q只有几百。排查顺序如下第一步看网格。BIC是高局域模式场集中在介质内部网格粗了Q值会被“数字损耗”拉低。这是网格问题而不是物理问题。可以把网格局部加密再看结果如果Q成倍提高说明之前的网格不够。第二步看PML。PML层太薄或者离结构太近会有反射。试着把PML加厚、间隔拉大如果Q值变化明显就是PML吸收不透。第三步看周期边界方向。如果Floquet边界设置成波矢方向反了模式会在边界处非物理地“泄漏”Q值自然低。检查箭头方向是否一致。第四步看模式选择。COMSOL会输出很多模式你跟踪的那个模式可能根本不是BIC模式。用场图确认它是局域化场而不是延伸到PML里的辐射场。5.2 特征模式追踪失败、模式跳变扫描kx时COMSOL给出的模式排序可能会发生变化——某个模式在kx0.02是第2个在kx0.03变成了第3个。此时如果只是机械地跟踪“第2个模式”你看到的曲线就是乱的。解决办法是在特征频率研究设置里选择“寻找接近指定值的特征值”把目标频率设成你关注的模式频率附近比如200 THz约1500 nm然后勾选“模式跟踪”。有时候COMSOL还是会跳那就缩小扫描步长并把网格质量固定下来——网格在不同kx点如果变化太大也会导致模式判别不一致。5.3 网格收敛性测试怎么做这个我建议养成习惯。先做一套粗糙网格算出一个模式的Q值然后把结构区域的网格细化如果Q值变化超过5%继续加密直到Q值趋于稳定。收敛的Q值才是可信的。有一个技巧分别加密结构区和PML区看看哪个区域的加密对Q影响最大。如果PML区域加密影响巨大说明PML反射占据主导如果是结构区加密影响巨大说明模式还没有被充分解析。这种拆分排查方式能帮你精准定位问题。5.4 内存不足与计算时间优化三维仿真里这个问题尤其突出。二维模型一般没问题但如果扫参量很大几十个格点每个格点一次特征求解可能也要十几分钟到半小时。建议利用COMSOL的“参数化扫描”和“扫描时隐藏图形更新”功能关掉图形刷新能大幅提升速度。另外特征值求解器用“MUMPS”或“PARDISO”比默认的迭代求解器对这类大规模稀疏矩阵更稳。5.5 灵敏度与误差分析最后提醒一个容易被忽略的问题仿真结果对几何参数非常敏感。光栅齿宽变化几纳米Q值就可能下降一个数量级——这是准BIC的固有属性不是仿真误差。所以如果你准备做论文或者做流片一定要在仿真里加入制造误差分析扫描齿宽、刻蚀深度、倾角等参数看Q值的鲁棒性。我见过不少研究仿真做得很好看参数一扰动就崩了就是忽略了这一步。6. 进阶思考从BIC到实用器件的仿真路线6.1 把BIC概念迁移到其他物理场这次我用的是电磁波模块但BIC并不只属于光子学。COMSOL的声学模块、结构力学模块也都支持特征频率分析只要物理场方程形式允许连续谱中存在束缚态BIC就有机会出现。比如声波波导里也可以用对称保护的BIC做声学传感。这种方法论是可迁移的——核心思路都是“找连续谱里的零辐射模式”你在一个领域里学会了换物理场只是换方程系数和边界条件。6.2 BIC与远场发射、片上集成的结合准BIC模式在远场有定向辐射特性这在COMSOL里可以通过“远场计算”进一步分析。具体做法是把PML作为场提取的边界用远场变换得到角分辨辐射图。我自己做的一个有趣结果是准BIC的辐射方向对不对称微扰的方向非常敏感可以通过调节微扰方向“剪裁”辐射束。这种方向可调性在片上光互连、激光雷达里都很有用。如果你做的是片上集成后续往往还需要把波导BIC结构嵌入到完整器件里比如一端接输入波导、一端接输出波导观察透射谱中的Fano共振线型。这个在COMSOL里直接用“频域”研究扫描输入波长观察透过率即可基本框架和标准波导仿真一致。只是在端口设置时要注意边界条件选“端口”而不是简单的“散射边界”。6.3 我自己用下来最顺手的操作习惯最后分享几个我在实践里养成的操作习惯。第一所有结构参数都定义成全局参数不直接写数字。这样扫参和优化会方便得多改一个数字全模型同步更新省掉大量重复工作。第二每次求解之前先做一次“网格质量”检查。COMSOL里可以绘制网格单元质量分布红色的区域就是可能出问题的地方优先处理。第三把关键结果Q值、频率、场积分用全局评估表达式固定下来然后用“结果表格”功能批量导出扫描数据。后续用Python或Origin画图比在COMSOL里做复杂后处理更顺手。第四不要迷信默认边界条件。工程电磁仿真里的误差来源有一大半出在边界条件上。每换一个结构都要回头重新审视边界条件是否还成立。做波导BIC仿真给我最大的感受是这个概念看似很高冷但一旦用对工具、理顺边界条件它其实非常直观而且能直接量化评估。如果你正在做光子晶体、超表面或者集成光学方面的研究强烈建议拿一个简单的二维光栅波导复现一下BIC的Q值尖峰。这个从“无辐射”到“可控辐射”的过程既练操作也长物理直觉。