等离激元增强的石墨烯光吸收*

赵承祥 郄媛 余耀 马荣荣 秦俊飞 刘彦

1) (山西师范大学现代文理学院, 临汾 041000)

2) (山西师范大学物理与信息工程学院, 临汾 041004)

3) (南京信息工程大学滨江学院, 无锡214105)

(2019 年 10 月 28日收到; 2019 年 12 月 23日收到修改稿)

石墨烯中等离激元具有特殊的光电性质, 其和入射光的强烈耦合可以引起光吸收的增强. 本文基于时域有限差分法和多体自洽场理论研究了等离激元对处于光学谐振腔中的石墨烯光吸收的影响. 由于石墨烯中等离激元与入射光动量和能量不匹配而不能直接相互作用, 因此石墨烯上施加了金属光栅结构. 研究发现光栅结构能够对入射光进行动量补偿并且能够引起其下石墨烯中的电场强度产生很大程度增强, 从而导致在该石墨烯结构中太赫兹等离激元和入射光发生强烈耦合而产生太赫兹等离极化激元, 同时引起石墨烯光吸收的增强. 希望本文能够加深对石墨烯光电特性的理解以及可以为基于石墨烯的太赫兹光电装置提供一定的理论依据.

1 引言

石墨烯因具有二维单层、线性能谱、高的载流子迁移率等特性[1], 被认为是比较理想的等离激元载体[2−5]. 相比金属表面等离激元, 石墨烯中的等离激元具有更强局域性、场增强、低的欧姆损耗和高的传播程[5]、易调控性[6]. 此外石墨烯中等离激元的频率是从太赫兹到远红外频段[5], 这些对于石墨烯用于太赫兹(THz)等离激元光电装置提供了条件. 因此近年来, 基于石墨烯优异的等离激元光学性质, 石墨烯等离激元光电装置的研究备受关注[7−11].

本征石墨烯几乎是透明的, 其吸收系数约为2.3%[12]. 而利用石墨烯中等离激元和入射光的强烈耦合可以实现对入射光吸收的增强. 然而, 由于石墨烯中等离激元和自由空间的光场的动量和能量不匹配[13], 因此很难实现石墨烯中等离激元和光场的直接耦合. 要实现等离激元和入射光的耦合需要一定的耦合机制. 在光学谐振腔中可以实现等离激元和入射光的耦合, 但是由于要实现太赫兹波段的等离激元和光的耦合, 谐振腔就必须具有较大的模式体积, 导致等离激元和光子的耦合较弱[14].文献[15]指出, 由于光栅的调制作用, 在光栅条的边缘附近电磁波的电场分量能够产生很大程度的增强, 等离激元和光子的耦合是正比于电场强度的平方的[14], 所以光栅的施加能够进一步增强其下石墨烯中等离激元和光的耦合从而增强光吸收. 而且金属光栅可以对入射光的动量进行弥补, 使入射光子的动量和石墨烯中的等离激元匹配从而使两者可以发生耦合[13]. 因此本文理论研究了施加了光栅的石墨烯光学谐振腔结构中等离激元对石墨烯光吸收的影响.

2 石墨烯结构模型

2.1 理论模型

图1显示的等离激元装置是由施加了光栅的石墨烯层和光学谐振腔两个基本部分组成. 在石墨烯和光栅之间存在一层厚度为20 nm的Al2O3薄膜, 作为上门电压的介电层. 石墨烯下面是一层厚度为300 nm的SiO2层并且放置在厚度为200 µm的Si基底上. 利用光在Al2O3薄膜上界面以及Si下界面的反射作用, 在虚线框内的Al2O3/SiO2/Si结构可以充当光学谐振腔. 而石墨烯作为等离激元的载体. 光学谐振腔内的光子在光学谐振腔中经过多次反射形成驻波, 对于垂直入射的光场, 腔内光子的模式可以由表示, 这里为腔模光子的角频率,是不同谐振腔光子模式的量子数,c是真空中的光速,L是谐振腔的长度,是谐振腔中介质的介电常数, 由于Al2O3和SiO2层的厚度远小于THz光的波长以及Si层的厚度, 所以我们取Si的介电常数

