激發(fā)光線寬對(duì)原子光致漂移速率的影響?

楊佳琦 劉加?xùn)| 劉濤? 張志忠

1 引 言

光致漂移現(xiàn)象是存在于光場(chǎng)中的一種氣體動(dòng)力學(xué)效應(yīng),其原理是:基于多普勒效應(yīng),利用激光對(duì)某種氣體粒子進(jìn)行“速度選擇性激發(fā)”,使得在激光傳播方向上具有特定速度分量的該種氣體粒子有一部分布居到激發(fā)態(tài)上,而具有其他速度分量的粒子保持基態(tài)不變.該種氣體粒子處于緩沖氣體的環(huán)境中,由于兩種能態(tài)上的粒子與緩沖氣體的碰撞輸運(yùn)特性有差異,使得氣體粒子在選擇的速度方向上產(chǎn)生宏觀的定向移動(dòng)[1,2].光致漂移現(xiàn)象在天體物理[3?6]、半導(dǎo)體物理[7]和同位素分離[8?12]等領(lǐng)域均具有廣泛的應(yīng)用.

光致漂移速率是表征光致漂移效應(yīng)的一個(gè)關(guān)鍵參數(shù),定義為分布在各個(gè)能級(jí)上的吸收氣體粒子沿激光傳播方向投影速度分量的平均值[13].光致漂移現(xiàn)象一經(jīng)發(fā)現(xiàn),研究者們便圍繞漂移速率的理論研究開展了大量的工作,其中大多數(shù)理論工作把堿金屬原子作為研究對(duì)象[13?21].對(duì)于光致漂移而言,激發(fā)光線寬會(huì)影響原子激發(fā)的速度選擇性,進(jìn)而影響漂移速率的大小,因此是光致漂移研究中不可忽視的一個(gè)影響因素.在國(guó)外已有的理論研究中,通常把激發(fā)光作為單色光進(jìn)行處理,很少考慮激發(fā)光線寬對(duì)漂移速率的影響[13?19].1981年,Popov等[20]首次發(fā)表與激發(fā)光線寬相關(guān)的理論研究,但該研究是采用二能級(jí)原子模型對(duì)原子光致漂移過(guò)程進(jìn)行描述,沒有考慮能級(jí)簡(jiǎn)并與超精細(xì)能級(jí)對(duì)光致漂移的影響.柴俊杰等[21]的研究中對(duì)線寬因素的影響稍有涉及,但針對(duì)碰撞問(wèn)題采用的是碰撞弛豫時(shí)間近似的數(shù)學(xué)處理方法,這種方法非常粗略,只適合對(duì)漂移速率問(wèn)題做一些簡(jiǎn)單的定性分析.

本文在理論建模時(shí),統(tǒng)一考慮能級(jí)簡(jiǎn)并、超精細(xì)結(jié)構(gòu)以及碰撞模型等各個(gè)因素,以文獻(xiàn)[13]中堿金屬原子光致漂移速率方程理論為基礎(chǔ),引入能級(jí)簡(jiǎn)并與超精細(xì)能級(jí)的影響,利用強(qiáng)碰撞模型描述原子與緩沖氣體的碰撞作用,建立四能級(jí)速率方程模型描述光致漂移過(guò)程中原子的能級(jí)躍遷.并運(yùn)用數(shù)值方法對(duì)速率方程進(jìn)行計(jì)算,根據(jù)計(jì)算結(jié)果分析了激發(fā)光線寬對(duì)堿金屬原子漂移速率的影響.

2 原子光致漂移速率方程理論

通常堿金屬原子的光譜結(jié)構(gòu)較為簡(jiǎn)單,基于堿金屬原子基態(tài)超精細(xì)分裂寬度與多普勒展寬基本在一個(gè)量級(jí),所以計(jì)算時(shí)要考慮基態(tài)超精細(xì)結(jié)構(gòu)能級(jí)的影響.而激發(fā)態(tài)的超精細(xì)分裂寬度遠(yuǎn)小于激發(fā)態(tài)精細(xì)結(jié)構(gòu)的能級(jí)差以及原子的多普勒展寬,故在建立速率方程時(shí)可以忽略激發(fā)態(tài)的超精細(xì)結(jié)構(gòu).基于上述原因,在建立理論模型時(shí)主要圍繞圖1所示的4個(gè)能級(jí)之間原子躍遷進(jìn)行討論.其中能級(jí)1和2代表原子基態(tài)的超精細(xì)能級(jí)結(jié)構(gòu),能級(jí)3和4代表原子激發(fā)態(tài)的精細(xì)結(jié)構(gòu).Aij(i=3,4;j=1,2)為自發(fā)輻射速率,hij為速度選擇性激發(fā)的速率,ζij為碰撞作用因子.

