基于狄拉克半金属纳米线的太赫兹可调七波段完美吸收器的模拟仿真

上一篇

下一篇

卢文强, 易颖婷, 宋前举, 周自刚, 易有根, 曾庆栋, 易早. 基于狄拉克半金属纳米线的太赫兹可调七波段完美吸收器的模拟仿真[J]. 物理学报, 2025, 74(3): 034101-1. doi: 10.7498/aps.74.20241516
引用本文: 卢文强, 易颖婷, 宋前举, 周自刚, 易有根, 曾庆栋, 易早. 基于狄拉克半金属纳米线的太赫兹可调七波段完美吸收器的模拟仿真[J]. 物理学报, 2025, 74(3): 034101-1. doi: 10.7498/aps.74.20241516
Wenqiang LU, Yingting YI, Qianju SONG, Zigang ZHOU, Yougen YI, Qingdong ZENG, Zao YI. Simulation of terahertz tunable seven-band perfect absorber based on high frequency detection function of Dirac semi-metallic nanowires[J]. Acta Physica Sinica, 2025, 74(3): 034101-1. doi: 10.7498/aps.74.20241516
Citation: Wenqiang LU, Yingting YI, Qianju SONG, Zigang ZHOU, Yougen YI, Qingdong ZENG, Zao YI. Simulation of terahertz tunable seven-band perfect absorber based on high frequency detection function of Dirac semi-metallic nanowires[J]. Acta Physica Sinica, 2025, 74(3): 034101-1. doi: 10.7498/aps.74.20241516

基于狄拉克半金属纳米线的太赫兹可调七波段完美吸收器的模拟仿真

    作者简介: 卢文强: 2434691674@qq.com .
    通讯作者: E-mail: qjsong@swust.edu.cn).;  E-mail: yizaomy@swust.edu.cn.
  • 中图分类号: 41.20.Jb, 78.20.Ci, 42.70.Qs

