李愛仙

(北京應(yīng)用物理與計算數(shù)學(xué)研究所,北京 100088)

0 引言

非線性光學(xué)的研究,包括三波混頻、四波混頻、高次諧波的產(chǎn)生,以及拉曼及超拉曼散射等非線性過程的研究,仍是當(dāng)今超強光物質(zhì)耦合方面的研究熱點[1-2]。分子/原子和固體,在強場的驅(qū)動下,都表現(xiàn)出了豐富有趣的非線性光學(xué)特性。特別是半導(dǎo)體納米結(jié)構(gòu),比如半導(dǎo)體量子點體系,由于其可操控性比較好,為我們提供了許多探索線性光學(xué)和非線性光學(xué)的可能性。半導(dǎo)體量子點在基礎(chǔ)量子力學(xué)的理解上起著十分重要的作用,在功能器件方面也有很好的應(yīng)用前景。為了能夠更好地實現(xiàn)對量子點的非線性光學(xué)調(diào)控,本文主要研究半導(dǎo)體量子點的非線性光學(xué)特性。由于半導(dǎo)體量子點在三個維度都受到了很強的限制,因而將表現(xiàn)出十分顯著的量子受限效應(yīng),在納米尺度范圍,量子受限效應(yīng)會導(dǎo)致許多有趣的物理現(xiàn)象,比如:與尺度密切相關(guān)的吸收譜和熒光譜。半導(dǎo)體量子點還經(jīng)常被用于非線性光學(xué)新特性的研究,本文研究的體系就是半導(dǎo)體量子點,但所獲得的主要物理圖像在普通的量子體系中均有效。

半導(dǎo)體量子點體系,在弱場的驅(qū)動下,光學(xué)響應(yīng)是線性的,輻射光的頻率主要由體系的能級差來決定,輻射光譜的強度主要由振蕩強度決定,并且與驅(qū)動場的強度是成正比的[3]。特別是,輻射模式是由基于反演對稱性的選擇定則決定的。強場驅(qū)動下的非線性光學(xué)響應(yīng)則包含了更加豐富的物理。一般來說,半導(dǎo)體量子點在強場的驅(qū)動下,表現(xiàn)出來的輻射譜主要由兩部分組成:高次諧波部分和超拉曼線部分。而高次諧波和超拉曼線在太赫茲的產(chǎn)生上都有豐富的應(yīng)用[3-4]。近幾十年來,人們對輻射譜中的高次諧波部分已經(jīng)投入了特別多的關(guān)注,研究的體系包括原子、分子和納米體系[3-12]。在重新制備好的態(tài)下,對大塊固體的高次諧波已經(jīng)有了一定的研究[6]。相對而言,目前對超拉曼線的研究還不太多,主要是由于超拉曼線的輻射強度在絕大多數(shù)情況下都非常弱,在實驗上很難被觀測到。

在本文中,我們主要研究輻射譜中的超拉曼組分,以及其與高次諧波組分之間的關(guān)系,對超拉曼譜的調(diào)控主要是通過動力學(xué)手段來實現(xiàn)。我們通過將一個光學(xué)過程映射到了一個輸運過程中去,從而將量子輸運模型中的能帶寬度與光學(xué)過程中的輻射強度聯(lián)系起來。本文還指出了,為什么在一般情況下,超拉曼線都比高次諧波要弱得多的原因,以至于我們在實驗上很難觀測到超拉曼線。找到原因則可以指導(dǎo)我們調(diào)控出很強的超拉曼線。

1 模型和計算方法

我們研究的半導(dǎo)體量子點體系,在強場的驅(qū)動下可以簡化成一個二能級模型,如圖1(a)所示，采用如下的哈密頓來描述:

(1)

(2)

Fig.1 Schematic representation of our system: (a) related optical process, (b) mapping to transport圖1 系統(tǒng)示意圖:(a)相關(guān)的光學(xué)過程,(b)輸運映射示意圖

我們在電子數(shù)和光子數(shù)表象下構(gòu)建了一個映射,并獲得了如下的哈密頓量:

(3)

我們用基|i,m〉來描述電子在態(tài)|i〉上,并伴有m個模式為q頻率為νq的光子。態(tài)|1,m〉 (|2,m〉)被映射到鏈I(鏈II)的m格點上的粒子態(tài)上,如圖1(b)所示。這里,bm(dm)是描述態(tài)|1,m〉 (|2,m〉)的湮滅算符。在兩條耦合鏈的粒子輸運表象下,我們可以對輻射譜有更多更好的定性上的理解。實質(zhì)上,粒子沿鏈傳輸?shù)迷竭h(yuǎn),伴隨產(chǎn)生的光子數(shù)就越多,在輻射譜上表現(xiàn)出的輻射強度就越強。因此,一個量子光學(xué)問題就被映射到一個含時的量子輸運問題上去了。周期場驅(qū)動下的量子體系可以借助Floquet理論來求解[14]。而準(zhǔn)能譜的能帶寬度在量子輸運中的影響至關(guān)重要,也就是說,對輻射譜的分析是至關(guān)重要的。準(zhǔn)能譜的能帶寬度越寬,量子輸運就越有效,就會伴隨有越多的光子數(shù)產(chǎn)生。而準(zhǔn)能譜可以通過對時間演化算符U(T=2π/ω0)的精確對角化獲得。

