从理论到仿真：用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，这是硅光子的标准厚度；
- 光栅周期p：650 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.48，k在低损耗区可以设成0；对SiN，n大概在2.0左右。周围介质是空气，n=1，不做色散。
实际操作时可以在“全局定义”里建立解析函数，然后材料属性引用这些函数。听起来复杂，但现在这些色散数据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域：在最外层加一个矩形框，和模型区域之间保持适当的间隙（我通常留300～500 nm空气间隔），然后在材料里设成“空气”，再在物理场里添加一个“完美匹配层”特征，选择这层矩形域。注意PML域不能和结构域重叠。
第三，周期边界。COMSOL里有“周期条件”设置，需要选中左右两条边界，并设置周期类型为“Floquet”或“自动”。这时候出现的是布洛赫矢的实部和虚部（kx, ky），在二维模型里只需要设置kx = 0（对应Γ点，即k=0），然后后续做参数扫描时把这个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。需要用实部取倒数换算到波长，虚部计算Q：
Q ≈ f_real / (2 * |f_imag|)这里面有个容易混淆的点：辐射损耗对应的是虚部为负值，Q只跟虚部的绝对值有关。有些模式虚部很大，说明损耗大，是辐射模式，不在关注范围内。真正要追踪的是虚部很小的那些点，它们往往是BIC或准BIC。
如何判断找对了模式？一个简单的方法：打开电场分布图（2D），BIC模式通常在结构内部形成高度局域化的场分布，光场并没有漏到PML里去。如果你看到的场图里光在上下空间有大幅传播的样子，那基本就是辐射模式。此外，BIC对称性很明确——在中心对称的几何里，模式要么是偶对称，要么是奇对称，场图会体现出明显的对称性。
还可以把kx作为参数扫描，观察模式的频率实部如何随kx变化。如果有一条平带（群速度低）或者带边处频率极值，并且Q在某个kx值附近爆炸式增长，基本可以断定这就是BIC点。
4. 参数扫描与BIC点验证：如何“逼出”准BIC
4.1 参数化扫描设计
在COMSOL里做参数扫描有两种常用方式：一种是把kx设置成全局参数，然后在研究里用“参数扫描”步骤扫kx；另一种是用“辅助扫描”，直接扫描几何参数。我建议先做kx扫描，因为BIC的位置主要由波矢决定，几何参数先固定。
扫描范围从Γ点（kx=0）开始，向第一布里渊区边界移动一步，然后逐渐加精细步长。一个典型的步长序列可以是：0, 0.01(2π/p), 0.02, ..., 0.1，然后把最接近BIC的区域细化。注意，如果kx离BIC太远，模式可能从束缚态过渡到漏模，Q值急剧下降，这时追踪模式的计算机逻辑可能会跳号——COMSOL的“特征值跟踪”功能可以帮你在扫描中持续追踪同一模式，但需要在“特征频率”设置里打开“寻找接近指定值的特征值”。
4.2 用Q值曲线定位准BIC
理论上，BIC处Q值趋近无穷大，但仿真中是有限的。你观察到的趋势应该是：随着kx趋近BIC点，Q值持续上升，比如从几百到几千再到几万。它并不会在某一个点突然变成无穷大，而是表现为“尖峰型”曲线。
我在自己的仿真里得到过比较漂亮的趋势：从kx=0.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给出的模式排序可能会发生变化——某个模式在kx=0.02是第2个，在kx=0.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值尖峰。这个从“无辐射”到“可控辐射”的过程，既练操作，也长物理直觉。
