企业官网建设流程全解析

1. 先说清楚“off-gamma BIC”到底算什么物理

1.1 BIC不是玄学：连续谱里的“隐身态”

我最早听到Bound State in the Continuum，第一反应是“这不就是一个永远不辐射的共振吗，跟理想模式有什么区别”。后来做超表面高Q谐振设计才发现，事情没那么简单。BIC指的是这样一个态：它的频率落在辐射连续谱内，理论上可以和所有平面波通道耦合，但因为对称性或参数配平的原因，它对外辐射的分量刚好为零，能量就老老实实待在结构里。换个说法，它是一间漏风的房子里的某个房间，外面再怎么刮风下雨，这个房间里的能量就是出不去。

在COMSOL里算BIC，核心工作其实不是“找不辐射的模式”，而是把辐射通道在数值上完整地呈现出来，然后找到辐射阻尼消失的那个参数点。BIC在倒空间中通常是一个孤立点：只有一个布洛赫波矢kx对应无限高Q。实际样品总有尺寸、入射角展宽，你不可能精确落在那个k点上，所以近几年大家都在做“Merging BIC”，也就是把多个孤立BIC通过调节几何参数合并到一起，让高Q区域从“一个尖峰”变成“一片平台”。

1.2 Γ点BIC和off-Γ BIC是两种完全不同的存在

Γ点（kx=0, ky=0）出现的BIC一般是对称性保护的。结构里有C2或C4旋转对称性，某些模式在这些对称操作下的变换性质和所有辐射通道都正交，对应模式就无法辐射。这种BIC只要结构对称性不被破坏，它就一直存在，不太需要专门调参数。

off-Γ的BIC就麻烦得多。它在kx不为零的某个位置出现，没有严格的对称性保护，纯粹靠结构的几何参数在某组特定数值下让偶然辐射通道的系数碰巧为零。这类BIC也被称为偶然BIC或参数型BIC，它的位置、存在性和Q值对几何参数极其敏感。椭圆柱的椭圆度、柱高、周期，只要微调一点点，off-Γ BIC的位置就会在k空间里移动，甚至消失。

“Merging off-gamma BIC”就是利用这种敏感性：把某个几何参数扫过去，让两个原本分开的off-Γ BIC在k空间里逐渐靠近，最后合并。合并成功之后，Q值在很大的k区间内都保持高数值。这个概念对做传感、窄带滤波、低阈值激光的人特别有用，因为实际光斑尺寸有限，角谱展宽不可避免，合并后的BIC对入射角不再苛刻。

1.3 Merging解决的真实痛点：Q值对环境误差的敏感

我见过不少刚接触这个方向的人，跑COMSOL找到某个kx处的Q峰，兴奋得不行，结果稍微改一下周期或者厚度，那个峰就消失了。这不是算法的问题，而是单一off-Γ BIC本质上就是“参数悬崖”上的平衡：存在条件非常脆弱。Merging BIC的意义在于，通过把BIC合并，让极端高Q区域从零维变成一个有一定宽度的区域。这样样品即便存在几十纳米的尺寸误差，或者激发光有一个小的角度偏差，Q值仍然能保持较高水平。

所以我建议所有打算在COMSOL里复现Merging BIC计算的人，动手之前先把这三个概念分清：对称保护BIC在Γ点、偶然BIC在off-Γ、通过参数扫描让多个偶然BIC汇合。章节后面讲的所有建模和扫描策略，都是围绕这三个对象的辨别与追踪展开的。

2. COMSOL建模之前的选型决定了后面所有计算的成败

2.1 用RF模块还是波动光学模块

很多人在这一步犹豫。我的习惯是直接用RF模块的“电磁波，频域”接口做特征频率研究。波动光学模块也能做特征频率，但它的边界条件和端口设置更偏向波导器件，对周期性超表面来说，RF模块的Floquet周期性边界条件用起来更顺手，材料折射率设置也直白。尤其在近红外或可见光波段，用RF模块求解本征模问题时，介质柱的折射率用实数，先不考虑吸收损耗，这样算出来的Q值反映的是纯辐射Q，方便你判断BIC是否真的“不辐射”了。

如果你最终关心的是Q_total，也就是实验上直接测到的品质因子，那要把材料吸收加回来的话，就采用1/Q_total = 1/Q_rad + 1/Q_abs这种思路，分两步算。这个拆分极其重要，后文踩坑部分会详细说为什么不能一开始就开吸收损耗。

2.2 几何参数化：为什么椭圆柱是主流选择

如果只是为了得到某个高Q模式，随便一个矩形柱或者圆柱都行。但你要跟踪的是off-Γ BIC随参数移动的过程，几何参数就必须选得均匀、单调、不引起模式身份突变。

