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

第 55 卷 第 7 期 2006 年 7 月 物 　理 　学 　报 Vol . 55 ,No . 7 ,July ,2006 ( ) 10003290200655 07 359007 ACTA PHYSICA SINICA 2006 Chin . Phys. Soc . 二维金属光子晶体的带结构研究 路志刚 　宫玉彬 　魏彦玉 　王文祥 ( 电子科技大学物理电子学院 ,成都 　610054) (2005 年 7 月 8 日收到 ;2005 年 12 月 14 日收到修改稿) 　　采用一种新的平面波展开法研究金属光子晶体的带结构 ,即在传统平面波展开法的基础上 ,将“原问题”拓展 , 引入一个“新问题”,通过求解“新问题”得到“原问题”的带结构 ,并论证了它们之间的关系. 为了准确求解“新问 题”,引入辅助函数 ,将其色散关系等价为一个积分微分本征方程 ,求解这个本征方程得到“新问题”的带结构 ,从而 由此导出“原问题”的带结构. 最后 , 以正方晶格二维金属光子晶体为例 ,进行数值计算 ,得到了满意的结果. 关键词 : 金属光子晶体 , 平面波展开法 , 带结构 , 正方晶格 PACC : 7125P , 4270Q 强烈波动 ,为了用平面波展开法比较精确的计算光 1 引 言 子晶体的带结构 ,需要大量的展开项 ,这将给计算带 来硬件和软件上的高要求. 对于金属光子晶体 , 电磁 自Yablonovitch 提出光子晶体概念以来[ 1] ,它已 场的不连续性更为突出 ,而且金属柱的介电常数随 成为人们研究的热点 ,众所周知 ,它的最大特征是存 频率变化 ,使得传统的平面波方法失效. ( ) 为了从物理上理解金属光子晶体 ,本文采用了 在完全光子带隙 PB G . 电磁波在周期结构中的传 播是用能带理论来描述的 ,一般来讲 ,采用能带理论 一种新的方法 ,即所谓的以积分方程为基础的混合 精确计算带结构是很困难的 ,尤其对于金属结构更 平面波展开法 ,求解二维金属光子晶体的带结构. 下 是如此. 目前 , 已经提出了几种分析方法 , 如 N 阶 面首先给出这种方法的数学推导 ,然后以正方晶格 法 、传输矩阵法和平面波展开法等. N 阶法是基于 为例进行数值计算 ,验证所述方法的正确性 。 时域有限差分原理的 ,它的优点是计算量与系统的 大小成正比,并且可以得到场分布 ;传输矩阵法是一 2 平面波展开法 种多散射方法 ,适合于计算复杂晶体结构的透射 、反 射系数 ,对光子晶体物理概念的理解帮助不大. 而平 　　考虑二维晶格金属柱 ,如图 1 所示 ,晶格基矢为 a 和 a ,矩形框区域为 WS 晶胞 , 用 Ω= {αa + 面波展开法则能够从物理概念上解释光子晶体的带 1 2 1 1 αa : α ≤12 }表示 ,在单位晶胞内,金属柱区域 结构 ,是一种重要的方法 ,它要求电磁波在倒格矢空 2 2 i 间展开 ,将麦克斯韦方程转化为一个矩阵本征值问 用 D 表示 ,边界用 D 表示 ,相应的外单位法向向量 题 ,其中