1. 项目概述：六角晶格高阶BIC的合并行为仿真

在光子晶体和超材料研究中，边界连续条件（Boundary Continuity Conditions, BIC）因其特殊的非辐射特性成为热点课题。本项目使用COMSOL Multiphysics对六角晶格结构中的高阶BIC合并行为进行全流程仿真，重点复现三个关键指标：能带结构、品质因子（Q因子）和远场偏振特性。这种仿真对设计高Q值光学谐振器、传感器和激光腔体具有直接指导意义。

六角晶格因其C6v对称性常表现出丰富的BIC模式，而高阶BIC的合并行为更涉及复杂的模式耦合机制。通过COMSOL的RF模块和波动光学模块组合仿真，我们可以完整捕捉从模式形成到远场辐射的全过程物理现象。下面将分步骤详解实现方法。

2. 核心模型构建与参数设置

2.1 六角晶格单元设计

采用硅基介质柱（折射率n=3.5）在空气背景中排列形成六角晶格，晶格常数a=500nm，圆柱半径r=0.3a。在COMSOL中通过"几何>参数"定义这些基础参数，便于后续参数化扫描。

关键技巧：

• 使用"阵列"功能生成晶格时，选择"旋转对称"而非简单线性阵列，确保严格的六角对称性
• 对圆柱边缘进行微米级倒角（0.01a），避免数值计算中的场奇异性

2.2 物理场选择与边界条件

添加"电磁波，频域"接口，设置：

• 研究频率范围：150-250THz（对应近红外波段）
• 边界条件：
  • 晶格单元侧边：Floquet周期性边界
  • 上下边界：完美匹配层（PML）厚度设为λ/2
  • 介质柱表面：连续边界条件

特别注意：

• Floquet边界需要正确定义k矢量方向，六角晶格需设置两个基本矢量方向（60°夹角）
• PML的坐标拉伸系数建议设为2-3，过大会导致数值不稳定

3. 能带结构计算方法

3.1 布里渊区路径选取

对六角晶格采用Γ-M-K-Γ路径扫描：

1. 在"研究1"中添加"频域"研究
2. 创建参数化扫描，k矢量沿：
   • Γ→M：kx从0到2π/(a√3), ky=0
   • M→K：kx固定2π/(a√3), ky从0到2π/(3a)
   • K→Γ：沿对角线返回原点

3.2 模式分析技巧

• 使用"模式分析"研究类型，设置搜索频率接近预期BIC位置
• 网格设置：在介质柱边缘进行至少5层边界层网格划分
• 求解器配置：启用"矢量元"计算，阶数设为二次

典型问题排查：

• 若出现模式交叉，需在"频域"研究中勾选"搜索附近频率"
• 能带出现断裂时，检查周期性边界是否正确定义了相位延迟

4. 品质因子精确计算

4.1 时域衰减法

1. 在已求解的频域模型上添加"电磁波，瞬态"研究
2. 设置高斯脉冲激励（中心频率对应BIC频率）
3. 添加"探针"记录场能量随时间衰减
4. 对衰减曲线做傅里叶变换得到Q=ω₀/Δω

4.2 本征频率法

通过频域本征值研究直接获取：

Q = real(ω)/(2*imag(ω))

关键参数：

• 本征值搜索方法：Arnoldi迭代
• 偏移频率设为预期BIC频率的95%
• 最大模式数设为10（确保覆盖所有简并模式）

实测对比：

• 六角晶格Γ点Q因子通常可达10⁶-10⁸量级
• 高阶BIC合并时可能出现Q因子突变（提升或降低2-3个数量级）

5. 远场偏振特性提取

5.1 远场计算设置

1. 在"电磁波"接口添加"远场域"节点
2. 定义半球形远场表面（半径≥10λ）
3. 设置偏振分析：
   • 添加"远场偏振椭圆"计算
   • 定义观测角度θ=0-90°, φ=0-360°

5.2 数据后处理技巧

• 使用"辐射方向图"绘图类型显示偏振态
• 导出Stokes参数进行偏振度计算：

  DOP = √(S1²+S2²+S3²)/S0

• 对合并BIC区域特别注意：
  • 偏振奇点位置变化
  • 拓扑荷数守恒验证

6. 高阶BIC合并行为分析

6.1 参数化扫描策略

通过改变晶格对称性（如引入椭圆度）触发BIC合并：

1. 定义变形参数δ=1-(b/a)，b为y方向半径
2. 从0到0.2进行参数扫描，步长0.02
3. 监测以下指标变化：
   • 本征频率实部/虚部
   • 远场偏振奇点位置
   • 模式场分布重叠积分

6.2 合并判据验证

当两个BIC模式满足：

• 频率差Δf < 本征线宽
• 模式重叠积分>0.7
• 偏振奇点发生湮灭 即判定发生合并行为

7. 常见问题解决方案

7.1 收敛性问题

现象：Q因子随网格加密不收敛 解决：

• 检查PML吸收效率（反射率应<-60dB）
• 改用矢量元+曲边单元组合
• 启用"几何形状修正"选项

7.2 模式识别错误

现象：自动模式编号混乱 解决：

• 手动指定初始场分布
• 通过"场重叠"功能追踪模式演化
• 导出模式场数据用MATLAB做相关性分析

7.3 远场计算内存不足

优化方案：

• 改用远场近似公式替代全波计算
• 降低远场网格密度（θ/φ步长增至5°）
• 启用"分布式计算"选项

8. 模型验证与实验对比

8.1 收敛性测试标准

• 能量误差：<1e-4
• 频率偏移：<0.1%
• Q因子变化：相邻加密<5%

8.2 典型实验结果

• 能带对比：Γ点频率误差<0.3THz
• Q因子差异：仿真值通常为实测值2-3倍（考虑材料损耗）
• 偏振角偏差：<5°

实际操作中发现，当圆柱高度达到0.6a时，会出现新的面外偏振BIC模式，这为设计多波段BIC器件提供了新思路。建议在参数扫描时同步监控z方向场分量，可能发现意外的模式耦合现象。