周曉霞 陳英 蔡力

1) (長(zhǎng)沙學(xué)院電子信息與電氣工程學(xué)院,長(zhǎng)沙 410000)

2) (國(guó)防科技大學(xué),智能科學(xué)學(xué)院裝備綜合保障技術(shù)重點(diǎn)實(shí)驗(yàn)室,長(zhǎng)沙 410073)

本文提出一種基于零折射率介質(zhì)的超窄帶光學(xué)濾波器.在線缺陷光學(xué)濾波器中引入具有類狄拉克點(diǎn)的光學(xué)超構(gòu)材料,利用其零折射效應(yīng)來實(shí)現(xiàn)濾波帶寬的壓縮.基于COMSOL Multiphysics 軟件的傳輸特性分析表明,當(dāng)超構(gòu)材料的類狄拉克點(diǎn)頻率與缺陷態(tài)的諧振頻率相符合時(shí),光學(xué)濾波器的透射峰能夠顯著壓縮;光場(chǎng)分布及等效介質(zhì)分析表明,零折射率介質(zhì)的零相位延遲效應(yīng)與強(qiáng)色散特性能夠增強(qiáng)缺陷態(tài)透射電磁響應(yīng)隨頻率變化的靈敏度,提高濾波器的品質(zhì)因子,并能在壓縮濾波帶寬的同時(shí)保持高的峰值透射率,實(shí)現(xiàn)高耦合、超窄帶的濾波設(shè)計(jì).該結(jié)果為基于光學(xué)超構(gòu)材料的波分復(fù)用系統(tǒng)設(shè)計(jì)與應(yīng)用提供了新的技術(shù)思路.

1 引言

窄帶光學(xué)濾波器件在密集波分復(fù)用(DWDM)[1]、激光雷達(dá)[2]等光通信與光信息處理系統(tǒng)及熒光顯微鏡、流式細(xì)胞儀、DNA 測(cè)序等精密激光儀器設(shè)備[3]中具有廣泛的應(yīng)用前景.利用光學(xué)超構(gòu)材料(photonic metamaterials,PMs)的帶隙特性[4,5]和缺陷態(tài)光子局域特性[6],開發(fā)在相對(duì)較寬的光子禁帶中具有高透過率和狹窄傳輸頻段的光濾波器,是超構(gòu)材料應(yīng)用探索研究廣為關(guān)注的問題[7-9].這些濾波器通常都采用微腔(諧振腔)結(jié)構(gòu),其品質(zhì)因子取決于諧振腔的結(jié)構(gòu)設(shè)計(jì).研究者發(fā)展了環(huán)型腔[10,11]、多諧振腔[12]、諧振腔+反射壁[13]、相移波導(dǎo)光柵結(jié)構(gòu)[14]等多種超構(gòu)材料濾波器模型,濾波帶寬能夠?qū)崿F(xiàn)20 nm 的信道間隔,達(dá)到粗波分復(fù)用系統(tǒng)的要求.但隨著5G、6G、物聯(lián)網(wǎng)、社區(qū)智能化等技術(shù)對(duì)大容量、高密度信息傳輸需求的迅速增長(zhǎng),需要DWDM 系統(tǒng)的復(fù)用信道數(shù)量進(jìn)一步增加以提高對(duì)光纖的利用率,使得信道間隔變得越來越窄.作為提高信噪比、保證光通信質(zhì)量的關(guān)鍵器件,超窄帶濾波器(半功率帶寬＜0.05 nm)的研究引起越來越多的重視[15].同時(shí),它對(duì)激光雷達(dá)、流式細(xì)胞儀等精密光學(xué)設(shè)備抑制背景干擾,提高分辨率和探測(cè)距離也具有重要意義.目前的光學(xué)超構(gòu)材料濾波器還難以同時(shí)滿足高透過率、超窄帶寬的要求.

