2026/4/7 15:00:09
网站建设
项目流程
找网站做q币,织梦网站如何播放mp4,一般的企业网站开发价格,it培训网站磁偶极子贡献准BIC#xff0c;多极子分析#xff0c;斜入射反射相位计算。在超表面设计里常遇到个头疼问题——结构参数调麻了也找不到BIC点。这时候试试磁偶极子可能就有惊喜#xff0c;最近发现用环形电流结构激发磁偶极矩#xff0c;配合特定晶格排列能稳定触发准BIC现象…磁偶极子贡献准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(thetatheta) 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, levels20, cmaptwilight)注意看色卡图中出现的涡旋状相位分布那个涡旋中心对应的就是鲁棒性最强的准BIC工作点。调试时遇到过坑当结构对称性破缺超过临界值这个涡旋会分裂成两个这时候得回炉重新优化环形结构的开口角度。