在非線性介質(zhì)中的光傳播行為

1、在非線性介質(zhì)中的光傳播行為文石明豐一、簡介1834年時，一位蘇格蘭科學(xué)家John S. Russell1沿著一條窄而淺的運河旁騎馬，他發(fā)現(xiàn)了如下景象，當(dāng)一艘小船突然停下後，原本在船艏前方被推動的水因慣性的關(guān)係繼續(xù)向前進，奇特的是，在這單一“突起的水”的前方和後方，運河的水面是非常平靜而沒有任何波動，他好奇地騎馬跟隨這單一“突起的水”走了好幾哩後，發(fā)現(xiàn)這“突起的水”的形狀、大小仍不見有任何改變，Russell 於是在他的筆記本記下了他觀察到這樣一個“孤立升起”(solitary elevation)的水波現(xiàn)象。50年後，兩位荷蘭科學(xué)家Korteweg 和de Vries 發(fā)現(xiàn)要觀察上述“孤立升起

2、”的水波現(xiàn)象，這升起水面的振幅必須非常的大，如果這升起水面的振幅不夠大時，就只能引起一般的水波，而且在傳播不久後就會消散。他們了解到，這表示著大的振幅和小的振幅會使水波有不同的行為，也可以說水-這個傳播水波的介質(zhì)-對水波振幅的反應(yīng)是非線性的?，F(xiàn)在回到我們的主題-光在非線性材料中傳播的行為，雖然我們在大學(xué)時學(xué)電磁學(xué)及電磁波的時候，教科書總會在例題或習(xí)題中假設(shè)，電磁波傳播的介質(zhì)是線性的，意即不管電磁波的振輻大小，折射率(介電係數(shù)和介磁係數(shù)的函數(shù))總是不會改變，這樣在解題時，線性的微分方程式是比較容易解的，然而介電係數(shù)代表的是電場(E)和極化強度(P)之間的比例關(guān)係，如同彈簧的虎克定律一樣，當(dāng)外力大

3、到某一個程度時，彈簧的形變量和外力就不再是成正比，相同地，當(dāng)電場強度大到某一個程度時，極化強度和電場就不再是正比關(guān)係，如此介電係數(shù)在電磁波的電場振幅較大時，和電場振幅較小時的值必然不同，這就表示這個傳播電磁波的介質(zhì)對電磁波振幅的反應(yīng)是非線性的，換句話說，折射率將會是電磁波振幅的一個函數(shù)。在光學(xué)上，我們也可以說，折射率不是一個定值而是一個光強度的函數(shù)。但如此一來，要解一個電磁波或光波在介質(zhì)中的行為就變得複雜的多了，而且在不同材料中，折射率和電場關(guān)係也不儘相同，所以我們無法找到一個電磁波或光波在非線性材料中傳播的通解，甚至在很多的情況下，找不到解析解。和折射率受到光強影響而改變的物理現(xiàn)象有許多種，

4、如相位共軛光 (如圖一)，當(dāng)Signal 和Pump 1入射到非線性材料時，因為干涉的原因而產(chǎn)生週期性的亮暗條紋，並因此在非線性材料中產(chǎn)生週期性的折射率的變化，這週期性折射率變化就如同光柵一般，另一道光Pump 2入射後受到這個光柵的影響會以如圖1中PC光的方向繞射出去，相同的Signal 和Pump 2所產(chǎn)生的光柵也會使Pump 1以PC的方向繞射出去，這道以PC方向繞射出去的光其相位和Signal光在空間中的相位會圖一：相位共軛光的形成示意圖。互為共軛複數(shù)，故我們稱PC為Signal的相位共軛光(Phase Conjugation)。此外，二光混合的作用原理也類似，Signal 和 Pum