作為光學(xué)超構(gòu)材料的一個(gè)重要分支,零折射率超構(gòu)材料(zero-index media,ZIM)由于具有高傳播相速度、零相位延遲及高透過率等特性,成為近年來研究的熱點(diǎn)[16-21].基于這些特性,ZIM 能夠產(chǎn)生能量隧穿[22,23]、高指向性輻射[24,25]、無反射彎曲波導(dǎo)[26,27]等獨(dú)特的電磁波操控效應(yīng),為新型光學(xué)器件設(shè)計(jì)提供了新穎的思路.由于ZIM 的等效折射率通常在狹窄的頻段內(nèi)出現(xiàn)零折射率特性,并在鄰近頻率范圍劇烈變化,表現(xiàn)出強(qiáng)色散特性[28,29].這一特性有利于實(shí)現(xiàn)狹窄頻段的波調(diào)控,基于此本文設(shè)想將ZIM 引入超窄帶濾波器設(shè)計(jì)中,在含缺陷態(tài)的超構(gòu)材料設(shè)計(jì)模型中嵌入具有狄拉克椎光子能帶特性的ZIM,通過調(diào)節(jié)尺寸和材料參數(shù)使得ZIM 的狄拉克點(diǎn)頻率位于光子禁帶中,且與缺陷態(tài)的諧振頻率相吻合.COMSOL 仿真結(jié)果表明,提出的設(shè)計(jì)模型在光子禁帶中產(chǎn)生半峰高寬遠(yuǎn)小于僅僅利用缺陷態(tài)局域效應(yīng)實(shí)現(xiàn)的窄帶透射峰,本文進(jìn)而研究了ZIM 對(duì)光學(xué)濾波器品質(zhì)因子及峰值透射率的調(diào)節(jié)規(guī)律和機(jī)理,結(jié)果表明濾波器能夠在壓縮濾波帶寬的同時(shí)保持高的峰值透過率.本研究工作為光學(xué)超構(gòu)材料的超窄帶光濾波器設(shè)計(jì)提供了一種新的思路.

2 理論模型

設(shè)想將基于類狄拉克點(diǎn)的零折射超構(gòu)材料與基于缺陷態(tài)能量局域效應(yīng)的光學(xué)超構(gòu)材料濾波器設(shè)計(jì)相結(jié)合.圖1(a)為二維光學(xué)超構(gòu)材料濾波器原理模型,超構(gòu)材料PMs1 由介電圓柱按正方晶格周期性排列于空氣中構(gòu)成,通過調(diào)節(jié)圓柱半徑和晶格常數(shù)在需要的濾波頻段產(chǎn)生光子禁帶[10].將其內(nèi)部第n列和n+1 列單元的間距拉開L1引入缺陷模態(tài),當(dāng)缺陷態(tài)的特征頻率ω0位于禁帶中時(shí),能夠?qū)⒐鈭?chǎng)局域在缺陷區(qū)域內(nèi)形成光學(xué)諧振腔,產(chǎn)生窄帶濾波效果[30].提高品質(zhì)因子,實(shí)現(xiàn)極窄帶寬和峰值高透射率是濾波器應(yīng)用于光通訊、傳感等工程領(lǐng)域的發(fā)展方向.本文設(shè)想在上述濾波器模型中引入ZIM,在缺陷中嵌入厚度為L(zhǎng)2的ZIM,如圖1(b)所示.ZIM 也由正方晶格周期性排列于空氣中的介電圓柱陣列構(gòu)成,設(shè)計(jì)中保持晶格常數(shù)a為圓柱半徑R的5 倍,即a=5R.對(duì)于TM 波(橫磁極化波,電場(chǎng)沿著柱子的軸向方向)極化,這時(shí)該超構(gòu)材料在Bloch 波矢k=0 處的Γ點(diǎn)位置存在一個(gè)由單極子模態(tài)和偶極子模態(tài)偶然簡(jiǎn)并形成的三重簡(jiǎn)并態(tài)[17,31,32],該模態(tài)附近的光子能帶曲線具有線性特性并在Γ點(diǎn)相交[33],形成狄拉克錐形能帶結(jié)構(gòu),Γ點(diǎn)為類狄拉克點(diǎn).