Simulation of terahertz tunable seven-band perfect absorber based on high frequency detection function of Dirac semi-metallic nanowires

    Corresponding authors: E-mail: qjsong@swust.edu.cn.;  E-mail: yizaomy@swust.edu.cn.
  • MSC: 41.20.Jb, 78.20.Ci, 42.70.Qs

  • 摘要: 设计了一种高灵敏度、高品质因子、高品质因数、高频探测、双固定功能的太赫兹可调完美吸收器. 该吸收器可实现4—14.5 THz范围内7个波段的完美吸收. 在进行结构设计时将线阵结构的参数与周期进行了关联. 通过计算吸收器的相对阻抗来对器件宏观层面的电磁进行解释, 并通过分析共振频率点的表面电场和磁场分布, 来分析该器件的物理机制. 计算了7个共振频点的品质因子Q, 其中最大Q值为219.41. 通过改变外部折射率, 该吸收器的灵敏度和品质因数值最大可达5421.43 GHz/RIU和35.204 RIU–1. 通过讨论关键参数对器件的影响, 得出该器件可实现双固定性能的选择、七波段吸收以及全波段反射. 通过改变狄拉克半金属的费米能级, 证明该吸收器具有良好的动态调节能力. 通过改变外部电磁波的入射角发现该器件在中低频段具有良好的稳定性, 但在高频段受外部入射角影响较大. 本文所提出的吸收器在成像、探测、检测等领域具有巨大的应用潜力, 相关工作对光电器件的设计提供了思路.
  • 加载中
  • 图 1  吸收器的单元结构组成及其结构参数

    Figure 1.  Unit structure of the absorber and its structural parameters.

    图 2  不同费米能级下BDS介电常数的实部(a)和虚部(b)随频率的变化

    Figure 2.  Variations of real part (a) and imaginary part (b) of the permittivity of BDS with frequency at different Fermi levels.

    图 3  (a)吸收器在4—14.5 THz范围内的特性曲线; (b)吸收器在其工作区间的相对阻抗(实部和虚部)示意图

    Figure 3.  (a) Characteristic curves of the absorber in the range of 4–14.5 THz; (b) diagram of the relative impedance (real and imaginary parts) of the absorber in its operating interval.

    图 4  吸收器在不同频率处的电场分布 (a) f 1 = 5.032 THz; (b) f 2 = 5.859 THz; (c) f 3 = 7.674 THz; (d) f 4 = 9.654 THz; (e) f 5 = 11.656 THz; (f) f 6 = 12.514 THz; (g) f 7 = 14.01 THz

    Figure 4.  Electric field distribution of absorber at different frequencies: (a) f 1 = 5.032 THz; (b) f 2 = 5.859 THz; (c) f 3 = 7.674 THz; (d) f 4 = 9.654 THz; (e) f 5 = 11.656 THz; (f) f 6 = 12.514 THz; (g) f 7 =14.01 THz.

    图 5  吸收器在不同频率处的磁场分布 (a) f 1 = 5.032 THz; (b) f 2 = 5.859 THz; (c) f 3 = 7.674 THz; (d) f 4 = 9.654 THz; (e) f 5 = 11.656 THz; (f) f 6 = 12.514 THz; (g) f 7 = 14.01 THz

    Figure 5.  Magnetic field distribution of absorber at different frequencies: (a) f 1 = 5.032 THz; (b) f 2 = 5.859 THz; (c) f 3 = 7.674 THz; (d) f 4 = 9.654 THz; (e) f 5 = 11.656 THz; (f) f 6 = 12.514 THz; (g) f 7 =14.01 THz.

    图 6  (a)不同折射率下吸收器的吸收光谱; (b)谐振频率点随折射率的变化; (c) 7种模式的吸收率与折射率的对应关系

    Figure 6.  (a) Absorption spectra of absorbers with different refractive indices; (b) the change of resonant frequency points with refractive index; (c) the corresponding relationship between absorptivity and refractive index of 7 modes.

    图 7  (a)不同BDS的费米能对吸收率的影响; (b) 7个共振频率点与费米能的关系; (c) 7个共振频率点处吸收率与费米能的关系

    Figure 7.  (a) Effect of Fermi energy of different Dirac semi-metals on absorption efficiency; (b) the relationship between seven resonance frequency points and Fermi energy; (c) the relationship between the absorption rate and Fermi energy at seven resonance frequency points.

    表 1  本文所提出的吸收器的Q值与近年来类似吸收器之间的比较

    Table 1.  Comparison of the Q value of the proposed absorber with similar absorbers in recent years.

    参考文献[48][49][50][51][52]本文
    Q55.597377.89106154219.41
    下载: 导出CSV

    表 2  本文所提出的吸收器的灵敏度与近年来类似吸收器之间的比较

    Table 2.  Comparison of the sensitivity of the proposed absorber with similar absorbers in recent years.

    参考文献 [51] [54] [55] [56] [57] 本文
    S/(GHZ·RIU–1) 23.8 74.43 96.2 560 2475 5421.43
    下载: 导出CSV
  • [1] Li J T, Wang G C, Yue Z, Liu J Y, Li J, Zheng C L, Zhang Y T, Zhang Y, Yao J Q 2022 Opto-Electron Adv. 5 210062 doi: 10.29026/oea.2022.210062
    [2] 张学进, 陆延青, 陈延峰, 朱永元, 祝世宁 2017 物理学报 66 148705 doi: 10.7498/aps.66.148705 Zhang X J, Lu Y Q, Chen Y F, Zhu Y Y, Zhu S N 2017 Acta Phys. Sin. 66 148705 doi: 10.7498/aps.66.148705
    [3] 黄若彤, 李九生 2023 物理学报 72 054203 doi: 10.7498/aps.72.20221962 Huang R T, Li J S 2023 Acta Phys. Sin. 72 054203 doi: 10.7498/aps.72.20221962
    [4] Yue Z, Li J T, Li J, Zheng C L, Liu J Y, Wang G C, Xu H, Chen M Y, Zhang Y T, Zhang Y, Yao J Q 2022 Opto-Electron Sci. 1 210014 doi: 10.29026/oes.2022.210014
    [5] Li W X, Zhao W C, Cheng S B, Yang W X, Yi Z, Li G F, Zeng L C, Li H L, Wu P H, Cai S S 2023 Surf. Interfaces 40 103042 doi: 10.1016/j.surfin.2023.103042
    [6] Sun W F, Wang X K, Zhang Y 2022 Opto-Electron. Sci. 1 220003 doi: 10.29026/oes.2022.220003
    [7] Li W X, Yi Y T, Yang H, Cheng S B, Yang W X, Zhang H F, Yi Z, Yi Y G, Li H L 2023 Commun. Theor. Phys. 75 045503 doi: 10.1088/1572-9494/acbe2d
    [8] 陈俊, 杨茂生, 李亚迪, 程登科, 郭耿亮, 蒋林, 张海婷, 宋效先, 叶云霞, 任云鹏, 任旭东, 张雅婷, 姚建铨 2019 物理学报 68 247802 doi: 10.7498/aps.68.20191216 Chen J, Yang M S, Li Y D, Cheng D K, Guo G L, Jiang L, Zhang H T, Song X X, Ye Y X, Ren Y P, Ren X D, Zhang Y T, Yao J Q 2019 Acta Phys. Sin. 68 247802 doi: 10.7498/aps.68.20191216
    [9] Zhao H, Wang X K, Liu S T, Zhang Y 2023 Opto-Electron Adv. 6 220012 doi: 10.29026/oea.2023.220012
    [10] Song H J, Nagatsuma T 2011 IEEE T. Thz. Sci. Techn. 1 256 doi: 10.1109/TTHZ.2011.2159552
    [11] Gigli C, Leo G 2022 Opto-Electron Adv. 5 210093 doi: 10.29026/oea.2022.210093
    [12] Li F Y, Li Y X, Tang T T, Liao Y L, Lu Y C, Liu X Y, Wen Q Y 2022 J. Alloys Compd. 928 167232 doi: 10.1016/j.jallcom.2022.167232
    [13] Zhu J, Xiong J Y 2023 Measurement 220 113302 doi: 10.1016/j.measurement.2023.113302
    [14] Cheng Y Z, Withayachumnankul W, Upadhyay A, Headland D, Nie Y, Gong R Z, Bhaskaran M, Sriram S, Abbott D 2015 Adv. Opt. Mater. 3 376 doi: 10.1002/adom.201400368
    [15] Li W X, Cheng S B, Zhang H F, Yi Z, Tang B, Ma C, Wu P H, Zeng Q D, Raza R 2024 Commun. Theor. Phys. 76 065701. doi: 10.1088/1572-9494/ad3b8f
    [16] Chen H T, Padilla W J, Zide J M O, Gossard A C, Taylor A J, Averitt R D 2006 Nature 444 597 doi: 10.1038/nature05343
    [17] Cao T, Lian M, Chen X Y, Mao L B, Liu K, Jia J Y, Su Y, Ren H N, Zhang S J, Xu Y H, Chen J J, Tian Z, Guo D M 2022 Opto-Electron Sci. 1 210010 doi: 10.29026/oes.2022.210010
    [18] 沈晓鹏, 崔铁军, 叶建祥 2012 物理学报 61 058101 doi: 10.7498/aps.61.058101 Shen X P, Cui T J, Ye J X 2012 Acta Phys. Sin. 61 058101 doi: 10.7498/aps.61.058101
    [19] He M Y, Wang Q Q, Zhang H, Xiong J, Liu X P, Wang J Q 2024 Phys. Scr. 99 035506 doi: 10.1088/1402-4896/ad203e
    [20] Landy N I, Sajuyigbe S, Mock J J, Smith D R, Padilla W J 2008 Phys. Rev. Lett. 100 207402 doi: 10.1103/PhysRevLett.100.207402
    [21] Zhan Y, Yin H Y, Wang J H, Yao H W, Fan C Z 2022 Res. Opt. 8 100255 doi: 10.1016/j.rio.2022.100255
    [22] Mao Y, Zhang H, Xiong J, Liu X P, Wang Q Q, Wang J Q 2024 J. Phys. D: Appl. Phys. 57 255111 doi: 10.1088/1361-6463/ad3376
    [23] Ullah K, Meng Y F, Sun Y, Yang Y K, Wang X J, Wang A R, Wang X R, Xiu F X, Shi Y, Wang F Q 2020 Appl. Phys. Lett. 117 011102 doi: 10.1063/5.0010707
    [24] Moll P J W, Nair N L, Helm T, Potter A C, Kimchi I, Vishwanath A, Analytis J G 2016 Nature 535 266 doi: 10.1038/nature18276
    [25] Borisenko S, Gibson Q, Evtushinsky D, Zabolotnyy V, Büchner B, Cava R J 2014 Phys. Rev. Lett. 113 027603 doi: 10.1103/PhysRevLett.113.027603
    [26] Cheng S B, Li W X, Zhang H F, Akhtar M N, Yi Z, Zeng Q D, Ma C, Sun T Y, Wu P H, Ahmad S 2024 Opt. Commun. 569 130816 doi: 10.1016/j.optcom.2024.130816
    [27] Chen Z Y, Cheng S B, Zhang H F, Yi Z, Tang B, Chen J, Zhang J G, Tang C J 2024 Phys. Lett. A 517 129675 doi: 10.1016/j.physleta.2024.129675
    [28] Li W X, Ma J, Zhang H F, Cheng S B, Yang W X, Yi Z, Yang H, Zhang J G, Wu X W, Wu P H 2023 Phys. Chem. Chem. Phys. 25 8489 doi: 10.1039/D2CP05562G
    [29] Liang S R, Cheng S B, Zhang H F, Yang W X, Yi Z, Zeng Q D, Tang B, Wu P H, Ahmad S, Sun T Y 2024 Ceram. Int. 50 23611 doi: 10.1016/j.ceramint.2024.04.085
    [30] Li W X, Zhao W C, Cheng S B, Zhang H F, Yi Z, Sun T Y, Wu P H, Zeng Q D, Raza R 2024 Opt. Lasers Eng. 181 108368 doi: 10.1016/j.optlaseng.2024.108368
    [31] Ma J, Wu P H, Li W X, Liang S R, Shangguan Q Y, Cheng S B, Tian Y H, Fu J Q, Zhang L B 2023 Diam. Relat. Mater. 136 109960 doi: 10.1016/j.diamond.2023.109960
    [32] Zhu J, Xiong J Y 2023 Opt. Express 31 36677 doi: 10.1364/OE.501242
    [33] Huang S L, Chen Y, Yu C C, Chen S J, Zhou Z K, Liang J, Dai W 2024 Chinese J. Phys. 89 740 doi: 10.1016/j.cjph.2023.09.005
    [34] Wang J L, Hassan M, Liu J W, Yu S H 2018 Adv. Mater. 30 1803430 doi: 10.1002/adma.201803430
    [35] Cheng C, Gonela R K, Gu Q, Haynie D T 2005 Nano Lett. 5 175 doi: 10.1021/nl048240q
    [36] Qu T, Zhao Y B, Li Z B, Wang P P, Cao S B, Xu Y W, Li Y Y, Chen A H 2016 Nanoscale 8 3268 doi: 10.1039/C5NR09140C
    [37] Fu R, Chen K X, Li Z L, Yu S H, Zheng G X 2022 Opto-Electron Sci. 1 220011 doi: 10.29026/oes.2022.220011
    [38] Shangguan Q Y, Zhao Y, Song Z J, Wang J, Yang H, Chen J, Liu C, Cheng S B, Yang W X, Yi Z 2022 Diam. Relat. Mater. 128 109273 doi: 10.1016/j.diamond.2022.109273
    [39] Li W X, Liu Y H, Ling L, Sheng Z X, Cheng S B, Yi Z, Wu P H, Zeng Q D, Tang B, Ahmad S 2024 Surf. Interfaces 48 104248 doi: 10.1016/j.surfin.2024.104248
    [40] Liang S R, Xu F, Li W X, Yang W X, Cheng S B, Yang H, Chen J, Yi Z, Jiang P P 2023 Appl. Therm. Eng. 232 121074 doi: 10.1016/j.applthermaleng.2023.121074
    [41] Zhang Y J, Yi Y T, Li W X, Liang S R, Ma J, Cheng S B, Yang W X, Yi Y G 2023 Coatings 13 531 doi: 10.3390/coatings13030531
    [42] Huang X M, Chen Y, Chen S J, Yang K, Liang J, Zhou Z K, Dai W 2023 Res. Phys. 47 106364 doi: 10.1016/j.rinp.2023.106364
    [43] Zeng C, Lu H, Mao D, Du Y Q, Hua H, Zhao W, Zhao J L 2022 Opto-Electron. Adv. 5 200098 doi: 10.29026/oea.2022.200098
    [44] Li W X, Liu M S, Cheng S B, Zhang H F, Yang W X, Yi Z, Zeng Q D, Tang B, Ahmad S, Sun T Y 2024 Diam. Relat. Mater. 142 110793 doi: 10.1016/j.diamond.2024.110793
    [45] Liang S R, Xu F, Yang H, Cheng S B, Yang W X, Yi Z, Song Q J, Wu P H, Chen J, Tang C J 2023 Opt. Laser Technol. 158 108789 doi: 10.1016/j.optlastec.2022.108789
    [46] Shangguan Q Y, Chen Z, Yang H, Cheng S B, Yang W X, Yi Z, Wu X, Wang S, Yi Y, Wu P H 2022 Sensors 22 6483 doi: 10.3390/s22176483
    [47] Shangguan Q Y, Chen H, Yang H, Liang S R, Zhang Y J, Cheng S B, Yang W X, Yi Z, Luo Y, Wu P H 2022 Diam. Relat. Mater. 125 108973 doi: 10.1016/j.diamond.2022.108973
    [48] Sourav A, Li Z, Huang Z, Botcha V D, Hu C, Ao J P, Peng Y, Kuo H C, Wu J, Liu X, Ang K W 2018 Adv. Opt. Mater. 6 1800461 doi: 10.1002/adom.201800461
    [49] Piper J R, Liu V, Fan S 2014 Appl. Phys. Lett. 104 251110 doi: 10.1063/1.4885517
    [50] Li H J, Qin M, Wang L L, Zhai X, Ren R Z, Hu J G 2017 Opt. Express 25 31612 doi: 10.1364/OE.25.031612
    [51] Wang Y L, Cheng W, Qin J Y, Han Z H 2019 Opt. Commun. 434 163 doi: 10.1016/j.optcom.2018.10.061
    [52] Zhou K, Cheng Q, Lu L, Li B W, Song J L, Luo Z X 2020 Opt. Express 28 1647 doi: 10.1364/OE.383988
    [53] Li W X, Xu F, Cheng S B, Yang W X, Liu B, Liu M S, Yi Z, Tang B, Chen J, Sun T Y 2024 Opt. Laser Technol. 169 110186 doi: 10.1016/j.optlastec.2023.110186
    [54] Zheng W, Fan F, Chen M, Bai J J, Chang S J 2017 Infrared Laser Eng. 46 420003 doi: 10.3788/IRLA201746.0420003
    [55] Pan W, Yan Y J, Ma Y, Shen D J 2019 Opt. Commun. 431 115 doi: 10.1016/j.optcom.2018.09.014
    [56] Lu W Q, Wu P H, Bian L, Yan J Q, Yi Z, Liu M S, Tang B, Li G F, Liu C 2024 Opt. Laser Technol. 174 110650 doi: 10.1016/j.optlastec.2024.110650
    [57] Chen T, Jiang W J, Yin X H 2021 Micro. Nanostructures 154 106898 doi: 10.1016/j.spmi.2021.106898
  • 加载中