5、p 1 因干涉所產(chǎn)生的光柵使得Pump 1繞射到Signal 的路徑而加強Signal的強度，(反過來亦可能發(fā)生，端視所產(chǎn)生的光柵和Signal 與 Pump 1所生干涉條紋之間的位置關(guān)係)。二、自聚焦效應(yīng)、空間光孤子和光調(diào)變不穩(wěn)定性因筆者主要的研究方向是非線性介質(zhì)中光傳播的行為，所以以下我會在這方面作較多的介紹。如前所述，要詳細了解光波在非線性材料中的傳播行為的微分方程式是一個非線性方程式，並不容易去解，目前有許多專門的方法(如Inverse Scattering Method)6是設(shè)計來解決這樣的問題，但在這裏我們只從原理上定性地來看這樣的現(xiàn)象，並且只須知道兩個原則：(一)海更斯原理告訴我

6、們，波的傳播是垂直於波前的；(二)光波前進的速度是反比於折射率的；就可以討論光波在非線性介質(zhì)中傳播的行為。圖二：光的繞射，自聚焦，和光孤子等現(xiàn)象，其光強度分布和相位波前的關(guān)係示意圖。關(guān)於這方面的研究大約起始於1960年代中期，其中最早開始的研究是R. Chiao2所建議的一道光束應(yīng)該可以完全避免繞射，而在介質(zhì)中形成我們現(xiàn)在所稱的“空間光孤子”的現(xiàn)象(optical spatial soliton)。光束在自由空間傳播時，因為其波的特性，會自然的繞射(diffraction)，若由海更斯原理來解釋，我們可以知道，這道光束的波前一定是中間部分凸出，而造成光束的能量由中間部分向旁邊擴散開來。若要這道

7、光束不產(chǎn)生繞射的話，光的波前必須作適當(dāng)?shù)奶幚?，因此這項“光孤子”現(xiàn)象的想法(如圖二)能夠成立的前提是光傳播介質(zhì)對光強的反應(yīng)必須是非線性的且光越強的地方折射率越大(稱為自聚焦self-focusing材料，反之，光越強的地方折射率越小為自散焦self-defocusing材料)；因為在光束中心光較強的部分可以導(dǎo)致較大的折射係數(shù)，折射率較大的部分光波前進的速度較慢，而改變了光波的波前，如果有“適當(dāng)”的光強分佈，波前是可以被調(diào)變成一平面，因為海更斯原理，波的傳播是垂直於波前的，如此這光孤子的光束在傳播時，就可完全的消除繞射，但是當(dāng)光強分佈和光束大小不適當(dāng)時，這光束就會繞射或自聚焦，而空間光孤子的形成

8、其實就是依靠自聚焦和繞射的完全平衡，這種光孤子又稱為亮孤子(bright soliton)。另外還有一種稱為暗孤子(dark soliton)的現(xiàn)象，和剛才所敘述的亮孤子完全相反，其所存在的非線性材料必須為自散焦類型的，而其光強分布恰和亮孤子相反，中間暗，旁邊亮，其傳播行為比亮孤子更為複雜，我們在此就不多作介紹。再回到亮孤子的介紹，在1965年時，Kelley卻提出二維的亮孤子(像是一根針而不繞射的光束)是不穩(wěn)定的3，因為在當(dāng)時，唯一已知可以產(chǎn)生自聚焦現(xiàn)象的介質(zhì)都是Kerr介質(zhì)，其折射係數(shù)的變化正比於光的強度，當(dāng)自聚焦現(xiàn)象發(fā)生時，光的強度會因聚焦而增加，並進一步的提高折射係數(shù)，產(chǎn)生更強的自聚焦

