搜索

x

留言板

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

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

太赫兹多波束调控反射编码超表面

黄若彤 李九生

引用本文:
Citation:

太赫兹多波束调控反射编码超表面

黄若彤, 李九生

Terahertz multibeam modulation reflection-coded metasurface

Huang Ruo-Tong, Li Jiu-Sheng
PDF
HTML
导出引用
  • 已报道的大多数编码超表面仅利用相位或幅度编码进行电磁波调控, 限制了太赫兹波调控灵活性. 本文提出了一种反射超表面单元, 通过相位编码构造超表面, 在圆极化波入射下获得反射波束分裂和偏转功能, 实现对圆极化波束的灵活调控; 同一超表面单元结构利用幅度编码构造超表面在线极化太赫兹波入射下, 实现空间成像功能. 通过相位编码和幅度编码结合构造超表面, 提高了对太赫兹波操控的灵活性, 该编码超表面构造思路可以为太赫兹器件设计提供一种全新思路.
    Most of reported coding metasurfaces only use phase encoding or amplitude encoding to regulate electromagnetic waves, which limits the flexibility of terahertz wave regulation. In this work, a metasurface element structure is proposed. The metasurface element is composed of three layers, i.e. metal pattern structure layer, intermediate medium layer, and metal base layer. According to the geometric phase principle, the phase coverage in the 2π range can be achieved by rotating the metal pattern structure layer under the incidence of the circular-polarized terahertz wave. The metasurface element structure is arranged reasonably by using the phase coding, and the 1-bit and 2-bit phase coding metasurface are designed. First of all, the coding metasurface with interlacing “0” and “1” is designed to generate a double beam reflection under the vertical incidence of circular polarized terahertz waves, while the two-dimensional checkerboard coding metasurface with “0” and “1” generates a symmetrical four-beam reflection. In addition, the metasurface is designed to deflect the reflected beam, and the coding period is changed to design the metasurface to deflect the reflected beam to the specified angle, showing good flexibility. Finally, the convolutional operation is introduced to flexibly regulate the circular polarized beam, and the functions of beam splitting and reflection beam deflection are obtained. The amplitude coded metasurface is designed under theincidence of the online polarized terahertz wave, and the near-field imaging effect can be realized by the amplitude differentiation of polarization reflection. The designed amplitude coded metasurface realizes the function of imaging in space, presenting the designed “CJLU” pattern, which has different imaging effects at different observation locations. When the observation plane distance is 80 μm at the observation frequency of 1.22 THz, the near-field imaging effect is best. In conclusion, we propose a terahertz multibeam modulation reflection-coded metasurface, which combines geometric phase and amplitude variation to achieve different terahertz wave modulation functions under different polarization incident terahertz waves. The results from the simulated near-field radiation model and the far-field radiation model are both in agreement with the theoretical calculation predictions. The designed metasurface provides a degree of freedom method for terahertz wave polarization and phase manipulation, which greatly improves the efficiency of terahertz wave manipulation and has potential applications in terahertz systems.
      通信作者: 李九生, lijsh2008@126.com
    • 基金项目: 国家自然科学基金(批准号: 61871355, 61831012)、浙江省科技厅人才工程 (批准号: 2018R52043)和浙江省重点研发项目 (批准号: 2021C03153, 2022C03166)资助的课题.
      Corresponding author: Li Jiu-Sheng, lijsh2008@126.com
    • Funds: Project supported by the National Natural Science Foundation of China (Grant Nos. 61871355, 61831012), the Talent Project of Zhejiang Provincial Department of Science and Technology, China (Grant No. 2018R52043), and the Zhejiang Key R & D Project of China (Grant Nos. 2021C03153, 2022C03166).
    [1]

    Zhang C, Deng L, Wang L, Chen X, Li S 2021 Appl. Sci. 11 7128Google Scholar

    [2]

    Huang J, Yin X, Xu M, Liu M, Zhang Y, Zhang H 2022 Res. Phys. 33 105204Google Scholar

    [3]

    Sun X, Xu M, Wang G, Song Q, Li Y, Gao X 2022 Appl. Opt. 61 34Google Scholar

    [4]

    Xu Z, Sheng H, Wang Q, Zhou L, Shen Y 2021 SN Appl. Sci. 3 1Google Scholar

    [5]

    Dash S, Liaskos C, Akyildiz IF, Pitsillides A 2020 Mater. Sci. Forum 1009 63Google Scholar

    [6]

    Wang L, Lan F, Zeng H, et al. 2021 IEEE. 1 93Google Scholar

    [7]

    Gong Y, Quan B, Hu F, Wang H, Zhang L, Jiang M 2022 E Low dimens. Syst. Nanostruct. 143 115334Google Scholar

    [8]

    Niu J, Li C, Mo W, Yao Q, Zhu A 2022 J. Phys. D-Appl. Phys. 55 395105Google Scholar

    [9]

    Cui T J, Qi M Q, Wan X, Zhao J, Cheng Q 2014 Light Sci. Appl. 3 218Google Scholar

    [10]

    Li Z, Wang W, Deng S, Qu J, Li Y, Lü B 2022 Opt. Lett. 47 441Google Scholar

    [11]

    Yang D, Wang W, Lü E, Wang H, Liu B, Hou Y 2022 Iscience 25 104824Google Scholar

    [12]

    He C, Song Z 2022 Opt. Express 30 25498Google Scholar

    [13]

    Liu C X, Yang F, Fu X J, Wu J W, Zhang L, Yang J 2021 Adv. Opt. Mater. 9 2100932Google Scholar

    [14]

    Pan W M, Li J S 2021 Opt. Express 29 12918Google Scholar

    [15]

    Zhao D, Tan Z, Zhao H, Fan F, Chang S 2022 Opt. Lett. 47 818Google Scholar

    [16]

    Zheng S, Li C, Fang G 2021 Opt. Express 29 43403Google Scholar

    [17]

    Wang T, Chen B, Wu J, Yang S, Shen Z, Cai J 2021 Appl. Phys. Lett. 118 081101Google Scholar

    [18]

    Li W, Hu X, Wu J, Fan K, Chen B, Zhang C 2022 Light Sci. Appl. 11 1Google Scholar

    [19]

    Lin Q W, Wong H, Huitema L, Crunteanu A 2022 Adv. Opt. Mater. 10 2101699Google Scholar

    [20]

    Ren B, Feng Y, Tang S, Wu J L, Liu B, Song J 2022 Opt. Express 30 16229Google Scholar

    [21]

    Ren B, Feng Y, Tang S, Wang L, Jiang H, Jiang Y 2021 Opt. Express 29 17258Google Scholar

    [22]

    Li J, Cheng Y, Fan J, Chen F, Luo H, Li X 2022 Phys. Lett. 428 127932Google Scholar

    [23]

    Wei J, Qi Y, Zhang B, Ding J, Liu W, Wang X 2022 Opt. Commun. 502 127425Google Scholar

    [24]

    Zang X, Yao B, Chen L, Xie J, Guo X 2021 Light: Adv. Manuf. 2 148Google Scholar

    [25]

    Liu S, Ouyang C, Yao Z, Zhao J, Li Y, Feng L 2022 Opt. Express 30 28158Google Scholar

    [26]

    Qi Y, Zhang B, Ding J, Zhang T, Wang X, Yi Z 2021 Chin. Phys. B. 30 024211Google Scholar

    [27]

    He J, Chen R, Li Y, Chen S, Liu Z, Zhang Q 2021 Appl. Opt. 60 5752Google Scholar

    [28]

    Zhong M, Li J S 2022 Opt. Commun. 511 127997Google Scholar

    [29]

    Kou W, Shi W, Zhang Y, Yang Z, Chen T, Gu J 2021 IEEE. 12 13Google Scholar

    [30]

    Chen D C, Zhu X F, Wei Q, Yao J, Wu D J 2020 J. Phys. D Appl. Phys 53 255501Google Scholar

    [31]

    Saifullah Y, Yang G, Feng X U 2021 J. Radars 10 382Google Scholar

    [32]

    Li S, Li Z, Han B, Huang G, Liu X, Yang H 2022 Front Mater 9 854062Google Scholar

  • 图 1  编码超表面结构示意图 (a) 幅度编码超表面; (b) 相位编码超表面; (c) 超表面单元

    Fig. 1.  Schematic diagram of the proposed coding metasurface: (a) Amplitude coding metasurface; (b) phase coding metasurface; (c) unit cell.

    图 2  圆极化波入射时超表面单元的反射振幅和相位 (a)反射振幅; (b)反射相位

    Fig. 2.  Reflection amplitude and phase of the unit cell under the circularly polarized wave incidence: (a) Reflection amplitude; (b) reflection phase.

    图 3  产生两分束太赫兹反射波的1-bit编码超表面 (a) 编码排布示意图; (b) 三维远场散射图; (c) 归一化反射振幅图

    Fig. 3.  1-bit coding metasurface for generating two-reflected beam: (a) Layout diagram of coding metasurface; (b) three-dimensional far field scattering diagram; (c) normalized reflection amplitude diagram.

    图 4  产生四分束的1-bit编码超表面 (a) 编码排布示意图; (b) 三维远场散射图; (c) 归一化反射振幅图

    Fig. 4.  1-bit coding metasurface for generating quadrant beam: (a) Layout diagram of coding metasurface; (b) three-dimensional far field scattering diagram; (c) normalized reflection amplitude diagram.

    图 5  编码超表面卷积过程示意图 (a)—(c) S1, S2和S3的超表面单元排布示意图; (d)—(f) 0.48 THz处左圆偏振入射时, S1, S2和S3的3D远场散射图; (g)—(i) 0.48 THz处左圆偏振入射时, S1, S2和S3的2D散射图

    Fig. 5.  Convolution process schematic diagram of the coding metasurface: (a)–(c) Layout diagram of S1, S2 and S3 metasurfaces; (d)–(f) 3D far-field scattering diagram of S1, S2 and S3 under left circularly polarized incidence at 0.48 THz; (g)–(i) 2D scattering diagram of S1, S2 and S3 under left circularly polarized incident at 0.48 THz.

    图 6  2-bit反射编码超表面(“00-01-10-11···”) (a) 超表面排布示意图; (b) 超表面结构; (c) 三维远场散射图; (d) 归一化反射振幅图

    Fig. 6.  2-bit reflected coding metasurface (“00-01-10-11···”): (a) Layout diagram of coding metasurface; (b) metasurface structure; (c) three-dimensional far field scattering diagram; (d) normalized reflection amplitude diagram.

    图 7  2-bit反射编码超表面(“0000-0101-1010-1111”) (a) 超表面排布示意图; (b) 超表面结构; (c) 三维远场散射图; (d) 归一化反射振幅图

    Fig. 7.  2-bit reflection encoding metasurface (“0000-0101-1010-1111”): (a) Layout diagram of coding metasurface; (b) metasurface structure; (c) three-dimensional far field scattering diagram; (d) normalized reflection amplitude diagram.

    图 8  线极化波入射时超表面单元产生的反射幅度和反射相位 (a) 反射幅度; (b)反射相位

    Fig. 8.  Reflected amplitude and phase of the unit cell under linearly polarized wave incidence: (a) Reflection amplitude; (b) reflection phase.

    图 9  “CJLU”超表面排布示意图

    Fig. 9.  Layout diagram of "CJLU" metasurface.

    图 10  观测频率1.22 THz处, 观测平面距离为60 μm的近场图像

    Fig. 10.  Near field image on an observation plane of 60 μm distance at 1.22 THz.

    图 11  观测频率1.22 THz处, 观测平面距离为80 μm的近场图像

    Fig. 11.  Near field image on an observation plane of 80 μm distance at 1.22 THz.

    图 12  观测频率1.22 THz处, 观测平面距离为100 μm的近场图像

    Fig. 12.  Near field image on an observation plane of 100 μm distance at 1.22 THz.

    表 1  编码超表面单元

    Table 1.  Unit cell of the proposed coding metasurface.

    1-bit相位编码01
    2-bit相位编码00 011011
    旋转角度/(°)0 4590135
    相位/(°)–137.21 –49.0943.99131.99
    俯视图
    下载: 导出CSV
  • [1]

    Zhang C, Deng L, Wang L, Chen X, Li S 2021 Appl. Sci. 11 7128Google Scholar

    [2]

    Huang J, Yin X, Xu M, Liu M, Zhang Y, Zhang H 2022 Res. Phys. 33 105204Google Scholar

    [3]

    Sun X, Xu M, Wang G, Song Q, Li Y, Gao X 2022 Appl. Opt. 61 34Google Scholar

    [4]

    Xu Z, Sheng H, Wang Q, Zhou L, Shen Y 2021 SN Appl. Sci. 3 1Google Scholar

    [5]

    Dash S, Liaskos C, Akyildiz IF, Pitsillides A 2020 Mater. Sci. Forum 1009 63Google Scholar

    [6]

    Wang L, Lan F, Zeng H, et al. 2021 IEEE. 1 93Google Scholar

    [7]

    Gong Y, Quan B, Hu F, Wang H, Zhang L, Jiang M 2022 E Low dimens. Syst. Nanostruct. 143 115334Google Scholar

    [8]

    Niu J, Li C, Mo W, Yao Q, Zhu A 2022 J. Phys. D-Appl. Phys. 55 395105Google Scholar

    [9]

    Cui T J, Qi M Q, Wan X, Zhao J, Cheng Q 2014 Light Sci. Appl. 3 218Google Scholar

    [10]

    Li Z, Wang W, Deng S, Qu J, Li Y, Lü B 2022 Opt. Lett. 47 441Google Scholar

    [11]

    Yang D, Wang W, Lü E, Wang H, Liu B, Hou Y 2022 Iscience 25 104824Google Scholar

    [12]

    He C, Song Z 2022 Opt. Express 30 25498Google Scholar

    [13]

    Liu C X, Yang F, Fu X J, Wu J W, Zhang L, Yang J 2021 Adv. Opt. Mater. 9 2100932Google Scholar

    [14]

    Pan W M, Li J S 2021 Opt. Express 29 12918Google Scholar

    [15]

    Zhao D, Tan Z, Zhao H, Fan F, Chang S 2022 Opt. Lett. 47 818Google Scholar

    [16]

    Zheng S, Li C, Fang G 2021 Opt. Express 29 43403Google Scholar

    [17]

    Wang T, Chen B, Wu J, Yang S, Shen Z, Cai J 2021 Appl. Phys. Lett. 118 081101Google Scholar

    [18]

    Li W, Hu X, Wu J, Fan K, Chen B, Zhang C 2022 Light Sci. Appl. 11 1Google Scholar

    [19]

    Lin Q W, Wong H, Huitema L, Crunteanu A 2022 Adv. Opt. Mater. 10 2101699Google Scholar

    [20]

    Ren B, Feng Y, Tang S, Wu J L, Liu B, Song J 2022 Opt. Express 30 16229Google Scholar

    [21]

    Ren B, Feng Y, Tang S, Wang L, Jiang H, Jiang Y 2021 Opt. Express 29 17258Google Scholar

    [22]

    Li J, Cheng Y, Fan J, Chen F, Luo H, Li X 2022 Phys. Lett. 428 127932Google Scholar

    [23]

    Wei J, Qi Y, Zhang B, Ding J, Liu W, Wang X 2022 Opt. Commun. 502 127425Google Scholar

    [24]

    Zang X, Yao B, Chen L, Xie J, Guo X 2021 Light: Adv. Manuf. 2 148Google Scholar

    [25]

    Liu S, Ouyang C, Yao Z, Zhao J, Li Y, Feng L 2022 Opt. Express 30 28158Google Scholar

    [26]

    Qi Y, Zhang B, Ding J, Zhang T, Wang X, Yi Z 2021 Chin. Phys. B. 30 024211Google Scholar

    [27]

    He J, Chen R, Li Y, Chen S, Liu Z, Zhang Q 2021 Appl. Opt. 60 5752Google Scholar

    [28]

    Zhong M, Li J S 2022 Opt. Commun. 511 127997Google Scholar

    [29]

    Kou W, Shi W, Zhang Y, Yang Z, Chen T, Gu J 2021 IEEE. 12 13Google Scholar

    [30]

    Chen D C, Zhu X F, Wei Q, Yao J, Wu D J 2020 J. Phys. D Appl. Phys 53 255501Google Scholar

    [31]

    Saifullah Y, Yang G, Feng X U 2021 J. Radars 10 382Google Scholar

    [32]

    Li S, Li Z, Han B, Huang G, Liu X, Yang H 2022 Front Mater 9 854062Google Scholar

  • [1] 陈闻博, 陈鹤鸣. 基于超材料复合结构的太赫兹液晶移相器. 物理学报, 2022, 71(17): 178701. doi: 10.7498/aps.71.20212400
    [2] 惠战强, 高黎明, 刘瑞华, 韩冬冬, 汪伟. 低损耗大带宽双芯负曲率太赫兹光纤偏振分束器. 物理学报, 2022, 71(4): 048702. doi: 10.7498/aps.71.20211650
    [3] 冯龙呈, 杜琛, 杨圣新, 张彩虹, 吴敬波, 范克彬, 金飚兵, 陈健, 吴培亨. 太赫兹实时近场光谱成像研究. 物理学报, 2022, 71(16): 164201. doi: 10.7498/aps.71.20220131
    [4] 刘紫玉, 亓丽梅, 道日娜, 戴林林, 武利勤. 基于VO2的波束可调太赫兹天线. 物理学报, 2022, 71(18): 188703. doi: 10.7498/aps.71.20220817
    [5] 闫志巾, 施卫. 太赫兹GaAs光电导天线阵列辐射特性. 物理学报, 2021, 70(24): 248704. doi: 10.7498/aps.70.20211210
    [6] 惠战强. 低损耗大带宽双芯负曲率太赫兹光纤偏振分束器. 物理学报, 2021, (): . doi: 10.7498/aps.70.20211650
    [7] 陈志文, 佘圳跃, 廖开宇, 黄巍, 颜辉, 朱诗亮. 基于Rydberg原子天线的太赫兹测量. 物理学报, 2021, 70(6): 060702. doi: 10.7498/aps.70.20201870
    [8] 冯正, 王大承, 孙松, 谭为. 自旋太赫兹源:性能、调控及其应用. 物理学报, 2020, 69(20): 208705. doi: 10.7498/aps.69.20200757
    [9] 王大勇, 李兵, 戎路, 赵洁, 王云新, 翟长超. 连续太赫兹波双物距叠层定量相衬成像. 物理学报, 2020, 69(2): 028701. doi: 10.7498/aps.69.20191310
    [10] 李晓楠, 周璐, 赵国忠. 基于反射超表面产生太赫兹涡旋波束. 物理学报, 2019, 68(23): 238101. doi: 10.7498/aps.68.20191055
    [11] 张真真, 黎华, 曹俊诚. 高速太赫兹探测器. 物理学报, 2018, 67(9): 090702. doi: 10.7498/aps.67.20180226
    [12] 闫昕, 梁兰菊, 张璋, 杨茂生, 韦德泉, 王猛, 李院平, 吕依颖, 张兴坊, 丁欣, 姚建铨. 基于石墨烯编码超构材料的太赫兹波束多功能动态调控. 物理学报, 2018, 67(11): 118102. doi: 10.7498/aps.67.20180125
    [13] 汪静丽, 刘洋, 钟凯. 基于领结型多孔光纤的双芯太赫兹偏振分束器. 物理学报, 2017, 66(2): 024209. doi: 10.7498/aps.66.024209
    [14] 张学进, 陆延青, 陈延峰, 朱永元, 祝世宁. 太赫兹表面极化激元. 物理学报, 2017, 66(14): 148705. doi: 10.7498/aps.66.148705
    [15] 冯伟, 张戎, 曹俊诚. 基于石墨烯的太赫兹器件研究进展. 物理学报, 2015, 64(22): 229501. doi: 10.7498/aps.64.229501
    [16] 鲍迪, 沈晓鹏, 崔铁军. 太赫兹人工电磁媒质研究进展. 物理学报, 2015, 64(22): 228701. doi: 10.7498/aps.64.228701
    [17] 鹿文亮, 娄淑琴, 王鑫, 申艳, 盛新志. 基于太赫兹时域光谱技术的伪色彩太赫兹成像的实验研究. 物理学报, 2015, 64(11): 114206. doi: 10.7498/aps.64.114206
    [18] 戴雨涵, 陈小浪, 赵强, 张继华, 陈宏伟, 杨传仁. 太赫兹波段谐振频率可调的开口谐振环结构. 物理学报, 2013, 62(6): 064101. doi: 10.7498/aps.62.064101
    [19] 韩煜, 袁学松, 马春燕, 鄢扬. 波瓣波导谐振腔太赫兹回旋管的研究. 物理学报, 2012, 61(6): 064102. doi: 10.7498/aps.61.064102
    [20] 张显斌, 施 卫. 基于可调谐准高斯波束太赫兹源的成像系统研究. 物理学报, 2008, 57(8): 4984-4990. doi: 10.7498/aps.57.4984
