陈旭生 李九生
(中国计量大学,太赫兹研究所,杭州 310018)
提出一种多缺陷组合嵌入VO2薄膜结构的可调太赫兹吸收器,它由上表面金属图案层、基体和底层金属板三层结构组成,在上表面和基体之间嵌入二氧化钒介质.计算结果表明在 f=4.08 THz和f=4.33 THz两频点吸收率分别为99.8%和99.9%.通过改变外界环境温度可控制二氧化钒相变,从而使两个频点吸收率从99.8%变化到1.0%.改变入射角和偏振态,计算结果表明在入射角0°—40°,吸收器在TE和TM两种极化波下吸收率都能在98%以上.该太赫兹波吸收器具有高吸收、动态调谐、极化不敏感等特性,本文所设计的可调太赫兹吸收器在太赫兹波相关领域,例如探测器、开关、动态调制器、隐身技术等方面具有很好的应用前景.
超材料吸收器具有厚度薄、尺寸小、结构简单和吸收率高的优点,可广泛应用于电磁隐身[1,2]、传感器[3,4]、热成像[5]、太阳能光伏[6,7]、光催化剂[8−10]等领域.Landy等[11]研究出一种超材料吸收器,通过对所设计的吸收器结构参数扫描优化,使得电磁波照射到吸收器中后不发生反射也不透射,实现完全吸收.在众多学者持续努力下,多频点与宽带超材料吸收器都有了很大的发展.Li等[12]利用闭合三圆环在太赫兹波段实现了三频段近完美吸收.Grzeskiewicz等[13]利用微结构金属与SU-8组合在太赫兹波段实现了单频点吸收.Bakshi等[14]利用简单的H形金属镂空结构实现了太赫兹多频带吸收.Mohanty等[15]将简易的PI形和U形金属条组合实现了太赫兹波段多频吸收.Karimi等[16]利用矩形石墨烯阵列实现了太赫兹波段宽带吸收,并利用等效电路模型进行了阐述说明.Biabanifard和Abrishamian[17]将多层条状石墨烯进行组合,不同层石墨烯对应不同费米能级,实现了太赫兹波段宽带吸收.Wang等[18]利用石墨烯与多缺陷金属环相结合在太赫兹波段实现了宽带吸收.Daraei等[19]运用等效电路模型阐述了两层尺寸不同的石墨烯圆盘阵列在太赫兹波段的宽带吸收原理.在固定频点吸收不可调谐的缺点日渐凸显,实现吸收率等可调谐的太赫兹吸收器在太赫兹波相关领域,例如探测器、开关、动态调制器、隐身技术等方面具有很好的应用前景.
本文设计了多缺陷组合嵌入VO2薄膜结构可调太赫兹吸收器,在顶层超材料和基体之间嵌入了一层VO2层,通过控制外界环境温度促使VO2从绝缘态相变成金属态,可实现不改变器件结构对吸收器吸收率的动态调谐.研究结果表明该吸收器在 f=4.08 THz 与 f=4.33 THz 处吸收率分别为99.8%和99.9%,通过改变工作环境温度可实现吸收率从99.9%变化为1%,且在入射角0°—40°对TE和TM两种极化不敏感.
图1所示为本文设计的多缺陷组合嵌入VO2薄膜结构可调太赫兹吸收器,它由上表面金属图案层、二氧化钒层、基体和金属底板组成.顶层金属图案由缺陷圆环、缺陷十字架以及中心圆环三部分构成,顶层金属图案与底层金属板所用材料都为金,电导率为 σgold=4.56×107S/m,厚度为 0.5 μm.二氧化钒层材料为 VO2薄膜,其厚度为 0.2 μm.基体所用材料为 Si,相对介电常数ε=11.9,厚度为 35 μm.图1(b)为吸收器的单元结构俯视图.
为了研究所设计吸收器的双频吸收机理,在室温下,分别对不同组合结构进行了仿真计算,得到三种结构的吸收率曲线(图2).为了分析不同组合结构的吸收机理,对每个结构各自高吸收频率点在TE模式下进行电场分析(图3).组合一结构只有中心圆环和缺陷圆环,低频率点的吸收峰在4.06 THz处吸收率达到93%,此时的能量主要集中在中心圆环的上下部位;高频率点的吸收峰在4.27 THz处吸收率只有73.8%,能量主要集中在缺陷圆环的缺口处和圆环的边缘处.组合二结构只有缺陷圆环和缺陷十字架,低频率点的吸收峰在4.03 THz处吸收率只有76.3%,大部分能量出现在左右两个十字架枝节和带缺陷圆环的间隙处,在缺陷圆环的外边缘也出现了部分能量团;高频率点的吸收峰在4.33 THz处吸收率达到99.5%,与低频率点相比,此处频率点在缺陷处和缺陷圆环边缘处能量都有所加强;此外这种模式还出现了第三个吸收峰,在 4.38 THz 处吸收率为 93.8%.组合三结构为中心圆环和缺陷十字架,低频率点吸收峰在4.06 THz处的吸收率为88.6%,能量主要集中在中心圆环和缺陷十字架的间隙处;高频率点吸收峰在4.39 THz处吸收率达到98.7%,上下两个十字架枝节和中心圆环的外边缘出现了能量团.综上所述可见,低频率点吸收主要受中心圆环影响,其次还受到缺陷圆环影响;高频率点吸收主要受缺陷圆环的影响,受中心圆环的影响较小.将上述三种多缺陷组合成一个整体,该结构可实现在太赫兹波f=4.08 THz和f=4.33 THz 处吸收率达到 99.8%和99.9%.将三种模式的组合进行整合后在TE模式下吸收器电场分布与电流分布如图4,图4(a)为 f=4.08 THz 处电场分布,吸收能量集中在中心圆环上下处,在缺陷圆环边缘也出现了少许能量;图4(b)为 f=4.33 THz 处电场分布,能量主要集中在缺陷圆环的边缘和缺口间隙处;图4(c)为f=4.08 THz处的电流分布,可观察到电流集中在带缺陷圆环、十字架与中心圆环;图4(d)是高频率点f=4.33 THz 处电流分布,电流主要分布在外围缺陷圆环和十字架组合上,中心圆环几乎没有电流分布.可见低频率处的吸收峰受到了中心圆环的影响,同时外围缺陷圆环和十字架的组合对低频率处的吸收率也有影响;高频率点的吸收峰主要受到外围缺陷圆环和十字架的影响.
图1 多缺陷组合结构可调太赫兹吸收器 (a) 单元结构三维立体图;(b) 单元俯视图Fig.1.Adjustable terahertz absorber with multiple defects:(a) Three-dimensional of unit structure;(b) top view of unit structure.
图2 不同组合形成的吸收器吸收曲线Fig.2.Absorption curves of the absorber formed by different combinations.
图3 TE 模式下电场 E 分布图 (a) 组合一,f=4.06 THz;(d) 组合一,f=4.27 THz;(b) 组合二,f=4.03 THz;(e) 组合二,f=4.33 THz;(c) 组合三,f=4.06 THz;(f) 组合三,f=4.39 THzFig.3.Electric field distribution in TE mode:(a) Combination 1,f=4.06 THz;(d) combination 2,f=4.27 THz;(b) combination 2,f=4.03 THz;(e) combination 2,f=4.33 THz;(c) combination 3,f=4.06 THz;(f) combination 3,f=4.39 THz.
对影响吸收频率点的重要参数进行扫描分析,计算结果如图5所示.从图5(a)中看出,随着R增大,低频率点发生蓝移,吸收率也稍微降低.当 R=8.5 μm 时,吸收器的吸收率最高,达到99.2%;当 R=9.0 μm 时,吸收器吸收率为 98.9%;当 R=9.5 μm 时,吸收器的吸收率为 98.9%;当R=10 μm 时,吸收器的吸收率为 98.8%.在这个过程中发现吸收器的高频率点并没有发生明显的移动,吸收率保持不变.可以得出这样的结论,缺陷十字架与中心距离R主要对吸收器的低频率点产生影响,调节R大小,可以实现吸收器的低频率点吸收率的优化.此外,本文还对缺陷圆环缺陷宽度W和中心圆环宽度D进行了仿真计算.从图5(b)中可以看出,缺陷宽度W对吸收器的低频率点几乎没有影响,主要对吸收器的高频率点产生影响.W=4.5 μm,f=4.31 THz 处吸收率为 98.6%;W=5.0 μm,f=4.33 THz 处吸收率为 99.9%;W=5.5 μm,f=4.34 THz 处的吸收率为 99.0%;W=6.0 μm,f=4.35 THz 处吸收率为 97.6%.可见随着W增加,高频率点发生了蓝移.中心圆环宽度D与吸收器吸收率的关系如图5(c)所示.从图5(c)中可以观察到,保持中心圆环圆心不变,改变圆环宽度D从2 μm变化到 4 μm,吸收器在低频率点产生蓝移,而高频率点产生红移,且两个频率点上的吸收率均下降,这是因为中心圆环宽度变化实际上改变了缺陷十字架与圆环之间间隙的宽度,金属之间电磁共振强度也随着距离的增大而减弱,当宽度 D=2.5 μm 时,共振强度最大,吸收效果最好.该吸收器优化得到的基本尺寸参数如下:周期为40 μm,顶层金属图案中缺陷圆环缺口宽度为W=5 μm,缺陷圆环宽度为 2 μm,其外半径为 16 μm,中心圆环宽度为 D=2.5 μm,其外半径为 8 μm,缺陷十字架与单元中心位置距离为R=8.5 μm,缺陷十字架宽度为4 μm.
图4 多缺陷组合TE模式下电场E和电流A分布图(a) f=4.08 THz,电场分布;(b) f=4.33 THz,电场分布 ;(c) f=4.08 THz,电流分布;(d) f=4.33 THz,电流分布Fig.4.Electric field and current distribution in TE mode with multiple defects:(a) f=4.08 THz,electric field distribution;(b) f=4.33 THz,electric field distribution;(c) f=4.08 THz,current distribution;(d) f=4.33 THz,current distribution.
VO2薄膜是一种比较典型的相变材料,VO2由绝缘态向金属态转变时,电导率会有3—5个数量级的变化[20].根据材料介电函数与电导率的关系[21]可以得到相变过程中不同温度所对应的VO2薄膜的电导率关系为
其中,σ 为复合体系的电导率,ε0为真空介电常数,复合体系的介电函数 εC与温度相关.图6所示为VO2电导率随温度变化的曲线,随着温度升高到达相变点 68 ℃,VO2电导率发生骤变;温度继续上升至 80 ℃,VO2电导率成稳定状态,此时相变完成,并且过程可逆.利用前文优化后的结构参数结合器件工作温度变化,计算得到不同温度下吸收器的吸收率曲线如图7(a)所示,当室温 T=25 ℃时,吸收器在 f=4.08 THz和f=4.33 THz 处吸收率分别为99.8%和99.9%.随着温度升高,吸收器的两处频点的吸收率逐渐降低.当 T=40 ℃时,f=4.08 THz 与 f=4.33 THz 处吸收率分别为 79.8%和65.1%,当温度升高到 80 ℃ 时,VO2已相变成金属相,吸收器吸收率变为1.0%.在这渐变的过程中,吸收器吸收峰发生轻微蓝移,这是由于VO2在绝缘态往金属态发生相变过程中电导率在不断变大导致的.根据吸收器理论,吸收率A=1 – R – T,反射率 R=|S11|2,透过率 T=|S21|2,S11和S21分别为吸收器的反射系数和透射系数.其中S11是由吸收器材料与自由空间的阻抗匹配程度决定的,二者之间匹配达到完美状态,也就是电磁波没有反射现象,完全进入到了吸收器中.为了减小S21,通常根据金属底板的趋肤深度,在吸收器底部加一层一定厚度的金属板,这样可以有效地阻止电
磁波从器件中透射出去.图7(b)是利用前文优化后的结构参数并在室温下吸收器归一化表面阻抗后的实部与虚部,其中,归一化阻抗可以表达为
图5 结构参数改变对应的吸收率曲线 (a) 缺陷十字架距离中心距离 R;(b) 缺陷圆环缺陷宽度 W;(c) 中心圆环宽度DFig.5.Absorption curves corresponding to the change of structural parameters:(a) Defect cross distance from center R;(b) defect ring width W;(c) center ring width D.
图6 VO2 的电导率随温度的变化Fig.6.Changes of VO2 conductivity with temperature.
图7 (a) 不同温度下吸收器的吸收率;(b) 吸收器的归一化表面阻抗Fig.7.(a) Absorption of the absorber at different temperatures;(b) normalized surface impedance of the absorber.
在 f=4.08 THz与 f=4.33 THz 处表面阻抗实部与虚部分别为 0.995 与 0.048,1.05 与–0.016.这表明所设计吸收器的归一化表面阻抗与自由空间的阻抗形成良好匹配,吸收器达到了完美吸收.
在设计整个吸收器时将偏振敏感性也列入了考虑范围,所以吸收器为中心对称结构.为了验证这一效果,对TE和TM两种模式下太赫兹波不同入射角照射吸收器的吸收光谱进行仿真分析,结果图8 所示.由图8 可见,入射角在 0°—40°吸收器的吸收率可达到98%以上,当太赫兹波入射角为55°时,吸收率下降至90%以下.所设计的吸收器具有偏振不敏感特性,该吸收器在隐身材料、传感等领域都有很好的应用价值.
图8 太赫兹吸收谱 (a) TE 模式;(b) TM 模式Fig.8.Terahertz absorption spectrum:(a) TE mode;(b) TM mode.
提出一种多缺陷组合嵌入VO2薄膜结构可调太赫兹吸收器,分析了吸收器的吸收机理并优化了结构参数,在 0°—40°入射角范围内 f=4.08 THz和f=4.33 THz频率处的吸收率达到99.8%以上.通过改变外界环境温度可以调控吸收器的吸收效率.计算结果表明该吸收器对TE和TM两种极化波不敏感.所设计的可调太赫兹吸收器在太赫兹波相关领域,例如探测器、开关、动态调制器、隐身技术等方面具有很好的应用前景.