传输矩阵法分析有限周期光子晶体的反射谱特征.docx

A ct a S i nica Q u ant u m O pt ica文章编号 : 1007 - 6654 (2004) 01 - 0038 - 04传输矩阵法分析有限周期光子晶体的反射谱特征张会云1 ,张玉萍2 ,3 ,苏希玉1 ,刘灿德1(1 . 曲阜师范大学物理工程学院 ,山东 曲阜 273165 ; 2 . 天津大学激光与光电子研究所 , 天津 300071 ; 3 . 郑州大学激光与光电信息技术重点实验室 , 河南 郑州 450052)摘要 : 利用传输矩阵法计算了一维光子晶体的光子带隙结构 ,分析了有限周期一维光子晶体的反射率及有缺陷的有限周期一维光子晶体的反射率 ,并讨论了它们的特征 。关键词 : 传输矩阵 ; 光子晶体 ; 光子能带 ; 反射率 ;缺陷中图分类号 : O431文献标识码 : A法4光子晶体属于光学与凝聚态物理交叉的新领域 ,是近年来应用物理学的一个热门研究课题 ,是一 种折射率周期变化的人工微结构材料 ,其典型结构 为一个折射率周期变化的三维物体 ,周期为光波长 量级 。光子晶体的特征是具有光子带隙 ,频率落在 光子带隙内的光在光子晶体中沿任何方向都不能传, 计算了有限周期一维光子晶体的反射率及有缺陷的有限周期一维光子晶体反射率 ,并分析讨论了它们的特征 。1物理模型考虑由两 种 不 同 介 质 构 成 的 一 维 光 子 晶 体 结播1 ,2,因此它具有能够抑制原子分子的自发辐射构 ,厚度分别为 a 和 b , 周期 a + b = d 。折射率分别为 n a 和 nb , 光线如图 1 所示入射到光子晶体上 。 根据文献4 ,电磁波在分层介质中的传输特性可用其传输矩阵表示 。在均匀介质膜的内部 ,假定等诱人的光电子学特性 ,这在基础研究和实际应用上都有着巨大的潜力 。利用光子晶体可以控制原子 的自发辐射 ,可以制作宽带 、低损耗的光反射镜 ,可 以制作零阈值的激光器 、高效率的发光二极管 ,可以 制作光滤波器 、光开关 ,光混频器和光存储器等3 。 光子晶体的理论计算已趋于成熟 ,常用的方法 有平面波展开法 ,时域有限差分法 ,多重散射法 ,传 输矩阵法等 。但如何利用较为简单易行的方法估算 制造参数 、选取介质仍有重要的研究意义 。相对而 言 ,一维光子晶体在结构上最为简单 、易于制备 ,在 计算上也最易于提出模型进行模拟计算 。本文利用 一种相对简单但实用的理论计算方法 ———传输矩阵介电常量ε、磁导率 μ 和折射率 n =εμ, 研究的所有介质都是非磁性的 (μ = 1) 。如果以 θ代表电磁波的传播方向和介质表面法线方向的夹角 ,考虑T E 波 ,其传输矩阵为 :M ( z )- i sin ( k0 n z co sθ)co s ( k0 n z co sθ)pco s ( k0 n z co sθ)=( 1)- i p sin ( k0 n z co sθ)收稿日期 : 2003 - 12 - 15作者简介 : 张会云 (1974 - ) 男 ,山东临沂沂水人 ,山东曲阜师范大学物理工程学院在读研究生.张会云等 : 传输矩阵法分析有限周期光子晶体的反射谱特征·39 ·图 1 一维周期性光子晶体2反射率式中 p1R = |=ε1 / μ1 co sθ1 ,r | ,( 4)ε o sθ。对于一维周期性结构的情况 ,其中 p =μcpl = εl / μl co sθl , ε1 、μ1逐层应用 (1) 式 ,可以写出多层介质的M ( z N )和 εl 、μl 分别为与光子晶体接触的左边和右边介质的介电常量和磁导率 ; θ和θl 分别为光子晶体的入 射波和透射波方向与介质表面法线之间的夹角 。由 (3) 、(4) 式利用 Matlab 数学工程软件编程得到周期 性结构光子晶体的反射率 R 随入射光角频率 ω的变化 关系 ,从而可推知周期性结构光子晶体的带隙结构 。= M 1 ( z 1 ) M 2 ( z 2 -z 1 )M N ( z N-z N - 1 )m 11m 21m 12m 22=( 2)从而进一步得到整个结构的反射系数及反射率分别为 :反射系数2结果与分析2 . 1光子晶体周期性结构情况在计算中假定组成光子晶体的媒质层 a 和 b 分 m 11 p 1 + m 12 p l p1 - m 21 - m 22 pl( 3)r=m 11 p1 + m 12 p p - m- m pl 1 21 22 l(a)( b)N = 1 , na = 1 . 378 , nb = 2 . 356N = 2 , na = 1 . 378 , nb = 2 . 356(c)N