计量
  • 文章访问数:  1093
  • PDF下载量:  61
  • 被引次数: 0
出版历程
  • 收稿日期:  2022-10-14
  • 修回日期:  2022-12-17
  • 上网日期:  2022-12-26
  • 刊出日期:  2023-03-05

太赫兹多波束调控反射编码超表面

  • 中国计量大学太赫兹技术与应用研究所, 杭州 310018
  • 通信作者: 李九生, lijsh2008@126.com
    基金项目: 国家自然科学基金(批准号: 61871355, 61831012)、浙江省科技厅人才工程 (批准号: 2018R52043)和浙江省重点研发项目 (批准号: 2021C03153, 2022C03166)资助的课题.

摘要: 已报道的大多数编码超表面仅利用相位或幅度编码进行电磁波调控, 限制了太赫兹波调控灵活性. 本文提出了一种反射超表面单元, 通过相位编码构造超表面, 在圆极化波入射下获得反射波束分裂和偏转功能, 实现对圆极化波束的灵活调控; 同一超表面单元结构利用幅度编码构造超表面在线极化太赫兹波入射下, 实现空间成像功能. 通过相位编码和幅度编码结合构造超表面, 提高了对太赫兹波操控的灵活性, 该编码超表面构造思路可以为太赫兹器件设计提供一种全新思路.

