COMSOL模拟下的六角晶格光子晶体四重简并狄拉克点与零折射率复现研究

comsol能带复现六角晶格光子晶体四重简并狄拉克点零折射率打开COMSOL新建模型时的手总是微微发抖——这破软件动不动就吃内存到32GB封顶。六角晶格光子晶体建模的关键在于准确构建蜂窝状排列的介质柱先打开几何面板用参数化建模更靠谱% 六角晶格基础参数 a 1e-6; % 晶格常数 r 0.3*a; % 介质柱半径 theta 0:60:300; % 六边形角度分布这里用极坐标生成六个方位角注意COMSOL的几何内核对极坐标转换存在约3%的坐标漂移误差得手动修正循环偏移量。建完单胞后别急着复制阵列先检查原胞边界的周期性——在边界条件里勾选Floquet周期边界时发现TM波的电场连续性需要特殊处理得手动设置电场切向分量连续。comsol能带复现六角晶格光子晶体四重简并狄拉克点零折射率频率扫描设置藏着魔鬼细节study model.study.create(study1); study.feature.create(param, Parametric); study.feature(param).set(plist, linspace(0.4,0.6,50));扫频范围得精确覆盖狄拉克点出现的归一化频率0.5附近。求解后发现特征频率在0.497a/λ处出现四个模式交汇这四重简并的能带结构像突然塌陷的沙堆——典型的狄拉克锥特征。查看电场分布时按住Ctrl鼠标滚轮放大十倍能清晰看见四个模式分别对应着不同方向的相位涡旋。零折射率现象验证需要魔改后处理模块% 计算等效折射率 epsilon_eff (E_total*epsilon*E_total)/norm(E_total)^2; mu_eff (H_total*mu*H_total)/norm(H_total)^2; n_eff sqrt(epsilon_eff*mu_eff);当频率对准狄拉克点时算出来的n_eff实部突然跳水到1e-3量级。这时候打开参数化扫描结果动画能看见电磁波像在迷宫里瞬间完成空间跳跃——相前传播延迟几乎消失这验证了零折射率材料的时空穿透特性。不过要注意COMSOL的频域求解器在接近奇异点时容易数值发散记得把残差容差调到1e-6以下才能稳住计算。