简单六方结构二维光子晶体能带的COMSOL模拟

1、 简单六方结构二维光子晶体能带的COMSOL模拟 北京东之星应用物理研究所 伍勇1.引言 COMSOL携带的案例库里，其中一篇<Bandgap Analysis of a Photonic Crystal>(以下简称< Bandgap >)对砷化镓简单正方格子2D光子能带进行了完整计算和研究。本文将程序用于简单六方结构，并将结果在此做一介绍。2. 关于 Floquet (弗洛盖） 波矢 这是入门COMSOL光子晶体能带模拟的重要概念，在另一案例<Porous Absorber>中，在Floquet周期性边界条件一段写明： 由此我判断Floquet 波矢就是B

2、loch（布洛赫）波矢,但“帮助”文档中有： ，以正格子基矢表示（其文没有任何几何插图和物理说明），使我决定必须在六方格子中选择矩形单胞作为周期单元，以使计算机程序能够运行我的几何方案。3.几何建模 图1作为试探选择的几何模型，圆形柱代表以GaAs作为格点材料，在空气介质中周期性排列，形成二维六方结构人造晶体。 是晶格常数。 是z方向的单位矢量以上根据倒格子基矢定义计算出，及其分量。由倒格子基矢 ， ，构建长方格子的布里渊区也是长方结构如图2： 4.二维光子晶体主方程 COMSOL在< Bandgap > “模型开发器” 电磁波，频域 写出方程形式如下： ， 在< Bandg

3、ap >中，下面目录 波方程，电 中直接简化为， 电磁波在光子晶体中的传播遵从麦克斯韦方程，上述方程可由麦克斯韦方程组出发导出 介质中的麦克斯韦方程组 ， 在电介质中一般认为自由电荷，自由电流密度（电导率）为零。本文档不考虑磁性质， 传播模态电场函数COMSOL表达为：， 在周期结构中，它应具有Bloch波的性质，不考虑衰减损耗。注意这里次上标的符号与我们习惯的教科书里正负符号相反。由(2)两端取旋度，将（4）代入得： ， ， ， ， ， 为真空波数，应用（5），如期得到前面提到的< Bandgap >中的假设波方程： 利用 偏振时（又称横磁场模式）对应的二维光子晶体主方程是

4、： 这里电场方向平行于z轴，而磁场方向只分布在x-y平面内。 二维光子晶体主方程 是计算光子能带的出发点和源依据。5.折射率色散函数 是COMSOL计算能带的第二个关键点，描述折射率随频率参量的变化。格点材料砷化镓的，在< Bandgap >有如下定义n_GaAs=3.3285e5s/ m*c_const/f+3.5031，参变量f,单位Hz；下限：0，上限：1表示真空光速，我在软件Mathematica 上化简这个函数，给出 ，并绘制图3如下： 经COMSOL计算机程序进行求解，得到E偏振情况下2D介质（GaAs)光子晶体能带结构。6.光子能带图（选取能量最低的5个能带） (1)

5、沿布里渊区 10方向() (图2): (2)沿布里渊区 01方向(): (1) (2)7. 电场传播模态Z分量的平面图(1)和高度图(2):A 10,band=1 （平移组合了6个单胞图） (1) (2)B 01,band=5 (1) (2) 8.讨论作者又在文献“Bandgap And Field Propagation Analysis Using Photonic Crystal”（International Journal of Innovative Research in Computer and Communication EngineeringAn ISO 3297: 2007 Certified Organization Vol.3, Special Issue 1, February 2015）中，看到晶体周期单元适合计算机运行的相同本文的图示: 受到启发：计算机看自然可能不会拘于人类目光所决定的一种格式，晶体状态也会有它本身具有的所有高低不同对称性和不同的能量范围的本征态。AlphaGo人机大战挑战人类的思维。如果以往在试题中出现选择矩形为