二维金属光子晶体的带结构的研究.pdfVIP

二维金属光子晶体的带结构研究1 路志刚2宫玉彬魏彦玉王文祥 (电子科技大学物理电子学院，四川成都610054) 摘要：为了研究金属光子晶体的带结构，本文给出了一种新的求解方法，即在传统平面波展开法的基础上， 引入一个新问题，首先论证了“新问题”与“原问题”之间的关系，然后引入辅助函数，将“新问题”的色散关 系等价为一个积分微分本征方程，求解这个本征方程得到“新问题”的带结构，通过“新闯题”的带结构就可以 导出了“原问题”的带结构。最后以正方晶格二维正方金属柱光子晶体为例进行了数值计算，得NT满意的结果。 关健字：带结构，平面波展开法，金属光子晶体 1引言 既Yablonovitch提出光子晶体概念以来【11，近几年它已成为人们研究的热点。众 所周知，它的最大特征是存在完全光子带隙(PBG)，电磁波在周期结构中的传播是 用能带理论来描述的，一般来讲，带结构的精确计算是很困难的，尤其对于金属结构 更是如此。平面波展开法是求解光子晶体带结构的重要方法，对于金属光子晶体，由 于电场的法向分量和磁场的切向分量在交界面上是不连续的，而且金属柱的介电常数 随频率变化，使得传统的平面波方法失效。为了从物理上理解金属光子晶体的带结构， 本文提出了一种新的求解方法，即所谓的以积分方程为基础的混合平面波展开法。下 面给出这种方法的数学推导。 2平面波展开法 考虑二维正方晶格金属柱，如图1所 示，背景媒质为空气，晶格基矢为面和磊， 矩形框区域为WS晶胞，表示为 Q=(a．蟊+％a2：圳≤112}，在单位晶胞 一 内，金属柱区域用D表示，边界用aD表示， 相应的外单位法向向量用元表示。为方便 起见，引入倒格矢量云和反，它们与晶格 Z=x3 基矢的关系为：am·％=2碱∥其中，瓯。。 图1正方晶格单位WS晶胞示意图 是Kronecker符号。为了计算电磁本征模 式，设平面波垂直金属柱所在的轴入射， 对于二维光子晶体，类似于波导的情况可以定义TE模式和TM模式，即磁场平行于 介质柱的TE模式和电场平行于介质柱的TM模式。二维平面内的点用尹=(五，屯)表 示，在E极化时令lf，=Ez，H极化时令y=H：，lf厂满足： V2lf，+卢V=0 (1) a，)． rF糨t卜、 f2a) 妙+：0…．．．(z：，T二、F’ … 丫 ，”。’’、’’I 7。。Z，一一 、 。j OD， (2b) ay+／a疗=0，．．．(五，x2)e (圄吸化) 1基金项目教育部新世纪优秀人才支持项目资助；重点实验室基金(51440020104DZ0202) 399 波数，∞是辐射电磁波的角频率，C是真空光速， Igexp(j／云．尹)具有周期性，在E极 化时，lf，在介质柱的边界上满足狄利克雷边界条件(2a)，在H极化时，则满足诺埃 曼边界条件(2b)，上标“+”表示相关的场量位于边界aD之外。 对于在倒格矢空间中给定的后，方程(1)存在无穷平凡解，表示为p=口(云)，这是 一个多值函数，而这个多值函数p(七)就是所要求解的带结构的数学表示。 在(1)式