基于動力論的塵埃等離子體充電方程研究

陳惠敏 石雁祥 蘇金善 孫 毅 張保花

(1.昌吉學院物理系,新疆 昌吉831100;2.伊犁師范學院物理與電子信息學院,新疆伊寧835000)

1.引 言

塵埃等離子體可簡單定義為在一般等離子體中考慮塵埃粒子作用的等離子體。

塵埃等離子體理論是20世紀80～90年代在國際上出現的一門新型理論。它建立在一般等離子體理論基礎上但考慮了塵埃粒子的作用。由于在天體物理、航空航天、大氣環境、微波通信、微電子器件制造等一些重要的科學技術領域存在著與塵埃等離子體有關的問題,塵埃等離子體理論及相關的應用研究開始受到越來越普遍的重視。

塵埃等離子體與一般等離子體最大的不同在于塵埃粒子的充放電作用。并由此導致塵埃等離子體的很多物理性質與一般等離子體不同。塵埃粒子充電方程是描述塵埃等離子體充放過程的基本方程。而充電頻率(電荷馳豫速率)是描述這一物理過程的基本物理量。有關塵埃粒子充電問題的研究已經發表了多篇文獻[1-11]。其中,Jana和 Tsytovich分別根據各自不同的方法給出了充電頻率的兩種不同表達式[1,2]。在研究塵埃等離子體微波散射及塵埃等離子體中各種縱波問題時,不同的文獻引入了不同的充電頻率表達式。然而,對兩種表達式之間的關系至今尚沒有給出明確的說明。本文在動力論基礎上研究塵埃粒子充電問題,導出塵埃粒子充電方程,給出了充電頻率表達式。證明了兩種充電頻率表達式是考慮不同近似條件得到的。

2.理論分析

2.1 塵埃粒子充電電流的OLM理論

我們考慮浸在非磁化電子—離子等離子體中的有限尺寸的塵埃粒子。因為電子的熱速度遠大于離子的熱速度,電子到達塵埃粒子表面比離子快。因此,塵埃粒子獲得的電子比離子多很多,結果它的表面電勢變成負值。當然,原來的電子電流和離子電流會受塵埃粒子表面電勢的影響,因為它們依賴等離子體和塵埃粒子的相對速度。當表面電勢是負的時候,電子就會被排斥而離子被吸引,即由電子攜帶的塵埃粒子電流降低而由離子攜帶的塵埃粒子電流增加。

采用軌道受限運動(the orbit-limited motion,OLM)方法來計算等離子體粒子j(電子或離子)對塵埃粒子的充電電流Ij。考慮等離子j體粒子從無限遠的距離接近半徑為rd帶電量為qd的塵埃粒子的情況,如圖1。當帶電粒子進入德拜球時,它會受到塵埃粒子的影響,由于庫侖電場的作用,它的路徑會發生變化。假設Vj和Vjg分別是等離子體粒子與塵埃粒子剛好擦過時前后的速度,顯然,對于一個確定的速度,小于碰撞參數bj時,等離子j體粒子j就會與塵埃粒子相碰。其充電碰撞的碰撞截面是=。

圖1 等離子體粒子j與帶電塵埃粒子之間的摩擦碰撞(qj＜0,qd＜0)

由于帶電粒子受有心力作用,角動量和能量守恒定律要求

塵埃粒子的帶電量qd=Cφd,φd=φg-φp是塵埃粒子勢φg與等離子體勢φp之間的電勢差。在等離子體中,球形塵埃粒子的電容C=rdexp(-rd/λD)。對于德拜半徑λd?rd的情況,C≈rd。聯立求解(1)式和(2)式并考慮到關系式qd=rdφd可以得到

若距塵埃粒子為無限遠的等離子體粒子j的速度分布是fj(vj),則由j類等離子體粒子產生的塵埃充電電流Ij是

式中,Vjmin是與塵埃粒子碰撞的等離子體粒子的最小速度。Vjmin的值可分兩種情況討論：當qiφd＜0時,等離子體粒子和塵埃粒子相互吸引,Vjmin=0;而當qjφd＞0時,等離子體粒子和塵埃粒子由于庫侖力作用而相互排斥,只有那些動能大于相互作用勢能的充電粒子才能接近塵埃粒子而對其充電,在這種情況下Vjmin可表示為[6]