圖1 原子的簡(jiǎn)化能級(jí)結(jié)構(gòu)Fig.1.Simplifi ed energy level structure of atoms.

在光致漂移過(guò)程中,堿金屬原子除了經(jīng)歷基態(tài)原子的激發(fā)、激發(fā)態(tài)原子的自發(fā)輻射和受激輻射等過(guò)程,還會(huì)與其他原子產(chǎn)生碰撞.由于堿金屬原子密度遠(yuǎn)小于緩沖氣體密度,此處忽略堿金屬原子之間的碰撞,僅考慮堿金屬原子同緩沖氣體之間的碰撞.堿金屬原子與緩沖氣體之間的碰撞過(guò)程會(huì)造成速度的改變和能級(jí)的躍遷,主要包括:1)碰撞前后能級(jí)不變,僅速度發(fā)生改變;2)激發(fā)態(tài)精細(xì)結(jié)構(gòu)能級(jí)上的原子由于碰撞引起的無(wú)輻射躍遷.由于基態(tài)超精細(xì)結(jié)構(gòu)能級(jí)之間的碰撞躍遷截面較小,此處忽略其碰撞躍遷過(guò)程.

綜合上述分析并結(jié)合圖1,可以得到光致漂移過(guò)程中堿金屬原子各能級(jí)粒子數(shù)布居的速率方程如下[13]:

式中ρi(v)(i=1,2,3,4)為能級(jí)i上原子的速度分布函數(shù)(其中v指激光傳播方向上的粒子速度);t為時(shí)間;Ai=Ai1+Ai2(i=3,4);速度選擇性激發(fā)的速率hij表達(dá)式為[13]

式中Bij為與受激輻射相關(guān)的愛因斯坦系數(shù),c為光速,Γ為均勻展寬,f為激發(fā)光頻率,fij為i能級(jí)躍遷至j能級(jí)時(shí)吸收光子的頻率,k為波數(shù),I(f)為激發(fā)光功率密度的譜分布,計(jì)算中認(rèn)為譜分布是高斯線型[22],

式中I0=I(f)d f,?f為激發(fā)光功率密度譜分布的半高寬;碰撞作用因子ζij表示由于碰撞引起的原子i→j能級(jí)的躍遷以及速度大小分布的變化.

由于碰撞散射角度的不同,碰撞作用中又包含大角度散射(large angle scattering,LAS)和小角度散射(small angle scattering,SAS)兩項(xiàng)[13,14].強(qiáng)碰撞模型理論認(rèn)為,LAS項(xiàng)中原子經(jīng)碰撞后的速度滿足麥克斯韋分布,SAS項(xiàng)中原子碰撞前后速度不發(fā)生變化[13].對(duì)于原子的兩個(gè)基態(tài)超精細(xì)能級(jí),由于忽略其碰撞躍遷過(guò)程,故只有LAS項(xiàng)起作用.對(duì)于原子的兩個(gè)激發(fā)態(tài)精細(xì)能級(jí),碰撞作用項(xiàng)既包含LAS項(xiàng)又包含SAS項(xiàng).故強(qiáng)碰撞模型下的速率方程可表示為[13]:

式中γij為L(zhǎng)AS碰撞頻率(γi= γii),Γij為SAS碰撞頻率.ni為能級(jí)i上的總粒子數(shù),ni=∫ρi(v)d v,各能級(jí)的粒子數(shù)密度滿足歸一化公式:

W(v)為一維麥克斯韋分布函數(shù):

式中v0為原子的最概然速度.

強(qiáng)碰撞模型認(rèn)為兩個(gè)基態(tài)超精細(xì)結(jié)構(gòu)能級(jí)上的原子具有相同的LAS碰撞頻率.LAS碰撞頻率γij計(jì)算公式如下:

