搜索

x

留言板

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

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

基于反射超表面产生太赫兹涡旋波束

李晓楠 周璐 赵国忠

引用本文:
Citation:

基于反射超表面产生太赫兹涡旋波束

李晓楠, 周璐, 赵国忠

Terahertz vortex beam generation based on reflective metasurface

Li Xiao-Nan, Zhou Lu, Zhao Guo-Zhong
PDF
HTML
导出引用
  • 具有螺旋波前的电磁波是携带轨道角动量的涡旋波束, 涡旋波束存在的相位奇点, 使其在微粒操控和通讯等领域有特殊的应用. 本文提出了一种基于反射型超表面的太赫兹宽带涡旋波束产生器, 该器件由超表面-电介质-金属三层结构构成, 顶层为两个正交I形金属结构单元组成的超表面, 中间层是聚酰亚胺介质, 最底层为金属. 通过对超表面单元结构参数的优化设计, 可以实现在不同旋转角度下反射波振幅高达90%以上, 同时反射波相位随旋转角线性变化的目的. 进一步利用这些单元结构按照Pancharatnam-Berry相位原理进行超表面布阵, 在0.8—1.4 THz频率范围内, 可以将圆偏振太赫兹波束转换为具有轨道角动量的涡旋波束, 这一器件的工作带宽相对较宽, 结构简单, 转换效率高, 在太赫兹涡旋波束产生方面具有潜在的应用价值.
    The electromagnetic wave with spiral wavefront is a vortex beam carrying orbital angular momentum. The phase singularity of the vortex beam has special applications in the fields of particle manipulation and communication. In this paper, a terahertz (THz) wide-band vortex beam generator based on reflective metasurface is proposed and simulated. The device consists of a metasurface-dielectric-metal three-layer structure, and the top layer is a metasurface composed of two orthogonal I-shaped metal structural units. The intermediate layer of polyimide medium, and the bottom layer is of metal as a reflecting plate. The CST microwave studio is used to simulate the reflection performance of unit cell. The structure parameters are optimized to obtain the better performance. A set of optimed structure parameters is determined. According to the phase principle of Pancharatnam-Berry (P-B), by rotating the angle of the top-layer I-type metal structure, the reflection amplitudes of the unit cell structure at different rotation angles are required to approximately equal while the phase changes linearly with rotation angle and reaches a range of 2lπ for the topological charge number l. These cell structures are arranged according to the phase principle mentioned above. The metasurfaces of different topological charge numbers are designed to generate the corresponding vortex beams. In this paper, the metasurfaces with topological charge numbers 1 and 2 are designed. The reflection amplitude and phase of the circularly polarized THz beam incident vertically on the metasurface are simulated by using CST microwave studio. The simulation results show that in a frequency range of 0.8−1.4 THz, the metasurface can convert the circularly polarized terahertz beam into a vortex beam with a different topological charge number. In addition, in order to illustrate that the metasurface designed can produce a higher topological charge number of vortex beam, a metasurface with a topological charge number of 3 is designed as an example. The reflection amplitude and phase of the circularly polarized THz beam at a frequency of 1.1 THz is simulated. The results show that the designed metasurface can produce a vortex beam with a topological charge number of 3. The higher topological charges of vortex beam can also be generated according to the corresponding phase arrangement. The device has a relatively wide operating bandwidth, simple structure, high conversion efficiency, and has the potential application in terahertz vortex beam generation.
      通信作者: 赵国忠, guozhong-zhao@126.com
    • 基金项目: 国家自然科学基金(批准号: 5307625130)资助的课题
      Corresponding author: Zhao Guo-Zhong, guozhong-zhao@126.com
    • Funds: Project supported by the National Natural Science Foundation of China (Grant No. 5307625130)
    [1]

    Beard M C, Turner G M, Schmuttenmaer C A 2002 J. Phys. Chem. B 106 7146Google Scholar

    [2]

    Vieweg N, Fischer B M, Reuter M, Kula P, Dabrowsk R, Celik M A, Frenking G, Koch M, Jepsen P U 2012 Opt. Express 20 28249Google Scholar

    [3]

    Janek M, Zich D, Naftaly M 2014 Mater. Chem. Phys. 145 278Google Scholar

    [4]

    Hu B B, Nuss M C 1995 Opt. Lett. 20 1716Google Scholar

    [5]

    Mittleman D M, Gupta M, Neelamani R, Baraniuk R G, Rudd J V, Koch M 1999 Appl. Phys. B 68 1085Google Scholar

    [6]

    Kohler R, Tredicucci A, Beltram F, Beere H E, Linfield E H, Davies A G, Ritchie D A, Lotti R C, Rossi F 2002 Nature 417 156Google Scholar

    [7]

    Zhang Z W, Wang K J, Lei Y, Zhang Z Y, Zhao Y M, Li C Y, Gu A, Shi N C, Zhao K, Zhan H L, Zhang C L 2015 Science China 58 124202Google Scholar

    [8]

    Heljo V P, Nordberg A, Tenho M, Virtanen T, Jouppila K, Salonen J, Maunu S L, Juppo A M 2012 Pharm. Res. 29 2684Google Scholar

    [9]

    Kirilenko M S, Khonina S N 2013 Optical Memory and Neural Networks 22 81Google Scholar

    [10]

    Chavez-Cerda S, Padgett M J, Allison I, New G H C, Gutierrez-Vega J C, Neil A T O, Vicar I M, Courtial J 2002 J. Optics B: Quantum Semiclass Opt. 4 S52Google Scholar

    [11]

    Genevet P, Lin J, Kats M A, Capasso F 2012 Nature Comm. 3 1Google Scholar

    [12]

    Mohammadi S M, Daldorff L K S, Bergman J E S, Karlsson R L, Thidé B, Forozesh K, Carozzi T D, Isham B 2010 IEEE Trans. Antenn. Propag. 58 565Google Scholar

    [13]

    Thide B, Then H, SjoHolm J, Palmer K, Bergman J, Carozzi T D, Istomin Y N, Ibragimov N H, Khamitova R 2007 Phys. Rev. Lett. 99 087701Google Scholar

    [14]

    付亚男, 张新群, 赵国忠, 李永花, 于佳怡 2017 物理学报 66 180701

    Fu Y N, Zhang X Q, Zhao G Z, Li Y H, Yu J Y 2017 Acta Phys. Sin. 66 180701

    [15]

    李永花, 周璐, 赵国忠 2018 中国激光 45 0314001Google Scholar

    Li Y H, Zhou L, Zhao G Z 2018 Chin. J. Lasers 45 0314001Google Scholar

    [16]

    Li H, Xiao B Y, Huang X J, Yang H L 2015 Phys. Scr. 90 1Google Scholar

    [17]

    Enoch S, Tayeb G, Sabouroux P, Guerin N, Vincent P 2002 Phys. Rev. Lett. 89 213902Google Scholar

    [18]

    Huang J, Pogorzelski R J 1998 IEEE Trans. Antenn. Propag. 46 650Google Scholar

    [19]

    Martynyuk A E, Martinez-Lopez J I, Martynyuk N A 2004 IEEE Trans. Antenn. Propag. 52 142Google Scholar

    [20]

    周璐, 赵国忠, 李晓楠 2019 物理学报 68 108701

    Zhou L, Zhao G Z, Li X N 2019 Acta Phys. Sin. 68 108701

    [21]

    Genevet P, Yu N F, Aieta F, Lin J, Kats M A, et al. 2012 Appl. Phys. Lett. 100 013101Google Scholar

    [22]

    Yu N, Genevet P, Kats M A, Aieta F, Tetienne J P, Capasso F, Gaburro Z 2011 Science 334 333Google Scholar

    [23]

    Zhang K, Yuan Y Y, Zhang D W, Ding X M, Rstni B, Burokur S N, Lu M J, Tang J, Wu Q 2018 Opt. Express 26 1351Google Scholar

    [24]

    Luo W J, Sun S L, Xu H X, He Q, Zhou L 2017 Phys. Rev. Appl. 7 044033Google Scholar

    [25]

    Xu H H, Wang G M, Cai T, Xiao J, Zhuang Y Q 2016 Opt. Express 24 27836Google Scholar

  • 图 1  超表面单元结构示意图 (a)顶视图; (b)侧视图

    Fig. 1.  Schematic diagram of the unit cell of metasurface: (a) Top view; (b) side view.

    图 2  线偏振光入射单元结构产生的反射相位和振幅谱 (a)相位谱; (b)振幅谱

    Fig. 2.  Reflected phase and amplitude spectrum produced by linearly polarized light incident unit call: (a) Phase; (b) amplitude.

    图 3  仿真得到的不同转角的相位和反射振幅

    Fig. 3.  The phase and reflection amplitude of different corners simulated.

    图 4  不同转角单元结构的反射相位和振幅谱 (a) LCP入射的相位谱; (b) RCP入射的相位谱; (c) LCP入射的振幅谱; (d) RCP入射的振幅谱

    Fig. 4.  Reflective phase and amplitude spectra of unit cell structure under different rotation angle: (a) Phase spectra at LCP incident; (b) phase spectra at RCP incident; (c) amplitude spectra at LCP incident; (d) amplitude spectra at RCP incident.

    图 5  两种用于产生拓扑荷数分别为 (a) l = 1和(b) l = 2的涡旋波束超表面结构

    Fig. 5.  Two kinds of metasurface structures for generating the vortex beam with topological charge numbers (a) l = 1 and (b) l = 2.

    图 6  超表面产生l = 1和l = 2的涡旋波束反射振幅和相位分布 LCP入射l = 1超表面的(a)振幅分布和(b)相位分布; RCP入射l = 2超表面的(c)振幅分布和(d)相位分布

    Fig. 6.  Reflective amplitude and phase distributions of vortex beams with l = 1 and l = 2 generated by metasurface. LCP incident l = 1 metasurface: (a) amplitude distribution and (b) phase distribution; RCP incident l = 2 metasurface: (c) amplitude distribution and (d) phase distribution.

    图 7  不同频率下l = 1和l = 3超表面产生的反射涡旋波束振幅、相位分布 l = 1超表面(a) 0.8 THz频率下振幅分布, (b) 0.8 THz频率下相位分布, (c) 1.4 THz频率下振幅分布, (d) 1.4 THz频率下相位分布; l = 3超表面(e) 1.1 THz频率下振幅分布, (f) 1.1 THz频率下相位分布

    Fig. 7.  The amplitude and phase distribution of reflective vortex beam generated by the LCP incident l = 1 and l = 3 metasurface at different frequencies. l = 1: (a) amplitude distribution at 0.8 THz, (b) phase distribution at 0.8 THz, (c) amplitude distribution at 1.4 THz, (d) phase distribution at 1.4 THz. l = 3: (e) amplitude distribution at 1.1 THz, (f) phase distribution at 1.1 THz.

  • [1]

    Beard M C, Turner G M, Schmuttenmaer C A 2002 J. Phys. Chem. B 106 7146Google Scholar

    [2]

    Vieweg N, Fischer B M, Reuter M, Kula P, Dabrowsk R, Celik M A, Frenking G, Koch M, Jepsen P U 2012 Opt. Express 20 28249Google Scholar

    [3]

    Janek M, Zich D, Naftaly M 2014 Mater. Chem. Phys. 145 278Google Scholar

    [4]

    Hu B B, Nuss M C 1995 Opt. Lett. 20 1716Google Scholar

    [5]

    Mittleman D M, Gupta M, Neelamani R, Baraniuk R G, Rudd J V, Koch M 1999 Appl. Phys. B 68 1085Google Scholar

    [6]

    Kohler R, Tredicucci A, Beltram F, Beere H E, Linfield E H, Davies A G, Ritchie D A, Lotti R C, Rossi F 2002 Nature 417 156Google Scholar

    [7]

    Zhang Z W, Wang K J, Lei Y, Zhang Z Y, Zhao Y M, Li C Y, Gu A, Shi N C, Zhao K, Zhan H L, Zhang C L 2015 Science China 58 124202Google Scholar

    [8]

    Heljo V P, Nordberg A, Tenho M, Virtanen T, Jouppila K, Salonen J, Maunu S L, Juppo A M 2012 Pharm. Res. 29 2684Google Scholar

    [9]

    Kirilenko M S, Khonina S N 2013 Optical Memory and Neural Networks 22 81Google Scholar

    [10]

    Chavez-Cerda S, Padgett M J, Allison I, New G H C, Gutierrez-Vega J C, Neil A T O, Vicar I M, Courtial J 2002 J. Optics B: Quantum Semiclass Opt. 4 S52Google Scholar

    [11]

    Genevet P, Lin J, Kats M A, Capasso F 2012 Nature Comm. 3 1Google Scholar

    [12]

    Mohammadi S M, Daldorff L K S, Bergman J E S, Karlsson R L, Thidé B, Forozesh K, Carozzi T D, Isham B 2010 IEEE Trans. Antenn. Propag. 58 565Google Scholar

    [13]

    Thide B, Then H, SjoHolm J, Palmer K, Bergman J, Carozzi T D, Istomin Y N, Ibragimov N H, Khamitova R 2007 Phys. Rev. Lett. 99 087701Google Scholar

    [14]

    付亚男, 张新群, 赵国忠, 李永花, 于佳怡 2017 物理学报 66 180701

    Fu Y N, Zhang X Q, Zhao G Z, Li Y H, Yu J Y 2017 Acta Phys. Sin. 66 180701

    [15]

    李永花, 周璐, 赵国忠 2018 中国激光 45 0314001Google Scholar

    Li Y H, Zhou L, Zhao G Z 2018 Chin. J. Lasers 45 0314001Google Scholar

    [16]

    Li H, Xiao B Y, Huang X J, Yang H L 2015 Phys. Scr. 90 1Google Scholar

    [17]

    Enoch S, Tayeb G, Sabouroux P, Guerin N, Vincent P 2002 Phys. Rev. Lett. 89 213902Google Scholar

    [18]

    Huang J, Pogorzelski R J 1998 IEEE Trans. Antenn. Propag. 46 650Google Scholar

    [19]

    Martynyuk A E, Martinez-Lopez J I, Martynyuk N A 2004 IEEE Trans. Antenn. Propag. 52 142Google Scholar

    [20]

    周璐, 赵国忠, 李晓楠 2019 物理学报 68 108701

    Zhou L, Zhao G Z, Li X N 2019 Acta Phys. Sin. 68 108701

    [21]

    Genevet P, Yu N F, Aieta F, Lin J, Kats M A, et al. 2012 Appl. Phys. Lett. 100 013101Google Scholar

    [22]

    Yu N, Genevet P, Kats M A, Aieta F, Tetienne J P, Capasso F, Gaburro Z 2011 Science 334 333Google Scholar

    [23]

    Zhang K, Yuan Y Y, Zhang D W, Ding X M, Rstni B, Burokur S N, Lu M J, Tang J, Wu Q 2018 Opt. Express 26 1351Google Scholar

    [24]

    Luo W J, Sun S L, Xu H X, He Q, Zhou L 2017 Phys. Rev. Appl. 7 044033Google Scholar

    [25]

    Xu H H, Wang G M, Cai T, Xiao J, Zhuang Y Q 2016 Opt. Express 24 27836Google Scholar

  • [1] 于博, 庄书磊, 王正心, 王曼诗, 郭兰军, 李鑫煜, 郭文瑞, 苏文明, 龚诚, 刘伟伟. 基于纳米印刷技术的双螺旋太赫兹可调超表面. 物理学报, 2022, 71(11): 117801. doi: 10.7498/aps.71.20212408
    [2] 高喜, 唐李光. 基于双层超表面的宽带、高效透射型轨道角动量发生器. 物理学报, 2021, 70(3): 038101. doi: 10.7498/aps.70.20200975
    [3] 孙胜, 阳棂均, 沙威. 基于反射超表面的偏馈式涡旋波产生装置. 物理学报, 2021, 70(19): 198401. doi: 10.7498/aps.70.20210681
    [4] 李国强, 施宏宇, 刘康, 李博林, 衣建甲, 张安学, 徐卓. 基于超表面的多波束多模态太赫兹涡旋波产生. 物理学报, 2021, 70(18): 188701. doi: 10.7498/aps.70.20210897
    [5] 龙洁, 李九生. 相变材料与超表面复合结构太赫兹移相器. 物理学报, 2021, 70(7): 074201. doi: 10.7498/aps.70.20201495
    [6] 周璐, 赵国忠, 李晓楠. 基于双开口谐振环超表面的宽带太赫兹涡旋光束产生. 物理学报, 2019, 68(10): 108701. doi: 10.7498/aps.68.20182147
    [7] 高强, 王晓华, 王秉中. 基于宽带立体超透镜的远场超分辨率成像. 物理学报, 2018, 67(9): 094101. doi: 10.7498/aps.67.20172608
    [8] 李小兵, 陆卫兵, 刘震国, 陈昊. 基于可调石墨烯超表面的宽角度动态波束控制. 物理学报, 2018, 67(18): 184101. doi: 10.7498/aps.67.20180592
    [9] 闫昕, 梁兰菊, 张璋, 杨茂生, 韦德泉, 王猛, 李院平, 吕依颖, 张兴坊, 丁欣, 姚建铨. 基于石墨烯编码超构材料的太赫兹波束多功能动态调控. 物理学报, 2018, 67(11): 118102. doi: 10.7498/aps.67.20180125
    [10] 宁仁霞, 鲍婕, 焦铮. 基于石墨烯超表面的宽带电磁诱导透明研究. 物理学报, 2017, 66(10): 100202. doi: 10.7498/aps.66.100202
    [11] 付亚男, 张新群, 赵国忠, 李永花, 于佳怡. 基于谐振环的太赫兹宽带偏振转换器件研究. 物理学报, 2017, 66(18): 180701. doi: 10.7498/aps.66.180701
    [12] 张银, 冯一军, 姜田, 曹杰, 赵俊明, 朱博. 基于石墨烯的太赫兹波散射可调谐超表面. 物理学报, 2017, 66(20): 204101. doi: 10.7498/aps.66.204101
    [13] 李唐景, 梁建刚, 李海鹏, 牛雪彬, 刘亚峤. 基于单层线-圆极化转换聚焦超表面的宽带高增益圆极化天线设计. 物理学报, 2017, 66(6): 064102. doi: 10.7498/aps.66.064102
    [14] 侯海生, 王光明, 李海鹏, 蔡通, 郭文龙. 超薄宽带平面聚焦超表面及其在高增益天线中的应用. 物理学报, 2016, 65(2): 027701. doi: 10.7498/aps.65.027701
    [15] 李唐景, 梁建刚, 李海鹏. 基于单层反射超表面的宽带圆极化高增益天线设计. 物理学报, 2016, 65(10): 104101. doi: 10.7498/aps.65.104101
    [16] 李勇峰, 张介秋, 屈绍波, 王甲富, 吴翔, 徐卓, 张安学. 二维宽带相位梯度超表面设计及实验验证. 物理学报, 2015, 64(9): 094101. doi: 10.7498/aps.64.094101
    [17] 范亚, 屈绍波, 王甲富, 张介秋, 冯明德, 张安学. 基于交叉极化旋转相位梯度超表面的宽带异常反射. 物理学报, 2015, 64(18): 184101. doi: 10.7498/aps.64.184101
    [18] 余积宝, 马华, 王甲富, 冯明德, 李勇峰, 屈绍波. 基于开口椭圆环的高效超宽带极化旋转超表面. 物理学报, 2015, 64(17): 178101. doi: 10.7498/aps.64.178101
    [19] 陈吴玉婷, 韩鹏昱, Kuo Mei-Ling, Lin Shawn-Yu, 张希成. 具有缓变折射率的太赫兹宽带增透器件. 物理学报, 2012, 61(8): 088401. doi: 10.7498/aps.61.088401
    [20] 张庆斌, 兰鹏飞, 洪伟毅, 廖青, 杨振宇, 陆培祥. 控制场对宽带超连续谱产生的影响. 物理学报, 2009, 58(7): 4908-4913. doi: 10.7498/aps.58.4908