進一步假設等離子體粒子的速度分布服從Maxwell分布,即[12]

式中,nj是等離子體粒子濃度。把式(5)、(6)代入式(4),用球坐標表示并進行積分可以得到充電電流為

當qjφd＜0 時 ,當qjφd＞0 時 ,

2.2 充電平衡條件和電中性條件

當粒子停止充電并達到穩定狀態時的塵埃粒子表面電勢稱為平衡電勢。塵埃粒子的充電平衡電勢可由充電平衡條件和準中性條件得到。若已知塵埃粒子的半徑,進一步還可得到塵埃粒子的帶電量。

由于離子質量遠大于電子質量,離子的熱運動速度要遠小于電子的熱運動速度,因而,在充電過程初期,離子的充電電流Ii要遠小于電子的充電電流Ie,充電使塵埃粒子帶上了負電。由于塵埃粒子表面帶負電,對于電子qjφd＞0;對于離子qiφd＜0,由式(7)和式(8)可知,電子和離子的充電電流分別可寫為

假定了離子是一次電離的,即對于離子qj=e(e=1.6×10-19C)。隨著充電的繼續,帶負電的塵埃粒子開始排斥電子同時吸引離子,這使得離子的充電電流開始增加,電子的充電電流開始減小。當Ii=Ie時,塵埃粒子處于充電平衡狀態。在多種帶電粒子對塵埃粒子充電的情形,充電平衡條件可以寫為

由于塵埃粒子的帶電,塵埃等離子體的電中性條件與一般等離子體不同,其表現形式為ne+Zd nd=ni,即

式中：Zd是塵埃粒子電荷數;nd是塵埃粒子濃度。當Zd nd/ni?1,塵埃粒子是孤立的,此時ni≈ne。隨著塵埃粒子濃度增大,塵埃粒子間距減小,塵埃粒子之間的庫侖相互作用將會逐漸增大,當Zdnd/ni與1可以比擬時,塵埃粒子不能被看作是孤立的。對于非孤立塵埃粒子,其表面電勢φg與等離子體電勢φp之差小于孤立塵埃粒子時的情形,因此其塵埃粒子的電量也小于孤立塵埃時的情形。這是由于當塵埃粒子增多時,雖然塵埃粒子總體會占有更多的電子,但每個塵埃粒子分得的電子卻會減少。

3.理論推導

3.1 塵埃粒子充電方程

在電子和離子對塵埃粒子充電的過程中,由于受到塵埃等離子體中縱波的影響,充電電流也會產生波動,結果會引起塵埃粒子表面電量的起伏。假定不存在波動時,塵埃粒子帶電量的平衡值為qd0,則縱波的影響使塵埃粒子帶電量圍繞qd0起伏,表示為qd0+qd1,qd1是塵埃粒子受縱波影響產生的電量起伏值,qd1遠小于qd0。充電頻率νch描述了塵埃粒子表面電量恢復到平衡值時的快慢程度。對于由電子﹑一次電離的離子以及塵埃粒子組成的塵埃等離子體,充電平衡時,由充電平衡條件式可得

用|Ie0|表示充電平衡時電子和離子的充電電流,則

有縱波存在時等離子體j粒子的分布函數為

fj0是未受擾動時j粒子的分布函數。fj1表示分布函數的微擾。

縱波存在時的充電截面為

式中：φd0是充電平衡時的塵埃粒子電勢與等離子體背景電勢之差;?φd是圍繞平衡電勢變化的浮動電勢。于是,有擾動時的充電電流為

式中

為沒有擾動時的充電電流。式(17)其余三項分別代表充電截面變化、分布函數變化對擾動電流的貢獻,以及充電截面變化及分布函數變化共同引起的擾動電流。

3.2 塵埃粒子的充電頻率

考慮到fα1?fα0,φd?φd0,式(17)等式右邊最后一項為二階小量,可以忽略。于是

