1、文章编号：1002-2082(2024)02-0282-10基于超宽带反射超表面产生太赫兹轨道角动量韩晓冰1，陈江侠1，王焱2,3,4，周远国1（1.西安科技大学通信与信息工程学院，陕西西安710054；2.中国航空研究院研究生院，江苏扬州225000；3.电磁环境效应航空科技重点实验室，辽宁沈阳110035；4.电磁信息控制与效应全国重点实验室，辽宁沈阳110035）摘摘 要：要：由于轨道角动量（orbitalangularmomentum，OAM）有望成为新的通信复用自由度，在拓展信道容量、提高频谱资源利用率方面有巨大潜力而受到越来越多研究人员的关注。目前的太赫兹涡旋波产生器件存在只能工作
2、在单频点、带宽窄、转换效率低等问题，如何在太赫兹频段高效产生 OAM 成为关键问题之一。该文设计了一种超宽带反射超表面单元，结合几何相位原理和相位叠加原理设计了单层反射超表面。仿真结果表明：设计的超表面实现了在 0.82THz2.09THz（相对带宽 87.3%）的宽频范围内，将圆极化太赫兹波转换为携带轨道角动量的太赫兹涡旋波，同极化反射谱的幅度大于 0.97，转换效率高于 94.7%，反射相位覆盖 0360。利用傅里叶变换分解反射场中各个 OAM 模式，定量分析了 OAM 模式纯度，在不同频率下涡旋波中均为主模式l=2 能量占比最高，进一步对设计的超表面进行优化，主模式能量占比明显提升。该文
3、设计的超表面具有转换效率高、工作带宽大、主模式强度高等优势，为宽带太赫兹涡旋波的高效产生提供了一种参考。关键词：关键词：轨道角动量；太赫兹；超宽带；超表面中图分类号：TN82;O43文献标志码：ADOI：10.5768/JAO202445.0201002Terahertz orbital angular momentum generation based onultra-wideband reflective metasurfaceHANXiaobing1，CHENJiangxia1，WANGYan2,3,4，ZHOUYuanguo1（1.CollegeofCommunicationandIn
4、formationEngineering,XianUniversityofScienceandTechnology,Xian710054,China；2.GraduateSchoolofChineseAeronauticalEstablishment,Yangzhou225000,China；3.AviationKeyLaboratoryofScienceandTechnologyonElectromagneticEnvironmentEffects,Shenyang110035,China；4.NationalKeyLaboratoryofElectromagneticInformation
5、ControlandEffects,Shenyang110035,China）Abstract：Orbitalangularmomentum(OAM)hasbeenincreasinglyinvestigatedbyresearchersbecauseitisexpectedtobeanewphysicalquantityforcommunicationmultiplexingandhasgreatpotentialforexpandingchannel capacity and improving spectrum resource utilization.Current terahertz
6、 vortex wave generationdevices are limited by operating at only a single frequency,having a narrow bandwidth and having lowconversionefficiency,sohowtoefficientlygenerateOAMintheterahertzbandhasbecomeoneofthekeyissues.Anultra-widebandreflectivemeta-atomwasproposed,andasingle-layerreflectivemetasurfa
7、cewasdesignedbycombiningthePancharatnam-Berryphaseconceptandthephasesuperpositionprinciple.Thesimulationresultsshowthatitachievestheconversionofcircularlypolarizedterahertzwavesintoterahertzvortexwavescarryingorbitalangularmomentuminawidefrequencyrangefrom0.82THzto2.09THz(withtherelativebandwidthof8
8、7.3%).Theamplitudeoftheco-polarizedreflectionspectrumishigherthan0.97,theconversionefficiencyismorethan94.7%,andthereflectionphasecovers0360.TheFouriertransform收稿日期：2023-06-21；修回日期：2023-08-26基金项目：国家自然科学基金（61801009）；陕西省自然科学基金面上项目（2024JC-YBMS-556）作者简介：韩晓冰（1965），男，教授，主要从事无线通信和超表面阵列天线等研究。E-mail：通信作者：周远国
9、（1981），男，博士，副教授，主要从事计算电磁学和计算光学等研究。E-mail：第45卷第2期应用光学Vol.45No.22024年3月JournalofAppliedOpticsMar.2024wasusedtodecomposeeachOAMmodeinthereflectedfield,andtheOAMmodepuritywasquantitativelyanalyzed,withthehighestenergyweightshareofthedominantOAMmodel=2inallvortexwavesatdifferentfrequencies,andthedesigned
10、metasurfacewasfurtheroptimizedtoincreasetheenergyshareofthedominantmode.Thedesignedmetasurfacehastheadvantagesofhighconversionefficiency,largeoperatingbandwidth,andhighmainmodeenergy,whichprovidesareferencefortheefficientgenerationofwidebandterahertzvortexwaves.Key words：orbitalangularmomentum；terah
11、ertz；ultra-wideband；metasurface引言太赫兹（Terahertz，THz）波是一种波长介于微波与红外之间的电磁波，波长范围在 0.03mm3mm，频率范围一般为 0.3THz10THz。与微波相比，太赫兹波波束较窄，具有通信带宽大、抗干扰能力强、分辨率高、光子能量低等优势，并且兼具毫米波的穿透性与光波的操控性，太赫兹在光谱检测1、生物成像2-3、全息投影4、高速通信5-6等领域已经展现出广阔的应用前景。电磁波携带线动量和角动量，其中角动量包括自旋角动量（spinangularmomentum，SAM）和轨道角动量（orb-italangularmomentum，OA
12、M），OAM 作为电磁波的基本物理量，在 1992 年由 ALLENL 等7人在拉盖尔-高斯（Laguerre-Gaussian，LG）光束中发现并定义后受到了国内外研究人员的广泛关注。携带 OAM的电磁波束具有中心为 0 的环形强度分布和螺旋的空间相位分布 exp（il），所以又将其称为涡旋波束，其中 是空间方位角，l 为 OAM 拓扑电荷数，又被称为模式数，可以取小数或整数。由于具有不同模式的涡旋波相互正交互不干扰，使得 OAM有望成为全新的传输复用自由度，在光学、通信领域中有着广泛的应用。目前 OAM 已经在光学微操控8、超分辨率成像9和高速率信息传输10等领域成功应用。因此，如何在宽带
13、高效地产生太赫兹涡旋波成为相关领域研究的关键问题之一。传统的涡旋波产生方法有环形阵列天线11、螺旋相位板12和全息衍射光栅等。其中，环形阵列天线需要设计复杂的馈电网络，系统成本高；其他光学器件尺寸大、造价高，不利于集成化和小型化。超表面13-14作为一种新颖的由亚波长单元组成的人工结构，能够灵活地调控电磁波的振幅、相位和极化，并且避免了复杂的馈电设计，具有体积小、重量轻、制造成本低和易于加工等优势而受到研究人员的青睐，超表面成为产生太赫兹涡旋波束的一种理想方法。近年来，越来越多的研究人员尝试通过超表面来产生涡旋波，如 2011 年，YUN 等15提出的“V”字单元结构的超表面，实现了在红外波段
14、产生 l=1 的涡旋光；2016 年，CHANGZ等16利用在超表面加载可重配置的石墨烯，在 1.6THz处产生了涡旋波束；2017 年，SHIY 等17提出了一种基于多层石墨烯结构的反射式超表面，在1.8THz2.8THz 频率范围内产生了涡旋光束；2019年，李晓楠等18提出十字架结构单元超表面在 0.8THz1.4THz 产生了涡旋波束，转换效率为 86%；2022 年，仲敏等19提出了频率可切换超表面，分别在 0.7THz 和 1.23THz 产生了不同模式的涡旋波。上述超表面在产生太赫兹涡旋波时存在工作频率单一、带宽较小、转换效率低等问题。本文提出了一种利用单层几何相位单元组成的宽带
15、太赫兹涡旋波超表面，几何相位又名panch-aratnam-berry（P-B）相位，利用 P-B 相位与频率无关，仅与单元旋转角度相关的相位调控相关，实现了在太赫兹频段高效率产生涡旋波。通过旋转超表面单元获得 0360反射相位范围，根据补偿相位分布构建了 l=2 的超表面，并对其进行了数值模拟。在圆极化平面波激励下的仿真结果表明，本文提出的超宽带太赫兹涡旋波超表面具有以下优势：1）实现了在 0.82THz2.09THz（相对带宽87.3%）宽带内产生涡旋波；2）超表面的转换效率高于 94.7%；3）OAM 模式纯度高，为太赫兹涡旋波的高效产生提供了一种参考。1 理论分析与单元设计反射型超表面
16、产生涡旋波束的原理是，将具有 SAM 的平面波反射为携带 OAM 的结构波。为了简化设计过程，以及更加方便地设计出高性能单元，本文使用琼斯矩阵对超表面的入射场和散射场进行分析。如图 1 所示，在直角坐标系中，将超表面单元旋转角度 后，在圆极化波入射下SAM-OAM 的转换过程如下：应用光学2024，45（2）韩晓冰，等：基于超宽带反射超表面产生太赫兹轨道角动量283(a)直角坐标(b)旋转坐标yyvuxx图 1 超表面单元旋转示意图Fig.1 Schematic diagram of reflective meta-atom with ro-tation anglerll=12(rxxryy)
17、+j(rxy+ryx)exp(j2)（1）rlr=12(rxx+ryy)+j(rxyryx)（2）rrr=12(rxxryy)j(rxy+ryx)exp(j2)（3）rrl=12(rxx+ryy)j(rxyryx)（4）式中：rxx，ryy分别为 x-，y-极化波垂直入射下的同极化反射系数；rll和 rrr分别为左圆和右圆极化波垂直入射下的同极化反射系数；rxy，ryx，rlr和 rrl为对应的交叉极化反射系数；为超表面单元相对 x 轴的旋转角度。从式（1）式（4）中可以看出，当反射场的圆极化状态相对于入射场不变时，反射场将携带一个额外的、与单元旋转角 相关的相位项 exp（j2）和 exp（
18、j2），当反射场的圆极化状态被转换时，在入射场和散射场之间将不存在相位差。根据 P-B 相位理论，可以通过旋转超表面单元的方法获得 0360的任意相位差，将螺旋相位波前调制成反射波，这是超表面能够产生涡旋波束的必要条件。利用具有不同旋转角的单元构造出设计的超表面补偿相位面，可以实现对反射电磁波的任意预期波前调控。在圆极化波入射下，为了实现较高的 SAM-OAM 转换效率，超表面单元的反射系数 rll和 rrr的幅度应该尽可能地接近于1。因此，实现高效同极化转换只需要设计的超表面单元满足如下条件：|rxx|?ryy?1（5）?arg(rxx)arg(ryy)?=180（6）为了实现更大的系统工作
19、带宽，在宽带宽中高效地产生涡旋波束，要求超表面单元在不同旋转角度下同极化反射系数的相位在（fmin，fmax）带宽内具有相似的斜率，可用公式表示为(l,l)fmin(l,l)fi(l,l)fmax（7）式中：fmin和 fmax分别为较低和较高的截止频率；fi为（fmin，fmax）频带内的任意频率；（l，l）为左旋圆极化（left-handcircularpolarization，LCP）波入射下的反射相位。在右旋圆极化（right-handcircularpol-arization，RCP）波入射下时与之相同。基于上述讨论，为了实现高性能的 SAM-OAM转换，本文设计了一个各向异性的超表
20、面单元，如图 2(a)所示。该单元由 3 层结构组成，如图 2(b)所示。上层为中心挖口，上下对称的“凹”字型挖槽矩形金属；中间层为介质基板，其材料为聚四氟乙烯（F4B，r=2.65，tan=0.001），高度 h=30m；下层为金属地。该单元为镜像对称结构，周期为 p=100m。为了实现更大的带宽，结合遗传算法对单元结构参数进行优化，最终得出满足上述理论要求的超表面单元结构，如图 2(c)所示。其具体参数为 a=28m，b=80m，a1=5m，b1=20m，a2=10m，b2=20m。通过对上层金属结构进行旋转，使超表面单元具有连续的 0360反射相位范围，提出的超表面单元的电尺寸约为中心频
21、率 fc=1.5THz的 0.50（0为中心频率波长）。(a)三维侧视图(c)单元结构示意图hpxzyPECPECF4Babb1a1a2b2(b)3层结构示意图图 2 超表面单元示意图Fig.2 Schematic diagram of reflective meta-atom284应用光学第45卷第2期2 宽带涡旋波超表面设计考虑到本文设计的超表面中心频率为 1.5THz，由 MN 个单元组成，设计中选取 M=N=16。超表面需要具有涡旋相位分布才能使反射波携带涡旋相位，因此，为了在超表面反射场中产生 OAM 涡旋波束，需要在超表面上调制 exp(il)的相位分布。超表面产生涡旋波原理如图
22、3 所示。x(xi,yi)y反射波入射波图 3 超表面产生涡旋波原理示意图Fig.3 Schematic diagram of OAM vortex wave generationbased on metasurface将超表面设置在图 3 所示 x-o-y 平面上，则超表面中的任意一个单元的位置可以用坐标表示为（xi，yi）。涡旋相位分布可表示为1(xi,yi)=larctan(yi/xi)（8）式中 l 表示相位奇点的 OAM 模式数，也称拓扑电荷数，可以为任意整数或小数。本文将以 l=2 为例说明其原理，如图 4（a）所示。由于馈源天线辐射球面波不满足涡旋超表面的平面波入射条件，需要结合
23、聚焦相位将球面波转化为平面波。本文将馈源天线放置在焦点 F处，其聚焦相位分布函数表示为2(xi,yi)=k0(F2+x2i+y2iF)（9）式中：k0（k0=2/0）为自由空间波数；F 为焦距，大小可以在设计中自由选择，本文设为 F=1000m。聚焦相位分布如图 4（b）所示。根据相位叠加原理，每个位置的超表面单元所需要的补偿相位可以由下式计算：tot(xi,yi)=1(xi,yi)+2(xi,yi)（10）根据补偿相位可以得到任意位置处单元的旋转角度为 tot/2。图 4（c）所示的数值结果表明，反射波的相位分布覆盖 0360并呈现螺旋状分布。理论上，通过补偿相位计算并排列不同旋转角度的超表
24、面单元，可以得到产生具有任意拓扑电荷数的连续空间相变 OAM 涡旋波超表面。148121614812160180360 x/m(a)涡旋相位1148121614812160180360 x/m(b)聚焦相位2148121614812160180360 x/m(c)总补偿相位toty/my/my/m相位/()相位/()相位/()图 4 超表面补偿相位计算示意图Fig.4 Schematic diagram of calculating phase distributionof reflective metasurface3 仿真结果与分析3.1 单元仿真结果与分析前面我们给出了超表面单元的结构及
25、参数，超表面单元的仿真设置和仿真结果如图5 所示。如图5（b）所示，在 x-和 y-极化波入射下反射谱的幅度都接近于1，满足式（5）；在宽频带内具有 180的相位差，满足式（6），这两点是在宽带内产生涡旋波的必要条件。应用光学2024，45（2）韩晓冰，等：基于超宽带反射超表面产生太赫兹轨道角动量285圆极化波入射下单元的转换效率可通过下式计算：Eefficiency=r2ll/(r2ll+r2lr)（11）yxz0.61.01.41.82.2频率/THz(b)线极化波激励下的反射谱和相位差00.20.40.60.81.0|rxx|&|ryy|8006004002000200400|rxx|r
26、yy|arg(rxx)arg(ryy)arg(rxx)arg(ryy)0.61.01.41.82.2频率/THz(c)圆极化波激励下的反射谱及转换效率00.20.40.60.81.000.20.40.60.81.0|rll|&|rlr|效率|rll|rlr|效率0.81.21.62.0频率/THz(d)圆极化波激励下不同旋转角的反射谱00.20.40.60.81.0|rll|60040020002000306090120150(a)单元仿真设置示意图反射相位/()反射相位/()模式2模式1入射端口图 5 超表面单元仿真设置和仿真结果Fig.5 Simulation settings and r
27、esults of reflective meta-atom如图 5（c）所示，在圆极化波入射下，设计的宽带单元在 0.82THz2.09THz（相对带宽为 87.3%）频率范围内，其同极化反射系数幅度均大于 0.97，转换效率均高于 94.7%，交叉极化反射系数幅度均小于 0.23。更重要的是，从图 5（d）中可以观察到，对于不同旋转角度下的超表面单元，其反射谱的相位如预期的那样在带宽频率范围内保持平行，满足式（7），并且同极化反射谱的幅度均高于 0.95，这两点对于构建具有高性能、高纯度 OAM 特性的超表面十分重要。上面超表面单元的电磁响应仿真结果，表明了本文提出的超表面单元具有高性能，
28、能够用于构造高效的宽带涡旋波超表面。表 1 给出了我们设计的单元与其他在太赫兹频段产生 OAM 涡旋波束的单元性能比较。表 1 中采用的数据均为相应文章中公开的仿真数据，结果表明，我们设计的单元可以在更大的带宽内实现高效工作，同时在带宽和效率之间获得了良好的平衡。表 1 OAM 超表面单元比较Table 1 Comparison with meta-atoms for OAM vortexwave generation参考文献工作频率/THz 相对带宽/%效率/%调相方式数据来源171.82.843.5石墨烯化学势 仿真180.81.454.592几何相位仿真200.91.866.790几何相
29、位仿真210.30.4540几何相位仿真225.17.436.880传输相位仿真本文0.822.0987.394.7几何相位仿真3.2 超表面仿真结果与讨论根据上面分析的 SAM-OAM 转换概念，通过将上述超表面单元以一定的旋转角度在不同位置进行排列，就可以构造出具有期望涡旋相位分布exp(il)的超表面。本文设计的 l=2 的超表面如图 6 所示。图 6 中 l 的正负仅与入射波的极化方图 6 l=2 的超表面拓扑结构Fig.6 Topology of proposed metasurface with OAM modeof l=2286应用光学第45卷第2期向有关，具体拓扑结构由 161
30、6 个旋转的超表面单元阵列组成。为了说明提出的超表面能够在宽带上产生涡旋波，取 1THz2THz 频率范围对其进行数值模拟，并用近场和远场电磁性能验证它的涡旋特性。在 RCP 波激励下，本文提出的超表面产生了模式数为 l=2 的太赫兹涡旋波。对不同频率下反射场中的右旋分量进行采样和分解，采样平面尺寸为 1600m1600m，如图 7 所示。从采样场中右旋分量可以发现电场幅度表现为空心环形，相位分布呈现为围绕奇点的对称双螺旋状，相位变化范围在 0360，这与 OAM 特征一致，证明了本文设计的超表面在宽带上可使反射电磁波携带 OAM。图 8 给出了阿基米德螺旋天线激励的反射场中右旋分量的几个典型
31、频率下三维远场辐射方向图，所有频率下的辐射增益均大于 10dB。MaxMin36001.00 THz2.00 THz1.25 THz1.75 THz1.50 THz图 7 不同频率下 l=2 涡旋波近场仿真结果的幅度与相位Fig.7 Amplitudes and phases of simulated near-field vortexwaves for l=2 at different frequencies(a)1.00 THz(e)2.00 THz(b)1.25 THz(d)1.75 THz(c)1.50 THz10505101520 xyz1510505101520 xyz151050
32、5101520 xyz1510505101520 xyz1510505101520 xyz增益/dB增益/dB增益/dB增益/dB增益/dB(a)1.00 THz(e)2.00 THz(b)1.25 THz(d)1.75 THz(c)1.50 THz10505101520 xyz1510505101520 xyz1510505101520 xyz1510505101520 xyz1510505101520 xyz增益/dB增益/dB增益/dB增益/dB增益/dB图 8 不同频率下 l=2 涡旋波三维远场辐射方向图Fig.8 3D far-field radiation diagram of v
33、ortex waves forl=2 at different frequenciesOAM 模式纯度可以用来描述反射场中各种不同模式数的涡旋波的能量占比分布情况，能够很好地说明涡旋波质量。为了定量分析 OAM 模式纯度，本文利用傅里叶变换来分解反射场中各个OAM 模式。以涡旋波束的相位奇点为圆心，沿主波束选取一个环形电场数据，并对其进行傅里叶变换。相应的计算公式如下：Al=12w20E()exp(i l)d（12）Eenergy weight=Al7l=7Al（13）式中：E()为选取的环形电场；Al为各个模态的幅度。考虑 l=7l=7 的 OAM 模式，傅里叶谱分析结果如图（9）所示。从图
34、（9）可以看出，在不同频率下的 OAM 模式占比中，主模式 l=2 所占的能量都最高，表明超表面产生的涡旋波主模式与超表面设计的模式相同。3.3 超表面优化仿真与讨论上面我们对由 1616 个单元组成的超表面进行了仿真及模式强度占比分析，从图 9 中可以看出，在各个频率下主模式 l=2 的强度占比均在 0.6左右。为了提升主模式强度占比，我们首先对超表面尺寸进行优化，增加单元数量，设计了由 2020个超表面单元组成的矩形超表面，如图 10（a）所示。矩形超表面在不同频率下反射的涡旋波和对应谱如图 11 所示。应用光学2024，45（2）韩晓冰，等：基于超宽带反射超表面产生太赫兹轨道角动量287
35、8*p(a)l=2矩形超表面结构(b)l=2圆形超表面结构图 10 优化的超表面拓扑结构Fig.10 Topology of optimized metasurface with OAMmode of l=282从图 11 可以看出，各个频率下主模式 l=2 的能量占比在 0.76 左右，比 1616 个单元组成的矩形超表面明显增加。矩形是超表面设计中最常用的形状，我们将其优化为圆形研究其对主模式强度的影响，设计了半径为p 的圆形超表面，如图 10（b）所示。图 12 给出了圆形超表面的仿真结7654321 0 1 2 3 4 5 6 700.20.40.60.81.0(a)1.00 THzM
36、axMin3600场幅度|E|场相位E7654321 0 1 2 3 4 5 6 700.20.40.60.81.0(b)1.25 THzMaxMin3600场幅度|E|场相位E权重权重模式数 模式数 7654321 0 1 2 3 4 5 6 700.20.40.60.81.0(a)1.00 THz7654321 0 1 2 3 4 5 6 700.20.40.60.81.0(b)1.25 THz7654321 0 1 2 3 4 5 6 700.20.40.60.81.0(c)1.50 THz7654321 0 1 2 3 4 5 6 700.20.40.60.81.0(d)1.75 TH
37、z7654321 0 1 2 3 4 5 6 700.20.40.60.81.0(e)2.00 THz权重权重权重权重权重模式数 模式数 模式数 模式数 模式数 图 9 不同频率下 OAM 谱分析Fig.9 OAM spectrum weight for simulated results atdifferent frequencies288应用光学第45卷第2期7654321 0 1 2 3 4 5 6 700.20.40.60.81.0(c)1.50 THzMaxMin3600场幅度|E|场相位E7654321 0 1 2 3 4 5 6 700.20.40.60.81.0(d)1.75
38、THzMaxMin3600场幅度|E|场相位E7654321 0 1 2 3 4 5 6 700.20.40.60.81.0(e)2.00 THzMaxMin3600场幅度|E|场相位E权重权重权重模式数 模式数 模式数 图 11 矩形超表面在不同频率下反射的涡旋波和对应谱分析Fig.11 Near-field observation and corresponding spectralanalyses of vortex wave at different frequencies re-flected by rectangular metasurface7 6 5 4 3 2 1 0 1 2
39、 3 4 5 6 700.20.40.60.81.0(b)1.25 THzMaxMin3600场幅度|E|场相位E7 6 5 4 3 2 1 0 1 2 3 4 5 6 700.20.40.60.81.0(c)1.50 THzMaxMin3600场幅度|E|场相位E7 6 5 4 3 2 1 0 1 2 3 4 5 6 700.20.40.60.81.0(d)1.75 THzMaxMin3600场幅度|E|场相位E7 6 5 4 3 2 1 0 1 2 3 4 5 6 700.20.40.60.81.0(e)2.00 THzMaxMin3600场幅度|E|场相位E权重权重权重权重模式数 模式数
40、模式数模式数图 12 圆形超表面在不同频率下反射的涡旋波和对应谱分析Fig.12 Near-field observation and corresponding spectralanalyses of vortex wave at different frequencies re-flected by circular metasurface7 6 5 4 3 2 1 0 1 2 3 4 5 6 700.20.40.60.81.0(a)1.00 THzMaxMin3600场幅度|E|场相位E权重模式数 应用光学2024，45（2）韩晓冰，等：基于超宽带反射超表面产生太赫兹轨道角动量289果。
41、从图 12 可以看出，圆形超表面反射的涡旋波主模式 l=2 的强度在各频率下均提升到了 0.88 左右，验证了圆形超表面设计能够有效地提升主模式强度。仿真计算的意义是为实践提供参考，超表面在实际加工中不可避免地存在误差，因此本文将可能产生的误差考虑进仿真，以研究其对反射涡旋波模场分布的影响。随机选取尺寸误差为2m，旋转角度误差分布如图 13（a）所示。构造圆形超表面并以中心频率进行仿真，仿真结果如图 13（b）所示。从图 13 可以看出，在误差影响下，主模式强度明显降低，其他串扰模式强度有所增大，但就整体而言，在一定范围内误差对反射涡旋波的模场分布影响不大，主模式强度仍然占据主要部分。2101
42、2角度误差/()(a)旋转角度误差分布7654321 0 1 2 3 4 5 6 700.20.40.60.81.0(b)有误差超表面仿真结果(1.5 THz)MaxMin3600场幅度|E|场相位E无误差有误差权重模式数 图 13 旋转角度误差分布及有误差超表面在中心频率下反射的涡旋波和对应谱Fig.13 Rotation angle error distribution,near-field obser-vation and corresponding spectral analyses of vor-tex wave at center frequencies4 结论本文提出了一种高性能
43、的几何相位单元，能够对电磁波的相位进行自由调控。结合几何相位原理和相位叠加原理设计了 l=2 的宽带太赫兹涡旋波超表面。对超表面的近场和远场性能进行了仿真验证。仿真结果表明，设计的超表面可以在 0.82THz2.09THz 频率范围内产生携带 OAM 的涡旋波束，转换效率高于 94.7%，辐射增益大于 10dB，OAM 主模式能量占比最高。对设计的超表面进行了优化仿真，通过增加单元数量以及将超表面设计为圆形均能够有效地提升主模式强度占比，仿真分析了误差对反射涡旋波模场分布的影响。本文设计的超表面具有结构简单、转换效率高、工作带宽大、主模式强度高等优势，对宽带太赫兹涡旋波的产生有一定参考价值，在
44、太赫兹通信、成像等方面具有潜在的应用前景。参考文献：PENGY,SHIC,ZHUY,etal.Terahertzspectroscopyinbiomedicalfield:a review on signal-to-noise ratio im-provementJ.PhotoniX，2020（1）：1-18.1杨书涛,于淼,张孟庭,等.以高斯波束为探针的太赫兹叠层成像系统J.应用光学，2021，42（4）：571-576.YANGShutao,YUMiao,ZHANGMengting,etal.Tera-hertz ptychography system using Gaussian bea
45、m asprobeJ.JournalofAppliedOptics，2021，42（4）：571-576.2WANM,HEALYJJ,SHERIDANJT.TerahertzphaseimagingandbiomedicalapplicationsJ.Optics&LaserTechnology，2020，122：105859.3WANZ,GAOZ,GAOF,etal.TerahertzmassiveMIMOwith holographic reconfigurable intelligent surfacesJ.IEEE Transactions on Communications，2021
46、，69（7）：4732-4750.4杨鸿儒,李宏光.太赫兹波通信技术研究进展J.应用光学，2018，39（1）：12-21.YANG Hongru,LI Hongguang.Research progress onterahertz communication technologyJ.Journal of Ap-pliedOptics，2018，39（1）：12-21.5CHENZ,MAX,ZHANGB,etal.AsurveyonterahertzcommunicationsJ.ChinaCommunications，2019，16（2）：1-35.6ALLENL,BEIJERSBERGEN
47、MW,SPREEUWRJC,etal.Orbitalangularmomentumoflightandthetrans-formationofLaguerre-GaussianlasermodesJ.PhysicalReviewA，1992，45（11）：8185.7HERRINGRA.AnewtwistforelectronbeamsJ.Sci-ence，2011，331：155-156.8290应用光学第45卷第2期YANG T,SHI H,GUO J,et al.Orbital-angular-mo-mentum-basedsuper-resolution ISAR imaging fo
48、r man-euveringtargets:modelingandperformanceanalysisJ.DigitalSignalProcessing，2021，117：103197.9GONGL,ZHAOQ,ZHANGH,etal.Opticalorbital-an-gular-momentum-multiplexed data transmission underhigh scatteringJ.Light:Science&Applications，2019，8（1）：27.10WUJ,HUANGZ,RENX,etal.Widebandmillimeter-wavedual-moded
49、ualcircularlypolarizedOAMantennausing sequentially rotated feeding techniqueJ.IEEEAntennasandWirelessPropagationLetters，2020，19（8）：1296-1300.11CHEN Y,ZHENG S,LI Y,et al.A flat-lensed spiralphaseplatebasedonphase-shiftingsurfaceforgenerationofmillimeter-waveOAMbeamJ.IEEEAntennasandWirelessPropagation
50、Letters，2015，15：1156-1158.12孙胜,阳棂均,沙威.基于反射超表面的偏馈式涡旋波产生装置J.物理学报，2021，70（19）：263-269.SUNSheng,YANGLingjun,SHAWei.Offset-fedvor-tex wave generator based on reflective metasurfaceJ.ActaPhysicaSinica，2021，70（19）：263-269.13樊俊鹏,王小召,张建新,等.一种光可调的多功能太赫兹几何相位超表面J.应用光学，2022，43（5）：986-993.FANJunpeng,WANGXiaozhao,