式中g(shù)3與g4分別為3,4能級(jí)的簡(jiǎn)并度,m為堿金屬原子的質(zhì)量,kB為玻爾茲曼常數(shù),T為溫度,Dg為基態(tài)原子的擴(kuò)散系數(shù),De為激發(fā)態(tài)原子的擴(kuò)散系數(shù).

SAS碰撞頻率Γij計(jì)算公式如下:式中Γitj為碰撞躍遷頻率,Γitj=Nvσij;N 為緩沖氣體的原子數(shù)密度;σij為由i能級(jí)至j能級(jí)的碰撞躍遷截面;v為堿金屬原子與緩沖氣體原子的平均運(yùn)動(dòng)速度.

根據(jù)漂移速率的定義,漂移速率vdr的表達(dá)式為[13]:

在穩(wěn)態(tài)的情形下,令方程(4)中的

運(yùn)用數(shù)值方法求解穩(wěn)態(tài)條件下的速率方程(4),并結(jié)合原子數(shù)密度的歸一化公式(5),可解得原子的速度分布函數(shù)ρi(v),進(jìn)而可根據(jù)漂移速率的定義式(9)獲得光致漂移速率vdr.

3 數(shù)值計(jì)算結(jié)果與討論

在進(jìn)行數(shù)值計(jì)算時(shí),選擇堿金屬原子7Li作為研究對(duì)象,氙氣作為緩沖氣體.計(jì)算所采用的參數(shù)如下:中心吸收波長(zhǎng)λ23=670.791 nm,λ24=670.776 nm,λ13=670.790 nm,λ14=670.775 nm;能級(jí)簡(jiǎn)并度g1=3,g2=5,g3=8,g4=16;A3=A4=3.77× 107s?1;σ34=1.2×10?18m2,σ43=6.6× 10?19m2. 在溫度T=600 K,壓強(qiáng)p=1.01×105Pa條件下Dg=0.77 cm2/s,De=0.52 cm2/s.

不同激發(fā)光線寬條件下光致漂移速率隨激發(fā)光波長(zhǎng)的變化曲線如圖2,靠近每條曲線處的文字表示激發(fā)光線寬.由圖2可以看出,無(wú)論激發(fā)光線寬如何,漂移速率隨波長(zhǎng)的變化趨勢(shì)是一致的,隨著波長(zhǎng)的變化,在正方向和反方向上分別有兩個(gè)漂移速率的極大值.在波長(zhǎng)分別逐漸偏離原子的中心吸收波長(zhǎng)670.776 nm和670.791 nm時(shí),漂移速率的值均呈現(xiàn)先增大后減小的趨勢(shì),且漂移速率為極大值時(shí)對(duì)應(yīng)的激發(fā)光波長(zhǎng)也近乎一致.

光致漂移速率隨激發(fā)光線寬的變化曲線如圖3.由圖3可以看出,隨著線寬的增大,漂移速率的值呈現(xiàn)先增大后減小的趨勢(shì).存在一個(gè)最佳的激發(fā)光線寬,在該線寬條件下原子的漂移速率達(dá)到最大值.從光致漂移的原理出發(fā),給出的解釋如下[1,2]:單頻激光只能有效激發(fā)某一特定速度附近的少量原子,具有其他速度的大部分原子則不能被有效激發(fā),導(dǎo)致原子的漂移速率較低;當(dāng)線寬逐漸增大時(shí),更大速度范圍內(nèi)的原子會(huì)被激發(fā),激發(fā)效率的提高使原子的漂移速率增大;但是隨著線寬持續(xù)增大,在激光功率密度一定的條件下,線寬的增加會(huì)導(dǎo)致譜功率密度下降,使得處于激光線寬內(nèi)的原子受到的激發(fā)作用減弱;此外,線寬過(guò)大也會(huì)導(dǎo)致反向運(yùn)動(dòng)的原子被激發(fā).因此,線寬增大到一定值后,原子漂移速率開始呈現(xiàn)下降的趨勢(shì).由圖3還可看出,最佳線寬點(diǎn)右側(cè)曲線變化較左側(cè)緩慢.該現(xiàn)象告訴我們,在進(jìn)行相關(guān)實(shí)驗(yàn)時(shí),如果激光線寬在最佳線寬附近出現(xiàn)不可控的波動(dòng),可以設(shè)置激發(fā)光線寬使之略大于最佳線寬,這樣可以減少線寬的隨機(jī)波動(dòng)對(duì)漂移速率的影響.

