极化无关连续束缚态BICs(多极子分析) COMSOL光子晶体超表面模拟光在周期性结构中跳舞的时候偶尔会出现任性的行为——明明应该存在的泄漏模式突然自我封闭这就是传说中的连续束缚态BICs。最近实验室里的小张拿着COMSOL模型折腾两周后突然发现某些特殊排列的纳米柱居然能无视光的偏振方向稳定锁住能量这事儿有点意思。先看个实际的模型边长500nm的钛酸钡方形纳米柱按六方晶格排列在二氧化硅基底上。我们在COMSOL里建模时边界条件设置得像三明治——上下用完美匹配层左右周期性边界。重点来了当晶格常数调整到720nm时参数扫描结果显示无论线偏振光怎么旋转1550nm附近总有个尖锐的共振峰死守位置不动。% 参数扫描脚本片段 for a 700:10:750 model.param.set(period, num2str(a)); model.study(std1).run; extract_resonance_peak(); % 自定义峰位提取函数 end plot(peak_positions, LineWidth,2);这段代码循环改变晶格常数并记录共振峰位跑完数据就能明显看到某个特定参数下峰位突然钉在固定波长。有意思的是这时候在频域监视器里看到的场分布既不是单纯的偶极子也不是四极子更像是多种模式在玩叠叠乐。多极子分解这时候就该上场了。我们通过后处理计算各阶电磁多极矩的贡献占比% 多极矩计算核心部分 Qe sum(eps*r.*E_field); % 电四极矩近似计算 Qm sum(mu*r.*H_field); % 磁四极矩分量 scattering abs(Qe)^2 abs(Qm)^2; % 散射功率估算当电偶极、磁偶极和环形偶极的散射功率曲线出现三重交叉时系统进入非辐射状态。这时候在动量空间画能带图会看到原本应该存在的辐射通道神奇地关闭了就像高速公路突然多出个隐形收费站。极化无关连续束缚态BICs(多极子分析) COMSOL光子晶体超表面模拟不过实际操作时有个坑——网格剖分不够细的话环形偶极项会被数值噪声淹没。有次我偷懒用默认网格结果多极子分解后的数据像心电图一样乱跳。后来把纳米柱边缘的网格尺寸压到1nm以下数据才老实下来。这告诉我们搞光子晶体模拟咖啡可以省计算资源不能省。最后验证偏振无关性时可以暴力旋转入射光偏振角for theta 0:15:180 model.param.set(polarization_angle, theta); model.study(std2).run; track_quality_factor(); % 自定义Q因子追踪 end当Q因子曲线平得像青藏高原时就知道这个BIC确实对偏振免疫了。不过这种理想情况只在严格对称时成立实际制备时纳米柱的尺寸误差超5%就会破功所以文章里的漂亮数据都是烧了无数机时后的幸存者。这种偏振无关BIC在传感领域挺有潜力。去年有组在类似结构上做到折射率传感灵敏度超过30000 nm/RIU不过他们没公开COMSOL模型文件——懂的都懂这行当的模型文件比祖传秘方还金贵。反正下次组会我要建议老板买个128核服务器现在的笔记本跑三维参数优化已经快冒烟了。