基于COMSOL弱形式方程求解色散光子晶体能带

摘要

为了研究光子晶体复杂色散关系,基于有限元数值仿真软件COMSOL Multiphysics, 我们发展了一套基于COMSOL弱形式方程光子晶体能带的计算方法。通过频率求解波矢k的本征值,将传统方法很难求解的色散能带问题转化为简单线性的本征值问题。利用这一新的数值工具,我们首先对比分析了无色散简单正方晶格光子晶体的体能带结构,弱形式方法与现有方法得到完全一致的能带结果,验证了这种计算方法的正确性。利用这种求解方法能快速求解光子晶体的等频率曲线(equi-frequency contour), 我们也验证了我们研究的正方晶格光子晶体具有类狄拉克锥能带结构。我们将这种方法应用于具有色散负介电常数背景下类石墨烯蜂巢晶格光子晶体的能带结构计算,这是传统光子晶体能带计算方法难以求解的结构。我们发现,在这种光子晶体的胡须状(Bearded)和锯齿状(Zigzag)界面都存在光子界面态。