图1 基于石墨烯的装置示意图(从上至下依次是金光栅/Al2O3薄膜/石墨烯/SiO2/Si; d和w分别代表光栅的周期和光栅条的宽度; 虚线框内的结构可以充当谐振腔的作用;L是谐振腔的厚度)Fig. 1. Schematic illustration of the device based on graphene. From top to bottom, there are the gold grating layer, Al2O3 dielectric medium, graphene sheet, and SiO2/Si layer. Here, d and w are respectively the period and the width of the gold strips. The structure sketched by the dotted line can be served as cavity and L is the cavity length.

考虑一束TM极化的远场平面波从光栅入射到硅基底, 利用 CST Microwave Studio软件的时域有限差分法(FDTD)求解器模拟计算得到了仿真模型的透射曲线和吸收曲线. 由于模型的对称性和周期性, 模拟中模拟区域在(x-z)平面的宽度是光栅的一个周期 (6 µm), 深度对应模型的厚度. 模拟区域在y方向的厚度为6 µm. 模拟区域的顶部和底部采用吸收边界条件, 而x方向和y方向应用周期性边界条件. 模拟区域划分为六面体网格, 在网格的全局设置中, 以 1.5 THz 为中心波长, 每个波长范围内划分为25条线, 最长与最短网格线比极限值为50. 在石墨烯层的局部网格处, 将石墨烯层沿厚度方向加密, 网格设为 0.02 nm. 模拟时温度设置为室温300 K, 金光栅的介电函数采用Drude介电函数的形式[16,17]:其中rad/s为金中的等离激元频率,Hz 为衰减项.

2.2 模拟结果

图2 (a)存在 (黑色实线)和不存在 (红色实线)石墨烯时的结构透射谱; (b)存在(黑色实线)和不存在(红色实线)石墨烯时的结构吸收谱Fig. 2. (a) Frequency dependence of the light transmission of the structure with (black solid line) and without graphene (red solid line); (b) frequency dependence of the light absorption of the structure with (black solid line) and without graphene (red solid line).

在图2 (a)和图2 (b)中, 红色实线分别代表模型中不存在石墨烯时的透射谱和吸收谱. 由于只有当入射光在谐振腔中满足共振条件时才会出现透射的峰值. 所以图2(a)中红色曲线的各峰值所对应的频率满足公式即这些频率模式就是谐振腔中的光子模式. 这些模式也可以通过解麦克斯韦方程得到[18]. 另外可以看到, 不存在石墨烯时整个模型几乎是透明的. 而当石墨烯存在时, 通过介电函数来考虑石墨烯对THz入射光的响应, 即这里分别代表高频介电常数、真空介电常数、石墨烯的厚度. 另外

是石墨烯的面电导率[19−21], 式中是入射太赫兹光的角频率,为等离激元的角频率[22], 其中费米能级是石墨烯的费米速度. 由于石墨烯不完全被金属光栅覆盖, 其有效介电常数为[23]其 中是石墨烯完全被金属光栅覆盖情况下的有效介电常数[24], 这里分别是 SiO2[25]和 Al2O3[26]层的介电常数,h是Al2O3层的厚度;是石墨烯完全不被覆盖时的有效介电常数[27]. 这里取石墨烯中电子浓度时等离激元的频率此频率也对应n=5 的腔模光子频率是等离激元的谱宽,为石墨烯中电子的弛豫时间.对于石墨烯, 电子浓度时, 石墨烯中的电子弛豫时间因此这里取从图2的黑色实线可以看到当谐振腔中存在石墨烯时, 透射率在等离激元共振频率的附近, 即的腔模频率附近, 有很明显的下降. 同时, 可以看到透射峰也分裂为两个峰和, 这是由于等离激元和腔模光子的耦合形成两个等离极化激元模式所致. 另外, 从图2(b)黑色实线可以看到, 当存在石墨烯时, 由于等离激元和腔模光子的耦合导致在共振频率附近石墨烯的光吸收产生很大程度增强.

3 等离激元和腔模光子的耦合

3.1 理论和方法