圖1 濾波器結(jié)構(gòu)示意圖 (a)線缺陷濾波器;(b)線缺陷+零折射超構(gòu)材料濾波器Fig.1.Schematic diagram of the filter: (a) Structure of the filter with line defect;(b) structure of the filter with line defect and zero index metamaterial.

電磁波在該模型諧振腔中傳播的相位延遲包括在背景介質(zhì)中傳播L1的相位延遲Δδ1和ZIM 中傳播L2的相位延遲Δδ2兩部分.在類狄拉克點(diǎn)頻率ZIM 可等效為介電常數(shù)和磁導(dǎo)率均為零的均勻介質(zhì),其中的電磁波傳播具有零相位延遲特性.基于此,考慮通過能帶結(jié)構(gòu)分析,設(shè)計(jì)ZIM 的類狄拉克點(diǎn)頻率ωd與濾波器的諧振頻率ω0相等,即ω0=ωd.則濾波器在ω0諧振時(shí),電磁波在ZIM中傳播L2的相位延遲Δδ2≈0,ZIM 將左側(cè)的電磁波場(chǎng)直接傳遞到右側(cè).這時(shí)濾波器整體的諧振特性與寬度為L(zhǎng)1的線缺陷濾波器類似.當(dāng)頻率偏離ω0時(shí),諧振腔中的相位延遲偏離諧振條件使得透射響應(yīng)迅速降低.不含有ZIM 時(shí),濾波器中只有L1部分的相位延遲Δδ1隨頻率發(fā)生變化;而對(duì)含ZIM的濾波器,ZIM 部分的相位延遲Δδ2由近似為零轉(zhuǎn)為一個(gè)非零的有限值,濾波器的響應(yīng)延遲為Δδ1+Δδ2.狄拉克錐陡峭的色散曲線意味著ZIM 的等效相速度(υp=c/n=c/)在偏離ω0時(shí)迅速減小,相位延遲Δδ2則迅速增大,因而利用ZIM在類狄拉克點(diǎn)附加的電場(chǎng)響應(yīng)特性能夠提高濾波器對(duì)頻率變化的敏感度,從而使得透射峰帶寬被有效壓縮.

3 特性仿真與機(jī)理分析

選取ZIM 的晶格常數(shù)a為838 nm,圓柱半徑R為167.6 nm,圓柱體選用GaAs,相對(duì)介電常數(shù)和磁導(dǎo)率分布為12.5 和1,背景介質(zhì)為空氣,基于COMSOL 的光子能帶計(jì)算結(jié)果如圖2(a)所示,在193.4 THz(波長(zhǎng)1552 nm)出現(xiàn)三重簡(jiǎn)并的Γ點(diǎn).兩側(cè)PMs1 的晶格常數(shù)為700 nm,圓柱半徑為60 nm,材料參數(shù)與ZIM 相同,能帶計(jì)算結(jié)果表明在150—236 THz 之間出現(xiàn)第一光子禁帶.利用COMSOL軟件的波動(dòng)光學(xué)模塊分布建立圖1(a)和圖1(b)結(jié)構(gòu)的仿真模型,計(jì)算區(qū)域的左右兩端采用散射邊界條件,在左端用1 V/m 的振幅表征平面TM 模電磁波從左向右的入射,上下兩端則采用周期性邊界,各區(qū)域均采用自由三角形網(wǎng)格,單元最大尺寸小于計(jì)算頻段最小波長(zhǎng)的1/10.

圖2 (a)零折射超構(gòu)材料的光子能帶結(jié)構(gòu);(b),(c)含零折射超構(gòu)材料濾波器的傳輸特性Fig.2.(a)The band structure of a metamaterial with zero refractive index;(b),(c) transmission characteristics of filter with zero index metamaterials.

