Optoelectronics
Vol. 08  No. 03 ( 2018 ), Article ID: 27013 , 10 pages
10.12677/OE.2018.83016

A New Method for Obtaining Multiband Band Gap between Two-Dimensional Photonic Crystals

Yiwan Chen1, Haixia Du2, Zhaorong Chen3

1School of Science, Hubei University of Technology, Wuhan Hubei

2The Rockets Army Engineering University, Xi’an Shaanxi

3Computer and Electronic Information School, Guangxi University, Nanning Guangxi

Received: Sep. 7th, 2018; accepted: Sep. 22nd, 2018; published: Sep. 29th, 2018

ABSTRACT

The two-dimensional tetragonal structure, triangular structure, and hexagonal structure composed of cylindrical dielectric were studied. The refractive index is n = 1.5 (glass), 1.8 (graphite), 2.5 (graphite) respectively. Between the cylinder is tangent or not tangent, the paper got the complete band gap distribution. Through the two different square structures combination, and triangular and hexagonal structure combination, the paper got band gap after combination, and found that after the combination of two different square structures, the gap was fully reproduced. After triangular and hexagonal structure combination, the band gap in triangular structure was completely reappeared; band gap of hexagonal structure was partly reappeared.

Keywords:Photonic Crystals, FDTD, Energy Band, Photonic Band Gap, Crystal Combination

一种获得二维光子晶体多个宽带隙的新方法

陈义万1,杜海霞2,陈昭蓉3

1湖北工业大学理学院,湖北 武汉

2火箭军工程大学,陕西 西安

3广西大学计算机与电子信息学院,广西 南宁

收稿日期:2018年9月7日;录用日期:2018年9月22日;发布日期:2018年9月29日

摘 要

本文系统的研究了由圆柱介质构成的二维正方结构、三角形结构、六角形结构,在折射率分别为n = 1.5 (玻璃),1.8 (石墨),2.5 (石墨),圆柱之间相切和不相切的条件下的带隙,得到了完整的带隙分布;通过把两种不同的正方形结构组合,以及三角形结构和六角形结构组合,得到组合后的带隙,发现两种不同的正方形结构组合后,这两种结构的带隙在组合后都得到完全的再现;三角形和六角形结构组合后,三角形结构的带隙得到完全再现,六角形结构的带隙得到部分再现。

关键词 :光子晶体,FDTD,能隙,光子带隙,晶体组合

Copyright © 2018 by authors and Hans Publishers Inc.

This work is licensed under the Creative Commons Attribution International License (CC BY).

http://creativecommons.org/licenses/by/4.0/

1. 引言

完整的晶体中的原子,由于周期性的排列,使晶体中运动的电子具有能带结构。由介质周期性排列形成的光子晶体对光的作用,类似于晶体对电子的作用,也让光在光子晶体中传播时,在一些波长范围内能够透射,而在另外的波长范围内不能透射。把不能透射的波长范围(也可以折算成对应的能量范围)叫带隙。光子晶体由于其优异的性能,成为光子学、光电材料研究的热门领域 [1] [2] 。研究光在光子晶体中传播性质的方法有平面波展开法,转移矩阵法,有限时域差分法(Finite Difference Time Domain,简称FDTD)。二维光子晶体是由介质柱在平面内周期性排列,在第三维无限延长构成的结构,它相对于三维光子晶体,计算量比较小,便于研究。如何利用二维光子晶体获得多个宽的带隙,使光子晶体的带隙满足工程的需要,是一个令人感兴趣的问题。Xiao-Jing Liu综述了函数光子晶体的带隙性质 [3] ,D. R. Sollia,J. M. Hickmann 研究了带隙作为填充率和折射率函数的性质 [4] ,Xiao-Jing Liu,Jing-Bin Lu,Ji Ma,Zhi-Guo Wang a研究了折射率作为半径函数的光子晶体的带隙的性质 [5] [6] 。本文希望用两种不同的光子晶体获得比较大的带隙的光子晶体。

2. 光子晶体的带隙理论和有限时域差分法

FDTD方法,是由K. S. Yee在1966年提出的 [7] ,其基本思想是把数学中的微分用差分代替,用中心差分法表达导数为 [8] :

d f ( x ) d x = f ( x + Δ x ) f ( x Δ x ) 2 Δ x ,

由于计算机的储存容量的限制,计算的物理区间是有限的,用FDTD方法计算时,被广泛使用的是完全匹配吸收层 [9] (PML),通过选择边界内的电导率和磁导率的比值:

