介质排序及入射角对一维光子晶体禁带特性影响.docVIP

介质排序及入射角对一维光子晶体禁带特性影响 　　摘 要 利用传输矩阵的方法，研究一维光子晶体基本周期的介质排列顺序对光子禁带的影响以及光波入射角度对禁带宽度和禁带中心波长的影响。通过MATLAB仿真，结果表明，在入射角一定的情况下，光子晶体的介质排列顺序对光子禁带的宽度和中心波长没有影响。在入射角从0°-85°的变化过程中，禁带的宽度由514 nm增加到614 nm，禁带的中心位置由1593 nm减小到1411 nm。 　　关键词 一维光子晶体；传输矩阵；基本周期；禁带宽度；禁带中心波长 　　中图分类号：O734 文献标识码：A 文章编号：1671-7597（2013）21-0049-02 　　光子晶体由S.John和E.Yablonovitch独立提出，它是由不同折射率的介质周期性排列而成的人工微结构。在一定频率范围内的光波能通过光子晶体，形成光子导带。另一些频率光波不能通过，形成光子禁带。自从提出光子晶体概念以来，光子晶体的理论得到迅速发展。王旭东等人研究了带隙相对宽度随介质折射率比值变化的规律，认为在结构及折射率一定的情况下，调整介质的大小，不能改变带隙的相对宽度。刘启能研究了一维光子晶体禁带的全貌结构，为设计一光子晶体全能反射器提供了理论依据。舒森等人对一维光子晶体滤波器进行了研究，设计出了水下激光通信的窄带滤光片。 　　笔者是在光子晶体折射率取定的情况下，改变基本周期内介质的排列顺序，进而分析其对一维光子晶体禁带的影响。利用传输矩阵的方法，分析并讨论了光子禁带宽度及禁带中心波长随入射光波角度的变化。 　　1 传输矩阵 　　折射率n1、n2，厚度为a、b的两种介质交替周期排列即可组成一维光子晶体结构，周期厚度即为d=a+b，这一结构与多层介质膜相同，空气折射率为n0=1。 　　根据薄膜光学理论，光在每层介质中的传输特性可用一个的特征矩阵表示，对于第j层介质，其特征矩阵为： 　　（1） 　　式中， 　　（2） 　　nj、dj为该介质层的光学厚度，θj为光线在该介质层中与界面法线方向的夹角，λ为入射光的波长。由两层介质组成的一个基本周期单元，其特征矩阵为： 　　（3） 　　= 　　N个基本周期组成的光子晶体总特征矩阵为： 　　（4） 　　则光子晶体对光波的反射系数为： 　　（5） 　　透射率为： （6） 　　反射率T=O时形成光子禁带。 　　2 数值计算结果与讨论 　　2.1 基本周期内介质排列顺序对光子禁带的影响 　　每个基本由两层介质组成，介质的折射率分别为n1=2.822，n2=1.655，中心波长λ=1550nm。每层介质的光学厚度n1a=n2b=λ/4，周期层数N=20，入射角θ=0°作出两层介质排列顺序不同时一维光子晶体的透射谱（图1（a）、（b））。 　　（a） （b） 　　图1 入射角为0°时两种介质不同排列顺序的透射谱： 　　（a）AB排列；（b）BA排列 　　图1（a）图和（b）中的禁带宽度都是533 nm，禁带的中心位置都在1595.5 nm处。图1表明，在入射角为0°时，改变基本周期内介质的排列顺序，禁带宽度和禁带中心波长并没有发生变化。把入射角改为45°。经过仿真，禁带宽度都是569 nm，禁带的中心位置都在1507.5 nm处。这表明入射角为45°时，改变基本周期内介质的排列顺序不影响禁带宽度和禁带中心波长。经过多次仿真，只要角度取定，改变基本周期内两层介质的排列顺序不影响禁带宽度和禁带中心波长。 　　如果基本周期由A、B、C三层介质组成，介质C的折射率为n3，光学厚度n1a=n2b=n3c=λ/4，在入射角为0°时，分别做出三层介质排列顺序不同时的一维光子晶体的透射谱（图2）。 　　（a） （b） 　　（c） （d） 　　图2 入射谱角为0°时三种介质不同排列顺序的透射谱： 　　（a）ABC排列；（b） BAC排列； （c）CBA排列；（d） CAB排列 　　图2表明，在入射角为0°时，图（a）的禁带上限是2828 nm，下限是1998 nm。图（b）的禁带上限是2830 nm，下限是2000nm。图（b）的禁带上限是2829 nm，下限是1999 nm。图（d）的禁带上限是2830 nm，下限是2000 nm。由于取数存在误差，我们可以认为在入射角一定情况下改变基本周期内三层介质的排列顺序不影响禁带宽度和禁带中心波长。经过多次仿真，无论基本周期的内的介质是多少层，改变其排列顺序，不影响禁带宽度和禁带中心波长。对比图1与图2可以发现，三元光子晶体与二元光子晶体的禁带特征有所不同，波长为300 nm至2000 nm的范围内三元光子晶体具有4个禁带，最宽的禁带是840 nm。而二元光子晶体具有两个禁带，较宽的禁带是530 nm。所