左右手材料光子晶体的双重光学滤波功能

许江勇，苏　安，周丽萍，高英俊，谭福奎，唐秀福

(1.兴义民族师范学院 物理与工程技术学院，兴义 562400；2.河池学院 物理与机电工程学院，宜州 546300；3.广西大学 物理科学与工程技术学院，南宁 530004)

引　言

2001年,美国SMITH等人成功研制具有负折射特性的左手材料，BERRIER等人实验验证有效负折射率光子晶体中光通讯波段的负折射现象[10]；之后,LU等人在微波段实现了有效负折射率光子晶体平板的次波长成像实验[11]；MOUSSA等人设计出可发生负折射的三角形晶格微波波段2维光子晶体[12]；ZHANG等人设计出近红外波段的金属双折射率材料光子晶体[13]。LIU等人根据TE模式和TM模式的折射效应设计出双折射偏振光分器2维光子晶体[14];德国的DOLLING等人研制出负折射率的780nm光波段金属材料2维光子晶体[15]。随着研究进展，通过特殊设计的完全由电介质材料组成的左手材料光子晶体也相继实现[16]，如土耳其比尔肯特大学研究小组的白宝石介质棒组成的正方形晶格2维光子晶体，实现了左手材料特性[17-18]。国内先后有浙江大学的HE团队也设计出具有有效负折射率光子晶体的偏振分离器和开口谐振腔[19-20];北京师范大学物理系的ZHANG设计出对TE和TM偏振光具有相同的有效负折射率的金属棒结构2维光子晶体[21];中国科学院物理研究所FENG等人设计出具有有效负折射率的光子晶体[22]。可见，左手材料光子晶体已经从理论研究走向了实验制备，随着制备工艺的不断完善以及对左手材料特殊物理特性研究的深入，左手材料光子晶体应用研究和设计将进入一个新的热潮[23]。

基于此，本文中在构造左右手材料光子晶体模型(HL)mDl(LH)m的基础上，利用传输矩阵法理论，通过计算机数值计算和可视化仿真，对比分析H为左手材料和右手材料时光子晶体传输特性的异同，找出光传输特性及其调制规律，为新型光学滤波器件的研究和设计提供指导。

1　研究模型和方法

构造和研究的左右手材料光子晶体模型为(HL)mDl(LH)m，不考虑色散，其中H介质层可为左手材料或右手材料，当其为右手材料时是氟化镁，折射率nH=1.38，光学厚度nHdH=λ0/4，当其为左手材料时，nH=-1.38，磁导率μL=-1，nHdH=-λ0/4。L和D介质层均为右手材料，其中L介质层为砷化镓，nL=2.35，nLdL=λ0/4，D介质层为碲化铅nD=4.1，nDdD=λ0/2，D介质亦即插入到镜像对称结构光子晶体(HL)m(LH)m中间的缺陷。λ0是与光子晶体禁带中心频率ω0对应中心波长，m和l是光子晶体(HL)m(LH)m与缺陷Dl的排列周期数，在计算中取正整数。

鉴于计算和研究的主要任务是左右手材料光子晶体的透射能带谱，因此研究方法采用相对成熟且比较简便、直观的传输矩阵法[2-7,24-25]。传输矩阵法在很多文献已经有详细的报道，在此不再赘述。

2　计算结果与分析

2.1　m对双重滤波功能的调制

固定光子晶体(HL)mDl(LH)m的其它参量不变，取缺陷介质层D的排列周期数l=1，基本介质单元(HL)m(LH)m的排列周期数m=1，2，3，4，5。由计算软件MATLAB编程，通过计算机计算模拟，即可绘制出H为右手材料(nH=1.38)和左手材料(nH=-1.38)时光子晶体(HL)mD(LH)m的透射能带谱，如图1和图2所示。图中横坐标单位为归一化频率ω/ω0。

Fig.1Tranmission spectrum for photonic crystals (HL)mD(LH)mwithnH=1.38

Fig.2Tranmission spectrum for photonic crystals (HL)mD(LH)mwithnH=-1.38

从图1可见，当H为右手材料(nH=1.38)时，光子晶体(HL)mD(LH)m的透射能带谱由禁带和通带交替排列形成，而且禁带的中心出现一条精细的窄透射峰(缺陷模)，即频率ω/ω0奇数倍处出现窄透射峰，两条窄透射峰之间的频率ω/ω0偶数倍数处则出现通带。考虑图形周期性及文章篇幅，文中各图只绘制两个周期。随着周期数m增大，各禁带中心的透射峰和透射峰之间的通带位置不变，但它们的带宽变窄，而且当m增大到一定数值时通带出现严重的劈裂现象而发展成透射率比较低的透射峰。若以透射峰的半峰全宽(full width at half maximun，FWHM)计量带宽[5-7]，则频率1.00ω/ω0处的透射峰在m为1,2,3,4,5时的带宽ΔW1分别为0.266×10-1ω/ω0,0.025×10-1ω/ω0,5.281×10-4ω/ω0,5.005×10-4ω/ω0,5.000×10-4ω/ω0，而频率2.00ω/ω0处的中心通带对应各m值的带宽ΔW2分别为4.137×10-1ω/ω0,3.106×10-1ω/ω0,2.543×10-1ω/ω0,2.197×10-1ω/ω0,1.967×10-1ω/ω0。可见，频率ω/ω0处的透射峰或通带的带宽均变窄，如图1a～图1e所示。光子晶体的这种光传输特性对窄、宽带双重光学滤波器件的设计具有一定的指导价值。

