寬絕對禁帶的一維磁性光子晶體結構?

陳敏 萬婷 王征 羅朝明 劉靖

(湖南理工學院信息與通信工程學院,復雜工業物流系統智能控制與優化湖南省重點實驗室,岳陽 414006)(2016年4月22日收到;2016年10月16日收到修改稿)

寬絕對禁帶的一維磁性光子晶體結構?

陳敏 萬婷 王征 羅朝明?劉靖

(湖南理工學院信息與通信工程學院,復雜工業物流系統智能控制與優化湖南省重點實驗室,岳陽 414006)(2016年4月22日收到;2016年10月16日收到修改稿)

提出了一種具有寬絕對禁帶的一維磁性光子晶體結構,該結構由相同的折射率和物理厚度以及不同的波阻抗的兩種磁性材料交替組合而成.通過傳輸矩陣法分析可得,相比于非磁性光子晶體,該光子晶體的禁帶對入射角和偏振都不敏感,從而具有更寬的絕對禁帶.合適地調節兩種磁性材料的參數,增加兩者波阻抗的差值,該光子晶體的絕對禁帶寬度也相應地增加;調節兩種磁性材料的物理厚度,其絕對禁帶中心也會隨之調整;最后,將兩個滿足上述條件的一維磁性光子晶體組成異質結構,其第一禁帶寬度與禁帶中心之間的比值可達到1.41以上.

絕對禁帶,光子晶體,磁性材料

1引 言

近年來,光子晶體的許多優良特性和潛在應用不斷被發現,其研究已經引起越來越多的研究者的關注[1,2].光子晶體又稱“光學半導體”,它具有的最為顯著的性質是光子禁帶[3,4].光子晶體能夠阻止頻率落在光子禁帶中的電磁波的傳播,但是通常情況下光子禁帶對入射角和偏振態都較為敏感.為了實現某一頻帶的電磁波在任意偏振和角度入射下都能被阻止,可以通過合理設計光子晶體結構使得其具有絕對禁帶[5?8].這種絕對禁帶在全方位高反射鏡[9,10]、低損耗波導[11]、光開關[12]和空間濾波器等[13?15]方面都具有許多潛在的應用.

非磁性光子晶體的絕對禁帶通常是比較窄的,有必要對其禁帶進行擴展.目前擴展光子晶體禁帶的方法主要有三種,第一種是通過調整光子晶體的材料和旋轉對稱性來擴展禁帶[8,16];第二種是將幾個光子晶體組成異質結構來擴展禁帶[17?20].但是前面兩種方法對非磁性光子晶體的禁帶擴展都是在某一種偏振下進行的,同時該禁帶還對入射角較為敏感,從而絕對禁帶擴展比較小.隨著超常介質材料的出現,人們提出了第三種擴展禁帶的方法,通過將超常介質材料引入到光子晶體中來擴展禁帶[21?23].由于超常介質材料具有對入射角和偏振的不敏感性,從而能夠擴寬光子的絕對禁帶,但這種超常介質材料在自然界中是不存在的且其制備相當困難.為了解決這一不足,Ouyang等[24]將磁性材料引入光子晶體中,使得其絕對禁帶得到了一定的擴展.

本文提出了一種寬絕對禁帶的一維磁性光子晶體結構,組成該結構的兩種磁性材料具有相同的折射率和物理厚度以及不同的波阻抗.該結構禁帶對光的入射角和偏振態都具有不敏感性.通過改變該結構的波阻抗、材料的物理厚度以及構成的異質結構來研究其對光子禁帶的影響.

2結構模型與理論

圖1(a)是由兩種磁性材料構成的一維磁性光子晶體結構(AB)N示意圖,組成光子晶體的A,B兩種材料的折射率和磁導率分別為:nA,μA,nB,μB,對應的物理厚度分別為dA,dB,N表示周期數.設平面電磁波由空氣入射到光子晶體中,θ為其入射角.為了準確地描述這種結構的傳輸特性,需要確定任意波矢分量的反射和透射,該結構可由下面的2×2矩陣來表示[25,26]:

這里的Tm?1,m為第m?1層到第m層的轉換矩陣,

其中rm?1,m和tm?1,m則分別表示從第m?1層到第m層的反射系數和透射系數;Pm為第m層的傳輸矩陣

上式中的dm和kmz分別表示第m層介質的厚度和沿z方向的波矢分量.在本文的討論中,rAB,tAB分別表示圖1(b)所示的介質A到介質B界面的反射系數和透射系數;rBA,tBA則分別表示圖1(c)所示的介質B到介質A界面的反射系數和透射系數.具體表示為[25,26]

