高级检索

ISSN1001-3806CN51-1125/TN 网站地图

留言板

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

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

基于π型结构双折射超表面的设计与应用

韩晓晓 童元伟

引用本文:
Citation:

基于π型结构双折射超表面的设计与应用

    作者简介: 韩晓晓(1990-), 女, 硕士研究生, 现主要从事超表面材料的研究.
    通讯作者: 童元伟, tyw0991@usst.edu.cn
  • 基金项目:

    国家自然科学基金资助项目 11504236

    上海理工大学科技发展资助项目 16KJFZ078

  • 中图分类号: O436

Design and application of birefringent metasurface based on π-shape structure

    Corresponding author: TONG Yuanwei, tyw0991@usst.edu.cn
  • CLC number: O436

  • 摘要: 为了实现结构简单、透射率高的双折射超表面, 采用广义薄板跃迁条件分析了该超表面结构与其周围入射场、反射场和透射场的关系, 并利用介质空间域的表面极化率、磁化率等描述了相应超表面的等效特性, 设计出一种基于π型金属结构单元的双折射超表面。通过将具有梯度透射相位的7个单元按照顺序排列, 形成具有对垂直入射x, y极化电磁波双折射性能的超表面。结果表明, 在波束折射和偏振分束超表面中, 损失均低于-8dB; λ/4波片中达到全透射; 所设计的双折射超表面对垂直入射的电磁波具有高透射特性, 并且能够分离x极化电磁波和y极化电磁波的波束, 实现双折射。该研究结果对高性能超表面的设计与实现具有一定的指导意义。
  • Figure 1.  Schematic diagram of electromagnetic waves refracting through matasurface

    Figure 2.  Polarization and magnetization coefficients of metasurface for normal incident plane wave and refraction angle of 45°

    a, c—χe, e and χm, m for equation (8) b, d—χe, e and χm, m for equation (11)

    Figure 3.  Scattering particles of metasurface

    a—generic unit cell consisted by cascaded sheet admittances b—工-shaped metallic particle c—jerusalem-cross-shaped metallic particle d—π-shaped metallic particle e—generic unit cell perspective view with dielectric substrate

    Figure 4.  Quarter wave-plate

    Figure 5.  Transmission coefficient S21 of quarter wave-plate

    a—transmission magnitude of quarter wave-plate b—transmission phase of quarter wave-plate

    Figure 6.  Design shapes of outer and inner layers of beam refractive metasurface metalattices

    a—middle layers b—outer layers

    Figure 7.  Simulation results of beam refraction metasurface

    a—transmission magnitude of transmission array b—electric field of x-polarized wave with refraction angle of 45°

    Figure 8.  Metasurface of the polarized beam splitting

    Figure 9.  Supercell with 7×7 unit cells

    a—the simulated metasurface of polarization beam splitting b—transmission magnitude of PBS c—φx=45° d—φy=135°

    Table 1.  Dimensions of quarter wave-plate

    parameters Ay/mm Bx/mm Lx/mm Sx/mm Wy/mm
    outer layers 0.75 5.5 0.6 1.55 4.25
    middle layer 0.75 5 1 1.25 4.25
    下载: 导出CSV

    Table 2.  Transmission phase shifts for x and y polarization of 7 unit cells forming the metacell

    No. 1 2 3 4 5 6 7
    φx/(°) 0 51.5 103 154.5 206 257.5 309
    φy/(°) 309 257.5 206 154.5 103 51.5 0
    下载: 导出CSV

    Table 3.  Design dimensions of beam refraction metasurface inner layer of π-shaped metal

    No. Ay/mm Bx/mm Lx/mm Sx/mm Wy/mm
    1 0.75 5 1 1.25 4.25
    2 0.75 5 1 1.25 4.25
    3 0.75 4.5 1 1.25 4.25
    4 0.5 4 1 0.5 4.5
    5 0.25 3.75 1.5 1.25 4.5
    6 0.25 3.75 1.5 1.5 3.5
    7 0.75 4.25 0.75 1 4.5
    下载: 导出CSV

    Table 4.  Design dimensions of beam refraction metasurface outer layer of π-shaped metal

    No. Ay/mm Bx/mm Lx/mm Sx/mm Wy/mm
    1 0.75 5.5 1.025 1.175 4.25
    2 0.75 5.5 0.6 1.55 4.25
    3 0.75 5.5 0.45 1.45 4
    4 0.5 5 1 0.5 4
    5 0.75 4.25 1.25 0.75 3.5
    6 0.5 4.25 1 0.25 3.25
    7 0.25 4.25 2 1 3.75
    下载: 导出CSV
  • [1]

    VESELAGO V G. The electrodynamics of substances with simultaneously negative values of ε and μ[J]. Soviet Physics Uspekhi, 1968, 10(4):509-514. doi: 10.1070/PU1968v010n04ABEH003699
    [2]

    GAO X, HAN X, CAO W P, et al. Ultrawideband and high-efficiency linear polarization converter based on double V-shaped metasurface[J]. IEEE Transactions on Antennas & Propagation, 2015, 63(8):3522-3530. 
    [3]

    LI Y, ZHANG J, QU S, et al. Achieving wide-band linear-to-circular polarization conversion using ultra-thin bi-layered metasurfaces[J]. Journal of Applied Physics, 2015, 117(4):011129. 
    [4]

    LIU X B, ZHANG J S, LI W, et al. Three-band polarization converter based on reflective metasurface[J]. IEEE Antennas & Wireless Propagation Letters, 2016, 16:924-927. 
    [5]

    ACHOURI K, SALEM M A, CALOZ C. General metasurface synthesis based on susceptibility tensors[J]. IEEE Transactions on Antennas and Propagation, 2015, 63(7): 2977-2991. doi: 10.1109/TAP.2015.2423700
    [6]

    NIEMI T, KARILAINEN A O, TRETYAKOY S A. Synthesis of polarization transformers[J]. IEEE Transactions on Antennas and Pro-pagation, 2013, 61(6): 3102-3111. doi: 10.1109/TAP.2013.2252136
    [7]

    KODERA T, SOUNAS D L, CALOZ C. Artificial Faraday rotation using a ring metamaterial structure without static magnetic field[J]. Applied Physics Letters, 2011, 99(3): 031114. doi: 10.1063/1.3615688
    [8]

    YU N, CAPASSO F. Flat optics with designer metasurfaces[J]. Nature Materials, 2014, 13(2):139-50. doi: 10.1038/nmat3839
    [9]

    SHI H Y, LI J X, ZHANG A X. Gradient metasurface with both polarization-controlled directional surface wave coupling and anomalous reflection[J]. IEEE Antennas & Wireless Propagation Letters, 2015, 14:104-107. 
    [10]

    SUN S, YANG K Y, WANG C M, et al. High-efficiency broadband anomalous reflection by gradient meta-surfaces[J]. Nano Letters, 2012, 12(12):6223. doi: 10.1021/nl3032668
    [11]

    SALEM M A, ACHOURI K, CALOZ C. Metasurface synthesis for time-harmonic waves: Exact spectral and spatial methods[J]. Progress in Electromagnetics Research, 2014, 149: 205-216. doi: 10.2528/PIER14100505
    [12]

    SALEM M A, CALOZ C. Manipulating light at distance by a metasurface using momentum transformation[J]. Optics Express, 2014, 22(12): 14530-14543. doi: 10.1364/OE.22.014530
    [13]

    CHEN P Y, ALU A. Subwavelength imaging using phaseconjugating nonlinear nanoantenna arrays[J]. Nano Letters, 2011, 11(12): 5514-5518. doi: 10.1021/nl203354b
    [14]

    ROY T, ROGERS E T, ZHELUDEV N I. Sub-wavelength focusing meta-lens[J]. Optics Express, 2013, 21(6): 7577-7582. doi: 10.1364/OE.21.007577
    [15]

    ACHOURI K, LAVIGNE G, SALEM M A, et al. Metasurface spatial processor for electromagnetic remote control[J]. IEEE Transactions on Antennas and Propagation, 2016, 64(5): 1759-1767. doi: 10.1109/TAP.2016.2537369
    [16]

    SOUNAS D L, KODERA T, CALOZ C. Electromagnetic modeling of a magnetless nonreciprocal gyrotropic metasurface[J]. IEEE Transactions on Antennas and Propagation, 2013, 61(1): 221-231. doi: 10.1109/TAP.2012.2214997
    [17]

    HADAD Y, SOUNAS D L, ALU A. Space-time gradient metasurfaces[J]. Physical Review, 2015, B92(10): 100304. 
    [18]

    KHAN B A, GUPTA S, CALOZ C. Spatial nonreciprocal and nongyrotropic structure[C]//International Symposium on Antennas and Propagation (ISAP). New York, USA: IEEE, 2015: 1-4.
    [19]

    ASADCHY V, ALBOOYEH M, TCVETKOVA S, et al. Metasurfa-ces for perfect and full control of refraction and reflection[C]//10th International Congress on Advanced Electromagnetic Materials in Microwaves and Optics (METAMATERIALS).New York, USA: IEEE, 2016: 364-366.
    [20]

    ZHU B O, FENG Y. Passive Metasurface for reflectionless and arbi-tary control of electromagnetic wave transmission[J]. IEEE Transactions on Antennas & Propagation, 2015, 63(12): 5500-5511. 
    [21]

    PFEIFFER C, GRBIC A. Bianisotropic metasurfaces for optimal polarization control: Analysis and synthesis[J]. Physical Review A-pplied, 2014, 2(4): 044011. 
  • [1] 冯可豪李白萍蔡艺军周远国 . 基于AdaBelief残差神经网络的超材料结构逆设计. 激光技术, 2022, 46(3): 307-311. doi: 10.7510/jgjs.issn.1001-3806.2022.03.003
    [2] 高丹盛荔金亓吴轲王梦樱陶音孔勇韩华 . 入射角度不敏感型石英光双折射滤波器. 激光技术, 2015, 39(3): 329-332. doi: 10.7510/jgjs.issn.1001-3806.2015.03.010
    [3] 李文煜章海锋刘婷马宇 . 一种波束扫描固态等离子体超表面的设计. 激光技术, 2018, 42(6): 822-826. doi: 10.7510/jgjs.issn.1001-3806.2018.06.018
    [4] 张雪妍郁步昭王吉明吴彤赫崇君刘友文路元刚 . 基于几何相位超表面的Ince-Gaussian矢量涡旋光场聚焦. 激光技术, 2022, 46(1): 85-93. doi: 10.7510/jgjs.issn.1001-3806.2022.01.008
    [5] 郑盛梅江孝伟江达飞王琳 . 基于黑磷的双频带超材料吸收体及其传感特性. 激光技术, 2023, 47(6): 846-853. doi: 10.7510/jgjs.issn.1001-3806.2023.06.017
    [6] 蒋孝鑫王吉明黄鑫赫崇君刘友文 . 基于狭槽天线超表面的复振幅光瞳滤波器. 激光技术, 2017, 41(6): 807-811. doi: 10.7510/jgjs.issn.1001-3806.2017.06.008
    [7] 李瑶莫伟成杨振刚刘劲松王可嘉 . 利用超表面天线阵列产生太赫兹涡旋光束. 激光技术, 2017, 41(5): 644-648. doi: 10.7510/jgjs.issn.1001-3806.2017.05.005
    [8] 张剑李国华 . 一种新的Lyot型双折射调谐滤光片. 激光技术, 2005, 29(6): 639-640,648.
    [9] 王建国贺洪波邵建达范正修 . 光波在具有柱状结构的双折射薄膜中的传输特性. 激光技术, 2005, 29(6): 601-603,607.
    [10] 苑秋红谢康韩艳芬 . 含负折射率材料的1维光子晶体掺杂缺陷模研究. 激光技术, 2010, 34(2): 232-235. doi: 10.3969/j.issn.1001-3806.2010.02.024
    [11] 赵培李国华 . 90°转向复合延迟器件的研制. 激光技术, 2004, 28(4): 424-426.
    [12] 胡晨曦王吉明吴彤赫崇君顾晓蓉刘友文 . 基于超表面的宽波带光束聚焦研究. 激光技术, 2018, 42(5): 681-686. doi: 10.7510/jgjs.issn.1001-3806.2018.05.018
    [13] 张大勇刘仓理李剑峰罗飞骆永全沈志学刘海涛李建明 . 楔型液晶盒对入射激光偏振态的影响. 激光技术, 2006, 30(6): 636-637,660.
    [14] 张浩马宇章海锋杨靖刘佳轩 . 一种带宽展宽的等离子体超材料吸波体的设计. 激光技术, 2019, 43(2): 256-262. doi: 10.7510/jgjs.issn.1001-3806.2019.02.020
    [15] 王鹏飞贺风艳刘建军井绪峰洪治 . 基于连续谱束缚态的高Q太赫兹全介质超表面. 激光技术, 2022, 46(5): 630-635. doi: 10.7510/jgjs.issn.1001-3806.2022.05.008
    [16] 邓之兵李岩韩艳梅张书练 . 用双折射激光频率分裂原理实现的静态与动态测量方法. 激光技术, 1999, 23(1): 29-33.
    [17] 张剑李国华郝殿中 . 离子注入铌酸锂平面波导双折射特性的研究. 激光技术, 2006, 30(4): 360-362,365.
    [18] 甄海龙 . 一种金填充高双折射光子晶体光纤偏振滤波器. 激光技术, 2016, 40(1): 1-4. doi: 10.7510/jgjs.issn.1001-3806.2016.01.001
    [19] 夏源谢卉孙莉萍胡强高 . 超窄带宽可调光滤波器的研究. 激光技术, 2013, 37(4): 493-497. doi: 10.7510/jgjs.issn.1001-3806.2013.04.017
    [20] 王海峰吴福全郝殿中王涛张姗王晶 . 光学胶折射率对渥拉斯顿棱镜透射性能的影响. 激光技术, 2007, 31(6): 668-670.
  • 加载中
