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", R+w/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的恢复点功能能救命。最后提醒：环形结构的方位角对称性直接影响结果精度，做参数优化时一定要锁定旋转对称边界条件。