仿真結(jié)果表明,對(duì)2n列散射體(n=6)構(gòu)成的PMs1,在ZIM 厚度L2為4a,間隙L1為300 nm 時(shí),在193.4 THz 產(chǎn)生窄帶的透射峰O1,如圖2(c)實(shí)線所示.作為對(duì)比,分析了圖1(a)所示僅由PMs1構(gòu)成的線缺陷濾波器特性,分析中圓柱半徑、晶格常數(shù)、材料參數(shù)和散射體列數(shù)均保持不變.結(jié)果表明在L1=288 nm 時(shí)能夠同樣在193.4 THz 形成透射峰,如圖2(c)虛線所示.圖2(b)為圖2(c)中透射峰頻段的局部放大圖,可以看到,含ZIM 濾波器的半峰高寬為70 GHz,不含ZIM 濾波器的半峰高寬為400 GHz.在兩側(cè)PMs1 的材料參數(shù)、厚度和結(jié)構(gòu)參數(shù)不變時(shí),引入ZIM 能夠顯著壓縮濾波器的透射帶寬.

進(jìn)一步,分析了ZIM 厚度L2變化時(shí)該集成結(jié)構(gòu)濾波特性的變化.分別取L2為2a,4a,8a,12a,兩側(cè)PMs1 保持不變,調(diào)節(jié)L1使得濾波器的透射峰保持在類狄拉克頻率193.4 THz 處.透射系數(shù)仿真結(jié)果如圖3(a)所示.可以看到,隨ZIM 厚度增大,濾波器的透射半峰高寬由147 GHz 減小至37 GHz,品質(zhì)因子明顯提高,ZIM 厚度能夠有效提高濾波器的窄帶濾波特性.圖3(b)是在上述分析中,間隙L1隨ZIM 厚度的變化.在ZIM 的厚度從2a增大到12a的過程中,L1的變化小于10%.說明在類狄拉克頻率處ZIM 的厚度對(duì)濾波器諧振特性的影響很小.圖4(a)和圖4(b)分別為ZIM 厚度L2=4a和L2=12a時(shí),峰值頻率處含ZIM 濾波器中的波場(chǎng)分布仿真結(jié)果,可以很明顯地看到,透射峰頻率處濾波器中產(chǎn)生零階諧振模式,在ZIM內(nèi)部,電場(chǎng)強(qiáng)度的分布幾乎是均勻的,ZIM 出射界面相位與入射界面相同,波傳播過程中相位幾乎不發(fā)生任何變化,即相位延遲Δδ2=0.進(jìn)一步,利用等效介質(zhì)理論計(jì)算類狄拉克頻率附近ZIM 超構(gòu)材料的等效參數(shù)[34-36],結(jié)果如圖4(c)所示,超構(gòu)材料的等效介電常數(shù)和磁導(dǎo)率在透射峰值處均為零,這與圖4(a)和圖4(b)的波場(chǎng)分布相符合.正是超構(gòu)材料在類狄拉克點(diǎn)的零折射特性,使得濾波峰值頻率受ZIM 厚度的影響很小,這時(shí)濾波器與寬度為L(zhǎng)1的不含ZIM 的線缺陷濾波器類似.為保持透射峰頻率不變,L2大幅增大時(shí)需L1略微增大的現(xiàn)象則說明在周期規(guī)模比較小時(shí),介質(zhì)的折射特性與理想周期結(jié)構(gòu)的零折射特性有所偏離,導(dǎo)致了微弱的相位延遲.而隨L2的增大,有限結(jié)構(gòu)介質(zhì)更趨近于理想零折射特性,相位延遲進(jìn)一步減小,所以需要略微增大L1來保持總的相位延遲不變.

圖3 (a)傳輸譜和(b)缺陷間隙L1 隨零折射超構(gòu)材料厚度的變化Fig.3.The variations of (a) the transmission spectrum and (b) the width of the defect with the thickness of the zero index metamaterials.

