维普资讯 http://www.cqvip.com 第29卷第6期 哈尔滨工程大学学报 Vo1.29№.6 2008年6月 Journal of Harbin Engineering University Jun 2008 负折射率介质的电磁特性参数研究 谢德权,关春颖,苑立波 (哈尔滨工程大学理学院,黑龙江哈尔滨150001) 摘要:为了研究负折射率介质的负折射率效应,构造了2种H型的单元结构.基于等效均匀介质理论,对2种H 型的单元结构进行了分析.得到了这2种结构呈现负折射率效应时的电磁特性参数,包括相对波阻抗、折射率、相对 介电常数和相对磁导率.比较得到的2组参数,发现透射率和负折射率频率区间的大小是在负折射率介质构造过程 中需要同时考虑的2个方面. 关键词:均匀介质;H型单元结构;负折射率;电磁特性参数 中图分类号:O441.6文献标识码:A文章编号:1006—7043(2008)06—0622—06 A study on electromagnetic characteristic parameters of mediums with negative refractive index XIE De—quan,GUAN Chun—ying,YUAN Li.bo (College of Science,Harbin Engineering University,Harbin 150001,China) Abstract:In order to study the effect of negative refractive effect of mediums with negative refractive index, two kinds of H—type element structures were constructed.An analysis was conducted on the structures of these tWO elements based on equivalent homogeneous medium theory,obtaining the electromagnetic char— acteristic parameters of these two structures in the case of negative refractive index,including relative wave impedance,refractive index,relative permittivity and relative permeability.Comparing the various parameters of these two elements,we found that the transmissivity and the size of frequency interval are two facets to be simultaneously considered in constructing a negative refractive index medium. Keywords:homogeneous medium;H—type element structures;negative refractive index;electromagnetic characteristic parameters 前苏联科学家Veselago在1968年从理论上提 问题,负折射率值可以通过测量楔形介质的偏转角 出当介质的介电常数和磁导率都为负值时,该介质 得到『3 引,也可以通过波束横移实验来获得E17],还 能够呈现负折射率效应『1],称之为负折射率介质。时 可以引入等效均匀介质的概念来推导获得[1 _21]. 隔30年后,负折射率介质在物理学界引起了广泛关 该文基于文献[193提供的方法,结合有限积分 注.首先从实验上验证了该介质的负折射率效 算法软件CST Microwave Studio(MWS),对构造 应[2_ ,同时对它的研究延伸到许多领域,包括超级 的2种负折射率介质的单元结构进行仿真计算,在 透镜成像的研究[6。。。,微带传输线负折射率效应的 等效均匀介质概念的基础上,获得了这2种结构的 实现[1l-la],负折射率介质波导的导波特性L1 .负 电磁特性参数,包括相对波阻抗、折射率、相对介电 折射率介质的应用是基于它的负折射率,如何获得 常数和相对磁导率的值. 相应频段的负折射率值一直是广大研究者所关心的 等效均匀介质理论 收稿日期:2006—09—04. 作者简介:谢德权(1982一),男,硕士; 苑立波(1962一),男,教授,博士生导师,E-mail:lbyuan@ 当负折射率介质呈现周期性结构特性,与其作 rip.sina.COITL 用的平面电磁波的波长大于6倍周期值时,该介质 维普资讯 http://www.cqvip.com 第6期 谢德权,等:负折射率介质的电磁特性参数研究 ・623・ 可以等效为均匀介质.如果入射平面电磁波在真空 定下来,下面就结合相对介电常数er和相对磁导率 与Z和 的关系式来获得这个条件,e 和 可以表示 为e ===n/Z, 一nZ,用e 、 、 来表示e 和 中的波矢值为k。,均匀介质的厚度为d,折射率和 相对波阻抗设为 和Z,在已知该介质两界面处的 反射和透射系数S 和S。 的情形下,可以得到 ,Z 的实部和虚部,易得 关于k。、d、S 和S。 的关系表达式,文献[2o一21]对 此已有阐述.结合该文,有限积分算法软件MWS采 用的平面波的指数表达式为expEi( ̄一忌・r)],得 到 c一(Ze一1)EL—exp(--i2nk o )] ,1、 一百 sz 一 悉 . 依据文献[18]的推导方法,得到相对波阻抗z的表 为 (3) 由式(1)、(2)得 eXp(一i磁。 )一F .(4) 从式(3)可以看出,相对波阻抗Z有2个根,可以表 示为Z—Z 一if,而它的实部需满足条件Z ≥0, 依此能够获得z的正确值.将z的值代人式(4),得 到表达式exp(一ink。 )的值,为了获得折射率 的 值,将 表示为 — 一in",因为介质有吸收损耗, ≥0,得 exp(一i o )一exp[--i(n ̄ko +2 ̄)]exp(- o ). (5) 式中,m是整数,对式(5)两边求对数得 Ira===一—{lnEexp( - ink)]———广od}+2azn+ —一i—Re{lnEexp(—inkod) ̄}‘ kod . ‘ (6)、 从式(6)可以看出,只要确定了exp(-ink。 )的值,就能 得到∥的值,加上条件 ≥0,可以获得 的合理值. 而 的值由于整数m的存在还需要其他条件才能确 (7) 一 2 z. (8) 因为 ≥0, ≥0 Bg],由式(7)、(8),得到关于 的 不等式: l≤d'z /I I. (9) 关于 的不等式也能得到 ≤ ≤ 丝: j l二 !2 二 n .。(1o) 、… 由不等式(10),整数m的取值范围可以确定,但不 能得到m的准确值,为了解决这个问题,采用文献 [21]的方法,若m在频率厂 处只能取单值,此时的 折射率 (fs)就能被确定,设 ( )和 ( )是 上一个频率^和下一个频率厂 处对应的折射率 值,则 exp[-in( )志。( ) ]≈ exp[一i ( )志。( ) ](1+△1+ 1凸21). (11) 其中,A 一in( )志。( ) —in( )志。( ) ,得到 ( )的表达式: ( )一f,n(is)/f.一 (12) ic(1干、 日)/(27【 ). 同理, ( )的表达式如下 ic(1干 可)/(2 ̄fbd). (13) 其中,△2===in( )志。( ) —in( )志。( ) .由式 维普资讯 http://www.cqvip.com 哈尔滨工程大学学报 第29卷 (4)、(12)和(13),应用迭代法可以求出各个频率处 对应的折射率值.因为计算结果包含2个分支,所以 必须以式(6)得到的 值为参考,取出合适的折射 率值,表示为 。= 。E1 。n a,。 。和 分别是 。的 实部和虚部.令 一 。,得 一 .(14) 通过式(4)、(14),单一的m值被求得,由于m是整 数,所以对式(14)的结果求整,得到各个频率处对应 的单一的m的整数值.因为获得m值的过程都是近 似的,所以有必要对它进行检验,将m值代人表达 式(6),求出 的值,然后代人不等式(9),看是否满 ’足,从而得到正确的m值. 下面就用上面阐述的等效均匀介质理论对构造 的2种负折射率介质的单元结构来进行分析,并获 取这2种结构的电磁特性参数,包括Z、 、e, 和 ,. 2 负折射率介质的电磁特性参数 2.1 H型结构 , 构造的H型的单元结构如图1所示,在电路板 -4 r 的两面镀上铜的涂覆层,电路板采用FR4型,相对 介电常数为4.4,损耗正切值为0.02,铜的电导率是 5.8×10 S/m,涂覆层的厚度是0.02 mm,宽度是 0.2 mm.图1中的黑色部分表示H型的磁谐振结 构,高度为h 一2.2 mm,h。一4.8 mm,宽度为 w一2.6 mm,开口大小是g一0.4 mm;灰黑色部 分是条状的铜导体,高度为10 mm;中间半透明部 分是FR4型的电路板,沿5C轴方向的高度为 10 mm,沿Y轴方向的厚度为0.25 mm,沿 轴方向 的长度为3.1 mm.应用MWS软件对该结构进行仿 真时,设人射平面电磁波的频率区间为3 OHz到 20 GHz,沿 轴方向传播,电场偏振方向沿z轴,磁 场偏振方向沿Y轴,单元结构沿z、Y、 方向的尺寸 分别为10 mm、2.5 mm、3.1 mm.此时将该结构等 效为均匀介质,厚度d一3.1 mm,在边界处依次满 足电场,磁场和开放边界条件.图2(a)所示为仿真 后的反射和透射系数的振幅部分,用l S。。l和 I s。 I表示;图2(b)所示为仿真后的反射和透射系 数的位相部分,用 和 z表示. 图1 H型的单元结构图 Fig.1 H—patterned unit structure (a)幅值』S 』和』52 』关于频率的曲线 --- 一92 -: : l l0 20 lfOHz (b)相位 和 z关于频率的曲线 图2 H型单元结构的反射系数和透射系数 Fig.2 Reflectance and transmissivity of H—patterned unit structure 反射系数的复振幅表示为S 一 l S l exp(i ̄  ̄/180),透射系数的复振幅表示为 S2 一l S2 l exp( 7/180),代人表达式(3),应用条 件Z ≥0,求出相对波阻抗Z的正确值,代人表达 式(4),得到exp(一ink。 )的值.由不等式(10),整 维普资讯 http://www.cqvip.com 第6期 谢德权,等:负折射率介质的电磁特性参数研究 数m的取值范围可以确定,发现m在频率 等于 5.04 GHz处只能取值0,代入表达式(12)和(13)得 到 左右两侧各个频率点对应的折射率值,与式 (6)得到的 值比较,取出合适的折射率值7/" ,由 式(14),单一的m值被求得,求整后代入式(6),将 求得的7"/值代入不等式(9)进行验证,去除不合理的 m值,得到满足式(9)的m值对应的频率区间为 [5.04 GHz,7.981 GHZ],E9.715 GHz,12.146 GHZ], 分析该区间内的7/" 、e, 和 , 值,得到它们同时取负 值时的负折射率效应频率区间:[9.715 GHz, 12.146 GHz-].此时12.146 GHz处的波长和均匀 介质厚度d的比值大于7,说明该结构等效为均匀 介质是合理的.图3所示为H型单元结构呈现负折 射率效应时的电磁特性参数Z、,z、e 和 ,关于频率 的曲线. (a)相对波阻抗Z (b)折射率 (c)相对介电常数£ (d)相对磁导率 图3 H型单元结构的电磁特性参数关于频率的曲线 Fig.3 Curves of electromagnetic characterized parameters of H—patterned unit structure versus frequency 2.2改进的H型结构 改进的H型结构如图4所示,与图1不同的地方 在于磁谐振结构由4个开口铜环构成,参数分别是: h1===3.575 inrn,h2—4.5 inrn,h3—1.125 mm,w一 2.45 mm,g一0.4 mm.此时单元结构沿z方向的尺 寸为3.25 mm,即d一3.25 mm. 与H型的单元结构的仿真过程类似,得到反射 和透射系数的振幅和位相部分关于频率的曲线如图 5所示,呈现负折射率效应时的电磁特性参数Z、7"/、 e,和 ,关于频率的曲线如图6所示,获得的负折射 率效应频率区间为E13.098 Gt-lz,14.322 GHZ], 14.322 GI-Iz ̄的波长和均匀介质厚度d的比值大于6, 说明该结构等效为均匀介质也是合理的. 图4改进的H型的单元结构图 Fig.4 Improved H—patterned unit structure 对图2和图5,图3和图6进行比较,发现改进 的H型结构的透射率要大,而负折射率效应频率区 间要小,所以在负折射率介质构造的过程中,透射率 维普资讯 http://www.cqvip.com 哈尔滨工程大学学报 第29卷 和负折射率效应频率区间要综合考虑,依据不同的 应用可以有所侧重. (a)幅值l S l和l S: l关于频率的曲线 e 、 (b)相位1和 z关于频率的曲线 图5改进的H型单元结构的反射系数和透射系数 Fig.5 Reflectance and transmissivity of improved H—patterned unit structure (a)相对波阻抗Z (b)折射率 lfGHz (c)相对介电常数e lfGHz (d)相对磁导率 图6改进的H型单元结构的电磁 特性参数关于频率的曲线 Fig.6 Curves of electromagnetic characterized parameters of impmved H-patterned unit structure versus frequency 3 结束语 该文在文献[19]的基础上,详细阐述了等效均 匀介质理论的分析方法,对如何获得公式(6)中整数 m的合理值作了详尽的推导,并应用该方法对构造 的2种H型的单元结构进行了分析,得到了这2种 结构呈现负折射率效应时的电磁特性参数Z、 、e, 和 ,关于频率的曲线,确定了负折射率效应的频率 区间.在实际构造负折射率介质的过程中,以图3和 图6的参数为基础,在CST Design Studio(Ds)软 件环境下,指定单元结构排列的周期数,就能得到该 周期介质的负折射率特性参数,为制作实物和实验 测量提供了参考依据. 参考文献: [1]VESELAGO V G.The electrodynamics of substances with simultaneously negative values of£and [J].Soy Phys Usp,1968,10(4):509—514.  ̄2]SHELBY R A,SMITH D R,SCHULTZ S.Experimen— 维普资讯 http://www.cqvip.com 第6期 谢德权,等:负折射率介质的电磁特性参数研究 ・ 627・ tal verification of a negative index of refraction[J].SCi— ence,2001,292:77—79. [3]PARAZZOLI C G,GREEC ̄R R B,LI K,et a1.Exper— imental verification and simulation of negative index of refraction using snell’s law[J].Phys Rev Lett,2003,90 (10):107401. [4]HOUCK A A。BROCK J B,CHUANG I L.Experimen— tal observation of a left-handed material that obeys snell’s law[J].Phys Rev Lett,2003,90(13):137401. [5]GREEGOR R B,PARAZZOLI C G,LI K,et a1.Exper— imental determination and numerical simulation of the properties of negative index of refraction materials[J]. Opt Express,2003,11(7):688—695. [61 PENDRY J B Negative refraction makes a perfect lens [J].Phys Rev Lett,2000,85(18):3966—3969. [7]FANG N。LEE H,SUN C,et a1.Sub-diffraction-limit— ed optical imaging with a silver superlens[J].Science, 2005,308:534—537. r8]BRATKOVSKY A M,CANO A,LEVANYUK A P. Strong effect of surfaces on resolution limit of negative- index“superlens”[J].Appl Phys Lett,2005,87: 103507. [9]PODOLSKIY V A,KUHTA N A。MILT0N G w.Op— timizing the superlens.manipulating geometry to en— hance the resolution[J].Appl Phys Lett,2005,87: 231113. [10]WEI Jingsong,GAN Fuxi.Dynamic readout of subdif— fraction-limited pit arrays with a silver superlens[J]. Appl Phys Lett,2005,87:211101. [11]GRBIC A,ELEFTHERIADES G v.Growing evanes— cent waves in negative-refractive-index transmission-line media[J].Appl Phys Lett,2003,82(12):1815—1817. [12]IYER A K。KREMER P C,ELEFTHERIADES G v. Experimental and theoretical verification of focusing in a large,periodically loaded transmission line negative re— fractive index metamaterial[J].Opt Express,2003,11 (7):696—708. [131 GRBIC A.ELEFTHERIADEs G V.Overcoming the diffraction limit with a planar left-handed transmission- line lens[J].Phys Rev Lett,2004,92(11):117403. [14]WU B L.GRZECRDRCzYK T M,zHANG Y,et a1. Guided modes with imaginary transverse wave number in a slab waveguide with negative permittivity and permea— bility[J].J Appl Phys,2003,93(11):9386—9388. [15]AGUANNO G D。MATTIUCCI N。SCALORA M,et al_TE and TM guided modes in an air waveguide with negative-index-material cladding[J].Phys Rev E,2005, 71:046603. [16]DONG Z G。ZHU S N,LIU H,et a1.Numerical simu— lation of negative-index refraction in wedge-shaped meta— materials[J].Phys Rev E,2005,72:016607. [17]RAN Lixin,HUANGFU Jiangtao,CHEN Hongsheng, et a1. Beam shifting experiment for the characterization of left-handed properties[J].J Appl Phys,2004,95(5): 2238—2241. r18]SMITH D R。SCHULTz S,MARKOS P,et a1.Deter— mination of effective permittivity and permeability of metamaterials from reflection and transmission coeffi— cients[J].Phys Rev B,2002,65:195104. [19]CHEN Xudong,GRZEG0RCzYK T M,wU B,et a1. Robust method to retrieve the constitutive effective pa— rameters of metamaterials[J].Phys Rev E,2004,70: 016608. [20]SMITH D R。VIER D C,KOSCHNY T,et a1.Elec— tromagnetic parameter retrieval from inhomogeneous metamaterials[J].Phys Rev E,2005,71:036617. [21]POPA B I.CUMMER S八Determination the effective electromagnetic properties of negative-refractive-index metamaterials from internal fields[J].Phys Rev B, 2005,72:165102.