图(9) / 表(4)
计量
  • 文章访问数:  6960
  • HTML全文浏览量:  5372
  • PDF下载量:  27
  • 被引次数: 0
出版历程
  • 收稿日期:  2019-03-19
  • 录用日期:  2019-04-15
  • 刊出日期:  2020-01-25

基于π型结构双折射超表面的设计与应用

    通讯作者: 童元伟, tyw0991@usst.edu.cn
    作者简介: 韩晓晓(1990-), 女, 硕士研究生, 现主要从事超表面材料的研究
  • 上海理工大学 理学院, 上海 200093
基金项目:  国家自然科学基金资助项目 11504236上海理工大学科技发展资助项目 16KJFZ078

摘要: 为了实现结构简单、透射率高的双折射超表面, 采用广义薄板跃迁条件分析了该超表面结构与其周围入射场、反射场和透射场的关系, 并利用介质空间域的表面极化率、磁化率等描述了相应超表面的等效特性, 设计出一种基于π型金属结构单元的双折射超表面。通过将具有梯度透射相位的7个单元按照顺序排列, 形成具有对垂直入射x, y极化电磁波双折射性能的超表面。结果表明, 在波束折射和偏振分束超表面中, 损失均低于-8dB; λ/4波片中达到全透射; 所设计的双折射超表面对垂直入射的电磁波具有高透射特性, 并且能够分离x极化电磁波和y极化电磁波的波束, 实现双折射。该研究结果对高性能超表面的设计与实现具有一定的指导意义。

