吴柔兰 李九生

(中国计量大学太赫兹技术与应用研究所,杭州 310018)

利用VO2 嵌入超表面设计了一种实现不同频率,且线极化和圆极化两种模式入射下均产生高效率吸收的太赫兹超表面.当VO2 为绝缘态时,设计的超表面对圆极化波的旋向产生选择性吸收,在1.30 THz 处对左旋圆极化波产生的吸收率大于95%,对右旋圆极化波不吸收,圆二色性为0.85.当VO2 为金属态时,在1.95 THz处,该超表面对TE 线极化入射波吸收率达到98.5%.结果表明,在线极化和圆极化波入射下,所设计的超表面结构具有良好的广角吸收性能.由于它具有形态简单、易于加工等特点,在太赫兹波传感、成像和通信领域具有广阔的应用前景.

1 引言

近年来超快激光技术的发展为太赫兹波脉冲的产生提供了稳定、可靠的激光发源,使太赫兹波研究获得长足发展[1].超表面作为一种新型二维人工超材料,可以有效控制电磁波的透射、反射与极化等参量,因此备受各国研究人员关注.另外,利用VO2[2,3]、石墨烯[4,5]和锑化铟[6,7]等构建复合超表面结构可以实现对微波、太赫兹波和光波的动态调控[8-10].最近,可调太赫兹吸收器作为太赫兹波调控器件之一,引起研究人员广泛关注[11,12].2018 年,Hu 等[13]利用石墨烯构建不对称分裂环图案层获得多重共振,实现太赫兹波双宽带吸收器.Yang 等[14]利用8 个尺寸渐变的单元组成光栅结构实现对TM 和TE 两种偏振波在不同频率的高效吸收.2019 年,Yan 等[15]采用不同尺寸石墨烯方环设计了可调谐的双频点太赫兹吸收超表面.2020 年,Divdel 等[16]利用希尔伯特分形结构设计出拥有多个吸收波段的太赫兹超表面吸收器,吸收率大于90%.上述报道的超表面主要针对线极化太赫兹波入射实现吸收,没有对圆极化入射到该类型超表面结构产生何种效果进行研究分析.最近,Zhao 等[17]提出基于阿基米德螺旋结构设计超表面实现对左旋圆极化(left-handed circularly polarized,LCP)波吸收的同时也对右旋圆极化(right-handed circularly polarized,RCP)波产生完全反射.2022 年,Liang 等[18]利用VO2和石墨烯构建复合超表面结构通过改变VO2相变状态,实现对圆极化波的非选择性完美吸收与选择性完美吸收功能可切换.VO2是一种典型的温控相变材料,其电导率在相变过程中会发生巨大突变.当低于临界温度(68 ℃)时,VO2具有高电阻率和良好的绝缘性能,而高于临界温度时,电阻由高阻态变为低阻态且表现出金属特性.因此,将VO2嵌入超表面中组成复合超表面结构,实现工作频点可切换,非常有必要研究VO2复合型超表面的太赫兹特性.

本文提出一种VO2复合型超表面实现对线极化和圆极化太赫兹波高效吸收.当VO2绝缘态(σ=0 S/m)时,在频率1.3 THz 处,所设计超表面对LCP 波吸收率达到95%,对RCP 波产生完全反射,圆二色性(CD)为0.85.当VO2金属态(σ=2×105S/m)时,在频率1.95 THz 处,所设计超表面对TE 线极化入射波的吸收率达到98.5%.研究发现该结构对线极化和圆极化入射波均能实现广角吸收,有望为太赫兹的新型探测器和传感器设计提供多样性的选择.

2 结构设计与理论分析

本文提出的太赫兹吸收器结构,如图1 所示,其工作频率可切换,对线极化和圆极化入射的太赫兹波均产生完美吸收.该复合超表面由嵌入VO2的金属图案层,石英基体和金属底板组成.利用CST Microwave Studio 软件优化后复合超表面结构尺寸: 周期为60 µm,w1=4 µm,w2=1 µm,方环长度为44 µm,圆环外半径为15 µm,θ=70°,石英基体厚度为18.5 µm.VO2的相对介电常数利用Drude 模型可表示为[19]