除了可以在光栅上施加门电压调节石墨烯中载流子浓度之外, 光栅还可以对入射光弥补动量另外, 光栅的施加也决定了石墨烯中被激发的等离激元波矢的大小为因此光栅的施加使得入射光可以和石墨烯中的等离激元耦合形成新的元激发模式-等离极化激元.由于光栅的调制作用, 入射光场的电场在光栅条下边缘附近会产生很大程度的增强, 由于等离激元和光子的耦合强度正比于入射光场电场强度的平方[14],所以光栅能够增强入射光和等离激元的耦合.

图3显示的是x方向偏振光垂直入射的情况下由FDTD模拟得到的结果. 黑线(0.75)和红线(0.5)分别显示不同光栅占空比情况下,光栅条下石墨烯层上的电场x和z分量, 其中电场分量按照入射电磁场电场分量做了归一化. 坐标区间对应宽度 为周 期 为µm的光栅条位置. 其他区间是光栅条之间的空间.从图3可以看到, 电场分量在光栅条下边缘附近有很大程度的增强, 这是由于光栅的调制作用, 入射光场和散射光场在光栅下叠加的结果. 另外可以看到, 增加光栅的占空比, 电场分量的增强程度也随之增加. 这意味着电磁场在石墨烯层上的强度和分布可以通过调节光栅占空比来调节.

图3 不同的光栅占空比下, 在石墨烯层中入射电磁波的电场分量 (a) 和 (b)沿着 x 方向的分布Fig. 3. Spatial distribution of the electric field (a) and (b) components of light field along the x direction in graphene sheet for different factor .

得到石墨烯层上的电场分布之后, 可以计算石墨烯中等离激元和腔模光子耦合形成的等离极化激元的模式. 石墨烯中由电子-电子(e-e)相互作用的导带内跃迁通道产生的低频等离激元模式可以避免朗道阻尼而衰减为电子空穴对[22]. 因此本文仅考虑导带内的电子-电子, 电子-光子(e-p)相互作用. 同时考虑光栅上施加正门电压, 这样和光栅下的电磁场耦合的等离激元是由导带中电子集体激发产生. 石墨烯中无质量电子的哈密顿量可以写成

电子和腔模光子之间的耦合势可以从电子-光子相互作用矩阵元的傅里叶变化得到:

在多体自洽场理论下, 存在电子-电子, 电子-光子相互作用的动力学介电常数可以写为[14]

式中

3.2 结果和讨论

图4 的腔模光子 (红色虚线) 和石墨烯中等离激元(红色实线)耦合形成等离极化激元模式(黑色实线)随着电子浓度的变化Fig. 4. The electron density dependence of the plasmon polariton modes (black solid lines) induced by coupling between cavity photons ( , red dotted lines) and plasmon modes in graphene (red solid line).

耦合共振的地方由腔模光子和等离激元的频率决定, 因此可通过调节谐振腔、光栅或者石墨烯的参数来调节, 例如谐振腔的长度、光栅的周期、石墨烯中电子的浓度等. 同样耦合强度也可以通过调节光栅占空比调节石墨烯层电场分量的强度来调节. 这些结果表明此基于石墨烯的装置可以用于太赫兹场效应装置. 另外, 目前还没有发现相关的实验数据验证, 但是我们注意到了关于GaN/AlGaN二维电子气的实验研究[23], 并且本文理论计算结果和此实验研究能够定性地符合. 因此, 希望本文理论结果能够为石墨烯的相关实验提供一定的理论预测和支持.

4 结论

本文理论研究了加光栅的石墨烯谐振腔结构中等离激元对石墨烯光吸收的影响, 并且得到了石墨烯中腔模光子和等离激元耦合形成的等离极化激元模式. 研究发现入射光的电场分量在石墨烯层中光栅条边缘附近产生很大的增强, 增强程度随着光栅占空比的增大变大. 这可以导致入射光和石墨烯中的太赫兹等离激元模式发生强烈的耦合, 形成两支太赫兹等离极化激元模式. 模拟发现在入射光和石墨烯中的太赫兹等离激元共振频率处系统透射谱产生很大程度的下降, 而共振频率处的光吸收大幅度增强. 这些结果意味着石墨烯中等离激元和光的耦合能够引起对太赫兹入射光吸收的大幅度增强, 同时这也为基于石墨烯的太赫兹光电装置提供了一定的理论依据.