当nH=-1.38，即H为左手材料时，如图2所示。光子晶体(HL)mD(LH)m的透射能带谱由精细的窄透射峰和发育很完整的通带交替排列形成，整个能带谱可看成是一个带宽很宽的大禁带里交替出现透射峰和透射带，窄透射峰处于禁带的频率ω/ω0奇数倍处，能带处于禁带的频率ω/ω0偶数倍处。随着周期数m增大，大禁带中的透射峰和通带频率位置也保持不变，但它们的带宽变窄。m=1,2,3,4,5时，频率1.00ω/ω0处的透射峰的带宽ΔW1分别为0.273×10-1ω/ω0,0.026×10-1ω/ω0,3.111×10-4ω/ω0,2.510×10-4ω/ω0,2.500×10-4ω/ω0，而频率2.00ω/ω0处通带对应的带宽ΔW2分别为3.508×10-1ω/ω0,3.064×10-1ω/ω0,2.745×10-1ω/ω0,2.530×10-1ω/ω0,2.374×10-1ω/ω0。对比图2和图1，特别值得注意的是，图2中大禁带中的通带随m增大带宽变窄的同时顶端劈裂比较微弱，即相比H为右手材料，当H为左手材料时光子晶体的通带更加完整。因此，由本光子晶体模型设计宽、窄带双重光学滤波器件时，可通过周期数m来调制滤波器各通道的品质，但当模型结构中的H介质层为左手材料时将会得到更好的滤波性能。

2.2　l对双重滤波功能的调制

接着，固定光子晶体(HL)mDl(LH)m的其它参量不变，并固定基本单元介质(HL)m(LH)m的排列周期数m=5，取缺陷介质层D的排列周期数l=1,2,3,4,5，也分别在介质H为右手材料(nH=1.38)和左手材料(nH=-1.38)时模拟出光子晶体(HL)5Dl(LH)5的透射能带谱，如图3和图4所示。

Fig.3Tranmission spectrum for photonic crystals (HL)5Dl(LH)5withnH=1.38

由图3可见，当H为右手材料(nH=1.38)时，光子晶体(HL)5Dl(LH)5透射能带谱中的频率ω/ω0奇数倍处仍然出现单条窄透射峰，同时频率ω/ω0偶数倍处亦出现带宽很宽的通带，但通带发育不完整，2.0ω/ω0频率处两侧很宽的频率范围内通带的透射率很低。当周期数l增大时，频率ω/ω0奇数倍处的窄透射峰透射率不变且保持精细带宽，但频率ω/ω0偶数倍处的通带则快速产生分裂，形成带宽大小不一的小通带，且l越大分裂越严重，当l=3时分裂成3条，当l=5时分裂成7条，由于分裂成的这些通带之间距离很短，因此带宽也越来越小，但这些透射峰仍然发育不完整，没有形成完整分立的窄透射峰，因为它们的底部相连的部分的透射率还不等于零。可见，当H为右手材料时，光子晶体(HL)5Dl(LH)5也可实现宽、窄带双重滤波功能，只是实现带宽滤波的通带发育不全且参差不齐。

Fig.4Tranmission spectrum for photonic crystals (HL)5Dl(LH)5withnH=-1.38

当H为左手材料(nH=-1.38)时，由图4可见，随着周期数l增大，频率ω/ω0偶数倍处的通带带宽变窄，但通带的透射率保持100%不变，而且通带形态比较规整。如l为1,2,3,4,5时，频率2.00ω/ω0处的通带对应各m值的带宽ΔW2分别为2.376×10-1ω/ω0,1.545×10-1ω/ω0,1.177×10-1ω/ω0,0.955×10-1ω/ω0,0.800×10-1ω/ω0。而频率ω/ω0奇数倍处的窄透射峰左右两侧则新增多条精细透射峰，以频率1.00ω/ω0处为例，当l=1时出现1条，当l为2和3时在频率1.00ω/ω0处窄透射峰的两侧新增了2条透射率很高的窄透射峰，l=2时透射率84.63%，l=4时透射率100%；当l=4时频率1.00ω/ω0处两侧又新增了2条，但此时新增的透射峰透射率比较低，只有11.90%，当l=5时透射率达到67.10%，经计算，当l=6时，这2条新增透射峰透射为98%，此时，频率1.00ω/ω0处及两侧即存在5条透射率很高的高品质窄透射峰。

从图4和图3可得，相对H为右手材料(nH=1.38)，当H为左手材料(nH=-1.38)时，通过缺陷自身周期数l调制，光子晶体(HL)5Dl(LH)5可更好地实现宽、窄带双重滤波功能，而且还可以实现多通道窄带滤波功能，这个特性可对设计光学滤波器件具有指导意义。