图1 所设计的工作频率可切换,对线极化和圆极化入射波均产生完美吸收的太赫兹吸收器结构示意图 (a) 三维结构图;(b) 单元俯视图Fig.1.Schematic diagram of the designed terahertz absorber with switchable operating frequency and perfect absorption for both linearly polarized and circularly polarized waves incidence: (a) 3D structure diagram;(b) top view of the unit.

式中,ε∞=12,碰撞频率γ=5.75×1013rad/s,ωp(σ)为等离子体频率.

式中,r++和r—+分别为右圆极化波入射下的共极化和交叉极化反射系数,r——和r+—为左圆极化波入射下的共极化和交叉极化反射系数.Rcirc为圆极化的反射矩阵,该矩阵由线极化反射系数表示为

式中,rxx(ryy)和rxy(ryx)分别为线性波入射下的共极化和交叉极化反射系数.若所设计结构对LCP入射波完全吸收并对RCP 入射波完全反射,此时反射系数满足条件为r——=r——=r+—=0,r—+=1,则线极化反射系数可表示为

式中,α为任意相.若所设计的复合超表面结构要实现旋向吸波,则该结构需满足

式中,φ代表单元结构旋转一圈后出现与原始结构重合的角度.φ=mπ (m=0,±1,···),说明只有二重对称结构才能实现对旋向入射波吸波.除旋转对称结构外,镜像对称结构若想实现旋向吸收性能,需满足如下条件:

由(4)式可知,若(rxx—ryy)/rxy为虚数,将无法满足(6)式中的镜像对称条件.因此,实现旋向选择吸收结构的同时必须打破旋转对称与镜像对称.

3 仿真结果分析

3.1 VO2 为绝缘态

图2 表示VO2为绝缘态(σ=0 S/m)时,LCP波和RCP 波入射到复合超表面结构时产生的电磁响应下曲线.从图2(a)可以看出,LCP 波的共极化反射曲线(r——)与RCP 波的共极化反射曲线(r++)一致,LCP 波的交叉极化反射(r+—)与RCP 波的交叉极化反射曲线(r—+)有较大差别.在频率1.30 THz 处,r—+＞ 0.9,而r+—＜ 0.1,结果表明在该频点处所设计超表面反射了大部分RCP波,对LCP 波不反射.为了更加清楚地解释不同圆极化波之间的吸收差,引入CD 进行描述,定义圆极化吸收率为

则CD 表示为

从图2(b)可以明显观察到,在频率1.30 THz处所设计的超表面结构对LCP 入射波吸收率达到95%,但是对RCP 入射波吸收率仅为10%,它的CD 值为0.85,表明几乎所有入射LCP 太赫兹波被吸收,同时大部分入射RCP 太赫兹波被反射.图3 是在LCP 波和RCP 波入射下,顶层金属结构在1.30 THz 处的电场和电流分布图.由图3(a)和图3(b)可以看出,当RCP 太赫兹波入射时,在方形缺口环和圆形缺口环的一侧存在单极共振,但是当LCP 太赫兹波入射时,圆形缺口环产生明显的偶极共振.同样地,从图3(c)和图3(d)可以看出,RCP 太赫兹波入射时圆环表面电流非常微弱,而LCP 太赫兹波入射时表面电流得到显著增强.

图2 LCP 波和RCP 波入射下所设计超表面结构的电磁响应曲线 (a)反射系数;(b)吸收率及圆二色性Fig.2.Electromagnetic response curves of the designed metasurface structures under the incident of LCP and RCP waves: (a) Reflection coefficient;(b) absorption and circular dichroism.

图3 LCP 和RCP 太赫兹波入射下的电场分布图与电流分布图 (a) RCP 入射下的电场分布;(b) LCP 入射下的电场分布;(c) RCP 入射下的电流分布;(d) LCP 入射下的电流分布Fig.3.Electric field distribution diagram and current distribution diagram under the LCP and RCP waves incidence: (a) Electric field distribution under RCP incidence;(b) electric field distribution under LCP incidence;(c) current distribution under RCP incidence;(d) current distribution under LCP incidence.