我建议用椭圆柱结构：短半轴设为r1，长半轴设为r2，椭圆度定义为γ = r2 / r1。γ从1变到0.6，结构对称性从圆对称逐渐过渡到椭圆对称。圆对称时，有些模式在Γ点简并；打破圆对称后简并消除，模式沿k空间的色散会重新排列，而off-Γ BIC在这种“由对称操作导致模式重组”的过程中出现或消失。用椭圆度当扫描参数，比直接改柱高或周期更平滑。

拿TiO2介质超表面举例，常见初始量级是：周期P=700 nm，柱高H=300 nm，长半轴r2=200 nm，短半轴r1=150 nm，衬底用熔石英，折射率约1.45，TiO2在近红外折射率约2.4。这些数值只是起跑点，真正的合并条件要靠扫描来确定，但以这个量级开始扫描，计算量比较友好，物理画面也比较清晰。

2.3 Floquet边界和倒空间扫描路径的设置思路

要在COMSOL里扫k向量，常见做法是给周期性超表面单元设置“周期性边界条件”，并启用Floquet子类型，然后指定布洛赫波矢分量kx和ky。很多教程会让你直接在“Floquet周期”里填具体的k_x和k_y数值。但在实际扫描时，应该把kx定义成全局参数：

kx = s * (2*pi / P)

其中s是辅助参数，从0扫到0.5。这样做的好处是，扫描的每一步都对应倒空间里的规则路径，不会因为手动填数出现k点分布不均的问题。扫描路径我一般沿着Γ-X方向，也就是固定ky=0，扫kx；如果要找更精细的合并行为，再沿Γ-M方向补扫。

这方面的设置还有一个坑：不要把“布洛赫波矢”和“入射角”混在一起。特征频率研究里求解的是本征频率ω，同一个本征频率对应自由空间波矢k0 = ω/c；而布洛赫波矢kx是倒空间里的独立坐标。只要你在边界条件里指定kx，求解器就能算出对应的模式频率，不需要再通过k0 = kx/sinθ去反推入射角，否则会给自己制造不必要的非线性。

3. 特征频率研究里如何把Q因子“逼”出来

3.1 从复本征频率到Q值

COMSOL特征频率研究输出的本征频率是复数。物理上，实部代表模式的谐振频率，虚部代表能量衰减速率。Q值按下面这个公式从本征频率提取：

Q = Re(f) / (2 * |Im(f)|)

如果是无损耗介质，Im(f)的来源就是辐射。BIC点对应的理想情况是Im(f) → 0，但数值上永远不可能得到无穷大的Q，因为离散化网格、求解器误差、边界条件近似都会引入虚假泄漏。所以判断BIC，不能等Q变成无穷大，而是看Q在某个k点是不是出现尖锐的极大值，同时这个极大值对网格加密是否稳定。

我在实际操作时习惯用log10(Q)来画图，而不是线性Q。原因很直接：普通共振的Q可能在几千到几万，BIC附近的Q可能到10^5甚至10^7，线性坐标会把低Q部分压成一条直线，根本看不出off-Γ峰的位置。取对数后，Q峰的形状和位置一目了然。

3.2 用辅助参数扫描k空间的两种做法

在COMSOL里让kx变化有两种常用方式。第一种：直接在研究设置的参数列表里把kx设为参数，选择“参数扫描”，把所有想算的k点列出来。这种方式直观、结果清晰，但每个k点都要重新做一次特征频率求解，比较费时间。

第二种：用辅助扫描。在研究设置的“扫描”里加一个辅助参数s，从0以固定步长扫到0.5，同时设定kx = s*(2π/P)。这两种方式本质没区别，我推荐先粗扫。第一步用20个k点、较粗网格把Q曲线趋势摸清楚，确认哪些区域有可疑的Q峰，再在峰附近加密取点。一次特征频率扫描不要贪多，宁可多跑几轮。

3.3 模式追踪：怎么确保你盯的就是同一个模式

特征频率问题最大的麻烦在于，随着kx变化，本征值的排序会跳。同一个物理模式在不同k点可能排到第几个完全不一样，尤其是频率接近的模式密集区，COMSOL列出的第5个特征值可能在下一个k点变成了第3个。

我有两个办法应对。一是把频率范围缩窄，在“特征频率”设置里选择“搜索特征频率附近”，给定一个预期频率，把无关的低频模式和高频伪模式挡在外面。二是每次算完一个k点就输出电场分布，确认模式剖面和你要跟踪的模式一致。我一般用两个判据：谐振频率的连续性，以及模式在z方向的电场剖面或模式偏振特征是否保持一致。如果看到频率出现不连续跳变，或者同一支色散上突然冒出一个Q完全异常的数值，大概率是模式跟踪跳线了，需要回到那个k点检查场分布。

