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

简单六方结构二维光子晶体能带的COMSOL模拟北京东之星应用物理研究所伍勇1.引言COMSOL携带的案例库里，其中一篇BandgapAnalysisofaPhotonicCrystal(以下简称Bandgap)对砷化镓简单正方格子2D光子能带进行了完整计算和研究。本文将程序用于简单六方结构，并将结果在此做一介绍。2.关于Floquet(弗洛盖）波矢Fk这是入门COMSOL光子晶体能带模拟的重要概念，在另一案例PorousAbsorber中，在Floquet周期性边界条件一段写明：)dk(ie)dx(p)x(p由此我判断Floquet波矢就是Bloch（布洛赫）波矢,但“帮助”文档中有：)sinancosa(sinkk21211F，以正格子基矢21a,a表示（其文没有任何几何插图和物理说明），使我决定必须在六方格子中选择矩形单胞作为周期单元，以使计算机程序能够运行我的几何方案。3.几何建模图1作为试探选择的几何模型，圆形柱代表以GaAs作为格点材料，在空气介质中周期性排列，形成二维六方结构人造晶体。a是晶格常数。zˆ是z方向的单位矢量形单胞六方格子光子晶体的矩图.1以上根据倒格子基矢定义计算出1b，2b及其分量。由倒格子基矢1b，2b，构建长方格子的布里渊区也是长方结构如图2：a3aKMxkykaa1aa32Saakji)aa(aaabyx10002222321321)aa(aaab3211322里渊区六方结构光子晶体的布图2.4.二维光子晶体主方程COMSOL在Bandgap“模型开发器”[电磁波，频域]写出方程形式如下：0)()(0201EjkErr，在Bandgap中，下面目录[波方程，电]中直接简化为，020Ek)E(r电磁波在光子晶体中的传播遵从麦克斯韦方程，上述方程可由麦克斯韦方程组出发导出介质中的麦克斯韦方程组)(D1)(B30)(tBE2)(tDJH4ED，HB，EJ在电介质中一般认为自由电荷，自由电流密度（电导率）为零。本文档不考虑磁性质，0，0J，1r传播模态电场函数COMSOL表达为：)(tiezzike)y,x(E)t,z,y,x(E5，在周期结构中，它应具有Bloch波的性质，不考虑衰减损耗。注意这里次上标的符号与我们习惯的教科书里正负符号相反。由(2)两端取旋度，将（4）代入得：22tE)E(，22k绝缘介质，2202020022nkckrrrr，rrvcn00折射率，1r2nr，0000ck为真空波数，应用（5），如期得到前面提到的Bandgap中的假设波方程：0220Enk)E(利用EE)E()E(22E偏振时（又称横磁场模式TM）对应的二维光子晶体主方程是：)()y,x(zEkn)y,x(zE)yx(602022222这里电场方向平行于z轴，而磁场方向只分布在x-y平面内。二维光子晶体主方程)(6是计算光子能带的出发点和源依据。5.折射率色散函数)f(n)f(n是COMSOL计算能带的第二个关键点，描述折射率随频率参量f的变化。格点材料砷化镓的)f(n，在Bandgap有如下定义n_GaAs=3.3285e5[s/m]*c_const/f+3.5031，参变量f,单位Hz；下限：0，上限：1]s/m[const_c299792458表示真空光速，我在软件Mathematica上化简这个函数，给出f/^*..GaAs_n131097884950313，并绘制图3如下：0.00.20.40.60.81.02.010144.010146.010148.010141.010151.21015fn_GaAs的折射率随频率的变化光子晶体图GaAs3.经COMSOL计算机程序进行求解，得到E偏振情况下2D介质（GaAs)光子晶体能带结构。6.光子能带图（选取能量最低的5个能带）(1)沿布里渊区[10]方向(K)(图2):(2)沿布里渊区[01]方向(M):(1)(2)7.电场传播模态Z分量)y,x(zE的平面图(1)和高度图(2):A[10],band=1（平移组合了6个单胞图）(1)(2)B[01],band=5(1)(2)8.讨论作者又在文献“BandgapAndFieldPropagationAnalysisUsingPhotonicCrystal”（InternationalJournalofInnovativeResearchinComputerandCommunicationEngineeringAnISO3297:2007CertifiedOrganizationVol.3,SpecialIssue1,February2015）中，看到晶体周期单元适合计算机运行的相同本文的图示:受到启发：计算机看自然可能不会拘于人类目光所决定的一种格式，晶体状态也会有它本身具有的所有高低不同对称性和不同的能量范围的本征态。AlphaGo人机大战挑战人类的思维。如果以往在试题中出现选择矩形为六方晶格的布拉菲格子单胞的答案，一定不给分数的，但自然本身却有更大的包容和可能。致谢百度文库，百度学术，道客巴巴的文献支持及小木虫生动深入有益的学术讨论。作者仅在此抛砖引玉。