计量
  • 文章访问数:  6450
  • PDF下载量:  294
  • 被引次数: 0
出版历程
  • 收稿日期:  2019-07-10
  • 修回日期:  2019-08-28
  • 上网日期:  2019-11-27
  • 刊出日期:  2019-12-05

基于反射超表面产生太赫兹涡旋波束

  • 1. 首都师范大学物理系, 北京 100048
  • 2. 北京市成像理论与技术高精尖创新中心, 北京 100048
  • 3. 太赫兹光电子学教育部重点实验室, 北京 100048
  • 通信作者: 赵国忠, guozhong-zhao@126.com
    基金项目: 国家自然科学基金(批准号: 5307625130)资助的课题

摘要: 具有螺旋波前的电磁波是携带轨道角动量的涡旋波束, 涡旋波束存在的相位奇点, 使其在微粒操控和通讯等领域有特殊的应用. 本文提出了一种基于反射型超表面的太赫兹宽带涡旋波束产生器, 该器件由超表面-电介质-金属三层结构构成, 顶层为两个正交I形金属结构单元组成的超表面, 中间层是聚酰亚胺介质, 最底层为金属. 通过对超表面单元结构参数的优化设计, 可以实现在不同旋转角度下反射波振幅高达90%以上, 同时反射波相位随旋转角线性变化的目的. 进一步利用这些单元结构按照Pancharatnam-Berry相位原理进行超表面布阵, 在0.8—1.4 THz频率范围内, 可以将圆偏振太赫兹波束转换为具有轨道角动量的涡旋波束, 这一器件的工作带宽相对较宽, 结构简单, 转换效率高, 在太赫兹涡旋波束产生方面具有潜在的应用价值.

