-
本文中采用的数学模型为严格耦合波分析法(rigorous coupled-wave analysis,RCWA),该方法是麦克斯韦方程应用于电磁散射结构的精确解。该方法直接、非迭代、具有确定性,已经广泛用于衍射结构的分析和设计中。利用RCWA方法,透射率T和反射率R可以通过累积前向和后向衍射效率分别计算[13],吸波率A=1-R-T。其中将石墨烯视为介电常数为ε的3维材料,电导率使用σ表示[14]:
$ \begin{gathered} \sigma=\frac{2 \mathrm{j} e^2 k_{\mathrm{B}} \hat{T}}{\pi \hbar^2 \omega} \ln \left[2 \cosh \left(\frac{\mu_{\mathrm{c}}}{2 k_{\mathrm{B}} \hat{T}}\right)\right]+ \\ \frac{e^2}{4 \hbar}\left\{\frac{1}{2}+\frac{1}{\pi} \arctan \left[\frac{\left(\hbar \omega-2\left|\mu_{\mathrm{c}}\right|\right)}{2 k_{\mathrm{B}} \hat{T}}\right]+\right. \\ \left.\frac{\mathrm{j}}{2 \pi} \ln \left[\frac{\hbar \omega-2\left|\mu_{\mathrm{c}}\right|^2+\left(2 k_{\mathrm{B}} \hat{T}\right)^2}{\left(\hbar \omega+2\left|\mu_{\mathrm{c}}\right|\right)^2}\right]\right\} \end{gathered} $
(1) 式中,e为电子电荷;$\hbar $为约化普朗克常数;kB为玻尔兹曼常数;$\hat{T}$为绝对温度;μc为石墨烯化学势,可以通过门电压调节进行改变;ω为入射光角频率。
由电导率可以计算出石墨烯的介电常数,具体表达式如下[15]:
$\varepsilon=1-\frac{\operatorname{Im}[\sigma]}{\omega \varepsilon_0 d_{\mathrm{g}}}+\mathrm{j} \cdot \frac{\operatorname{Re}[\sigma]}{\omega \varepsilon_0 d_{\mathrm{g}}} $
(2) 式中,ε0为空气的绝对介电常数,dg=0.34 nm为单层石墨烯的厚度。
-
图 1为两种F-P谐振完美吸波结构,分别记为结构1和结构2。其中构成光子晶体的1/4波长介质A和B分别为二氧化硅(SiO2)和硅(Si),SiO2的介电常数εA=3.9,厚度$d_{\mathrm{A}}=\lambda_0 / 4 \sqrt{\varepsilon_{\mathrm{A}}}$,Si的介电常数εB=11.9,厚度$d_{\mathrm{B}}=\lambda_0 / 4 \sqrt{\varepsilon_{\mathrm{B}}}$,λ0为中心频率对应波长。结构中的顶部镜子和底部镜子为A和B的堆栈结构,n和m分别为顶部镜子和底部镜子的堆栈的数目。假设中心频率f0=2.5 THz,λ0=c/f0,c为光速,取m=3,n=10,将各参数代入(3)式可得到光子禁带范围Δf为2.062 THz~2.938 THz。
$ \Delta f=\frac{4 f_0}{\pi}\left[\arcsin \left(\frac{\left|\sqrt{\varepsilon_{\mathrm{A}}}-\sqrt{\varepsilon_{\mathrm{B}}}\right|}{\sqrt{\varepsilon_{\mathrm{A}}}+\sqrt{\varepsilon_{\mathrm{B}}}}\right)\right] $
(3) 结构1为双层石墨烯F-P谐振腔结构,如图 1a所示。结构1以SiO2材料、厚度为d1的F-P谐振腔为第1层谐振腔,在其中心嵌入厚度为d2、材料为Si的缺陷腔(第2层谐振腔),石墨烯单层在Si材料缺陷腔中心。结构1的顶部镜子采用Si和SiO2堆栈结构(A/B)n,底部镜子采用(A/B)m结构。
结构2为3层谐振腔吸波体结构,在结构1的基础上,结构2在结构1材料为Si的第2层谐振腔内再次嵌入厚度为d3、材料为SiO2的第3层谐振腔,石墨烯单层保持在第3层SiO2缺陷腔中心位置,第1层腔总厚度为d1,第2层腔总厚度为d2。再进一步优化设计,将结构2的第2层腔、第3层腔在上级腔内的位置设为待优化参数,其中第1层腔顶部厚度为t1d1/2、底部厚度为d1(1-t1/2),第2层腔的顶部厚度为t2d2/2、底部厚度为d2(1-t2/2),如图 1b所示,t1和t2为谐振腔位置的比例因子,且保持3层腔的总厚度为d1+d2+d3不变。最终3层谐振腔结构如图 1b所示。
-
为实现完美吸收,需要对结构的参数进行优化设计。对于结构1设计的双层谐振腔,其优化设计参数如表 1所示。
表 1 结构1的双层F-P谐振腔的优化参数表
Table 1. Optimization parameters of double-layer F-P resonator cavity of structure 1
parameter value μc 0.25 eV τ 0.1 ps d1 0.135λ0 d2 0.275λ0 图 2为结构1的吸收谱、反射谱和传输谱。从图中可以看出,结构1由于双层谐振腔产生的F-P谐振形成了两个完美吸收峰,它们的吸收率分别达到了99.27%和99.57%,谐振点频率分别为2.2114 THz和2.8366 THz,将它们分别标记为模式1和模式2。
图 2 结构1的双层F-P谐振腔的吸收谱、反射谱和传输谱
Figure 2. Absorbance, reflectance and transmittance of double-layer F-P resonance cavity of structure 1
模式1和模式2这两个光子带隙内完美吸波峰的产生得益于谐振腔的加入,为了分析两层谐振腔对两个吸波峰的具体贡献,进一步对其两个谐振频率的电场进行分析。图 3为结构1谐振频率的电场分布图。其中图 3a为模式1,谐振频率为2.2114 THz;图 3b为模式2,谐振频率为2.8366 THz,黑色虚线和白色实线分别标记了石墨烯和F-P谐振腔的位置。从图 3中可以看出,模式1中电场局域在临近石墨烯单层的第2层Si腔的上半腔中,而模式2的电场则局域在第2层谐振腔的下半腔中。两个谐振模式的最强电场并不分别分布于两个谐振腔中,而是分布在同一谐振腔中,验证了多层谐振腔的工作原理并不是从数目上增加从而实现更多的F-P谐振模式。但是多层谐振腔的加入使得入射波在光子禁带范围内达到了两个临界耦合条件,因此谐振模式数量增加。对于结构2设计的3层F-P谐振腔,其优化设计参数如表 2所示。
图 3 结构1的双层F-P谐振腔谐振频率的电场分布
Figure 3. Electric field distribution at resonant frequency of double-layer F-P resonance cavity of structure 1
表 2 结构2的3层F-P谐振腔的优化参数表
Table 2. Optimization parameters of three-layer F-P resonator cavity of structure 2
parameter value μc 0.2 eV τ 0.2 ps d1 0.6λ0 d2 0.135λ0 d3 0.35λ0 t1 0.5 t2 0.5 图 4为结构2的吸收谱、反射谱和传输谱。从图中可以看出,结构2设计的3层F-P谐振腔产生了3个完美吸收峰,它们的吸收率分别达到了97.8%、96.87%和96.54%,谐振点频率分别为2.0782 THz、2.4103 THz和2.8100 THz,将它们分别标记为模式1、模式2和模式3。
图 4 结构2的3层F-P谐振腔的吸收谱、反射谱和传输谱
Figure 4. Absorbance, reflectance and transmittance of three-layer F-P resonance cavity of structure 2
图 5为结构2谐振频率的电场分布图。其中图 5a为模式1,谐振频率为2.0782 THz;图 5b为模式2,谐振频率为2.4103 THz;图 5c为模式3,谐振频率为2.8100 THz。可以看出,模式1的电场主要局域在第1层和第3层腔体上半部分,模式2主要为第1层腔体上部以及石墨烯薄膜附近,模式3则为第1层腔体下部。测量3个模式的半峰全宽(full width at half maxima,FWHM),可求得模式1、模式2、模式3的Q因子分别为60、50和312,该结果与模式3局域电场最强对应,表明局域电场强度与Q值相关,局域电场越强,Q值越高。模式3的局域电场出现位置并不在石墨烯临近电场附近,而是出现在距离石墨烯最远的第1层腔内。
图 5 结构2的3层F-P谐振腔谐振频率的电场分布
Figure 5. Electric field distribution at resonant frequency of three-layer F-P resonance cavity of structure 2
综合分析两种结构,它们在达到近100%吸波时的电场分布(见图 3和图 5)说明了系统吸收并不依赖于石墨烯处的电场。完美吸收产生的本质取决于临界耦合条件,临界耦合发生时最强的局域电场并不一定发生在石墨烯周围。因此在系统中引入单层石墨烯,在禁带范围内并不只有一个谐振模式,谐振模式的个数与在谐振腔内满足临界耦合条件的模式个数有关。通过多层F-P谐振腔结构设计,可以在系统中引入多个满足临界耦合条件的模式,越丰富的谐振腔结构能够在禁带模式内创造越多样的近100%吸收模式。
-
石墨烯作为结构中唯一的吸波材料,对系统的吸波特性有重要的影响,图 6和表 3分别显示了随着石墨烯化学势μc的变化,结构吸波模式和3个模式的中心频率的改变。可以看出,随着μc在0.15 eV~0.25 eV增加,模式1和模式2的谐振频率蓝移,模式3的谐振频率红移,说明多层F-P谐振腔系统中谐振模式均受石墨烯化学势改变的影响。
图 6 结构2的3层F-P谐振腔中,不同石墨烯μc下,吸收率与中心频率的关系
Figure 6. Relationship between absorbance and center frequency under different μc of three-layer F-P resonance cavity of structure 2
表 3 结构2的3层F-P谐振腔中,不同石墨烯μc下,不同吸波模式的中心频率
Table 3. Center frequency of absorption mode with different μc of three-layer F-P resonance cavity of structure 2
μc/eV model 1/THz model 2/THz model 3/THz 0.15 2.0724 2.4033 2.8124 0.20 2.0782 2.4105 2.8100 0.25 2.0862 2.4179 2.8076 光子晶体缺陷的厚度影响吸波模式的偏移,F-P谐振条件可表述为[16]:
$ \operatorname{Re}(k d)=i \pi $
(4) 式中,k为波数,d为谐振腔长度,i为F-P谐振模式的数目。当i一定时,增加d将直接导致谐振波长的增加。
通过改变每层F-P谐振腔的厚度来探究3层谐振腔分别对吸波模式的影响,图 7a、图 7b、图 7c分别表示调整d1、d2、d3对系统吸收谱产生的影响。如图 7所示,3个近100%吸波模式由3个谐振腔厚度共同调节,每个腔体厚度的改变都会改变所有吸波频率,保持其余两个谐振腔厚度不变,依次减少d1、d2、d3的长度,对应的3个吸波模式均发生蓝移,出现了多个吸收模式的联动调节效应。此外,通过改变d2,可以对模式1的吸收率进行调节。d2的改变对3个吸波模式有不同的调节效果,模式1的吸收率随着d2的增加出现明显的衰退,模式2的吸收率对于d2的改变具有鲁棒性,根据这一特性,可以通过调节d2对系统吸收模式的数量进行控制。
基于多层F-P谐振腔的石墨烯吸波体设计
Design of graphene absorber based on multilayer F-P resonator
-
摘要: 为了研究当单层石墨烯位于多层Fabry-Pérot(F-P)谐振腔中时,系统近全吸波模式与多层F-P谐振腔数目之间的关系,同时提高系统对吸波模式的调控能力,采用严格耦合波分析法,对石墨烯多层F-P谐振腔系统的吸波响应进行了研究,分析了临界耦合条件下双层和3层F-P谐振腔结构的光谱响应特征。结果表明,双层和3层F-P谐振腔可调谐近全吸波体,分别形成了两个99%以上和3个96%以上的近全吸波模式; 通过对石墨烯掺杂可以实现对3层F-P谐振腔系统吸波特性的调节,通过改变3层谐振腔的结构可以控制系统吸收模式的数量和相对位置。该研究为可调近全吸波系统引入更丰富的吸波线型。
-
关键词:
- 光学器件 /
- 吸波体 /
- 严格耦合波分析 /
- Fabry-Pérot谐振 /
- 石墨烯
Abstract: To find the connection between the number of near-full absorbing mode and the Fabry-Pérot (F-P) multilayer F-P cavity when single-layer graphene is located in the multilayer F-P cavity and to improve the system's controlling ability to the diversity of absorbing modes, rigorous coupled-wave analysis was adopted. The response characteristics of two-layer and three-layer F-P resonator under critical coupling conditions were analyzed. The results show that perfect absorbing modes of more than 99% and 96% can be respectively formed in two-layer and three-layer F-P cavity systems. Among them, the absorption characteristics of three-layer F-P nested cavity system can be adjusted by doping graphene, and the number and relative positions of absorption modes of three-layer nested cavity can be controlled by changing the structure of three-layer F-P cavity. The system introduces a richer absorbing line type.-
Key words:
- optical devices /
- absorber /
- rigorous coupled-wave analysis /
- Fabry-Pérot resonance /
- graphene
-
表 1 结构1的双层F-P谐振腔的优化参数表
Table 1. Optimization parameters of double-layer F-P resonator cavity of structure 1
parameter value μc 0.25 eV τ 0.1 ps d1 0.135λ0 d2 0.275λ0 表 2 结构2的3层F-P谐振腔的优化参数表
Table 2. Optimization parameters of three-layer F-P resonator cavity of structure 2
parameter value μc 0.2 eV τ 0.2 ps d1 0.6λ0 d2 0.135λ0 d3 0.35λ0 t1 0.5 t2 0.5 表 3 结构2的3层F-P谐振腔中,不同石墨烯μc下,不同吸波模式的中心频率
Table 3. Center frequency of absorption mode with different μc of three-layer F-P resonance cavity of structure 2
μc/eV model 1/THz model 2/THz model 3/THz 0.15 2.0724 2.4033 2.8124 0.20 2.0782 2.4105 2.8100 0.25 2.0862 2.4179 2.8076 -
[1] ANDRYIEUSKI A, LAVRINENKO A V. Graphene metamaterials based tunable terahertz absorber: Effective surface conductivity approach[J]. Optics Express, 2013, 21(7): 9144-9155. doi: 10.1364/OE.21.009144 [2] WU J. Tunable ultranarrow spectrum selective absorption in a graphene monolayer at terahertz frequency[J]. Journal of Physics, 2016, D49(21): 215108. [3] FURCHI M, URICH A, POSPISCHIL A, et al. Microcavity-integra-ted graphene photodetector[J]. Nano Letters, 2012, 12(6): 2773-2777. doi: 10.1021/nl204512x [4] 卞立安, 刘培国, 陈雨薇, 等. 石墨烯介质堆栈提高系统调控Fano共振能力[J]. 激光技术, 2018, 42(2): 187-191. BIAN L A, LIU P G, CHEN Y W, et al. Improvement of system tu-nability for Fano resonance by graphene-dielectric stack[J]. Laser Technology, 2018, 42(2): 187-191(in Chinese). [5] NAIR R R, BLAKE P, GRIGORENKO A N, et al. Fine structure constant defines visual transparency of graphene[J]. Science, 2008, 320(5881): 1308. doi: 10.1126/science.1156965 [6] GRANDE M, VINCENTI M A, STOMEO T, et al. Absorption and losses in one-dimensional photonic-crystal-based absorbers incorporating graphene[J]. IEEE Photonics Journal, 2014, 6(6): 1-8. [7] HSU H T, WU C J. Design rules for a Fabry-Perot narrow band transmission filter containing a metamaterial negative-index defect[J]. Progress in Electromagnetics Research Letters, 2009, 9: 101-107. doi: 10.2528/PIERL09032803 [8] GRANDE M, VINCENTI M A, STOMEO T, et al. Graphene-based perfect optical absorbers harnessing guided mode resonances[J]. Optics Express, 2015, 23(16): 21032-21042. doi: 10.1364/OE.23.021032 [9] SCHNEIDER G J, WATSON G H. Nonlinear optical spectroscopy in one-dimensional photonic crystals[J]. Applied Physics Letters, 2003, 83(26): 5350-5352. doi: 10.1063/1.1636249 [10] SHI B, JIAN Z M, WANG X. Defective photonic crystals with greatly enhanced second-harmonic generation[J]. Optics Letters, 2001, 26(15): 1194-1196. doi: 10.1364/OL.26.001194 [11] BIAN L A, LIU P, LI G, et al. Multi-mode absorption in multi-ca-vity photonic crystal with two graphene monolayers[J]. Superlattices and Microstructures, 2017, 112: 303-310. doi: 10.1016/j.spmi.2017.09.037 [12] BIAN L A, YANG L, LIU P, et al. Controllable perfect absorption in a double-cavity photonic crystal with one graphene monolayer[J]. Journal of Physics, 2017, D51(2): 025106. [13] MOHARAM M G, POMMET D A, GRANN E B, et al. Stable implementation of the rigorous coupled-wave analysis for surface-relief gratings: Enhanced transmittance matrix approach[J]. Journal of the Optical Society of America, 1995, A12(5): 1077-1086. [14] FALKOVSKY L A, VARLAMOV A A. Space-time dispersion of graphene conductivity[J]. The European Physical Journal, 2007, B56(4): 281-284. [15] VAKIL A, ENGHETA N. Transformation optics using graphene[J]. Science, 2011, 332(6035): 1291-1294. doi: 10.1126/science.1202691 [16] LIU J T, LIU N H, LI J, et al. Enhanced absorption of graphene with one-dimensional photonic crystal[J]. Applied Physics Letters, 2012, 101(5): 052104. doi: 10.1063/1.4740261