PT對稱晶格勢中的光孤子傳輸特性的研究

程家曦，倪書進，徐四六*

(1.湖北科技學院 電子與信息工程學院，湖北 咸寧 437100；2.咸寧實驗中小學，湖北 咸寧 437100 )

PT對稱晶格勢中的光孤子傳輸特性的研究

程家曦1，倪書進2，徐四六1*

(1.湖北科技學院 電子與信息工程學院，湖北 咸寧 437100；2.咸寧實驗中小學，湖北 咸寧 437100 )

利用改進的平方算符法和傅里葉配點法以及分步傅里葉算法研究了二維非線性PT對稱晶格勢中的基極和二級光孤子。研究表明，兩種孤子功率都隨著傳播常數的增加而增大，隨著調制深度的增加而減小。二極孤子內部的排斥力作用導致光束在PT勢中分裂成四極孤子。而對于PT 晶格勢，當實部的晶格調制系數V0增大時，PT勢的相對強度的帶隙結構發生變化，PT強度呈波浪式增大，而孤子在低能量區域是穩定的，高能量區域不穩定。

改進的平方算符法；分步傅里葉算法；多極光孤子；PT對稱勢

孤子又是一種特殊的超短脈沖，孤立波是一種在傳播過程中保持形狀、速度、幅度不變的脈沖狀行波[1]。鑒于孤子具有的這些特性，孤子系統在現實世界有著廣泛的應用潛力[2]。因此，孤子的形成及其傳輸特性的研究在目前是一個十分誘人的課題。

近年來, PT對稱晶格介質中的空間孤子一直是研究的熱點,人們對它們的各種獨特的性質例如相互作用、穩定性進行了大量研究[3-5]。 PT對稱即哈密頓量宇稱-時間對稱，在量子力學中，哈密頓量是要對所有物理可觀察量都必須是厄米共軛[6]。1998年，C.M.Bender提出了一個非厄米-哈密頓量，即PT對稱哈密頓量，它可以完全展現出實譜，要滿足PT對稱性的要求，則該哈密頓量要與算子共有一系列相同的本征值[7]。則勢函數,即PT對稱勢的實部必須是偶函數,虛部必須是奇函數。當超過某一閾值時，這種哈密頓量會發生相變，這時的本征值譜不再全為實數，而是變成部分為復數[8]。由于首次引入PT對稱勢的概念 ，自此，線性與非線性系統中的PT對稱哈密頓量相關的物理現象被廣泛關注。光在PT對稱哈密頓量介質中具有獨特的傳輸特性，2008年，Z.H.Musslimani等人對一維和二維PT周期晶格中的光孤子傳輸特性做出了跨越性的研究[9]，使得該領域又迎來了新的研究熱潮。最近幾年，PT對稱晶格中的非局域帶隙孤子[10]，渦旋多峰孤子[11,12]等；盡管如此，對PT對稱光空間中光孤子的傳輸特性研究仍有待深入，還有很多孤子類型與特殊的性質等著我們去發現。

本文基于二維非局域非線性薛定諤方程, 在前期研究[13-26]的基礎上，通過數值模擬得到非局域暗孤子解, 然后提出了暗孤子穩定性分析理論, 并對其數值求解得到了非局域暗孤子的穩定性分析圖, 最后利用加噪聲的傳輸驗證了穩定性分析理論的正確。

一、理論模型

光束在(2+1)維PT對稱晶格勢非線性介質中沿著z軸傳播，那么歸一化的非線性薛定諤方程

(1)

其中，φ為包絡函數，T是調制深度，V(x,y),W(x,y)分別為PT對稱晶格勢的實部和虛部,其非線性項為三五次競爭。

假設方程(1)中孤子解的形式為φ=f(x,y)eiμz，μ為傳播常數，模f(x,y)滿足該方程

(2)

求解孤子，我們利用平方算符法對(2)進行數值求解，而對于所得到的孤子穩定性，我們對方程(1)采用微擾解進行穩定性分析