选取超表面结构参数线宽w1,w2和间隙θ研究其对所设计超表面器件CD 值的影响情况,结果如图4 所示.图4(a)表示w1从2 µm 增加至4 µm时,所设计超表面结构的CD 只有微小变化,这是由于LCP 入射超表面时方形缺口环的电场能量较弱,对CD 值影响较小.图4(b)表示随着w2增加,所设计超表面的CD 峰值频率出现蓝移,且峰值迅速减小,在w2=5 µm 处CD 值接近于0.这是由于所设计的超表面结构谐振频率与缺口环等效的电容和电感成反比,增加线宽会减小表面电流密度和等效电感,从而使峰值频率发生蓝移.图4(c)表示所设计超表面的CD 峰值对应的频率同样随着θ增加发生明显蓝移,这是因为圆环缺口可等效成电容,当缺口角度不断增加时等效电容不断减小,使得峰值频率向更高频方向移动.

图4 不同结构参数对所设计超表面的CD 影响情况 (a) 线宽w1;(b) 线宽w2;(c) 间隙θFig.4.Influence of different structural parameters on the CD of the designed metasurface: (a) Line width w1;(b) line width w2;(c) angle θ.

图5 讨论了太赫兹波斜入射角度对所设计超表面结构的圆极化吸收性能影响.当太赫兹波以入射角70°入射时,在频率1.30 THz 处入射的LCP太赫兹波吸收率大于80%,如图5(a)所示,结果表明所设计的超表面具有良好的广角圆极化吸收性能.图5(b)表示当入射角为70°时,在频率1.30 THz处,所设计超表面的CD 大于0.7,表明该超表面在大入射角情况下依旧能够保持着良好的旋向吸收性能.

图5 入射角对所设计超表面的吸收特性影响 (a) LCP;(b) CDFig.5.Incident angle vs.absorption properties of the designed metasurface: (a) LCP;(b) CD.

研究分析在频率0.8—1.6 THz 内,线性太赫兹波(TE 极化波与TM 极化波)入射到所设计超表面结构产生的吸收特性如图6 所示.从图6(a)可看出,TE 波入射下产生的共极化反射曲线(ryy)与TM 波入射下产生的共极化反射曲线(rxx)变化趋势一致,TE 波入射下产生的交叉极化反射(rxy)与TM 波入射下产生的交叉极化反射(ryx)重叠.在频率1.30 THz 处,rxx=0.35,ryy＞ 0.53.如图6(b)所示,所设计超表面结构在频率1.30 THz处对TE 和TM 入射波的吸收率分别为43.86%和58.40%.

图6 TE 波和TM 波入射下所设计结构的电磁响应曲线 (a) 反射系数;(b) 吸收率Fig.6.Electromagnetic response curves of the designed structures under the incident of TE and TM waves: (a) Reflection coefficient;(b) absorption.

3.2 VO2 为金属态

当VO2为金属态(σ=2×105S/m)时,LCP和RCP 太赫兹波入射到所设计的超表面结构时产生的电磁响应曲线见图7.图7(a)表示在LCP和RCP 太赫兹波入射下所设计超表面结构产生的交叉极化反射和共极化反射曲线.可见,在频段1.6—2.4 THz 内,两种圆极化波对应的交叉极化和共极化反射曲线完全一致.图7(b)为LCP 与RCP太赫兹波入射到所设计超表面结构获得的电磁波吸收率,在频率1.95 THz 处吸收率大于90%,圆二色性CD 值为0,此时所设计超表面结构表现为对不同偏振态的太赫兹波没有选择性吸收.

图7 LCP 和RCP 太赫兹波入射到所设计超表面结构产生的电磁响应曲线 (a) 反射系数;(b) 吸收率及圆二色性Fig.7.Electromagnetic response curves of the designed structures under LCP and RCP waves incidence: (a) Reflection coefficient;(b) absorption and circular dichroism.