通過密度矩陣來計算輻射譜也十分有效和方便,其中密度矩陣滿足的演化方程如下(下面取?=1):

(4)

其中

(5)

這里Δ=E2-E1是兩個態(tài)之間的能級間距。方程(4)中的最后一項用來描述可能的耗散效應(yīng),如自發(fā)的聲子散射等。而偶極矩的平均值則通過下式來計算[15]:

(6)

輻射譜

(7)

則可以通過快速傅里葉變換獲得。

2 理論結(jié)果與討論

在數(shù)值計算時,二能級量子體系的偶極矩的含時平均是通過解密度矩陣的方程(4)獲得的。我們采用龍格-庫塔方法來求解,選取的時間步長為0.002/ω0,總步數(shù)為3.0×106.然后,就可以通過對偶極矩做傅里葉變換來獲得輻射譜。準(zhǔn)能譜的數(shù)值計算則是通過對時間演化算符U(T)做精確對角化獲得的。我們選取驅(qū)動場的頻率?ω0作為能量的單位。并取gq=0.005,Δ/?ω0=6.0.取弛豫常數(shù)Γ=1.0×10-6ω0,弛豫時間足夠長,以保證輻射譜的主要相干特性不受破壞。接下來我們將主要介紹強場驅(qū)動下,如何實現(xiàn)超拉曼譜的調(diào)控。

弱場驅(qū)動下的線性輻射譜相對比較簡單,其光學(xué)響應(yīng)是線性的?；旧?我們只能看到兩個峰:一個與共振頻率有關(guān),出現(xiàn)在ν=Δ的位置;另一個來源于瑞利散射,散射頻率為ω0.一般情況下,強場驅(qū)動下的非線性輻射譜主要由兩部分組成,既有高次諧波,又有相應(yīng)的超拉曼線,這是由于Floquet態(tài)具有相同的對稱性造成的。通常情況下,超拉曼線的強度要比高次諧波弱得多,在實驗上很難觀測到[3],這給超拉曼譜的研究帶來了巨大的阻礙。

Fig.2 Quasi-energy spectrum of the symmetric configuration(without time-dependent diagonal term, i.e. G11=G22=0). G12/ω0=30.93圖2 對稱結(jié)構(gòu)下的準(zhǔn)能譜(沒有含時的對角元項,即G11=G22=0)。G12/ω0=30.93

Fig.3 Corresponding emission spectra of the symmetric configuration.The black line is emission spectrum for the strong driving field.The red dots are obtained by the formula log10CB2/p, for ν=p/q,where B is the bandwidth from Fig.2 and C is a constant.G12/ω0=30.93. The initial condition is ρ11(t=0)=1.黑線是強場驅(qū)動下,輻射譜的數(shù)值計算結(jié)果。帶有圓點的線是由公式log10CB2/p計算得來的,其中ν=p/q。B是從圖2中得出的能帶寬度,C是常數(shù)。G12/ω0=30.93。初始條件是ρ11(t=0)=1。圖3 對稱結(jié)構(gòu)下對應(yīng)的輻射譜

圖3并不是一個特例,在其它的參數(shù)下,能帶的寬度與輻射譜的這種關(guān)系依然是存在的,如圖4和圖5所示。其中超拉曼線的位置移動也是通過動力學(xué)方法來實現(xiàn)的。而準(zhǔn)能譜的能帶寬度隨拉比頻率的變化顯示,一般情況下,高次諧波譜對應(yīng)的能帶寬度要比超拉曼譜對應(yīng)的能帶寬度大得多。正因如此,高次諧波的強度要比超拉曼線的強度大得多。因此,在非線性光學(xué)實驗中,超拉曼線很難被觀測到,很多情況下都忽略不計。