σ σ * = ε 0 μ 0 ,使电磁波从介质入射到边界后,被完全吸收,不被反射。

用FDTD方法推导的2维TE模式的电磁波递推表达式 [10] :

;

E x n + 1 ( i , j ) = ε 0 Δ t σ y ε 0 Δ t + σ y E x n 1 ( i , j ) + 1 Δ y ( ε 0 Δ t + σ y ) [ H z x n ( i , j + 1 ) + H z y n ( i , j + 1 ) H z x n ( i , j 1 ) H z y n ( i , j 1 ) ] ;

H z x n + 1 ( i , j ) = μ 0 Δ t σ x * μ 0 Δ t + σ x * H z x n 1 ( i , j ) 1 Δ x ( μ 0 Δ t + σ x * ) [ E y n ( i + 1 , j ) E y n ( i 1 , j ) ] ;

H z y n + 1 ( i , j ) = μ 0 Δ t σ y * μ 0 Δ t + σ y * H z y n 1 ( i , j ) + 1 Δ y ( μ 0 Δ t + σ y * ) [ E x n ( i , j + 1 ) E x n ( i , j 1 ) ] ;

在完全绝缘介质中,取 σ x = σ y = σ x * = σ y * = 0

3. 由圆柱介质构成的二维光子晶体的透射谱

对于二维正方格子光子晶体的带隙已经有很多研究工作 [11] - [19] ,但是,不同结构的带隙的对比,以及三角形结构和六角形结构的带隙并没有系统的结果,下面分别进行研究。

3.1. 圆柱介质构成的正方形二维光子晶体的带隙

圆柱介质折射率分别设置为n = 1.5,1.8,2.5,介质柱半径设置为 r = 0.2 μm ,介质柱构成正方形结构。研究两种不同的结构:介质柱相切排列,晶格常数 a = 0.4 μm ;介质柱不相切排列, a = 0.6 μm ,整个结构的背景介质是空气,等距离排列成 20 × 30 30 × 30 方阵,整个结构在Z方向没有长度,构成二维晶体。信号源是高斯型,波长范围设置为0.4~2.0 μm。使用PML边界,PML层数为14层,反射率0.0001。结构如图1

(1) a = 0.4 μm (2) a = 0.6 μm

Figure 1. The structure of (1): 20 × 20 and (2): two-dimensional square photonic crystals and light source

图1. (1): 20 × 20 和(2): 30 × 30 方阵二维光子晶体的结构和光波源

光波穿过光子晶体,在检测的一条线上被接收,检测抽样频率点500,计算区域网格最小步长为 2.5 × 10 10 m

后面计算例子的PML层的参数、计算网格参数的设置与此相同,不再赘述。位于光源正前方的接收处,坡印廷矢量与波长之间的关系如图2所示。

把计算结果列为表1

Figure 2. Transmission properties of (a) 20 × 20 and (b) 30 × 30 two-dimensional square photonic crystal to light waves

图2. (a) 20 × 20 和(b) 30 × 30 方阵二维光子晶体对光波的传播性质

Table 1. Transmission properties of two-dimensional square photonic crystals

表1. 二维正方形光子晶体的透射性质

最大带隙分别为0.30 μm和0.26 μm,已经在表中用黑体字表示。

3.2. 圆柱介质构成的三角形二维光子晶体的带隙

圆柱介质折射率分别设置为n = 1.5,1.8,2.5,介质柱半径设置为 r = 0.2 μm ,介质柱构成三角形结构。研究两种不同的结构:介质柱相切排列,晶格常数 a = 0.4 μm ;介质柱不相切排列, a = 0.6 μm ,整个结构的背景介质是空气,等距离排列成 16 × 30 方阵,整个结构在Z方向没有长度,构成2维晶体。信号源是高斯型,波长范围设置为0.4~2.0 μm,使用PML边界,结构如图3

高斯型光波穿过光子晶体,在检测的一条线上被接收。位于光源正前方的接收处,坡印廷矢量与波长之间的关系如图4所示。

紧密排列的三角形晶体,与不相切的三角形晶体,在圆柱介质具有不同的折射率时,对光的透射性质不同,把图4中的带隙数值,总结为表2

3.3. 圆柱介质构成的六角形二维光子晶体(石墨烯结构)的带隙

为了研究石墨烯结构的二维光子晶体的带隙,圆柱介质折射率分别设置为n = 1.5,1.8,2.5,介质柱半径设置为 r = 0.2 μm ,介质柱构成六角形(石墨烯)结构。研究两种不同的结构:介质柱相切排列,晶格常数 a = 0.4 μm ;介质柱不相切排列, a = 0.6 μm ,整个结构的背景介质是空气,等距离排列成 15 × 30