圖4 (a),(b)類狄拉克點(diǎn)頻率處濾波器中的波場(chǎng)分布;(c)類狄拉克點(diǎn)附近ZIM 的等效參數(shù);(d),(e)濾波器非零折射率透射峰O2,O3 的波場(chǎng)分布Fig.4.(a),(b) The field distribution in the filter near the Dirac-like point;(c) the effective permittivity and permeability of the ZIM as a function of frequency near the Dirac-like point;(d),(e) the field distribution in the filter at transmission peak O2 and O3 with nonzero refractive index.

當(dāng)頻率偏離ω0時(shí),濾波器中相位延遲偏離諧振響應(yīng)條件,使得透射系數(shù)下降.對(duì)不含ZIM 的濾波器,只有間隙L1導(dǎo)致的相位延遲Δδ1發(fā)生變化;而對(duì)含ZIM 的濾波器,圖4(c)中狄拉克錐陡峭的色散曲線意味著ZIM 的等效相速度在偏離ω0時(shí)由趨于無窮大轉(zhuǎn)為迅速減小,相位延遲Δδ2則由近似為零轉(zhuǎn)為迅速增大,濾波器的相位延遲Δδ1+Δδ2遠(yuǎn)大于Δδ1,即利用ZIM 在類狄拉克點(diǎn)附加的強(qiáng)色散效應(yīng)能夠提高濾波器頻率響應(yīng)的敏感度,使得透射峰帶寬被有效壓縮,且Δδ2的變化隨ZIM 厚度的增大而增大,因此透射帶寬隨ZIM 厚度的增大而減小.同時(shí),在頻率偏離狄拉克點(diǎn)頻率時(shí),εeff和μeff的不同的變化趨勢(shì)會(huì)使得ZIM 的阻抗迅速變化,與空腔中空氣層間產(chǎn)生阻抗失配,會(huì)使部分光在空腔中產(chǎn)生多次反射及相位延遲,也能提高透射峰對(duì)頻率變化的敏感度.

含ZIM 濾波器中相位延遲Δδ1+Δδ2遠(yuǎn)大于Δδ1還導(dǎo)致在光子禁帶中產(chǎn)生更多的高階透射峰.對(duì)圖2(c)中含ZIM 濾波器在O2(181.2 THz)和O3(205.5 THz)透射峰處的波場(chǎng)分布分析表明,這兩個(gè)透射峰都是由濾波器的一階諧振模式導(dǎo)致的,如圖4(d)和圖4(e)所示,同樣的分析還表明更多的透射峰是由更高階的諧振模式導(dǎo)致的.通常諧振腔的一階模式出現(xiàn)在遠(yuǎn)高于零階諧振頻率的位置,ZIM 的零折射和強(qiáng)色散特性使得濾波器在零階諧振兩側(cè)都產(chǎn)生更高階的諧振模式.

性能良好的濾波器需要在實(shí)現(xiàn)超窄帶濾波的同時(shí)保持高耦合效率.對(duì)圖1(a)所示的不含ZIM的超構(gòu)材料濾波器,增大缺陷兩側(cè)超構(gòu)材料PMs1的厚度將使得濾波器趨于理想缺陷態(tài),因而能夠有效壓縮透射峰帶寬.分析表明,當(dāng)缺陷兩側(cè)PMs1的2n列散射體從n為10 增大到16 時(shí),濾波帶寬可由30 GHz 壓縮到3 GHz 以下,如圖5(a)所示.但兩側(cè)結(jié)構(gòu)厚度的增大會(huì)使得諧振腔的耦合效率降低,因而峰值透過率降低至0.3.圖5(b)為含ZIM濾波器的帶寬壓縮結(jié)果,在兩側(cè)PMs1 的散射體保持n為10 時(shí),將ZIM 的厚度從L2=8a增大到L2=30a,能夠?qū)⑼干鋷捰?2 GHz 壓縮到3 GHz 以下,同時(shí)峰值透過率能夠保持在0.9 以上.零相位延遲特性使得ZIM 在厚度大幅增大時(shí)也能與電磁場(chǎng)保持高的耦合效率,而增大ZIM 的厚度能顯著壓縮濾波帶寬,因而在超構(gòu)材料濾波器中引入ZIM 能夠?yàn)榧骖櫝瓗捄透叻逯低高^率設(shè)計(jì)提供思路.

