搜索

x

留言板

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

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

含石墨烯分界面有耗分层介质的传播矩阵

王飞 魏兵

含石墨烯分界面有耗分层介质的传播矩阵

王飞, 魏兵
PDF
HTML
导出引用
  • 给出一种适用于含导电界面的有耗分层介质的传播矩阵方法. 利用相位匹配原理给出斜入射时有耗介质波矢量的实部和虚部, 二者方向不同使得在介质中传播非均匀平面波. 根据边界条件, 推导了跨越石墨烯界面的传播矩阵, 以及“无限薄”石墨烯层的反透射系数解析式. 最终将传播矩阵方法推广应用于含石墨烯界面的有耗分层介质情形, 可用于快速解析分析分层介质与导电界面复合结构的反透射和电波传播特性.
      通信作者: 王飞, wfei79@163.com
    • 基金项目: 国家自然科学基金(批准号: 61401344, 61571348)和高等学校学科创新引智计划(批准号: B17035)资助的课题
    [1]

    Novoselov K S, Geim A K, Morozov S V, Jiang D, Zhang Y, Dubonos S V, Grigorieva I V, Firsove A A 2004 Science 306 666

    [2]

    Geim A K 2009 Science 324 1530

    [3]

    Sensale-Rodriguez B, Yan R, Kelly M, Fang T, Tahy K, Hwang W S, Jena D, Liu L, Xing H G 2012 Nature Commun. 3 780

    [4]

    Alaee R, Farhat M, Rockstuhl C, Lederer F 2012 Opt. Express 20 28017

    [5]

    Fallahi A, Perruisseau-Carrier J 2012 Phys. Rev. B 86 195408

    [6]

    Sensale-Rodriguez B, Yan R, Rafique S, Zhu M, Li W, Liang X, Gundlach D, Protasenko V, Kelly M M, Jena D, Liu L, Xing H G 2012 Nano Lett. 12 4518

    [7]

    Fu M X, Zhang Y 2013 JECT 11 352

    [8]

    Lee S H, Choi M, Kim T T, Lee S, Liu M, Yin X, Choi H K, Lee S S, Choi C G, Choi S Y, Zhang X, Min B 2012 Nat. Mater. 11 936

    [9]

    Zuo Z G, Wang P, Ling F R, Liu J S, Yao J Q 2013 Chin. Phys. B 22 097304

    [10]

    张玉萍, 张洪艳, 尹贻恒, 刘陵玉, 张晓, 高营, 张会云 2012 物理学报 61 047803

    Zhang Y P, Zhang H Y, Yin Y H, Liu L Y, Zhang X, Gao Y, Zhang H Y 2012 Acta Phys. Sin. 61 047803

    [11]

    Zhu X L, Yan W, Jepsen P U, Hansen O, Mortensen N A Xiao S S 2013 Appl. Phys. Lett. 102 131101

    [12]

    Pomar J L G, Alexey Y N, Luis M M 2013 ACS Nano 7 4988

    [13]

    Thongrattanasiri S, Koppens F H L, de Abajo F J G 2012 Phys. Rev. Lett. 108 047401

    [14]

    Ferreira A, Peres N M R, Ribeiro R M, Stauber T 2012 Phys. Rev. B 85 115438

    [15]

    Tian Y C, Jia W, Ren P W, Fan C Z 2018 Chin. Phys. B 27 124205

    [16]

    Jia W, Ren P W, Fan C Z, Tian Y C 2019 Chin. Phys. B 28 026102

    [17]

    Liu J T, Liu N H, Li J, Li X J, Huang J H 2012 Appl. Phys. Lett. 101 052104

    [18]

    Peres N M R, Bludov Y V 2013 EPL 101 58002

    [19]

    谢凌云, 肖文波, 黄国庆, 胡爱荣, 刘江涛 2014 物理学报 63 057803

    Xie L Y, Xiao W B, Huang G Q, Hu A R, Liu J T 2014 Acta Phys. Sin. 63 057803

    [20]

    Zhang H J, Zheng G G, Chen Y Y 2018 Chin. Phys. Lett. 35 038102

    [21]

    Lovat G 2012 IEEE Trans. Electromagn. Compat. 54 101

    [22]

    葛德彪, 魏兵 2011 电磁波理论 (北京: 科学出版社) 第32, 56−65页

    Ge D B, Wei B 2011 Electromagnetic Wave Theory (Beijing: Science Press) pp32, 56−65 (in Chinese)

    [23]

    George W H 2008 J. Appl. Phys. 103 064302

    [24]

    孙旺, 李粮生, 张景, 殷红成 2018 雷达学报 7 67

    Sun W, Li L S, Zhang J, Yin H C 2018 J. Radars 7 67

  • 图 1  分层介质

    Fig. 1.  Stratified medium.

    图 2  石墨烯电导率和SE (a) 电导率实虚部; (b) SE

    Fig. 2.  Graphene conductivity and shielding effectiveness: (a) Real and imaginary parts of conductivity; (b) SE.

    图 3  石墨烯单层的反射和透射系数模值 (a) TE模; (b) TM模

    Fig. 3.  Modulus of reflective and transmittance coefficients of a graphene sheet:(a) TE mode; (b) TM mode.

    图 4  含石墨烯涂层CdTe半空间的反透射系数模值(TE模) (a) 反射系数; (b) 透射系数

    Fig. 4.  Modulus of reflective and transmittance coefficients of CdTe half-space containing graphene coating (TE mode): (a) Reflective coefficient; (b) transmittance coefficient.

    图 5  含石墨烯涂层CdTe半空间的反透射系数模值(TM模) (a) 反射系数; (b) 透射系数

    Fig. 5.  Modulus of reflective and transmittance coefficients of CdTe half-space containing graphene coating (TM mode): (a) Reflective coefficient; (b) transmittance coefficient.

    图 6  CdTe半空间的反透射光场(TE模) (a) 无石墨烯涂层; (b) 含石墨烯涂层

    Fig. 6.  Optical field of reflection and transmission coefficients of CdTe half-space (TE mode): (a) Without graphene coating; (b) with graphene coating.

    图 7  CdTe半空间的反透射光场(TM模) (a) 无石墨烯涂层; (b) 含石墨烯涂层

    Fig. 7.  Optical field of reflection and transmission coefficients of CdTe half-space (TM mode): (a) Without graphene coating; (b) with graphene coating.

    图 8  Si/SiO2周期结构型1DPC

    Fig. 8.  Si/SiO2 1DPC with periodic structure.

    图 9  含石墨烯涂层Si/SiO2周期结构1DPC的反透射系数 (a) 反射系数; (b) 透射系数

    Fig. 9.  Modulus of reflective and transmittance coefficients of Si/SiO2 1DPC containing graphene sheet: (a) Reflective coefficient; (b) transmittance coefficient.

    图 10  含石墨烯界面Si/SiO2周期结构1DPC的吸收率

    Fig. 10.  Absorbance of Si/SiO2 1DPC containing graphene sheet.

    图 11  含石墨烯涂层Si/SiO2周期结构1DPC的吸收率(TE模) (a) 无涂层; (b) 表面涂层; (c) 底层涂层

    Fig. 11.  Contour plots of the absorbance of the Si/SiO2 1DPC as a function of the light frequency and the incident angles for the TE mode: (a) Without graphene sheet; (b) graphene sheet on the top; (c) graphene sheet on the bottom.

    图 12  含石墨烯涂层Si/SiO2周期结构1DPC的吸收率(TM模) (a) 无涂层; (b) 表面涂层; (c) 底层涂层

    Fig. 12.  Contour plots of the absorbance of the Si/SiO2 1DPC as a function of the light frequency and the incident angles for the TE mode: (a) Without graphene sheet; (b) graphene sheet on the top; (c) graphene sheet on the bottom.

  • [1]

    Novoselov K S, Geim A K, Morozov S V, Jiang D, Zhang Y, Dubonos S V, Grigorieva I V, Firsove A A 2004 Science 306 666

    [2]

    Geim A K 2009 Science 324 1530

    [3]

    Sensale-Rodriguez B, Yan R, Kelly M, Fang T, Tahy K, Hwang W S, Jena D, Liu L, Xing H G 2012 Nature Commun. 3 780

    [4]

    Alaee R, Farhat M, Rockstuhl C, Lederer F 2012 Opt. Express 20 28017

    [5]

    Fallahi A, Perruisseau-Carrier J 2012 Phys. Rev. B 86 195408

    [6]

    Sensale-Rodriguez B, Yan R, Rafique S, Zhu M, Li W, Liang X, Gundlach D, Protasenko V, Kelly M M, Jena D, Liu L, Xing H G 2012 Nano Lett. 12 4518

    [7]

    Fu M X, Zhang Y 2013 JECT 11 352

    [8]

    Lee S H, Choi M, Kim T T, Lee S, Liu M, Yin X, Choi H K, Lee S S, Choi C G, Choi S Y, Zhang X, Min B 2012 Nat. Mater. 11 936

    [9]

    Zuo Z G, Wang P, Ling F R, Liu J S, Yao J Q 2013 Chin. Phys. B 22 097304

    [10]

    张玉萍, 张洪艳, 尹贻恒, 刘陵玉, 张晓, 高营, 张会云 2012 物理学报 61 047803

    Zhang Y P, Zhang H Y, Yin Y H, Liu L Y, Zhang X, Gao Y, Zhang H Y 2012 Acta Phys. Sin. 61 047803

    [11]

    Zhu X L, Yan W, Jepsen P U, Hansen O, Mortensen N A Xiao S S 2013 Appl. Phys. Lett. 102 131101

    [12]

    Pomar J L G, Alexey Y N, Luis M M 2013 ACS Nano 7 4988

    [13]

    Thongrattanasiri S, Koppens F H L, de Abajo F J G 2012 Phys. Rev. Lett. 108 047401

    [14]

    Ferreira A, Peres N M R, Ribeiro R M, Stauber T 2012 Phys. Rev. B 85 115438

    [15]

    Tian Y C, Jia W, Ren P W, Fan C Z 2018 Chin. Phys. B 27 124205

    [16]

    Jia W, Ren P W, Fan C Z, Tian Y C 2019 Chin. Phys. B 28 026102

    [17]

    Liu J T, Liu N H, Li J, Li X J, Huang J H 2012 Appl. Phys. Lett. 101 052104

    [18]

    Peres N M R, Bludov Y V 2013 EPL 101 58002

    [19]

    谢凌云, 肖文波, 黄国庆, 胡爱荣, 刘江涛 2014 物理学报 63 057803

    Xie L Y, Xiao W B, Huang G Q, Hu A R, Liu J T 2014 Acta Phys. Sin. 63 057803

    [20]

    Zhang H J, Zheng G G, Chen Y Y 2018 Chin. Phys. Lett. 35 038102

    [21]

    Lovat G 2012 IEEE Trans. Electromagn. Compat. 54 101

    [22]

    葛德彪, 魏兵 2011 电磁波理论 (北京: 科学出版社) 第32, 56−65页

    Ge D B, Wei B 2011 Electromagnetic Wave Theory (Beijing: Science Press) pp32, 56−65 (in Chinese)

    [23]

    George W H 2008 J. Appl. Phys. 103 064302

    [24]

    孙旺, 李粮生, 张景, 殷红成 2018 雷达学报 7 67

    Sun W, Li L S, Zhang J, Yin H C 2018 J. Radars 7 67

  • [1] 王飞, 魏兵. 电磁偏置各向异性石墨烯界面的传播矩阵. 物理学报, 2020, (): . doi: 10.7498/aps.69.20201089
    [2] 王飞, 魏兵. 分层有耗手征介质中斜入射电磁波的传播矩阵. 物理学报, 2017, 66(6): 064101. doi: 10.7498/aps.66.064101
    [3] 邓伟胤, 朱瑞, 邓文基. 有限尺寸石墨烯的电子态. 物理学报, 2013, 62(8): 087301. doi: 10.7498/aps.62.087301
    [4] 彭艳玲, 薛文瑞, 卫壮志, 李昌勇. 涂覆石墨烯的非对称并行电介质纳米线波导的模式特性分析. 物理学报, 2018, 67(3): 038102. doi: 10.7498/aps.67.20172016
    [5] 卫壮志, 薛文瑞, 彭艳玲, 程鑫, 李昌勇. 基于涂覆石墨烯的三根电介质纳米线的THz波导的模式特性分析. 物理学报, 2018, 67(10): 108101. doi: 10.7498/aps.67.20180036
    [6] 程鑫, 薛文瑞, 卫壮志, 董慧莹, 李昌勇. 涂覆石墨烯的椭圆形电介质纳米线光波导的模式特性分析. 物理学报, 2019, 68(5): 058101. doi: 10.7498/aps.68.20182090
    [7] 张保磊, 王家序, 肖科, 李俊阳. 石墨烯-纳米探针相互作用有限元准静态计算. 物理学报, 2014, 63(15): 154601. doi: 10.7498/aps.63.154601
    [8] 张玉萍, 刘陵玉, 陈琦, 冯志红, 王俊龙, 张晓, 张洪艳, 张会云. 具有分离门电抽运石墨烯中电子-空穴等离子体的冷却效应. 物理学报, 2013, 62(9): 097202. doi: 10.7498/aps.62.097202
    [9] 俎凤霞, 张盼盼, 熊伦, 殷勇, 刘敏敏, 高国营. 以石墨烯为电极的有机噻吩分子整流器的设计及电输运特性研究. 物理学报, 2017, 66(9): 098501. doi: 10.7498/aps.66.098501
    [10] 郭伟玲, 邓杰, 王嘉露, 王乐, 邰建鹏. 具有石墨烯/铟锑氧化物复合透明电极的GaN发光二极管. 物理学报, 2019, 68(24): 247303. doi: 10.7498/aps.68.20190983
    [11] 乔文涛, 龚健, 张利伟, 王勤, 王国东, 廉书鹏, 陈鹏辉, 孟威威. 梳状波导结构中石墨烯表面等离子体的传播性质. 物理学报, 2015, 64(23): 237301. doi: 10.7498/aps.64.237301
    [12] 邓新华, 刘江涛, 袁吉仁, 王同标. 全新的电导率特征矩阵方法及其在石墨烯THz频率光学特性上的应用. 物理学报, 2015, 64(5): 057801. doi: 10.7498/aps.64.057801
    [13] 卢晓波, 张广宇. 石墨烯莫尔超晶格. 物理学报, 2015, 64(7): 077305. doi: 10.7498/aps.64.077305
    [14] 秦志辉. 类石墨烯锗烯研究进展. 物理学报, 2017, 66(21): 216802. doi: 10.7498/aps.66.216802
    [15] 谢凌云, 肖文波, 黄国庆, 胡爱荣, 刘江涛. 光子晶体增强石墨烯THz吸收. 物理学报, 2014, 63(5): 057803. doi: 10.7498/aps.63.057803
    [16] 叶振强, 曹炳阳, 过增元. 石墨烯的声子热学性质研究. 物理学报, 2014, 63(15): 154704. doi: 10.7498/aps.63.154704
    [17] 金芹, 董海明, 韩奎, 王雪峰. 石墨烯超快动态光学性质. 物理学报, 2015, 64(23): 237801. doi: 10.7498/aps.64.237801
    [18] 卢琪, 吕宏鸣, 伍晓明, 吴华强, 钱鹤. 石墨烯射频器件研究进展. 物理学报, 2017, 66(21): 218502. doi: 10.7498/aps.66.218502
    [19] 张晓波, 青芳竹, 李雪松. 化学气相沉积石墨烯薄膜的洁净转移. 物理学报, 2019, 68(9): 096801. doi: 10.7498/aps.68.20190279
    [20] 韩同伟, 贺鹏飞. 石墨烯弛豫性能的分子动力学模拟. 物理学报, 2010, 59(5): 3408-3413. doi: 10.7498/aps.59.3408
  • 引用本文:
    Citation:
计量
  • 文章访问数:  863
  • PDF下载量:  16
  • 被引次数: 0
出版历程
  • 收稿日期:  2019-05-28
  • 修回日期:  2019-09-29
  • 上网日期:  2019-11-27
  • 刊出日期:  2019-12-01

含石墨烯分界面有耗分层介质的传播矩阵

  • 1. 西安电子科技大学物理与光电工程学院, 西安 710071
  • 2. 西安电子科技大学信息感知技术协同创新中心, 西安 710071
  • 通信作者: 王飞, wfei79@163.com
    基金项目: 国家自然科学基金(批准号: 61401344, 61571348)和高等学校学科创新引智计划(批准号: B17035)资助的课题

摘要: 给出一种适用于含导电界面的有耗分层介质的传播矩阵方法. 利用相位匹配原理给出斜入射时有耗介质波矢量的实部和虚部, 二者方向不同使得在介质中传播非均匀平面波. 根据边界条件, 推导了跨越石墨烯界面的传播矩阵, 以及“无限薄”石墨烯层的反透射系数解析式. 最终将传播矩阵方法推广应用于含石墨烯界面的有耗分层介质情形, 可用于快速解析分析分层介质与导电界面复合结构的反透射和电波传播特性.

English Abstract

    • 2004年, 世界上最薄的半导体材料—石墨烯, 由Novoselov和Geim的团队采用微机械剥离法制备成功[1]. 石墨烯是由单层碳原子构成的二维晶体结构, 具有高电子迁移率、低电阻率、光学透明性、高灵活性及结构稳定等特性[2], 这使得其在光学和光电传感等方面具有广泛的应用前景[3-6], 同时也被认为非常适合于太赫兹光学设备[7-10]. 石墨烯非常薄, 厚度只有约0.34 nm, 这导致其光吸收率很低, 应用受到限制. 为此, 常采用周期掺杂、表面等离子体、光学微腔等方法来增强光学吸收[11-16], 其中利用光子晶体增强石墨烯光吸收是一种常用方法. 2012年, Liu等[17]利用石墨烯和间隔层在一维光子晶体(one-dimensional photonic crystal, 1DPC)表面构造表面缺陷, 导致光的局域化, 使可见光光吸收被增强约4倍. 2013年, Peres和Bludov[18]利用1DPC使石墨烯的THz吸收增强了超过3倍. 2014年, Xie等[19]分析了1DPC表面石墨烯在应力赝磁场作用下的THz吸收. 2018年, Zhang等[20]设计了一种高对比光栅耦合波导结构, 提高了石墨烯在近红外和可见光频段的光吸收.

      1DPC本质上是一个周期分层介质, 而单层的石墨烯可以看做是“无限薄”的[21], 因此可以把石墨烯和1DPC构造的结构看成一个含石墨烯分界面的分层介质. 传播矩阵方法[22]是一种分析分层介质中光或电磁波传播特性的解析理论方法, 可以直接计算分层介质对平面电磁波的反透射等. 根据电磁理论, 当分层介质中某层为有耗介质时, 该层中的波矢量为复数, 而在斜入射情形下, 考虑相位匹配, 该层中波矢量的实部和虚部并不平行, 即相位和振幅的传播方向不同, 此时该层中传播的应为非均匀平面波. 本文首先根据相位匹配条件, 推导了有耗分层介质中复数波矢量的实部和虚部, 给出传播矩阵中所需要的纵向波矢量计算式; 然后考虑石墨烯表面电导率, 根据边界条件推导了跨越石墨烯分界面的传播矩阵; 最终, 使得传播矩阵方法可用于解析分析任意入射角情形下含导电界面有耗分层介质的反透射和电波传播特性. 本文中时谐因子取$\exp \left( { - {\rm{j}}\omega t} \right)$.

    • 设平面波从真空入射到$n$层分层介质, 如图1所示, 入射角为$\theta $. 考虑一般情形, 设各层介质的相对介电系数和磁导系数${\varepsilon _{{\rm{r}}l}}, {\mu _{{\rm{r}}l}}\;(l = 1, 2, \cdots, n)$皆为复数,

      图  1  分层介质

      Figure 1.  Stratified medium.

      ${\varepsilon _{{\rm{r}}l}} = {\varepsilon '_{{\rm{r}}l}} + {\rm{j}}{\varepsilon ''_{{\rm{r}}l}},\;{\mu _{{\rm{r}}l}} = {\mu '_{{\rm{r}}l}} + {\rm{j}}{\mu ''_{{\rm{r}}l}}.$

      此时, 介质有耗, 介质层中的波矢量${{{k}}_l}$也为复数. 下面来具体分析介质层中的复数波矢量${{{k}}_l}$.

      各分层中的复数波矢量${{{k}}_l}$可以写为

      ${{{k}}_l} = {{{k}}_{{\rm{R}}l}} + {\rm{j}}{{{k}}_{{\rm{I}}l}},$

      其中${{{k}}_{{\rm{R}}l}}$${{{k}}_{{\rm{I}}l}}$分别表示${{{k}}_l}$的实部和虚部. 如图1所示, 入射面为$xoz$平面, 故有

      $\begin{split} & {{{k}}_{{\rm{R}}l}} = \hat z{k_{{\rm{R}}lz}} + \hat x{k_{{\rm{R}}lx}}, \\ & {{{k}}_{{\rm{I}}l}} = \hat z{k_{{\rm{I}}lz}} + \hat x{k_{{\rm{I}}lx}}. \end{split} $

      根据界面处的相位匹配(波矢量切向分量连续)条件可得

      ${k_{{\rm{R}}lx}} = {k_{0x}} = {k_0}\sin \theta,{\rm{ }}{k_{{\rm{I}}lx}} = 0.$

      ${{{k}}_l}$的虚部${{{k}}_{{\rm{I}}l}}$只有$z$分量, 实部${{{k}}_{{\rm{R}}l}}$$x$分量在每个区域都相同, 都等于区域0中入射波矢量${{{k}}_0}$$x$分量${k_{0 x}}$. 由此可以看到, 平面波斜入射有耗分层介质时, 介质中波矢量的实部和虚部方向不平行, 虚部只有纵向分量. 由于波矢量实部代表相位传播方向, 虚部代表振幅传播方向, 可知此时平面波是非均匀的, 波矢量为

      $\begin{split} {{{k}}_l} ={}& \hat z\left( {{k_{{\rm{R}}lz}} + {\rm{j}}{k_{{\rm{I}}l}}} \right) + \hat x{k_{{\rm{R}}lx}} \\ ={}& \hat z\left( {{k_{{\rm{R}}lz}} + {\rm{j}}{k_{{\rm{I}}l}}} \right) + \hat x{k_{0x}}. \end{split} $

      由电磁波理论可知在有耗介质中[22]

      $\begin{split} {{{k}}_l} \cdot {{{k}}_l} ={}& \omega {\varepsilon _0}{\mu _0}{\varepsilon _{{\rm{r}}l}}{\mu _{{\rm{r}}l}} \\ ={}& k_0^2\left( {{{\varepsilon '}_{{\rm{r}}l}} + {\rm{j}}{{\varepsilon ''}_{{\rm{r}}l}}} \right)\left( {{{\mu '}_{{\rm{r}}l}} + {\rm{j}}{{\mu ''}_{{\rm{r}}l}}} \right). \end{split}$

      而由(2)和(5)式可得

      $\begin{split} {{{k}}_l} \cdot {{{k}}_l}\, &= \left( {{{{k}}_{{\rm{R}}l}} + {\rm{j}}{{{k}}_{{\rm{I}}l}}} \right) \cdot \left( {{{{k}}_{{\rm{R}}l}} + {\rm{j}}{{{k}}_{{\rm{I}}l}}} \right) \\ &= {{{k}}_{{\rm{R}}l}} \cdot {{{k}}_{{\rm{R}}l}} - {{{k}}_{{\rm{I}}l}} \cdot {{{k}}_{{\rm{I}}l}} + 2{\rm{j}}\left( {{{{k}}_{{\rm{R}}l}} \cdot {{{k}}_{{\rm{I}}l}}} \right) \\ & = k_{{\rm{R}}lz}^2 - k_{{\rm{I}}l}^2 + 2{\rm{j}}\left( {{k_{{\rm{R}}lz}} \cdot {k_{{\rm{I}}l}}} \right) + k_{0x}^2. \end{split} $

      则根据(6)和(7)式得

      $\begin{split} & k_{{\rm{R}}lz}^2 - k_{{\rm{I}}l}^2 = k_0^2\left( {{{\varepsilon '}_{{\rm{r}}l}}{{\mu '}_{{\rm{r}}l}} - {{\varepsilon ''}_{{\rm{r}}l}}{{\mu ''}_{{\rm{r}}l}}} \right) - k_0^2{\sin ^2}\theta, \\ & 2{k_{{\rm{R}}lz}} \cdot {k_{{\rm{I}}l}} = k_0^2\left( {{{\varepsilon '}_{{\rm{r}}l}}{\mu ''}_{\rm rl} + {{\varepsilon ''}_{{\rm{r}}l}}{{\mu '}_{{\rm{r}}l}}} \right). \end{split} $

      $\begin{split} & AA = ({{\varepsilon '}_{{\rm{r}}l}}{{\mu '}_{{\rm{r}}l}} - {{\varepsilon ''}_{{\rm{r}}l}}{{\mu ''}_{{\rm{r}}l}}) - {\sin ^2}\theta,\\ & BB = \left( {{{\varepsilon '}_{{\rm{r}}l}}{{\mu ''}_{{\rm{r}}l}} + {{\varepsilon ''}_{{\rm{r}}l}}{{\mu '}_{{\rm{r}}l}}} \right), \\ & CC = {{BB} / {AA}}. \end{split} $

      求解(8)式得到

      $\left\{ \begin{split} & {k_{{\rm{I}} l}} = {k_{{\rm{R}}lz}} = {k_0}\sqrt {\frac{{BB}}{2}},\quad\quad\quad\quad\quad\quad\quad\;\; AA = 0; \\ & {k_{{\rm{I}} l}} = {k_0}\sqrt {AA \cdot \frac{{\sqrt {1 + {{(CC)}^2}} - 1}}{2}},\\ & {k_{{\rm{R}}lz}} = {k_0}\sqrt {AA \cdot \frac{{\sqrt {1 + {{(CC)}^2}} + 1}}{2}},\quad\;\;\;\; AA > 0; \\ & {k_{{\rm{I}} l}} = {k_0}\sqrt { - AA \cdot \frac{{\sqrt {1 + {{(CC)}^2}} + 1}}{2}},\\ & {k_{{\rm{R}}lz}} = {k_0}\sqrt { - AA \cdot \frac{{\sqrt {1 + {{(CC)}^2}} - 1}}{2}}, \quad\; AA < 0. \end{split} \right.$

      综上所述, 考虑到相位匹配, 有耗分层介质中平面波的复数波矢量可通过(10)式计算. 计算得到波矢量虚部${{{k}}_{{\rm{I}}l}}$沿$z$方向, 所以波在介质中沿深度方向衰减. 波矢量实部${{{k}}_{{\rm{R}}l}}$的方向代表波传播方向, 介质中的折射角可以依此计算.

    • 对整个分层介质, 在图1中从入射点${M_0}$到透射点${N_{\rm{t}}}$的传播用前向传播矩阵表示为

      $\left[ {\begin{array}{*{20}{c}} 0 \\ {T{B_0}\exp ( - {\rm{j}}{k_{{\rm{t}}z}}{z_n})} \end{array}} \right] = {{{V}}_{{\rm{t}}0}}\left[ {\begin{array}{*{20}{c}} {R{B_0}\exp ({\rm{j}}{k_{0z}}{z_0})} \\ {{B_0}\exp ( - {\rm{j}}{k_{0z}}{z_0})} \end{array}} \right],$

      其中${B_0}$${M_0}$点处下行波幅值. 整个分层介质总的前向传播矩阵为[22]

      ${{{V}}_{{\rm{t0}}}} = {{{U}}_n} \cdot {{{V}}_n} \cdots {{{U}}_1} \cdot {{{V}}_1} \cdot {{{U}}_0},$

      其中${{{U}}_l}$${{{V}}_l}$分别是跨越$l-l + 1$层界面和第$l$层中的传播矩阵. 对有耗介质, 传播矩阵计算时${\varepsilon _l}, {\mu _l}$取复数, 另外要特别注意${k_{lz}}$, 应该按照前述方法计算.

    • 单层的石墨烯可以看做是“无限薄”的[21], 石墨烯薄层两个侧面的表面电导率为${\sigma _{\rm{g}}}\left( {\omega, {\mu _{\rm c}}, \varGamma, T} \right)$, 其中$\omega $, ${\mu _{\rm c}}$, $\varGamma $$T$分别是角频率、化学势、散射率和室温. 由Kobu公式得到石墨烯表面电导率表达式为${\sigma _{\rm{g}}} = {\sigma _{{\rm{intra}}}} + {\sigma _{{\rm{inter}}}}$, ${\sigma _{{\rm{intra}}}}$取决于能带内的能量, 而${\sigma _{{\rm{inter}}}}$取决于能带间的能量. 因为能带内的能量占能量主要部分, 所以只考虑${\sigma _{{\rm{intra}}}}$这一项, 而忽略${\sigma _{{\rm{inter}}}}$项.

      ${\sigma _{{\rm{intra}}}}$可用Drude表达式表示为

      ${\sigma _{{\rm{intra}}}}(\omega,{\mu _{\rm{c}}},\varGamma,T) = \frac{{{\sigma _0}}}{{1 - {\rm{j}}\omega \tau }},$

      式中

      ${\sigma _0} = \frac{{{e^2}\tau {k_{\rm{B}}}T}}{{{\text{π}}{\hbar ^2}}}\left\{ {\frac{{{\mu _{\rm{c}}}}}{{{k_{\rm{B}}}T}} + 2\ln \left[ {\exp \left( { - \frac{{{\mu _{\rm{c}}}}}{{{k_{\rm{B}}}T}}} \right) + 1} \right]} \right\},$

      为直流电导率, 其中e是电子电量, ${k_{\rm{B}}}$是玻尔兹曼常数, $\hbar $是约化普朗克常量; $\tau = {1 / {\left( {2\varGamma } \right)}}$是电子弛豫时间.

      设在分层介质区域$l$$l + 1$的分界面$z = {z_l}$处有一层“无限薄”石墨烯层, 如图1所示, 根据电磁理论, 在介质分界面两侧的电磁场满足边界条件, 所以在分界面$z = {z_l}$处有

      $\begin{split} & \hat z \times \left( {{{{E}}_l} - {{{E}}_{l + 1}}} \right) = - {{J}}_{\rm{m}}^{\rm{s}}, \\ & \hat z \times \left( {{{{H}}_l} - {{{H}}_{l + 1}}} \right) = {{J}}_{\rm{e}}^{\rm{s}}, \end{split} $

      其中${{J}}_{\rm{e}}^{\rm{s}}$${{J}}_{\rm{m}}^{\rm{s}}$表示表面电流和磁流. 由于该界面为“无限薄”石墨烯层, 则此处[23]

      $\begin{split} & {{J}}_{\rm{m}}^{\rm{s}} = 0, \\ & J_{{\rm{e}},x}^{\rm{s}} = {\sigma _{\rm{g}}}{E_x},~~ J_{{\rm{e}},y}^{\rm{s}} = {\sigma _{\rm{g}}}{E_y}. \end{split} $

      联合(15)式可得

      ${E_{l\alpha }}\left( {z = {z_l}^ + } \right) = {E_{\left( {l + 1} \right)\alpha }}\left( {z = {z_l}^ - } \right),{\rm{ }}\alpha = x,y,$

      ${H_{\left( {l + 1} \right)y}}\left( {z = {z_l}^ - } \right) - {H_{ly}}\left( {z = {z_l}^ + } \right) = {\sigma _{\rm{g}}}{E_x}\left( {z = {z_l}} \right),$

      ${H_{( {l + 1} )x}}( {z = {z_l}^ - } ) - {H_{lx}}( {z = {z_l}^ + } ) = - {\sigma _{\rm{g}}}{E_y}( {z = {z_l}} ).$

      首先考虑TE波情形, 区域l中的总场可写为[22]

      $\begin{split} &{E_{ly}} = [ {{A_l}\exp ( {{\rm{j}}{k_{lz}}z} ) + {B_l}\exp ( { - {\rm{j}}{k_{lz}}z} )} ]\exp ( {{\rm{j}}{k_x}x} ), \\ & {H_{lx}} = \frac{{ - {k_{lz}}}}{{\omega {\mu _l}}}[ {{A_l}\exp ( {{\rm{j}}{k_{lz}}z} ) \!-\! {B_l}\exp ( {\! - {\rm{j}}{k_{lz}}z} )} ]\exp ( {{\rm{j}}{k_x}x} ), \\ & {H_{lz}} = \frac{{{k_x}}}{{\omega {\mu _l}}}[ {{A_l}\exp ( {{\rm{j}}{k_{lz}}z} ) \!+\! {B_l}\exp ( {\! - {\rm{j}}{k_{lz}}z} )} ]\exp ( {{\rm{j}}{k_x}x} ), \end{split} $

      式中${A_l}$${B_l}$分别代表区域$l$中的上行波和下行波幅值. 在区域$l$$l + 1$的分界面$z = {z_l}$处, 根据边界条件(17)和(19)式可得

      $\begin{split} & {A_{\left( {l + 1} \right)}}\exp \left( {{\rm{j}}{k_{\left( {l + 1} \right)z}}z} \right) - {B_{\left( {l + 1} \right)}}\exp \left( { - {\rm{j}}{k_{\left( {l + 1} \right)z}}z} \right) \\ =\, & \left( {p_{\left( {l + 1} \right)l}^{{\rm{TE}}} + s_{\left( {l + 1} \right)}^{{\rm{TE}}}} \right){A_l}\exp \left( {{\rm{j}}{k_{lz}}z} \right) \\ &+ \left( { - p_{\left( {l + 1} \right)l}^{{\rm{TE}}} + s_{\left( {l + 1} \right)}^{{\rm{TE}}}} \right){B_l}\exp \left( { - {\rm{j}}{k_{lz}}z} \right), \\[-16pt] \end{split} $

      $\begin{split} & \left( {1 - s_{\left( {l + 1} \right)}^{{\rm{TE}}}} \right){A_{\left( {l + 1} \right)}}\exp \left( {{\rm{j}}{k_{\left( {l + 1} \right)z}}z} \right) \\ &- \left( {1 + s_{\left( {l + 1} \right)}^{{\rm{TE}}}} \right){B_{\left( {l + 1} \right)}}\exp \left( { - {\rm{j}}{k_{\left( {l + 1} \right)z}}z} \right) \\ = \; &p_{\left( {l + 1} \right)l}^{{\rm{TE}}}\left[ {{A_l}\exp \left( {{\rm{j}}{k_{lz}}z} \right) - {B_l}\exp \left( { - {\rm{j}}{k_{lz}}z} \right)} \right], \end{split} $

      其中

      $s_l^{{\rm{TE}}} = \frac{{\omega {\mu _l}{\sigma _{\rm{g}}}}}{{{k_{lz}}}}.$

      求解并整理(21)和(22)式可得

      $\left[\!\!\!{\begin{array}{*{20}{c}} {{A_{\left( {l + 1} \right)}}\exp \left( {{\rm{j}}{k_{\left( {l + 1} \right)z}}{z_l}} \right)} \\ {{B_{\left( {l + 1} \right)}}\exp \left( { - {\rm{j}}{k_{\left( {l + 1} \right)z}}{z_l}} \right)} \end{array}}\!\!\!\right] \!=\! {{U}}_{\rm{g}}^{{\rm{TE}}} \cdot \left[\!\!\!{\begin{array}{*{20}{c}} {{A_l}\exp \left( {{\rm{j}}{k_{lz}}{z_l}} \right)} \\ {{B_l}\exp \left( { - {\rm{j}}{k_{lz}}{z_l}} \right)} \end{array}}\!\!\!\right],$

      其中

      ${{U}}_{\rm{g}}^{{\rm{TE}}} = \frac{1}{2}\left[\!\!\!{\begin{array}{*{20}{c}} {1 + p_{\left( {l + 1} \right)l}^{{\rm{TE}}} + s_{l + 1}^{{\rm{TE}}}}&{1 - p_{\left( {l + 1} \right)l}^{{\rm{TE}}} + s_{l + 1}^{{\rm{TE}}}} \\ {1 - p_{\left( {l + 1} \right)l}^{{\rm{TE}}} - s_{l + 1}^{{\rm{TE}}}}&{1 + p_{\left( {l + 1} \right)l}^{{\rm{TE}}} - s_{l + 1}^{{\rm{TE}}}} \end{array}}\!\!\! \right]$

      是穿越$z = {z_l}$处石墨烯分界面从${M_l}$点到${N_{l + 1}}$点的TE波前向传播矩阵.

      对TM波情形, 区域$l$中的总场可写为

      $\begin{split} &{H_{ly}} = \left[ {{A_l}\exp \left( {{\rm{j}}{k_{lz}}z} \right) + {B_l}\exp \left( { - {\rm{j}}{k_{lz}}z} \right)} \right]\exp \left( {{\rm{j}}{k_x}x} \right), \\ &{E_{lx}} = \frac{{{k_{lz}}}}{{\omega {\varepsilon _l}}}\left[ {{A_l}\exp \left( {{\rm{j}}{k_{lz}}z} \right) \!- \!{B_l}\exp \left( {\! - {\rm{j}}{k_{lz}}z} \right)} \right]\exp \left( {{\rm{j}}{k_x}x} \right), \\ & {E_{lz}} = \frac{{ - {k_x}}}{{\omega {\varepsilon _l}}}\left[ {{A_l}\exp \left( {{\rm{j}}{k_{lz}}z} \right) \!+\! {B_l}\exp \left( {\! - {\rm{j}}{k_{lz}}z} \right)} \right]\exp \left( {{\rm{j}}{k_x}x} \right). \end{split} $

      在区域$l$$l + 1$的分界面$z = {z_l}$处, 根据边界条件(17)和(18)式可得

      $\begin{split} &{A_{\left( {l + 1} \right)}}\exp \left( {{\rm{j}}{k_{\left( {l + 1} \right)z}}z} \right) + {B_{\left( {l + 1} \right)}}\exp \left( { - {\rm{j}}{k_{\left( {l + 1} \right)z}}z} \right) \\ =\, & \left( {1 + s_l^{{\rm{TM}}}} \right){A_l}\exp \left( {{\rm{j}}{k_{lz}}z} \right) \\ & + \left( {1 - s_l^{{\rm{TM}}}} \right){B_l}\exp \left( { - {\rm{j}}{k_{lz}}z} \right), \\[-16pt] \end{split} $

      $\begin{split} & \left( {1 - s_{\left( {l + 1} \right)}^{{\rm{TM}}}} \right){A_{\left( {l + 1} \right)}}\exp \left( {{\rm{j}}{k_{\left( {l + 1} \right)z}}z} \right) \\ &+ \left( {1 + s_{\left( {l + 1} \right)}^{{\rm{TM}}}} \right){B_{\left( {l + 1} \right)}}\exp \left( { - {\rm{j}}{k_{\left( {l + 1} \right)z}}z} \right) \\=\,& {A_l}\exp \left( {{\rm{j}}{k_{lz}}z} \right) + {B_l}\exp \left( { - {\rm{j}}{k_{lz}}z} \right), \end{split}$

      其中

      $s_l^{{\rm{TM}}} = \frac{{{k_{lz}}{\sigma _{\rm{g}}}}}{{\omega {\varepsilon _l}}}.$

      求解并整理(27)和(28)式可得

      $\left[\!\!\!{\begin{array}{*{20}{c}} {{A_{( {l + 1} )}}\exp ( {{\rm{j}}{k_{( {l + 1} )z}}z} )} \\ {{B_{( {l + 1} )}}\exp ( { - {\rm{j}}{k_{( {l + 1} )z}}z} )} \end{array}} \!\!\!\right] = {{U}}_{\rm{g}}^{{\rm{TM}}} \cdot \left[\!\!\! {\begin{array}{*{20}{c}} {{A_l}\exp ( {{\rm{j}}{k_{lz}}z} )} \\ {{B_l}\exp ( { - {\rm{j}}{k_{lz}}z} )} \end{array}}\!\!\!\right],$

      其中

      ${{U}}_{\rm{g}}^{{\rm{TM}}} = \frac{1}{2}\left[\!\!\!{\begin{array}{*{20}{c}} {1 + p_{\left( {l + 1} \right)l}^{{\rm{TM}}} + s_l^{{\rm{TM}}}}&{1 - p_{\left( {l + 1} \right)l}^{{\rm{TM}}} - s_l^{{\rm{TM}}}} \\ {1 - p_{\left( {l + 1} \right)l}^{{\rm{TM}}} + s_l^{{\rm{TM}}}}&{1 - p_{\left( {l + 1} \right)l}^{{\rm{TM}}} - s_l^{{\rm{TM}}}} \end{array}}\!\!\!\right]$

      是穿越$z = {z_l}$处石墨烯分界面从${M_l}$点到${N_{l + 1}}$点的TM波前向传播矩阵.

      实际上, 上述方法中所涉及的石墨烯界面可拓展到表面电导率为${\sigma _{\rm{g}}}$的导电界面.

      综上, 若分层介质中某界面含“无限薄”导电层, 则可以通过将总传播矩阵${{{V}}_{{\rm{t0}}}}$((12)式)中该界面处的跨界面传播矩阵${{{U}}_{l = 0 - n}}$替换为${{U}}_{\rm{g}}^{{\rm{TE/TM}}}$((25)和(31)式), 从而得到含导电界面的分层介质传播矩阵. 相反, 若令${\sigma _{\rm{g}}} = 0$, 则${{U}}_{\rm{g}}^{{\rm{TE/TM}}}$又可以退化为${{{U}}_{l = 0 - n}}$, 即为一般电介质分界面.

      注意到, 由于石墨烯界面的导电特性, 此处TE和TM波的边界条件不对偶, 因而其跨界面传播矩阵在TE和TM波情形下并不对偶.

    • 设“无限薄”石墨烯上、下半空间分别为均匀介质1和2, 在(24)和(30)式中令$l = 1$${A_2} = 0$, 结合(25)式、(31)式和反透射系数定义可以得到“无限薄”石墨烯界面的反透射系数

      $\begin{split} {R^{{\rm{TE}}}} = & - \frac{{s_1^{{\rm{TE}}} - s_2^{{\rm{TE}}} + s_1^{{\rm{TE}}}s_2^{{\rm{TE}}}}}{{s_1^{{\rm{TE}}} + s_2^{{\rm{TE}}} + s_1^{{\rm{TE}}}s_2^{{\rm{TE}}}}} \\ = & - \frac{{{\mu _1}{k_{2z}} - {\mu _2}{k_{1z}} + {\mu _1}{\mu _2}\omega {\sigma _{\rm{g}}}}}{{{\mu _1}{k_{2z}} + {\mu _2}{k_{1z}} + {\mu _1}{\mu _2}\omega {\sigma _{\rm{g}}}}}, \\ {\rm{ }}{T^{{\rm{TE}}}} =\; & \frac{{s_1^{{\rm{TE}}} + s_2^{{\rm{TE}}} - s_1^{{\rm{TE}}}s_2^{{\rm{TE}}}}}{{2s_1^{{\rm{TE}}}}} \\ &+ {R^{{\rm{TE}}}} \cdot \left( {\frac{{s_1^{{\rm{TE}}} - s_2^{{\rm{TE}}} - s_1^{{\rm{TE}}}s_2^{{\rm{TE}}}}}{{2s_1^{{\rm{TE}}}}}} \right)\\ = \;& \frac{{{\mu _1}{k_{2z}} + {\mu _2}{k_{1z}} - {\mu _1}{\mu _2}\omega {\sigma _{\rm{g}}}}}{{2{\mu _1}{k_{2z}}}} \\ & + {R^{{\rm{TE}}}} \cdot \left( {\frac{{{\mu _1}{k_{2z}} - {\mu _2}{k_{1z}} - {\mu _1}{\mu _2}\omega {\sigma _{\rm{g}}}}}{{2{\mu _1}{k_{2z}}}}} \right), \end{split} $

      $\begin{split} {R^{{\rm{TM}}}} =\, & - \frac{{s_2^{{\rm{TM}}} - s_1^{{\rm{TM}}} - s_1^{{\rm{TM}}}s_2^{{\rm{TM}}}}}{{s_2^{{\rm{TM}}} + s_1^{{\rm{TM}}} + s_1^{{\rm{TM}}}s_2^{{\rm{TM}}}}}\\ =\, & - \frac{{{\varepsilon _1}{k_{2z}} - {\varepsilon _2}{k_{1z}} - {k_{1z}}{k_{2z}}\left( {{{{\sigma _{\rm{g}}}} / \omega }} \right)}}{{{\varepsilon _1}{k_{2z}} + {\varepsilon _2}{k_{1z}} + {k_{1z}}{k_{2z}}\left( {{{{\sigma _{\rm{g}}}} / \omega }} \right)}}, \\ {T^{{\rm{TM}}}} =\, & \frac{{s_2^{{\rm{TM}}} + s_1^{{\rm{TM}}} - s_1^{{\rm{TM}}}s_2^{{\rm{TM}}}}}{{2s_2^{{\rm{TM}}}}} \\ & + {R^{{\rm{TM}}}} \cdot \left( {\frac{{s_2^{{\rm{TM}}} - s_1^{{\rm{TM}}} + s_1^{{\rm{TM}}}s_2^{{\rm{TM}}}}}{{2s_2^{{\rm{TM}}}}}} \right) \\ =\, & \frac{{{\varepsilon _1}{k_{2z}} + {\varepsilon _2}{k_{1z}} - {k_{1z}}{k_{2z}}\left( {{{{\sigma _{\rm{g}}}} / \omega }} \right)}}{{2{\varepsilon _1}{k_{2z}}}} \\ & + {R^{{\rm{TM}}}} \cdot \left( {\frac{{{\varepsilon _1}{k_{2z}}\omega - {\varepsilon _2}{k_{1z}}\omega + {k_{1z}}{k_{2z}}\left( {{{{\sigma _{\rm{g}}}} / \omega }} \right)}}{{2{\varepsilon _1}{k_{2z}}\omega }}} \right). \end{split} $

      当垂直入射时, ${k_{1 z}} = {k_1} = {k_0}\sqrt {{\varepsilon _{{\rm{r1}}}}{\mu _{{\rm{r1}}}}}, {k_{2 z}} = {k_2} = $${k_0}\sqrt {{\varepsilon _{{\rm{r2}}}}{\mu _{{\rm{r2}}}}} $, 由(32)和(33)式可得“无限薄”石墨烯界面的垂直反透射系数

      $\begin{split} & R = \frac{{{\eta _2} - {\eta _1} - {\sigma _{\rm{g}}}{\eta _1}{\eta _2}}}{{{\eta _2} + {\eta _1} + {\sigma _{\rm{g}}}{\eta _1}{\eta _2}}},\\ & {\eta _n} = \sqrt {\frac{{{\mu _n}}}{{{\varepsilon _n}}}}\;\;\left( {n = 1,2} \right), \\ & T = \frac{{2{\eta _2}}}{{{\eta _2} + {\eta _1} + {\sigma _{\rm{g}}}{\eta _1}{\eta _2}}} = 1 + R. \end{split} $

      结果和文献[23]一致. 注意到, 在(32), (33)和(34)式中若令${\sigma _{\rm{g}}} = 0$, 则可以自然过渡为一般介质分界面的Fresnel公式.

    • 例1 石墨烯单层的屏蔽效率

      设各向同性石墨烯层的参数为T = 300 K, $ {\mu _{\rm{c}}} = 0, \varGamma = 0.11$ MeV, 应用传播矩阵计算其反透射系数, 进而得到屏蔽效率(shielding effectiveness, SE)[21]随频率变化曲线如图2(b)所示(图2(a)为该石墨烯层电导率实虚部与频率的关系曲线). SE的计算结果和文献[21]结果相符, 同时由图2可见, 石墨烯层对TM模的SE低于对TE模的.

      图  2  石墨烯电导率和SE (a) 电导率实虚部; (b) SE

      Figure 2.  Graphene conductivity and shielding effectiveness: (a) Real and imaginary parts of conductivity; (b) SE.

      以下算例中石墨烯层的参数取

      $T = 300\;{\rm{ K}},{\mu _{\rm{c}}} = 0.5\;{\rm{ eV}},\tau = 0.5\;{\rm{ ps}}.$

      例2 空气中石墨烯单层的反透射

      设平面电磁波入射真空中一单层石墨烯层, 根据(32)和(33)式可计算得到该石墨烯层的反透射系数, 计算结果如图3所示. 图3(a)图3(b)分别是TE和TM波的反透射系数模值, 其中实线、短划线、点线分别代表入射角为${0^ \circ }$, ${30^ \circ }$, ${60^ \circ }$时的反射系数, 方形、圆形、星形分别代表入射角为${0^ \circ }$, ${30^ \circ }$, ${60^ \circ }$时的透射系数. 可以看到, 随着频率增大, 该石墨烯层的反射系数减小而透射系数增大; 随着入射角增大, TE波反射增强透射减弱, 而TM波相反.

      图  3  石墨烯单层的反射和透射系数模值 (a) TE模; (b) TM模

      Figure 3.  Modulus of reflective and transmittance coefficients of a graphene sheet:(a) TE mode; (b) TM mode.

      例3 涂覆石墨烯层的半空间有耗介质(碲化镉CdTe)界面

      设在碲化镉(CdTe)表面上涂覆着一层石墨烯, CdTe的相对介电系数为[24]

      ${\varepsilon _{\rm{r}}} = {\varepsilon _\infty } + {\varepsilon _\infty }\left( {\frac{{\omega _{{\rm{LO}}}^2 - \omega _{{\rm{TO}}}^2}}{{\omega _{{\rm{TO}}}^2 - {\omega ^2} - {\rm{j}}\gamma \omega }}} \right),$

      其中${\varepsilon _\infty }$, ${\omega _{{\rm{LO}}}}$, ${\omega _{{\rm{TO}}}}$, $\gamma $分别是CdTe的高频极限相对介电常数、纵波光学声子角频率、横波光学声子角频率、碰撞角频率. 本算例取文献[24]给出的测量值

      $\begin{split} &{\varepsilon _\infty } = 6.93,\;\;{\omega _{{\rm{TO}}}} = 2{\text{π}} \times 4.25\;{\rm{ THz}},\\ &{\omega _{{\rm{LO}}}} = 2{\text{π}} \times 5.01\;{\rm{ THz}},\;\;\gamma = 2{\text{π}} \times 0.125\;{\rm{ THz}}.\end{split}$

      计算得到CdTe的禁带为

      ${f_{\rm{L}}} < f < {f_{\rm{H}}},~{f_{\rm{L}}} = 4.2548\;{\rm{ THz}},~ {f_{\rm{H}}} = 5.0044\;{\rm{ THz}}.$

      根据(32)和(33)式可计算该CdTe界面的反透射系数, 计算结果如图4图5所示. 图4图5分别是TE波和TM波情形, 图(a)和(b)分别是反射、透射系数模值, 各图中实线、短划线、点线分别代表入射角为${0^ \circ }$, ${30^ \circ }$, ${60^ \circ }$时不含石墨烯涂层时的情形, 方形、圆形、星形分别代表入射角为${0^ \circ }$, ${30^ \circ }$, ${60^ \circ }$时含石墨烯涂层时的情形. 由图4图5可见, 石墨烯的高导电特性使得其涂覆CdTe表面后, 反射增强而透射减弱, 并且在禁带以下的低频段尤为明显.

      图  4  含石墨烯涂层CdTe半空间的反透射系数模值(TE模) (a) 反射系数; (b) 透射系数

      Figure 4.  Modulus of reflective and transmittance coefficients of CdTe half-space containing graphene coating (TE mode): (a) Reflective coefficient; (b) transmittance coefficient.

      图  5  含石墨烯涂层CdTe半空间的反透射系数模值(TM模) (a) 反射系数; (b) 透射系数

      Figure 5.  Modulus of reflective and transmittance coefficients of CdTe half-space containing graphene coating (TM mode): (a) Reflective coefficient; (b) transmittance coefficient.

      图6图7是该CdTe半空间界面TE和TM入射波的反透射光场分布, 图(a)和(b)分别是界面不含和含有石墨烯涂层时的情形, 从图6图7可见, 石墨烯涂层不改变反透射角, 但增强了表面反射, 起到了一定屏蔽作用.

      图  6  CdTe半空间的反透射光场(TE模) (a) 无石墨烯涂层; (b) 含石墨烯涂层

      Figure 6.  Optical field of reflection and transmission coefficients of CdTe half-space (TE mode): (a) Without graphene coating; (b) with graphene coating.

      图  7  CdTe半空间的反透射光场(TM模) (a) 无石墨烯涂层; (b) 含石墨烯涂层

      Figure 7.  Optical field of reflection and transmission coefficients of CdTe half-space (TM mode): (a) Without graphene coating; (b) with graphene coating.

      例4 涂覆石墨烯层的Si/SiO2周期结构型一维光子晶体

      有一含SiO2缺陷层的Si/SiO2周期结构型1DPC, 如图8所示. 该周期结构周期数$N = 8$, 介质a (Si)和介质b (SiO2)参数为

      图  8  Si/SiO2周期结构型1DPC

      Figure 8.  Si/SiO2 1DPC with periodic structure.

      $\begin{split} & {\varepsilon _{{\rm{ra}}}} = {{\varepsilon }'_{{\rm{ra}}}} + {\rm{j}}{\varepsilon ''_{{\rm{ra}}}},\;\; {{\varepsilon }'_{{\rm{ra}}}} = 12,\;\;{{\varepsilon }''_{{\rm{ra}}}} = 0, \\ & {\varepsilon _{{\rm{rb}}}} = {{\varepsilon }'_{{\rm{rb}}}} + {\rm{j}}{{\varepsilon}''_{{\rm{rb}}}},\;\; {{\varepsilon }'_{{\rm{rb}}}} = 4.5,\;\; {{\varepsilon }''_{{\rm{rb}}}} = 0.07. \end{split} $

      Si, SiO2层和缺陷层的厚度分别为${d_{\rm{a}}} = {{\lambda _0}} / {\sqrt {{{\varepsilon}'_{{\rm{ra}}}}} },$$ {d_{\rm{b}}} = {{{\lambda _0}} / {\sqrt {{{\varepsilon}'_{{\rm{rb}}}}} }}, \;{d_{\rm{s}}} = {{{d_{\rm{b}}}}/ {0.85}}$, 其中${\lambda _0} = $ 20 μm.

      图9是平面波垂直入射时该光子晶体的反透射系数, 其中圆圈表示无石墨烯涂层时的反透射系数, 作为对比用虚线给出了无耗情形(${\varepsilon ''_{{\rm{rb}}}} = 0$)时的反透射系数, 可以看到, 有耗时反透射都有所减弱. 图9中实心和空心星形分别代表石墨烯涂层位于光子晶体表面(第1个界面位置)和底层(最后一个界面位置)时的情形, 可以看到, 在禁带下方低频段, 石墨烯涂层的加入使反射增强, 透射减弱, 在禁带内位于底层的石墨烯涂层对反射系数几乎没有影响而位于表面的石墨烯涂层对反射系数影响较大, 使之减弱.

      图  9  含石墨烯涂层Si/SiO2周期结构1DPC的反透射系数 (a) 反射系数; (b) 透射系数

      Figure 9.  Modulus of reflective and transmittance coefficients of Si/SiO2 1DPC containing graphene sheet: (a) Reflective coefficient; (b) transmittance coefficient.

      图10是平面波垂直入射时该光子晶体的吸收率曲线, 石墨烯涂层增强了光子晶体的吸收, 而且当石墨烯涂层位于表面时, 石墨烯和间隔层构成的表面缺陷导致光的局域化, 使得THz光吸收被增强[18], 在禁带及以上频段, 吸收率显著提高; 当石墨烯涂层位于底层时, 吸收率的提高不明显, 而且主要集中在禁带下方低频段.

      图  10  含石墨烯界面Si/SiO2周期结构1DPC的吸收率

      Figure 10.  Absorbance of Si/SiO2 1DPC containing graphene sheet.

      图11图12分别是TE和TM波情形下, 该光子晶体的吸收率随频率和入射角变化的伪色彩图, 其中子图(a), (b)和(c)分别是无石墨烯涂层、涂层位于表面和底层时的情形, 整体上看, 位于表面的涂层对光吸收的增强效应更为明显, 而且对TE波的吸收要强于TM波.

      图  11  含石墨烯涂层Si/SiO2周期结构1DPC的吸收率(TE模) (a) 无涂层; (b) 表面涂层; (c) 底层涂层

      Figure 11.  Contour plots of the absorbance of the Si/SiO2 1DPC as a function of the light frequency and the incident angles for the TE mode: (a) Without graphene sheet; (b) graphene sheet on the top; (c) graphene sheet on the bottom.

      图  12  含石墨烯涂层Si/SiO2周期结构1DPC的吸收率(TM模) (a) 无涂层; (b) 表面涂层; (c) 底层涂层

      Figure 12.  Contour plots of the absorbance of the Si/SiO2 1DPC as a function of the light frequency and the incident angles for the TE mode: (a) Without graphene sheet; (b) graphene sheet on the top; (c) graphene sheet on the bottom.

    • 根据相位匹配, 讨论给出了平面波斜入射分层介质时, 有耗层中波矢量实部和虚部的计算方法. 有耗层中复数波矢量的实部和虚部不平行, 传播的平面波非均匀. 代表波振幅传播方向的波矢量虚部沿纵向方向, 即波沿深度方向衰减. 波矢量实部代表相位传播方向, 折射角依此计算. 依据边界条件, 推导了跨越石墨烯界面的传播矩阵, 以及“无限薄”石墨烯层的反透射系数解析式. 推广了传播矩阵方法, 使之可应用于包括含导电界面的有耗分层介质反透射和电波传播特性的计算和分析, 为分层介质与导电界面复合结构的设计和验证提供了支撑.

参考文献 (24)

目录

    /

    返回文章
    返回