Table 2. Transmission properties of two-dimensional triangular photonic crystals

表2. 二维三角形光子晶体的透射性质

最大带隙分别为0.07 μm和0.30 μm,已经在表中用黑体字表示。

(1) a = 0.4 μm (2) a = 0.6 μm

Figure 3. The structure of 16 × 30 two-dimensional triangle photonic crystals and light source

图3. 16 × 30 二维三角形光子晶体的结构和光波源

Figure 4. Transmission properties of 16 × 30 2-dimensional triangle photonic crystal to light waves

图4. 16 × 30 三角形二维光子晶体对光波的传播性质

阵,整个结构在Z方向没有长度,构成2维六角形晶体。信号源是高斯型,波长范围设置为0.4~2.0 μm,使用PML边界,结构如图5

高斯型光波穿过光子晶体,在检测的一条线上被接收。位于光源正前方的接收处,坡印廷矢量与波长之间的关系如图6所示。

把六角形结构的二维光子晶体的对光的透射结果,总结在表3中。把最大的带隙用黑体字表示。

4. 由两种结构得到二维光子晶体的带隙叠加

在科研和工程实践中,我们往往需要宽度更大,范围更广的带隙,考虑到前面几种二维光子晶体的带隙,我们把具有比较宽大的带隙的结构进行组合,希望得到新的带隙。

4.1. 相切的正方形和不相切的正方形结构的组合

表1数据可以看到,正方形二维结构中,相切的结构带隙最大的是n = 2.5,正方形不相切的结构,带隙最大的也是n = 2.5的结构,现在把这两种结构组合,观察组合后的光子晶体带隙(图7)。

图8可以观察到,两种正方形结构的组合的带隙为λ = 0.5~0.53 μm,0.66~0.70 μm,0.88~0.97 μm,

(1) a = 0.4 μm (2) a = 0.6 μm

Figure 5. The structure of 16 × 30 two-dimensional triangle photonic crystals and light source

图5. 16 × 30 二维三角形光子晶体的结构和光波源

Figure 6. Transmission properties of 15 × 30 2-dimensional hexagonal photonic crystal to light waves

图6. 15 × 30 六角形二维光子晶体(石墨烯结构)对光波的传播性质

Figure 7. A combination of two square structures

图7. 两种正方形结构的组合

Figure 8. Transmission-wavelength relation

图8. 结构光的透射与波长关系

Table 3. Transmission properties of two-dimensional hexagonal photonic crystals (graphene structure)

表3. 二维六角形光子晶体(石墨烯结构)的透射性质

最大带隙分别为0.07 μm和0.30 μm,已经在表中用黑体字表示。

1.11~1.18 μm,1.56~1.75 μm,组合结构的带隙,再现了正方形结构的两种成份的全部的带隙。

4.2. 六角形结构和三角形结构的组合

表2中可以看到,三角形结构中,带隙最大的是折射率n = 2.5,不相切的结构的带隙,带隙宽度

Figure 9. The combination of a triangle and a hexaganal structure

图9. 三角形和六角形结构的组合

Figure 10. The band gap of a two-dimensional photonic crystal with triangular and hexagonal structures

图10. 三角形和六角形结构组合的二维光子晶体的带隙

为0.3 μm;

表3中的数据显示,六角形结构中,带隙最大的是折射率n = 1.8,相切的结构,带隙宽度为0.2 μm,但波长位置不同。现在这把两种结构组合,结构如图9

图10的透射图中可以看到,这种组合结构的带隙在:λ = 0.65~0.72 μm,0.80~0.85 μm,1.02~1.13 μm,1.49~1.61 μm,把组合结构的带隙和原结构的带隙相比较,三角形不相切结构的带隙完全再现,六角形相切结构的带隙λ = 1.50~1.70 μm也得到部分再现。

5. 结论

本文系统的研究了由圆柱介质构成的二维正方结构、三角形结构、六角形结构,在折射率分别为n = 1.5 (玻璃),1.8 (石墨),2.5 (石墨),圆柱之间相切和不相切的条件下的带隙,得到了完整的带隙分布;通过把两种不同的正方形结构组合后,这两种结构的带隙在组合后都得到完全的再现;三角形结构和六角形结构组合,得到组合后的带隙,三角形结构的带隙得到完全再现,六角形结构的带隙得到部分再现。六角形结构的带隙部分再现的原因,还需要进一步研究。通过分别设计两种不同带隙的材料,然后把这两种材料组合,可以得到两种材料的带隙的再现,为寻找多个宽带隙材料提供了一个新思路。