上述公式中θA和θB分別是入射電磁波在介質A和介質B中的傳輸角,表示介質A對應的波阻抗,ηB則表示介質B對應的波阻抗.整個光子晶體的透射系數和反射系數可由下面的公式來描述[25]:

圖1 (網刊彩色)光子晶體結構示意圖及兩個不同界面的傳輸圖Fig.1.(color online)Schematic of photonic crystal structure and two transmission diagrams of di ff erent interfaces.

其中字母p,s分別對應平行和垂直偏振情形.整個介質的反射率和透射率用下式來表示[25]:上式中kmz,θm和k0z,θ0分別表示出射介質m和入射介質中沿z方向的波矢分量和傳輸角.

根據(1)式,圖1(a)所給出的一維光子晶體結構可用下面的矩陣來表示:

M=T0A(PATABPBTBA)N?1PATABPBTB0.(10)為了獲得絕對禁帶,必須考慮其偏振特性,只有當禁帶對偏振不敏感時才能獲得更寬的絕對禁帶.由(1)—(9)式可知,物理厚度相同的各層對應相同的傳輸矩陣Pm,因而Pm是與偏振無關的;而轉換矩陣Tm,m?1由于rm?1,m和tm?1,m的偏振敏感性而偏振相關.在圖1結構中,很明顯影響該結構偏振特性的是介面A到界面B的反射系數rAB和透射系數tAB,以及界面B到界面A的反射系數rBA和透射系數tBA.結合(4)—(7)式分析可知,只有當θA和θB相等時(即組成光子晶體結構的兩種材料的折射率相等時),不同偏振對應的透射系數和反射系數才相等,從而使得TE和TM是偏振不敏感的,以達到擴展光子晶體絕對禁帶的目的.因此,我們可得出滿足下面條件的磁性光子晶體結構能擴展其絕對禁帶寬度,即組成光子晶體結構的兩種材料具有相同的折射率和物理厚度以及不同的波阻抗.

3數值模擬與分析

首先,我們比較研究一維磁性和非磁性光子晶體兩種結構的禁帶特性.為了研究的方便,我們參照文獻[24,27,28]忽略材料損耗來選取相關參數,理論探索一維磁性光子晶體的禁帶特性.根據第二部分得出的擴展禁帶的條件(相同的折射率和物理厚度以及不同的波阻抗),我們選取的兩種磁性材料A和B的結構參數如下:折射率nA=nB=n=3.4;磁導率μA=2.5,μB=1.5;每層介質的物理厚度為dA=dB=λ0/(4n),其中λ0為真空中禁帶的中心波長;周期數N=6.通過采用傳輸矩陣法[13?15],得到了在不同偏振情況下以不同角度入射的一維磁性光子晶體(AB)6的傳輸譜,如圖2(a)—(d)所示.從該圖中可以看出,在禁帶處TE和TM偏振的傳輸曲線幾乎重合,從而可認為該禁帶是偏振不敏感的;而且禁帶出現的位置幾乎不隨角度的改變而移動,這些都有利于絕對禁帶的擴展,該結構絕對禁帶的歸一化頻率帶寬(禁帶寬度與禁帶中心之間的比值)達到了0.41.為了與之進行對比,我們選取了具有相同波阻抗比值的兩種非磁性材料C和D(即ηA/ηB= ηC/ηD=5/3),其結構參數為nC=1.5,nD=2.5,dC= λ0/(4nC),dD= λ0/(4nD),該非磁性光子晶體 (CD)6的傳輸譜如圖2(a′)—(d′)所示.從該圖中可以看出,非磁性光子晶體的TE和TM偏振的傳輸曲線隨著入射角的增加差別越來越大,從而該禁帶是偏振相關的;并且隨著角度的增加TE偏振的禁帶左右邊緣都向右移動,但TM偏振的禁帶左邊緣向右移動而右邊緣先向右移后又向左偏移;其絕對禁帶的歸一化頻率帶寬只有0.14,相比于前面的磁性光子晶體而言,其絕對禁帶寬度明顯窄得多.

圖2 (網刊彩色)一維磁性光子晶體和非磁性光子晶體的傳輸譜Fig.2.(color online)Transmission spectrum of one dimensional magnetic photonic crystals and non-magnetic photonic crystals.