English Abstract

    • 超表面是由二维元素微结构周期性或非周期性排列的人工电磁材料, 能够操纵电磁波的振幅和相位, 实现偏振转换[1-4] 、完美吸收[5]振幅和相位调制[6-8] 等特殊功能. 2014年, Cui等[9]提出了微波的数字编码超表面概念. 近年来, 编码超表面在太赫兹频率波段也受到了广泛关注, 研究人员通过使用编码超表面实现了太赫兹波束偏转[10-15]、光谱成像[16-18]、偏振操纵[19-26]、聚焦[27-29]等功能. 2020年, Chen等[30]提出用于涡旋太赫兹波束操纵和聚焦的多功能反射编码超表面. 2021年, Saifullah等[31]提出一种结构为“四叶草形”的编码超表面, 并实现了超宽带漫散射. 2022年, Li等[32]设计了一种左右圆极化太赫兹多功能编码超表面. 上述报道的超表面都是在单一偏振态太赫兹波入射下通过相位、幅度等参数的差异编码构造超表面实现波束分裂、波束偏转和涡旋波束操纵等功能, 由于只利用电磁波单一偏振态, 所以对太赫兹波调控灵活性降低.

      本文提出太赫兹多波束调控反射编码超表面, 它将几何相位和幅度变化相结合, 在不同偏振的太赫兹波入射下实现不同的太赫兹波调控功能. 模拟的近场和远场辐射模式, 结果与理论计算预测一致, 设计的超表面提供了一种太赫兹波偏振和相位操控自由度方法, 极大地提高了太赫兹波操纵的效率.

    • 本文提出的编码超表面结构如图1所示, 其中图1(a)为所设计的幅度编码超表面, 利用入射线偏振 (LP) 波产生反射波的幅度差值进行编码, 图1(b)为相位编码超表面, 通过旋转顶层金属图案利用入射圆偏振 (CP) 波产生PB相位 (几何相位)进行编码. 超表面单元结构如图1(c)所示, 单元周期P = 65 μm, 由3层介质组成, 分别为金属图案(外径r1 = 30 μm, 内径r2 = 25 μm, 厚度为d1 = 0.125 μm), 中间硅介质层(εr = 11.9, 厚度h = 25 μm), 金属板薄膜层(厚度d2 = 2 μm). α是单元结构金属图案沿逆时针方向旋转角度.

      图  1  编码超表面结构示意图 (a) 幅度编码超表面; (b) 相位编码超表面; (c) 超表面单元

      Figure 1.  Schematic diagram of the proposed coding metasurface: (a) Amplitude coding metasurface; (b) phase coding metasurface; (c) unit cell.

      利用几何相位原理构造编码超表面, 编码单元所对应的俯视图和相位如表1所示. 对于1-bit数字编码超表面, 数字“0”和“1”代表反射波相对于入射波的相位差为180°的两种编码单元, 对应单元的反射相位为 –137.21°和 43.99°; 对于2-bit数字编码超表面, 它由4个编码单元组成, 分别对应 –137.21°, –49.09°, 43.99°和 131.99°的反射相位响应, 对应的数字 “00”, “01”, “10”和“11”代表反射波相对于入射波的相位差为90°的4个编码超表面单元. 编码单元对应的反射幅度和反射相位如图2(a), (b)所示, 可以看出, 在0.48—0.68 THz之间, 太赫兹波反射幅度均大于0.8, 相邻编码单元反射波相对于入射波的相位差接近90°, 满足2-bit编码超表面设计要求.

      1-bit相位编码01
      2-bit相位编码00 011011
      旋转角度/(°)0 4590135
      相位/(°)–137.21 –49.0943.99131.99
      俯视图

      表 1  编码超表面单元

      Table 1.  Unit cell of the proposed coding metasurface.

      图  2  圆极化波入射时超表面单元的反射振幅和相位 (a)反射振幅; (b)反射相位

      Figure 2.  Reflection amplitude and phase of the unit cell under the circularly polarized wave incidence: (a) Reflection amplitude; (b) reflection phase.

      编码超表面一般包含有限种单元, 凭借离散化的数字编码序列来调控电磁波, 核心思想是将数字化的二进制编码融入超表面的结构、电磁参数、功能等设计的各个方面. 编码超表面由N × N个相同尺寸为D的方形栅格构成, 每个栅格填充由“0”和“1”单元构成的子阵列, 其散射相位设为φ = (m, n). 在平面波垂直入射的情形下, 超表面的远场函数可表示为

      $ \begin{split} f(\theta ,\varphi ) =\;& {f_{\text{e}}}(\theta ,\varphi )\sum\limits_{m = 1}^N \sum\limits_{n = 1}^N \exp \Bigr\{ - {\text{i}} \big\{ \varphi (m,n) \\ & + KD\sin \theta [ (m - 1/2)\cos \varphi \\ & + (n - 1/2)\sin \varphi \big] \big\} \Bigr\}, \end{split} $

      式中θφ 为任意方向上的俯仰角和方位角, fe(θ, φ)为栅格方向函数, 其中θ = arcsin (λ/Γ), λ是自由空间波长, Γ 表示编码序列周期长度. 根据广义斯涅耳定律和超表面散射的远场函数, 反射波束的方位角(φ)和单元“0”或“1”编码粒子的长度(Dx)和宽度(Dy)满足以下关系:

      $ \begin{split} & \varphi = \pm \arctan ({D_x}/{D_y}) \text{, } \\ & \varphi = {\text{π }} \pm \arctan ({D_x}/{D_y}). \end{split} $

    • 圆极化太赫兹波的垂直入射下, “0”和“1”交错的编码超表面将产生双波束反射, 而用“0”和“1”二维棋盘格编码时, 编码超表面会产生对称的四波束反射. 首先, 设计了如图3(a)所示的编码超表面, 每个“0”和“1”单元格都是由6×6个编码粒子组成, 编码序列“000000-111111···”沿 x 方向呈周期性排列, 沿y方向排布, 编码周期Γ = 780 µm. 当圆极化太赫兹波垂直入射到超表面上时, 在0.5 THz处, 入射太赫兹波沿着x轴被分为两束对称的反射波, 其三维远场如图3(b)所示. 图3(c)描绘了在0.5 THz处反射模式下两分束的归一化幅度曲线, 两个反射波束的峰值分别在–50°和50°出现, 即偏转角为50°. 利用公式计算θ = arcsin (λ/Γ) = 50.2°, 结合归一化反射幅度曲线和三维远场图可以得出模拟与理论计算结果符合.

      图  3  产生两分束太赫兹反射波的1-bit编码超表面 (a) 编码排布示意图; (b) 三维远场散射图; (c) 归一化反射振幅图

      Figure 3.  1-bit coding metasurface for generating two-reflected beam: (a) Layout diagram of coding metasurface; (b) three-dimensional far field scattering diagram; (c) normalized reflection amplitude diagram.

      此外, 本文还排布了一个棋盘编码超表面, 该结构周期性地使用数字序列“000000-111111/111111-000000··· ” 沿 x 方向呈周期性排布, 沿 y 方向排列, 编码周期Γ = 1103 µm, 如图4(a)所示. 当圆极化太赫兹波垂直入射到编码超表面上时, 在0.55 THz处的三维远场和归一化反射振幅图如图4(b), (c)所示. 可以观察到, 入射太赫兹波被分为四束对称的反射波.

      图  4  产生四分束的1-bit编码超表面 (a) 编码排布示意图; (b) 三维远场散射图; (c) 归一化反射振幅图

      Figure 4.  1-bit coding metasurface for generating quadrant beam: (a) Layout diagram of coding metasurface; (b) three-dimensional far field scattering diagram; (c) normalized reflection amplitude diagram.

      为了更加灵活地获得任意散射角调控功能, 本文将编码序列进行基本的卷积运算, 不同编码序列的加法卷积运算可以实现不同的编码周期, 从而获得对入射太赫兹光波不同的散射角度或任意角度. 对于“00”, “01”, “10”, “11”的编码粒子, 每一个编码单元由6×6个编码粒子组成, 四位卷积计算可以执行为 “00 + 00 = 00”, “00 + 11 = 11”, “10 + 10 = 00”, “11+ 10 = 01” 和 “11 + 11 = 10”, 编码超表面的新卷积序列“S3 = S1 + S2”. 结合数字编码原理和卷积定理, 所设计的编码超表面序列S1, S2和S3如图5(a)(c)所示. 编码超表面序列S1(“00-01-10-11…”)沿x方向呈周期性排列, 沿y方向排布, 当左圆偏振太赫兹波入射到超表面S1时, 在0.48 THz处产生了反射波束的偏转, 相应的三维(3D)和二维(2D)散射图如图5(d), (g)所示, 入射太赫兹波以θ = arcsin (λ/Γ) = 24°的反射波束角度偏转. 编码超表面序列S2 (“00-10”···)是沿y方向呈周期性排列、沿x方向排布的编码超表面, 其3D和2D散射图如图5(e), (h)所示. 可以看出, 太赫兹波被分成两个对称的反射波束, 偏转角为θ = arcsin (λ/Γ) = 53.25°. 编码超表面序列S3 (“00-01-10-11···\10-11-00-01···”) 由S1和S2叠加组成, 根据傅里叶变换的卷积定理, 卷积编码超表面S3的远场散射图由两个编码序列的远场散射图叠加组成, 编码超表面序列S3的3D和2D散射图如图5(f), (i)所示, 可见当左圆极化太赫兹波照射到超表面S3上时, 其反射波束分离成了4个偏转光束, 卷积编码超表面可以实现对反射波束的偏转和分裂.

      图  5  编码超表面卷积过程示意图 (a)—(c) S1, S2和S3的超表面单元排布示意图; (d)—(f) 0.48 THz处左圆偏振入射时, S1, S2和S3的3D远场散射图; (g)—(i) 0.48 THz处左圆偏振入射时, S1, S2和S3的2D散射图

      Figure 5.  Convolution process schematic diagram of the coding metasurface: (a)–(c) Layout diagram of S1, S2 and S3 metasurfaces; (d)–(f) 3D far-field scattering diagram of S1, S2 and S3 under left circularly polarized incidence at 0.48 THz; (g)–(i) 2D scattering diagram of S1, S2 and S3 under left circularly polarized incident at 0.48 THz.

    • 2-bit编码超表面设计是利用圆极化太赫兹波入射控制反射波束, 实现波束散射角的偏转. 预设的2-bit编码超表面“00-01-10-11···”沿x 方向周期性排布, 沿 y 方向排列, 每个“00” “01” “10” “11”单元格都是由6×6个编码粒子组成, 周期为 Γ = 1560 µm, 超表面单元排布如图6(a), (b)所示. 在0.48 THz处, 当左圆极化太赫兹波垂直入射到编码超表面上时, 太赫兹反射波束发生了偏转, 其三维远场图及对应的归一化幅度曲线分别如图6(c), (d)所示. 从图中可以看出, 相对于z轴的正方向, 反射波束的偏转角接近24°, 偏转角可以利用公式计算出θ = arcsin (λ/Γ) = 23.62°, 仿真结果与计算结果一致. 类似地, 设计的其他 2-bit编码超表面编码序列“0000-0101-1010-1111”沿 x 方向周期性排布, 沿y方向排列, 每个“00” “01” “10” “11”单元格都是由6×6个编码粒子组成, 周期 Γ = 3120 µm, 其超表面单元排布如图7(a), (b)所示. 在0.48 THz处, 当左圆极化太赫兹波垂直入射到编码超表面, 太赫兹反射波束发生了偏转. 编码超表面在反射模式下太赫兹波散射主瓣的三维远场图和对应的归一化反射幅度曲线图如图7(c), (d)所示, 可以看出相对于z轴正方向出现了12°偏转, 理论计算的θ = 11.56°, 这与图7(c), (d)所示的模拟结果非常一致. 结果表明, 可以通过改变编码周期来设计超表面使其反射波束偏转到指定的角度.

      图  6  2-bit反射编码超表面(“00-01-10-11···”) (a) 超表面排布示意图; (b) 超表面结构; (c) 三维远场散射图; (d) 归一化反射振幅图

      Figure 6.  2-bit reflected coding metasurface (“00-01-10-11···”): (a) Layout diagram of coding metasurface; (b) metasurface structure; (c) three-dimensional far field scattering diagram; (d) normalized reflection amplitude diagram.

      图  7  2-bit反射编码超表面(“0000-0101-1010-1111”) (a) 超表面排布示意图; (b) 超表面结构; (c) 三维远场散射图; (d) 归一化反射振幅图

      Figure 7.  2-bit reflection encoding metasurface (“0000-0101-1010-1111”): (a) Layout diagram of coding metasurface; (b) metasurface structure; (c) three-dimensional far field scattering diagram; (d) normalized reflection amplitude diagram.

    • 当线极化太赫兹波入射到超表面单元时, 反射幅度和反射相位的曲线如图8(a), (b)所示. 线极化波入射到超表面单元时, 两个不同超表面单元的反射幅度差异较大, 利用偏振反射幅度差异化特性可以实现近场成像效应. 编码超表面设计为具有两个灰度级(即亮和暗)对应于“1”和“0”的编码单元, 利用振幅编码显示超表面. 幅度高的编码为“1”, 幅度低的编码为“0”, 其编码单元结构如图8(a)所示.

      图  8  线极化波入射时超表面单元产生的反射幅度和反射相位 (a) 反射幅度; (b)反射相位

      Figure 8.  Reflected amplitude and phase of the unit cell under linearly polarized wave incidence: (a) Reflection amplitude; (b) reflection phase.

      设计 “CJLU” 图案由两种不同类型的超表面单元结构分别排布在字母方框内外两块区域, 字母部分选择用幅度编码为“1”的单元排布, 字母以外的其余部分用幅度编码为“0”的单元排布, 编码超表面由32×32个单元组成, 如图9(a)(d)所示, 超表面设置的 “CJLU” 编码图案的轮廓在电场中明显显示出来. 红色部分对应高幅度的编码单元, 蓝色部分对应低幅度的编码单元. 在观测频率为1.22 THz时, 观测平面距离编码超表面60, 80和100 μm的近场图像显示分别如图10图12所示. 编码图案的近场图像显示随着编码图案单元的增加而增大, “CJLU”编码图案轮廓边缘存在些许的不规则, 是由于两个编码单元幅度差异分布不均匀引起的电场能量分布不均匀, 观测距离为80 μm时可得到最佳的近场图像. 总体而言, 仿真得到的成像效果与预设图像大小、位置、轮廓方面的模拟结果较符合, 验证了利用偏振反射幅度差异化特性可以实现近场成像效应.

      图  9  “CJLU”超表面排布示意图

      Figure 9.  Layout diagram of "CJLU" metasurface.

      图  10  观测频率1.22 THz处, 观测平面距离为60 μm的近场图像

      Figure 10.  Near field image on an observation plane of 60 μm distance at 1.22 THz.

      图  11  观测频率1.22 THz处, 观测平面距离为80 μm的近场图像

      Figure 11.  Near field image on an observation plane of 80 μm distance at 1.22 THz.

      图  12  观测频率1.22 THz处, 观测平面距离为100 μm的近场图像

      Figure 12.  Near field image on an observation plane of 100 μm distance at 1.22 THz.

    • 本文结合几何相位和幅度两种参量的变化, 提出了一种太赫兹多波束调控反射编码超表面, 在圆偏振入射太赫兹光波作用下, 所设计的超表面结构在0.48—0.68 THz频段内可以实现太赫兹波束分裂和反射模式下波束偏转, 利用卷积运算可以对圆偏振入射太赫兹光波产生反射多波束调控. 在线偏振态太赫兹光波入射下, 所设计的幅度编码超表面可以在1.22 THz处实现成像功能. 研究结果表明, 所设计的超表面对太赫兹波偏振和相位操控自由度提供了一种新思路, 为实现太赫兹波前的动态操作开辟了一条新途径, 在太赫兹系统具有广阔的应用潜力.

参考文献 (32)

目录

    /

    返回文章
    返回