基金项目

湖北省教育科学“十二五”规划重大项目,编号2012A016。

文章引用

陈义万,杜海霞,陈昭蓉. 一种获得二维光子晶体多个宽带隙的新方法
A New Method for Obtaining Multiband Band Gap between Two-Dimensional Photonic Crystals[J]. 光电子, 2018, 08(03): 113-122. https://doi.org/10.12677/OE.2018.83016

参考文献

  1. 1. John, S. (1987) Strong Localization of Photons in Certain Disordered Dielectric Lattice. Physical Review Letters, 58, 2486-2489.

  2. 2. Yablonovitch, E. (1987) Inhibited Spontaneous Emission in Solid State Physics and Electronics. Physical Review Letters, 58, 2059-2061.

  3. 3. Liu, X.J., Liang, Y., et al. (2017) Two-Dimensional Function Photonic Crystals. Physica E, 85, 227-237.

  4. 4. Solli, D.R. and Hickmann, J.M. (2011) Study of the Properties of 2D Photonic Crystal Structures as a Function of the Air-Filling Fraction and Refractive Index Contrast. Optical Materials, 33, 523-526.

  5. 5. Liu, X.J., Lu, J.B., et al. (2017) Adjustable Band Gaps Structures and Semi-Dirac Points in Two-Dimensional Linear Function Photonic Crystals. Optik, 149, 252-260.

  6. 6. Ma, J. and Wang, Z.G., et al. (2017) Photonic Band Gaps Structure Properties of Two-Dimensional Function Photonic Crystals. Physica E, 89, 61-66.

  7. 7. Yee, K.S. (1966) Numerical Solution of Initial Boundary Value Problems involving Maxwell Equations in Isotropic Medie. IEEE Transactions on Antennas Propagate, 14, 302-307.

  8. 8. 葛德彪, 阎玉波. 电磁波时域有限差分方法[M]. 第3版. 西安: 西安电子科技大学出版社, 2011: 10.

  9. 9. 崔俊海, 钟顺时. 采用PML吸收边界条件的FDTD法在分析平面微带结构中的应用[J]. 微波学报, 2000, 16(20): 537-541.

  10. 10. 陈义万, 杜海霞, 陈昭蓉, 闵锐. 二维光子晶体带隙与结构的关系[J]. 湖北大学学报(自然科学版), 2016, 33(3): 220-223.

  11. 11. Inci, H.D. and Ozsoy, S. (2012) A Theoretical Study of Large Solid-Core Square-Lattice Silica Photonic Crystalfibers with Square Air-Holes. Optical Materials, 35, 205-210.

  12. 12. Fan, Q.B., Li, C.Q., et al. (2016) Polarization-Independent Waveguides in Air Holes Photonic Crystals and Its Slow Light. Optics Communications, 380, 227-232.

  13. 13. Inci, H.D. and Ozsoy, S. (2014) Birefringence, Dispersion and Loss Properties for PCFs with Rectangular Air-Holes. Infrared Physics & Technology, 67, 354-358.

  14. 14. Wu, F.G., Wu, T.W., Liu, Z.Y. and Liu, Y.Y. (2006) Effect of Point Defect Geometry on Localized Defect Modes in Two-Dimensional Photonic Crystals. Physics Letters A, 349, 285-290.

  15. 15. 胡学美, 沈昱明, 杨毅彪. FDTD法对二维正方格子光子晶体带隙的数值计算[J]. 光学技术, 2013, 39(3): 256-262.

  16. 16. Chen, H.B., Li, Z.F., Liu, W., et al. (2006) Line Defect Splitters for Self-Collimated Beams in Photonic Crystals. Optics Communications, 262, 120-124.

  17. 17. Shi, P., Huang, K. and Li, Y.P. (2012) Photonic Crystal with Complex Unit Cell for Large Complete Band Gap. Optics Communications, 285, 3128-3132.

  18. 18. Fang, Y.T. and Shen, T.G. (2004) Three Absolute Photonic Gaps from 2d Square Compound Lattices. Crystal Research and Technology, 39, 434-438.

  19. 19. Wang, J.J., Ouyang, Z.B., Wen, G.H. et al. (2014) Wide Ab-solute-Photonic-Bandgap 2D Square-Lattice Photonic Crystal Based on Hollow Cylinders and Cross Connecting Plates. Applied Mechanics and Materials, 670-671, 105-108.

期刊菜单