需要指出的是,通常情況下的磁性材料都是有損耗的,下面研究損耗對本文中磁性光子晶體禁帶特性的影響.眾所周知,損耗參數可用折射率的虛部來描述[29,30],在這里將磁性材料A和B的折射率調整為:nA=nB=3.4+0.1i,其他參數與圖2(a)—(d)保持一致,相應的含損耗磁性材料的一維光子晶體(AB)6的傳輸譜如圖3(a)—(d)所示.由圖中可看出其透過率相對于圖2(a)—(d)中不考慮損耗時的情形略有降低,但絕對禁帶寬度變化很小,幾乎不變.因此在本文的討論中,為了研究的方便,都忽略了磁性材料的損耗參數對禁帶特性的影響,即假定磁性材料是沒有損耗的.

根據一維磁性光子晶體絕對禁帶擴展的條件,將進一步研究兩種材料的波阻抗對其絕對禁帶寬度的影響.假定折射率和物理厚度都與前面的分析一致,即nA=nB=n=3.4,dA=dB=λ0/(4n).在保持兩種材料介電常數不變的情況下調節磁導率使得波阻抗之間的差值依次增大,圖4(a)—(c)分別是兩種磁性材料比值為ηA/ηB=2,4,6的光子晶體結構(AB)6的禁帶.從圖中可以得出,ηA/ηB=2,4,6對應的絕對禁帶的歸一化頻率范圍分別為0.78—1.25,0.6—1.44,0.50—1.53,其歸一化的頻率帶寬分別為0.47,0.84,1.03.由此可以總結出,隨著兩種材料波阻抗差值的增加,該光子晶體結構的絕對禁帶的頻率帶寬得到了顯著的擴寬.

然后,將研究探討每層材料的物理厚度對該光子晶體結構傳輸特性的影響.我們選取了圖4(b)中除每層物理厚度外相同的其他結構參數,即折射率nA=nB=n=3.4,波阻抗比值為ηA/ηB=4,磁導率μA=6.0,μB=1.5;同時將每層物理厚度分別調整為dA=dB=λ0/(8n)和dA=dB=λ0/(16n),其對應的不同角度下的傳輸譜如圖5(a)—(d) 和圖5(a′)—(d′)所示. 對比兩者可以發現,當材料的物理厚度變小時,禁帶中心依次向右移動,其歸一化頻率帶寬分別為0.78和0.80.由此我們可以總結得出,隨著材料每層的物理厚度d的縮小,絕對禁帶的禁帶中心依次右移,而對禁帶的歸一化頻率帶寬的影響很小.

圖3 (網刊彩色)含損耗磁性材料一維光子晶體的傳輸譜Fig.3.(color online)Transmission spectrum of one dimensional magnetic photonic crystals containing the lossy magnetic materials.

圖4 (網刊彩色)不同波阻抗比的一維磁性光子晶體禁帶結構Fig.4.(color online)Bandgap structures of one-dimensional magnetic photonic crystals with di ff erent wave impedances.

圖5 (網刊彩色)不同厚度的一維磁性光子晶體傳輸譜Fig.5.(color online)Transmission spectra of one-dimensional magnetic photonic crystals with di ff erent thicknesses.

圖6 (網刊彩色)三種不同結構的光子晶體的禁帶結構Fig.6.(color online)Bandgap structures of three kinds of photonic crystals with di ff erent structures.

最后研究利用該一維磁性光子晶體組成的異質結構來進一步擴展其絕對禁帶.異質結構的絕對禁帶的擴展應滿足疊加原理,即要求組合成異質結構的兩個光子晶體結構中的一個結構所產生的導帶被包含在另一個結構的絕對禁帶中.由此,我們選擇了圖5中的兩個磁性光子晶體,并重畫了它們的禁帶結構圖,即每層材料的物理厚度為dA=dB= λ0/(8n)和dA′=dB′= λ0/(16n)時的光子晶體的禁帶結構圖,具體如圖6(a)—(b)所示.從圖6(a)中可以看出,第一個光子晶體結構(AB)6的第一和第二絕對禁帶的歸一化頻率范圍依次為 1.18—2.85和5.37—6.85;圖6(b)為第二個光子晶體結構(A′B′)6的禁帶結構,其絕對禁帶的歸一化頻率范圍為2.37—5.68.很明顯,第一個光子晶體的第一和第二絕對禁帶之間的導帶被包含在第二個光子晶體的絕對禁帶中,從而有利于禁帶的擴展.圖6(c)為這兩個光子晶體組成的異質結構(AB)6(A′B′)6的禁帶結構,其絕對禁帶的歸一化頻率范圍為1.18—6.85,對應的歸一化頻率帶寬達1.41.結合圖6(a)—(c)可以總結出,滿足上述疊加條件的兩個光子晶體結構所構成的異質結構的絕對禁帶就是這兩個光子晶體絕對禁帶的疊加,從而實現了絕對禁帶的擴展.

