孔猛

（北京師范大學(xué)珠海分校工程技術(shù)學(xué)院北京519087）

含左手材料與Kerr 非線性不定介質(zhì)缺陷層一維光子晶體雙穩(wěn)態(tài)位移的動(dòng)態(tài)調(diào)制

孔猛

（北京師范大學(xué)珠海分校工程技術(shù)學(xué)院北京519087）

根據(jù)傳輸矩陣法和穩(wěn)態(tài)相位理論，在五個(gè)不同的頻率范圍內(nèi)研究光束在含左手材料與Kerr非線性不定介質(zhì)缺陷層一維光子晶體中的非線性傳播特性。結(jié)果表明：當(dāng)非線性不定介質(zhì)缺陷層為cut-off和never cut-off介質(zhì)時(shí)，透射光強(qiáng)隨入射光強(qiáng)的變化呈S-型曲線關(guān)系，也就是說(shuō)透射光強(qiáng)與入射光強(qiáng)存在光學(xué)雙穩(wěn)態(tài)特性。同時(shí)，透射光束的橫向位移隨入射強(qiáng)度也呈現(xiàn)出回滯曲線響應(yīng)關(guān)系，并且通過(guò)改變?nèi)肷涔馐念l率，可以實(shí)現(xiàn)一維光子晶體橫向透射位移的動(dòng)態(tài)調(diào)制。

一維光子晶體非線性不定介質(zhì)雙穩(wěn)態(tài)位移

一堯引言

隨著人造媒介的發(fā)展，特別是超常材料（metamaterials）[1]的出現(xiàn)大大的拓展了電磁材料的特性研究。典型的超常材料如介電常數(shù)ε和磁導(dǎo)率μ同時(shí)為負(fù)值的負(fù)折射率材料，因?yàn)殡姶挪ㄔ谶@種材料中傳輸時(shí)其波矢→，電場(chǎng)→和磁場(chǎng)→形成了左手關(guān)系[5]，故又稱為左手材料（left-handed materials）。電磁波在左手材料中傳輸時(shí)會(huì)出現(xiàn)許多反常物理現(xiàn)象，如負(fù)折射率、逆Doppler效應(yīng)、逆Cherenkov效應(yīng)和完美透鏡等[2]。值得關(guān)注的是，D. R. Smith等[3]首次發(fā)現(xiàn)各向異性的特異材料——不定介質(zhì)（indefinite metamaterials），不定介質(zhì)這一新型材料立即成為物理學(xué)，微電子學(xué)以及信息技術(shù)等領(lǐng)域的研究熱點(diǎn)。

近年來(lái)，隨著對(duì)GH（Goos-H?nchen）位移[4]研究的深入，非線性GH位移亦引起了廣泛關(guān)注，如O. Emile等[5]在實(shí)驗(yàn)上探測(cè)到自聚焦Kerr非線性介質(zhì)界面處的非線性GH位移，Y. S. Kivshar等[6]研究了空間光孤子在向列液晶表面處的非線性GH位移特性，P. Hou等[7]研究了含通常Kerr非線性缺陷一維光子晶體雙穩(wěn)態(tài)位移的特性等。本文研究光束在含左手材料與Kerr非線性不定介質(zhì)缺陷層一維光子晶體中的非線性傳播特性，分析入射頻率對(duì)雙穩(wěn)態(tài)位移的影響。

二、物理模型與推導(dǎo)

考慮一個(gè)對(duì)稱的多層膜結(jié)構(gòu)(AB)mANA( BA)m作為一維有限維光子晶體結(jié)構(gòu)，晶體兩側(cè)的對(duì)稱介質(zhì)C(εC,μC)為右手線性材料。A(εA,μA)為右手線性材料，B(εB,μB)為左手線性材料，N(?N,?N)為非線性不定介質(zhì)，m為周期數(shù)。介質(zhì)層A，B和N的幾何厚度分別為dA，dB和dN。由于非線性各向異性超常材料的介電常數(shù)ε?N和磁導(dǎo)率μ?N是二階張量，為使討論的問(wèn)題簡(jiǎn)化，D. R. Smith等[8]假定材料的結(jié)構(gòu)參數(shù)ε?N和μ?N可以同時(shí)對(duì)角化。在這種情況下，介電張量和磁導(dǎo)率張量可分別表示為：

其對(duì)角分量εxx可以表示為，ε和χ分別為該分量的線性介電

x3常數(shù)和Kerr非線性系數(shù)。為簡(jiǎn)單起見，我們忽略非線性材料的吸收特性。εx，μyy和μzz的色散關(guān)系可通過(guò)Drude模型來(lái)描述，即：（2）

其中Ω=ω/ ωp為歸一化的入射光頻率，Ωry=ωry/ωp和Ωrz=ωrz/ωp分別為y方向和z方向上的歸一化磁共振頻率，ωp為電等離子頻率，F(xiàn)y和Fz分別表示開口環(huán)諧振腔在y方向和z方向所占區(qū)域的比例系數(shù)。

