资源描述：
《用传输矩阵法计算一维光子晶体的带隙特性研究.doc》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。
1、QINGDAOUNIVERSITY俐新实验结题报告题目：一维光子晶体带隙特性研究学院：自动化工程学院专业：自动化姓名：阚枢指导教师：李长红2010年10月150关键词：光子带隙特征矩阵规律1引言光子晶体是一种折射率周期变化的人工微结构材料，其典型结构为一个折射率周期变化的三维物体，周期为光波长量级.光子在光子晶体中传播存在光子带隙.，频率落在光子带隙内的电磁波不能在光子晶体中传播，光子晶体的这种特性具有极大的理论价值和潜在的应用前景。在光子晶体中掺杂后，会在光子能隙中引入局域模式，这将给激光技术和非线性光学等带来全新的
2、应用，如制作零阈值激光器、光滤波器、慢光缓存器、慢光传感器等。理论研究发现，对于含有缺陷的一维光子晶体，在光子禁带(PBG：PhotonicBandGap)的带边和缺陷模对应的频率位置，光的传输具有极低的群速度，Scalorta等人发现在带边处，光脉冲传输速度可以降低到c/17(c为真空中光速)，大约为1.76x107m/so光子晶体的理论计算已相对成熟，本文旨在应用现有的计算方法，建立一维光子晶体模型并讨论一维光子晶体在不同结构参数和参数下的光学传输特性。与平而波算法相比，传输矩阵法的计算量大为降低。传输矩阵方法可以
3、计算一个有限尺寸光子晶体的反射系数及透射系数；根据光的群速度定义，基于此可以计算光在光子晶体中的传输速度；另外，在光子晶体透射特性分析基础上，还可以推导计算一维光子晶体中光传输时电场在各层中的分布情况，因此本文采用传输矩阵法进行光子晶体规律的探讨。2—维光子晶体2.1模型的建立一维光子晶体由两种不同相对介电常量(£a,£b)、厚度(a,b)的薄介质层交替排列构成的一维周期性结构材料.如图1所示，空间周期为d=a+b,—朿频率为3的光从左向右正入射到图中所示的一维周期性结构材料中.fF_FQ-'d+F)S—nnoe_~
4、10匕lz=bz=0恥一心式屮，心=-Q%=-纟®/cos0(仇=色=為)表示，波矢为k的平而波在介质C~层中垂直横跨两个界而时的相位差，同样，E/E/,根据E,由以上各式得:HA(3-4)E}=cos3hEusin^?Hn%H}=-ir/bsin3bEu+cos3bHugo他式中，%=cos2eb,将式(3-4)写成矩阵形式为:(3-5)(3-6)(3-7)sin^%cos氏则对第j层单层介质的传输矩阵为(TE模):cosJ-sin/一ir/jsin巧cos8}式中,以上计算关系是对于TE模式的，对TM模式可以得到类
5、似丁上式的结果，只是在TM模式情况下，77，匕国三」，若考虑正入射情况(0=0),则TE模也cosp与TM模传输矩阵相同。对于一维周期性结构(图3-1),可逐层应用式(3-6)的传输矩阵方程，对具有N层介质的一维结构：En+i=MABEgH,H亦CDH亦=MaMbMa^MaMbMlM2---MvE/v+iHz式中M=M(MbMa--MaMb就是一维光子晶体的传输矩阵。分别在第I个界面及第N+1个界面上列电磁场分量方程，并将求得的E/,川,Ez叭代入式(3-8),展开可求得透射系数和反射系数：反射系数：E,.、K—91
6、_的0+Bryig-C-Dg(3-9-1)r—A%+Bryig+C+Dt]n+[反射率：R=rr^(3-9-2)透射系数：t_®v+i2%(3-9-3)1—Ejg++C+Dm透射率：T=t・f(3-9-4)若IDPC结构置于空气中，则有加+尸"，利用式(3-8)和式(3-9)编程计算可以得到一维光子晶体的反射谱或者透射谱，从而得到具带隙结构。在计算所得反射系数基础上，1DPC结构内部各层电场分布可由式(3-10)计算：E&,劲=[1+r(0)]M,-〃()[1-心)]妙；(3-10)式中，是从第一层到第k层的传输矩阵
7、元素，根据式(3・10)就可以计算出电场在光子晶体内部各层的电场强度。对丁一维光子晶体(1DPC)结构，采用含有缺陷层的两种高低折射率材料A和B交替排列的周期结构，基本周期层(AB)的光学厚度取四分Z—波长，缺陷层光学厚度为二分Z—波长。在有限一维光子晶体结构中，群折射率定义为：吩c/g,可以采用有效折射率方法进行计算：复折射率：n(,lf(0)=n{a))+ik(co)(3-11)其中的实部和虚部由透射系数定义：透射系数：有效折射率实部可写为:虚部为:(3-13)(3-14)(3-15)k(co)=-—InIt{co
8、)ILeo其中，0(劲为透射率r(効的复角，L为有限一维光子晶体结构的总体厚度；这样，群折射率可以写为：从而，群速度可以通过下式计算:(3-16)利用获得的群速度，对丁有限长度的一维光子晶体慢光结构，就可以计算光脉冲在其中传输产生的时延。4数据计算与规律研究4・1一维光子晶体周期性情况通过MATLAB编程，采用传输矩阵法可以方便的
