留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

复合石墨烯/硅半球的宽带太赫兹超材料吸收器

孟宪睿 张铭 席宇鹏 王如志 王长昊 王波

孟宪睿, 张铭, 席宇鹏, 王如志, 王长昊, 王波. 复合石墨烯/硅半球的宽带太赫兹超材料吸收器[J]. 红外与激光工程, 2022, 51(6): 20210648. doi: 10.3788/IRLA20210648
引用本文: 孟宪睿, 张铭, 席宇鹏, 王如志, 王长昊, 王波. 复合石墨烯/硅半球的宽带太赫兹超材料吸收器[J]. 红外与激光工程, 2022, 51(6): 20210648. doi: 10.3788/IRLA20210648
Meng Xianrui, Zhang Ming, Xi Yupeng, Wang Ruzhi, Wang Changhao, Wang Bo. Wideband terahertz metamaterial absorber for composite graphene/silicon hemispheres[J]. Infrared and Laser Engineering, 2022, 51(6): 20210648. doi: 10.3788/IRLA20210648
Citation: Meng Xianrui, Zhang Ming, Xi Yupeng, Wang Ruzhi, Wang Changhao, Wang Bo. Wideband terahertz metamaterial absorber for composite graphene/silicon hemispheres[J]. Infrared and Laser Engineering, 2022, 51(6): 20210648. doi: 10.3788/IRLA20210648

复合石墨烯/硅半球的宽带太赫兹超材料吸收器

doi: 10.3788/IRLA20210648
详细信息
    作者简介:

    孟宪睿,男,硕士生,主要从事光电信息与能源材料方面的研究

  • 中图分类号: O436

