Comsol连续体中的束缚态BIC。 涉及能带计算与Q因子计算，包含一维光栅和二维光子晶体板。 注: 不包含拓扑荷计算。

在光子晶体和超表面设计中，连续体中的束缚态（BIC）就像是被光场意外捕获的"幽灵模式"。这种特殊状态下的光既不被材料吸收也不向外辐射，理论上Q因子能达到无限大。今天咱们用Comsol扒一扒一维光栅和二维光子晶体板里的BIC实现套路。

一维光栅的对称陷阱

先搞个550nm周期的硅光栅，厚度取220nm。核心在于破坏结构对称性来触发BIC泄露。建模时注意设置周期性边界条件：

% COMSOL LiveLink脚本片段
model.param.set('period', 550e-9);
model.geom('geom1').feature.create('grating', 'Rectangle');
model.geom('geom1').feature('grating').set('size', {'period' '220e-9'});
model.physics('wave').feature.create('pw', 'PlaneWave');
model.physics('wave').feature('pw').set('WaveType', 'TE');

这段代码建立了TE偏振的平面波入射环境。当我们将光栅倾斜角从0°扫到5°时，神奇的事情发生了——在特定角度下Gamma点的Q值突然飙升三个数量级。这种对称性保护的BIC可以通过监测散射损耗的突变来捕捉：

% 扫参脚本
sweep = model.study.create('sweep');
sweep.set('plist', 'linspace(0,5,50)');
sweep.set('pname', 'theta');
model.sol('sol1').feature.create('s1', 'StoreSolution');

记得在频域求解后提取S参数，用洛伦兹拟合估算Q因子。实际操作中发现，当结构误差超过10nm时Q值会断崖式下跌，这解释了实验中BIC难捕捉的原因。

二维光子晶体的暗波魔术

六方晶格光子晶体板更带劲。构建三角晶格时，孔的旋转对称性至关重要：

model.geom.create('hex', 'Hexagon');
model.geom('hex').set('side', 'a'); % a为晶格常数
model.geom.create('holes', 'Cylinder');
model.geom('holes').set('pos', {'0' '0' '-h/2'}); % h为板厚

设置Floquet边界时要注意波矢扫描方向。有个骚操作：在Gamma点附近设置泄漏通道监测器，当模式的有效折射率接近环境介质时，辐射损耗突然归零。这对应着参数空间中的奇异点。

Q因子计算的三个野路子

1. 时域衰减法：跑个脉冲仿真，记录场衰减包络。用半衰期算Q时记得排除数值耗散影响
2. 本征频移：在频域求解器中，虚部频率直接关联Q值（Q = Re(ω)/(2*Im(ω))）
3. 线宽反推：散射谱的3dB带宽计算，适合弱耦合情况

实测发现二维结构的Q值对网格密度极其敏感。当采用自动网格时，结果可能飘忽不定。建议在孔边缘做边界层网格：

model.mesh('mesh1').feature.create('boundLyr', 'BoundaryLayer');
model.mesh('mesh1').feature('boundLyr').set('thickness', 'lambda/20');

最后丢个冷知识：BIC的Q值理论上随系统尺寸呈指数增长，但实际受限于材料损耗。某次仿真中把硅换成GaAs后，Q值直接从1e4跌到5e3，可见材料选择的玄学。

玩转BIC的关键在于理解模式泄露的开关机制。下次设计时不妨先找个对称破缺参数当调节旋钮，配合Comsol的参数化扫参，说不定就能钓到那条Q值暴涨的"大鱼模式"。记住，仿真里的完美BIC就像理想气体——知道它存在就好，实际应用还得留足容差余量。