假設(shè)一角頻率為ω的平面電磁波從左側(cè)各向同性均勻右手材料中以入射角θ0入射到光子晶體中。我們假定入射電磁波的波矢位于y-z平面，入射電磁波的電場(chǎng)可以表示為，其中β和k分別是入射電磁波Cz波矢在分界面y方向和z方向的分量，。這里我們只考慮TE偏振波，對(duì)于TM偏振波的處理方法類似。

根據(jù)Maxwell方程組及邊界條件，相鄰任意兩層的電場(chǎng)分量和磁場(chǎng)分量可以通過(guò)傳輸矩陣相連接[9] 。對(duì)于Kerr非線性不定介質(zhì)層，傳輸矩陣可表示為：（3）

式中，kz+和kz-分別為前向波和后向波波矢在z方向上的分量，且他們均依賴于前向波和后向波的電場(chǎng)強(qiáng)度U±，即：（4）

+-Nz不考慮非線性效應(yīng)時(shí)缺陷層中電磁波波矢在z方向上的分量，滿足以下形式：（5）

這里要注意的一點(diǎn)是kNz的符號(hào)必須與μyy保持一致[10]。

對(duì)于一給定參數(shù)的透射強(qiáng)度，我們可采用定點(diǎn)迭代法[11]求解一組關(guān)于U±的非線性耦合方程，便可以得到Kerr非線性不定介質(zhì)層傳輸矩陣MN的具體形式。利用傳輸矩陣法，可以得到透射系數(shù)：（6）

三、數(shù)值計(jì)算與分析

由式（2）可知，通過(guò)調(diào)整相應(yīng)的結(jié)構(gòu)參數(shù)與入射光頻率，可以在不考慮非線性效應(yīng)的情況下實(shí)現(xiàn)對(duì)不定介質(zhì)張量元εx，μyy和μzz正負(fù)符號(hào)的轉(zhuǎn)換。從式（5）可看出，正是由于張量元εx，μyy和μzz不同符號(hào)的組合，使得這種材料的色散關(guān)系或者為橢圓關(guān)系，或者為雙曲線關(guān)系。因此不同類型的Kerr非線性不定介質(zhì)將會(huì)對(duì)一維光子晶體的電磁傳播特性產(chǎn)生不同的影響。在計(jì)算中，各物理參數(shù)選擇如下：對(duì)于Kerr非線性不定介質(zhì)N[13]，F(xiàn)y=0.56，F(xiàn)z=0.78，ωp/2π=18.4GHz ，ωry/2π=ωrz/2π=15.3GHz，χ3=0.01，對(duì)于線性材料A，B和C[14]，εA=5.29，εB=-2.46，εC=6.25，μA=μC=1.00，μB=-1.00，m=3。Kerr非線性不定介質(zhì)的張量元εx，μyy和μzz均與入射頻率有關(guān)，在不同的入射頻率范圍內(nèi)這些物理量有著不同的符號(hào)，為了更明確地研究光束在含左手材料與Kerr非線性不定介質(zhì)缺陷層一維光子晶體中的電磁傳播特性，接下來(lái)我們將在五個(gè)不同的入射頻率范圍分別進(jìn)行討論。

（i）0＜Ω＜0.8315，此時(shí)εx＜0，μyy＞0，μzz＞0，即always cut-off 介質(zhì)的一種情況。在這種情況下，kN2z總是負(fù)值，任何電磁波都不會(huì)透過(guò)這種介質(zhì)板，而是被全部反射回來(lái)。

（ii）0.8315＜Ω＜1，這里εx＜0，μyy＜0，μzz＜0，此時(shí)當(dāng)入射角θ0滿足θ0＜θc=arcsin時(shí)，由式（5）可知，波矢kNz為實(shí)數(shù)。那么，電磁波便可以透過(guò)含此類型的非線性不定介質(zhì)的一維光子晶體。這種材料是cut-off 介質(zhì)的一種情況，在不考慮Kerr非線性效應(yīng)時(shí)，具有負(fù)折射效應(yīng)。

（iii）1＜Ω＜1.2536，這里εx＞0，μyy＜0，μzz＜0，即always cut-off 介質(zhì)的另一種情況，電磁傳播特性與（i）類似。

（iv）1.2536＜Ω＜1.7728，這時(shí)εx＞0，μyy＞0，μzz＜0，在這種情況下，kN2

z總是正值，任何電磁波都可以在含這種never cut-off非線性不定介質(zhì)的一維光子晶體中傳播。因此，當(dāng)光束以任何角度入射到這樣的一維光子晶體時(shí)，均可以產(chǎn)生透射位移。

（v）1.7728＜Ω，此時(shí)εx＞0，μyy＞0，μzz＞0，即cut-off 介質(zhì)的另一種情況。當(dāng)入射角θ0滿足條件θ0＜θc=arcsin時(shí)，由式（5）可知，電磁波同樣也可以透過(guò)含此類型的非線性不定介質(zhì)的一維光子晶體。