9、現(xiàn)象，然而繞射的平衡作用卻不能以同一速率來增加，而使得光束不斷的縮小，直到光束的強度過高把介質(zhì)破壞為止，稱為“災(zāi)難性自聚焦”(catastrophic self-focusing)，也因如此，在此問題解決之前，並沒有任何實驗可觀察到穩(wěn)定的二維空間光孤子。附帶一提的是，在此之後，於1970年代，Hasegawa4等提出在光纖中的時間光孤子，並且由Mollenauer5等首先觀察到，時間光孤子在80年代遂成為光孤子的研究主流，許多有關(guān)時間光孤子的研究已經(jīng)成熟並且已進入實用階段。圖三：光調(diào)變不穩(wěn)定性圖案生成的示意圖。另一個和空間光孤子類似的光傳播行為就是光調(diào)變不穩(wěn)定性(modulation inst

10、ability)6，它所須要的條件和亮孤子的條件一樣，就是介質(zhì)必須是自聚焦類型的，不同的是，光孤子在入射到非線性光學(xué)材料之前，我們必須將其光束大小依其光的強度作一個調(diào)整，使得自聚焦和繞射之間可以完全的平衡，而光調(diào)變不穩(wěn)定性則不須要，光調(diào)變不穩(wěn)定性可以從平面波開始。以下我們來看看，光調(diào)變不穩(wěn)定性如何生成，如圖三，假設(shè)有一同調(diào)的平面光(coherent plane wave)入射到非線性自聚焦材料，因任何的光無論如何都會帶有雜訊(這雜訊可能來自光源的自發(fā)輻射spontaneous emission，光在空間行進時受到空氣中微粒的影響而造成，或者是光在介質(zhì)中傳播時，受到介質(zhì)中缺陷的影響)，那麼光在進

11、入自聚焦材料時，就會產(chǎn)生有的地方較亮，有的地方較暗的情形，雖然這亮暗的差別非常非常的小，但在進入自聚焦材料後，仍然會造成亮的地方折射率稍微大一點點，暗的地方折射率小一點點，如同我們在電磁學(xué)或光學(xué)學(xué)到的，光在前進時，會往折射率高的地方靠過去(就像波導(dǎo)，光纖內(nèi)的全反射一樣)，如此光在行進時，較亮的地方因其折射率高的緣故，會吸引旁邊較暗部分的光過來，如此一來，會造成亮的地方更亮，折射率變化更大，而能吸引更多的光，所以一開始即使只有很小的光強度的不同，在傳播一段距離後，光強的分佈就會有很大的不同，而能形成圖案，這樣的現(xiàn)象就稱為調(diào)變不穩(wěn)定性。三、最近十年的研究進展二維的空間光孤子無法被觀察到的情形一直持

12、續(xù)到90年代初，直到某些可飽合(saturable)非線性介質(zhì)被人發(fā)現(xiàn)可以用來支持二維的空間光孤子 ( 早在1969年, Daws和Marburger已經(jīng)以理論證明可飽合非線性介質(zhì)可以抵消“災(zāi)難性自聚焦”而導(dǎo)致穩(wěn)定的二維空間光孤子，此乃因為其折射係數(shù)的最大變化量會飽合而無法超過某個極限值，使得自聚焦作用有一個煞車的機制)。這些包括二階型光孤子7(quadratic solitons)，在可飽合非線性原子蒸氣中的光孤子8(solitons in saturable nonlinear atomic vapor)，和光折變空間光孤子9 (photorefractive optical spatia

13、l solitons)。 其中，二階型光孤子的發(fā)生是在非線性晶體中產(chǎn)生二階頻(second harmonic generation)時，基頻(fundamental frequency)和二階頻互相調(diào)制相位，而使得繞射作用消失；而在可飽合非線性原子蒸氣中的光孤子，則是光的頻率非常接近原子的躍遷能量，而引致非常大的非線性折射率變化，進而調(diào)制光孤子光束的波前相位；光折變空間光孤子則是首先由Segev於1992年提出，它們是由光折變所產(chǎn)生的自聚焦效應(yīng)來平衡繞射現(xiàn)象進而支持空間光孤子。大約在同一時候，澳大利亞科學(xué)家Snyder繼續(xù)Ashotovich10的想法提出以線性波導(dǎo)的方法做為研究非線性光孤子的

