极化无关BIC光子晶体超表面设计与COMSOL仿真实践

1. 极化无关BIC现象与光子晶体超表面基础

连续束缚态（Bound states in the continuum, BICs）是光子晶体和超表面研究中的特殊现象，表现为在辐射连续谱中存在的局域态。与传统束缚态不同，BICs理论上具有无限大的品质因子（Q因子），其能量被完美束缚在结构中不向外辐射。极化无关特性则指该现象对入射光偏振方向不敏感，这在实用化器件设计中具有重要价值。

钛酸钡（BaTiO₃）纳米柱选择源于其显著的电光效应和非线性光学特性，相对介电常数在可见至近红外波段可达30-40。六方晶格排列提供了C₆v对称性，这是实现极化无关BIC的关键。二氧化硅基底（SiO₂）不仅提供机械支撑，其低折射率（~1.44@1550nm）还能增强场局域效果。

在COMSOL Multiphysics中建模时，电磁波频域（emw）模块通过求解麦克斯韦方程组处理光与纳米结构的相互作用。完美匹配层（PML）的厚度通常设为工作波长的1/2，吸收率需达到99.9%以上以避免虚假反射。周期性边界条件通过Floquet端口实现，其相位匹配关系直接影响能带结构计算精度。

2. 多极子分解理论与数值实现

多极子分解是理解BIC物理本质的核心工具。对于亚波长纳米结构，电偶极（ED）、磁偶极（MD）、电四极（EQ）、磁四极（MQ）和环形偶极（TD）是主要贡献者。散射消逝条件可表示为：

$$
\sum_{i} |a_i|^2 = 0 \quad (i = ED, MD, TD,...)
$$

其中$a_i$为各多极矩的散射系数。在COMSOL中实现需通过后处理积分计算：

matlab复制% 电偶极矩计算
P = sum(eps*(r - r0).*E_field)*dV; 
% 磁偶极矩计算
M = 0.5*sum(cross((r - r0),conj(H_field)))*dV;
% 环形偶极矩需计算位移电流环
Jd = 1i*omega*eps*E_field;
T = 0.1*sum(cross((r - r0),Jd))*dV;

网格敏感性测试表明，当纳米柱边缘网格尺寸大于5nm时，环形偶极矩计算误差可达40%。建议采用：

1. 曲率自适应网格（最大尺寸λ/20）
2. 边界层网格（3层，厚度比1:2:4）
3. 纳米柱顶部增强剖分（最小尺寸1