本文系统的研究了由圆柱介质构成的二维正方结构、三角形结构、六角形结构，在折射率分别为n = 1.5 (玻璃)，1.8 (石墨)，2.5 (石墨)，圆柱之间相切和不相切的条件下的带隙，得到了完整的带隙分布；通过把两种不同的正方形结构组合，以及三角形结构和六角形结构组合，得到组合后的带隙，发现两种不同的正方形结构组合后，这两种结构的带隙在组合后都得到完全的再现；三角形和六角形结构组合后，三角形结构的带隙得到完全再现，六角形结构的带隙得到部分再现。 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.
陈义万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 beplay安卓登录
This work is licensed under the Creative Commons Attribution International License (CC BY).
http://creativecommons.org/licenses/by/4.0/
完整的晶体中的原子，由于周期性的排列，使晶体中运动的电子具有能带结构。由介质周期性排列形成的光子晶体对光的作用，类似于晶体对电子的作用，也让光在光子晶体中传播时，在一些波长范围内能够透射，而在另外的波长范围内不能透射。把不能透射的波长范围(也可以折算成对应的能量范围)叫带隙。光子晶体由于其优异的性能，成为光子学、光电材料研究的热门领域 [
FDTD方法，是由K. S. Yee在1966年提出的 [
d f ( x ) d x = f ( x + Δ x ) − f ( x − Δ x ) 2 Δ x ,
由于计算机的储存容量的限制，计算的物理区间是有限的，用FDTD方法计算时，被广泛使用的是完全匹配吸收层 [
σ σ * = ε 0 μ 0 ，使电磁波从介质入射到边界后，被完全吸收，不被反射。
用FDTD方法推导的2维TE模式的电磁波递推表达式 [
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 。
对于二维正方格子光子晶体的带隙已经有很多研究工作 [
圆柱介质折射率分别设置为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. (1)： 20 × 20 和(2)： 30 × 30 方阵二维光子晶体的结构和光波源
光波穿过光子晶体，在检测的一条线上被接收，检测抽样频率点500，计算区域网格最小步长为 2.5 × 10 − 10 m 。
后面计算例子的PML层的参数、计算网格参数的设置与此相同，不再赘述。位于光源正前方的接收处，坡印廷矢量与波长之间的关系如图2所示。
把计算结果列为表1。
图2. (a) 20 × 20 和(b) 30 × 30 方阵二维光子晶体对光波的传播性质
| 圆柱相切( r = 0.2 μm , a = 0.4 μm ) | 圆柱不相切( r = 0.2 μm , a = 0.6 μm ) | |
|---|---|---|
| n = 1.5 (窄带的中心波长，宽带隙的范围) | 0.55, 0.58, 1.05~1.15 | 0.42, 0.48, 0.68~0.75, 1.30~1.46 |
| n = 1.8 (窄带的中心波长，宽带隙的范围) | 0.65, 0.68, 1.2~1.37 | 0.54, 0.58, 0.78~0.88, 1.38~1.58 |
| n = 2.5 (窄带的中心波长，宽带隙的范围) | 0.52, 0.68~0.7, 0.88~0.94, 1.5~1.8 | 0.52~0.56, 0.65~0.8, 0.8~0.9, 0.92~1.18, 1.58~1.75 |
| 最宽的的带隙 | 0.30 μm | 0.26 μm |
表1. 二维正方形光子晶体的透射性质
最大带隙分别为0.30 μm和0.26 μm，已经在表中用黑体字表示。
圆柱介质折射率分别设置为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。
为了研究石墨烯结构的二维光子晶体的带隙，圆柱介质折射率分别设置为n = 1.5，1.8，2.5，介质柱半径设置为 r = 0.2 μm ，介质柱构成六角形(石墨烯)结构。研究两种不同的结构：介质柱相切排列，晶格常数 a = 0.4 μm ；介质柱不相切排列， a = 0.6 μm ，整个结构的背景介质是空气，等距离排列成 15 × 30 方
| 圆柱相切( r = 0.2 μm , a = 0.4 μm ) | 圆柱不相切( r = 0.2 μm , a = 0.6 μm ) | |
|---|---|---|
| n = 1.5 (窄带的中心波长，宽带隙的范围) | 0.5, 1.0 | 0.62, 1.18~1.28 |
| n = 1.8 (窄带的中心波长，宽带隙的范围) | 0.55~0.60, 1.2 | 0.65~0.84, 1.25~1.38 |
| n = 2.5 (窄带的中心波长，宽带隙的范围) | 0.75, 0.85, 1.58~1.65 | 0.65~0.72, 0.8~1.1, 1.52 |
| 最宽的的带隙宽度(μm) | 0.07 μm | 0.30 μm |
表2. 二维三角形光子晶体的透射性质
最大带隙分别为0.07 μm和0.30 μm，已经在表中用黑体字表示。
图3. 16 × 30 二维三角形光子晶体的结构和光波源
图4. 16 × 30 三角形二维光子晶体对光波的传播性质
阵，整个结构在Z方向没有长度，构成2维六角形晶体。信号源是高斯型，波长范围设置为0.4~2.0 μm，使用PML边界，结构如图5。
高斯型光波穿过光子晶体，在检测的一条线上被接收。位于光源正前方的接收处，坡印廷矢量与波长之间的关系如图6所示。
把六角形结构的二维光子晶体的对光的透射结果，总结在表3中。把最大的带隙用黑体字表示。
在科研和工程实践中，我们往往需要宽度更大，范围更广的带隙，考虑到前面几种二维光子晶体的带隙，我们把具有比较宽大的带隙的结构进行组合，希望得到新的带隙。
从表1数据可以看到，正方形二维结构中，相切的结构带隙最大的是n = 2.5，正方形不相切的结构，带隙最大的也是n = 2.5的结构，现在把这两种结构组合，观察组合后的光子晶体带隙(图7)。
从图8可以观察到，两种正方形结构的组合的带隙为λ = 0.5~0.53 μm，0.66~0.70 μm，0.88~0.97 μm，
图5. 16 × 30 二维三角形光子晶体的结构和光波源
图6. 15 × 30 六角形二维光子晶体(石墨烯结构)对光波的传播性质
图7. 两种正方形结构的组合
图8. 结构光的透射与波长关系
| 圆柱相切( r = 0.2 μm , a = 0.4 μm ) | 圆柱不相切( r = 0.2 μm , a = 0.6 μm ) | |
|---|---|---|
| n = 1.5 (窄带的中心波长，宽带隙的范围) | 0.58, 1.39~1.45 | 0.58~0.65, 0.75, 1.14~1.20, 1.82 |
| n = 1.8 (窄带的中心波长，宽带隙的范围) | 0.65, 1.00, 1.50~1.70 | 0.78~0.84, 1.21~1.25, 1.88~1.94 |
| n = 2.5 (窄带的中心波长，宽带隙的范围) | 0.73, 1.32, 1.74~1.90 | 0.90~0.97, 1.10~1.17, 2.05~2.12 |
| 最宽的的带隙宽度(μm) | 0.20 μm | 0.07 μm |
表3. 二维六角形光子晶体(石墨烯结构)的透射性质
最大带隙分别为0.07 μm和0.30 μm，已经在表中用黑体字表示。
1.11~1.18 μm，1.56~1.75 μm，组合结构的带隙，再现了正方形结构的两种成份的全部的带隙。
从表2中可以看到，三角形结构中，带隙最大的是折射率n = 2.5，不相切的结构的带隙，带隙宽度
图9. 三角形和六角形结构的组合
图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也得到部分再现。
本文系统的研究了由圆柱介质构成的二维正方结构、三角形结构、六角形结构，在折射率分别为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