可以做一个小建议：在参数扫描时开启COMSOL的模式跟踪功能，它会尝试从上一个参数点的解来找当前点的对应模式。并不是万能，但在小步长扫描下很好用；一旦步长太大，模式跟踪照样会跟丢。

4. Merging off-γ BIC的完整跟踪流程

4.1 第一步：固定几何参数，摸清Q-kx曲线

拿到模型之后，不要急着扫椭圆度γ，先把一个固定几何下的Q曲线做出来。选择γ=1.0或接近1的值，沿Γ-X方向扫描kx从0到π/P。此时输出模式的Q值曲线，重点关注两个位置：Γ点是否有对称保护的尖峰，以及kx≠0的非Γ区域是否有第二个尖峰。

这一步能帮你确认几件事：Floquet边界设置是否正确、目标模式在k空间里的连续色散长什么样、哪些区域有off-Γ BIC的雏形。如果Γ点没有出现预想中的高Q尖峰，先别急着调物理参数，检查边界条件和网格；Γ点的对称保护BIC在理想条件下应该非常明显。

4.2 第二步：把椭圆度γ加进来，观察两个峰之间的相对运动

确定一个基础几何能出off-Γ峰后，把γ变成扫描参数。一般来说，γ偏离1越远，结构的面内对称性调整越剧烈，off-Γ BIC在倒空间的位置会移动，Q峰值也会变化。扫描时可以先做γ从1.0到0.6的粗扫，步长0.05，每个γ下用20个k点跑趋势。你会看到Q峰的位置随γ连续移动，有些γ下峰朝Γ点方向走，有些γ下两个峰之间靠拢。

这个扫描阶段的关键是同时记录三个量：kx、频率、log10(Q)。不要只记Q峰值，因为Merging BIC的判据是整个k空间里高Q区域的宽度，不是单个孤峰。建议在二维扫描结束时直接画出以γ为纵轴、kx为横轴的Q热力图，能很直观地看到两个高Q轨迹如何相交和汇合。

4.3 第三步：合并判定的客观标准

很多人以为看到两个Q峰扎到一块就算Merging BIC成立。我的经验是，至少再加两个判据。

第一个判据是Q^{-1}的展宽。合并后，在kx区域中间的一段范围内，Q^{-1}应该比合并前低一个量级以上。实际操作中，你可以在某个特定γ参数下，统计Q^{-1}低于某个阈值（比如10^-5）的kx区间宽度，再对比单BIC情况下的宽度。如果合并后宽度明显扩展，说明高Q区域确实变大了。

第二个判据是远场或PB分量检查。Merging BIC在物理上对应两个辐射通道同时关闭，因此模式下某个远场方向的辐射应显著下降。如果有条件，在候选参数下加一层PML并计算模式在结构周围的漏泄场，看看漏泄波是不是真的被抑制了。COMSOL里加PML之后，Im(f)会变成一个有限的较小值，此时Q值才是可信的“开放结构下的数值Q”。

4.4 计算资源不够时的扫描策略

完整二维扫描（γ×kx）非常耗时。我通常分两轮：第一轮用二维模型、粗网格、20个k点×8个γ值，找出靶心区域；第二轮只在靶心区域加密，比如γ扫描步长收窄到0.01，kx点加密到100个。如果机器内存吃紧，可以把频率求解范围进一步收紧，并且把COMSOL的并行线程设置成和核心数一致。计算能跑通是一回事，跑通之后会不会判断结果是另一回事，省下的时间要用在验证上。

5. 那些会让你误判BIC结果的坑，我都替你踩过

5.1 网格疏密直接决定Q值“假峰”的位置

BIC计算对网格极其敏感，这点必须放在最前面说。我遇到过不止一次：粗网格下在某个kx处扫出了一个Q很高的峰，表面上看是off-Γ BIC，加密网格之后峰位置居然移动了，甚至消失。原因很简单：网格离散引入的人工泄漏在k空间分布并不均匀，它可能在局部与真实模式的能量分布“耦合”出虚假的高Q尖峰。

应对方法很笨但有效：对每个候选峰位置，至少做三套网格密度对比，最大单元尺寸分别取λ/10、λ/20、λ/30。如果Q峰值随着网格加密持续上升且峰位几乎不动，说明是物理的BIC；如果峰位漂移或Q值先升后降，那基本是数值假象。这一条是所有BIC计算的第一守则。

