-
具有螺旋波前的电磁波是携带轨道角动量的涡旋波束, 涡旋波束存在的相位奇点, 使其在微粒操控和通讯等领域有特殊的应用. 本文提出了一种基于反射型超表面的太赫兹宽带涡旋波束产生器, 该器件由超表面-电介质-金属三层结构构成, 顶层为两个正交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.
-
Keywords:
- terahertz /
- vortex beam /
- broadband /
- metasurface
1. 引 言
近年来, 非晶氧化锌镓铟薄膜晶体管(amorphous InGaZnO thin-film transistor, InGaZnO TFT)的研究及其应用广受关注. 由于具有高迁移率、低关态电流、高透光率、均匀性良好以及低温可制备性等优点[1-5], InGaZnO TFT有望取代硅基TFT, 在平板显示、光学图像传感器、触控和指纹传感等领域拥有良好的应用前景[6-9]. 在传统的应用中, 例如薄膜晶体管液晶显示(TFT-LCD), 或者有源矩阵有机发光二极管显示(active-matrix organic light emitting display, AM-OLED), 主要利用的是InGaZnO TFT导通态的电学性能[10-12], 但在InGaZnO TFT集成的众多新兴领域中, 例如光学图像传感器, InGaZnO TFT的泄漏电流特性至关重要.
以往的研究表明, TFT中常见的泄漏电流产生机制包括陷阱辅助热发射[13,14]、陷阱辅助场发射[13]、带间隧穿[15,16]和包含Poole-Frenkel (PF)效应的陷阱辅助热电子场致发射[17]等. Lui和Migliorato [18]基于PF效应和陷阱辅助隧穿效应, 提出了载流子产生复合速率模型, 其适合于器件的数值计算. Wu等[19]在此基础上提出了多晶硅TFT泄漏电流紧凑模型, 其中包含了PF效应, 无需繁琐的数值计算, 并且可以通过实验数据提取得到模型参数. Servati和Nathan [20]研究了非晶硅(a-Si:H) TFT的泄漏电流的产生来源, 基于欧姆传导模型、反向亚阈传导模型和前沟道传导模型建立了分段式泄漏电流模型. 但是迄今为止, 几乎未见有关InGaZnO TFT的泄漏电流模型的研究报道, 这不利于InGaZnO TFT集成图像传感器等新兴领域的研究.
本文基于载流子的产生-复合机制建立了InGaZnO TFT的泄漏电流模型. 通过与TCAD模拟以及测试结果的对比来检验所提出的模型. 基于所提出的泄漏电流模型, 分析了InGaZnO TFT的沟道宽度、沟道长度和栅介质层厚度对泄漏电流的影响.
2. 基于TCAD的InGaZnO TFT泄漏电流的数值分析
图1示意了底栅顶接触型InGaZnO TFT的剖面图. 表1列出了InGaZnO TFT器件结构的几何参数, 其中 InGaZnO TFT的沟道宽长比W/L = 300
μm /50μm , a-InGaZnO半导体层厚度tIGZO = 40 nm, 栅介质层厚度tSiOx=150nm [21].表 1 InGaZnO TFT器件结构的几何参数Table 1. Geometric parameters of InGaZnO TFT device structure.参数 数值 栅介质层厚度/nm 150 a-InGaZnO半导体层厚度/nm 40 沟道宽度/μm 300 沟道长度/μm 50 图2(a)给出了在实验[21]与TCAD数值模拟中得到的InGaZnO TFT在关断区和亚阈区的IDS-VGS转移特性的曲线对比, 其中VDS取值从2 V增加到10 V. 由图可知, InGaZnO TFT具有较低的泄漏电流, 且IDS随着VDS增大而增加, IDS的变化范围为10–13—10–12 A. 泄漏电流增加的主要原因是随着TFT上所施加偏压的增加, 陷阱态中的电子获得能量跃迁至导带的几率就越大. 实验值和TCAD模拟值的变化趋势及范围大致相同, 在亚阈区内泄漏电流均呈数量级的增长.
图2(b)为同一InGaZnO TFT IDS-VDS转移特性曲线的实验值和TCAD模拟值的对比, 结果表明IDS和VGS, VDS均呈正相关, 在器件的饱和区(VDS ≥ VGS-VTH), IDS的值趋于稳定且不随VDS发生变化. 从图中可以看出实验值与TCAD模拟值吻合度较高, 说明TCAD模拟较为可靠. 因此, 后续将以TCAD模拟值作为参考, 验证所建立InGaZnO TFT泄漏电流解析模型的可靠性.
在InGaZnO TFT的TCAD电学仿真中, 主要通过金属功函数和UST (universal Schottky tunnling model)模型描述肖特基势垒; 有源层材料InGaZnO的特性主要由缺陷态(defects)密度模型决定. InGaZnO的缺陷态模型与泄漏电流密切相关, 具体模型参数见表2.
表 2 InGaZnO缺陷态密度模型参数Table 2. Density of states model parameters for InGaZnO.参数 描述 数值 单位 nta 导带尾类受主能态密度 1.04 × 1019 cm–3·eV–1 ntd 价带尾类施主能态密度 5.0 × 1020 cm–3·eV–1 wta 类受主态特征能量 0.04 eV wtd 类施主态特征能量 0.1 eV nga 高斯分布的受主态密度 0 cm–3·eV–1 ngd 高斯分布的施主态密度 2.0 × 1016 cm–3·eV–1 egd 高斯分布施主能态峰值能量 2.9 eV wgd 高斯分布施主能态特征能量 0.1 eV 3. 泄漏电流的紧凑模型建立与分析
根据陷阱辅助热电子场致发射机制, 在较强的栅-漏电场作用下, 大量电子将从局域态中隧穿到导带[18]. 载流子的产生率与Poole-Frenkel效应和声子辅助隧穿效应密切相关. InGaZnO层带隙内的局域态可表示为
RT=RA+RD, (1) 其中RT是陷阱态中的载流子产生复合率, RA和RD分别是受主态和施主态的载流子产生复合率. 作为n沟道器件, InGaZnO TFT中类施主态呈中性, 对自由载流子数量的影响几乎可以忽略. 但是, InGaZnO的类受主态呈正电性, 可捕获自由电子, 对泄漏电流的大小的影响较大. 在建模过程中, 为了简化计算, 忽略了施主态, 只考虑了受主态的载流子产生. 因此(1)式可简化为:
RT=RA, (2) RA=np−n2in+n1cp(χF+ΓCoulp)+p+p1cn(1+ΓDiracn)NA(Et), (3) (3)式中NA(Et)是类受主态在导(价)带附近处深能态和带尾态的密度, 受主态包括尾态和深能态级, 由于InGaZnO TFT中深能态级密度远低于带尾态, 因此, 计算过程中主要考虑带尾态密度;
χF 是Poole-Frenkel增强效应因子; ni是本征载流子浓度; n1和p1分别代表陷阱能级上的电子浓度和空穴浓度, cn,p代表俘获系数, 且n1=niexp(Et−EikT), (4) p1=niexp(Ei−EtkT), (5) χF=exp(ΔECkT), (6) 式中k为玻尔兹曼常数, T为温度; 其中
ΔEC 可表征为ΔEC=q√qFπεIGZO, (7) 其中q代表电荷量, F为垂直方向的电场强度,
εIGZO 为InGaZnO材料的介电常数. (7)式来源于描述Poole-Frenkel效应下泄漏电流与电场强度关系的公式[22], 在Poole-Frenkel效应作用下, 泄漏电流与电场强度的指数正相关,ΔEC 表示在电场强度F作用下势垒的降低量. (3)式中ΓpCoul与ΓDiracn 分别是Dirac阱与Coulombic阱中载流子的隧穿比率, 可表征为:ΓDiracn=ΔEnkT∫10exp(ΔEnkT−knu32)du, (8) ΓCoulp=ΔEpkT∫1ΔECΔEpexp{ΔEpkTu−kpu32[1−(ΔECuΔEp)53]}du, (9) ΔEn=EC−Et, (10) ΔEp=Et−EV, (11) kn,p=4√2m∗e,h(ΔEn,p)33qℏF. (12) (12)式中
m∗e,h 表示隧穿效应中电子(空穴)的有效质量,ℏ 表示普朗克常数.NA(Et)=NTailAexp[Et−ECETailA], (13) U=∫ECEVRT(Et)dEt=∫ECEVRA(Et)dEt. (14) 将(1)—(13)式代入(14)式中, 可以求解得到InGaZnO TFT的载流子产生率. 在上述模型中, 包含了InGaZnO TFT可能的泄漏电流产生机制: 场致电子发射、带间隧穿以及Poole-Frenkel效应[23-25]. 由于InGaZnO TFT器件一般工作于常温, 热电子发射产生的泄漏电流极小, 这里忽略了热电子发射效应. 同时, (14)式也是一双积分表达式, 计算过程繁琐且需要消耗大量的内存资源. 因此, 本文采用分区段计算的方法来简化建模过程.
分析式(14)可知, 第一重积分计算的是场效应积分即
ΓCouln,p和ΓDiracn,p 的值, 积分对象是u, 积分过程与Poole-Frenkel增强因子χF 无关; 第二重积分的对象是Et, 积分过程同样与Poole-Frenkel增强因子χF 无关; 在整个积分过程中Poole-Frenkel增强因子χF 始终保持指数形式, 在最终积分结果中也保持原状.在低负电场区域(F ≤ FC), Poole-Frenkel效应增强下的热发射是泄漏电流产生的主要机制, 在载流子产生复合率U的表达式里Poole-Frenkel增强因子
χF 的值占主导地位, 因此低负电场区载流子产生复合率U可简化成U=a1exp(b1√F),F⩽FC. (15) 在高负电场区域(F ≥ FC), 陷阱中的载流子首先热激发到一个中间能级Et, 然后再从该中间能级隧穿到导带或者价带; 根据文献[26]可知在高负电场区域热离子场助发射是产生泄漏电流的主要原因. 文献[26]中计算得出, 在高负电场区域lnU与电场强度F呈线性关系, 故在高负电场区载流子产生复合率U可表示为
U=a2exp(b2⋅F),F>FC. (16) U={a1exp(b1√F),F⩽FC,a2exp(b2F),F>FC, (17) 其中a1, a2, b1, b2均为拟合参数. 泄漏电流ILeak可由(18)式描述[27]:
ILeak=qVolU, (18) 其中
Vol=W⋅Lov⋅Xd, (19) 式中q是电荷量, Vol是源漏极交叠耗尽区的体积, W是InGaZnO TFT沟道宽度, Lov是栅漏交叠区长度, Xd为耗尽区厚度. 将(19)式代入(18)式中得到关断区泄漏电流模型:
ILeak-off={qVol⋅a1exp(b1√F),F⩽FC,qVol⋅a2exp(b2F),F>FC, (20) 其中
F=VDS−VGStSiOx, (21) tSiOx 为栅氧化层厚度.由图2(a)可知在亚阈区后半段lnIDS与VGS呈线性关系:
ILeak-sub=a3exp(b3VGS). (22) 使用平滑函数tanh(x)[28]得到InGaZnO TFT泄漏电流的统一模型为
ILeak=qWLovXd⋅[a1exp[b1(VDS−VGStSiOx)]×{1−tanh[β(VGS−VC)]2}+a2exp(b2√VDS−VGStSiOx)×(1+tanh[β(VGS−VC)]2)]×(1−tanh(βVGS)2)+a3exp(b3VGS)⋅[1+tanh(βVGS)2], (23) 其中a1, a2, a3, b1, b2, b3均为拟合参数.
根据上述模型可得到InGaZnO TFT泄漏电流的计算值, 此外还可利用TCAD模拟的方式得到InGaZnO TFT泄漏电流的模拟值. 图3是模型计算结果和TCAD模拟结果的对比图. 从模拟结果来看, VDS从2 V增大到10 V的过程中, InGaZnO TFT泄漏电流也随之增大. 这与(20)式和(21)式描述的情况相符, VDS增大导致垂直方向电场F增大, 在强电场的作用力下, Poole-Frenkel效应愈加明显, 在陷阱态的电子只需要获得较小的热能便能脱离陷阱态变成自由载流子, 沟道内自由载流子增多使得泄漏电流增大.
图4给出了不同沟道宽度情况下的InGaZnO TFT的TCAD模拟结果和模型计算结果的对比. TFT的沟道宽度W从200
μm 增大至500μm , 泄漏电流随TFT的沟道宽度线性增大. 这与(19)式描述的情况较为符合. 这是因为TFT的沟道宽度W增大导致耗尽区的面积和体积增大, 从而沟道内的感应载流子增多, 致使泄漏电流增大. 由于(19)式中的Lov是指栅漏交叠区域的长度, 因此, TFT的沟道长度与泄漏电流的大小并无关系[29]. 由图5可知, TFT的沟道长度L从50μm 增大至90μm , 在关断区域内漏电流大小增长幅度仅为0.6%. 然而, 在导通区域, 漏电流随着TFT沟道长度L的增大而减小. 该模拟结果表明TFT的沟道长度L对关断区泄漏电流几乎无影响.图6为InGaZnO TFT不同栅氧化层厚度
tSiOx 下模型计算和TCAD模拟的泄漏电流对比图. 可见, 泄漏电流随着栅氧化层厚度tSiOx 的增大而减小, 栅氧化层厚度tSiOx 从150 nm增大至250 nm过程中, 泄漏电流的值减小了20%. 这与(21)式和(23)式描述的情况一致, 栅氧化层厚度tSiOx 增加导致电场强度F减弱, 陷阱态内电子无法获得足够的动力挣脱束缚, 沟道内自由载流子减少, 载流子复合产生率减小. 且在高电场作用下, 栅氧化层厚度tSiOx 的增加也会导致器件内隧穿电流减小, 最终亦会导致泄漏电流减小.4. 结 论
本文基于载流子产生-复合机制提出了一种InGaZnO TFT的泄漏电流模型. 基于低负电场区的Poole-Frenkel增强热发射、高负电场的场辅助热发射的泄漏电流产生机制, 在不同电场条件下提出了双积分载流子产生-复合率的近似模型, 并利用平滑函数将其统一成适用于关断区和亚阈区的泄漏电流模型. 将模型计算结果与TCAD模拟结果进行了对比, 在器件的亚阈区和关断区, 泄漏电流IDS的模型计算值和TCAD模拟值吻合程度较高. 基于上述模型, 本文讨论了InGaZnO TFT的沟道宽度、沟道长度和栅介质层厚度等关键参数对泄漏电流的影响. 模型计算值与TCAD模拟值的对比结果证明该模型在预测InGaZnO TFT沟道宽度、沟道长度和栅介质层厚度等结构参数对泄漏电流的影响方面较为可靠.
[1] Beard M C, Turner G M, Schmuttenmaer C A 2002 J. Phys. Chem. B 106 7146
Google 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 28249
Google Scholar
[3] Janek M, Zich D, Naftaly M 2014 Mater. Chem. Phys. 145 278
Google Scholar
[4] Hu B B, Nuss M C 1995 Opt. Lett. 20 1716
Google Scholar
[5] Mittleman D M, Gupta M, Neelamani R, Baraniuk R G, Rudd J V, Koch M 1999 Appl. Phys. B 68 1085
Google 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 156
Google 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 124202
Google 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 2684
Google Scholar
[9] Kirilenko M S, Khonina S N 2013 Optical Memory and Neural Networks 22 81
Google 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 S52
Google Scholar
[11] Genevet P, Lin J, Kats M A, Capasso F 2012 Nature Comm. 3 1
Google 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 565
Google 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 087701
Google 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 0314001
Google Scholar
Li Y H, Zhou L, Zhao G Z 2018 Chin. J. Lasers 45 0314001
Google Scholar
[16] Li H, Xiao B Y, Huang X J, Yang H L 2015 Phys. Scr. 90 1
Google Scholar
[17] Enoch S, Tayeb G, Sabouroux P, Guerin N, Vincent P 2002 Phys. Rev. Lett. 89 213902
Google Scholar
[18] Huang J, Pogorzelski R J 1998 IEEE Trans. Antenn. Propag. 46 650
Google Scholar
[19] Martynyuk A E, Martinez-Lopez J I, Martynyuk N A 2004 IEEE Trans. Antenn. Propag. 52 142
Google 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 013101
Google Scholar
[22] Yu N, Genevet P, Kats M A, Aieta F, Tetienne J P, Capasso F, Gaburro Z 2011 Science 334 333
Google 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 1351
Google Scholar
[24] Luo W J, Sun S L, Xu H X, He Q, Zhou L 2017 Phys. Rev. Appl. 7 044033
Google Scholar
[25] Xu H H, Wang G M, Cai T, Xiao J, Zhuang Y Q 2016 Opt. Express 24 27836
Google Scholar
-
图 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.
图 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 7146
Google 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 28249
Google Scholar
[3] Janek M, Zich D, Naftaly M 2014 Mater. Chem. Phys. 145 278
Google Scholar
[4] Hu B B, Nuss M C 1995 Opt. Lett. 20 1716
Google Scholar
[5] Mittleman D M, Gupta M, Neelamani R, Baraniuk R G, Rudd J V, Koch M 1999 Appl. Phys. B 68 1085
Google 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 156
Google 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 124202
Google 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 2684
Google Scholar
[9] Kirilenko M S, Khonina S N 2013 Optical Memory and Neural Networks 22 81
Google 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 S52
Google Scholar
[11] Genevet P, Lin J, Kats M A, Capasso F 2012 Nature Comm. 3 1
Google 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 565
Google 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 087701
Google 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 0314001
Google Scholar
Li Y H, Zhou L, Zhao G Z 2018 Chin. J. Lasers 45 0314001
Google Scholar
[16] Li H, Xiao B Y, Huang X J, Yang H L 2015 Phys. Scr. 90 1
Google Scholar
[17] Enoch S, Tayeb G, Sabouroux P, Guerin N, Vincent P 2002 Phys. Rev. Lett. 89 213902
Google Scholar
[18] Huang J, Pogorzelski R J 1998 IEEE Trans. Antenn. Propag. 46 650
Google Scholar
[19] Martynyuk A E, Martinez-Lopez J I, Martynyuk N A 2004 IEEE Trans. Antenn. Propag. 52 142
Google 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 013101
Google Scholar
[22] Yu N, Genevet P, Kats M A, Aieta F, Tetienne J P, Capasso F, Gaburro Z 2011 Science 334 333
Google 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 1351
Google Scholar
[24] Luo W J, Sun S L, Xu H X, He Q, Zhou L 2017 Phys. Rev. Appl. 7 044033
Google Scholar
[25] Xu H H, Wang G M, Cai T, Xiao J, Zhuang Y Q 2016 Opt. Express 24 27836
Google Scholar
计量
- 文章访问数: 14595
- PDF下载量: 504