English Abstract

    • 新型人工电磁材料是一种由周期或者非周期的人工单元结构排列构成的复合电磁材料[1],其单元尺寸一般远小于工作波长,主要通过调节单元结构满足所需的电磁参量,进而实现媒质对电磁波的灵活调控。超材料通常会有一些超常的物理特性如负折射率、负介电常数等。在过去的20多年里,新型人工电磁材料呈现出飞速发展的趋势,在各个领域都得到了深入的研究。但是这些设计仍然存在体积过大以及工作带宽太窄的问题。人工电磁超表面[2-4]可代替这些体积较大的超材料,设计出更加轻薄、微小和高效的器件,这些优势使得它在工程领域得到了广泛的应用。

      人工电磁超表面[5]是由亚波长尺度的散射粒子按特定顺序组成的2维结构的超材料,它主要利用电磁波在超表面面上的相位突变来控制透射波或反射波,与传统的光学器件依靠传播路径来积累相位相比,其厚度一般小于操控电磁波波长的1/10,可以有效避免入射波的能量损失,具有损失小、易集成等优点。电磁波的传播在超表面两侧遵循广义斯涅耳定律,可以灵活控制电磁波传播方向、极化方式、透射反射强度等。由各向异性单元结构构成的超表面在极化调制[6-10]方面的应用主要包括实现电磁波的圆极化转换、宽带极化转换、多频点极化转换等。超表面还可以对电磁波的波束方向[11-12]进行灵活调控,如双折射。目前已经被应用到几种主要光学元件中,如晶体偏振器和波片等,本文中研究的双折射是指超表面对入射的x偏振和y偏振电磁波的透射系数(包括幅值和相位)进行独立、同步控制。除此之外,超表面在高分辨率成像[13-14]、电磁干涉仪[15]和非互易超表面等[16-18]方面也取得了显著的发展。

      在传统的超表面设计中,“工”字型单元结构是一种典型的各向异性结构单元,由于它的调节对象比较单一,只能调节与“工”字型中间金属线平行方向上电场的相位,所以一般只能用来实现线极化到圆极化的转换或者圆极化波的异常反射和透射,限制了它在极化调控中的应用。因此,研究人员对“工”字型单元结构进行了改进,提出了两个彼此垂直的“工”字型组成的正交“工”字型单元结构。它继承了“工”字型单元结构的优点,并且能够对两个正交方向电场的相位进行独立调节。但这种各向异性结构单元通常至少含有8个参量,可调参量的增多也加大了设计难度。本文中提出的π型金属结构单元仅包含5个参量,既可以实现对电磁波极化的灵活调控,也可以有效减少散射粒子的设计参量。通过λ/4波片、波束折射超表面和偏振分束超表面的设计与仿真,实现了极化转换和xy极化波的偏振分离,在覆盖全相位的情况下,可达到高效传输的目的,对于未来高效、轻薄光学器件的设计具有重要意义。

    • 根据Idemen推导出的通过广义薄板跃迁条件分析超表面结构与其周围入射场、反射场和透射场的关系,采用介质空间域的表面极化率、磁化率等描述超表面的等效特性。超表面的等效表面极化率、磁化率与超表面两侧的电磁场不连续性相关,超表面的等效极化、磁化张量χe, eχm, m与超表面两侧的电磁场具有如下的关系:

      $ \left\{ {\begin{array}{*{20}{l}} {\hat z \times \Delta \mathit{\boldsymbol{H}} = {\rm{j}}\omega {\varepsilon _0}{\mathit{\boldsymbol{\chi }}_{e,{\rm{e}}}} \cdot {\mathit{\boldsymbol{E}}_{{\rm{av}}}}}\\ {\Delta \mathit{\boldsymbol{E}} \times \hat z = {\rm{j}}\omega {\mu _0}{\mathit{\boldsymbol{\chi }}_{{\rm{m}},{\rm{m}}}} \cdot {\mathit{\boldsymbol{H}}_{{\rm{av}}}}} \end{array}} \right. $

      (1)

      式中,χe, e的第1个下标表示电场响应,第2个下标表示电场激励,χm, m的第1个下标表示磁场响应,第2个下标表示磁场激励,ΔH=Ht-(Hi+Hr),Eav=(Et+Ei+Er)/2,ΔH和ΔE分别代表超表面两侧的磁场和电场的差,HavEav分别代表超表面两侧的平均磁场和平均电场,Ei, Hi分别表示入射场的电场和磁场,Er, Hr分别表示反射场的电场和磁场,Et, Ht分别表示透射场的电场和磁场,ω表示单位时间内的场量相位的变化, $\hat z$表示方向向量,ε0表示真空介电常数,μ0表示真空磁导率。将厚度极薄的超表面(其厚度远小于入射波波长,可视为0)放置在x-y平面,取z=0,电磁波经过超表面异常折射的示意图,如图 1所示。

      Figure 1.  Schematic diagram of electromagnetic waves refracting through matasurface

      图 1中,假设电磁波入射场的入射角ϕi=θi、反射场的反射角ϕr=θr、透射场的折射角ϕt=θt,且入射场幅值为1,则对应的p偏振的电场和磁场分别描述为(省略与时间相关相位项ωt):

      $ \left\{ \begin{array}{l} {\mathit{\boldsymbol{E}}_{\rm{i}}} = \left( {\cos {\theta _{\rm{i}}}\mathit{\boldsymbol{\hat x}} + \sin {\theta _{\rm{i}}}\mathit{\boldsymbol{\hat z}}} \right)\exp \left( {{\rm{j}}{k_0}x\sin {\theta _{\rm{i}}} - {\rm{j}}{k_0}z\cos {\theta _{\rm{i}}}} \right)\\ {\mathit{\boldsymbol{E}}_{\rm{t}}} = \mathit{\boldsymbol{T}}\left( {\cos {\theta _{\rm{t}}}\mathit{\boldsymbol{\hat x}} + \sin {\theta _{\rm{t}}}\mathit{\boldsymbol{\hat z}}} \right)\exp \left( {{\rm{j}}{k_0}x\sin {\theta _{\rm{t}}} - {\rm{j}}{k_0}z\cos {\theta _{\rm{t}}}} \right)\\ {\mathit{\boldsymbol{E}}_{\rm{r}}} = \mathit{\boldsymbol{R}}\left( {\cos {\theta _{\rm{r}}}\mathit{\boldsymbol{\hat x}} - \sin {\theta _{\rm{r}}}\mathit{\boldsymbol{\hat z}}} \right)\exp \left( {{\rm{j}}{k_0}x\sin {\theta _{\rm{r}}} + {\rm{j}}{k_0}z\cos {\theta _{\rm{r}}}} \right)\\ {\mathit{\boldsymbol{H}}_{\rm{i}}} = \mathit{\boldsymbol{\hat y}}\exp \left( {{\rm{j}}{k_0}x\sin {\theta _{\rm{i}}} - {\rm{j}}{k_0}z\cos {\theta _{\rm{i}}}} \right)/{\eta _0}\\ {\mathit{\boldsymbol{H}}_{\rm{t}}} = \mathit{\boldsymbol{T\hat y}}\exp \left( {{\rm{j}}{k_0}x\sin {\theta _{\rm{t}}} - {\rm{j}}{k_0}z\cos {\theta _{\rm{t}}}} \right)/{\eta _0}\\ {\mathit{\boldsymbol{H}}_{\rm{r}}} = - \mathit{\boldsymbol{R\hat y}}\exp \left( {{\rm{j}}{k_0}x\sin {\theta _{\rm{r}}} + {\rm{j}}{k_0}z\cos {\theta _{\rm{r}}}} \right)/{\eta _0} \end{array} \right. $

      (2)

      s偏振的电场和磁场分别描述为:

      $ \left\{ \begin{array}{l} {\mathit{\boldsymbol{E}}_{\rm{i}}} = \mathit{\boldsymbol{\hat y}}\exp \left( {{\rm{j}}{k_0}x\sin {\theta _{\rm{i}}} - {\rm{j}}{k_0}z\cos {\theta _{\rm{i}}}} \right)\\ {\mathit{\boldsymbol{E}}_{\rm{t}}} = \mathit{\boldsymbol{T\hat y}}\exp \left( {{\rm{j}}{k_0}x\sin {\theta _{\rm{t}}} - {\rm{j}}{k_0}z\cos {\theta _{\rm{t}}}} \right)\\ {\mathit{\boldsymbol{E}}_{\rm{r}}} = \mathit{\boldsymbol{R\hat y}}\exp \left( {{\rm{j}}{k_0}x\sin {\theta _{\rm{r}}} + {\rm{j}}{k_0}z\cos {\theta _{\rm{r}}}} \right)\\ {\mathit{\boldsymbol{H}}_{\rm{i}}} = \left( { - \cos {\theta _{\rm{i}}}\mathit{\boldsymbol{\hat x}} - \sin {\theta _{\rm{i}}}\mathit{\boldsymbol{\hat z}}} \right) \times \\ \;\;\;\;\;\;\;\frac{{\exp \left( {{\rm{j}}{k_0}x\sin {\theta _{\rm{i}}} - {\rm{j}}{k_0}z\cos {\theta _{\rm{i}}}} \right)}}{{{\eta _0}}}\\ {\mathit{\boldsymbol{H}}_{\rm{t}}} = \mathit{\boldsymbol{T}}\left( { - \cos {\theta _{\rm{t}}}\mathit{\boldsymbol{\hat x}} - \sin {\theta _{\rm{t}}}\mathit{\boldsymbol{\hat z}}} \right) \times \\ \;\;\;\;\;\;\;\frac{{\exp \left( {{\rm{j}}{k_0}x\sin {\theta _{\rm{t}}} - {\rm{j}}{k_0}z\cos {\theta _{\rm{t}}}} \right)}}{{{\eta _0}}}\\ {\mathit{\boldsymbol{H}}_{\rm{r}}} = \mathit{\boldsymbol{R}}\left( {\cos {\theta _{\rm{r}}}\mathit{\boldsymbol{\hat x}} - \sin {\theta _{\rm{r}}}\mathit{\boldsymbol{\hat z}}} \right) \times \\ \;\;\;\;\;\;\;\frac{{\exp \left( {{\rm{j}}{k_0}x\sin {\theta _{\rm{r}}} - {\rm{j}}{k_0}z\cos {\theta _{\rm{r}}}} \right)}}{{{\eta _0}}} \end{array} \right. $

      (3)

      式中,k0表示电磁波自由空间波数,η0表示真空中的波阻抗,R, T分别表示反射系数和透射系数。

      根据x偏振和y偏振的正交性,可将(1)式简化为:

      $ \left\{ {\begin{array}{*{20}{l}} {X_{{\rm{e}},{\rm{e}}}^{xx} = - \Delta {H_y}/\left( {{\rm{j}}\omega {\varepsilon _0}{E_{x,{\rm{av}}}}} \right)}\\ {\chi _{{\rm{e}},{\rm{e}}}^{yy} = \Delta {H_x}/\left( {{\rm{j}}\omega {\varepsilon _0}{E_{y,{\rm{av}}}}} \right)}\\ {\chi _{{\rm{m}},{\rm{m}}}^{yy} = - \Delta {E_x}/\left( {{\rm{j}}\omega {\mu _0}{H_{y,{\rm{av}}}}} \right)}\\ {\chi _{{\rm{m}},{\rm{m}}}^{xx} = \Delta {E_y}/\left( {{\rm{j}}\omega {\mu _0}{H_{x,{\rm{av}}}}} \right)} \end{array}} \right. $

      (4)

      式中,χe, exx的第1个上标表示场响应的方向,第2个上标表示场激励的方向,第1个下标表示电场响应,第2个下标表示电场激励,χm, myy, χe, eyy, χm, mxx同理。x极化波入射时,χe, exx表示x方向的电场激励引起x方向的电场变化;χm, myy表示y方向的磁场激励引起y方向的磁场变化;ΔHy表示y方向超表面两侧磁场的差;ΔEx表示x方向超表面两侧电场的差;Ex, av表示x方向超表面两侧的平均电场; Hy, av表示y方向超表面两侧的平均磁场。y极化波入射时,ΔHx表示x方向超表面两侧磁场的差;ΔEy表示y方向超表面两侧电场的差; Ey, av表示y方向超表面两侧的平均电场; Hx, av表示x方向超表面两侧的平均磁场。

      当平面电磁波垂直照射到超表面上(θi=0),且透射波的折射角为α时,结合(2)式、(3)式、(4)式得到的反射和透射系数分别为:

      $ \left\{ \begin{array}{l} \mathit{\boldsymbol{R}} = \frac{{2{\rm{j}}{k_0}\left( {{\mathit{\boldsymbol{\chi }}_{{\rm{m}},{\rm{m}}}} - {\mathit{\boldsymbol{\chi }}_{{\rm{e}},{\rm{e}}}}} \right)}}{{\left( {2 + {\rm{j}}{k_0}{\mathit{\boldsymbol{\chi }}_{{\rm{e}},{\rm{e}}}}} \right)\left( {2 + {\rm{j}}{k_0}{\mathit{\boldsymbol{\chi }}_{{\rm{m}},{\rm{m}}}}} \right)}}\\ \mathit{\boldsymbol{T}} = \frac{{4 + {\mathit{\boldsymbol{\chi }}_{{\rm{e}},{\rm{e}}}}{\mathit{\boldsymbol{\chi }}_{{\rm{m}},{\rm{m}}}}k_0^2}}{{\left( {2 + {\rm{j}}{k_0}{\mathit{\boldsymbol{\chi }}_{{\rm{e}},{\rm{e}}}}} \right)\left( {2 + {\rm{j}}{k_0}{\mathit{\boldsymbol{\chi }}_{{\rm{m}},{\rm{m}}}}} \right)}} \end{array} \right. $

      (5)

      即通过介质空间域的表面极化率、磁化率描述了相应超表面的反射和透射系数。

    • 基于全透射的超表面模型,设计将垂直入射电磁波转化为沿相对于z轴有45°折射的超表面,按照以下步骤实现。首先,利用(4)式计算超表面的极化、磁化系数,并验证设计方法的可行性。将超表面放置在x-y平面,取z=0。此时,电磁波的入射场(z < 0)和透射场(z>0)对应的x偏振的电场和磁场分别描述为:

      $ \left\{ \begin{array}{l} {\mathit{\boldsymbol{E}}_{\rm{i}}} = \mathit{\boldsymbol{\hat x}}\exp \left( { - {\rm{j}}{k_0}z} \right)\\ {\mathit{\boldsymbol{E}}_{\rm{t}}} = \mathit{\boldsymbol{T}}\left[ {(\sqrt 2 /2)\mathit{\boldsymbol{\hat x}} + (\sqrt 2 /2)\mathit{\boldsymbol{\hat z}}} \right] \times \\ \;\;\;\;\;\;\;\exp \left[ {{\rm{j}}{k_0}(\sqrt 2 /2)x - {\rm{j}}{k_0}(\sqrt 2 /2)z} \right]\\ {\mathit{\boldsymbol{H}}_{\rm{i}}} = \mathit{\boldsymbol{\hat y}}\exp \left( { - {\rm{j}}{k_0}z} \right)/{\eta _0}\\ {\mathit{\boldsymbol{H}}_{\rm{t}}} = \mathit{\boldsymbol{T\hat y}}\exp \left[ {{\rm{j}}{k_0}(\sqrt 2 /2)x - {\rm{j}}{k_0}(\sqrt 2 /2)z} \right]/{\eta _0} \end{array} \right. $

      (6)

      y偏振的电场和磁场分别描述为:

      $ \left\{ \begin{array}{l} {\mathit{\boldsymbol{E}}_{\rm{i}}} = \mathit{\boldsymbol{\hat y}}\exp \left( { - {\rm{j}}{k_0}z} \right)\\ {\mathit{\boldsymbol{E}}_{\rm{t}}} = \mathit{\boldsymbol{T\hat y}}\exp \left[ {{\rm{j}}{k_0}(\sqrt 2 /2)x - {\rm{j}}{k_0}(\sqrt 2 /2)z} \right]\\ {\mathit{\boldsymbol{H}}_{\rm{i}}} = \mathit{\boldsymbol{\hat x}}\exp \left( { - {\rm{j}}{k_0}z} \right)/{\eta _0}\\ {\mathit{\boldsymbol{H}}_{\rm{t}}} = - \mathit{\boldsymbol{T}}\left[ {(\sqrt 2 /2)\mathit{\boldsymbol{\hat x}} + (\sqrt 2 /2)\mathit{\boldsymbol{\hat z}}} \right] \times \\ \;\;\;\;\;\;\;\exp \left[ {{\rm{j}}{k_0}(\sqrt 2 /2)x - {\rm{j}}{k_0}(\sqrt 2 /2)z} \right]/{\eta _0} \end{array} \right. $

      (7)

      z=0处的x偏振和y偏振平面波的电磁场参量分别代入(4)式,得到的超表面等效极化、磁化系数可表示为:

      $ \left\{ {\begin{array}{*{20}{l}} {\chi _{{\rm{e}},{\rm{e}}}^{xx} = \frac{{4{\rm{j}}}}{{{k_0}}}\left\{ {\frac{{\exp \left[ {{\rm{j}}{k_0}(\sqrt 2 /2)x} \right] - 1}}{{2 + \sqrt 2 \exp \left[ {{\rm{j}}{k_0}(\sqrt 2 /2)x} \right]}}} \right\}}\\ {\chi _{{\rm{m}},{\rm{m}}}^{yy} = \frac{{\rm{j}}}{{{k_0}}}\left\{ {\frac{{\sqrt 2 \exp \left[ {{\rm{j}}{k_0}(\sqrt 2 /2)x} \right] - 2}}{{1 + \exp \left[ {{\rm{j}}{k_0}(\sqrt 2 /2)x} \right]}}} \right\}} \end{array}} \right. $

      (8)

      $ \left\{ {\begin{array}{*{20}{l}} {\chi _{{\rm{e}},{\rm{e}}}^{yy} = \frac{{\rm{j}}}{{{k_0}}}\left\{ {\frac{{\sqrt 2 \exp \left[ {{\rm{j}}{k_0}(\sqrt 2 /2)x} \right] - 2}}{{1 + \exp \left[ {{\rm{j}}{k_0}(\sqrt 2 /2)x} \right]}}} \right\}}\\ {\chi _{{\rm{m}},{\rm{m}}}^{xx} = \frac{{4{\rm{j}}}}{{{k_0}}}\left\{ {\frac{{\exp \left[ {{\rm{j}}{k_0}(\sqrt 2 /2)x} \right] - 1}}{{2 + \sqrt 2 \exp \left[ {{\rm{j}}{k_0}(\sqrt 2 /2)x} \right]}}} \right\}} \end{array}} \right. $

      (9)

      当全透射超表面被垂直入射时,由(5)式中的R=0,可得:χe, e=χm, m, 令χe, e=χm, m=χ,由(5)式可推出:

      $ \mathit{\boldsymbol{\chi }} = \frac{{2{\rm{j}}\left( {\mathit{\boldsymbol{T}} - 1} \right)}}{{{k_0}\left( {\mathit{\boldsymbol{T}} + 1} \right)}} $

      (10)

      基于全透射超表面,可将透射系数T表示为:T=Ψt/Ψi,其中,ΨiΨt分别表示入射波和透射波投射到超表面上的相位信息。由于入射波为垂直入射,入射波相位信息为Ψi= 1,则复振幅透射率T=Ψt=$\left[ { - {\rm{j}}{k_0}\left( {\sqrt 2 /2} \right)x} \right]$,结合(10)式可推出:

      $ \mathit{\boldsymbol{\chi }} = \frac{2}{{{k_0}}}\tan \left( {\frac{{{k_0}x}}{{2\sqrt 2 }}} \right) $

      (11)

      根据(8)式和(11)式,可得到将垂直入射电磁波转化为相对于z轴有45°折射的超表面的极化、磁化系数分布曲线, 如图 2所示。

      Figure 2.  Polarization and magnetization coefficients of metasurface for normal incident plane wave and refraction angle of 45°

      图 2中,λ0表示真空中的波长,实线和虚线分别表示极化、磁化系数的实部和虚部,图 2a图 2c中的极化、磁化系数由(8)式得到,图 2b图 2d中的极化、磁化系数由(11)式得到。(8)式中的极化、磁化系数为复数,虚部为负,其相应的超表面具有损耗,产生损耗的原因是入射波与透射波之间的法向功率流不相等[19-20];而(11)式中极化、磁化系数为纯实数,相应的超表面为无损结构。

      图 2还可看出,在自由空间中,极化、磁化系数的空间周期大于一个波长,此超表面容易由散射粒子实现。本文中将超表面的一个周期长度离散为7个具有不同透射相位的单元,单元结构采用以两层介质板为连接的3层金属结构,其中外层为具有相同结构的金属板,结构如图 3a所示。含两层介质的3层金属结构是实现全透射功能且具有2π相位覆盖所需的最少层数超表面结构单元[21]。图中,E1, H1, η1分别表示区域1处的电场、磁场、波阻抗;E2, H2, η2分别表示区域2处的电场、磁场、波阻抗。通常利用传输矩阵ABCD表示单元结构两侧的入射场和透射场的关系,图 3a中的3层级联导纳金属板用传输矩阵可描述为:

      Figure 3.  Scattering particles of metasurface

      $ \begin{array}{*{20}{c}} {\left[ {\begin{array}{*{20}{c}} \mathit{\boldsymbol{A}}&\mathit{\boldsymbol{B}}\\ \mathit{\boldsymbol{C}}&\mathit{\boldsymbol{D}} \end{array}} \right] = \left[ {\begin{array}{*{20}{c}} 1&0\\ {{Y_{{\rm{s}},1}}}&1 \end{array}} \right]\left[ {\begin{array}{*{20}{c}} {\cos (\beta d)}&{{\rm{j}}{\eta _{\rm{d}}}\sin (\beta d)}\\ {\frac{{{\rm{j}}\sin (\beta d)}}{{{\eta _{\rm{d}}}}}}&{\cos (\beta d)} \end{array}} \right] \cdot }\\ {\left. {\left[ {\begin{array}{*{20}{c}} 1&0\\ {{Y_{{\rm{s}},2}}}&1 \end{array}} \right]\left[ {\begin{array}{*{20}{c}} {\cos (\beta d)}&{{\rm{j}}{\eta _{\rm{d}}}\sin (\beta d)}\\ {\sin (\beta d)}&{\cos (\beta d)} \end{array}} \right]\left[ {\begin{array}{*{20}{c}} 1&0\\ {{Y_{{\rm{s}},3}}}&1 \end{array}} \right]} \right]} \end{array} $

      (12)

      式中,β是相移常数,d是介质厚度,ηd是介质的波阻抗,Ys, 1Ys, 2Ys, 3分别对应于图 3a单元中的3层金属板导纳。该单元结构的S参量(S11为反射系数,S21为透射系数)与用传输矩阵ABCD表示的级联导纳金属板的关系如下:

      $ \begin{array}{*{20}{c}} {\left[ {\begin{array}{*{20}{c}} {{\mathit{\boldsymbol{S}}_{11}}}&{{\mathit{\boldsymbol{S}}_{12}}}\\ {{\mathit{\boldsymbol{S}}_{21}}}&{{\mathit{\boldsymbol{S}}_{22}}} \end{array}} \right] = }\\ {\left( {\begin{array}{*{20}{c}} {\frac{{\mathit{\boldsymbol{A}} + \frac{\mathit{\boldsymbol{B}}}{{{\eta _0}}} - \mathit{\boldsymbol{C}}{\eta _0} - \mathit{\boldsymbol{D}}}}{{\mathit{\boldsymbol{A}} + \frac{\mathit{\boldsymbol{B}}}{{{\eta _0}}} + \mathit{\boldsymbol{C}}{\eta _0} + \mathit{\boldsymbol{D}}}}}&{\frac{{2\left( {\mathit{\boldsymbol{AD}} - \mathit{\boldsymbol{BC}}} \right)}}{{\mathit{\boldsymbol{A}} + \frac{\mathit{\boldsymbol{B}}}{{{\eta _0}}} + \mathit{\boldsymbol{C}}{\eta _0} + \mathit{\boldsymbol{D}}}}}\\ {\frac{2}{{\mathit{\boldsymbol{A}} + \frac{\mathit{\boldsymbol{B}}}{{{\eta _0}}} + \mathit{\boldsymbol{C}}{\eta _0} + \mathit{\boldsymbol{D}}}}}&{\frac{{ - \mathit{\boldsymbol{A}} + \frac{\mathit{\boldsymbol{B}}}{{{\eta _0}}} - \mathit{\boldsymbol{C}}{\eta _0} + \mathit{\boldsymbol{D}}}}{{\mathit{\boldsymbol{A}} + \frac{\mathit{\boldsymbol{B}}}{{{\eta _0}}} + \mathit{\boldsymbol{C}}{\eta _0} + \mathit{\boldsymbol{D}}}}} \end{array}} \right)} \end{array} $

      (13)

      本文中利用计算机仿真技术(computer simulation technology,CST)微波仿真软件模拟计算3层金属结构单元的传输系数S21,所有的仿真均采用周期性边界条件。

      通常用“工”字型单元结构实现超表面,如图 3b所示,但其调节对象比较单一、应用有限,研究人员对“工”字型单元结构进行了改进,提出的正交“工”字型结构单元,在设计超表面时涉及至少8个参量(AxAyBxByLxLyWxWy),如图 3c所示。本文中提出具有π型图案的金属单元仅包含5个参量(AyBxLxSxWy)。散射粒子参量的缩减可以有效地降低设计过程的复杂程度并缩短设计周期。

      本文中提出的π型结构单元的整体厚度t=3.04mm(约为λ0/10),单元尺寸为6mm×6mm(等于λ0/5),π型金属片结构如图 3d所示,金属层厚度为0.018mm,连接3层金属的两层介质采用介电常数εr=3的Rogers RO3003基底材料,介质层厚度为1.52mm,设计的仿真单元如图 3e所示。

      将7个具有梯度透射相位的单元按照顺序排列,这7个单元分别对x偏振和y偏振具有梯度透射相位,形成一个超晶格。在优化设计每一单元时,通过不断地进行迭代计算,优化各个单元的几何参量,使超晶格的透射能量达到最大值。优化过程如下:(1)通过CST微波仿真软件的迭代计算,优化超晶格的第1个单元;(2)达到预期性能后,更新超晶格;(3)移动到下一个单元,重复步骤(2),直到7个单元全部优化。

      通过上述步骤,本文中实现了全透射超表面,并将垂直入射的线偏振平面波转化为与z轴呈折射角为45°的线偏振平面波。通过不断优化,使超晶格的透射能量达到最大值,整个设计过程可调参量少,实现过程简单。

    • 为了验证π型金属结构高性能超表面的优势,分别设计了λ/4波片、波束折射超表面以及偏振分束器(polarization beam splitter, PBS)超表面。

    • 电磁波片是使相互垂直的两振动电场间产生特定相位差(Δφ=|φx-φy|)的一种光学器件,其中φxφy分别为x偏振和y偏振的透射相位。λ/4波片可以使沿与x轴,y轴分别成45°角的线偏振光波在通过λ/4波片后,使相互垂直的xy极化产生Δφ=90°的相位差,从而将线偏振光转换为圆偏振光或反向转换,如图 4所示。

      Figure 4.  Quarter wave-plate

      假设入射电场的振幅为1,入射场ϕi=0°、透射场ϕt=0°,x, y偏振分量分别为${\sqrt 2 /2}$, 出射后,xy偏振产生了相位差Δφ=|φx-φy|=90°,入射波、透射波的电场表达式分别为:

      $ \left\{ \begin{array}{l} {\mathit{\boldsymbol{E}}_{\rm{i}}} = [(\sqrt 2 /2)\mathit{\boldsymbol{\hat x}} + (\sqrt 2 /2)\mathit{\boldsymbol{\hat y}}]\exp \left( { - {\rm{j}}{k_0}z} \right)\\ \mathit{\boldsymbol{E}}_{\rm{t}}^x = (\sqrt 2 /2)\mathit{\boldsymbol{T}}\exp \left( { - {\rm{j}}{k_0}z + {\varphi _x}} \right)\\ \mathit{\boldsymbol{E}}_{\rm{t}}^y = (\sqrt 2 /2)\mathit{\boldsymbol{T}}\exp \left( { - {\rm{j}}{k_0}z + {\varphi _y}} \right) \end{array} \right. $

      (14)

      式中,Etx表示x偏振透射波的电场;Ety表示y偏振透射波的电场。

      由于本文中设计的λ/4波片是将垂直入射的线偏振波转化为垂直出射的圆偏振波,因此设计的超表面是均一的,并不需要梯度相位排列,从而只需将满足条件的单元结构周期性排列即可。具体的仿真单元如图 3e所示、单元设计尺寸如表 1所示。

      Table 1.  Dimensions of quarter wave-plate

      parameters Ay/mm Bx/mm Lx/mm Sx/mm Wy/mm
      outer layers 0.75 5.5 0.6 1.55 4.25
      middle layer 0.75 5 1 1.25 4.25

      图 5a图 5b中分别为电磁波垂直入射到λ/4波片的透射系数曲线(S21)的幅度和相位。从图中可看到, 在10GHz时,S21幅度基本达到全透射,且φx=32.9°,φy=128.2°,Δφ=|φx-φy|接近90°,实现了将线偏振转换为圆偏振的功能。

      Figure 5.  Transmission coefficient S21 of quarter wave-plate

    • 在10GHz频率处设计一个将垂直入射的线偏振平面波转换为折射ϕt=45°的全透射超表面。设计的超表面的性能参量如下:折射角ϕt=45°,入射角ϕi=0°,波长λ0=30mm,单元尺寸d=6mm,厚度t=3.04mm,单元个数N=7。

      假设电磁波的入射场ϕi=0°、透射场ϕt=45°的x偏振入射波、透射波的电场、磁场表达式如(6)式所示。

      7个单元都被分别设计为对x偏振和y偏振具有特定的梯度透射相位φxφy,各个单元的传输相移如表 2所示。

      Table 2.  Transmission phase shifts for x and y polarization of 7 unit cells forming the metacell

      No. 1 2 3 4 5 6 7
      φx/(°) 0 51.5 103 154.5 206 257.5 309
      φy/(°) 309 257.5 206 154.5 103 51.5 0

      一个构成全透射超表面的超晶格结构优化后的内外层π型金属层设计如图 6a图 6b所示,各个单元设计尺寸如表 3表 4所示。

      Figure 6.  Design shapes of outer and inner layers of beam refractive metasurface metalattices

      Table 3.  Design dimensions of beam refraction metasurface inner layer of π-shaped metal

      No. Ay/mm Bx/mm Lx/mm Sx/mm Wy/mm
      1 0.75 5 1 1.25 4.25
      2 0.75 5 1 1.25 4.25
      3 0.75 4.5 1 1.25 4.25
      4 0.5 4 1 0.5 4.5
      5 0.25 3.75 1.5 1.25 4.5
      6 0.25 3.75 1.5 1.5 3.5
      7 0.75 4.25 0.75 1 4.5

      Table 4.  Design dimensions of beam refraction metasurface outer layer of π-shaped metal

      No. Ay/mm Bx/mm Lx/mm Sx/mm Wy/mm
      1 0.75 5.5 1.025 1.175 4.25
      2 0.75 5.5 0.6 1.55 4.25
      3 0.75 5.5 0.45 1.45 4
      4 0.5 5 1 0.5 4
      5 0.75 4.25 1.25 0.75 3.5
      6 0.5 4.25 1 0.25 3.25
      7 0.25 4.25 2 1 3.75

      电磁波垂直入射到波束折射超表面的透射系数曲线(S21)和平面波折射的电场等相位面情形图分别如图 7a图 7b所示。由于7个相邻单元之间的相互耦合,全透射频率向上漂移到10.26GHz,此时透射系数的幅值为入射电场的1/3。

      Figure 7.  Simulation results of beam refraction metasurface

      图 7b中的仿真结果可见,沿z轴正方向入射的x偏振电磁波在照射到超表面后,沿与z轴成45°从超表面的右上侧出射,实现了垂直入射电磁波被超表面折射45°的结果。

    • 广义的双折射是指独立控制正交极化波的反射角、折射角及其幅度。全透射的偏振分束(PBS)超表面通常使垂直入射的x偏振和y偏振平面波的透射波产生不同方向的偏折,从而实现不同极化的电磁波的波束分离,如图 8所示。

      Figure 8.  Metasurface of the polarized beam splitting

      假设电磁波的入射场ϕi=0°、透射场φx=45°,φy=135°,出射后,xy偏振分离,其x偏振和y偏振的入射波电场(EixEiy)、透射波电场(EtxEty)表达式分别为:

      $ \left\{ \begin{array}{l} \mathit{\boldsymbol{E}}_{\rm{i}}^x = \mathit{\boldsymbol{\hat x}}\exp \left( { - {\rm{j}}{k_0}z} \right)\\ \mathit{\boldsymbol{E}}_{\rm{i}}^y = \mathit{\boldsymbol{\hat y}}\exp \left( { - {\rm{j}}{k_0}z} \right)\\ \mathit{\boldsymbol{E}}_{\rm{t}}^x = T[(\sqrt 2 /2)\mathit{\boldsymbol{\hat x}} + (\sqrt 2 /2)]\mathit{\boldsymbol{\hat z}} \times \\ \;\;\;\;\;\;\;\exp \left[ {{\rm{j}}{k_0}(\sqrt 2 /2)x - {\rm{j}}{k_0}(\sqrt 2 /2)z} \right]\\ \mathit{\boldsymbol{E}}_{\rm{t}}^y = \mathit{\boldsymbol{T}}[(\sqrt 2 /2)\mathit{\boldsymbol{\hat y}} + (\sqrt 2 /2)\mathit{\boldsymbol{\hat z}}] \times \\ \;\;\;\;\;\;\;\exp \left[ {{\rm{j}}{k_0}(\sqrt 2 /2)y - {\rm{j}}{k_0}(\sqrt 2 /2)z} \right] \end{array} \right. $

      (15)

      本文中设计的超表面在x方向和y方向透射系数梯度相位分布具有非对称性,如图 9a所示,在x方向实现对x极化电磁波透射系数相位梯度减小,在y方向实现对y极化电磁波透射系数相位梯度增加,每个单元的几何尺寸如表 3表 4所示。

      Figure 9.  Supercell with 7×7 unit cells

      全透射偏振分束超表面的透射曲线(S21)如图 9b所示,由于相邻单元之间的相互耦合,全透射频率略有偏移,此时x偏振和y偏振的透射系数S21幅值分别为-8.1dB和-5.6dB。考虑到7×7的超晶格中涉及了大量单元,导致全波仿真由于对计算机的性能要求太高而无法模拟,所以仍采用7×1的超晶格代替了7×7的超晶格,选择7×7超晶格的第1行和第1列分别进行仿真。仿真结果分别如图 9c图 9d所示。从图中可发现,x偏振电磁波和y偏振电磁波被分别折射为φx=45°和φy=135°。在仿真过程中,x偏振波和y偏振波均为沿z轴正方向垂直入射到超表面,被超表面折射后沿两个不同的方向出射,仿真结果达到了预期效果。

    • 本文中根据极化、磁化系数张量,阐述了超表面的介质等效性质以及物理实现过程。采用了π型图案的级联金属单元,设计了具备高透射率的超表面模型,并将其应用到λ/4波片、波束折射超表面和偏振分束超表面的设计中。结果表明,λ/4波片的透射系数S21幅度达到全透射,实现了线偏振转圆偏振;波束折射超表面的损失为-8.3dB,实现了将垂直入射电磁波被超表面折射45°;偏振分束超表面的x偏振和y偏振的透射系数损失分别为-8.1dB和-5.6dB,实现了xy偏振分离,仿真均达到了预期效果,充分验证了基于π型金属结构所设计的超表面模型具有设计参量少,易于实现高透射率等特点,对未来高性能超表面的设计实现具有积极意义。

参考文献 (21)

目录

    /

    返回文章
    返回