磁偶极子贡献准BIC，多极子分析，斜入射反射相位计算。
在超表面设计里常遇到个头疼问题——结构参数调麻了也找不到BIC点。这时候试试磁偶极子可能就有惊喜，最近发现用环形电流结构激发磁偶极矩，配合特定晶格排列能稳定触发准BIC现象。看这段COMSOL后处理脚本：
extract_magnetic_dipole = lambda fields: np.sum([(r × J) for J in current_density]) phase_curve = [calc_phase(extract_magnetic_dipole(sim_data)) for wavelength in sweep_range] plt.plot(wavelengths, phase_curve, marker='*', color='#FF6B6B')关键在环形电流的旋度计算（那个叉乘操作），这玩意儿直接对应磁偶极矩强度。当相位曲线出现π突跳时，赶紧用多极子分解工具验证下：
[multipoles, params] = jcmpt_multipole('project.jcm', 'Frequency', linspace(200e12,300e12,50)); hold on; plot(params.frequency, abs(multipoles.magnetic_dipole).^2, 'LineWidth',2); plot(params.frequency, abs(multipoles.electric_quadrupole).^2, '--'); legend('Magnetic Dipole','Electric Quadrupole');这时候会看到磁偶极子分量突然暴增，而电四极子像被掐住脖子似的往下掉——典型的BIC特征。但实际做斜入射时相位响应容易跑偏，用这个算法锁定相位奇点：
def calc_oblique_phase(theta): rc, _ = fdfd_solver.solve(theta=theta) return np.angle(rc[:,1,1]) # 取TE极化分量 thetas = np.linspace(0, 40, 21) phase_maps = np.array([calc_oblique_phase(np.deg2rad(t)) for t in thetas]) plt.contourf(wavelengths, thetas, phase_maps.T, levels=20, cmap='twilight')注意看色卡图中出现的涡旋状相位分布，那个涡旋中心对应的就是鲁棒性最强的准BIC工作点。调试时遇到过坑：当结构对称性破缺超过临界值，这个涡旋会分裂成两个，这时候得回炉重新优化环形结构的开口角度。
)