圖2 不同激發(fā)光線寬條件下漂移速率隨激發(fā)光波長(zhǎng)變化曲線 計(jì)算參數(shù)為I0=0.63 W/cm2,p=300 Pa,T=600 KFig.2.Drift velocity as a function of wavelength at diff erent linewidths.Calculation parameters:I0=0.63 W/cm2,p=300 Pa,T=600 K.

不同激發(fā)光線寬條件下光致漂移速率隨激發(fā)光功率密度的變化曲線如圖4.靠近每條曲線處的數(shù)字表示激發(fā)光線寬的大小.由圖4可知,無(wú)論線寬如何,隨著功率密度的增大,漂移速率均呈現(xiàn)先增大后減小的趨勢(shì).當(dāng)線寬增加時(shí),獲得最大漂移速率所對(duì)應(yīng)的最佳功率密度逐漸減小.該現(xiàn)象告訴我們,在以后進(jìn)行相關(guān)實(shí)驗(yàn)時(shí),可采取適當(dāng)增大激光線寬的方法減小最佳激光功率密度,降低激光成本.由圖4還可看出,激發(fā)光線寬越小,在最佳功率密度附近漂移速率隨功率密度變化越平緩.該現(xiàn)象告訴我們,在激光線寬較小時(shí),激光功率密度在最佳值附近有一定偏差不會(huì)對(duì)漂移速率產(chǎn)生很大影響.激光線寬較大時(shí),則需嚴(yán)格設(shè)置激發(fā)光功率密度.

圖3 漂移速率隨激發(fā)光線寬變化曲線 計(jì)算參數(shù)為λL=670.78 nm,I0=0.63 W/cm2,p=300 Pa,T=600 KFig.3. Drift velocity as a function of linewidth.Calculation parameters:λL = 670.78 nm,I0=0.63 W/cm2,p=300 Pa,T=600 K.

圖4 不同激發(fā)光線寬條件下光致漂移速率隨激光功率密度的變化曲線 計(jì)算參數(shù)為λL=670.78 nm,p=300 Pa,T=600 KFig.4.Drift velocity as a function of laser power density at diff erent linewidths.Calculation parameters:λL=670.78 nm,p=300 Pa,T=600 K.

最佳線寬隨激發(fā)光功率密度的變化曲線如圖5.由圖5可知,在功率密度很小時(shí),最佳線寬隨功率密度的增大呈現(xiàn)急劇增大的趨勢(shì),待功率密度達(dá)到一定值后,最佳線寬又開始下降,且下降速度極為緩慢.

圖5 最佳線寬隨激發(fā)光功率變化曲線 計(jì)算參數(shù)為λL=670.78 nm,p=300 Pa,T=600 KFig.5.Best linewidth as a function of laser power density.Calculation parameters:λL=670.78 nm,p=300 Pa,T=600 K.

圖7 最佳線寬隨緩沖氣體壓強(qiáng)變化曲線 計(jì)算參數(shù)為λL=670.78 nm,I0=0.63 W/cm2,T=600 KFig.7. Best linewidth as a function of pressure.Calculation parameters:λL = 670.78 nm,I0=0.63 W/cm2,T=600 K.

最佳線寬隨溫度的變化曲線如圖6.由圖6可知,最佳線寬隨著溫度的升高近似線性增大.這是因?yàn)闇囟壬邥?huì)造成原子吸收譜線多普勒展寬的增大,在總功率不變的情況下,激光需要在更寬的頻率范圍內(nèi)分配功率才能獲得最佳漂移效應(yīng).

最佳線寬隨緩沖氣體壓強(qiáng)的變化曲線如圖7.由圖7可知,最佳線寬隨著緩沖氣體壓強(qiáng)的升高而逐漸減小.且壓強(qiáng)越高,最佳線寬下降越緩慢.該現(xiàn)象可以解釋為:壓強(qiáng)越低,緩沖氣體與堿金屬原子的碰撞頻率越低,碰撞弛豫作用越弱,在總功率不變的情況下,需要增大激光線寬提升激光的速度選擇性,從而獲得最佳漂移效應(yīng).