图( 7) 表( 2)
计量
  • 文章访问数:  96
  • HTML全文浏览数:  96
  • PDF下载数:  2
  • 施引文献:  0
出版历程
  • 收稿日期:  2024-10-29
  • 刊出日期:  2025-02-05

基于狄拉克半金属纳米线的太赫兹可调七波段完美吸收器的模拟仿真

    通讯作者: E-mail: qjsong@swust.edu.cn).; 
    通讯作者: E-mail: yizaomy@swust.edu.cn.
    作者简介: 卢文强: 2434691674@qq.com
  • 1. 西南科技大学数理学院, 绵阳 621010
  • 2. 中南大学物理学院, 长沙 410083
  • 3. 湖北工程学院物理与电子信息工程学院, 孝感 432000
  • 4. 吉首大学化学化工学院, 吉首 416000

摘要: 设计了一种高灵敏度、高品质因子、高品质因数、高频探测、双固定功能的太赫兹可调完美吸收器. 该吸收器可实现4—14.5 THz范围内7个波段的完美吸收. 在进行结构设计时将线阵结构的参数与周期进行了关联. 通过计算吸收器的相对阻抗来对器件宏观层面的电磁进行解释, 并通过分析共振频率点的表面电场和磁场分布, 来分析该器件的物理机制. 计算了7个共振频点的品质因子Q, 其中最大Q值为219.41. 通过改变外部折射率, 该吸收器的灵敏度和品质因数值最大可达5421.43 GHz/RIU和35.204 RIU–1. 通过讨论关键参数对器件的影响, 得出该器件可实现双固定性能的选择、七波段吸收以及全波段反射. 通过改变狄拉克半金属的费米能级, 证明该吸收器具有良好的动态调节能力. 通过改变外部电磁波的入射角发现该器件在中低频段具有良好的稳定性, 但在高频段受外部入射角影响较大. 本文所提出的吸收器在成像、探测、检测等领域具有巨大的应用潜力, 相关工作对光电器件的设计提供了思路.