因此,要實現(xiàn)超拉曼譜的調(diào)控,首先應(yīng)抑制高次諧波譜的強度,而這可以通過動力學(xué)調(diào)控的方法產(chǎn)生能帶塌縮來實現(xiàn)[3]。在能帶塌縮的條件下,高次諧波的輻射譜將被強烈抑制,從而使輻射譜上的超拉曼譜相對變強,便于我們觀察。我們可以通過粒子輸運圖像來理解相應(yīng)的光子輻射過程。基本上,能帶塌縮就意味著粒子的有效質(zhì)量變得無限大,變得很難移動,這將導(dǎo)致動力學(xué)局域化,其本質(zhì)原因是由隧穿的相干破壞造成的。能帶塌縮和其相關(guān)的動力學(xué)局域化效應(yīng)是周期場驅(qū)動下粒子輸運中產(chǎn)生的一種有趣現(xiàn)象,目前已經(jīng)有了許多相關(guān)的研究[16-18]。發(fā)生能帶塌縮時,粒子將無法移動到其他格點上,也就無法輻射出光子,導(dǎo)致了相應(yīng)輻射峰的消失。

Fig.4 Quasi-energy spectrum of the symmetric configuration(without time-dependent diagonal term, i.e.,G11=G22=0). G12/ω0=32.77.圖4 對稱結(jié)構(gòu)下的準(zhǔn)能譜(沒有含時的對角元項,即G11=G22=0)。G12/ω0=32.77

Fig.5 Corresponding emission spectra of the symmetric configuration.The black line is emission spectrum for the strong driving field.The red dots are obtained by the formula log10CB2/p, for ν=p/q,where B is the bandwidth from Fig.4 and C is a constant.G12/ω0=32.77. The initial condition is ρ11(t=0)=1.黑線是強場驅(qū)動下,輻射譜的數(shù)值計算結(jié)果。帶有黑色圓點的線是由公式log10CB2/p計算得來的,其中ν=p/q。B是從圖4中得出的能帶寬度,C是常數(shù)。G12/ω0=32.77。初始條件是ρ11(t=0)=1。圖5 對稱結(jié)構(gòu)下對應(yīng)的輻射譜

圖6(a)是準(zhǔn)能譜的能帶寬度隨拉比頻率G12的變化曲線,既包括輻射光頻率是整數(shù)的高次諧波,也包括分?jǐn)?shù)的超拉曼譜。從圖中可見,v=nω0對應(yīng)的能帶寬度一般都要比超拉曼線對應(yīng)的能帶寬度大得多。這也正好解釋了在非線性光學(xué)實驗中很難觀測到超拉曼譜的原因。然而,在一些特殊點,如G12/ω0≈30.35,32.95,34.55時,v=1,3的能帶寬度將趨近于零,這些點對應(yīng)于能帶塌縮的情況。在發(fā)生能帶塌縮條件下,高次諧波的輻射譜將被強烈抑制,圖6(b)和(c)分別對應(yīng)v=1和3的情況,我們看到,位于v=1和3處的高次諧波峰消失了,輻射譜上只有超拉曼線的峰。總之,輻射譜中的高次諧波組分和超拉曼線組分的位置和強度均是由準(zhǔn)能譜決定的。而選擇定則是基于時空對稱性的。基于增頻轉(zhuǎn)換和降頻轉(zhuǎn)換,高次諧波和超拉曼線都可以應(yīng)用于太赫茲波的產(chǎn)生。

Fig.6 (a)The quasi-energy band width for v=1,3,1/2. The lines with black squares,red dots and blue triangles corresponding to v=1,3,1/2, respectively.(b) and (c) are the corresponding emission spectra.(b) G12/ω0=31.37 for ν=1, (c) G12/ω0=31.35for ν=3. G11=G22=0. The initial condition is ρ11(t=0)=1.(a)中帶方塊的黑線、帶圓點的線和帶三角的線分別對應(yīng)于v=1,3,1/2的情況。(b)G12/ω0=31.37對應(yīng)v=1的情況,(c) G12/ω0=31.35對應(yīng)于v=3的情況。G11=G22=0。初始條件是ρ11(t=0)=1。圖6 v=1,3,1/2的準(zhǔn)能能帶寬度(a)和能帶塌縮時的輻射譜(b)、(c)

3 結(jié)論

通過將一個量子光學(xué)過程映射到一個量子輸運過程,找出了分?jǐn)?shù)拉曼譜的準(zhǔn)能譜的能帶寬度與輻射譜之間的對應(yīng)關(guān)系。進(jìn)一步推動我們找出采用動力學(xué)的方法,不僅可以實現(xiàn)對準(zhǔn)能的能帶寬度的調(diào)控,還可以調(diào)控分?jǐn)?shù)準(zhǔn)能的位置,從而調(diào)節(jié)超拉曼線峰的位置。即超拉曼線峰的位置也可以通過入射場來調(diào)節(jié)[4]。尤其是對于一些特殊的外場,可以實現(xiàn)高次諧波譜的能帶塌縮,從而使我們可以觀測到超拉曼譜。