下面，我們針對(duì)（ii）,（iv）,（v）三種情況下的非線性傳播特性分別進(jìn)行詳細(xì)討論與分析。

3.1 0 .8315＜Ω＜1（εx＜0，μyy＜0，μzz＜0）

當(dāng)Kerr非線性結(jié)構(gòu)作為缺陷嵌入到一維光子晶體的中間時(shí)，禁帶中相應(yīng)出現(xiàn)缺陷模頻率。由于缺陷層中強(qiáng)烈的非線性效應(yīng)，導(dǎo)致缺陷模頻率隨入射光強(qiáng)度

的改變發(fā)生移動(dòng)，當(dāng)缺陷模頻率與入射光頻率接近時(shí)，產(chǎn)生雙穩(wěn)態(tài)現(xiàn)象[15]。數(shù)值計(jì)算中，我們采用的參數(shù)為：θ0=20°，ω/2π=16.07GHz ，

根據(jù)穩(wěn)態(tài)相位法分析，曲線的斜率表明透射光束的雙穩(wěn)態(tài)GH位移為負(fù)值。在一個(gè)非常狹窄的β區(qū)間內(nèi)（β/k0=0.8551），當(dāng)入射光強(qiáng)度在[1.31500,1.40600]區(qū)間內(nèi)增大時(shí)，曲線斜率突然變陡，表示透射光束的GH位移S在高閾值點(diǎn)2處發(fā)生躍變，下跳到一個(gè)較大值；當(dāng)入射光強(qiáng)度在區(qū)間[0.04092,0.02669]內(nèi)減小時(shí)，曲線斜率發(fā)生躍變，即S在低閾值點(diǎn)1處達(dá)到負(fù)最大，然后上跳到一個(gè)較小值。

此時(shí)一頻率為ω/2π=25.17GHz的平面電磁波從左側(cè)各向同性均勻右手材料中以入射角θ0=30°入射到含這種never cut-off非線性不定介質(zhì)的一維光子晶體中，在這種情況下，可以得到負(fù)的雙穩(wěn)態(tài)GH位移。此時(shí)所得的結(jié)果與（ii）類似，但雙穩(wěn)態(tài)GH位移的峰值略小，僅為10.65個(gè)波長(zhǎng)。透射光束相位φt在高閾值與低閾值處發(fā)生跳變。

當(dāng)入射光強(qiáng)度從一個(gè)大于高閾值點(diǎn)處的值開始減小時(shí)，S隨之增大，當(dāng)入射光強(qiáng)度減小到低閾值點(diǎn)時(shí)，S在此處達(dá)到正最大值。透射光束相位tφ的斜率在高閾值與低閾值處發(fā)生躍變，這恰恰說(shuō)明了透射光束雙穩(wěn)態(tài)GH位移的合理性。

四、結(jié)論

本文利用傳輸矩陣法和穩(wěn)態(tài)相位理論，在五個(gè)不同的頻率范圍內(nèi)研究了光束在含左手材料與Kerr非線性不定介質(zhì)缺陷層一維光子晶體中的非線性傳播特性。研究表明當(dāng)非線性不定介質(zhì)缺陷層為cut-off和never cut-off介質(zhì)時(shí)，透射光強(qiáng)與入射光強(qiáng)存在光學(xué)雙穩(wěn)態(tài)特性。。研究成果為制作高性能的光電探測(cè)及光開關(guān)器件提供了有益的參考。

[1]R.A.Shelby,D.R.Smith,and S.Schultz,Science 292,77(2001).

[2]V.G.Veselago,Sov.Phys.Usp.10,509(1968).

[3]D.R.Smith,W.J.Padilla,D.C.Vier,S.C.Nemat-Nasser,S. Schultz,Phys.Rev.Lett.84,4184(2000).

[4]F.Goos,H.H?nchen,Ann.Phys.1,333(1947).

[5]O.Emile,T.Galstyan,A.Le Floch,F.Bretenaker,Phys.Rev.Lett. 75,1511(1995).

[6]Y.S.Kivshar,Nature Physics 2,729(2006).

[7]P.Hou,Y.Y.Chen,X.Chen,J.L.Shi,Q.Wang,Phys.Rev.A 75,045802(2007).

[8]D.R.Smith,D.Schurig,Phys.Rev.Lett.90,077405(2003).

[9]N.H.Liu,S.Y.Zhu,H.Chen,X.Wu,Phys.Rev.E 65,046607 (2002).

[10]X.Chen,C.F.Li,Phys.Rev.E 69,066617(2004).

[11]S.Dutta Gupta,D.S.Ray,Phys.Rev.B 38,3628(1988).

[12]T.E.Hartman,J.Appl.Phys.33,3427(1962).

[13]R.Ruppin,Phys.Lett.A 277,61(2000).

[14]S.D.Gupta and G.S.Agarwal,J.Opt.Soc.Am.B 4,691(1987).

[15]R.Wang,J.Dong,D.Y.Xing,Phys.Rev.E 55,6301(1997).