1、.创新实验结题报告题 目：一维光子晶体带隙特性研究 学 院： 自动化工程学院 专 业： 自动化 姓 名： 阚枢 指导教师： 李长红 2010 年 10 月 15 日.关键词： 光子带隙 特征矩阵 规律1 引言光子晶体是一种折射率周期变化的人工微结构材料 ,其典型结构为一个折射率周期变化的三维物体 ,周期为光波长量级. 光子在光子晶体中传播存在光子带隙.，频率落在光子带隙内的电磁波不能在光子晶体中传播,光子晶体的这种特性具有极大的理论价值和潜在的应用前景。在光子晶体中掺杂后 ,会在光子能隙中引入局域模式 ,这 将给激光技术和非线性光学等带来全新的应用 ,如 制作零阈值激光器 、光滤波器、慢光缓存
2、器、慢光传感器等。理论研究发现，对于含有缺陷的一维光子晶体，在光子禁带（PBG ：Photonic Band Gap）的带边和缺陷模对应的频率位置，光的传输具有极低的群速度，Scalorta 等人发现在带边处，光脉冲传输速度可以降低到c/17（c 为真空中光速），大约为 1.76107m/s。光子晶体的理论计算已相对成熟 ,本文旨在应用现有的计算方法，建立一维光子晶体模型并讨论一维光子晶体在不同结构参数和参数下的光学传输特性。与平面波算法相比，传输矩阵法的计算量大为降低。传输矩阵方法可以计算一个有限尺寸光子晶体的反射系数及透射系数；根据光的群速度定义，基于此可以计算光在光子晶体中的传输速度；另
3、外，在光子晶体透射特性分析基础上，还可以推导计算一维光子晶体中光传输时电场在各层中的分布情况，因此本文采用传输矩阵法进行光子晶体规律的探讨。2 一维光子晶体2.1 模型的建立一维光子晶体由两种不同相对介电常量 (a , b ) 、厚度( a , b) 的薄介质层交替排列构成的一维周期性结构 材料. 如图 1 所示 ,空间周期为 d = a + b ,一束频率 为 的光从左向右正入射到图中所示的一维周期 性结构材料中.dA B将光波在介质层中的行 进看作是正向行进电磁波 (下行波) 和反向行进电磁 波 (上行波 ) 的叠加. 介质交界面处的电磁场满足边 界条件. 每一介质层与光波的相互作用可由其
4、矩阵完全决定. 介质层两边的场矢量 E , H , E H 的模可用特征矩阵联系起来 :EMH2.2 数学理论推导为行进一步的计算，下面给出光子晶体的特殊矩阵表达式。如图 2 所示 , E0 , H0表示界面 的na一侧矢量 , E, H表示界面的 nb 一侧的场矢量 , EII，HII 表示界面一侧的场矢量 , 在界面上有入射光波 Ei1、反射光波Er1 、透射光波 Et1以及由介质 nb入射到界面上的光波Er2图 3-2 1D PC 结构中任意层的电磁场传播情况以 TE 模为例展开讨论，对于 TE 模，电场分量 E 垂直入射面，根据电磁场边界条件，界面处电场偏振 E 和磁场分量 H 的切向
5、分量连续。以 En、 Hn表示界面 n 处的电场偏振 E 和磁场偏振 H 的切向分量，因为在同一界面两侧，有En-1=En， Hn-1=Hn。对于界面 I：（3-2）112122,coscs.irtriitiricoss 界面 II 上 EII， HII有类似公式。考察界面 I 上的透射场 Et1与界面 II 上入射场.Ei2：（3-3）()1002 11|xz bziktt iikbi bttEeEe式中， 表示，波矢为 k 的平面波在介质12cos()bzbbbtikhn层中垂直横跨两个界面时的相位差，同样， ，根据bire，由以上各式得：0rHE（3-4）1cossinicobbEHH式
6、中， ，将式（3-4）写成矩阵形式为：20cosbb（3-5）sinsicobbEEHH 则对第 j 层单层介质的传输矩阵为（TE 模）：（3-6）cossinicoj jjjjjM式中，（3-7）20csojjjjjjh以上计算关系是对于 TE 模式的，对 TM 模式可以得到类似于上式的结果，只是在 TM 模式情况下， ，若考虑正入射情况（ =0） ，则 TE02cosii模与 TM 模传输矩阵相同。对于一维周期性结构（图 3-1） ，可逐层应用式（3-6）的传输矩阵方程，对.具有 N 层介质的一维结构：（3-8）112111NNNababNEEMHHEEABMHCD 式中 就是一维光子晶体
7、的传输矩阵。abab分别在第 I 个界面及第 N+1 个界面上列电磁场分量方程，并将求得的EI， HI，E N+1，H N+1代入式（3-8） ，展开可求得透射系数和反射系数：反射系数： （3-9-10111NNriABCDE）反射率： （3-9-2Rr）透射系数： （3-9-3100112tNi NNEABCD）透射率： （3-9-4Tt）若 1D PC 结构置于空气中，则有 N+1=0，利用式（ 3-8）和式（3-9）编程计算可以得到一维光子晶体的反射谱或者透射谱，从而得到其带隙结构。在计算所得反射系数基础上，1D PC 结构内部各层电场分布可由式（3-10）计算：（3-10）2012(,
8、)1()()kkkEzrMr式中， 是从第一层到第 k 层的传输矩阵元素，根据式（3-10）就可以计算kijM出电场在光子晶体内部各层的电场强度。对于一维光子晶体（1D PC）结构，采用含有缺陷层的两种高低折射率材料 A 和 B 交替排列的周期结构，基本周期层（AB）的光学厚度取四分之一波长，缺陷层光学厚度为二分之一波长。在有限一维光子晶体结构中，群折射率定义为：n g=c/vg，可以采用有效折射率方法进行计算：.复折射率：（3-11）()()efnik其中的实部和虚部由透射系数定义：透射系数：（3-12）()()|ite有效折射率实部可写为：（3-13）()cnL虚部为：（3-14）()l|
9、()|kt其中， 为透射率 的复角，L 为有限一维光子晶体结构的总体厚度；()t这样，群折射率可以写为：（3-15）()()gdnn从而，群速度可以通过下式计算：（3-16）gcvdn利用获得的群速度，对于有限长度的一维光子晶体慢光结构，就可以计算光脉冲在其中传输产生的时延。4 数据计算与规律研究 4.1 一维光子晶体周期性情况通过 MATLAB 编程，采用传输矩阵法可以方便的求出一位光子晶体的的特征矩阵 M 中的 A , B , C , D 4 个阵元 , ,即可求出一维光子晶体的带隙结构.通过改变各个参数来谈论光子晶体各参数变对光子晶体禁带宽度和群速度的影响规律。按照上文所述建立一维光子晶
10、体模型，折射率为 a，b 的两种介质交替构成光子晶体。折射率 Na=2.22，Nb=1.41；根据四分之一系统，介质层厚度a=abs(lamd0/4/Na)，b=abs(lamd0/4/Nb)，其中 lamd0=1550e-9，空间周期d=a+b；一维光子晶体的空间周期数为 10，加上插入缺陷层共 21 层介质。则通过编程得到一维光子晶体的带隙结构，如图.图 1 光子晶体带隙附近折射率的变化情况图 2 光子晶体带隙附近群速度的变化情况利用本文讨论光子晶体结构和参数对光子近体带隙特性的研究。4.2 层数 NN 对光子晶体性能的影响光子晶体插入缺陷层产生光子禁带，如果改变光子晶体的层数，如图3、图
11、 4 所示，可以发现随引入缺陷的层数增加，光子晶体禁带的带宽变窄，可以得到的更小的群速度。.图 4 一维光子晶体禁带宽度随 NN 的变化曲线4.3 各介质层折射率对光子晶体性能的影响由 4.2 的结论知光子晶体的带隙特性与总层数 NN 有关，但它们的关系不是线性的，于是我们对光子晶体各介质层对其特性的影响进行了研究。下图是研究引入缺陷层折射率对 Vg 和禁带带宽的影响得出的结果：图 3 一维光子晶体群速度随层数 NN 的变化.图5 一维光子晶体群速度Vg随缺陷层折射率Nd的变化曲线图6 一维光子晶体禁带宽度随Nd的变化曲线由图可知随引入缺陷折射率的增大，光子晶体可以得到更大的群速度，更大的禁带
12、宽度。同样的方法我们可以得到光子晶体介质 a,b 对光子晶体参数的影响情况图7 一维光子晶体禁带宽度随Na的变化曲线.为进一步寻找引起光子晶体带隙特性参数变化的因素，我们进行了|Na|-|Nb|和|Na|/|Nb| 的规律探寻：图 8 一维光子晶体禁带宽度随折射率 Na 的变化曲线图 10 一维光子晶体禁带宽度随折射率 Nb 的变化曲线图 11 一维光子晶体群速度随|Na|-|Nb|的变化曲线|图 9 一维光子晶体群速度随折射率 Nb 的变化曲线.图 12 一维光子晶体禁带宽度随|Na|-|Nb|的变化曲线|图 13 一维光子晶体禁带宽度随|Na|/|Nb|的变化曲线|图 14 一维光子晶体群
13、速度随|Na|/|Nb|的变化曲线|.4.3 介质厚度对光子晶体性能的影响（1）介质 a 的厚度对群速度和带宽的影响Na=1.81;Nb=1.41; NN=42.1 2.15 2.2 2.25 2.3 2.35 2.4 2.45x 10-70.3940.3960.3980.40.4020.4040.4060.4080.41与 与 a与 与 与与与与Vg与 与 a与 与 与 与 与 与 与 与 与 与图 14 一维光子晶体群速度随介质 a 的厚度的变化曲线2.1 2.15 2.2 2.25 2.3 2.35 2.4 2.45x 10-70.090.0920.0940.0960.0980.10.1
14、02与 与 a与 与 与与与与 与 a与 与 与 与 与 与 与 与 与图 14 一维光子晶体带宽随介质 a 的厚度 的变化曲线.由图 13、图 14 可得出，随着介质 a 的厚度的增加，群速度和带宽都减小，都近似成线性。由此可知，得到宽的带宽和得到小的群速度是相互矛盾的，所以，要选择合适的 a 的厚度，从而兼顾带宽和群速度的要求。（2）介质 b 的厚度对群速度和带宽的影响Na=1.81;Nb=1.41 NN=4;2.75 2.8 2.85 2.9 2.95 3 3.05 3.1x 10-70.4060.40650.4070.40750.4080.40850.4090.40950.410.41
15、05与 与 b与 与 与与与与Vg与 与 b与 与 与 与 与 与 与 与 与 与图 15 一维光子晶体群速度随介质 b 的厚度的变化曲线2.75 2.8 2.85 2.9 2.95 3 3.05 3.1x 10-70.0970.0980.0990.10.1010.1020.103与 与 b与 与 与与与与 与 b与 与 与 与 与 与 与 与 与图 16 一维光子晶体带宽随介质 b 的厚度的变化曲线.由图 15、图 16 可得出，随着介质 b 的厚度的增加，群速度呈对数增长，而带宽则近似成线性减小。要得到宽的带宽和小的群速度，就需要减小 b 的厚度。（2）介质 d 的厚度对群速度和带宽的影响
16、Na=1.81;Nb=1.41;NN=43.2 3.25 3.3 3.35 3.4 3.45 3.5 3.55x 10-70.3960.3980.40.4020.4040.4060.408与 与 d与 与 与与与与Vg与 与 d与 与 与 与 与 与 与 与 与 与图16 一维光子晶体群速度随介质d的厚度的变化曲线3.2 3.25 3.3 3.35 3.4 3.45 3.5 3.55x 10-70.09950.10.10050.1010.10150.1020.1025与 与 d与 与 与与与d与 与 与 与 与 与 与 与 与图17 一维光子晶体带宽随介质d的厚度的变化曲线.由图 16、图 17 可得出，随着介质 d 的厚度的增加，群速度和带宽都减小。由此可以得到与 a 厚度变化相同的规律：宽的带宽和小的群速度是相互矛盾的，所以，要选择合适的 d 的厚度，从而兼顾带宽和群速度的要求。5 结论本文介绍了特殊矩阵法，并利用特殊矩阵法对一维光子晶体的带隙特性的影响因素进行了探究，得出了一系列因素对光子晶体特性影响的变化曲线，为今后光子晶体带隙方面的实际应用提供一定的数据和基础。参考文献：1.王辉、李永平物理学报. 20012.赵绪新. . 量子电子学报.20043. 宁学峰 中国科学院上海冶金研究所.20004 太原师范学院学报.2007.35.张曼、潘炜.。2009.6