Wideband terahertz metamaterial absorber for composite graphene/silicon hemispheres

  • 摘要: 提出了一种性能可调的宽带、极化与入射角不敏感的超材料太赫兹吸收器,该吸收器自上而下分为四层结构,分别是:硅半椭球/半球体复合结构、连续石墨烯层、PDMS介质层和金属背板。通过在TE波垂直入射条件下仿真,在已知结果基础上,对不同石墨烯化学势和不同结构条件下的电场结果分析表明,在硅半椭球/半球体亚波长复合结构所形成的连续、多模法布里-珀罗共振,以及由石墨烯所激发的多个离散的等离子体共振的协同作用下,其吸收光谱得到平滑和扩展,使该结构可实现吸收率宽范围可调,以及接近100%吸收率的宽频带吸收特性。特别的,当石墨烯化学势分别为0.2与0.9 eV时,其分别可获得约5.7 THz与7 THz的宽带太赫兹波吸收(吸收率超过90%),且其最大吸收率接近完美吸收(约99.8%)。此外,该结构还具有360°极化不敏感和高于60°的入射角不敏感等优异特性,在以上角度范围内,吸收器吸收率仍可保持到90%以上。在太赫兹波探测、光谱成像以及隐身技术等方面具有潜在的应用前景。
  • 图  1  (a)宽带吸收器周期结构x-y平面示意图;(b)吸收器的三维单元示意图,其中P = 104 μm, R1= 11 μm, R2= 38 μm, R3= 26 μm, ts = 15 μm, tm = 2 μm

    Figure  1.  (a) x-y plane diagram of the periodic structure of the broadband absorber;(b) 3 D schematic diagram of the absorber, where P=104 μm, R1=11 μm, R2= 38 μm, R3=26 μm, ts =15 μm, tm =2 μm

    图  2  不同石墨烯化学势μ(0~0.9 eV)条件下吸收器吸收率随频率变化图

    Figure  2.  Absorption spectra of the proposed absorber with various values of the graphene chemical potential μ from 0 to 0.9 eV

    图  3  石墨烯层上方z = 50 nm处吸收器在不同石墨烯化学势μ下电场分布图:(a)μ = 0.1 eV;(b) μ = 0.2 eV;(c) μ = 0.4 eV

    Figure  3.  The electric field distributions (|E|) of the proposed absorber on the xy plane with z = 50 nm at the different graphene chemical potential μ of (a) 0.1, (b) 0.2 and (c) 0.4 eV

    图  4  μ = 0.2 eV时五种不同结构的吸收谱:(1)无石墨烯无谐振结构;(2)有石墨烯且只有半球结构(蓝线);(3)有石墨烯且只有半椭球结构(绿线);(4)无石墨烯但具有半椭球/半球谐振结构(粉线);(5)有石墨烯与半椭球/半球谐振结构(红线)

    Figure  4.  (a) Simulated absorption spectra of the proposed G-SemiEllip/SemiSphere absorber (red line) and four other absorbers with NN (black line), G-SemiSphere (blue line), G-SemiEllip (green line) and N-SemiEllip/SemiSphere (pink line)

    图  5  在 (a) 0.82 THz和(b) 4.95 THz的不同频率下,在z =50 nm的x-y平面上吸收器的电场分布 (|E|)

    Figure  5.  The electric field distributions (|E|) of the proposed absorber on the x-y plane with z = 50 nm at the different frequencies of (a) 0.82 THz, and (b) 4.95 THz

    图  6  (a) TE波斜入射时吸收率随极化角ϕ(电场方向与y轴的夹角)的变化率; (b)~(c) TE/TM波斜入射时吸收率随入射角θ(入射方向与z轴的夹角)的变化率

    Figure  6.  (a) The change rate of absorptivity with polarization angle (angle between electric field direction and y-axis) at oblique incident of TE wave; (b)-(c) The rate at which the absorptivity varies with the incidence angle θ (the angle between the incident direction and the z-axis) for oblique incidence of TE/TM waves

  • [1] O’Hara J F, Withayachumnankul W, Al-Naib I. A review on thin-film sensing with terahertz waves [J]. Journal of Infrared Millimeter & Terahertz Waves, 2012, 33(3): 245-291.
    [2] 刘欣缘, 曾昊旻, 田昕, 等. 太赫兹辐射在皮肤中的传输仿真与安全性分析[J]. 光学精密工程, 29(5): 9.

    Liu Xinyuan, Zeng Hanwen, Tian Xi, et al. Transmission simulation and safety analysis of terahertz radiation in skin [J]. Optics and Precision Engineering, 2021, 29(5): 9. (in Chinese)
    [3] Chen G, Lin X, Wang Z. Enhanced reflective dichroism from periodic graphene ribbons via total internal reflection [J]. Opt Express, 2019, 27(16): 22508-22521. doi:  10.1364/OE.27.022508
    [4] Kong X, Wang Z, Du L, et al. Optically transparent metamirror with broadband chiral absorption in the microwave region [J]. Opt Express, 2019, 27(26): 38029-38038. doi:  10.1364/OE.27.038029
    [5] 王 花, 孙晓红, 王 真等. 跟踪微分器在太赫兹波超材料吸波体的特性分析[J]. 红外与激光工程, 2016, 45(12): 1225003. doi:  10.3788/IRLA201645.1225003

    Wang Hua, Sun Xiaohong, Wang Zhen, et al. Characteristic analysis of metamaterial absorber in terahertz wavelength [J]. Infrared and Laser Engineering, 2016, 45(12): 1225003. (in Chinese) doi:  10.3788/IRLA201645.1225003
    [6] Forouzeshfard M R, Ghafari S, Vafapour Z. Solute concentration sensing in two aqueous solution using an optical metamaterial sensor [J]. Journal of Luminescence, 2021, 230: 117734. doi:  10.1016/j.jlumin.2020.117734
    [7] 徐进, 李荣强, 蒋晓云等. 基于方形开口环的超宽带线性极化转换器[J]. 物 理 学 报, 2019, 68(11): 117801. doi:  10.7498/aps.68.20190267

    Xu J, Li R-Q, Jiang X-P, et al. Ultra-wideband linear polarization converter based on square split ring [J]. Acta Physica Sinica, 2019, 68(11): 117801. (in Chinese) doi:  10.7498/aps.68.20190267
    [8] Shen G, Zhang M, Ji Y, et al. Broadband terahertz metamaterial absorber based on simple multi-ring structures [J]. AIP Advances, 2018, 8(7): 075206. doi:  10.1063/1.5024606
    [9] Landy N I. A perfect metamaterial absorber [J]. Physical Review Letters, 2008, 100: 207402. doi:  10.1103/PhysRevLett.100.207402
    [10] Baqir M A, Naqvi S A. Electrically tunable terahertz metamaterial absorber comprised Cu/Graphene strips [J]. Plasmonics, 2020, 15(6): 2205-2211. doi:  10.1007/s11468-020-01252-1
    [11] 王鑫, 孙立宁, 施云波等. 微型多参数水质传感器芯片的设计及应用[J]. 光学精密工程, 2020, 28(10): 2215-2226. doi:  10.37188/OPE.20202810.2215

    Wang Xin, Sun Lining, Shi Yunbo. Design and application of miniature multi-parameter water quality sensor chip [J]. Optics and Precision Engineering, 2020, 28(10): 2215-2226. (in Chinese) doi:  10.37188/OPE.20202810.2215
    [12] Phon R, Lim S. Design and analysis of active metamaterial modulated by RF power level [J]. Sci Rep, 2020, 10(1): 8703. doi:  10.1038/s41598-020-65318-0
    [13] 李泉, 刘姗姗, 路光达等. 利用石墨烯-金属复合结构实现太赫兹电磁诱导透明超表面主动调控[J]. 红外与激光工程, 2020, 49(9): 20210246. doi:  10.3788/IRLA20210246

    Li Quan, Liu Shanshan, Lu Guangda, et al. Active control of terahertz electromagnetically induced transparency metasurface using a graphene-metal hybrid structure [J]. Infrared and Laser Engineering, 2020, 49(9): 20210246. (in Chinese) doi:  10.3788/IRLA20210246
    [14] Zhong R, Yang L, Liang Z, et al. Ultrawideband terahertz absorber with a graphene-loaded dielectric hemi-ellipsoid [J]. Opt Express, 2020, 28(20): 28773-28781. doi:  10.1364/OE.401069
    [15] Andryieuski A, Lavrinenko A V. Graphene metamaterials based tunable terahertz absorber: effective surface conductivity approach [J]. Opt Express, 2013, 21(7): 9144-9155. doi:  10.1364/OE.21.009144
    [16] Zhang Y, Feng Y, Zhu B, et al. Graphene based tunable metamaterial absorber and polarization modulation in terahertz frequency [J]. Opt Express, 2014, 22(19): 22743-22752. doi:  10.1364/OE.22.022743
    [17] 曹暾, 刘宽, 李阳等. 可调谐光学超构材料及其应用[J]. 中国光学, 2021, 14(4): 968-985. doi:  10.37188/CO.2021-0080

    Cao Tun, Li Kuan, Li Yang, et al. Tunable optical metamaterials and their applications [J]. Chinese Optics, 2021, 14(4): 968-985. (in Chinese) doi:  10.37188/CO.2021-0080
    [18] 陈勰宇, 田震. 石墨烯太赫兹波动态调制的研究进展[J]. 中国光学, 2017(1): 86-97. doi:  10.3788/CO.20171001.0086

    Chen Xieyu, Tian Zhen. Recent progress in terahertz dynamic modulation based on graphene. [J]. Chinese Optics, 2017, 10(1): 86-97. (in Chinese) doi:  10.3788/CO.20171001.0086
    [19] Hwang E H, Das Sarma S. Dielectric function, screening, and plasmons in two-dimensional graphene [J]. Physical Review B, 2007, 75(20): 205418. doi:  10.1103/PhysRevB.75.205418
    [20] Yang J, Zhu Z, Zhang J, et al. Broadband terahertz absorber based on multi-band continuous plasmon resonances in geometrically gradient dielectric-loaded graphene plasmon structure [J]. Sci Rep, 2018, 8(1): 3239. doi:  10.1038/s41598-018-21705-2
    [21] Rahmanzadeh M, Rajabalipanah H, Abdolali A. Multilayer graphene-based metasurfaces: Robust design method for extremely broadband, wide-angle, and polarization-insensitive terahertz absorbers [J]. Appl Opt, 2018, 57(4): 959-968. doi:  10.1364/AO.57.000959
  • [1] 王茂旭, 于优, 汤振华, 肖永川, 王超, 高晖.  小型化高集成宽带光收发组件 . 红外与激光工程, 2021, 50(6): 20200513-1-20200513-7. doi: 10.3788/IRLA20210513
    [2] 吴俊政, 倪维平, 严卫东, 张晗.  圆周阵列太赫兹干涉成像中目标场景仿真 . 红外与激光工程, 2019, 48(1): 125004-0125004(8). doi: 10.3788/IRLA201948.0125004
    [3] 张钊, 陈勰宇, 田震.  使用太赫兹快速探测器测量硅少数载流子寿命 . 红外与激光工程, 2019, 48(9): 919003-0919003(6). doi: 10.3788/IRLA201948.0919003
    [4] 潘奕, 郑渚, 丁庆, 姚勇.  宽带太赫兹偶极子光电导接收天线研究 . 红外与激光工程, 2019, 48(1): 125002-0125002(6). doi: 10.3788/IRLA201948.0125002
    [5] 夏祖学, 刘发林, 邓琥, 陈俊学, 刘泉澄.  频率可调太赫兹微结构光电导天线 . 红外与激光工程, 2018, 47(5): 520002-0520002(7). doi: 10.3788/IRLA201847.0520002
    [6] 吴俊政, 严卫东, 倪维平, 张晗.  圆周阵列太赫兹干涉成像仿真 . 红外与激光工程, 2017, 46(8): 825002-0825002(8). doi: 10.3788/IRLA201746.0825002
    [7] 马学, 李琦, 鲁建业.  太赫兹高斯光束整形环形光束 . 红外与激光工程, 2017, 46(5): 525002-0525002(8). doi: 10.3788/IRLA201746.0525002
    [8] 王俊龙, 杨大宝, 邢东, 梁士雄, 张立森, 赵向阳, 冯志红.  0.2THz宽带非平衡式倍频电路研究 . 红外与激光工程, 2017, 46(1): 106003-0106003(4). doi: 10.3788/IRLA201746.0106003
    [9] 王俊龙, 杨大宝, 邢东, 梁士雄, 张立森, 赵向阳, 徐鹏, 冯志红.  0.22 THz宽带混频电路研究 . 红外与激光工程, 2017, 46(11): 1125003-1125003(4). doi: 10.3788/IRLA201746.1125003
    [10] 葛栋森, 许全, 魏明贵, 张学迁, 韩家广.  基于曲折线型介质超材料的宽带太赫兹四分之一波片 . 红外与激光工程, 2017, 46(9): 921002-0921002(5). doi: 10.3788/IRLA201746.0921002
    [11] 娄淑琴, 院楚君, 王鑫.  全光纤同带泵浦宽带掺镱超荧光光纤光源的实验研究 . 红外与激光工程, 2016, 45(8): 802001-0802001(6). doi: 10.3788/IRLA201645.0802001
    [12] 侯宇, 杨会静.  垂直双空芯宽带太赫兹偏振分离器 . 红外与激光工程, 2016, 45(12): 1225005-1225005(5). doi: 10.3788/IRLA201645.1225005
    [13] 王启超, 汪家春, 赵大鹏, 林志丹, 苗雷.  太赫兹波对烟幕的透射能力研究 . 红外与激光工程, 2015, 44(12): 3696-3700.
    [14] 许文忠, 钟凯, 梅嘉林, 徐德刚, 王与烨, 姚建铨.  太赫兹波在沙尘中衰减特性 . 红外与激光工程, 2015, 44(2): 523-527.
    [15] 何晓阳, 张屹遐, 杨春, 陈琦.  太赫兹光子晶体光纤与天线设计 . 红外与激光工程, 2015, 44(2): 534-538.
    [16] 王蓉蓉, 吴振森, 张艳艳, 王明军.  太赫兹波段信号在雾中的传输特性研究 . 红外与激光工程, 2014, 43(8): 2662-2667.
    [17] 鞠智鹏, 李德华, 周薇, 马建军, 李乾坤, 屈操.  相位阶跃变化型太赫兹波带片 . 红外与激光工程, 2013, 42(6): 1519-1522.
    [18] 刘晓旻, 李苏贵, 弓巧侠, 鲁旭, 马省, 梁二军, 李新建.  多孔硅薄膜对p 型单晶硅太赫兹波段透射特性的影响 . 红外与激光工程, 2013, 42(5): 1236-1240.
    [19] 周亚训, 徐星辰.  C+L 波段宽带增益平坦铋基掺铒光纤放大器的设计 . 红外与激光工程, 2012, 41(8): 2119-2124.
    [20] 林栩凌, 周峰, 张建兵, 戴志敏.  采用飞秒装置的高功率宽带太赫兹源 . 红外与激光工程, 2012, 41(1): 116-118.
  • 加载中