φ=eiμz{f(x,y)+[g(x,y)-h(x,y)]eλz+

[g(x,y)+h(x,y)]*eλ*z}

(3)

將上式代入到式(1)并對g(x,y),h(x,y)進行線性化，得到他們的本征方程為：

(4)

其中，

(5)

我們利用傅里葉配點法[22]對孤子進行穩定性求解，若本征值的實部為零，則孤子線性穩定，反之，線性不穩定。孤子的總功率P=∫-∞+∞∫-∞+∞|f(x,y)|2dxdy。

求解非線性方程的解析解的條件是滿足可積，但是在一般情況下可積性條件是不滿足的，為了探究方程(1)的孤子解，我們將使用改進的平方算符法求得其數值解，即近似解。在此，選用高斯光為入射光，取PT晶格為雙曲-三角對稱晶格勢

V(x,y)=V0(sech2(x)+sech2(y))(cos(-x2)+cos(-y2))

(6a)

W(x,y)=W0(tanh(2x)+tanh(2y))(cos(-x2)+cos(-y2))

(6b)

為了便于討論，取PT勢的實部強度V0=1。虛部強度W0=0.1。PT勢的分布如圖1所示，可以看到其實部V 關于原點偶對稱，虛部W關于原點奇對稱。PT勢強度隨著調制參數V0、W0的增大而波浪式增大。

圖1 PT勢分布(a)(b)(c)為PT勢實部分布；(d)(e)(f)為PT勢虛部分布；(g)為PT勢實部調制參數V0與PT勢的關系分布；(h)為PT勢虛部調制參數W0與PT勢的關系分布；(i)為PT晶格勢函數系數對PT勢強度影響的關系圖；

二、基極孤子

首先對方程(1)進行求解得到基極孤子，圖2(a)為μ與P的關系圖，可以看出功率P隨著μ的增加而增加，圖中實線表示穩定區域,虛線表示不穩定區域。圖2(b)是T與功率P的關系，隨著T的增加P減小，當μ增加時P隨著T減小的速度變緩。圖2(c)是μ與微擾實部Re(λ)之間的關系，當實部>0時，對應的基極孤子不穩定。孤子在光晶格傳輸有一個存在區域與穩定區域，調制深度與傳播常數μ的存在閾值μ1穩定閾值μ2是一一對應的，基極孤子存在時傳播常數μ>μ1(圖2(a))，而在穩定時μ>μ2(圖2(c))，所以基極孤子的穩定區域是δ(μ)=μ2-μ1，圖2(d)，2(e)，2(f)分別為T與μ1，μ2，δ(u)之間的關系?？梢钥闯龌鶚O孤子在調制深度T≥0時存在(圖2(d))，且隨著T的增大，μ1和μ2均增加，而穩定區域δ(μ)先增大后減小，在T=7時穩定區域最小。

圖2 (a)傳播常數μ與功率P的關系；(b) 調制深度T與P的關系；(c)傳播常數 μ與微擾增長率實部Re(λ)的關系；(d)調制深度T與傳播常數存在閾值μ1的關系；(e)調制深度T與傳播常數穩定閾值μ2的關系；(f)調制深度T與穩定區域δ(μ)的關系。

圖3(a)、(b)、(c)表示不同調制深度T和傳播常數μ時基極孤子的強度分布，圖3(a)、(b)、(c)分別對應圖2(a)功率曲線中所標出的三個點。當T=2，μ增加時，圖3(a)中孤子的光強小于圖3(b)；當μ=10，T增加時，圖3(b)中孤子的光強大于3(c)。圖3(d)、(e)、(f)是對應孤子實部的光場分布，圖3(g)、(h)、(i)為孤子的虛部分布?？梢姽鈭鰧嵅筷P于原點偶對稱，而虛部關于原點奇對稱，基極孤子與PT勢的對稱相關。圖3(j)、(k)、(l)是對應孤子的穩定性光譜，由于孤子圖3(d)的光譜僅僅包含虛部本征值，因此這個是孤子穩定的；而孤子圖3(a)、(b)的光譜包含兩對復數和兩對實數本征值，所以此孤子不穩定。對應圖2(a)孤子的穩定區域。