4 結(jié) 論

本文在考慮激發(fā)光線寬的情況下,采用速率方程理論對(duì)堿金屬原子的漂移速率進(jìn)行了研究.得到以下結(jié)論:存在一個(gè)最佳的激發(fā)光線寬,在該線寬條件下原子漂移速率達(dá)到最大值;激發(fā)光功率密度、溫度和緩沖氣體壓強(qiáng)會(huì)影響最佳線寬的值;在進(jìn)行相關(guān)實(shí)驗(yàn)時(shí),為了獲得最大漂移速率,激發(fā)光應(yīng)當(dāng)工作在最佳線寬的條件下;當(dāng)激光線寬在最佳線寬附近波動(dòng)時(shí),設(shè)置激光線寬略大于最佳線寬可減少線寬隨機(jī)波動(dòng)對(duì)漂移速率的影響.本文的研究結(jié)果為以后開展相關(guān)實(shí)驗(yàn)研究提供了理論基礎(chǔ),對(duì)于實(shí)驗(yàn)中激發(fā)光的選擇具有一定的參考價(jià)值.

[1]Gel’mukhanov F K,Shalagin A M 1979 Pis’ma Zh.Eksp.Teor.Fiz.29 773

[2]Antsigin V D,Atutov S N,Gel’mukhanov F K,Telegin G G,Shalagin A M 1979 Pis’ma Zh.Eksp.Teor.Fiz.30 262

[3]Shalagin A M 1989 Sov.Phys.Usp.32 281

[4]Leiblanc F,Michaud G 1993 Astron.J.408 251

[5]Aret A,Sapar A 2002 Astron.Nachr.323 21

[6]Sapar A,Aret A,Sapar L,Poolam?e R 2009 New Astron.Rev.53 240

[7]Shalaev V M,Douketis C,Moskovits M 1992 Phys.Lett.A 169 205

[8]Chapovsky P L,Shalagin A M 1981 Opt.Commun.40 129

[9]Streater A D,Mooibroek J,Woerdman J P 1987 Opt.Commun.64 137

[10]Streater A D,Mooibroek J,Woerdman J P 1988 Appl.Phys.Lett.52 602

[11]Gangrsky Y P,Hradecny C,Slovak J,Thethal T,Yermolayev I M 1992 Phys.Lett.A 168 230

[12]Atutov S N,Kolinko P V,Shalagin A M 1993 Laser Phys.3 855

[13]Streater A D,Woerdman J P 1989 J.Phys.B:At.Mol.Opt.Phys.22 677

[14]Kryszewski S,Nienhuis G 1987 J.Phys.B:At.Mol.Phys.20 3027

[15]Dubetsky B Y 1985 Zh.Eksp.Teor.Fiz.88 1586

[16]Nienhuis G 1985 Phys.Rev.A 31 1636

[17]Gel’mukhanov F K,Il’ichov L V,Shalagin A M 1986 J.Phys.A:Math.Gen.19 2201

[18]Werij H G C,Haverkort J E M,Planken P C M,Eliel E R,Woerdman J P,Atutov S N,Chapovskii P L,Gel’mukhanov F K 1987 Phys.Rev.Lett.58 2660

[19]Haverkort J E M,Werij H G C,Woerdman J P 1988 Phys.Rev.A 38 4054

[20]Popov A K,Shalagin A M,Shalaev V M,Yakhnin V Z 1981 Appl.Phys.25 347

[21]Chai J J,Chen R S,Xu W Q 2015 Acta Optica Sin.35 0102001(in Chinese)[柴俊杰,陳日升,許文強(qiáng) 2015光學(xué)學(xué)報(bào)35 0102001]

[22]Zhou B K,Gao Y Z,Chen T R,Chen J H 2014 Laser Principle(Beijing:National Defence Industry Press)pp134–137(in Chinese)[周炳琨,高以智,陳倜嶸,陳家驊2014激光原理(北京:國(guó)防工業(yè)出版社)第134—137頁(yè)]