4結 論

本文提出了一種具有寬絕對禁帶的一維磁性光子晶體結構,該結構由折射率相同和物理厚度相同但波阻抗不同的兩種磁性材料分別交替組合而成.首先,我們對比分析了磁性和非磁性一維光子晶體的禁帶特性,進一步論證了上述絕對禁帶擴展的條件.然后研究了波阻抗的改變對該絕對禁帶擴展的影響,結果發現:隨著兩種材料波阻抗差值的增大,絕對禁帶得到了很明顯的擴展.同時也得出了在保持所有波阻抗不變的情況下,通過磁性材料的物理厚度的調整,其禁帶中心隨之改變,但絕對禁帶寬度幾乎不變.最后我們用兩個一維磁性光子晶體構成異質結構,研究發現該異質結構的絕對禁帶寬度是這兩個光子晶體絕對禁帶的疊加.這種磁性光子晶體具有寬的絕對禁帶,在集成光學、光纖通信以及激光系統等領域有著潛在的應用,比如可實現具有偏振無關、全方向等良好性能的反射鏡、光開關、光學濾波器等.

[1]Yablonovitch E 1987 Phys.Rev.Lett.58 2059

[2]John S 1987 Phys.Rev.Lett.58 2486

[3]Joannopoulos J D,Meade R D,Winn J N 1995 Photonic Crystals:Molding the Flow of Light(Princeton:Princeton Univ.Press)

[4]Sakoda K 2001 Optical Properties of Photonic Crystals(Berlin:Springer-Verlag)

[5]Dowling J P 1998 Science 282 1841

[6]Liu H,Yao J Q,Li E B,Wen W Q,Zhang Q,Wang P 2006 Acta Phys.Sin.55 230(in Chinese)[劉歡,姚建銓,李恩邦,溫午麒,張強,王鵬2006物理學報55 230]

[7]Cheng X P,Cao Q X 2008 Acta Phys.Sin.57 3249(in Chinese)[程旭攀,曹全喜 2008物理學報 57 3249]

[8]Fink Y,Winn J N,Fan S,Chen C,Michel J,Joannopoulos J D,Thomas E L 1998 Science 282 1679

[9]Ibanescu M,Fink Y,Fan S,Thomas E L,Joannopoulos J D 2000 Science 289 415

[10]Jiang L,Zheng G,Shi L,Yuan J,Li X 2008 Opt.Commun.281 4882

[11]Hart S D,Maskaly G R,Temelkuran B,Prideaux P H,Joannopoulos J D,Fink Y 2002 Science 296 510

[12]Lu Y H,Huang M D,Park S Y,Kim P J,Nahm T U,Lee Y P,Rhee J Y 2007 J.Appl.Phys.101 036110

[13]Luo Z M,Tang Z,Xiang Y,Luo H,Wen S 2009 Appl.Phys.B 94 641

[14]Luo Z M,Qu S,Liu J,Tian P 2013 J.Mod.Opt.60 171

[15]Luo Z M,Chen M,Liu J,Lei D J 2016 Opt.Commun.365 120

[16]Winn J N,Fink Y,Fan S,Joannopoulos J D 1998 Opt.Lett.23 1573

[17]Zhang J,Benson T M 2013 J.Mod.Opt.60 1804

[18]Suthar B,Bhargava A 2012 Opt.Commun.285 1481

[19]Joseph S,Ha fi z A K 2014 Optik 125 2734

[20]Han P,Wang H 2005 J.Opt.Soc.Am.B 22 1571

[21]Feng X,Li H 2013 Eur.Phys.J.D 67 1

[22]Xiang Y,Dai X,Wen S,Fan D 2007 J.Opt.Soc.Am.A 24 A28

[23]Yin C P,Dong J W,Wang H Z 2009 Eur.Phys.J.B 67 221

[24]Ouyang Z B,Mao D,Liu C P,Wang J C 2008 J.Opt.Soc.Am.B 25 297

[25]Yariv A,Yeh P 2007 Optical Electronics in Modern Communications(New York:Oxford University Press)pp199–204

[26]Yeh P 1988 Optical Waves in Layered Media(New York:Wiley)pp58–67