English Abstract

    • 太赫兹波具有穿透性、吸水性、高分辨率等诸多优越特性, 在光谱[1-3]、成像[4-6]、无损检测[7,8]等领域已经展示了重要的应用前景. 携带轨道角动量的涡旋光束, 具有由exp(ilΦ)描述的螺旋相位结构, 其中Φ是方位角, l是拓扑荷数, 它在不增加带宽的情况下通过对光通讯信息进行编码, 可以极大地扩展通信容量, 在光通信[9]领域具有广阔的应用前景. 传统上人们利用螺旋相位板[10]、全息衍射光栅[11]以及天线阵列[12,13]等方法来产生涡旋波束, 然而这些方法存在体积大, 不利于器件集成化、成本高等问题. 作为一种人工合成材料, 超表面具有自然材料所不具备的性质, 因而受到广泛关注. 人们通过在金属表面设计微结构来实现对电磁波的相位和振幅的调制. 近年来应用超表面实现的功能器件有偏振转换器[14-16]、高方向辐射天线[17]、高增益透镜阵列[18,19]以及涡旋光束产生器[20,21]等.

      利用超表面产生涡旋波束不仅解决了体积大、不利于集成的问题, 还大大降低了制造成本. 近年来, 在可见光和微波波段基于超表面产生涡旋波束的研究较多. 2011年, 哈佛大学Yu等[22]在光学薄界面上设计V字形天线单元, 通过选择八个相位间隔45°的V字形单元结构排布超表面, 实现了在红外波段产生拓扑荷数为1的涡旋光. 2014年, Zhang等[23]在微波频段内基于开口谐振环提出了单层透射超表面, 圆偏振光入射超表面产生了拓扑荷数为1的涡旋光, 但是透过率仅有24.7%. 2017年, Lou等[24]基于三层超表面结构在10.1—10.9 GHz产生了涡旋光, 但是存在制备复杂等问题.

      目前为止, 关于涡旋波束的产生, 普遍存在效率不高、制备困难等问题, 且大部分工作是在微波、红外频段, 在太赫兹频段产生涡旋波束的研究甚少. 本文提出了一种新型反射式超表面结构, 实现了在太赫兹频段产生高效率涡旋光的目的.

    • 为了实现利用超表面产生涡旋光的目的, 需要设计能够独立地控制入射和反射波的单元结构. 根据Pancharatnam-Berry (P-B)相位原理, 超表面上可以进行单元结构布阵, 各个单元应该能够获得振幅相近、相位可以独立调节的反射波. 我们提出了一个由正交I形金属结构单元组成的结构, 可以满足上述控制反射波振幅和相位的设计要求.

      当一束平面波垂直照射到反射型超表面单元结构, 单元结构以其中心为旋转中心, 波束传播方向为旋转轴, 逆时针旋转α角时, 如图1所示, 反射矩阵可以写作[25]:

      图  1  超表面单元结构示意图 (a)顶视图; (b)侧视图

      Figure 1.  Schematic diagram of the unit cell of metasurface: (a) Top view; (b) side view.

      ${ R}_\alpha^{XY} \!=\! {\left[\!\!\!{\begin{array}{*{20}{c}} {\cos \alpha }\!\!&\!\!{ - \sin \alpha } \\ {\sin \alpha }\!\!&\!\!{\cos \alpha } \end{array}}\!\!\!\right]^{ - 1}}\left[\!\!\!{\begin{array}{*{20}{c}} {{r_{xx}}}\!\!&\!\!{{r_{xy}}} \\ {{r_{yx}}}\!\!&\!\!{{r_{yy}}} \end{array}}\!\!\!\right]\left[\!\!\!{\begin{array}{*{20}{c}} {\cos \alpha }\!\!&\!\!{ - \sin \alpha } \\ {\sin \alpha }\!\!&\!\!{\cos \alpha } \end{array}}\!\!\!\right],$

      其中${r_{xx}}$表示沿x轴正方向的线偏振光入射时, 沿x轴负方向出射的反射系数, 同理${r_{xy}}$, ${r_{yx}}$, ${r_{yy}}$有类似的意义.

      对于圆偏振光入射到超表面, 反射矩阵表示为

      ${ R}_\alpha ^{{\rm{LR}}} = \frac{1}{2}{\left[\!\!{\begin{array}{*{20}{c}} 1&{ - {\rm{j}}} \\ 1&{\rm{j}} \end{array}}\!\! \right]}{{R}}_\alpha ^{XY}{\left[\!\!{\begin{array}{*{20}{c}} 1&{ - {\rm{j}}} \\ 1&{\rm{j}} \end{array}}\!\!\right]^{ - 1}},$

      其中参数α为单元结构相对于其中心的转角, LR代表入射光的偏振态, L表示左旋偏振, R表示右旋偏振, ${{{R}}_\alpha }^{{\rm{LR}}}$表示圆偏振光入射时的反射矩阵.

      ${r_{{\rm{LL}}}} = \frac{1}{2}\left[ {({r_{xx}} - {r_{yy}}) - {\rm{j}}({r_{xy}} + {r_{yx}})} \right]{{\rm{e}}^{ - {\rm{j}}2\alpha }},\tag{3a}$

      ${r_{{\rm{RR}}}} = \frac{1}{2}\left[ {({r_{xx}} - {r_{yy}}) + {\rm{j}}({r_{xy}} + {r_{yx}})} \right]{{\rm{e}}^{{\rm{j}}2\alpha }},\tag{3b}$

      ${r_{{\rm{RL}}}} = \frac{1}{2}\left[ {({r_{xx}} + {r_{yy}}) + {\rm{j}}({r_{xy}} - {r_{yx}})} \right],\qquad \tag{3c}$

      ${r_{{\rm{LR}}}} = \frac{1}{2}\left[ {({r_{xx}} + {r_{yy}}) - {\rm{j}}({r_{xy}} - {r_{yx}})} \right].\quad \tag{3d}$

      从(3)式可以看出, 只有${r_{{\rm{LL}}}}$, ${r_{{\rm{RR}}}}$携带轨道角动量, 当单元结构满足$\left| {{r_{xx}}} \right| \approx \left| {{r_{yy}}} \right| \approx 1$或者$\left| {{r_{xy}}} \right| \approx $$ \left| {{r_{yx}}} \right| \approx 0$且相位差为π时, 两个交叉极化之间不存在相位差, 因此当单元结构旋转角α时, 相位φ = ± 2α, “+”和“–”表示右旋和左旋入射光.

      基于上述理论可知, 在设计P-B超表面时, 通过旋转单元结构找出具有不同相位的单元, 使得单元结构具有相同且近似为1的反射振幅, 同时具有180°的相位差, 可实现在圆偏振光入射时, 反射超表面可以使入射光成为携带轨道角动量的涡旋光. 但是, 目前涡旋光产生器件存在带宽窄的问题, 为了增加超表面的工作带宽, 需要确保在不同频率处具有相同的斜率:

      $ {{\partial \varphi ({\rm{R}},{\rm{R}})}}/{{\partial {f_i}}} \approx {{\partial \varphi ({\rm{R}},{\rm{R}})}}/{{\partial {f_j}}}, $

      其中$\partial \varphi ({\rm{R}}, {\rm{R}})$表示右旋入射光的反射相位变化, ${f_i}$${f_j}$表示任意两处不同的频率, 为了保证超表面有合适的工作带宽, 需要保证$({f_j} - {f_i})$足够大.

    • 在超表面结构设计过程中, 首先考虑的主要问题为如何设计单元结构、单元结构参数如何确定、衬底如何选择等. 我们设计了一个由超表面-电介质-金属三层结构组成的反射型超表面器件, 如图1所示. 超表面由三层组成. 其中顶层为正交I型金属结构, 厚度为d1; 中间层为介质层, 材料为聚酰亚胺, 厚度为d2; 底层为金膜, 厚度为d1. 考虑实验室材料、各个仿真结果以及样品制备, 最终得出单元结构具体参数如下: 周期p = 90 μm, 金属线宽w = 3 μm, l1 = 60 μm, l2 = 28 μm, s = 27 μm, 金属层厚度d1 = 0.2 μm, 衬底厚度d2 = 35 μm. 入射光为偏振方向沿xy轴方向的线偏光, 波矢沿z轴方向.

    • 利用CST microwave studio (CST)软件, 以周期性边界条件对单元结构反射振幅和反射相位进行仿真, 仿真频率范围设为0.5—2.0 THz, 结果如图2. 在线性x偏振和线性y偏振入射下, 在0.8—1.4 THz频率范围内共极化反射幅值均大于0.9且相位差接近180°, 在0.93 THz时共极化相位差刚好为180°且反射幅值近似为1. 因此可以通过改变转角实现对反射波振幅和相位的调制.

      图  2  线偏振光入射单元结构产生的反射相位和振幅谱 (a)相位谱; (b)振幅谱

      Figure 2.  Reflected phase and amplitude spectrum produced by linearly polarized light incident unit call: (a) Phase; (b) amplitude.

      以0.85 THz为例, 对单元结构反射振幅和相位进行数值仿真, 可得出不同转角与相位及振幅的一一对应关系, 如图3. 不同转角下单元结构的反射振幅近似相等, 转角与相位呈现线性关系, 与上述P-B相位原理符合得很好. 因此可得到不同相位值对应的不同转角的单元结构.

      图  3  仿真得到的不同转角的相位和反射振幅

      Figure 3.  The phase and reflection amplitude of different corners simulated.

      入射光为圆偏振光时, 对参数α进行扫描, 扫描范围为0°—157.5°, 扫描步长设置为22.5°, 仿真结果如图4所示.

      图  4  不同转角单元结构的反射相位和振幅谱 (a) LCP入射的相位谱; (b) RCP入射的相位谱; (c) LCP入射的振幅谱; (d) RCP入射的振幅谱

      Figure 4.  Reflective phase and amplitude spectra of unit cell structure under different rotation angle: (a) Phase spectra at LCP incident; (b) phase spectra at RCP incident; (c) amplitude spectra at LCP incident; (d) amplitude spectra at RCP incident.

      图4(a)图4(b)分别为左旋圆偏光(LCP)和右旋圆偏光(RCP)入射到超表面单元结构时, 转角α在0°—157.5°范围内以22.5°变化时得到的相位谱. 从图4(a)图4(b)可以看出, 相位以二倍转角(45°)间隔增大(减小), 不同转角对应的相位覆盖了2π. 图4(c)图4(d)分别为左旋圆偏光和右旋圆偏光入射到超表面单元结构时, 转角α在0°—157.5°范围内以22.5°的变化改变时得到的振幅谱, 可以看出, 反射振幅几乎相同且高达0.9, 满足控制反射波振幅和相位的设计要求.

    • 基于文中提出的反射型超表面单元, 可以通过相位分布设计产生不同拓扑荷数的涡旋波束. 以超表面中心为原点, 超表面每个位置(x, y)的相位需要满足

      $ \varphi (x,y) = l \cdot \arctan \left({y}/{x}\right), $

      其中$\varphi (x, y)$表示超表面$(x, y)$位置需要满足的相位, l表示拓扑荷数. 通过改变l的值可以产生不同拓扑荷数的涡旋波束.

      本文设计了拓扑荷数l = 1和l = 2的反射型超表面. 根据(4)式确定的相位-位置关系以及单元结构相位是转角α的2倍关系, 排布了拓扑荷数为1和2的超表面结构, 如图5所示.

      图  5  两种用于产生拓扑荷数分别为 (a) l = 1和(b) l = 2的涡旋波束超表面结构

      Figure 5.  Two kinds of metasurface structures for generating the vortex beam with topological charge numbers (a) l = 1 and (b) l = 2.

      l = 1和l = 2的两种超表面反射波束特性进行仿真. 不失一般性, 假设入射波束为左旋圆偏振波, 半径r1, r2, r3分别设为120, 200, 300 μm, 其中, 半径的选择考虑到各个单元结构之间不能够重叠, 以及实验测量中的尺寸、加工难度的因素. 以0.93 THz为例, 图6给出了反射圆偏振波的振幅和相位分布.

      图  6  超表面产生l = 1和l = 2的涡旋波束反射振幅和相位分布 LCP入射l = 1超表面的(a)振幅分布和(b)相位分布; RCP入射l = 2超表面的(c)振幅分布和(d)相位分布

      Figure 6.  Reflective amplitude and phase distributions of vortex beams with l = 1 and l = 2 generated by metasurface. LCP incident l = 1 metasurface: (a) amplitude distribution and (b) phase distribution; RCP incident l = 2 metasurface: (c) amplitude distribution and (d) phase distribution.

      图6(a)图6(b)是0.93 THz处左旋圆偏振波垂直入射到拓扑荷数为1的反射型超表面上, 反射光场距离超表面500 μm处的反射振幅和相位分布. 从图中可以看出振幅呈现中心为暗环且中心强度始终为零的分布, 相位围绕一周改变2π. 图6(c)图6(d)是0.93 THz的右旋圆偏振波垂直入射到拓扑荷数为2的反射型超表面上, 反射光场距离超表面500 μm处的反射振幅和相位分布. 可以看出振幅分布呈现中心场强为零, 外环场强相对较大, 相位覆盖4π. 表明了在圆偏振波束入射到超表面可以产生拓扑荷数为1和2的涡旋波束.

      图6(a)图6(c)可以看出拓扑荷数为1的涡旋光束的中心暗环半径相对于拓扑荷数为2的涡旋光中心暗环半径较小, 这是因为随着拓扑荷数的增大, 超表面相邻的相位梯度也随之增大导致的.

      在0.8—1.4 THz频段范围内, 圆偏光入射到超表面时均可产生涡旋光. 图7给出了在0.8 THz和1.4 THz两个频率下, 左旋圆偏振光入射在l = 1的超表面时, 反射波束在传播一定距离处的相位分布.

      图  7  不同频率下l = 1和l = 3超表面产生的反射涡旋波束振幅、相位分布 l = 1超表面(a) 0.8 THz频率下振幅分布, (b) 0.8 THz频率下相位分布, (c) 1.4 THz频率下振幅分布, (d) 1.4 THz频率下相位分布; l = 3超表面(e) 1.1 THz频率下振幅分布, (f) 1.1 THz频率下相位分布

      Figure 7.  The amplitude and phase distribution of reflective vortex beam generated by the LCP incident l = 1 and l = 3 metasurface at different frequencies. l = 1: (a) amplitude distribution at 0.8 THz, (b) phase distribution at 0.8 THz, (c) amplitude distribution at 1.4 THz, (d) phase distribution at 1.4 THz. l = 3: (e) amplitude distribution at 1.1 THz, (f) phase distribution at 1.1 THz.

      图7(a)图7(c)分别为左旋圆偏光在0.8, 1.4 THz下入射到l = 1的超表面上产生的反射光场距离超表面700 μm处的振幅分布, 可以看出, 振幅分布相对均匀, 呈现中间暗环, 外环场强相对较大的甜甜圈分布, 图7(b)图7(d)分别为左旋圆偏光在0.8, 1.4 THz下入射到l = 1的超表面上产生的反射光场距离超表面700 μm处的相位分布, 可以看出相位变化了2π, 综上符合涡旋光的相位性质, 并且反射率达到86%.

      此外, 为证明设计的超表面可以产生更高阶的涡旋光束, 设计并仿真了拓扑荷数为3的超表面振幅和相位图像. 以1.1 THz频率下为例, 图7(e)图7(f)给出了反射场距离超表面300 μm处的反射振幅和相位分布. 从图中可以看出振幅呈现中间暗环的甜甜圈形状, 相位覆盖了6π, 仍然符合涡旋光束的特性.

      结合图6图7的仿真结果, 在0.8, 0.93, 1.4 THz频率下, 圆偏振光垂直入射到l = 1和l = 2的超表面时, 反射波束具有螺旋相位, 可产生中心强度为零的涡旋光. 证明提出的反射型超表面可以在0.8—1.4 THz频率范围内产生涡旋光, 并且可以产生更高阶的涡旋光束.

    • 基于相位突变的结构单元, 利用P-B相位原理设计了一个简单而高效的反射型超表面, 该超表面由正交I型结构-介质-金属三层结构组成. 仿真结果表明, 我们设计的超表面结构可以在0.8—1.4 THz频段圆偏振光正入射条件下, 产生携带轨道角动量的涡旋波束, 反射率达到86%, 最高可达92%. 我们设计的超表面具有效率高、宽带并且结构简单等优点, 对基于超表面的太赫兹涡旋波束产生具有参考价值.

参考文献 (25)

目录

    /

    返回文章
    返回