14、方法，他是以前後一致性(self-consistent)的方式為邏輯，假設(shè)一道光束入射非線性物質(zhì)，經(jīng)由光和物質(zhì)的非線性作用，它可以產(chǎn)生一光波導(dǎo)，而如果這光束是它所引致光波導(dǎo)的一個模(mode)，那麼這光束將受到它引致光波導(dǎo)的限制而以同樣的波形繼續(xù)向前傳播，因為一個波導(dǎo)的模在傳播中是不會改變形狀及大小，因此它所形成的光波導(dǎo)的形狀也會保持固定，而繼續(xù)使得光孤子的光束向前傳播而不繞射，如此空間光孤子就可穩(wěn)定傳播。Snyder的方法可以應(yīng)用於一維或者是二維的光孤子，且可以延申做為研究光孤子的穩(wěn)定性和交互作用的研究工具。空間光孤子及光調(diào)變不穩(wěn)定性等研究的發(fā)展在1995年後，由於找到了上述三種合適的非線性

15、光學(xué)材料，且由於觀念上的突破，而有了重大的進步，在這裏筆者選取下列幾個重點的研究方向做介紹：(1) 二維光孤子在三維空間的交互作用在1991年，Mitchell11提出兩個靠的很近而前進方向互相歪斜的空間光孤子會因為“交互作用力”(有關(guān)這作用力的解釋請參照12，這裏我們可以如此簡單解釋，每個光孤子建立了一個折射係數(shù)的勢井，而另一個光孤子受到該勢井的影響而會被吸引)的關(guān)係互相環(huán)繞，以DNA的形式繼續(xù)向前傳播。筆者於1995年時正在進行同一平面上光折射空間光孤子互相碰撞和交互作用的實驗，受到此研究的啟示，了解到光折射空間光孤子將是最好的工具，可以作為試圖觀察Mitchell所描述的光孤子環(huán)繞現(xiàn)象，

16、於是展開藉由交互作用而相互環(huán)繞的空間光孤子對的研究13(另有一關(guān)於相互環(huán)繞空間光孤子的研究，描述的是由一“光旋”14(optical vortex)分裂為兩個互相旋轉(zhuǎn)的空間光孤子，它們互相旋轉(zhuǎn)是因為“光旋”本身已帶有角動量)。在這個實驗中, 兩個相互不同調(diào)(mutually incoherent)的二維光折變空間光孤子在入射時的起始軌跡是相互歪斜的，它們之間的交互作用是由它們之間的距離及所建立的折射係數(shù)勢井來決定。當(dāng)它們的軌跡距離相對遙遠時，這兩個光孤子因吸引力不足，而互相遠離。然而，當(dāng)它們的軌跡相距適當(dāng)時，它們之間的吸引力剛好能夠平衡互相環(huán)繞所需的離心力時，它們就互相環(huán)繞(如圖四)。當(dāng)它們的

17、軌跡相距圖四：(上)兩個光孤子在空間中互相環(huán)繞的示意圖(下)實驗時，在(a)、(b)、(c)三個斷面所觀察到的兩個光孤子。太近而使吸引力太大時，它們之間的吸引力就使它們?nèi)诤显谝黄稹Ｈ绻苁惯@兩個光孤子不受外界影響的話，以上實驗暗示著在光孤子的交互作用系統(tǒng)中，總角動量是守恆的，以上的論述，已在理論及數(shù)值模擬中獲得証實15。(2) 部分非同調(diào) (partially incoherent)空間光孤子在1996以前，所有的空間光孤子的光源都是屬於同調(diào) (coherent)的，也就是說，知道光孤子一點的相位， 我們都可以知道光孤子任何其它一點的相位，此外要觀察到光孤子，因其光源為同調(diào)光源，所以就必須在實