图(6)
计量
  • 文章访问数:  299
  • HTML全文浏览量:  89
  • PDF下载量:  43
  • 被引次数: 0
出版历程
  • 收稿日期:  2021-09-08
  • 修回日期:  2021-11-04
  • 刊出日期:  2022-07-05

复合石墨烯/硅半球的宽带太赫兹超材料吸收器

doi: 10.3788/IRLA20210648
    作者简介:

    孟宪睿,男,硕士生,主要从事光电信息与能源材料方面的研究

  • 中图分类号: O436

摘要: 提出了一种性能可调的宽带、极化与入射角不敏感的超材料太赫兹吸收器,该吸收器自上而下分为四层结构,分别是:硅半椭球/半球体复合结构、连续石墨烯层、PDMS介质层和金属背板。通过在TE波垂直入射条件下仿真,在已知结果基础上,对不同石墨烯化学势和不同结构条件下的电场结果分析表明,在硅半椭球/半球体亚波长复合结构所形成的连续、多模法布里-珀罗共振,以及由石墨烯所激发的多个离散的等离子体共振的协同作用下,其吸收光谱得到平滑和扩展,使该结构可实现吸收率宽范围可调,以及接近100%吸收率的宽频带吸收特性。特别的,当石墨烯化学势分别为0.2与0.9 eV时,其分别可获得约5.7 THz与7 THz的宽带太赫兹波吸收(吸收率超过90%),且其最大吸收率接近完美吸收(约99.8%)。此外,该结构还具有360°极化不敏感和高于60°的入射角不敏感等优异特性,在以上角度范围内,吸收器吸收率仍可保持到90%以上。在太赫兹波探测、光谱成像以及隐身技术等方面具有潜在的应用前景。