三、二極孤子

我們再對方程(1)求解得到二極孤子。圖4(a)是功率P與傳播常數μ之間的關系，可以看到P隨μ的增大而增大，圖中實線是表示穩定區域，虛線表示不穩定區域，其中，調制深度T增大，功率P隨μ的增加而增加的速率變緩。而圖4(b)為功率P與調制深度T之間的關系，隨著T的增加P減小，其中，傳播常數μ增大，P隨T的增加而減小的速率變緩。圖4(c)是傳播常數μ與微擾因子實部Re(λ)之間的關系，當實部>0時，對應的孤子不穩定。孤子傳輸的存在區域與穩定區域中，調制深度與傳播常數μ的存在閾值μ1穩定閾值μ2是一一對應，孤子存在時傳播常數 (圖4(a))，而孤子穩定時μ<μ2(圖4(c))，所以二極孤子的穩定區域是δ(μ)=μ2-μ1，圖4(d)-(f)分別為T與μ1，μ2，和δ(μ)之間的關系?？梢钥闯龆O孤子在調制深度T≥0時存在(圖4(d))，且隨著T的增大，μ1和μ2均增加，而穩定區域δ(μ)先增大后減小，在T=4穩定區域最小。

圖3 (a)(b)(c) 基極孤子強度分布圖，參數為(a)T = 2，μ =5，(b) T = 2 ，μ = 10，(c) T = 3，μ = 10，分別對應圖3(a)功率曲線所標出的3個點；(d)(e)(f) 基極孤子光場實部分布；(g)(h)(i) 基極孤子光場虛部分布；(j)(k)(l) 基極孤子線性穩定性光譜

圖4 (a) 傳播常數μ與功率P的關系；(b) 調制深度T與P的關系；(c)傳播常數 μ與微擾增長率實部Re(λ)的關系；(d)調制深度T與傳播常數存在閾值μ1的關系；(e)調制深度T與傳播常數穩定閾值μ2的關系；(f)調制深度T與穩定區域δ(μ)的關系。

圖5(a)-(c)表示不同調制深度T和傳播常數μ時二極孤子的光場分布，其分別對應圖4(a)中所標出的三個點。當T=6，傳播常數μ增加時，圖5(b)中孤子的光強大于5(a)；當μ=13.8時，圖5(c)中孤子的光強小于圖5(b)。圖5(d)-(f)時對應的二極孤子的光場實部，可以看出它是關于原點奇對稱。圖5(g)-(i)是二極孤子虛部分布。它是關于原點偶對稱。圖5(j)(k)(l)是對應孤子的線性穩定性光譜，由于孤子5(a)，5(b)的本征譜僅僅包含虛數特征值，因此這兩個孤子穩定；而孤子圖3(c)的本征譜包含兩對復數和兩對實特征值，因此該孤子線性不穩定。這是由于方程1中的非線性項是三五次競爭，它在與晶格的相互作用下限制了色散，非線性效應加強光束寬度變小抑制了孤子的能量轉移，使得在這種平衡狀態下得到穩定孤子。二極孤子之間存在著相互排斥的作用，這是因為光束在傳輸過程中基極孤子之間出現了π個相位差，所以在得到穩定的二極孤子時需要對基極孤子更大的調制深度，進而分裂得到二極孤子。

圖5 (a)(b)(c)二極孤子強度分布圖，參數為(a)T = 6，μ = 11.2，(b) T = 6，μ = 13.8，(c) T = 7，μ = 13.8，分別對應圖2(a)功率曲線所標出的3個點；(d)(e)(f) 二極孤子光場實部分布；(g)(h)(i) 二極孤子光場虛部分布；(j)(k)(l) 二極孤子線性穩定性光譜

