本文设计了圆形、六边形两种散射元并引用二维石墨烯的蜂窝状结构,基于对光子晶体禁带的研究,利用平面波展开法,对这二种散射元构建的蜂窝状光子晶体的禁带特性进行模拟计算。分别得到了在TE和TM两种电磁场模式下光子晶体的禁带宽度,对相对禁带宽度关于散射元结构参数变化的关系做了分析,依据相对禁带宽度的大小,得到了相对禁带最大的模型。 According to the honeycomb structure of the two-dimensional graphene, we designed circular and hexagonal scatterers. On the basis of a general analysis of the photonic crystal bandgaps, using the plane wave expansion method, the characteristics of the band gap of photonic crystal constituted by two kinds of scattering elements with honeycomb structure were calculated. The bandgap of photonic crystals in TE and TM electromagnetic field modes is obtained, respectively. The relationship between the relative bandgap and the structural parameters of scattering elements is analyzed. According to the value of the relative bandgap, the model with the maximum relative band gap is obtained.
于鲁辉1*,杨忠泽2,刘浩楠1,刘可歆1,万勇1#
1青岛大学物理科学学院,山东 青岛
2索斯兰国际学校青岛为明中学,山东 青岛
收稿日期:2019年8月27日;录用日期:2019年9月5日;发布日期:2019年9月12日
本文设计了圆形、六边形两种散射元并引用二维石墨烯的蜂窝状结构,基于对光子晶体禁带的研究,利用平面波展开法,对这二种散射元构建的蜂窝状光子晶体的禁带特性进行模拟计算。分别得到了在TE和TM两种电磁场模式下光子晶体的禁带宽度,对相对禁带宽度关于散射元结构参数变化的关系做了分析,依据相对禁带宽度的大小,得到了相对禁带最大的模型。
关键词 :光子晶体,石墨烯,散射元,禁带
Copyright © 2019 by author(s) and Hans Publishers Inc.
This work is licensed under the Creative Commons Attribution International License (CC BY).
http://creativecommons.org/licenses/by/4.0/
光子晶体 [
石墨烯是一种由碳原子呈蜂窝晶格构成的单层片状结构,光子晶体结构可以模仿石墨烯结构,并人为的进行设计 [
本文创新性地结合石墨烯蜂窝状结构,将其应用于二维光子晶体的禁带研究,分别以圆形和六边形为散射元进行构建,采用平面波展开法,计算了光子晶体在不同结构参数下禁带的变化规律,为光子晶体在通讯和太阳能利用等领域提供了方向,也为石墨烯蜂窝状结构的潜在应用提供了途径。
实验模型中参数设置取值一般与材料、晶格常数和散射元有关 [
由于石墨烯光子晶体中传播的电磁波模式可以分为TE和TM两种模式,因此模拟时可分为四种情况:1) 是TE模介质柱型;2) 是TE模空气孔型;3) 是TM模介质柱型;4) 是TM模空气孔型。四种情况下,都是通过单一的改变半径R或L的参数值,计算相对禁带宽度,分别得到相对禁带随参数变化的规律,找出相对禁带最大值对应的参数值。
实验模拟采用的是通用方法:通过改变参数分别对于圆形散射元和六边形散射元构建的石墨烯结构光子晶体进行模拟计算,得到不同的参数对应的绝对禁带相对禁带的大小,然后经过换算,获得最大的相对禁带宽度 [
图1. 石墨烯结构光子晶体示意图
经模拟计算,得到在TM模式下得到的相对禁带宽度不理想,在TE模式下得到了较大的相对禁带宽度,如图2所示,图2(a)、图2(b)分别表示了在TE模介质柱和TE模空气孔两种情况下相对禁带宽度关于参数R的变化规律。可见,二者变化规律完全不同:
1) TE模介质柱情况下,相对禁带宽度随R呈先增大后减小的趋势,R = 0.24 μm时,相对禁带最大;
2) TE模空气孔情况下,相对禁带宽度随R呈增大趋势,R < 0.28 μm时没有禁带出现,R = 0.50 μm时,相对禁带最大。
图2. (a) TE模介质柱禁带随R参数变化规律;(b) TE模空气孔柱禁带随R参数变化规律
经模拟计算,在TE模空气孔和TM模介质柱两种情况下得到的禁带,二者有相似的地方:在TE模介质柱和TM模空气孔两种情况下得到了较大的的相对禁带宽度,相对禁带宽度随L都呈先增大后减小的趋势。如图3所示,图3(a)、图3(b)分别表示了这两种情况下相对禁带宽度关于参数R的变化规律。但是,而且也有区别:1) TE模介质柱情况下,当L = 0.19 μm时,相对禁带较大;2) TM模空气孔情况下,当L = 0.47 μm时,相对禁带最大。
通过对以上进行的模拟进行总结,可以获得圆形和六边形散射元石墨烯结构四种情况下的相对禁带的最大值对应的结构,具体数据如表1所示。
图3. (a) TE模介质柱禁型带随L参数变化规律;(b) TM模空气孔型禁带随L参数的变化规律
散射元形状 | 圆形 | 六边形 | ||||
---|---|---|---|---|---|---|
散射元类型 | R (μm) | 禁带中心频率(a/λ) | 相对禁带宽度(%) | L (μm) | 禁带中心频率(a/λ) | 相对禁带宽度(%) |
TE模介质柱 | 0.24 | 0.29 | 29.7 | 0.19 | 0.42 | 49.0 |
TE模空气孔 | 50 | 0.14 | 26.7 | 0.57 | 0.55 | 14.6 |
TM模介质柱 | 0.50 | 0.16 | 29.6 | 0.39 | 0.31 | 11.5 |
TM模空气孔 | 0.48 | 0.36 | 14.56 | 0.47 | 0.36 | 53.6 |
表1. 四种情况下最大相对禁带对应的参数
分析模拟实验结果可得:以六边形为散射元的石墨烯结构光子晶体比以圆形为散射元的石墨烯结构光子晶体所得到的相对禁带宽度理想:1) 以六边形为散射元的石墨烯结构光子晶体结构中,TE模介质柱型情况下,L = 0.19 μm时,相对禁带最大,其能带情况如图4(a)所示,其中阴影部分代表禁带,其归一化中心频率为0.29,相对禁带宽度为49.0%;2) TM模空气孔情况下,L = 0.47μm时,相对禁带最大,其能带情况如图4(b)所示,其中阴影部分代表禁带,其归一化中心频率为0.36,相对禁带宽度为53.6%。
图4. 实验模拟的能带图
本文创新性设计了石墨烯蜂窝状结构的二维光子晶体,分别以圆形和六边形为散射元进行构建,采用平面波展开法,计算了光子晶体在结构参数的变化下相对禁带宽度的变化规律。最终得到六边形散射元所构建的更加理想的该光子晶体相对禁带宽度,且在TE模介质柱型和TM模空气孔型两种情况下,相对禁带宽度都呈先增大后减少的趋势,最后得到相对禁带最理想的情况:六边形散射元构建的光子晶体,TE模介质柱型情况下,L = 0.19 μm时,归一化中心频率为0.29,相对禁带宽度为49.0%;TM模空气孔型,L = 0.47 μm时,归一化中心频率为0.36,相对禁带宽度达到最大为53.6%。可见,散射元结构和周期性排列类似时,更容易形成相对宽的禁带,这为制作实际的光子晶体提供了相应的数据以及合理的模型,尤其对太阳能方面的应用具有实际的意义,也为石墨烯结构的应用提供了方向。
于鲁辉,杨忠泽,刘浩楠,刘可歆,万 勇. 石墨烯结构的二维光子晶体禁带特性研究 Research on the Characteristics of Bandgaps of Two-Dimensional Photonic Crystal with Graphene Structure[J]. 应用物理, 2019, 09(09): 379-383. https://doi.org/10.12677/APP.2019.99045