COMSOL 多极分解，分方向多级展开通用模型，电磁场，面上箭头，透射率光学 BIC 仿真

COMSOL 多极分解，分方向多级展开通用模型，电磁场，面上箭头，透射率光学 BIC 仿真。 COMSOL 准 BIC控制石墨烯临界耦合光吸收。 COMSOL 光学仿真，石墨烯，光吸收，费米能级可调

咱们今天来点硬核的——玩转COMSOL里的边界态控制。最近在实验室折腾光学BIC和石墨烯调控，发现多极分解这玩意儿真是神器。先给你们看看我的模型界面截图（此处假装有图），面上那些五彩斑斓的箭头可不是装饰，它们暴露了电磁场分布的玄机。

做BIC仿真最头疼的就是模式分解，直接在模型里敲这段代码：

mphglobalexpr('Q_scat', 'emw.Qsca_emw'); mphglobalexpr('Q_abs', 'emw.Qabs_emw');

这俩全局变量算的是散射和吸收效率。关键是要在频域研究中勾选"存储场解"，否则后续多极分解会报错。记得把计算域设置成比结构大1.5倍，不然边界反射会污染结果。

重点来了——分方向多级展开。在电磁波模块的派生值里新建"多极分解"，设置方位角展开阶数从-3到3足够用。有个骚操作是同时勾选笛卡尔和球坐标系输出，对比着看更容易发现模式特性。比如当某个模式在笛卡尔坐标系下呈现偶极-四极混合态，但在球坐标系里呈现纯环形分布，基本可以确定是准BIC。

石墨烯调参才是真功夫。这里有个自定义材料参数的模板：

sigma_drude = (e^2 * Ef)/(pi*hbar^2 * (1i*omega - gamma)) * 0.5;

注意这里的费米能级Ef要关联到全局参数，配合静电模块做双向耦合。实测当石墨烯电导率的实部等于介质腔的辐射损耗时，吸收率会突然飙升到90%+，这就是临界耦合点的魅力。

后处理时必看透射率曲线突变点。用参数化扫描连续调节费米能级，会观察到明显的反交叉现象。有意思的是当结构对称性被刻意破坏时，原本的BIC会转变成具有高品质因子的准BIC，这时候吸收峰宽度突然收窄，配合石墨烯的主动调控可以实现光开关效果。

最后说个坑：网格剖分一定要用曲边形单元，特别是处理微纳结构中的弧形边界时，直角剖分会导致多极分解结果出现鬼影。还有别忘了在求解器设置里把残差容差调到1e-6以下，否则模式纯度指标会飘得亲妈都不认识。