[27]Sigalas M M,Soukoulis C M,Biswas R,Ho K M 1997 Phys.Rev.B 56 959

[28]Teng C C,Zhou W,Zhuang Y Y,Chen H M 2016 Acta Phys.Sin.65 024210(in Chinese)[滕晨晨,周雯,莊煜陽,陳鶴鳴2016物理學報65 024210]

[29]Kong J A(translated by Wu J)2003 Electromagnetic Wave Theory(Beijing:Publishing House of Electronics Industry)pp81,82(in Chinese)[孔金甌著(吳季譯)2003電磁波理論(北京:電子工業出版社)第81,82頁]

[30]Wang L G,Chen H,Zhu S Y 2005 Opt.Lett.30 2936

PACS:42.70.Qs,42.55.Tv,75.50.–yDOI:10.7498/aps.66.014204

*Project supported by the National Natural Science Foundation of China(Grant No.61205126),the Science and Technology Program of Hunan Province,China(Grant No.2016TP1021),and the Experimental Project of College Students in Hunan Province and Hunan Institute of Science and Technology,China.

?Corresponding author.E-mail:zhaomingluo@hnu.edu.cn

One-dimensional magnetic photonic crystal structures with wide absolute bandgaps?

Chen Min Wan Ting Wang Zheng Luo Zhao-Ming?Liu Jing

(Key Laboratory of Hunan Province on Intelligent Control and Optimization of Complex Industrial Logistics System,College of Information and Telecommunications Engineering,Hunan Institute of Science and Technology,Yueyang 414006,China)(Received 22 April 2016;revised manuscript received 16 October 2016)

The photonic absolute bandgaps have many potential applications in speci fi c fi elds,and some methods to enlarge the absolute bandgaps,such as adjusting the material and the rotational symmetry,constituting a heterostructure have been explored.Recently,with the occurring of metamaterial,the photonic crystal based on metamaterial has also realized the wide absolute bandgaps.However,the metamaterial is an arti fi cially structured material of which the construction is more complicated.In this paper,one-dimensional magnetic photonic crystal structure with wide absolute bandgaps is proposed,which is composed of two kinds of magnetic materials with the same refractive index and physical thickness but di ff erent wave impedances.First of all,the transmission properties of one-dimensional magnetic and non-magnetic photonic crystals with the same wave impedance ratio are studied by using transfer matrix method.It is shown that the normalized frequency bandwidth of magnetic photonic crystal,i.e.the ratio of the band of bandgap to its center,is 0.41,while the normalized frequency bandwidth of the non-magnetic photonic crystal is 0.14.From the results,we can conclude that the absolute bandgap of the above magnetic photonic crystal is wider than that of non-magnetic photonic crystal because the former bandgap is not sensitive to the incident angle nor polarization.Secondly,we adjust the wave impedance ratios of the two kinds of magnetic materials and make them respectively reach 2,4and 6,with the refractive index and the physical thickness kept unchanged.By analyzing their transmission properties,it is found that the normalized frequency bandwidths of the absolute bandgaps are respectively 0.47,0.84and 1.03,and the greater the di ff erence between the two wave impedances,the wider the normalized frequency bandwidth is.Thirdly,we investigate the in fl uence of the per-layer physical thickness of the magnetic material on the bandgap,with the other parameters remaining unchanged.It is shown that the center of the absolute bandgap shifts toward high frequency with the decrease of the per-layer physical thickness.Finally,a kind of heterostructure is constructed by the above two one-dimensional magnetic photonic crystals.The normalized frequency ranges of the fi rst and the second absolute bandgap of one magnetic photonic crystal structure are respectively 1.18–2.85and 5.37–6.85.The normalized frequency range of the absolute bandgap of the other magnetic photonic crystal is 2.37–5.68.The normalized frequency range of the absolute bandgap of the heterostructure can be enlarged to 1.18–6.85and the corresponding normalized frequency bandwidth can reach more than 1.41.The wide absolute bandgaps can be applied to integrated optics,optical fi ber communication and high-power laser systems,according to which we may design the polarization-independent and omnidirectional devices such as re fl ectors,optical switchers and optical fi lters.

absolute bandgap,photonic crystals,magnetic material

10.7498/aps.66.014204

?國家自然科學基金(批準號:61205126)、湖南省科技計劃項目(批準號:2016TP1021)和湖南省和湖南理工學院大學生實驗項目(批準號:湘教通[2016]283號,校[2016]21號)資助的課題.

?通信作者.E-mail:zhaomingluo@hnu.edu.cn