COMSOL模拟环偶极子增强磁光克尔效应

comsol环偶极子增强磁光克尔效应最近在玩COMSOL模拟磁光克尔效应的时候发现环偶极子结构对增强效果特别有意思。这玩意儿就像给光波装了个磁力放大器咱们今天直接上干货看看怎么用COMSOL玩转这个现象。先搞明白环偶极子怎么在模型里构建。用App开发器写个参数化脚本最方便double R 800e-9; // 环半径 double w 100e-9; // 环宽度 model.geom(geom1).feature().create(r1, Ring); model.geom(geom1).feature(r1).set(radinner, R-w/2); model.geom(geom1).feature(r1).set(radouter, Rw/2);这个环形结构的关键在于宽度和半径的比例太宽了会破坏环形电流的对称性太窄了场增强不够。建议先用参数扫描找最佳值后面再固定参数做深入分析。材料设置要特别注意磁光特性。用YIG钇铁石榴石这类磁光材料时介电张量得这么写epsilon [eps_xx 0 0; 0 eps_yy i*eps_xy; 0 -i*eps_xy eps_yy];这里的非对角项eps_xy就是磁光效应的来源记得在物理场设置里勾选启用磁光效应。有个坑要注意COMSOL默认用的是高斯单位制实际输入得换算成国际单位制。边界条件配置是另一个重点。环形结构外围用完美匹配层PML顶部设置端口激发model.physics(emw).feature().create(port1, Port, 3); model.physics(emw).feature(port1).set(PortName, 入射端口); model.physics(emw).feature(port1).set(PortType, Circuit);这里有个骚操作把入射波的偏振方向设置成与环形结构对称轴成45度角可以同时激发TE和TM模式更容易观察到克尔旋转角的增强效果。comsol环偶极子增强磁光克尔效应求解器配置建议用频域扫频但要注意网格划分。环形边缘处要做至少三层边界层网格model.mesh(mesh1).feature().create(bndl1, BoundaryLayers); model.mesh(mesh1).feature(bndl1).set(nbrlayers, 3); model.mesh(mesh1).feature(bndl1).set(thickness, 20e-9);后处理提取克尔角时用这个公式最靠谱θK (1/2i)*ln(rp/r_s)用COMSOL的派生值直接计算rp emw.S21_1; // p偏振反射系数 rs emw.S21_2; // s偏振反射系数 theta_K 0.5*imag(log(rp/rs));跑完仿真你会发现在特定共振频率下克尔角能飙到普通结构的5倍以上。这是因为环形结构产生了强烈的环形偶极矩和入射光的轨道角动量发生耦合把原本微弱的磁光效应给放大了。调试时容易遇到收敛问题建议先关掉非线性材料特性等找到共振峰后再逐步开启高阶效应。记住保存每次的计算快照COMSOL的恢复点功能能救命。最后提醒环形结构的方位角对称性直接影响结果精度做参数优化时一定要锁定旋转对称边界条件。