圖5 (a)不含ZIM 濾波器的傳輸譜隨超構(gòu)材料厚度的變化;(b)含ZIM 濾波器的傳輸譜隨ZIM 厚度的變化Fig.5.The variations of the transmission spectrum (a) with the thickness of PMs1 for the filter without ZIM and (b) with the thickness of ZIM for the filter with ZIM.

對(duì)上述含有復(fù)雜結(jié)構(gòu)的光學(xué)器件,邊界輻射、材料損耗等都是影響其光學(xué)性能的重要因素.為此,進(jìn)一步將圖1 中的二維模型拓展為有限長(zhǎng)度圓柱陣列的三維模型,如圖6(a)所示.圖中結(jié)構(gòu)尺寸與圖1 相同,厚度為200 nm,上下表面設(shè)置為良導(dǎo)體(Ag).研究表明,當(dāng)上下邊界設(shè)置完為完美導(dǎo)體邊界時(shí),可以實(shí)現(xiàn)與二維模型類似的窄帶濾波效果.但將上下邊界層設(shè)置為實(shí)際良導(dǎo)體時(shí),透射峰減小,表明對(duì)含ZIM 的器件,邊界的輻射損耗對(duì)濾波器透射效果的影響更為明顯.同時(shí),考慮分別在缺陷兩側(cè)PMs1 中和ZIM 中的GaAs 介質(zhì)柱中引入介電損耗.不同損耗因子下,濾波器透射系數(shù)隨頻率的變化如圖6(b)和圖6(c)所示.可以看到,同樣的損耗因子,ZIM 介質(zhì)中材料損耗導(dǎo)致的透射峰降低要遠(yuǎn)大于兩側(cè)的PMs1.顯然,ZIM 中狄拉克點(diǎn)的電磁場(chǎng)模式為特殊的自旋模式,由于這樣的模式更利于電磁場(chǎng)能量的轉(zhuǎn)化吸收,使得濾波器對(duì)材料損耗更為敏感.因此,對(duì)含ZIM 材料的器件,其材料損耗和結(jié)構(gòu)輻射是設(shè)計(jì)時(shí)需要重視的問題.

圖6 (a)有限厚度的線缺陷+零折射超構(gòu)材料濾波器示意圖;(b)考慮PMs1 中介質(zhì)柱損耗的濾波器傳輸譜;(c)考慮ZIM 中材料損耗的濾波器傳輸譜Fig.6.(a) Schematic diagram of the finite thickness filter with line defect and zero index metamaterial;(b),(c) transmission spectrum considering different material losses of dielectric columns in defects and in ZIM.

4 結(jié)論

本文提出一種基于零折射超構(gòu)材料的高耦合超窄帶光學(xué)濾波器設(shè)計(jì)思路.研究表明,在缺陷態(tài)濾波器中嵌入具有類狄拉克錐光子能帶特性的超構(gòu)材料,通過調(diào)節(jié)尺寸和材料參數(shù)使得狄拉克點(diǎn)頻率與缺陷態(tài)的諧振頻率重合時(shí),零折射材料的零相位延遲和強(qiáng)色散特性能夠使濾波器在保持高峰值透過率的同時(shí)有效壓縮濾波帶寬,實(shí)現(xiàn)高耦合超窄帶設(shè)計(jì).該結(jié)果為光學(xué)超構(gòu)材料在波分復(fù)用等光集成系統(tǒng)的應(yīng)用提供了新的技術(shù)途徑,也為零折射率超構(gòu)材料的應(yīng)用探索提供了新的思路.