四、結語

我們研究了基極和二級光孤子在PT對稱晶格勢中的傳輸特性。他們具有共同的特點：功率都隨著傳播常數的增加而增大，隨著調制深度的增加而減小。二極孤子內部的排斥力作用導致光束在PT勢中分裂成四極孤子。而對于PT 晶格勢，當實部的晶格調制系數V0增大時，PT勢的相對強度的帶隙結構發生變化，PT強度呈波浪式增大，而孤子在低能量區域是穩定的，高能量區域不穩定。

[1] 黃景寧，徐濟仲，熊吟濤．孤子：概念、原理和應用[M]．北京：高等教育出版社，2004.

[2] 何影記．空間孤子傳輸的研究[D]．激光與光譜學研究所光電材料與技術國家重點實驗室,2007.

[6]R.El-Ganainy,K.G.Makris,D.N.Christodoulides,andZ.H.Musslimani.TheoryofcoupledopticalPT-symmetricstructures[J].Opt.Lett.2007,(32): 2632～2634 .

[7]C.M.BenderandS.Boettcher,PT-symmetricquantummechanics[J].Phys.Rev.Lett.,1998,(80):5243.

[8]A.Guo,G.J.Salamo,D.Duchesne,etal.ObservationofPT-symmetrybreakingincomplexopticalpotentials[J].Phys.Rev.Lett.2009,(103):093902.

[9]Z.H.Musslimani,K.G.Makris,R.El-Ganainy,etal.OpticalsolitonsinPTperiodicpotentials[J].Phys.Rev.Lett.,2008:100, 030402 .

[10] 周博臻，花春波，徐四六，等.PT對稱晶格勢中渦旋光孤子[J]. 中國激光,2015，(5):165～171.

[11]D.MIHALACHE,D.MAZILU.RomanianReportsinPhysics[J].2008,60,3,749～761.

[12] 周博臻,徐四六,程正則. 光彈在庫墨-高斯晶格中傳輸特性的研究[J]. 物理學,ActaPhys.Sin.2013, 62(8):084210.

[13]SiliuXu,JiangchuLiangandZhongmingLi.Self-similarspatiotemporalsolitarywavesinBesselHermiteopticallattices[J].Chin.Phys.B,2011,(20):11～17.

[14] 徐四六,劉會平, 易林, 強非局域非線性介質中的二維庫墨-高斯孤子簇[J].物理學報, 2010,59(2):1069～1075.

[15]SiliuXu,JiangchuLiangandLinYi.Self-similarsolitarywavesinBesselopticallattices[J].J.Opt.Soc.Am.B,2010, 27(1):99～106.

[16]SiliuXu,JiangchuLiangandLinYi.ExactsolitonsolutionstoageneralizednonlinearSchrodingerequation[J].Commun.Theo.Phys.,2010,531, 59～66.

[22]JiankeYang.NonlinearwavesinIntegrableandNon-integrableSystens[M].PA:SocietyforIndustrialandAppliedNathematicsphiladelpkia,2010.

[24]S.L.Xu,G.P.Zhou,N.Z.Petrovi,etal.NonautonomousvectormatterwavesintwocomponentBose-Einsteincondensateswithcombinedtime-dependentharmonic-latticepotential[J].J.Opt，17,2015,105605～105611.

[26]SiliuXu,JiangchuLiangandLinYi.SpatialSolitonSolutionsinaHarmonicPotential[J].CHIN.PHYS.LETT.,2009,(26): 074215～0742219.

2095-4654(2016)12-0083-06

2016-10-08 基金項目：國家大學生創新項目資助“可控光子晶格中光孤子傳輸特性及其應用研究”(201410927007) 通訊作者：徐四六(1969-)，男，湖北省通山縣人，博士，湖北科技學院電子與信息工程學院教授，E-mail：1158447603@qq.com。

TN929.11

A