2.3　m不对称变化对双重滤波功能的调制

在实际设计中，光子晶体模型结构往往不是对称的，当模型对称性遭到破坏时，光子晶体的内部局域电场的量子化程度也将受到破坏，体现在宏观上的透射特性(分立透射峰)将随之改变。因此，为找出不对称因素对该光子晶体双重滤波功能的影响规律，下面以(HL)m和(LH)m的排列周期数不等值变化为例进行研究，即固定(HL)m的排列周期数m=5和缺陷D的排列周期数l=1，而(LH)m的排列周期数取m为1,2,3,4,5依次递增，则可模拟出当介质H为右手材料(nH=1.38)和左手材料(nH=-1.38)时，(HL)5D(LH)m的透射能带谱，如图5和图6所示。

Fig.5Tranmission spectrum for photonic crystals (HL)5D(LH)mwithnH=1.38

如果以Δm=5-m表示光子晶体(HL)5D(LH)m对称结构的破坏程度。Δm称为不对称度[24-25]，则从图5可见，当介质H为右手材料(nH=1.38)时，随着不对称度Δm增大，频率ω/ω0奇数倍处的窄透射峰透射率快速下降，而且频率ω/ω0奇数倍处的各窄透射峰透射率下降的速度相等，如当m为1,2,3,4,5，亦即不对称度Δm为0,1,2,3,4时，频率ω/ω0奇数倍处的窄透射峰透射率为100%,31.72%,3.55%,0.03%和0%。显然，最终频率ω/ω0奇数倍处及其周围将形成透射率为0的全反射带。而随着Δm增大，频率ω/ω0偶数倍处的通带分布频率范围则变窄，即通带有向频率ω/ω0偶数倍处压缩形成规整通带的趋势。当Δm为0,1,2,3,4时，频率ω/ω0偶数倍处通带频率分布范围ΔW2分别为1.968×10-1ω/ω0,2.090×10-1ω/ω0,2.283×10-1ω/ω0,2.640×10-1ω/ω0,4.170×10-1ω/ω0。

Fig.6Tranmission spectrum for photonic crystals (HL)5D(LH)mwithnH=-1.38

如图6所示，当H为左手材料(nH=-1.38)时，随着Δm增大，频率ω/ω0奇数倍处的窄透射峰透射率也下降，而且各频率ω/ω0奇数倍处的窄透射峰透射率下降速度相同，当Δm为0,1,2,3,4时，窄透射峰透射率为100%,31.72%,3.55%,0.03%和0%。显然，最终频率ω/ω0奇数倍处及其周围也将形成透射率为0的全反射带。另外，对比右手材料的情况知，无论H是右手材料还是左手材料，随着不对称度Δm增大，频率ω/ω0奇数倍处的窄透射峰透射率均下降，而且下降的速度相等。对于频率ω/ω0偶数倍处通带，当Δm增大时，通带顶端的劈裂消失，且带宽变窄的同时发展成比较完整的通带，当Δm为0,1,2,3,4时，频率ω/ω0偶数倍处通带频率分布范围ΔW1分别为2.735×10-1ω/ω0,2.473×10-1ω/ω0,2.469×10-1ω/ω0,2.396×10-1ω/ω0,2.136×10-1ω/ω0。

综合图6和图5可得，随着周期不对称度Δm增大，无论H是右手材料还是左手材料，光子晶体(HL)5D(LH)m窄带滤波功能将慢慢弱化，但其保持宽带滤波功能，而且宽带滤波的频率范围随着Δm增大而变窄，同时H在左手材料情况下实现的滤波性能比右手材料情况下的好。即不对称度Δm对光子晶体(HL)5D(LH)m宽、窄带双重或单重滤波功能具有灵敏的调制作用。

3　结　论

通过对光子晶体(HL)mDl(LH)m的光传输特性进行数值计算、模拟，得出研究结论：无论H是右手材料还是左手材料，光子晶体结构模型(HL)mDl(LH)m均可实现窄带、宽带双重光学滤波功能，但在相同调制机制情况下，两者的滤波性能不一样。

(1)当H为左手材料时，周期数m可调制窄带和宽带通道的滤波品质，而当H为右手材料时，周期数m仅对窄带通道的滤波品质具有调制作用。

(2)当H为左手材料时，周期数l不仅可调制宽带通道的滤波品质，而且还可以调制窄带滤波通道的通道数目，而当H为右手材料时，周期数m仅对宽、窄带通道的滤波特性调制作用不明显。

(3)无论H为左手材料或右手材料，周期数不对称度Δm均对窄带滤波通道的透射率具有调制作用，而且两者的调制机制效果相同。而Δm对于宽带滤波通道性能的调制效果，则是左手材料优于右手材料。

可见，对于实现双重光学滤波功能的光子晶体(HL)mDl(LH)m模型，当H为左手材料时，可获得最好的滤波特性，而且这些特性对研究和设计多重、多通道的光学滤波器件具有指导意义。