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以下才能稳住计算。