1. 引言
近些年来,超材料因其具有负介电常数、负磁导率等独特的电磁响应而引起了广泛的关注,应用于隐形[1] [2] [3] [4] [5]、传感、[6]-[13]完美吸收器[14] [15] [16] [17]和能量收集等领域。超材料在结构上层或里面设计周期性排列的亚波长金属谐振器,亚波长谐振单元可以调控其电磁特性。超材料是具有可变介电常数ε和磁导率μ的人工周期性结构,可以达到
和
[18],有效磁导率
和有效介电常数
是独立指定,通过麦克斯韦方程组可以得到超材料的折射率为:
,而超材料结构的表面整阻抗为:
。
在2008年,Landy等人[19]首先提出并设计了一款金属–介质–金属的三层结构的太赫兹超材料吸收器,在1.125 THz处有一个谐振吸收峰,吸收率达到了98%。在2022年,Wang等人[20]提出了由两条水平金属条和两条垂直连接的金属条组成的双波段极化可控太赫兹超材料吸收体。为了实现太赫兹超材料吸收器的“完美吸收”,就是使吸收率最大化即透射率和反射率最小化,满足吸收率为1,透射率和反射率为零。为了使超材料吸收器透射率为零,只需要在结构底部加上一层金属反射层即可。入射电磁波在金属中传播具有趋肤深度,所以只要金属反射层的厚度远远大于趋肤深度,使得原本应该透射而出的电磁波被金属反射层所阻挡,因此便实现了透射率为零。为了使超材料吸收器反射率为零,需要超
材料结构的表面阻抗和空气阻抗相匹配,即超材料结构的阻抗Z = 1。根据公式
可知,当有效介电常数与磁导率相等时,超材料结构的阻抗Z = 1,即实现了超材料结构的表面阻抗和空气阻抗相匹配。设A为吸收率,R为反射率,T为透射率,则关系式为:
。若R为反射率
,T为透射率
,那么会得到以下关系式:
[21]。
本文设计了一款“川”字型的金属–介质–金属的“三明治”式的太赫兹超材料吸收器,由于该结构底部为金属反射层,,所以,透射率T即
为0。因此,
,即
,由此便可得到该太赫兹超材料吸收器的吸波特性。利用有限元仿真软件进行建模与仿真优化。该太赫兹超材料吸收器在1~3 THz之间有一个在2.84 THz谐振吸收峰,在此频率处实现了对太赫兹波99.93%的“完美吸收”,这是因为所设计的结构能够与入射的电磁波发生耦合现象[22]。为了进一步扩展太赫兹超材料吸收器的功能,在其金属条带表面随机生长一定尺寸范围内的金纳米颗粒,仿真结果表明,一定尺寸范围内的金纳米颗粒的存在,会引起太赫兹超材料吸收器在太赫兹波段谐振吸收峰的偏移,但对峰强度的响应影响较小。此外,选取金属条带上的一颗金纳米颗粒作为可见光波段仿真的结构,验证了该结构在可见光波段也有响应。所以,该超材料多级结构在太赫兹和可见光波段均有响应,这对于太赫兹超材料吸收器应用广度的拓展有着一定的影响。
2. 太赫兹超材料吸收器的建模与仿真
2.1. 太赫兹超材料吸收器的设计
太赫兹超材料的单元结构示意图1所示,(a)为单元结构,(b)为正视图,典型的金属–介质–金属结构。顶层是由简单金属条带组成,金属材料为金,金属条带宽度为0.6 μm,金属条带间距为5 μm,沿着X轴方向“川”型排列。中间介质层为FR-4,介质层长度为50 μm,厚度为20 μm。底层是全覆盖的金膜,厚度为0.1 μm。利用基于有限元的电磁仿真软件,分别对图1所示的太赫兹超材料进行仿真分析。为了模拟无限周期阵列,X和Y方向设置为unit cell边界,Z方向上设置开放性的边界条件。
Figure 1. Schematic diagram of the unit structure of terahertz metamaterial absorber, (a) is a three-dimensional view, and (b) is an elevation view
图1. 太赫兹超材料吸收器单元结构示意图,(a) 为立体图,(b) 为正视图
设定太赫兹波垂直入射到太赫兹超材料结构的表面,电场方向沿着y轴正方向,磁场方向沿着x轴正方向,由此得到太赫兹超材料结构的频率吸收谱,如图2所示。由图2可知,该太赫兹超材料在2.84 Thz处有一个吸收率为99.93%的谐振吸收峰。
Figure 2. Frequency absorber spectrum of terahertz metamaterial absorber
图2. 太赫兹超材料吸收器的频率吸收器谱
2.2. 太赫兹超材料吸收器的优化
采用控制变量法对金属条带间距、金属条带宽度、介质层长度、介质层厚度进行参数化扫描,得到如图3所示的频率吸收谱。太赫兹波垂直入射在太赫兹超材料上,如图3所示。因为该太赫兹超材料吸收器的结构由下到上分别为金基底、FR-4以及带有空气槽的金,是典型的“三明治”结构,所以形成了一个法布里-珀罗(Fabry- Perot, F-P)腔,腔吸收效应是由于腔内多次反射之间的干扰而产生的[23] [24] [25],这导致了显著的场捕获和增强。所以,腔长(FR-4的厚度)也会对该太赫兹超材料吸收器的透射曲线产生很大的影响。
对参数化仿真结果进行逐一对比,选定金属条带宽度为0.6 μm、金属条带间距为5 μm、介质层长度为50 μm、厚度为20 μm为最优结构。
3. 太赫兹超材料多级结构的建模与仿真
3.1. 太赫兹超材料多级结构的设计
基于金纳米颗粒的太赫兹超材料多级结构如图4所示,(a)为单元结构立体图,(b)为正视图,顶层金属条带上随机分布着一定尺寸的金纳米颗粒。
利用基于有限元的电磁仿真软件图4所示的超材料多级结构进行仿真分析。为了模拟无限周期阵列,X和Y方向设置为unit cell边界,Z方向上设置开放性的边界条件。
设定太赫兹波垂直入射到太赫兹超材料多级结构的表面,电场方向沿着y轴正方向,磁场方向沿着x轴正方向,对金纳米颗粒的半径R (30 nm、50 nm、60 nm、70 nm)进行频域扫描,分别得到太赫兹超材料多级结构的频率吸收谱。如图5所示,金纳米颗粒的存在会引起太赫兹超材料吸收器谐振峰微弱的偏移,但对强度基本没有影响。若金纳米颗粒尺寸更小,那么引起的吸收峰位置偏移量会更小,对太赫兹超材料吸收器的影响微乎其微。根据超材料吸收器结构设计的原理,如果金纳米颗粒较大,随机生长在超材料金属谐振层边缘位置,可能会造成超材料吸收器谐振吸收峰的位置及强度较大的变化。
Figure 3. Frequency absorber spectrum of terahertz metamaterial absorber. (a) Parameterization of metal strip width; (b) Parameterized scanning of the spacing between metal strips; (c) Parameterized scanning of the length of the dielectric layer; (d) Parameterized scanning of the thickness of the intermediate dielectric layer
图3. 太赫兹超材料吸收器的频率吸收器谱。(a) 对金属条带宽度进行参数化;(b) 对金属条带间距进行参数化扫描;(c) 对介质层长度进行参数化扫描;(d) 对中间介质层厚度进行参数化扫描
Figure 4. Shows a schematic diagram of the multi-level structure of metamaterials. (a) is a three-dimensional view of the unit structure, and (b) is a front view
图4. 超材料多级结构示意图。(a) 为单元结构立体图,(b) 为正视图
Figure 5. The effect of gold nanoparticle size on the frequency absorption of terahertz metamaterials with multi-level structures
图5. 金纳米颗粒尺寸大小对太赫兹超材料多级结构频率吸收的影响
3.2. 太赫兹超材料多级结构在可见光波段的响应特性
金纳米颗粒是具有纳米尺寸(通常指直径在1~100 nm之间)的金材料粒子。纳米金颗粒具有局部表面等离子体共振、大比表面积、较强的生物亲和力和优良的化学特性,使得其在很多领域具有广泛的应用和研究价值。拉曼增强效应是指金纳米颗粒表面等离激元共振引起的拉曼散射信号的增强。当激光光束与金纳米颗粒相互作用时,金纳米颗粒表面的电子受到激发,产生一种局域表面等离激元(LSP),这种激元的共振频率与激光波长相匹配。当分子与金纳米颗粒接近时,分子的振动模式会受到这种局域场的显著增强,从而产生更强的拉曼散射信号。这种效应可用于增强拉曼信号的灵敏度,从而实现对样品中微量成分的检测和表征。
在超材料多级结构上选取部分来研究研究金膜上金纳米颗粒在可见光波段的响应特性。选取的结构的正视图和侧视图,如图6所示。选定金纳米颗粒的半径R为30 nm、50 nm、60 nm、70 nm,可见光波段范围为400~700 nm。
Figure 6. Shows the multi-level structure of selected terahertz metamaterials, where (a) is the front view and (b) is the stereo view
图6. 选取的部分太赫兹超材料多级结构,(a) 为正视图,(b) 为立体图
首先,选定入射光波长为532 nm,可见光垂直入射在该结构上,金纳米颗粒的半径R为30 nm、50 nm、60 nm、70 nm,利用仿真软件对金纳米颗粒的半径频域扫描,分别得到相对应的XY平面电场图如图7所示。由图7可知,电场强度最强的位置在金纳米颗粒与金膜的间隙,也就是所谓的“热点”。在实际检测时,该超材料多级结构产生的拉曼增强效果会由于金纳米颗粒之间的相互作用而更好。所以,该超材料多级结构在可见光波段的响应可用于物质的检测与分析。
Figure 7. The effect of radius R variation of gold nanoparticles on the electric field strength of the selected structure. (a) When R is 30 nm; (b) When R is 30 nm; (c) When R is 30 nm; (d) When R is 60 nm
图7. 金纳米颗粒的半径R变化对选取结构的电场强度的影响。(a) 当R为30 nm时;(b) 当R为30 nm时;(c) 当R为30 nm时;(d) 当R为60 nm时
选定金纳米颗粒的半径R为30 nm,利用仿真软件对入射光波长在400~700 nm进行参数化扫描,步长为50 nm,得到XY平面电场图如图8所示。
从图8可以清楚的看出,此结构在几乎整个可见光波段皆有响应。所以,一定尺寸的金纳米颗粒结合柔性太赫兹超材料而成的超材料多级结构在太赫兹波段和可见光波段皆有响应,这种将大尺寸结构的太赫兹超材料吸收器与小尺寸的金纳米颗粒相结合的超材料多级结构,为拓展超材料吸收器的应用广度具有重要的意义。
4. 总结
本文设计了设一款“川”字型的金属–介质–金属的“三明治”式的太赫兹超材料吸收器,利用有
Figure 8. The influence of visible light wavelength variation on the electric field intensity of the selected structure. (a) When the visible light wavelength is 400 nm; (b) When the visible light wavelength is 500 nm; (c) When the visible light wavelength is 600 nm; (d) When the visible light wavelength is 700 nm
图8. 可见光波长变化对选取结构的电场强度的影响图。(a) 当可见光波长为400 nm时;(b) 当可见光波长为500 nm时;(c) 当可见光波长为600 nm时;(d) 当可见光波长为700 nm时
限元仿真软件进行建模与仿真优化,使得该太赫兹超材料吸收器在2.84 THz处实现了对太赫兹波99.93%的“完美吸收”。为了进一步扩展太赫兹超材料吸收器的功能,在其金属条带表面随机生长一定尺寸范围内的金纳米颗粒,仿真结果表明,一定尺寸范围内的金纳米颗粒的存在对太赫兹超材料吸收器在太赫兹波段的响应没有太大影响。此外,选取金属条带上的一颗金纳米颗粒作为可见光波段仿真的结构,验证了该结构在可见光波段也有响应。所以,该超材料多级结构在太赫兹和可见光波段均有响应,这对于太赫兹超材料吸收器应用广度的拓展有着一定的影响。
致 谢
感谢国家自然科学基金(No. 61875126)资助。感谢张玲老师在超表面研究中提供了建设性的意见,以及实验室成员在此过程中给予的帮助。
基金项目
国家自然科学基金(No. 61875126)资助。