式中

代表分布函數變化對擾動電流的貢獻。

在忽略塵埃粒子運動速度的情形下,塵埃電量qd與充電電流的關系為

將塵埃粒子電量寫為qd=qd0+qd1。qd0為充電平衡時的塵埃電量,它不隨時間發生變化。將式(19)代入式(21),并考慮式(13)得

式(22)右端第二項的求和包括電子e和離子i。以下討論式(22)右邊積分,將式(6)代入該積分可得

對于電子,將=代入上式

得

對于離子,qjφd＜0,將Vjmin=0代入式(23)得

考慮到式(14),可得

設塵埃粒子電容為C,則,將式(26)代入式(22)得

式中,νch是充電頻率[2],在一些文獻里也被稱為電荷弛豫速率[8-11]

通常也將式(27)寫為

式(29)即為充電方程的常見形式[8]。

假設塵埃粒子帶負電,并設電荷數為Z,則φd0=-Ze2/rd,將式(14)|Ie0|的具體形式代入式(28),得

這正是文獻[2]給出的結果。式(28)和式(31)常出現在以往的文獻中,但各種文獻對兩者的關系沒有做出說明,這里給出進一步的說明。

4.結 論

塵埃等離子體與等離子體最大的不同在于塵埃粒子的充電作用。在動力論基礎上推導塵埃粒子充電方程,給出的充電頻率具有統一形式表達式,文獻中給出的兩種不同形式的充電頻率是在不同近似條件下給出的。

[1] JANA M R,SEN A,KAW P K,et al.Collective effects due to charge-fluctuation dynamics in a dusty plasma[J].Phys.Rev.E,1993,48(4)：3930-3933.

[2] TSYTOVICH V N,HAVNES O.Charging processes,Dispersion properties and anomalous transport in dusty plasma[J].Comments Plasma Phys.Controlled Fusion,1993,15(5)：267-280.

[3] HAVNEs O.Charges on dust particles[J].Adv.Space Res.,1984,4(9)：75-83.

[4] GOREE J.Ion trapping by a charged dust grain in a plasma[J].Phys.Rev.Lett.,1992,69(2)：277-280.

[5] ALLEN J E.Probe theory-theorbital motion approach[J].Phys.Scr.,1992,45：497-503.

[6] MA J X,YU M Y.Langmuir wave instability in a dusty plasma[J].Phys.Rev.E,1994,50：2431-2433.

[7] BARKAN A,D＇ANGELO N,MERLINO R L,et al.Charging of dust grains in a plasma[J].Phys.Rev.Lett.,1994,73(23)：3093-3096.

[8] SHUKLA P K,MAMUN A A,Introduction to dusty plasma physics[M].London：IOP Publishing Ltd,2002.

[9] 石雁祥,葛德彪,吳 健.環境壓強對固體火箭噴焰電導率的影響[J].電波科學學報,2008,23(6)：1115-1118.SHI Yanxiang,GE Debiao,WU Jian.Influence of pressure of environment on conductivity of solid rocket exhaust plume[J].Chinese Journal of Radio Science,2008,23(6)：1115-1118.(in Chinese)

[10] 薛 昆,郭立新,吳 健,等.高緯極區電離層離子速度分布函數的多項式解[J].電波科學學報,2009,24(5)：813-819.XUE Kun,GUO Lixin,WU Jian,et al.The polynomial solutions of the ion velocity distributions in the high-latitude auroral ionosphere[J].Chinese Journal of Radio Science,2009,24(5)：813-819.(in Chinese)

[11] 石雁祥,王 菊,吳 健,等.對兩種弱電離塵埃等離子體特征參量的定量估計[J].電波科學學報,2008,23(1)：95-99.SHI Yanxiang,WANG Ju,WU Jian,et al.Characteristic parameters estimation of two weakly ionized dusty plasma[J].Chinese Journal of Radio Science,2009,24(5)：813-819.(in Chinese)

[12] ALEXANDROV A F,BOGDANKEVICH L S,RUKHADZE A A,et al.Principles of Plasma Electrodynamics[M].Berlin：Springer,1984：59.