图8 表示TE 和TM 波垂直入射所设计的超表面结构时产生的太赫兹波反射和吸收谱.图8(a)中TE 与TM 极化波入射所设计的超表面结构的交叉反射系数rxy=ryx=0,且在1.95 THz 处,TM极化波的共极化反射系数rxx＞ 0.42,TE 极化波的共极化系数ryy=0.12.图8(b)表示TE 和TM 波入射到所设计的超表面结构产生的吸收曲线,在1.95 THz 处,吸收率分别为98.5%和82.0%.

图8 TE 和TM 极化太赫兹波入射到复合超表面结构产生的电磁响应曲线 (a) 反射谱;(b) 吸收谱Fig.8.Electromagnetic response curves of TE and TM polarized terahertz wave is incident on the designed composite metasurface structure: (a) Reflection coefficient;(b) absorption coefficient.

所设计的复合超表面结构对太赫兹波吸收率为A=1—RT=1—|S11|2—|S21|2,式中R和T分别为反射率和透射率,S11与S21是超表面太赫兹波反射参数与透射系数.由于该复合超表面结构的底层金属板厚度为1 µm,这远大于太赫兹波的趋肤深度,因此太赫兹波的透射率为0.引入阻抗匹配理论[21]分析该复合超表面结构对太赫兹波的吸收特性:

式中,Z1表示复合超表面结构等效阻抗,Z0=120π Ω 表示自由空间阻抗.当等效阻抗与自由空间阻抗匹配时,Z≈ 1,此时该结构可以实现完美吸收.图9 表示TE 极化波入射到复合超表面时产生的吸收谱与等效阻抗特性.在频率1.95 THz 处,复合超表面的等效阻抗实部为1.29,趋近于1,表明复合超表面在该频率可以实现太赫兹波的完美吸收.由图10 可见,频率1.95 THz 处TE 极化波入射到所设计的复合超表面结构时产生的电场能量主要集中在圆环和方环上下侧间隙处,此频率处的太赫兹波吸收峰是由方环和圆环的偶极共振效应叠加产生.

图9 TE 波入射到复合超表面结构产生的等效阻抗实部与虚部Fig.9.Real and imaginary parts of the equivalent impedance of the designed composite metasurface structure under the TE wave incidence.

图10 TE 波入射下,复合超表面结构电场分布Fig.10.Electric field distribution of the designed composite metasurface structure under TE wave incidence.

图11 为线极化波入射下超表面复合结构的不同结构参数对太赫兹吸收性能影响情况.图11(a)显示当w1从2 µm 增加至6 µm 时,复合超表面结构的太赫兹吸收峰出现轻微蓝移,吸收率仍大于90%.图11(b)表示随着w2增加,复合超表面结构的太赫兹吸收峰出现蓝移,吸收率逐渐降低.由图11(c)可见,复合超表面结构的吸收峰随着R增加产生红移现象.研究表明当太赫兹波入射角从0°变化到70°时,该复合超表面结构的太赫兹吸收率均大于90%,具有完美的广角吸收性能.如图12所示,当入射角大于30°时,在更高频处产生了另外一个吸收峰,并且随着入射角度继续增加,该谐振频率也发生轻微红移,这主要由超表面结构产生高阶模式共振所致.

图11 不同结构参数对线极化吸收影响 (a) w1;(b) w2;(c) RFig.11.Influence of different structural parameters on linear polarization absorption: (a) w1;(b) w2;(c) R.

4 结论

本文设计了一种VO2嵌入式复合超表面结构,可以实现对多个频率点的线极化和圆极化入射波完美吸收.当VO2绝缘态时,该超表面能够在1.30 THz 处实现对入射LCP 太赫兹波的吸收,同时对入射的RCP 太赫兹波产生反射.当VO2为金属态时,线极化波入射到超表面,在频率1.95 THz处实现太赫兹波的完美吸收.该结构对TE 极化和LCP 太赫兹波均具有良好的广角吸收性能.该复合超表面结构能实现不同频率、不同偏振态的太赫兹波产生良好吸收效应,所以该复合超表面结构的设计构思可用于其他超表面太赫兹器件设计,而且由于对不同偏振信号的响应特征不同,可以用于太赫兹成像与传感系统.