COMSOL 光学 手性 BIC 仿真 光子晶体板中连续域束缚态 BIC 赋予的手性。 包含正入射斜入射琼斯矩阵透射谱，模式耦合各种透射谱分量，动量空间偏振图
最近在鼓捣手性光子晶体板里的连续域束缚态（BIC），发现这玩意儿跟琼斯矩阵透射谱的耦合效应特别有意思。举个栗子，当结构参数调对了，原本对称性保护的BIC会在动量空间出现手性特征——简单说就是左旋和右旋圆偏振光透射谱突然开始"分道扬镳"。
先上COMSOL建模的干货：光子晶体板选二维六方晶格，晶格常数a=500nm，孔半径r=0.3a，厚度h=220nm。关键操作是在参数化扫描里同时扫入射角θ和方位角φ，注意这里要开波动矢量模块处理斜入射。边界条件设置有个坑——上下端口用端口扫描代替传统散射边界，否则斜入射时端口模式会耦合出错。
% 从COMSOL导出S参数后的数据处理 theta = 0:2:30; % 入射角扫描范围 lambda = 1400:5:1600; % 波长扫描 for n=1:length(theta) S = load(sprintf('S_params_theta%d.mat',n)); JonesMatrix(:,:,n) = [S.s11 S.s12; S.s21 S.s22]; end这段代码处理透射谱的琼斯矩阵分量。重点观察s21的实部和虚部随入射角的变化——当接近BIC点时，会看到s21real² + s21imag²趋近于零，同时s12分量突然增强，这是模式耦合触发的典型特征。
COMSOL 光学 手性 BIC 仿真 光子晶体板中连续域束缚态 BIC 赋予的手性。 包含正入射斜入射琼斯矩阵透射谱，模式耦合各种透射谱分量，动量空间偏振图
动量空间的偏振演化更带感。用COMSOL的频域电场导出功能，配合后处理脚本做斯托克斯参数计算。有个骚操作是把x-y偏振基变换到圆偏振基：
def Jones2Circular(E_x, E_y): E_RCP = (E_x - 1j*E_y)/np.sqrt(2) E_LCP = (E_x + 1j*E_y)/np.sqrt(2) return np.abs(E_RCP)**2, np.abs(E_LCP)**2在Γ点附近扫动量kx,ky时会发现，当结构引入面内不对称时，原本简并的RCP和LCP透射谷开始分裂。这其实对应着BIC从对称保护态向参数可调的手性态转变。有个反直觉的现象：斜入射30度时，TM偏振的透射谱在1550nm附近会出现"双峰劈裂"，这是TE-TM耦合与BIC拓扑荷相互作用的结果。
最后给个实用技巧：处理动量空间偏振图时，用二维傅里叶变换把实空间电场转成动量分布后，记得做高斯滤波消除数值噪声。COMSOL的网格收敛性测试表明，当最大单元尺寸小于λ/8时，偏振椭圆度的计算误差能控制在5%以内。