English Abstract

    • 太赫兹波是介于微波与红外波之间, 频率范围是在0.1—10 THz波段的电磁波 [13]. 迄今为止, 相对于其他波段, 研究者们对于太赫兹波段的研究较为有限, 对应太赫兹器件研究也不如其他波段 [46]. 所以, 对于发展太赫兹波段的器件研究依旧是时代所需. 近年来, 诸多太赫兹器件应运而生, 包括太赫兹探测器、太赫兹传输器以及太赫兹吸收器等[79]. 这些太赫兹器件在通信、探测、医学成像、国家安全等领域都具有非常大的应用潜力[1012]. 但是, 传统材料在应用到太赫兹器件的时候, 其劣势非常明显. 尽管有少部分材料可以通过构筑谐振腔体来实现对太赫兹波的吸收, 但其作用效果却总不明显[13,14]. 太赫兹吸收器的工作原理基于不同材料与太赫兹辐射的相互作用. 常见的吸收机制包括: 电导吸收、带隙吸收、界面吸收与共振吸收. 对于共振吸收来说, 在特殊材料或结构中, 当太赫兹波的频率与材料或器件的固有频率相匹配时, 会发生共振吸收. 这种共振吸收通常出现在具有周期性结构(如表面等离激元、超表面)的材料中. 而传统金属材料往往面临不可调节性的问题, 因此所设计的太赫兹器件发挥的作用单一, 没有多种用途, 这极大地限制了太赫兹器件的发展.

      作为一种人工自主设计结构的材料, 超材料为解决上述问题提供了一种有效且可行的方案. 自2000年超材料被成功报道出来后, 目前, 超材料已被广泛地应用在电磁学领域 [1517]. 由于自身结构的特殊性, 超材料具备自然材料所不具备的诸多特性, 例如负折射率、可动态调控、完美吸收等 [18,19]. 2008年, Landy等[20]设计出第一个太赫兹波段的超材料吸收器, 其表现出良好的吸收性能. 目前, 采用超材料来构筑的太赫兹波吸收器成为大家的研究热点[21,22]. 其中, 由于其特殊的传感潜力, 窄带超材料吸收器得到了广泛研究.

      狄拉克半金属(BDS)的费米能级具有线性色散关系, 类似于狄拉克费米子的色散关系[23,24]. BDS的介电常数可以通过改变费米能级来控制, 而费米能级可以通过化学掺杂来实现. 此外, 由于晶体对称性保护作用, BDS具有很高的载流子迁移率. Borisenko等[25]在他们的Cd3As2实验中展示了BDS的优越性能. 基于BDS的吸收器器件, 具有优越的传感性能[2628].

      本文提出了一种基于BDS的七波段完美吸收器. 所设计的结构是BDS-介质-金属的三层结构. 该吸收器可实现4—14.5 THz范围内f1 = 5.032 THz(84.43%), f2 = 5.859 THz(96.23%), f3 = 7.674 THz(91.36%), f4 = 9.654 THz(99.02%), f5 = 11.656 THz(93.84%), f6 = 12.514 THz(98.47%), f7 = 14.01 THz(97.32%) 7个波段的完美吸收. 通过阻抗匹配理论证实了这一结果, 并通过分析7个频点的电场和磁场得到吸收器的物理机制. 探究了器件吸收峰的品质因子和灵敏度, 得到每个吸收峰的品质因数. 讨论了关键的结构参数对器件吸收率的影响. 通过外加偏置电压以及化学掺杂等手段改变BDS的费米能级, 实现吸收器的动态调节. 最后, 研究了外部入射的电磁波的角度对吸收率的影响, 发现中低频段的吸收可用稳定到50°之前, 但在高频段的吸收只能稳定到20°.

    • 本文所提出完美吸收器的结构及具体参数如图1所示. 从图1可知, 该器件有3层结构. 它的底部由金(Au)薄膜组成, 其厚度H1 = 0.5 μm. 中间采用了无损的二氧化硅(SiO2)作为介质层, 其介电常数$ {\varepsilon }_{{{\mathrm{S}}{\mathrm{i}}{\mathrm{O}}}_{2}} $ = 2.13, 厚度H2 = 17.0 μm. 顶部采用的是BDS材料, 我们将其设计成线阵结构. 纳米线BDS的半径为R = 2.1 μm, 线阵中心的差值W = 6.3 μm. RW是与周期有关的参数, 即R1 = P/10, W = 3R1. 周期Px = Py = 2P = 42 μm, P = 21 μm. 线阵在y方向的长度为L = 37.8 μm, 其值也与P相关, L = KP/10, 其中K表示权重, K = 9. 这样设置线阵的参数有利于我们结构的稳定, 通过改变半周期P与权重K的值的大小, 可以调整微结构之间的比例关系. 上述的这些参数的设定可以使得我们所提出的吸收器在4—14.5 THz波段实现良好的电磁响应.

      本文使用CST仿真软件, 并基于FDID算法进行了模拟仿真. 通过在XY方向设置周期性的边界条件, 在Z方向上设置开放的边界条件, 让电磁波沿Z轴的负方向正入射到该吸收器的表面 [2931]. 我们利用单个周期单元实现对全部周期的仿真, 在保证结果正确的情况下, 加快了计算速度, 得到S参数. 利用S参数, 分别计算了吸收率(A)和传输率(T ), 其中, A = 1 – R – T, 其中R为反射率, R = $ |{S}_{11}| $2, 而传输率T = $ |{S}_{21}| $2. 由于底部Au的厚度远大于THz波对Au的趋肤深度, 可以认为传输接近于0, 因此, 该吸收器的吸收率可以简化为A = 1 – R [32,33].

      本吸收器的金属膜衬底可以通过沉积工艺实现, 接着可利用化学蒸汽法形成SiO2层. 周期性阵列可以通过自组装方法实现. 之前关于单层阵列的金属纳米线的自组装的报道具有制造可行性 [34]. 例如, 导电纳米线可以通过施加交变电场进行自组装 [35]. 此外, 利用嵌段共聚物, 可以实现介电圆柱形谐振器从垂直方向到平行排列的方向转变[36]. 而金属辅助化学刻蚀和组装技术可以为纳米结构引入新的制造方法[37].

      利用随机相位近似理论(RPA)获得BDS的复导电率, BDS的动态导电率表示为 [23]

      式中, G(E ) = n(–E ) – n(E ), 其中n(E )为费米分布函数; $ \varOmega = {\hbar}\omega / {E}_{{\mathrm{F}}} +{\rm i}{\hbar} \tau -1/ {E}_{{\mathrm{F}}} $; $ {k}_{{\mathrm{F}}} = {E}_{{\mathrm{F}}} /({\hbar} {v}_{{\mathrm{F}}}) $为费米动量, $ {v}_{{\mathrm{F}}} \approx {10}^{6} \;{\rm m/s}$表示费米速度; $ \varepsilon = E/ {E}_{{\mathrm{F}}} $, $ {\varepsilon }_{{\mathrm{c}}} = {E}_{{\mathrm{c}}}/ {E}_{{\mathrm{F}}} $ ($ {E}_{{\mathrm{c}}} $为狄拉克谱, 不是线性的截止能量); 简并因子g = 40. 在带电子跃迁的情况下, 可以通过使用双带模型来获得BDS的介电常数[24]:

      式中, $ {\varepsilon }_{0} $为真空介电常数, $ {\varepsilon }_{{\mathrm{b}}} $ = 1为有效背景介电常数. 将(1)式和(2)式代入(3)式, 得到不同费米能级下BDS介电常数实部和虚部随频率的变化, 如图2所示.

    • 通过器件结构参数优化, 得到完美吸收器在4—14.5 THz范围内的特征光谱, 如图3(a)所示. 该吸收器可实现在f1 = 5.032 THz, f2 = 5.859 THz, f3 = 7.674 THz, f4 = 9.654 THz, f5 = 11.656 THz, f6 = 12.514 THz, f7 = 14.01 THz一共7个频率点的高吸收. 从图3(a)可以看到, 其中有6个频率点的吸收率都超过了90%. 这里f1处的吸收率为84.43%, f2处的吸收率为96.23%, f3处的吸收率为91.36%, f4处的吸收率为99.02%, f5处的吸收率93.84%, f6处的吸收率为98.47%, f7处的吸收率为97.32%. 该结果说明该吸收器对外部正入射的电磁波具有良好的特性响应. 此外, 图1表明该器件的结构并不是旋转对称的, 所以可以预测到对于TE和TM两种极化模式下的光谱图并不一致. 说明该器件具有极化相干性.

      为了解释上述结果, 先从宏观电磁学的层面, 利用阻抗匹配理论对其进行分析. 相对阻抗Z的表达式为[38,39]

      图3(b)所示, 通过参数反演给出了该吸收器在4—14.5 THz范围内相对阻抗的实部和虚部曲 线. 根据阻抗匹配理论, 当吸收器的等效阻抗与自由空间的阻抗Z0相匹配时, 吸收器可以表现出良好的电磁响应, 这里Z0 = 376.7 Ω. 此时吸收器的相对阻抗Z应等于1. 从图3(b)可以观察到, 在7个共振频率点处, 吸收器的相对阻抗的实部靠近1, 虚部靠近0. 这表明该吸收器在这7个频率点处达到临界耦合条件, 系统的反射被抑制, 从而显示出强电磁吸收特性[40]. 图3(b)的理论推导与图3(a)的仿真结果均表明在这7个频率点发生完美吸收.

      为了进一步说明该吸收器背后的物理机理, 如图4图5所示, 通过添加场监视器分别探究了该吸收器在7个频率点处的超材料表面电场和整体的磁场分布情况. 图4展示了电场能量聚集的区域, 表征共振强度的分布情况. 而图5则展示出磁能分布, 表征了能量的变化趋势. 图4(a)(g)给出了上述7个波段处顶部电场分布情况. 从f1f7的表面电场分布情况可以发现一个统一的规律, 电场主要分布在BDS线阵结构的表面, 尤其是在y方向上的BDS线阵的端面处, 电场更加集中. 这一现象可利用Fano共振的等离子体来进行理解[41,42]. 7个频率点的电磁波入射到器件的表面激发了等离子体. 并在此区域发生了强烈的Fano共振, 入射电磁波的能量被转化成微结构表面的等离子体振动能量, 等离子体Fano共振具有较窄的光谱线宽, 所以我们才能看到电场在该区域高度集中, 吸收光谱的线宽较窄[43]. 同时, 电场分布的特性也表明共振吸收是沿着纵向方向激发的. 此外, 还可以观察到部分电场集中于BDS线阵的间隔区域, 电场同样集中于相邻线阵的边缘处. 这很明显是偶极性等离子体共振和相邻的等离子体的近场耦合效应的表现形式. 特别地, f4相对于f1而言, 耦合效益最为强烈且集中, 这也是f4吸收率最高且峰值最窄的原因. 而f7则是偶极性等离子体共振最为强烈.

      接着又探究了该吸收器在7个共振频率点处的磁场分布, 如图5(a)(g)所示. 当观察y-z截面处的磁场分布时, 可以看出当电磁波正入射到吸收器上时, 电磁波将通过线阵结构之间的缝隙进入到吸收器中间的SiO2介质层中. 于是底部的Au和顶部的BDS就会形成法布里-珀罗谐振腔[44,45]. 电磁波在其中不断地来回反射, 电磁波的能量由此被消耗掉. 此外可以观察到在线阵的顶部存在磁能的聚集, 顶部磁场的强度分布也证实了吸收是由顶部的Fano共振引起的. 同时证实了共振所导致的吸收是线阵的纵向特性. 我们还可以观察到有一部分磁能被限制在BDS的内部. 这是因为线阵是存在间隙的, 腔内反射的能量与吸收器周期边缘入射的能量相互耦合使得能量不再被消耗于介质层内, 而部分被分散到BDS中. 如图5所示, 随着频率的增大, 电磁波的能量会随着变大, Fano共振以及耦合的程度也会加深.

      品质因子(Q值)是谐振模式电磁存储能力的一个重要评价指标, Q值越大往往代表着吸收带宽越小, 探测更加灵敏. Q值的表达式为 [46,47]

      式中FWHM表示半峰全宽, 7个频率点的FWHM值见补充材料表S1 (online), 对应的Q值也如表S1(online)所列. 其中f4频率点处吸收峰的品质因子最大为219.41, 这非常有利于高精度的信号聚集与探测. 本文吸收器的Q值与近年来类似吸收器之间的比较见表1, 可以看出, 我们所提出的吸收器对于能量的集中和探测具有优势[4852].

      此外, 灵敏度作为传感器的一个主要指标会影响吸收器的应用范围. 为此, 探究了吸收器的吸收率随外部折射率的变化情况, 如图6所示. 图6(a)是外部折射率分别为1, 1.02, 1.04, 1.06, 1.08, 1.1时吸收器的吸收峰值的变化情况. 图6(b)是吸收频点随外部折射率的变化情况. 当外部折射率增大时, 吸收器的7个吸收峰都发生了红移, 这使得吸收器具有跟随外部环境来动态调节吸收峰的能力. 图6(c)是7个吸收峰的峰值随外部折射率的变化情况. 可以看出, 在外部折射率变大时, f1频点处的吸收率表现出先减少后上升的趋势, 整体表现出下降的趋势; f2频点处的吸收率的变化趋势则与f1相反; f3频点处的吸收率则是先下降后上升, 但整体吸收率略微升高; f5频点处的吸收率整体表现出上升的趋势; f6频点适合于固定的检测. 而f4f7频点处的吸收率都表现出快速的下降; f7频点处吸收率的下降趋势最为强烈, 表明f7受外部折射率的影响最大.

      结合图6(c)的变化趋势和图6(b)拟合出的直线方程, 可以得到评判吸收器应用领域的另一个非常重要的指标——灵敏度, 其表达式如下[53]:

      式中$ {{\Delta }}f $表示频点的变化, $ {{\Delta }}n $表示外部折射率的变化. 根据图6(b)拟合直线的斜率可以得到7个频点的灵敏度, 见补充材料表S2 (online). 同时, 在频率增加的过程中, 伴随着Fano共振以及耦合程度的加深, 介质层中的能量逐渐转移到器件顶部的BDS结构中. 其中在f7频点(14.01 THz)处能量被分散的最为强烈, 所受影响最大. 同时当外部折射率增大时, 进入到介质层的光场减小, 这也会引起吸收率的下降.

      吸收器的另一个重要指标是FOM值(品质因数), 其表达式如下:

      7个频点的FOM值也列于补充材料表S2 (online), 可以看到7种吸收模式下FOM值最高的是f7, 达到了35.204 RIU–1. 上述结果表明, 该吸收器的特性与外部折射率紧密相关. 表2对比了近些年所提出的吸收器[51, 54-57], 可以看出我们所提出的工作具有更优异的灵敏度.

      以上所有讨论都是在理想结构参数下的吸收器对外部正入射的电磁波的响应状态. 但是在实际操作中, 往往会存在对结构大小和形状的制造误差, 以及对器件是否具有动态调节能力的需求等. 因此我们所提出的吸收器是否具有制造容忍性、使用稳定性和良好的调节性将决定了吸收器的应用范围. 在其他参数保持在最佳水平的条件下, 研究每个特定参数对器件性能的影响, 具体结果见补充材料图S1和图S2 (online).

      BDS的一个显著特性是可以通过调整费米能量EF来改变其电学特性[2528]. 可以通过在吸收器表面涂上等离子体凝胶来施加外部偏压或者化学掺杂实现对BDS费米能的改变. 图7给出了不同EF下器件的吸收率. 从图7(b)可看到, 随着BDS费米能的增大, 7个频点的位置都发生了蓝移. 同时从图7(c)可以看到, 在这一过程中, f4, f5, f7频点处的吸收率在提升, 在Ef = 90 meV时接近于100%, f1f3频点处的吸收率虽然有些下降, 但在EF = 90 meV时仍保持着79.4%和80.24%的吸收率. 整个过程中, f2频点处的吸收率一直保持在90%以上, f6频点处的吸收率先基本保持不变, 在EF = 90 meV时降至87.11%. 综合来看, 利用BDS的费米能可以很好地实现吸收器的化学调谐性能.

      除了对上述各种参数的讨论外, 还需要考虑在实际应用环境中, 往往不是垂直入射的平面波, 为此, 需要探究不同入射角的电磁波对吸收效应的影响, 具体结果见补充材料图S3 (online), 该吸收器在0—10 THz波段内的吸收率在50°之前都是很优异的.

    • 本文提出了一种基于BDS纳米线的七波段太赫兹可调完美吸收器. 该吸收器可实现对f1 = 5.032 THz, f2 = 5.859 THz, f3 = 7.674 THz, f4 = 9.654 THz, f5 = 11.656 THz, f6 = 12.514 THz, f7 = 14.01 THz一共7个频率点的高吸收. 先从阻抗匹配理论的角度出发证实了结果的正确性. 接着通过分析7个频率点处的表面电场和磁场, 利用Fano共振原理解释了吸收形成背后的物理机理. 计算了7个吸收峰的Q值, 最大品质因子为219.41. 通过改变外部折射率的变化探究了吸收器的最大灵敏度为5421.43 GHz/RIU, 最大品质因数FOM值为35.204. 然后又讨论了关键参数对吸收器吸收率的影响, 说明该吸收器具有双固定性能的选择能力. 通过外部施加偏置电压或是化学掺杂改变BDS的费米能以实现吸收器的动态协调. 最后为了使该器件可以应用在合适的场合, 探究了外部电磁波的入射角对吸收率的影响, 表明该吸收器在中低频对外部电磁波的入射角的改变可稳定到50°之前, 但在高频处却只能维持到20°. 综上所述, 本器件具备简单结构、完美吸收、适用于中低频并能够额外检测高频的特点, 在传感和光电相关器件领域具有很好的参考价值.

    参考文献 (57)

目录

/

返回文章
返回