English Abstract

    • 太赫兹波是一种频率范围位于微波与红外波段之间的电磁波,因其在军事、通信、生物与医学检测以及安保等领域的巨大应用潜力而日渐获得广泛的关注。但目前,针对太赫兹波的研究与应用仍然是所有电磁波中最不发达的频率范围之一,因而被称为太赫兹间隙[1-2]。电磁超材料的出现为填补这一空白开辟了新途径,它是一种由亚波长结构组成的周期性阵列,因其独特的电磁特性备受关注。与其他材料不同,超材料无需改变材料本身,只需通过设计不同结构和尺寸来达到操控电磁波,进而对THz波产生响应的目的[3-4]。利用超材料构造的各种太赫兹功能器件已越来越多地见诸于报道[5],例如超材料传感器[6]、极化转换器[7]与吸收器[8]等。2008年,Landy等提出了基于超材料的完美吸收器概念,他们利用开口环等构成的电磁谐振器可对特定频率电磁波实现完美吸收[9]。然而早期的超材料吸收器由于其共振结构具有强烈的色散,仅能够实现很窄的吸收带宽。为了解决这一问题,研究者们提出了例如多共振器级联以及色散控制等多种方法去扩展超材料吸收器的带宽[10]。不过,基于传统超材料的电磁波吸收器仍然存在着吸收光谱不可调节的局限性。近十年来,研究者们相继提出了利用微机电系统[11]和肖特基二极管结构[12]等方法实现动态可调电磁吸收,但是这些结构复杂的吸收器具有较大的加工难度。

      自石墨烯被证明可在自然界稳定存在以来,因其独特的结构与性质吸引了研究者们的大量关注[13]。在石墨烯的众多特性中,其可通过电掺杂改变光学特性最为引人注目:石墨烯的等离子体频率可以通过改变石墨烯的化学势而改变[14],在动态吸收等方面具有独特的优势。例如,2013年A. Andryieuski等人利用单层石墨烯结合亚波长金属结构在太赫兹波段实现了窄带和宽带动态可调吸收[15],但是,该吸收器的调节形式仅局限于吸收率的调整。2014年,Y. Zhang等人将双层石墨烯网格与十字金属亚波长结构相结合,在太赫兹波段实现了极化无关的窄带吸收频率动态可调吸收[16],可惜的是该吸收器的频率相对调节范围仅有15%,并且未能实现吸收带宽的扩展。纵观近五年来国内外对基于石墨烯的太赫兹动态吸收研究,实现带宽的扩展和在保证吸收率的前提下实现吸收频率的调节仍然是主要研究目标[17-18]

      文中利用单层均匀石墨烯与硅材料球面-椭球面复合法布里-珀罗(FP)腔结构相结合,设计了一种新型的复合结构的石墨烯基太赫兹超材料吸收器。通过将石墨烯层与硅半球层的梯度变化组合,同时实现偏振无关和宽带高吸收的性能。该吸收器具有高度对称性,偏振角在90°范围内呈现出广角无关性;且入射角在入射角低于60°范围内可保持90%左右的宽带吸收性能。

    • 石墨烯是一种由碳原子组成的单原子层六方晶格二维材料,其电导率可由Kubo方程推导出的表面电导率表示。当其化学势μC >> kBT时,其表面电导率可表示为[19]

      $$ \begin{split} {\sigma _{2D}} \approx & \frac{{i{e^2}{k_B}T}}{{{\textit{π}} {\hbar ^2}\left( {\omega} + i2 {{\textit{π}} \varGamma } \right)}}\left[ {\frac{{{\mu _C}}}{{{k_B}T}} + 2ln\left( {{e^{\frac{{{\mu _C}}}{{{k_B}T}}}} + 1} \right)} \right] + \frac{{i{e^2}}}{{4{\textit{π}} \hbar }}\\ & {\rm ln}\left[ {\frac{{2\left| {{\mu _C}} \right| - \left( {\omega + i2\varGamma } \right)\hbar }}{{2\left| {{\mu _C}} \right| + \left( {\omega + i2\varGamma } \right)\hbar }}} \right] \end{split}$$ (1)

      式中:$ \hbar $kBe分别为约化的普朗克常数、玻耳兹曼常数以及电子电量;kBT表示热能;$\varGamma $是带电粒子散射率,其值设定为0.43 eV;$\mu _C $为石墨烯化学势,可由载流子浓度决定[20]

      $$ \begin{split} {n_0} = \frac{2}{{\pi {\hbar ^2}{v^2}}}\mathop \smallint \limits_0^{ + \infty } \varepsilon \left[ {{f_0}\left( {\varepsilon - {\mu _C}} \right) - {f_0}\left( {\varepsilon + {\mu _C}} \right)} \right]{\rm d}\varepsilon \end{split} $$ (2)

      式中:${f_0}=1+e^{\frac{{{(\varepsilon - \mu _C)}}}{{{k_B}T}}}$是费米-狄拉克分布函数,而费米速度v = 1.0×106 m/s。因为,载流子浓度可通过加载在石墨烯上的偏置电压控制。

      $$ \begin{split} {n_0} = \frac{{{\varepsilon _{dielectric}}{\varepsilon _0}{V_{bias}}}}{{ed}} \end{split}$$ (3)

      式中:d为介质层厚度;$\varepsilon _{dielectric} $为介质层相对介电常数。因此,可通过偏置电压实现对石墨烯化学势的有效调控。文中利用以下关系:

      $$ \begin{split} \sigma \left( \omega \right) = \frac{{{\sigma _{2D}}\left( \omega \right)}}{{{t_n}}} \end{split} $$ (4)

      式中:tn是石墨烯薄膜的膜厚,在文中的仿真计算中设定为1 nm。化学势具体调控手段为:如图1(b)所示,将吸收器石墨烯层外接电源作为正极,底层金属层作为负极,通过增大外加电压的方式增大石墨烯层的化学势。

      图  1  (a)宽带吸收器周期结构x-y平面示意图;(b)吸收器的三维单元示意图,其中P = 104 μm, R1= 11 μm, R2= 38 μm, R3= 26 μm, ts = 15 μm, tm = 2 μm

      Figure 1.  (a) x-y plane diagram of the periodic structure of the broadband absorber;(b) 3 D schematic diagram of the absorber, where P=104 μm, R1=11 μm, R2= 38 μm, R3=26 μm, ts =15 μm, tm =2 μm

    • 吸收器的结构示意图如图1所示,由下到上依次为连续金属接地层、PDMS(聚二甲基硅氧烷)薄膜、均匀单层石墨烯以及硅半球-半椭球复合结构层组成。该超材料结构单元的晶格周期P为104 μm,金属金作为接地层,其电导率为4.561×107 S/m,厚度tm为2 μm。其上为介质层,材料为低损耗的PDMS聚合物,相对介电常数ε2= 1.72,损耗角正切tanδ2= 0.15,介质厚度ts = 15 μm。介质层上为由半椭球与半圆球构成的硅周期结构,其单元结构分别由位于边缘的四个旋转对称半椭球和一个位于结构中心的半圆球组成。半椭球的长轴长R1 = 38 μm,短轴长R2 = 11 μm,且椭圆中心距离相邻侧壁距离分别为m = 40 μm和n = 12 μm,而中心半圆球半径R3= 26 μm。相对介电常数ε1和损耗角正切tanδ1分别为11.65和0.174。仿真基于有限元方法,利用商业计算软件CST Microwave Studio 2020中的频域求解器完成。仿真频率范围为0.1~10 THz,选择TE波作为入射波。其中,沿xy方向单胞采用周期性边界条件,则吸收器的吸收率可通过公式(5)计算:

      $$ A{\rm{ }} = 1{\rm{ }} - R-T$$ (5)

      式中:RT分别是反射率和透过率。这里,反射率R = | S11 | 2,其中S11是入射波的反射系数。由于底部金属反射层具有足够的厚度(ts = 2 μm),其远大于太赫兹波的趋肤深度而使得透过率T ≈ 0。因此,吸收率A仅与反射率有关,具体为:

      $$ A \approx {\rm{ }}1{\rm{ }} - {\rm{ }}{\left| {{S_{11}}} \right|^2} $$ (6)

      对于该结构上层的硅半椭球和半圆球,在TE波垂直入射时,波在半球内部沿x方向不断反射,该结构等效于几何轮廓尺寸不断随zy变化的法

      法布里-珀罗(FP)谐振腔,如图2(b)所示。所提出的吸收器中的激发石墨烯等离子体共振(GPRS)与具有波矢量的石墨烯等离子体的驻波模式的形成相关联。

      图  2  不同石墨烯化学势μ(0~0.9 eV)条件下吸收器吸收率随频率变化图

      Figure 2.  Absorption spectra of the proposed absorber with various values of the graphene chemical potential μ from 0 to 0.9 eV

      $$ \begin{split} {K_{GPRs}}\left( \omega \right) = \frac{{\pi {\hbar ^2}}}{{{e^2}{E_f}}}{\varepsilon _0}\left( {{\varepsilon _1} + {\varepsilon _2}} \right)\omega \left( {\omega + i{\tau ^{ - 1}}} \right) \end{split} $$ (7)

      式中:ε0是真空介电常数;ε1是硅半球的相对介电常数;ε2是PDMS的相对介电常数。半椭球和半圆球谐振腔对应的共振角频率ω可以用公式(8)来计算:

      $$ \begin{split} {\rm Re}({K_{GPRs}}\left( \omega \right)) = \frac{{\pi N + \varphi }}{{L\left( {y,z} \right)}} \end{split} $$ (8)

      式中:${\rm Re}({K_{GPRs}}\left( \omega \right))$${K_{GPRs}}\left( \omega \right) $的实部。在某个位置的半椭球F-P腔长L(y,z)可以近似为$L\left( {{y},{{z}}} \right)=$$2{R_2}\sqrt {1 - \dfrac{{{y^2} - {z^2}}}{{R_1^2}}} $,而半圆球的F-P腔长L(y,z)可以近似为$L\left( {{{y}},{{z}}} \right){\rm{ = }} 2{R_3}\sqrt {1 - \dfrac{{{y^2} - {z^2}}}{{R_3^2}}}$ N是正整数,是谐振模式的阶数;φ是FP腔中的反射系数的相位,文中研究可以忽略不计。因此,半椭球和半圆球第N阶模式的ωN1ωN2可以分别简化为:

      $$ \begin{split} {\omega _{N1}} = \frac{e}{\hbar }\sqrt {\frac{{{E_f}\left( {N\pi + \varphi } \right){R_2}}}{{2\pi {\varepsilon _0}\left( {{\varepsilon _1} + {\varepsilon _2}} \right)\sqrt {{R_2}^2{R_1}^2 - {R_1}^2{y^2} - {R_2}^2{z^2}} }}} \end{split} $$ (9)
      $$ \begin{split} {\omega _{N2}} = \frac{e}{\hbar }\sqrt {\frac{{{E_f}\left( {N\pi + \varphi } \right){R_3}}}{{2\pi {\varepsilon _0}\left( {{\varepsilon _1} + {\varepsilon _2}} \right)\sqrt {{R_3}^4 - {R_3}^2({y^2} - {z^2})} }}} \end{split} $$ (10)

      不同的整数N对应多个不同的离散GPR,然而,文中所示吸收器,半椭球和半圆球的几何轮廓不断变化,因此,通过连续改变L,相对于YZ形成一系列FP谐振腔;因此,将导出具有连续频率的谐振腔等离子体共振(FPRS)。连续FPRS不仅填充离散GPRS的频率间隔,还不会形成重叠的谐振频率。因此,该宽带吸收器的吸收谱较平滑[14]

    • 前期研究表明,石墨烯的表面电导率与其THz波段的化学势(费米能级)密切相关。而且,石墨烯的费米能级可以通过施加偏置电压来调节。因此,可将石墨烯材料优越的电磁特性应用于吸波器的振幅可调特性。为了验证利用偏置电压的吸收器的可调吸收性能,课题组模拟了TE极化法向入射下的吸收谱与频率和化学势的函数关系。

      图2所示,在室温下(T = 300 K),当调整石墨烯片的化学势时,吸收器的吸收幅值可发生显著变化。当化学势在0.0~0.9 eV之间变化时,其0.1 THz处吸收率可获得10~92%宽范围调节;而在1.9 THz处亦可获得约50~100%的吸收率变化。而且,当化学势从0 eV增加到0.2 eV时,约1~10 THz吸收曲线大体上在逐渐变得平坦,其平均吸收率获得显著得到提升;当化学势超过0.2 eV后,吸收光谱除原位于约0.82 THz共振峰发生明显蓝移外,其平均吸收率随化学势的增加变化不明显。此外,当化学势为0.2 eV时,吸收器可分别在0.1~3.1 THz和3.7~6.3 THz(总带宽约5.7 THz)的频率范围内获得高于90%的吸收性能,其最高吸收率可达到99.2%;而且,在0~10 THz的整个太赫兹波段均可获得超过70%的吸收率。类似地,当化学势调整为0.9 eV时,可分别在1.7~6.7 THz、6.9~7.9 THz,以及8.9~9.9 THz(总带宽约7 THz)三个频率范围获得超过90%的吸收率,最高吸收率可达到99.8%。

      为了更加清晰了解吸收光谱随石墨烯化学势变化的作用机理,计算了石墨烯化学势μ分别为0.1、0.2、0.4 eV时,在TE垂直入射情况下,在石墨烯层上方50 nm处的电场分布,结果如图3所示。

      图  3  石墨烯层上方z = 50 nm处吸收器在不同石墨烯化学势μ下电场分布图:(a)μ = 0.1 eV;(b) μ = 0.2 eV;(c) μ = 0.4 eV

      Figure 3.  The electric field distributions (|E|) of the proposed absorber on the xy plane with z = 50 nm at the different graphene chemical potential μ of (a) 0.1, (b) 0.2 and (c) 0.4 eV

      图3可知,石墨烯化学势不同时,吸收器谐振模式存在明显不同。当化学势为0.1 eV时,电场主要被强制约束在半椭球/半球体之间的间隙中,主要由相邻半椭球/半球体之间的耦合作用造成太赫兹波的捕获与吸收。

      当化学势增加至0.2 eV时,电场峰值由间隙移动至半椭球/半球体边缘处,说明其共振模式由相邻半椭球/半球体之间的耦合转变为半椭球或半球体自身的基频模式,且当化学势继续增加时,电场主要分布于半球体上,其太赫兹吸收主要来源于半球体的电偶极子共振。这说明该结构的优异吸收特性与其结构密切相关。由于石墨烯的存在,可以激发多个离散的离子体共振(GPRs)。再结合半椭球/半球体结构所具有的连续改变的几何轮廓,形成一系列法布里-珀罗谐振腔(F-P腔),因此可导出具有连续频率的法布里珀罗共振(FPRs)[21]。而且,连续的FPRs不仅会填充离散GPRs的频率间隔,而且还不会形成重叠的谐振频率。即,吸收器的宽带吸收光谱是F-P腔长度随yz的连续变化的结果,不同谐振频率的重叠可以平滑和扩展整个吸收光谱,从而获得一种具有平滑吸收光谱的宽带吸收器。如将该结构与介质模式的椭圆柱体或单纯的半球体结构相比,半椭球体的存在使其在yz方向上具有更温和的梯度分布,因此有可能使其存在更复杂、更紧凑的FPRs,从而获得更宽、更平滑的吸收曲线。

      为了探讨石墨烯以及亚波长F-P谐振腔的作用,分别计算了无石墨烯无谐振结构(NN)、有石墨烯且只有半球结构(G-SemiSphere)、有石墨烯且只有半椭球结构(G-SemiEllip)、无石墨烯但具有半椭球/半球谐振结构(N- SemiEllip SemiSphere)以及具有石墨烯与半椭球/半球谐振结构(G- SemiEllip/SemiSphere)的吸收率随频率变化关系,结果如图4所示。

      图  4  μ = 0.2 eV时五种不同结构的吸收谱:(1)无石墨烯无谐振结构;(2)有石墨烯且只有半球结构(蓝线);(3)有石墨烯且只有半椭球结构(绿线);(4)无石墨烯但具有半椭球/半球谐振结构(粉线);(5)有石墨烯与半椭球/半球谐振结构(红线)

      Figure 4.  (a) Simulated absorption spectra of the proposed G-SemiEllip/SemiSphere absorber (red line) and four other absorbers with NN (black line), G-SemiSphere (blue line), G-SemiEllip (green line) and N-SemiEllip/SemiSphere (pink line)

      总体而言,相对于NN结构(在0~10 THz频率范围内,其平均吸收率约为40.1%),表面F-P谐振结构与石墨烯的加入可显著提升其吸收性能。如吸收器均具有石墨烯层,随着在吸收器表面增加半圆球、半椭球以及半椭球/半圆球复合结构,吸收器的平均吸收率逐步递增,其分别为59.8%、84.6%与89.3%。

      此外,通过比较有、无石墨烯层的计算结果可知,石墨烯层的存在对吸收性能有显著影响。当石墨烯不存在时,特别是在5 THz以下的频率区域,整体吸收光谱性能显著恶化。在吸收器中增加均匀石墨烯层,可显著改善0~5 THz范围的吸收性能,从而将吸收器平均吸收率从76.2%提高到89.3%。上述结果表明,石墨烯层与表面半椭球/半圆球复合F-P谐振腔结构对于改善吸收器性能是非常有效的。

      此外,课题组通过分析吸收器在不同共振频率下的电场分布以了解其吸收器作用机理。图5为以0.82 THz与4.95 THz共振频率为例,分别取石墨烯化学势为0.2 eV,且在石墨烯层上方50 nm处的电场强度分布。计算均假设入射波为垂直入射的TE极化波。

      图  5  在 (a) 0.82 THz和(b) 4.95 THz的不同频率下,在z =50 nm的x-y平面上吸收器的电场分布 (|E|)

      Figure 5.  The electric field distributions (|E|) of the proposed absorber on the x-y plane with z = 50 nm at the different frequencies of (a) 0.82 THz, and (b) 4.95 THz

      由图可知,在这两个频率下,吸收体分别在半椭球/半球体边缘或其间隙形成强场约束,从而产生高效的太赫兹捕获和吸收。当共振频率为0.82 THz时,电场被限制在半椭球与半球体的边缘;而在共振频率4.95 THz时,电场则被强制约束在半椭球/半球体之间的间隙中。这说明在0.82 THz共振频率下的吸收主要是由于单独的半椭球或半球体的基频共振模式(电偶极子共振);而在4.95 THz共振频率下,则主要是由于相邻半椭球/半球体之间的耦合作用。再结合图1(a)可知,不同共振区可以产生重叠以实现宽带吸收。而且,据参考文献[14]报道,半椭球F-P谐振腔结构在在较高的谐振频率下,由于不同的谐振模式的严重重叠,即形成混合模式谐振,模式特性将变得更加复杂,可使吸收曲线变得更为平坦,从而有利于实现宽带高效吸收。

      进一步,笔者课题组研究了太赫兹吸收行为与极化和入射角的关系。为了简化分析,将化学势取为0.2 eV,并固定除极化角或入射角以外的其他所有吸收器参数。图6 (a)给出了垂直入射波偏振角在0°~90°之间变化时该吸收体的吸收光谱。吸收率甚至随入射角增大而获得一定提高,这可能是由于随着入射角θ的增大,结构某些部分的寄生共振增加,从而产生一定的广角吸收特性,该吸收特性应该主要与产生亚波长结构的谐振有关,而与入射角的依赖性较小。对于TM偏振,即使在70~80°的角度下,最大吸收率仍然很大。这是因为入射波的磁场方向可随入射角的变化而变化,这对于有效驱动表面亚波长谐振结构、石墨烯,以及底部金属层之间的循环电流是有利的,对于保持阻抗匹配也是非常重要的。上述结果表明,该吸收体的半椭球/半球体亚波长结构可以补偿石墨烯超材料表面与自由空间太赫兹波之间的动量不匹配,从而表现出优良的广角特性。

      图  6  (a) TE波斜入射时吸收率随极化角ϕ(电场方向与y轴的夹角)的变化率; (b)~(c) TE/TM波斜入射时吸收率随入射角θ(入射方向与z轴的夹角)的变化率

      Figure 6.  (a) The change rate of absorptivity with polarization angle (angle between electric field direction and y-axis) at oblique incident of TE wave; (b)-(c) The rate at which the absorptivity varies with the incidence angle θ (the angle between the incident direction and the z-axis) for oblique incidence of TE/TM waves

      最后,在制造方面,首先300 ℃下通过常规的真空镀膜工艺在石英基板表面镀上1~2 μm厚度的金属反射镜;PDMS溶液的预聚体与固化剂按照10 : 1的比例混合,用60 wt%甲苯稀释,将制备好的PDMS溶液放于Si基板上在100 ℃下用热板固烘烤两小时,得到12~18 μm厚度的PDMS薄膜;运用化学气相沉积法 (CVD)在铜箔衬底上生长2~5层石墨烯, 然后将石墨烯薄膜用湿化学腐蚀基底法从铜箔衬底上转移至PDMS表面;最后,采用光刻的方法在石墨烯表面刻蚀出硅半球层。由于许多应用需要廉价,大面积的等离子体吸收器,因此制造方法是至关重要的。课题组使用步进和闪光紫外线印记光刻制造等离子体吸收器,其和连续制造方法相比具有可重复使用的掩模和大面积制造的优点。

    • 文中设计了一种基于石墨烯的性能可调、宽带、极化与入射角不敏感的太赫兹超材料吸收器。该吸收器的吸收单元由表面半椭球/半球体亚波长谐振结构、石墨烯、介质层以及金属背板等组成。仿真工作由CST Studio Suite 2020 这款光学仿真工具来完成。

      通过石墨烯化学势的调控,可在其化学势分别为0.2 eV与0.9 eV时获得5.7 THz与7 THz左右的超宽带吸收,其吸收率均可高于90%,并可在一定谐振频率(如0.82、4.95 THz等)实现近乎完美的太赫兹波吸收。器件形成宽带吸收的机理源于由石墨烯激发的等离子体共振,以及源于表面亚波长结构所提供的连续可变的多模法布里-珀罗共振的协同作用。此外,该吸收器同时具备入射角度60°以下的高吸收特性和360°极化不敏感特性,最后,介绍了该吸收器的具体制备工艺。

参考文献 (21)

目录

    /

    返回文章
    返回