王飞 魏兵

1) (西安电子科技大学物理与光电工程学院, 西安 710071)

2) (西安电子科技大学信息感知技术协同创新中心, 西安 710071)

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

1 引 言

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]是一种分析分层介质中光或电磁波传播特性的解析理论方法, 可以直接计算分层介质对平面电磁波的反透射等.根据电磁理论, 当分层介质中某层为有耗介质时, 该层中的波矢量为复数, 而在斜入射情形下, 考虑相位匹配, 该层中波矢量的实部和虚部并不平行, 即相位和振幅的传播方向不同, 此时该层中传播的应为非均匀平面波.本文首先根据相位匹配条件, 推导了有耗分层介质中复数波矢量的实部和虚部, 给出传播矩阵中所需要的纵向波矢量计算式; 然后考虑石墨烯表面电导率, 根据边界条件推导了跨越石墨烯分界面的传播矩阵; 最终, 使得传播矩阵方法可用于解析分析任意入射角情形下含导电界面有耗分层介质的反透射和电波传播特性.本文中时谐因子取 e xp(-jωt) .

2 理论公式

2.1 有耗分层介质中的复数波矢量

设平面波从真空入射到n层分层介质, 如图1所示, 入射角为θ.考虑一般情形, 设各层介质的相对介电系数和磁导系数εrl,µrl(l=1,2,···,n) 皆为复数,

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

图1 分层介质Fig.1.Stratified medium.

各分层中的复数波矢量kl可以写为

其中kRl和kIl分别表示kl的实部和虚部.如图1 所示, 入射面为xoz平面, 故有

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

即kl的虚部kIl只有z分量, 实部kRl的x分量在每个区域都相同, 都等于区域0中入射波矢量k0的x分量k0x.由此可以看到, 平面波斜入射有耗分层介质时, 介质中波矢量的实部和虚部方向不平行,虚部只有纵向分量.由于波矢量实部代表相位传播方向, 虚部代表振幅传播方向, 可知此时平面波是非均匀的, 波矢量为

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

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

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

求解(8)式得到

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

2.2 有耗分层介质的传播矩阵

对整个分层介质, 在图1中从入射点M0到透射点Nt的传播用前向传播矩阵表示为

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

其中Ul和Vl分别是跨越l-l+1 层界面和第l层中的传播矩阵.对有耗介质, 传播矩阵计算时εl,µl取复数, 另外要特别注意klz, 应该按照前述方法计算.

2.3 跨越石墨烯界面的传播矩阵

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

σintra可用Drude表达式表示为

式中

为直流电导率, 其中e是电子电量,kB是玻尔兹曼常数, ℏ 是约化普朗克常量;τ=1/(2Γ) 是电子弛豫时间.

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

联合(15)式可得

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

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

其中

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

其中

是穿越z=zl处石墨烯分界面从Ml点到Nl+1点的TE波前向传播矩阵.

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

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

其中

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

其中

是穿越z=zl处石墨烯分界面从Ml点到Nl+1点的TM波前向传播矩阵.

实际上, 上述方法中所涉及的石墨烯界面可拓展到表面电导率为σg的导电界面.

综上, 若分层介质中某界面含“无限薄”导电层, 则可以通过将总传播矩阵Vt0((12)式)中该界面处的跨界面传播矩阵Ul=0-n替换为((25)和(31)式), 从而得到含导电界面的分层介质传播矩阵.相反, 若令σg=0 , 则又可以退化为Ul=0-n, 即为一般电介质分界面.

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

2.4 “无限薄”石墨烯层的反透射系数

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

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

3 算 例

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

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

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

图2 石墨烯电导率和 SE (a) 电导率实虚部; (b) SEFig.2.Graphene conductivity and shielding effectiveness：(a) Real and imaginary parts of conductivity; (b) SE.

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

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

图3 石墨烯单层的反射和透射系数模值 (a) TE 模;(b) TM 模Fig.3.Modulus of reflective and transmittance coefficients of a graphene sheet：(a) TE mode; (b) TM mode.

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

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

其中ε∞,ωLO,ωTO,γ分别是 CdTe 的高频极限相对介电常数、纵波光学声子角频率、横波光学声子角频率、碰撞角频率.本算例取文献[24]给出的测量值

计算得到CdTe的禁带为

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

图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.

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

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

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

Si, SiO2层和缺陷层的厚度分别为

图8 Si/SiO2 周期结构型 1DPCFig.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.

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

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

图10 含石墨烯界面 Si/SiO2周期结构 1DPC 的吸收率Fig.10.Absorbance of Si/SiO2 1DPC containing graphene sheet.

图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.

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

4 结 论

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