1. 引言
完整的晶体是由原子在空间周期性排列,电子在晶体中受到周期性的势能的作用,导致电子出现能带。类似的,光子晶体由介质球(三维)和介质柱周期性排列,光在该晶体中传播时,也能够得到类似电子能带的光子带隙 [1] [2],在光子带隙的范围内,光在晶体中不能传播。这样的材料由于对光的传播有频率选择作用,而得到广泛的应用。比如激光腔体对震荡光波需要有强反射,就可以使用光子晶体材料做成激光的振荡腔,把需要增强的光波段约束在激光腔体内部。隐身飞机对雷达的探测波需要有很低的反射,可以通过设计光子晶体结构,让雷达探测波从光子晶体的表面涂层透射而不反射,从而达到极高的隐身效果。光子晶体中研究比较成熟的是常规光子晶体,它的介质球或介质柱的折射率分布是均匀的;相应的,光子晶体的带隙位置和带隙宽度也会是固定的 [3] - [9]。为了按照要求灵活设计光子晶体的带隙,有人提出了函数光子晶体 [10] [11] [12]。所谓函数光子晶体,是指通过介质折射率和光的电场的非线性关系,由电场强度E在介质柱半径方向控制介质柱的折射率大小,即
,从而灵活调控光子晶体的带隙。但是,在实验中沿半径方向控制电场的强度实际上难以实现,笔者提出了沿半径方向阶梯增加折射率的方法来实现函数光子晶体的功能,这样的介质柱,在实际生产中通过晶体生长或者介质的分子束溅射很容易实现。同时,考察设置缺陷介质柱对光子晶体带隙的影响,也是科研中一个有兴趣的问题。
2. 折射率沿介质柱半径方向阶梯增加的二维光子晶体结构
在光子晶体的模拟计算中,使用较多的方法有转移矩阵方法,平面波展开法,有限时域差分方法(FDTD)。FDTD方法是把光在晶体中的传播过程 [13],用离散化方法表达出来,通过迭代计算,模拟光波在晶体中的传播过程。FDTD方法的优点是:对晶体的大小和边界条件要求不高,通过设定合适的计算网格,可以适用各种情况的光子晶体的模拟计算,而且可以实时显示光通过晶体时的能量分布。
由于计算机的储存容量的限制,计算的物理区间是有限的,用FDTD方法计算时,被广泛使用的是
完全匹配吸收层 [14] (PML),通过选择边界内的电导率和磁导率的比值:
,使电磁波从介质入射
到边界后,被完全吸收。
用FDTD方法推导的2维TE模式和TM模式的电磁波递推表达式 [15],经过计算机的逐步计算,模拟光波在光子晶体中的传播情况,从而得出结果。计算过程中,为了适当控制内存的需求和计算时间,需要恰当的划分计算网格,在计算开始前,对涉及到的各个参数进行适当的设置,既要保证计算精度,也要考虑计算时间和计算机内存大小。
介质柱为圆柱形,半径为r,沿半径方向分为均匀的3层到4层,折射率n阶梯增加,圆柱介质柱排列成正方形结构,晶格常数为a,背景材料为空气,折射率为1。缺陷柱有两种:一种是折射率n = 1.0,一种是折射率n = 4.0。光波源是波长范围为0.1~5 μm的平面波,计算时的频率抽样点为200,为了不影响光波的传播,包围计算区域的边界是完全匹配层(PML),通过光子晶体的光波被完全匹配层完全吸收,没有反射。在平面光波沿y方向通过光子晶体后,测量监测点的透过部分,得出光波的透射系数和透射强度。整个模拟实验使用专门的FDTD软件计算。二维光子晶体的结构图见图1。
Figure 1. Photonic crystal structure of square lattice of cylindrical medium with increasing refractive index step by step with radius
图1. 折射率沿半径阶梯增加的圆柱介质正方形晶格光子晶体结构图
3. 折射率沿介质柱半径分3层增加的光子晶体设置的光子晶体的带隙
图2所对应的光子晶体结构与图1中的相同,介质柱半径r = 0.6 μm,晶格常数a = 1.8 μm,折射率沿介质柱半径分别为n = 1.0,n = 2.5,n = 4.0的光子晶体。使用平面波光源,从结果中看到,在1.0~1.2 μm,1.50~1.62 μm,1.80~1.95 μm光的透射能量低,是晶体在Y方向的带隙。
Figure 2. Light wave transmission of photonic crystals with dielectric column radius r = 0.6 μm, lattice constant a = 1.8 μm and refractive index n = 1.5 and n = 2.5, n = 4.0 respectively
图2. 介质柱半径r = 0.6 μm,晶格常数a = 1.8 μm,折射率沿介质柱半径分别为n = 1.0,n = 2.5,n = 4.0的光子晶体的光波透射图
为了对比在相同的晶体结构条件下,引进空气柱缺陷时,光子晶体的带隙的变化,从而发现空气柱缺陷对光子晶体带隙的影响,我们与图1的结构和介质柱分布相同的情况下,在晶体中把一个介质柱换成空气柱,形成一个空气柱缺陷,结构类似于图1。使用相同的平面波光源,计算网格的划分与图1也完全相同,得到如图3的透射曲线。
Figure 3. Light wave transmission of photonic crystals with dielectric column radius r = 0.6 μm, lattice constant a = 1.8 μm and refractive index n = 1.0 and n = 2.5, n = 4.0 respectively, with an air column defect in the middle
图3. 介质柱半径r = 0.6 μm,晶格常数a = 1.8 μm,折射率沿介质柱半径分别为n = 1.0,n = 2.5,n = 4.0,中间设置一个空气柱缺陷的光子晶体的光波透射图
图3的光子晶体结构与图1中的类似,介质柱半径r = 0.60 μm,晶格常数a = 1.8 μm,折射率沿介质柱半径分别为n = 1.0,n = 2.5,n = 4.0,中间设置一个空气柱缺陷的光子晶体。使用平面波光源,从结果中看到,在1.0~1.2 μm,1.50~1.62 μm,1.80~1.95 μm光的透射系数低,是晶体在Y方向的带隙。图3的光子晶体结构与图2的光子晶体结构相同,不同的地方在于在晶体的中间设置了一个空气柱缺陷,但是它们的带隙位置完全相同。
为了对比在相同的晶体结构条件下,引进空气柱缺陷时,光子晶体的带隙的变化,从而发现空气柱缺陷对光子晶体带隙的影响,我们与图1的结构和介质柱分布相同的情况下,在晶体中把两个介质柱换成空气柱,形成两个间隔空气柱缺陷,结构如图4,类似于图1。使用相同的平面波光源,计算网格的划分与图1也完全相同,得到如图5的透射曲线。
Figure 4. The structure of photonic crystals with two separated air defect columns with dielectric column radius r = 0.6 μm, lattice constant a = 1.8 μm, refractive index n = 1.0, n = 2.5, and n = 4.0 along the dielectric column radius
图4. 介质柱半径r = 0.6 μm,晶格常数a = 1.8 μm,折射率沿介质柱半径分别为n = 1.0,n = 2.5,n = 4.0,设置两个间隔空气缺陷柱的光子晶体的结构图
Figure 5. The dielectric column radius r = 0.6 μm, the lattice constant a = 1.8 μm, the refractive index along the dielectric column radius is n = 1.0, n = 2.5, n = 4.0, respectively. The photonic wave transmission pattern of two separated air defect columns is set
图5. 介质柱半径r = 0.6 μm,晶格常数a = 1.8 μm,折射率沿介质柱半径分别为n = 1.0,n = 2.5,n = 4.0,设置两个间隔空气缺陷柱的光子晶体的光波透射图
图5中介质柱半径r = 0.6 μm,晶格常数a = 1.8 μm,折射率沿介质柱半径分别为n = 1.0,n = 2.5,n = 4.0,设置两个间隔空气缺陷柱的光子晶体的光波透射图,从图中可以看到,1.0~1.2 μm,1.50~1.62 μm,1.80~1.95 μm是3个明显的Y方向的带隙。图5的结构与图1的结构相同,不同之处在于图5的结构设置了两个间隔空气缺陷柱,但是带隙的位置没有发生变化,说明介质柱缺陷对光子晶体的带隙没有影响。
为了对比在相同的晶体结构条件下,引进空气柱缺陷时,光子晶体的带隙的变化,从而发现空气柱缺陷对光子晶体带隙的影响,我们与图4的结构和介质柱分布相同的情况下,在晶体中把两个介质柱换成空气柱,形成两个相邻空气柱缺陷,结构如图6,类似于图4。使用相同的平面波光源,计算网格的划分与图4也完全相同,得到如图7的透射曲线。
Figure 6. The structure diagram of the photonic crystals with dielectric column radius r = 0.6 μm, the lattice constant a = 1.8 μm, the refractive index along the dielectric column radius n = 1.0, n = 2.5, n = 4.0, and three connected air defect columns are set
图6. 介质柱半径r = 0.6 μm,晶格常数a = 1.8 μm,折射率沿介质柱半径分别为n = 1.0,n = 2.5,n = 4.0,设置三个相连空气缺陷柱的光子晶体的结构图
Figure 7. Transmission diagram of photonic crystals with three connected air defect columns with dielectric column radius r = 0.6 μm, lattice constant a = 1.8 μm, refractive index n = 1.0, n = 2.5, and n = 4.0, respectively
图7. 介质柱半径r = 0.6μm,晶格常数a = 1.8 μm,折射率沿介质柱半径分别为n = 1.0,n = 2.5,n = 4.0,设置三个相连空气缺陷柱的光子晶体的透射图
从图7中可以看到,在1.0~1.2 μm,1.50~1.62 μm,1.80~1.96 μm有三个明显的Y方向的带隙。尽管图7的结构中介质柱是相邻的,而图5的介质柱缺陷之间间隔一个折射率介质柱,但是它们的带隙完全相同,说明介质柱缺陷的形式对光子晶体带隙没有影响。
Table 1. Comparison of photonic crystal structure and light transmittance with radius r = 0.60 μm and lattice constant a = 1.80 μm
表1. 半径r = 0.60 μm,晶格常数a = 1.80 μm光子晶体结构和光透射率对比
表1中的几种光子晶体的结构,半径相同,r = 0.60 μm,晶格常数相同,a = 1.80 μm,填充率均为0.349,都是分3层折射率沿半径方向阶梯增加。从表1中几种结构的带隙的对比看出:尽管这几种结构的缺陷设置不同,但是它们的带隙位置和带隙宽度都相同。说明,光子晶体带隙的位置和宽度由晶体的结构参数决定,缺陷对带隙几乎没有影响。
4. 折射率沿介质柱半径分4层阶梯增加的光子晶体的带隙
为了对比在不同的晶体结构条件下,引进空气柱缺陷时,光子晶体的带隙的变化,我们改变晶体的晶格常数,设定a = 2.4 μm,介质柱半径r = 0.8 μm,结构类似于图1。使用相同的平面波光源,计算网格的划分与图1也完全相同,得到如图8的透射曲线。
Figure 8. Light wave transmission diagram of a photonic crystal with Radius of dielectric column r = 0.8 μm, lattice constant a = 2.4 μm, refractive index along dielectric column radius is n = 1.0, n = 2.0, n = 3.0, n = 4.0, respectively
图8. 介质柱半径r = 0.8 μm,晶格常数a = 2.4 μm,折射率沿介质柱半径分别为n = 1.0,n = 2.0,n = 3.0,n = 4.0的光子晶体的光波透射图
图8的带隙位置为1.34~1.404 μm,2.0~2.24 μm。
为了对比在相同的晶体结构条件下,引进空气柱缺陷时,光子晶体的带隙的变化,我们与图8对应的结构和介质柱分布相同的情况下,在晶体中把一个介质柱换成空气柱,形成一个空气柱缺陷,结构如图9。类似于图8的结构。使用相同的平面波光源,计算网格的划分与图8也完全相同,得到如图10的透射曲线。
Figure 9. Photonic crystal structure with radius of dielectric column r = 0.8 μm, lattice constant a = 2.4 μm, refractive index along dielectric column radius n = 1.0, n = 2.0, n = 3.0, n = 4.0, respectively. an air defect column is set up
图9. 介质柱半径r = 0.8 μm,晶格常数a = 2.4 μm,折射率沿介质柱半径分别为n = 1.0,n = 2.0,n = 3.0,n = 4.0,设置一个空气缺陷柱的光子晶体的结构图
Figure 10. The light transmission diagram of photonic crystal with an air defect column is set, and radius of dielectric column r = 0.8 μm, lattice constant a = 2.4 μm, refractive index along dielectric column radius n = 1.0, n = 2.0, n = 3.0, n = 4.0, respectively
图10. 介质柱半径r = 0.8 μm,晶格常数a = 2.4 μm,折射率沿介质柱半径分别为n = 1.0,n = 2.0,n = 3.0,n = 4.0,设置一个空气缺陷柱的光子晶体的光透射图
从图10中可以看出,在Y方向,带隙位置为1.34~1.50 μm,2.0~2.2 μm。
为了对比在相同的晶体结构条件下,引进不同空气柱缺陷时,光子晶体的带隙的变化,我们与图9的结构和介质柱分布相同的情况下,在晶体中把2个介质柱换成空气柱,形成两个间隔的空气柱缺陷,结构如图11。类似于图9,使用相同的平面波光源,计算网格的划分与图9也完全相同,得到如图12的透射曲线。
Figure 11. Photonic crystal structure of two air defect columns with dielectric column radius r = 0.8 μm, lattice constant a = 2.4 μm, refractive index n = 1.0, n = 2.0, n = 3.0, n = 4.0, respectively
图11. 介质柱半径r = 0.8 μm,晶格常数a = 2.4 μm,折射率沿介质柱半径分别为n = 1.0,n = 2.0,n = 3.0,n = 4.0,设置2个空气缺陷柱的光子晶体的结构图
Figure 12. Photonic crystal transmittance diagram of two air defect columns with dielectric column radius r = 0.8 μm, lattice constant a = 2.4 μm, refractive index n = 1.0, n = 2.0, n = 3.0, n = 4.0, respectively
图12. 介质柱半径r = 0.8 μm,晶格常数a = 2.4 μm,折射率沿介质柱半径分别为n = 1.0,n = 2.0,n = 3.0,n = 4.0,设置2个空气缺陷柱的光子晶体的光波透射图
从图12可以看出,Y方向带隙的位置是:1.34~1.50 μm,2.0~2.2 μm。
为了对比在相同的晶体结构条件下,引进介质柱缺陷时,光子晶体的带隙的变化,我们与图11的结构和介质柱分布相同的情况下,在晶体中把一个介质柱换成折射率n = 4.0的介质柱,形成一个介质柱缺陷,结构如图13。类似于图11,使用相同的平面波光源,计算网格的划分与图11也完全相同,得到如图14的透射曲线。
Figure 13. Photonic crystal structure diagram of dielectric column radius r = 0.8 μm, lattice constant a = 2.4 μm, refractive index along dielectric column radius n = 1.0, n = 2.0, n = 3.0, n = 4.0, setting 1 defect column n = 4.0
图13. 介质柱半径r = 0.8 μm,晶格常数a = 2.4 μm,折射率沿介质柱半径分别为n = 1.0,n = 2.0,n = 3.0,n = 4.0,设置1个n = 4.0缺陷柱的光子晶体的结构图
Figure 14. Light wave transmission diagram of photonic crystal with dielectric column radius r = 0.8 μm, lattice constant a = 2.4 μm, refractive index along dielectric column radius n = 1.0, n = 2.0, n = 3.0, n = 4.0, set 1 defect column n = 4.0
图14. 介质柱半径r = 0.8 μm,晶格常数a = 2.4 μm,折射率沿介质柱半径分别为n = 1.0,n = 2.0,n = 3.0,n = 4.0,设置1个n = 4.0缺陷柱的光子晶体的光波透射图
从图14可以看出,带隙的位置为:1.34~1.50 μm,1.70~1.72 μm,2.0~2.2 μm。
为了对比在相同的晶体结构条件下,引进不同介质柱缺陷时,光子晶体的带隙的变化,我们与图13的结构和介质柱分布相同的情况下,在晶体中把2个介质柱换成折射率n = 4.0的介质柱,形成间隔介质柱缺陷,结构如图15。类似于图13,使用相同的平面波光源,计算网格的划分与图13也完全相同,得到如图16的透射曲线。
Figure 15. Photonic crystal structure of dielectric column radius r = 0.8 μm, lattice constant a = 2.4 μm, refractive index along the dielectric column radius n = 1.0, n = 2.0, n = 3.0, n = 4.0, and set two n = 4.0 interval defect columns
图15. 介质柱半径r = 0.8 μm,晶格常数a = 2.4 μm,折射率沿介质柱半径分别为n = 1.0,n = 2.0,n = 3.0,n = 4.0,设置2个n = 4.0间隔缺陷柱的光子晶体的结构图
Figure 16. Light wave transmission diagram of a photonic crystal of radius of dielectric column r = 0.8 μm, lattice constant a = 2.4 μm, refractive index along dielectric column radius n =1.0, n = 2.0, n = 3.0, n = 4.0, set 2 columns with n = 4.0 interval defect
图16. 介质柱半径r = 0.8 μm,晶格常数a = 2.4 μm,折射率沿介质柱半径分别为n = 1.0,n = 2.0,n = 3.0,n = 4.0,设置2个n = 4.0间隔缺陷柱的光子晶体的光波透射图
从图16可以看出,沿Y方向的带隙位置为:1.34~1.50 μm,2.0~2.2 μm。
表2中的几种光子晶体的结构,半径相同,r = 0.80 μm,晶格常数相同,a = 1.80 μm,填充率均为0.62,都是分4层折射率沿半径方向阶梯增加。图8是无缺陷结构,其它都是设置缺陷柱的结构。在图14中,多出了1.70~1.72 μm的带隙,说明在折射率沿半径方向阶梯增加的光子晶体结构中,带隙主要是由晶体的结构和填充率决定,缺陷柱的折射率和缺陷柱的数量对带隙没有影响。
Table 2. Comparison of photonic crystal structure and light transmittance with radius r = 0.80 μm and lattice constant a = 2.4 μm
表2. 半径r = 0.80 μm,晶格常数a = 2.4 μm光子晶体结构和光透射率对比
5. 折射率沿介质柱半径分4层阶梯增加的光子晶体和设置缺陷柱的光子晶体的带隙
为了进一步对比在相同的晶体结构条件下,引进不同介质柱缺陷数量和位置时,光子晶体的带隙的变化,我们与图15的结构和介质柱分布相同的情况下,在晶体中把5个介质柱换成折射率n = 1.0的空气介质柱,形成5个间隔空气柱缺陷,结构如图17。类似于图15,使用相同的平面波光源,计算网格的划分与图15也完全相同,得到如图18的透射曲线。
Figure 17. Photonic crystal structure of medium column radius r = 0.8 μm, lattice constant a = 2.4 μm, refractive index along medium column radius n = 1.0, n = 2.0, n = 3.0, n = 4.0, setting 5 air defect columns n = 1.0
图17. 介质柱半径r = 0.8 μm,晶格常数a = 2.4 μm,折射率沿介质柱半径分别为n = 1.0,n = 2.0,n = 3.0,n = 4.0,设置5个n = 1.0空气缺陷柱的光子晶体的结构图
Figure 18. Light wave transmission diagram of photonic crystals with medium column radius r = 0.8 μm, lattice constant a = 2.4 μm, refractive index n = 1.0, n = 2.0, n = 3.0, n = 4.0 along medium column radius respectively, and 5 air defect columns with n = 1.0
图18. 介质柱半径r = 0.8 μm,晶格常数a = 2.4 μm,折射率沿介质柱半径分别为n = 1.0,n = 2.0,n = 3.0,n = 4.0,设置5个n = 1.0空气缺陷柱的光子晶体的光波透射图
从图18中可以看到,沿Y方向的带隙位置为:1.35~1.50 μm,2.0~2.2 μm,2.38~2.48 μm。
为了与图17对比在相同的晶体结构条件下,引进不同介质柱缺陷数量和位置时,光子晶体的带隙的变化,我们与图17的结构和介质柱分布相同的情况下,在晶体中把5个介质柱换成折射率n = 4.0的介质柱,形成5个间隔介质柱缺陷,这5个介质柱排列成5边形,结构如图19。类似于图17,使用相同的平面波光源,计算网格的划分与图17也完全相同,得到如图20的透射曲线。
Figure 19. The structure diagram of photonic crystals with dielectric column radius r = 0.8 μm, the lattice constant a = 2.4 μm, the refractive index along the dielectric column radius n = 1.0, n = 2.0, n = 3.0, n = 4.0, and setting five defect columns n = 4.0
图19. 介质柱半径r = 0.8 μm,晶格常数a = 2.4 μm,折射率沿介质柱半径分别为n = 1.0,n = 2.0,n = 3.0,n = 4.0,设置5个n = 4.0缺陷柱的光子晶体的结构图
Figure 20. Light wave transmission of photonic crystals with medium column radius r = 0.8 μm, lattice constant a = 2.4 μm, refractive index n = 1.0, n = 2.0, n = 3.0, n = 4.0 along medium column radius, respectively. Five defect columns with n = 4.0 are set
图20. 介质柱半径r = 0.8 μm,晶格常数a = 2.4 μm,折射率沿介质柱半径分别为n = 1.0,n = 2.0,n = 3.0,n = 4.0,设置5个n = 4.0缺陷柱的光子晶体的光波透射图
从图20可以看出,沿Y方向的带隙位置为:1.35~1.50 μm,2.0~2.2 μm,2.4~2.5 μm。
表3中的几种光子晶体的结构,半径相同,r = 0.80 μm,晶格常数相同,a = 2.40 μm,填充率均为0.349,都是分4层折射率沿半径方向阶梯增加。从表3中几种结构的带隙的对比看出:
图8是无缺陷结构,图18,图20是有缺陷的结构,它们都有共同的带隙:1.34~1.40,2.0~2.2 μm。图18,图20在无缺陷的结构基础上,增加了一个带隙2.38~2.48 μm,虽然图18和图20的缺陷柱折射率不同,但是它们的带隙位置和宽度很接近。说明缺陷介质柱的折射率对光子晶体带隙影响很小。
Table 3. Comparison of photonic crystal structure and light transmittance with radius r = 0.80 μm and lattice constant a = 2.4 μm
表3. 半径r = 0.80 μm,晶格常数a = 2.4 μm光子晶体结构和光透射率对比
6. 结论
本文对比研究了介质柱折射率沿半径方向阶梯增加的二维光子晶体的带隙,对比了完整光子晶体和设置缺陷柱的光子晶体的带隙,结果表明:光子晶体的带隙,主要由晶体的结构和填充率决定,各种类型的缺陷介质柱对光子晶体的带隙影响很小,这为光子晶体带隙的设计提供了重要参考。