18、驗上使用雷射光源，但在1996年時，非同調(diào)光也可以用來作光孤子首度被提出16， 部分非同調(diào)空間光孤子之所以引起注意是因為它們有豐富的統(tǒng)計光學(xué)性質(zhì)，從第一個部分非同調(diào)空間光孤子被發(fā)現(xiàn)到目前已有“白光光孤子”17(white light soliton)，“部分非同調(diào)性空間暗光孤子”18(partially incoherent dark soliton)，及相關(guān)的理論和這些光孤子的交互作用的分析和數(shù)值模擬被發(fā)表，這些研究除了致力了解部分非同調(diào)空間光孤子的交互作用的動態(tài)學(xué)外，並研究部分非同調(diào)空間光孤子的交互作用會對空間光孤子的同調(diào)性(即光的統(tǒng)計性質(zhì))有何影響。(3) 光孤子計算(computati

19、on with solitons)1998年，Princeton University的計算機科學(xué)家Steiglitz等19首先提出了以空間光孤子，可以作為計算的媒介，他的做法是將許多個空間光孤子以不同的角度入射到非線性光學(xué)介質(zhì)中，經(jīng)由光孤子的交互作用後，在非線性光學(xué)介質(zhì)的出射面，光孤子的排列和能量的分佈會是入射光孤子的排列和能量分佈的函數(shù)，這樣的對應(yīng)的關(guān)係可以看做是一個計算(computation)的結(jié)果。Steiglitz等的貢獻是提出了這樣的觀念，並且設(shè)計了一個演算法實現(xiàn)這樣一個光孤子計算。至於在實驗驗証方面，一個最近被觀察到的實驗結(jié)果20已為這個光孤子計算演算法打下了基礎(chǔ)，該實驗的做法

20、是以兩個向量光孤子(向量光孤子，Vector Solitons 是由Manakov於1974年提出21，光孤子可以由兩個成分以向量和的方式來組成，例如：兩個垂直極化方向的光以向量和的方式來組成向量光孤子)做碰撞，而這兩個光孤子在非線性交互作用下，能量能由一個向量光孤子傳遞到另一個光孤子，該演算法就是以這樣的能量傳遞來做為計算的媒介。(4) 非同調(diào)光的調(diào)變不穩(wěn)定性雖然光調(diào)變不穩(wěn)定性的研究始於1960年代，但是在2000年時，Segev等22提出，既然非同調(diào)光可以作為光孤子的光源，那麼也應(yīng)該可以用非同調(diào)光來觀察光調(diào)變不穩(wěn)定性，而且，因為非同調(diào)光的繞射角度較大，(註：同調(diào)性越差，繞射的角度越大)，所

21、以要觀察到週期一樣的光調(diào)變不穩(wěn)定性條紋所須要的自聚焦非線性的強度就必須越大(即同樣的光強度，可以得到較大的折射率變化)。他們以理論和實驗證實這樣的想法，還發(fā)現(xiàn)新的臨界現(xiàn)象，也就是說，在某個固定自聚焦非線性強度，當(dāng)光的同調(diào)性較小時，並不會有光調(diào)變不穩(wěn)定性的圖案條紋發(fā)生，而在光的同調(diào)性大於某個值時，光調(diào)變不穩(wěn)定性的圖案條紋會突然出現(xiàn)，而且這個臨界值還跟自聚焦非線性的強度有關(guān)。除此之外，這項研究也預(yù)測了具有相吸碰撞力(attractive collision force) 的原子氣(atomic gas) 在溫度略高於波色-愛因斯坦凝聚態(tài)(Bose-Einstein condensation)臨界溫

