刘芷宁1. 超表面与BIC：传统设计思路的局限性
在光学超表面领域，束缚态连续体（Bound states in the continuum, BIC）因其独特的品质因子和光场局域能力，近年来成为研究热点。传统超表面实现BIC的主流方法是通过打破结构对称性——比如在周期性结构中引入不对称的几何参数或材料分布。这种对称破缺会在原本连续的辐射谱中产生离散的束缚态。
但这种方法存在几个固有缺陷：首先，对称性破缺程度直接影响BIC的品质因子，设计时需要精确控制结构参数，工艺容差极小；其次，大多数BIC超表面对入射光的偏振态敏感，往往只能在特定偏振下工作；再者，通过结构不对称性实现的BIC通常伴随着复杂的近场耦合效应，给理论分析和优化带来挑战。
关键提示：传统BIC超表面的偏振依赖性主要源于非对称结构对电磁场不同偏振分量的差异化响应。例如，x方向的结构不对称会导致TE和TM偏振光产生不同的散射相位。
需要模型API调用？ 免费领10W Token，多模型网关一键接入 Claude、DeepSeek 等主流模型。
2. COMSOL在超表面设计中的独特优势
COMSOL Multiphysics作为一款多物理场仿真软件，在超表面设计中展现出独特价值。其核心优势在于：
2.1 全波电磁仿真能力
通过RF模块或波动光学模块，COMSOL可以精确求解麦克斯韦方程组，特别是对亚波长结构的近场耦合效应有出色表现。对于超表面这类包含复杂边界条件的系统，软件内置的有限元方法（FEM）能自动适应各种几何形状。
2.2 参数化扫描与优化
COMSOL LiveLink™ for MATLAB允许用户建立参数化模型，批量扫描结构参数（如单元尺寸、材料属性等）。结合内置的优化模块，可以实现自动化的超表面性能优化。例如，可以通过Nelder-Mead算法寻找使谐振峰最锐化的结构参数组合。
2.3 多物理场耦合分析
超表面往往涉及电磁-热-力等多场耦合效应。COMSOL能同步求解不同物理场的相互作用，这对分析BIC状态下的能量耗散机制尤为重要。软件中的"多物理场"节点可以直观地定义各场之间的耦合关系。
3. 偏振不依赖BIC的实现路径
要实现偏振无关的BIC超表面，COMSOL提供了几种创新设计思路：
3.1 旋转对称性结构设计
在模型中构建具有C4或C6旋转对称性的超表面单元。这类结构对任意线偏振光的响应相同，从根本上避免了偏振依赖性。在COMSOL中，可以通过"几何序列"功能快速生成旋转对称结构。
具体操作步骤：
- 在几何界面绘制基础结构（如十字形、圆环等）
- 右键点击"几何"→"变换"→"旋转"
- 设置旋转角度为90°（C4对称）或60°（C6对称）
- 勾选"复制"选项生成完整对称结构
3.2 双各向异性超表面
通过精心设计单元结构，使其在x和y方向具有等效但独立的谐振特性。在COMSOL中，可以通过参数化扫描验证不同偏振下的谐振频率一致性。典型的实现方法包括：
- 嵌套十字结构
- 正交排列的劈裂环谐振器
- 多层堆叠的各向异性结构
3.3 拓扑优化设计
利用COMSOL的密度法拓扑优化功能，以偏振无关性为目标函数自动生成最优结构。具体流程：
- 定义设计区域和材料分布变量
- 设置TE和TM偏振的仿真场景
- 建立目标函数：min(|f_TE - f_TM|)
- 运行优化并导出最佳结构
4. 仿真中的关键技术细节
4.1 周期性边界条件设置
超表面仿真必须正确设置周期性边界条件：
- 在"定义"中创建"周期"节点
- 选择相对的边界对
- 设置Bloch周期条件：Floquet周期类型
- 输入波矢分量（根据入射角度计算）
常见错误：忘记在端口设置中勾选"周期端口"，导致仿真结果不准确。
4.2 网格划分策略
BIC对网格精度极为敏感，推荐采用：
- 边界层网格：在金属-介质界面添加3-5层边界层
- 扫掠网格：对于规则结构域优先使用
- 自定义尺寸：在谐振区域设置更细的网格
当遇到"创建域的扫掠网格失败"错误时，可以：
- 检查几何是否存在微小间隙
- 尝试简化复杂几何特征
- 使用自由四面体网格作为备选方案
4.3 结果后处理技巧
- 使用"全局计算"提取谐振频率：es.Q因子计算
- 电场分布可视化：启用"箭头图"显示近场增强
- 参数提取：通过"派生值"计算散射参数
5. 实测案例：偏振无关BIC超表面设计
以一个实际设计案例说明完整流程：
5.1 模型构建
- 创建20×20μm²的单元结构
- 设计双层十字结构，上层硅（ε=12），下层金
- 设置C4旋转对称性
- 定义周期性边界条件
5.2 材料属性定义
matlab复制% COMSOL with MATLAB脚本示例
model.param.set('epsilon_Si', '12-0.01i');
model.param.set('gold_eps_inf', '1');
model.param.set('gold_omega_p', '2*pi*2.175e15');
model.param.set('gold_gamma', '2*pi*6.5e12');
5.3 仿真结果分析
通过参数扫描发现：
- 传统非对称结构：TE/TM谐振频率差达12%
- 旋转对称设计：TE/TM频率偏差<0.5%
- 品质因子Q均超过2000
5.4 实验验证要点
- 加工公差控制在±20nm以内
- 光学表征时需使用宽谱光源
- 注意消除基底背向反射的影响
6. 进阶应用与挑战
6.1 动态可调BIC超表面
通过引入相变材料或液晶层，实现BIC状态的电/光调控。在COMSOL中需要：
- 定义材料参数随控制变量的变化关系
- 设置多步骤研究：先稳态后频域
- 使用"事件"接口处理相变过程
6.2 制造工艺适配
仿真到实际的gap主要来自：
- 材料损耗（仿真中常被低估）
- 侧壁粗糙度效应
- 层间对准误差
建议在COMSOL中：
- 添加表面粗糙度模型
- 进行蒙特卡洛公差分析
- 建立等效电路模型辅助验证
6.3 多物理场耦合效应
BIC状态下的热-力效应不可忽视：
- 光热转换导致局部温升
- 热膨胀引起结构形变
- 形变反作用于光学性能
解决方法：
- 启用"热膨胀"多物理场耦合
- 设置双向耦合迭代
- 监测关键参数的温度系数
在实际项目中，我们发现旋转对称结构虽然解决了偏振依赖问题，但带来了新的挑战——模式简并度增加可能导致杂散谐振出现。这时需要在COMSOL中通过本征频率研究明确识别出真正的BIC模式，通常需要结合场分布和能量衰减分析进行判断。