5.2 绝对不要一开始就开材料吸收

材料吸收和辐射损耗都会让Im(f)非零，但它们的物理来源完全不同。BIC的判断依据是辐射损耗趋近于零。如果你在模型里用了复折射率，那算出来的Q是总Q，根本无法判断该模式的辐射是否关闭。更麻烦的是，在参数扫描里吸收项会掩盖BIC带来的Q尖峰，让你误以为没有BIC。

推荐流程是：无损耗折射率→找到辐射Q尖峰→加吸收项→对比Q_total。如果你想在无损耗条件下判断合并是否发生，就用无损耗折射率。这也是为什么我前面一直强调用实折射率起步的原因。

5.3 周期边界方向填反导致Γ点对称保护BIC消失

Floquet周期边界条件里的布洛赫波矢分量，顺序千万不能填错。对于二维晶格，第一个分量对应x方向的相移，第二个分量对应y方向的相移。如果你调换了方向，在kx=0时应该受对称保护的BIC会莫名其妙地出现有限辐射，而且Γ点的Q值会低得不正常。验证方法很简单：在kx=0、ky=0处算一个理想对称结构，理论上这个模式Q应该趋近于一个很大的数值。如果算出来只有几千，强烈提示边界条件或网格有问题，先解决这个问题再继续扫描。

另外，周期单元的相对边界之间相位要一致，否则结构的平移会引入额外相位误差，看似细微，实际会让off-Γ BIC位置偏移好几个k点。

5.4 伪模式、局域模式和真实模式的辨别

特征频率一次能吐出一堆本征值，里面混着很多没有物理意义的伪模式，常见的有两种：网格模式（场分布集中在几个网格节点附近，波动剧烈）和边角模式（能量集中在结构边缘或周期边界角落）。这类模式通常Q值极端（要么极低，要么异常高），频率也不按色散规律变化。

我的识别手段是：扫一幅k空间场剖面图，看电场在单元内部是否呈现规整的驻波形式。如果是，继续跟踪；如果场分布里出现强烈的高频波纹或边界尖点，直接放弃，不要在这种模式上浪费时间。判断模式“干净”与否的经验，是所有BIC计算里最容易被人忽略但又最关键的一环。

6. 用PML和远场结果验证Merging BIC是否成立

6.1 加PML对开放结构的必要性

周期性结构在水平方向用Floquet边界，垂直方向是开放的，能量可以从上下表面辐射进自由空间。如果上下边界直接用完美磁导体或完美电导体去截断，辐射会被反射回来，造成Q值虚高。这种虚高在BIC点附近尤其阴险，因为BIC本来辐射就小，反射回来会掩盖真实的有限辐射。

建议在结构上下加PML层，厚度一般取一个波长量级，材料用PML类型，并与结构之间留一段空气或衬底材料的间隔，间隔不小于半个波长。加PML之后，Im(f)会包含真实的辐射通道，Q值才代表开放结构里的实际Q。做Merging BIC判断时，至少在最终候选参数下用PML模型验证一次，不要只在无PML模型里看峰。

6.2 远场方向图和Q热力图怎么用来“下结论”

Merging BIC验证的最后一公里，是看远场辐射有没有在某个辐射方向被压掉。在COMSOL中，如果研究里包含远场域，可以在后处理里加一个远场辐射方向图。合并前后的对比非常明显：合并前，远场图在某几个角度上有明显的漏泄瓣；合并后，这些瓣变窄消失，说明辐射通道彻底关闭。

此外，把所有扫描结果整理成Q热力图时，强烈建议对Q值取对数，并在色标上封顶。比如把超过10^6的Q一律显示为同一颜色，否则数值噪声会造成彩色糊斑，你根本看不清合并区域在哪。这一手处理也是写论文出图时经常用到的，对读者、对合作者都更友好。

6.3 从仿真走向实验对照时要注意的现实差异

如果你后面打算做实验验证，COMSOL模型最好从一开始就包含衬底和上方覆盖介质，而不是悬空薄膜结构。加入衬底后，模式的等效折射率分布会变化，off-Γ BIC的合并参数也会跟着移动。实验上最常用的验证手段是角分辨反射光谱：固定波长，扫描入射角，观察窄带共振在角谱上的展宽行为。Merging BIC在实验上的标志，是在一个较宽的角度范围内都能看到高Q窄带反射谷，而不是只在某个特定角度出现。

就算只做仿真不发文章，我也建议你把衬底加进去。因为衬底的引入会让某些在悬空结构里看起来完美的BIC直接消失，提前发现这个问题，远比你投稿或被审稿人问到再补算要舒服得多。