22、度時，調(diào)變不穩(wěn)定性也可能會發(fā)生，因為：(1)原子氣的原子間具有相吸碰撞力時，描述波色-愛因斯坦凝聚態(tài)物質(zhì)波的薛丁格方程式(Gross-Pitaevski equation23)，和在自聚焦非線性材料中的電磁波的波方程式有一模一樣的形式，其波函數(shù)的解也應(yīng)該相同。(2) 波色-愛因斯坦凝聚態(tài)的亮孤子和暗孤子，也已被觀察報導(dǎo)，而不管在非線性光學(xué)中或其他非線性系統(tǒng)中，孤子和調(diào)變不穩(wěn)定性的形成原因是非常接近的，唯一的差別在於初始條件，既然波色-愛因斯坦凝聚態(tài)孤子的現(xiàn)象已被發(fā)現(xiàn)，那麼調(diào)變不穩(wěn)定性也應(yīng)該可以被觀察到。(3)在略高於波色-愛因斯坦凝聚態(tài)的臨界溫度時，原子氣的組成原子的波函數(shù)為基態(tài)(ground

23、 state) 和激發(fā)態(tài)的疊加，和部分非同調(diào)光類似。四、結(jié)語空間光孤子和光調(diào)變不穩(wěn)定性在60年代被提出後，因為未發(fā)現(xiàn)適合的非線性光學(xué)材料的關(guān)係，而使得其在70至80年代時的進展減緩了一陣子，但是這個領(lǐng)域在90年代因為新材料及新現(xiàn)象的發(fā)現(xiàn)而又有了非?？焖俚倪M步，雖然這裏我們主要介紹的是光在非線性材料中傳播的行為，但是在許多其他的非線性系統(tǒng)中，如水面波(water surface wave)24, 電漿中的離子聲波(ion-acoustic waves in plasma), 鐵磁性薄膜中傳遞的自旋微波(microwave spin waves in ferrite films)等25,26，或具

24、有相吸碰撞力的波色-愛因斯坦凝聚態(tài)系統(tǒng)，也都觀察到了孤子或調(diào)變不穩(wěn)定性的現(xiàn)象，也許不久的將來，我們也可以在這些系統(tǒng)裏觀察到類似本文所介紹的光在非線性材料中傳播的行為和現(xiàn)象，也可以在這些系統(tǒng)上出現(xiàn)。參考資料：1 J. S. Russell, in "14th meeting of the British Association Reports", York,（1844）.2 R. Y. Ciao, E. Garmire and C. H. Townes, Phys. Rev. Lett. 13, 479 (1964).3 P. L. Kelley, Phys. Rev. Le

25、tt. 15, 1005 (1965). 4 A. Hasegawa and F. Tappert, Appl. Phys. Lett. 23, 142 (1973). 5 L. F. Mollenauer, R. H. Stolen and J. P. Gordon, Phys. Rev. Lett. 45, 1095 (1980). 6 G. A. Agrawal, Nonlinear Fiber Optics, 2 ed., Academic Press, San Diego,（1995）.7 W. E. Torruellas, Z. Wang, D. J. Hagan, E. W. V

26、an Stryland, G. I. Stegeman, L. Torner and C. R. Menyuk, Phys. Rev. Lett. 74, 5036 (1995). 8 V. Tikhonenko, J. Christou, and B. Luther-Davies, Opt. Lett. 19, 1817 (1994). 9 M. Segev, B. Crosignani, A. Yariv and B. Fischer, Phys. Rev. Lett. 68, 923 (1992). 10 A. Snyder, D. Mitchell, and Y. Kivshar, M

27、od. Phys. Lett. B 9, 1479 (1995).11 D. J. Mitchell, A. Snyder, and L. Polodian, Opt. Comm. 85, 59 (1991).12 M. Segev and G. I. Stegeman, Physics Today 51(8), 42 (1998).13 M. Shih, M. Segev, and G. Salamo, Phys. Rev. Lett. 78, 2551 (1997). 14 V. Tikhonenko, J. Christou, and B. Luther-Davies, Phys. Re

28、v. Lett. 76, 2698 (1996); J. Opt. Soc. Am. 12, 2046 (1995).15 A. Buryak, V. Kivshar, M. Shih, and M. Segev, Phys. Rev. Lett. 82, 81 (1999). 16 M. Mitchell, A. Chen, M. Shih, and M. Segev, Phys. Rev. Lett. 77, 490 (1996).17 M. Mitchell and M. Segev, Nature 387, 880 (1997).18 Z. Chen, M. Mitchell, M. Segev, T. Coskun and D. Christodoulides, Science 280, 889 (1998)19 M. Jakubowski, K. Steiglitz, and R. Squier, Phys. Rev. E 58, 6752 (1998). 20 C. Anastassiou, M. Segev, K. Steiglitz, J. Giordmai