高溫等離子體的狀態(tài)方程及其熱力學(xué)性質(zhì)

湯文輝徐彬彬 冉憲文 徐志宏

(國(guó)防科學(xué)技術(shù)大學(xué)理學(xué)院,長(zhǎng)沙 410073)(2016年10月17日收到;2016年12月17日收到修改稿)

專題:高壓下物質(zhì)的新結(jié)構(gòu)與新性質(zhì)研究進(jìn)展

高溫等離子體的狀態(tài)方程及其熱力學(xué)性質(zhì)

湯文輝?徐彬彬 冉憲文 徐志宏

(國(guó)防科學(xué)技術(shù)大學(xué)理學(xué)院,長(zhǎng)沙 410073)(2016年10月17日收到;2016年12月17日收到修改稿)

高溫下等離子體的狀態(tài)方程及其熱力學(xué)性質(zhì)在天體物理、可控核聚變以及武器設(shè)計(jì)與破壞效應(yīng)等領(lǐng)域有著廣泛應(yīng)用.本文主要回顧了高溫等離子體在不同狀態(tài)區(qū)域的狀態(tài)方程的理論模型和處理方法.對(duì)于理想等離子體,離子之間的相互作用可以忽略,其狀態(tài)方程較簡(jiǎn)單,已趨于完善.在超高溫下,原子完全電離,離子和電子都可以采用理想氣體狀態(tài)方程描述;當(dāng)溫度不太高時(shí),離子部分電離,可以采用Saha方程及其修正模型描述;原子在高度壓縮狀態(tài)下,其狀態(tài)方程可以采用Thomas-Fermi模型及其改進(jìn)模型得到.對(duì)于非理想等離子體,離子之間存在強(qiáng)耦合,還沒(méi)有單一的理論模型能夠在任意密度和溫度范圍內(nèi)對(duì)離子之間的相互作用進(jìn)行統(tǒng)一描述.量子分子動(dòng)力學(xué)方法原則上可以在較大溫度密度范圍內(nèi)給出可靠結(jié)果,但由于計(jì)算量太大以及高溫下的計(jì)算存在收斂問(wèn)題,也較難應(yīng)用到溫度較高的稠密等離子體區(qū)域.半經(jīng)驗(yàn)的經(jīng)典分子動(dòng)力學(xué)方法雖然簡(jiǎn)單、計(jì)算量小,但只能在一定的區(qū)域范圍內(nèi)給出較精確的狀態(tài)方程結(jié)果.在不同溫度密度區(qū)域內(nèi)采用不同的計(jì)算模型,再在空白區(qū)域進(jìn)行插值從而得到全局狀態(tài)方程在目前不失為一種簡(jiǎn)單有效的方法.

高溫等離子體,狀態(tài)方程,熱力學(xué)性質(zhì)

1 引 言

自然界中大約90%—95%的可觀測(cè)物質(zhì)處于極端高溫高壓條件下[1].太平洋底部(海平面下11 km處)壓強(qiáng)為0.12 GPa;地球中心處的壓強(qiáng)為360 GPa,溫度為0.5 eV(1 eV=11604.5 K),密度約為10—20 g/cm3;木星中心壓強(qiáng)約為4—6 TPa,溫度約為2 eV,密度約為30 g/cm3;太陽(yáng)中心壓強(qiáng)為2.4×107GPa,溫度為1.6×103eV,密度約為150 g/cm3;白矮星壓強(qiáng)為1012—1018GPa,密度為106—109g/cm3,溫度為103eV.在慣性約束聚變(ICF)中,中心點(diǎn)火需要氘氚(DT)熱斑溫度達(dá)到5 keV以上,密度100 g/cm3以上,壓強(qiáng)達(dá)到1.5×107—2.0×107GPa.目前已知的自然界的物態(tài)范圍大致如圖1[1]所示(1 bar=105Pa).然而,在實(shí)驗(yàn)室條件下只能產(chǎn)生相對(duì)較低的壓強(qiáng)和溫度.靜高壓加載能達(dá)到550 GPa的壓強(qiáng)[2],結(jié)合激光加熱技術(shù)可使實(shí)驗(yàn)溫度能夠達(dá)到3000 K以上[3].動(dòng)高壓加載雖然能產(chǎn)生104GPa的瞬態(tài)高壓,但持續(xù)時(shí)間短,實(shí)驗(yàn)測(cè)量非常困難,測(cè)量結(jié)果誤差較大[4].事實(shí)上,人們真正了解的物質(zhì)狀態(tài)范圍只是圖1中很小的一個(gè)區(qū)域,研究更大溫度密度范圍的物質(zhì)狀態(tài)需要發(fā)展更先進(jìn)的實(shí)驗(yàn)手段和理論方法.

狀態(tài)方程是描述物質(zhì)系統(tǒng)中各狀態(tài)變量之間關(guān)系的函數(shù)表達(dá)式,反映物質(zhì)在一定熱力學(xué)條件下的宏觀性質(zhì)[1,4].近代自然科學(xué)和工程技術(shù)中的大量實(shí)際問(wèn)題,比如天體物理、恒星以及行星內(nèi)部構(gòu)造、慣性約束聚變、武器系統(tǒng)的設(shè)計(jì)及其破壞效應(yīng)等問(wèn)題,都涉及到物質(zhì)處于高溫、高壓、高密度等極端條件,要對(duì)這些問(wèn)題進(jìn)行研究,就需要了解物質(zhì)在各種極端條件下的狀態(tài)方程及其熱力學(xué)性質(zhì)[5].

眾所周知,高溫條件下的物質(zhì)以等離子體的形式存在.隨著科學(xué)技術(shù)的發(fā)展以及國(guó)防科技的驅(qū)動(dòng),等離子體的熱力學(xué)性質(zhì)及其描述方法已成為目前的研究熱點(diǎn)之一.然而,等離子體的熱力學(xué)性質(zhì)及其描述方法與其所處壓強(qiáng)和溫度狀態(tài)密切相關(guān),難度很大.為了便于人們對(duì)等離子體的熱力學(xué)性質(zhì)開展深入研究,本文在文獻(xiàn)報(bào)道的基礎(chǔ)上對(duì)高溫(1 eV以上)等離子體的狀態(tài)方程模型及其熱力學(xué)性質(zhì)進(jìn)行了綜述.

圖1 自然界和實(shí)驗(yàn)室中的物態(tài)范圍[1],圖中括號(hào)里的數(shù)字表示密度的對(duì)數(shù)值,密度單位為g/cm3;“Statics”區(qū)域?qū)?yīng)靜高壓實(shí)驗(yàn)所能達(dá)到的溫度壓強(qiáng)范圍,“Dynamics”區(qū)域表示動(dòng)高壓加載實(shí)驗(yàn)所能達(dá)到的溫度壓強(qiáng)范圍Fig.1.Extreme states in natu re and in the laboratory[1],the numbers in parentheses ind icate the logarithmof density(g/cm3).The “statics” domain corresponds tostatic techniques for high-pressure production,the“dynamics” domain tothe dynamic ones.

2 等離子體狀態(tài)區(qū)域的劃分

等離子體可以根據(jù)離子間相互作用的強(qiáng)弱劃分為理想等離子體和非理想等離子體[6].定義ΓD為離子之間相互作用能與離子動(dòng)能的比值,ΓD可用于評(píng)估離子之間耦合的強(qiáng)弱,稱為離子耦合系數(shù).ΓD?1表示離子之間的相互作用非常小,這時(shí)可以認(rèn)為等離子體是理想等離子體;ΓD～1表示離子之間的相互作用不容忽略,等離子體是非理想等離子體;ΓD?1時(shí),離子間存在較強(qiáng)的相互作用,離子間的相互作用能占主導(dǎo)地位.

另一方面,根據(jù)粒子量子效應(yīng)的強(qiáng)弱可以將等離子體分為強(qiáng)簡(jiǎn)并等離子體和弱簡(jiǎn)并等離子體,粒子簡(jiǎn)并度參數(shù)可用α表示如下:

式中,n是粒子數(shù)密度,λˉ是電子德布羅意波長(zhǎng),?是約化Planck常數(shù),m是粒子質(zhì)量,k是Boltzmann常數(shù),T是溫度.α=n3?1表示粒子間距遠(yuǎn)大于粒子的德布羅意波長(zhǎng),粒子處于非簡(jiǎn)并或者弱簡(jiǎn)并態(tài),量子效應(yīng)不明顯,這時(shí)粒子可以作為經(jīng)典粒子處理,遵循經(jīng)典的Boltzmann統(tǒng)計(jì).當(dāng)粒子間距較小時(shí),粒子處于強(qiáng)簡(jiǎn)并態(tài)(α?1),服從Fermi量子統(tǒng)計(jì).相對(duì)于電子,離子的質(zhì)量較大,德布羅意波長(zhǎng)非常小,因此離子的量子效應(yīng)只有在極低溫情況下才需要考慮.對(duì)于質(zhì)量最小的氫離子,在固體密度狀態(tài)下,溫度需要降低到幾十開爾文(K),氫離子的量子效應(yīng)才需要考慮.對(duì)于慣性約束聚變中心點(diǎn)火熱斑,氘氚燃料被壓縮到原來(lái)固體密度的1000倍以上,溫度到達(dá)5 keV以上[7],由于氘氚離子滿足α?1的條件,因此中心氘氚熱斑不需要考慮氘氚離子的量子效應(yīng),但對(duì)于電子,其量子效應(yīng)不能忽略.

為了方便,我們根據(jù)離子耦合系數(shù)以及電子簡(jiǎn)并度參數(shù)的大小,把高溫等離子體分成四個(gè)區(qū)域進(jìn)行討論,如圖2所示.I區(qū)和II區(qū)分界線是原子完全電離所需溫度,II區(qū)和IV區(qū)是根據(jù)離子相互作用耦合強(qiáng)弱劃分的,I區(qū)和III區(qū)是根據(jù)電子簡(jiǎn)并強(qiáng)弱劃分的.由于不同元素完全電離溫度不一致,圖2中的區(qū)域劃分主要適用于低Z元素,對(duì)于高Z元素,由于原子完全電離所需溫度增加,因此I區(qū)和II區(qū)的分界線需要適當(dāng)上移.I區(qū)是完全電離的等離子體,在溫度非常高的情況下,原子中所有電子電離,這時(shí)粒子之間的相互作用可以忽略不計(jì),電子也處于非簡(jiǎn)并狀態(tài)(ΓD?1,α?1),滿足經(jīng)典Boltzmann統(tǒng)計(jì),此區(qū)域的等離子體中電子和原子核都可以采用理想氣體狀態(tài)方程描述.II區(qū)是部分電離弱耦合等離子體,溫度不是特別高,部分電子處于自由態(tài),部分處于束縛態(tài).此區(qū)域密度較低,電子和電子之間以及電子和離子之間的相互作用較弱(ΓD? 1,α? 1),但在靠近 ΓD=1的區(qū)域,需要考慮粒子間的遠(yuǎn)程Coulomb相互作用.III區(qū)是強(qiáng)簡(jiǎn)并等離子體,物質(zhì)密度非常高,原子處于高度壓縮狀態(tài),電子分布在原子核周圍,對(duì)原子核產(chǎn)生屏蔽,原子核與原子核的相互作用可以忽略,電子處于強(qiáng)簡(jiǎn)并狀態(tài)(ΓD?1,α?1).IV區(qū)是強(qiáng)耦合等離子體,物質(zhì)密度較高,溫度相對(duì)較低,離子之間耦合作用非常強(qiáng),電子處于部分電離部分簡(jiǎn)并狀態(tài)(ΓD?1或者ΓD～1,α～1).I區(qū)、II區(qū)以及III區(qū)的等離子體由于離子間相互作用可以忽略不計(jì),因此可以當(dāng)作理想等離子體處理,現(xiàn)有的理論模型能夠較好地對(duì)這些區(qū)域的等離子體狀態(tài)方程進(jìn)行描述,而IV區(qū)是強(qiáng)耦合等離子體,離子之間的相互作用還沒(méi)有一個(gè)單一的理論模型能夠在任意密度和溫度范圍內(nèi)進(jìn)行統(tǒng)一的描述.

原子序數(shù)為Z的離子間Coulomb作用能為

式中

其中,M為原子摩爾質(zhì)量,ρ為密度,NA為Avogadro常數(shù),e為電子電荷量,r是粒子間的平均距離,r～ n?1/3.

單個(gè)粒子的平均動(dòng)能為

對(duì)于I區(qū)和II區(qū)的稀薄理想等離子體,ΓD?1.如果離子之間主要是Cou lomb相互作用,則耦合系數(shù)滿足條件

經(jīng)過(guò)變換,得到滿足稀薄理想等離子體的條件為

式中溫度單位為K.由(6)式可知,氣體在常態(tài)密度情況下,溫度只需要幾個(gè)eV就可以滿足稀薄氣體條件;對(duì)于鋰、鈹?shù)容^輕的元素,常態(tài)密度下,溫度需要達(dá)到幾十eV才能滿足稀薄氣體條件;對(duì)于重元素,則需要上百eV的溫度才能滿足稀薄氣體條件.對(duì)于完全電離的稀薄等離子體,電子和離子都滿足經(jīng)典Boltzmann條件,這時(shí)可以把電子和離子都作為理想氣體處理.對(duì)于溫度較低的稀薄理想等離子體,其狀態(tài)方程和電離度相關(guān),而電離狀態(tài)可以采用Saha方程[8]進(jìn)行描述.雖然帶電粒子間的長(zhǎng)程Coulomb作用能相對(duì)粒子動(dòng)能較小,但Cou lomb相互作用引起附加的自由能,會(huì)導(dǎo)致電離能減小,因此在密度稍高或者溫度較低的情況下,也就是ΓD?1條件滿足較弱的情況下,需要考慮粒子間長(zhǎng)程Coulomb作用的影響,這時(shí)的狀態(tài)方程可以采用Debye-Hückel模型[9]計(jì)算.

圖2 溫度密度平面上等離子體區(qū)域的劃分,其中,I區(qū)為完全電離等離子體,II區(qū)為部分電離區(qū)弱耦合等離子體,III區(qū)為強(qiáng)簡(jiǎn)并等離子體,IV區(qū)為強(qiáng)耦合等離子體;圖中粗實(shí)線是電子簡(jiǎn)并強(qiáng)弱分界線,即α=1,粗虛線表示離子耦合強(qiáng)弱分界線,即ΓD=1,細(xì)實(shí)線表示四個(gè)區(qū)域的分界線,Sun以及ICF分別表示太陽(yáng)和慣性約束聚變條件下氫的同位素所處的大致區(qū)域Fig.2.D ivision of plasma on the density and temperatu re region:Iis the fu lly ionized plasma,IIis partly ionized and week coupling plasma,IIIis the strong degenerate plasma and IV is the strong coupling plasma.The thick solid lines are the boundary ofweak and strong degenerate of electron(α=1)and the thick dashed are the boundary ofweak and strong coupling of ion(ΓD=1).The thin solid lines are the boundary of four regions.Sun and ICF mean the regions of hyd rogen’s isotopes under the condition of sun and inertial con finement fusion.

III區(qū)為強(qiáng)簡(jiǎn)并等離子體,是原子被高度壓縮所形成的.原子的內(nèi)層電子被均勻擠壓,電子處于強(qiáng)簡(jiǎn)并態(tài),滿足關(guān)系,電子的分布可以使用Thomas-Fermi(TF)模型[10,11]進(jìn)行描述.TF模型認(rèn)為原子中電子連續(xù)分布在原子核周圍[12,13],原子核被連續(xù)的電子云所包圍,電子是一團(tuán)簡(jiǎn)并氣體,并遵循Fermi量子統(tǒng)計(jì).由于分布在原子核周圍高密度電子的屏蔽作用,原子核與原子核之間的相互作用可以忽略.當(dāng)原子中電子處于強(qiáng)簡(jiǎn)并狀態(tài)時(shí),費(fèi)米能為Ek～2n2/3/(2m),于是ΓD正比于n?1/3.粒子數(shù)密度n越大時(shí),ΓD越小,因此電子強(qiáng)簡(jiǎn)并條件下的等離子體在高度壓縮條件下也趨于理想等離子體.

對(duì)于IV區(qū)的強(qiáng)耦合等離子體,離子之間存在較強(qiáng)的相互作用(ΓD～1或者ΓD?1).由于電子處于部分簡(jiǎn)并狀態(tài),這時(shí)單純采用經(jīng)典Boltzmann統(tǒng)計(jì)模型或者Fermi量子統(tǒng)計(jì)模型都不能正確地描述電子的統(tǒng)計(jì)狀態(tài).物質(zhì)處于稠密狀態(tài)時(shí),相鄰原子之間的距離非常小,離子間存在較強(qiáng)的相互作用,因此必須正確描述離子間的相互作用才能準(zhǔn)確描述物質(zhì)的熱力學(xué)性質(zhì).分子動(dòng)力學(xué)方法可以較好地描述粒子之間的相互作用,但前提是需要知道粒子之間的相互作用勢(shì).根據(jù)所使用的原子間相互作用勢(shì)的不同,分子動(dòng)力學(xué)方法可以分為兩大類:一類是經(jīng)典分子動(dòng)力學(xué)方法(MD),粒子之間相互作用采用經(jīng)驗(yàn)勢(shì)或半經(jīng)驗(yàn)勢(shì)進(jìn)行描述;而另一類是量子分子動(dòng)力學(xué)方法(QMD),原子間相互作用勢(shì)基于量子力學(xué)理論計(jì)算得到,比如基于密度泛函理論計(jì)算的分子動(dòng)力學(xué)方法[14].兩種分子動(dòng)力學(xué)方法各有優(yōu)缺點(diǎn).經(jīng)典分子動(dòng)力學(xué)方法由于直接使用經(jīng)驗(yàn)勢(shì),因此計(jì)算方法簡(jiǎn)單,計(jì)算量小;但由于不同物質(zhì)的相互作用勢(shì)差別較大,一種經(jīng)驗(yàn)勢(shì)一般只適用于一種原子系統(tǒng)或一部分系統(tǒng),而且經(jīng)驗(yàn)勢(shì)也局限于一定的密度溫度范圍,因此使用某種經(jīng)驗(yàn)勢(shì)的經(jīng)典分子動(dòng)力學(xué)方法只適用于特定溫度密度范圍內(nèi)的某些物質(zhì).量子分子動(dòng)力學(xué)方法不用假設(shè)粒子間的相互作用勢(shì),而是由第一性原理計(jì)算原子相互作用勢(shì),所以具有普適性,但由于在高溫時(shí)計(jì)算量較大,而且難以收斂,因此主要適用于溫度不太高的狀態(tài).在計(jì)算電子結(jié)構(gòu)時(shí)為了節(jié)約機(jī)時(shí),量子分子動(dòng)力學(xué)常采用贗勢(shì)的方法,因此在溫度較高時(shí),電子躍遷和電子激發(fā)的處理也存在困難,這時(shí)需要采用一些半經(jīng)驗(yàn)的分子動(dòng)力學(xué)或改進(jìn)的量子分子動(dòng)力學(xué)模型來(lái)描述稠密物質(zhì)的狀態(tài).總之,等離子體處于溫稠密和熱稠密狀態(tài)時(shí),難以采用單一的理論模型對(duì)其狀態(tài)方程進(jìn)行精確描述,一般需要采用多種理論模型相結(jié)合,各自描述一部分密度和溫度范圍,然后通過(guò)插值得到過(guò)渡段的狀態(tài)方程.

3 高溫稀薄等離子體

對(duì)于理想氣體,不考慮原子之間的相互作用力,其狀態(tài)方程可表示為

式中P為壓強(qiáng);ε為比內(nèi)能;n為粒子數(shù)密度;γ為氣體比熱比,與氣體的內(nèi)部自由度相關(guān),如果氣體有q個(gè)自由度,則γ=1+2/q.對(duì)于單原子分子q=3,γ=5/3;對(duì)于雙原子分子q=5,γ=7/5.

高溫條件下,分子先離解成原子,隨后原子電離成各級(jí)離子.分子的離解以及原子的電離都會(huì)導(dǎo)致氣體的粒子數(shù)密度發(fā)生明顯變化.當(dāng)溫度升高到原子完全電離(圖2中I區(qū)),等離子體的粒子(包括離子與電子)數(shù)密度由原來(lái)的n變?yōu)?1+Z)n,Z為原子的核電荷數(shù).因此,原子完全電離的高溫稀薄等離子體狀態(tài)方程可以表示為

對(duì)于部分電離的稀薄等離子體(圖2中II區(qū)),其狀態(tài)方程可采用Saha方程[8]描述.溫度不是很高時(shí),原子并未完全電離,帶電粒子的平均動(dòng)能遠(yuǎn)遠(yuǎn)大于粒子之間的Coulomb相互作用能,這時(shí)Saha方程可以很好地描述其電離行為.Saha方程適用的范圍是ΓD?1,對(duì)于原子序數(shù)為Z的粒子,稀薄氣體條件即是(6)式.在低密度或者高溫情況下都可以滿足Saha方程條件,但在密度較高或者溫度較低情況下(圖2中II區(qū)靠近IV區(qū)位置),電離氣體系統(tǒng)中帶電粒子間的長(zhǎng)程Coulomb作用不可忽視.Coulomb相互作用引起附加的自由能,使電離能減小.對(duì)于這樣的系統(tǒng),必須對(duì)Saha模型添加Debye-Hückel近似進(jìn)行修正.

3.1 稀薄電離氣體的Saha模型

大部分分子在1 eV以上時(shí)已經(jīng)離解為原子,有些原子已經(jīng)發(fā)生了電離,隨著溫度的進(jìn)一步升高,較重原子所組成的氣體會(huì)多次電離.考察某種原子所組成的氣體,單位質(zhì)量氣體中含有N個(gè)原子,單個(gè)原子的總能量為E=Ek+Q,其中Ek為單個(gè)原子的動(dòng)能,包括離子以及電子的動(dòng)能;Q為原子內(nèi)電子的勢(shì)能.在氣體足夠稀薄的情況下,電子遵循Boltzmann統(tǒng)計(jì)規(guī)律,這時(shí)氣體中每一個(gè)粒子(包括中性原子、離子以及電子)的平均熱能是(3/2)kT,因此N個(gè)原子系統(tǒng)總的熱能為

其中αe為平均電離度.

用Im表示第m次激發(fā)電子所必須消耗的能量,則從原子中激發(fā)m個(gè)電子需要消耗的總能量為

在給定溫度T以及密度ρ的情況下,單位質(zhì)量氣體中含有N0個(gè)中性原子,N1個(gè)一次電離的離子,Nm個(gè)電離m次的離子,以及Ne個(gè)電子.此外電離m次的離子還具有電子激發(fā)能Hm,因此氣體總的內(nèi)能為

(11)式中,αe是氣體的平均電離度,即Ne/N;αm是電離m次離子的濃度,表示為Nm/N.由于氣體比較稀薄,在粒子之間相互作用可以忽略的條件下,所有粒子都可以當(dāng)作理想氣體,因此可以得到電離氣體的壓強(qiáng)為

氣體的壓強(qiáng)與總的粒子數(shù)密度相關(guān),總的粒子數(shù)密度等于電子數(shù)密度與離子數(shù)密度之和.等離子體電離度越大,其壓力越高.電離度可以通過(guò)原子數(shù)守恒、電荷守恒以及Saha方程聯(lián)立求得

3.2 考慮輻射超高溫等離子氣體

當(dāng)?shù)入x子體氣體溫度非常高或者密度非常低時(shí),這時(shí)輻射壓以及輻射能相對(duì)于氣體的壓強(qiáng)以及能量并不能被忽略,需要把輻射的能量以及壓強(qiáng)加入到氣體的能量和壓強(qiáng)中.平衡輻射條件下,輻射的比能量等于輻射能量密度除以物質(zhì)的密度,即

輻射壓強(qiáng)可以表示為

輻射壓和溫度的四次方成正比,而根據(jù)(12)式,物質(zhì)壓和溫度的一次方成正比,因此當(dāng)溫度超過(guò)一定范圍之后,輻射壓將完全占據(jù)主導(dǎo)地位.令輻射壓等于物質(zhì)壓,可以得到兩種壓強(qiáng)平衡時(shí)溫度和密度之間的關(guān)系式為

式中密度單位取g/cm3.對(duì)標(biāo)準(zhǔn)密度的空氣來(lái)說(shuō),當(dāng)空氣溫度超過(guò)100 eV時(shí),輻射壓起主要作用;對(duì)于密度為1 g/cm3的固體來(lái)說(shuō),輻射壓要在溫度達(dá)到幾keV以上時(shí)才起主要作用.

如果輻射能流和物質(zhì)能流為同一量級(jí),這時(shí)輻射輸運(yùn)在能量輸運(yùn)過(guò)程中必須考慮.輻射能量以光速進(jìn)行傳播,而物質(zhì)能量以物質(zhì)的速度進(jìn)行傳播,在輻射溫度不是很高的情況下,也需要對(duì)輻射能流進(jìn)行考慮.同樣令物質(zhì)能流等于輻射能流,假設(shè)物質(zhì)速度為光速的千分之一,則在固體密度條件下,當(dāng)溫度在幾十eV時(shí)需要考慮輻射輸運(yùn)過(guò)程,這時(shí)不能單純使用流體動(dòng)力學(xué)過(guò)程描述物質(zhì)的運(yùn)動(dòng),而需要使用輻射流體動(dòng)力學(xué)模型來(lái)考察物質(zhì)的運(yùn)動(dòng).

3.3 沖擊壓縮下的等離子體氣體

在單次沖擊壓縮條件下,當(dāng)沖擊波強(qiáng)度足夠大時(shí),氣體的極限壓縮度可以表示為

對(duì)于單原子理想氣體,γ=5/3,h=4,即氣體密度最多只能被壓縮4倍;對(duì)于雙原子理想氣體,γ=7/5,h=6.可見,氣體比熱比γ越接近于1,其壓縮度也就越大.

在沖擊波壓縮下,波后物質(zhì)的壓強(qiáng)、密度與比內(nèi)能之間的關(guān)系由Hugoniot方程給出:

等離子體的內(nèi)能可以分為兩部分,即E=Ek+Q,其中Ek為粒子的平均動(dòng)能,Q包含了粒子的勢(shì)能以及粒子內(nèi)部自由度的能量.由于壓強(qiáng)只和粒子的平動(dòng)動(dòng)能相關(guān),即P=(2/3)ρEk,代入(18)式,在強(qiáng)沖擊波條件下,忽略沖擊波波前壓力P0和內(nèi)能E0,將ρ/ρ0用h代替,得到

從(19)式可以看出,勢(shì)能Q相對(duì)于熱運(yùn)動(dòng)動(dòng)能忽略不計(jì)(3Q/Ek?1)時(shí),單原子氣體的壓縮度最大值為4.在電離情況下,總的粒子數(shù)目增加,但(19)式表明,壓縮度和總的粒子數(shù)密度沒(méi)有關(guān)系,因此總粒子數(shù)目的增加并不會(huì)增加壓縮度.另一方面,電離度增加會(huì)導(dǎo)致用于電離的能量損耗增加,也就是勢(shì)能對(duì)總內(nèi)能的相對(duì)貢獻(xiàn)增加,因此原子在沖擊壓縮下發(fā)生電離時(shí),會(huì)使氣體變得更加容易壓縮.隨著原子最外層電子開始電離,到最外層電子全部電離,這個(gè)過(guò)程用于電離的能量所占總內(nèi)能的比例增加,因此這個(gè)過(guò)程壓縮度是增加的.在原子最外層電子(K層電子)全部電離而次外層電子(M層電子)還未開始電離時(shí),隨著沖擊波強(qiáng)度的增加,在相當(dāng)一段范圍內(nèi)離子不再發(fā)生電離(因?yàn)榇瓮鈱与娮与婋x能比最外層電子電離能大很多),勢(shì)能Q保持不變,而溫度升高成為內(nèi)能增加的主要形式,因此Q/Ek會(huì)減小,從而導(dǎo)致壓縮度降低.如果沖擊波強(qiáng)度足夠大,能夠?qū)е麓瓮鈱与娮与婋x,這時(shí)隨著沖擊波強(qiáng)度的進(jìn)一步增加,次外層電子會(huì)像最外層電子一樣一個(gè)個(gè)剝落,用于電離的勢(shì)能相對(duì)貢獻(xiàn)又增加,因此導(dǎo)致壓縮度再一次增大.如果是重原子氣體,原子含有多個(gè)殼層,則隨著各層電子的依次電離,沖擊壓縮條件下的壓縮度先增大再減小,然后又增大再減小.如果沖擊壓縮過(guò)后原子發(fā)生完全電離,隨著沖擊波強(qiáng)度進(jìn)一步增加,會(huì)導(dǎo)致溫度急劇上升,粒子的熱運(yùn)動(dòng)動(dòng)能急劇增加,而勢(shì)能Q保持不變,因此壓縮度最后會(huì)減小.在沖擊波強(qiáng)度非常大的情況下,如果不考慮熱輻射,根據(jù)(19)式,完全電離氣體的壓縮度也是趨近于4.

3.4 考慮粒子間Cou lomb作用的稀薄等離子體

對(duì)于滿足理想等離子體條件較弱的稀薄等離子體,粒子之間的Coulomb相互作用相對(duì)于熱運(yùn)動(dòng)能量雖然較小,但不能忽視.這時(shí)Coulomb相互作用可以采用Debye-Hückel模型進(jìn)行計(jì)算.Debye-Hückel模型假定所有帶電粒子的勢(shì)場(chǎng)都是球?qū)ΨQ的,帶電粒子之間只考慮Coulomb作用力,但仍然假定Coulomb作用能仍遠(yuǎn)小于粒子的熱運(yùn)動(dòng)動(dòng)能,因此有

或者表示為

式中?(r)表示r處的電勢(shì),ni(r)是r處第i種粒子的數(shù)密度,Zi是第i種粒子的電荷數(shù).第i種帶電粒子的勢(shì)能可以表示為

由于粒子滿足Boltzmann統(tǒng)計(jì)規(guī)律(稀薄電離氣體條件下,電子和離子都可以看作經(jīng)典氣體),粒子數(shù)密度的空間分布為

式中ni0是電荷為Zi的粒子的平均數(shù)密度.(23)式后面約等號(hào)成立的條件是粒子Coulomb作用能相對(duì)于粒子的平均動(dòng)能是小量.由于滿足電荷守恒,把(23)式代入泊松方程(21)得到

在粒子中心附近有r?λD,因此有

式中第一項(xiàng)為中心粒子本身所產(chǎn)生的勢(shì),第二項(xiàng)?′=Zie/λD是周圍其他粒子在r處所形成的勢(shì).

2015年11月底，《中共中央國(guó)務(wù)院關(guān)于進(jìn)一步推進(jìn)農(nóng)墾改革發(fā)展的意見》中要求，用3年左右時(shí)間，將國(guó)有農(nóng)場(chǎng)承擔(dān)的社會(huì)管理和公共服務(wù)職能納入地方政府統(tǒng)一管理，基本完成農(nóng)墾國(guó)有土地使用權(quán)確權(quán)登記發(fā)證任務(wù)。

在體積V中,等離子體的Coulomb相互作用能為所有粒子相互作用勢(shì)的一半,利用(22)式有

異性帶電粒子存在吸引力作用,每個(gè)離子周圍都有大量的電子,所以Coulomb作用能以及壓強(qiáng)都是負(fù)值.Cou lomb作用影響氣體的狀態(tài)主要體現(xiàn)在兩個(gè)方面,一個(gè)是減小氣體的內(nèi)能和壓強(qiáng),還有一個(gè)效應(yīng)是降低電離能,從而使得電離更加容易發(fā)生.Coulomb相互作用中的電子并不是完全自由的,由于電子還在離子的勢(shì)場(chǎng)中,其勢(shì)能并不為0,而是負(fù)值,因此把電子從原子或者離子中電離出來(lái)需要消耗的電離能會(huì)下降.電離能的減小可以表示為

在較低溫度下,?In相對(duì)于kT可以忽略不計(jì),但在上百eV高溫下,?In/In可以達(dá)到10%.在密度較高情況下,考慮粒子之間的遠(yuǎn)程Coulomb作用,壓強(qiáng)和內(nèi)能會(huì)減小,也會(huì)使得電子電離能下降,從而使電離更容易發(fā)生.

4 強(qiáng)簡(jiǎn)并等離子體

前面所討論的電離氣體,電子和離子都遵循經(jīng)典的Boltzmann統(tǒng)計(jì),但只有在溫度相當(dāng)高或者密度相當(dāng)?shù)偷那樾蜗?才能滿足經(jīng)典Boltzmann統(tǒng)計(jì)條件.物質(zhì)在高度壓縮條件下,原子內(nèi)電子高度簡(jiǎn)并,這時(shí)必須考慮電子的量子效應(yīng).對(duì)于物質(zhì)在高溫高壓條件下的物態(tài)性質(zhì),比較簡(jiǎn)潔方便的處理方法是采用TF模型[15].TF模型描述的是電子對(duì)壓強(qiáng)和內(nèi)能的貢獻(xiàn),原子核對(duì)壓強(qiáng)以及內(nèi)能的貢獻(xiàn)可以采用理想氣體狀態(tài)方程描述,

式中,na為原子核數(shù)密度.

電子對(duì)內(nèi)能的貢獻(xiàn)包括電子的動(dòng)能、電子與其他電子以及原子核之間的相互作用能.下面主要討論電子對(duì)系統(tǒng)壓強(qiáng)以及內(nèi)能的貢獻(xiàn).

4.1 費(fèi)米簡(jiǎn)并態(tài)自由電子氣體

自由電子在低密度或者溫度較高情況下服從經(jīng)典的Boltzmann統(tǒng)計(jì),但在溫度較低或者密度較高的情況下,由于電子處于簡(jiǎn)并狀態(tài),其量子態(tài)必須要考慮.零溫時(shí),自由電子處于完全簡(jiǎn)并狀態(tài),在dVdp相空間的量子數(shù)為4πp2dpdV/3,由于電子有兩個(gè)不同的自旋方向,因此在相空間內(nèi)量子數(shù)還需要乘以2.根據(jù)泡利不相容原理,每一個(gè)確定的自旋方向的量子狀態(tài)中只能有一個(gè)電子,假設(shè)體積V中N個(gè)電子緊湊地占據(jù)了動(dòng)量從0到p0的量子態(tài),因此有關(guān)系式

定義費(fèi)米溫度TF=ε0/k為

電子的簡(jiǎn)并程度也可以根據(jù)費(fèi)米溫度進(jìn)行判斷,如果電子溫度T?TF,則電子處于弱簡(jiǎn)并態(tài)或者非簡(jiǎn)并態(tài),電子可以作為經(jīng)典氣體處理;如果T?TF,則電子處于強(qiáng)簡(jiǎn)并狀態(tài),服從量子分布.比較(1)和(32)式有α∝(TF/T)3/2,可見α?1與T?TF是等價(jià)關(guān)系,都表示電子處于弱簡(jiǎn)并狀態(tài);而α?1與T?TF也是等價(jià)關(guān)系,表明電子處于強(qiáng)簡(jiǎn)并狀態(tài).因此,(1)式和(32)式都可以評(píng)判電子的簡(jiǎn)并程度,并且評(píng)判結(jié)果是一致的.

體積V中N個(gè)電子的動(dòng)能為

根據(jù)熱力學(xué)關(guān)系Td S=d E+P d V,零溫下完全簡(jiǎn)并電子氣體的壓強(qiáng)為

從(34)式可以看出,在極低溫條件下,自由簡(jiǎn)并電子的壓強(qiáng)正比于密度的5/3次方,而在超高溫情況下,電子服從經(jīng)典Boltzmann統(tǒng)計(jì),這時(shí)可以采用理想氣體狀態(tài)方程描述,壓強(qiáng)與密度的一次方成正比,這說(shuō)明電子氣體處于簡(jiǎn)并狀態(tài)時(shí)更難壓縮.高溫壓縮下原子中的電子并不是完全自由的,而是在原子核以及其他電子所形成的勢(shì)場(chǎng)中運(yùn)動(dòng),這時(shí)電子的行為可以使用TF模型進(jìn)行描述.

4.2 TF模型

單個(gè)原子中電子在原子核周圍服從Fermi-Dirac分布,分布函數(shù)為

式中,μ是電子的化學(xué)勢(shì),E是電子的能量.在TF模型中,電子在原子核以及其他電子所產(chǎn)生的勢(shì)場(chǎng)中運(yùn)動(dòng),其能量主要包括動(dòng)能以及在原子內(nèi)的靜電勢(shì)能,E=p2/(2me)?eU(r),U(r)是r處的電勢(shì),p為電子的動(dòng)量.利用電子的分布函數(shù)可以得到電子數(shù)密度為

原子內(nèi)電子總數(shù)等于核電荷數(shù),即

(37)式稱為電荷守恒條件,由此可得到化學(xué)勢(shì)μ.

在電子數(shù)密度已知條件下,原子內(nèi)的電勢(shì)可以通過(guò)泊松方程得到

(38)式就是TF方程.TF方程是二階微分方程,對(duì)其求解需要一定的邊界條件.在原子中心附近,也就是半徑趨向于0時(shí),由于電子均勻連續(xù)分布在原子核周圍,因此原子中心處的電勢(shì)全部由原子核提供,這時(shí)的邊界條件可以表示為

TF方程再加上兩個(gè)邊界條件可以封閉求解,從而得到TF模型下原子中的電子分布.由于電子的分布影響勢(shì)場(chǎng),而勢(shì)場(chǎng)又會(huì)影響電子的分布,勢(shì)場(chǎng)和電子分布耦合在一起,因此TF方程的解析解難以得到,一般采用數(shù)值計(jì)算或者近似方法求解.雖然TF方程沒(méi)有解析解,但TF方程的解具有普適性質(zhì),只要溫度和體積以TZ?4/3,VZ形式出現(xiàn),則壓強(qiáng)、內(nèi)能、比熱容以PZ?10/3,EZ?7/3,CVZ?1的形式出現(xiàn),這意味著只要知道某種元素的TF狀態(tài)方程,根據(jù)上面所得到的普適性質(zhì),可以得到任意元素的TF狀態(tài)方程.對(duì)于多種原子所組成的物質(zhì),可以利用平均原子的概念得到其TF狀態(tài)方程.平均原子序數(shù)和平均原子量可以表示為

式中xi是第i種原子數(shù)目所占總原子數(shù)目的百分比,Zi和Mi分別為第i種原子的原子序數(shù)和原子量.

TF模型不區(qū)分電子是否處于束縛狀態(tài),因此這也意味著物質(zhì)的熱力學(xué)性質(zhì)和原子的殼層結(jié)構(gòu)沒(méi)有任何關(guān)系,只和原子序數(shù)相關(guān).從TF模型的普適性也可以看出,在相同的密度條件下,隨著溫度升高,物質(zhì)的壓強(qiáng)和內(nèi)能是隨著原子序數(shù)的增加而單調(diào)增加的.由于不區(qū)分電子的殼層結(jié)構(gòu),溫度的升高只是改變?cè)又须娮用芏鹊闹匦路植?因此電子勢(shì)能的變化較緩慢,電子密度在原子體積中的分布更趨于均勻.在溫度非常高的情況下,也就是電子溫度遠(yuǎn)遠(yuǎn)大于費(fèi)米溫度(T?TF),這時(shí)電子的簡(jiǎn)并性將被消除,等離子體從圖2中的III區(qū)過(guò)渡到I區(qū),這時(shí)的狀態(tài)方程和完全電離理想氣體狀態(tài)方程一樣,壓強(qiáng)和密度成正比.

4.3 改進(jìn)的TF模型

簡(jiǎn)單的TF模型只考慮了電子之間的Cou lomb相互作用,而沒(méi)有考慮交換相互作用.電子是費(fèi)米子,服從泡利不相容原理,其所引起的自旋平行電子之間的相互作用在計(jì)算電子在原子核周圍的分布時(shí)需要考慮.Dirac[16]于1930年引入交換能對(duì)零溫TF模型進(jìn)行了修正,由此得到了Thomas-Fermi-Dirac(TFD)模型.Cowan和Ashkin[17]進(jìn)一步將TFD模型推廣到有限溫度情況.TF模型對(duì)電子的能量只考慮了電子的動(dòng)能、電子與原子核以及其他電子相互作用勢(shì)能,TFD模型還考慮了電子與所有自旋平行電子的交換能.考慮交換能后電子分布函數(shù)可以表示為

式中Eex是自由電子的交換能.由Bloch[18]估算得到的交換能表示為

如果Eex～kT或者Eex? kT,則電子交換能對(duì)TF模型的修正非常明顯,而在Eex? kT時(shí),電子交換能可以忽略,這時(shí)溫度需要滿足的條件為T? 1.37×10?15n4/3.當(dāng)物質(zhì)處于固體密度時(shí),溫度需要幾十甚至上百eV,電子交換項(xiàng)的影響才可以忽略不計(jì).與TF模型結(jié)果進(jìn)行比較,TFD模型的壓強(qiáng)在低溫下(幾十eV以下)比TF模型要低,隨著溫度的升高,兩者差別越來(lái)越小,最后會(huì)重合在一起.McCarthy[19]采用不同模型進(jìn)行計(jì)算給出了鉛在不同溫度下壓強(qiáng)隨密度的變化,如圖3所示[19].可以看出,零溫(kT=0 eV)時(shí),采用TFD模型得到的壓強(qiáng)明顯要比TF模型得到的壓強(qiáng)低,10 eV時(shí)兩者差別減小,在100 eV時(shí)已基本重合在一起.

圖3 鉛在100,10和0 eV下的壓強(qiáng)曲線:TF結(jié)果、TFD結(jié)果以及TFK結(jié)果比較[19]Fig.3.Pressu re profi le of lead at the temperatu re of 100,10 and 0 eV,the solid,dotted and dashed lines are the resu lts of TF,TFD and TFKmodels,respectively[19].

TFD模型雖然引入了交換修正,但并沒(méi)有考慮電子的量子效應(yīng).1957年,Kirzhnits[20]將電子分布函數(shù)展開,在Thomas-Fermi模型中同時(shí)引入交換修正和量子修正,形成了Thomas-Fermi-Kirzhnits(TFK)理論.TFK模型引入的修正項(xiàng)和TF模型一樣具有普適性,對(duì)于TF模型有

TFK模型對(duì)TF模型的修正項(xiàng)和TF模型的普適變化關(guān)系有所不同,其修正項(xiàng)可以表示為[21]

McCarthy[19]詳細(xì)地列出了較大溫度(TZ?4/3)和體積范圍內(nèi)(VZ)壓強(qiáng)和能量的計(jì)算結(jié)果.在計(jì)算精度要求不是很高的情況下,TFK方法的解析擬合表達(dá)式[12,22]可以為高能量密度條件下的流體模擬提供簡(jiǎn)單的狀態(tài)方程.

從圖3可以看出,在100 eV的高溫時(shí),三種模型(TF,TFD,TFK)基本重合在一起;10 eV時(shí),TF模型結(jié)果稍微偏高,但在零溫且密度較低時(shí),TF模型壓強(qiáng)明顯偏大,而TFD模型和TFK模型由于考慮了交換修正以及量子效應(yīng)修正,比TF模型要精確.

TF模型把所有電子作為準(zhǔn)自由電子處理,但實(shí)際上靠近原子核的電子受到原子核強(qiáng)烈的Coulomb作用,仍處于束縛態(tài),具有殼層結(jié)構(gòu).靠近原子核內(nèi)層電子處于束縛狀態(tài),遵從量子力學(xué)規(guī)律,而較外層的電子由于熱電離或者壓致電離較容易受激發(fā)成為自由電子,因此將原子內(nèi)的電子分成束縛電子和自由電子兩部分來(lái)處理更接近實(shí)際情況.考慮殼層結(jié)構(gòu)修正的TFS(Thomas-Fermi shell)方法[23]就是基于電子分成兩部分的設(shè)想而改進(jìn)的TF模型.

將原子體積分為兩個(gè)區(qū)域:內(nèi)層靠近原子核是電子密度較高的束縛區(qū)域;外層是密度較低的自由電子區(qū)域.外層可以采用TF模型方法處理,而內(nèi)層束縛態(tài)可以采用Wenzel-Kramers-Brillouin近似處理.TFL方法把束縛電子看成準(zhǔn)自由電子,雖然考慮了內(nèi)層束縛電子但并不求解薛定諤方程.Zink[23]分別使用TFL(Thomas-Fermi like)方法和TFS方法計(jì)算了鐵的壓強(qiáng),如圖4所示.殼層修正在低溫低密度條件下起作用,在高溫或者高密度條件下,兩者差不多完全重合.在低溫低密度時(shí),TFS結(jié)果反映了電子的凝聚現(xiàn)象,即?p/?ρ較小,因?yàn)橐恍┦`態(tài)能級(jí)處于空缺狀態(tài),外層自由電子較容易重新回到束縛態(tài).當(dāng)密度增加到一定程度時(shí),束縛態(tài)能級(jí)充滿電子,隨著壓縮繼續(xù)進(jìn)行,將會(huì)出現(xiàn)壓致電離,這時(shí)壓強(qiáng)會(huì)突然升高.于是,這種電子凝聚和壓致電離交替現(xiàn)象發(fā)生在TFS結(jié)果中.Lee和Thorsos[24]對(duì)TFS方法進(jìn)行了改進(jìn),認(rèn)為有些電子可能處于自由態(tài)和束縛態(tài)之間的共振狀態(tài),將電子分成自由電子以及束縛電子和共振態(tài)電子三項(xiàng),提出了TFSR方法.他們對(duì)Z=5到Z=30的各種元素進(jìn)行了計(jì)算,結(jié)果如圖5[24]所示.可以看出,壓強(qiáng)隨Z呈現(xiàn)周期性變化,這和TF模型壓強(qiáng)隨原子序數(shù)Z單調(diào)上升有明顯區(qū)別.

圖4 鐵在7.7和770 eV下,不同密度時(shí)的壓強(qiáng),圖中分別給出了TF模型、TFL模型以及TFS模型結(jié)果[23]Fig.4.Pressure profi le of iron with diff erent density at the temperature of 7.7 and 770 eV,the solid,dashed and dashed-dotted lines are the resu lts of TF,TFL and TFSmodels,respectively[23].

平均原子(average atom,AA)模型[25?27]也可用于高溫下計(jì)算等離子體中的電子結(jié)構(gòu).AA模型不區(qū)分離子的不同組態(tài),而是把等離子體中所有離子等效成一個(gè)平均原子.在這個(gè)原子中,電子按照Fermi-Dirac統(tǒng)計(jì)規(guī)律排布在各個(gè)單電子軌道上,忽略了電子之間不同的耦合方式所帶來(lái)的差異.AA模型電子軌道類似于氫原子軌道,軌道上的電子分布呈現(xiàn)出殼層結(jié)構(gòu).這是AA模型和TF模型的主要區(qū)別,也是其在低溫低密度下具有更高精度的原因.AA模型發(fā)展至今已有許多改進(jìn),比如不可分辨的躍遷系模型[28,29],超級(jí)躍遷系模型[30,31],細(xì)致譜模型[32]等,這些模型對(duì)電子殼層結(jié)構(gòu)考慮得越精細(xì),所計(jì)算的狀態(tài)方程越精確,但所需要的計(jì)算量也越大.

圖5 不同溫度及電子數(shù)密度下TF模型和TFSR模型壓強(qiáng)和原子序數(shù)的關(guān)系[24]Fig.5.The relations of pressure and atomic number at d iff erent density and diff erent temperature,the dotted and solid lines are the resultsof TFSR and TFmodels,respectively[24].

5 強(qiáng)耦合等離子體

前面所介紹的模型都是基于單原子假設(shè),即離子之間的相互作用能相對(duì)于離子的動(dòng)能可以忽略不計(jì),物質(zhì)的性質(zhì)可簡(jiǎn)化為求單個(gè)原子的性質(zhì).當(dāng)物質(zhì)處于密度較高,而溫度并不是特別高的溫稠密或者熱稠密狀態(tài)時(shí),滿足ΓD～1或者ΓD?1,α～1,即離子間的耦合較強(qiáng),電子處于部分電離部分簡(jiǎn)并態(tài),這時(shí)使用單原子的性質(zhì)來(lái)推導(dǎo)物質(zhì)的宏觀性質(zhì)將會(huì)帶來(lái)較大誤差.

把強(qiáng)耦合等離子體分為電子和離子兩部分進(jìn)行描述.對(duì)于電子部分的描述,可以采用TF模型、AA模型、無(wú)軌道(orbital-free,OF)模型以及量子力學(xué)理論.離子間相互作用可采用超網(wǎng)鏈(hypernetted-chain,HNC)模型[33?35]、蒙特卡羅(MC)方法以及MD方法進(jìn)行描述.結(jié)合電子的描述模型以及離子相互作用模型,涌現(xiàn)了一大批計(jì)算高溫稠密物質(zhì)物態(tài)性質(zhì)的模型方法:Ofer等[36]采用TF模型來(lái)描述電子結(jié)構(gòu),離子間相互作用采用HNC近似來(lái)處理.Furukawa和Nishihara[37]采用AA模型并考慮量子修正的 HNC近似處理離子相互作用.Zérah等[38?41]采用TF模型計(jì)算電子密度空間分布,采用MD來(lái)描述離子間相互作用,提出了TFMD模型.Dharma-wardana和Francois[42]在密度泛函框架下通過(guò)Kohn-Sham方程處理電子問(wèn)題,然后采用HNC近似來(lái)描述高溫稠密等離子體;Dharma-wardana和Murillo[43]采用雙溫模型解決高溫稠密物質(zhì)中的離子間相互作用問(wèn)題.Zérah等[44,45]結(jié)合無(wú)軌道密度泛函理論和MD方法,發(fā)展了無(wú)軌道分子動(dòng)力學(xué)方法(OFMD).侯永等[46,47]結(jié)合考慮能級(jí)展寬效應(yīng)的AA模型和MD方法提出了AAMD模型.Car和Parrinello[14]結(jié)合密度泛函理論和MD方法,提出了第一性原理分子動(dòng)力學(xué)方法(CPMD).Dai等[48?50]基于量子分子動(dòng)力學(xué)方法,考慮了高溫下電子碰撞離子對(duì)離子位置所形成的漲落,發(fā)展了能描述從低溫凝聚態(tài)到高溫理想等離子體的量子Langevin分子動(dòng)力學(xué)方法(QLMD).前面幾種基于TF模型、OF模型、AA模型的半經(jīng)驗(yàn)方法,適用于特定的溫度以及密度條件,而基于第一性原理的QMD,以及基于路徑積分的蒙特卡羅(PIMC)方法[51]被認(rèn)為是比較精確的模型.但它們所需計(jì)算量非常巨大,尤其是PIMC方法,計(jì)算量隨著電子個(gè)數(shù)急劇增加,其應(yīng)用受限于電子數(shù)很少的低Z原子.由于計(jì)算量大,QMD方法所能給出的狀態(tài)方程數(shù)據(jù)點(diǎn)較少,所以一般被用來(lái)作為評(píng)判其他半經(jīng)驗(yàn)方法是否適用的基準(zhǔn)[52?54].

隨著計(jì)算方法和超級(jí)計(jì)算機(jī)的發(fā)展,分子動(dòng)力學(xué)方法得到了迅速發(fā)展,分子動(dòng)力學(xué)方法可以有效地描述離子之間的相互作用,但前提是需要離子之間相互作用勢(shì).下面分別介紹使用經(jīng)驗(yàn)或者半經(jīng)驗(yàn)勢(shì)的經(jīng)典分子動(dòng)力學(xué)方法和利用量子力學(xué)原理直接計(jì)算相互作用勢(shì)的量子分子動(dòng)力學(xué)方法.

5.1 經(jīng)典分子動(dòng)力學(xué)方法

MD方法是一種用來(lái)計(jì)算經(jīng)典多體體系平衡和傳遞性質(zhì)的方法,其目的是模擬系統(tǒng)的微觀細(xì)致動(dòng)力學(xué)行為.在分子動(dòng)力學(xué)方法中,忽略原子核運(yùn)動(dòng)的量子效應(yīng),原子核的運(yùn)動(dòng)滿足經(jīng)典力學(xué)的牛頓定律.當(dāng)體系中所有粒子上的作用力已知時(shí),便可以通過(guò)牛頓運(yùn)動(dòng)方程得到系統(tǒng)中所有粒子的運(yùn)動(dòng)軌跡,再通過(guò)統(tǒng)計(jì)平均得到該系統(tǒng)的宏觀熱力學(xué)參數(shù).由于粒子在勢(shì)場(chǎng)中運(yùn)動(dòng),其受力等于勢(shì)函數(shù)的空間導(dǎo)數(shù),因此,從經(jīng)典力學(xué)出發(fā),求解原子間聯(lián)立牛頓運(yùn)動(dòng)方程的分子動(dòng)力學(xué)方法主要包括兩個(gè)問(wèn)題:一是粒子相互作用勢(shì)函數(shù)的選取,二是積分牛頓運(yùn)動(dòng)方程的求解.

對(duì)于一個(gè)含有N個(gè)粒子的系統(tǒng),每個(gè)粒子的坐標(biāo)、速度、動(dòng)量以及受力可以分別表示為ri(t),vi(t),pi(t),f(ri,t),系統(tǒng)的哈密頓量可以表示為

式中第一項(xiàng)是系統(tǒng)動(dòng)能,第二項(xiàng)U(rN)為系統(tǒng)勢(shì)能函數(shù),rN=(r1,r2,...,rN)代表N個(gè)原子的坐標(biāo).哈密頓方程可以表示為

(48)式說(shuō)明,離子運(yùn)動(dòng)所受到的力由離子所處的勢(shì)場(chǎng)所決定.已知離子所處的勢(shì)場(chǎng),通過(guò)(48)式加上初始條件以及周期性邊界條件,可以進(jìn)行數(shù)值求解,從而得到離子的運(yùn)動(dòng)軌跡,進(jìn)而得到各種宏觀物理量.根據(jù)統(tǒng)計(jì)物理思想,系統(tǒng)經(jīng)過(guò)足夠長(zhǎng)的時(shí)間之后會(huì)處于平衡狀態(tài),從而表征系統(tǒng)宏觀性質(zhì)的物理量A可以通過(guò)時(shí)間平均得到

在時(shí)間足夠長(zhǎng)、系統(tǒng)達(dá)到平衡態(tài)后,有

分子動(dòng)力學(xué)方法的實(shí)質(zhì)是利用計(jì)算機(jī)對(duì)系統(tǒng)中所有粒子的牛頓運(yùn)動(dòng)方程進(jìn)行數(shù)值求解,得到粒子坐標(biāo)、速度等動(dòng)力學(xué)量隨時(shí)間變化的函數(shù),然后通過(guò)統(tǒng)計(jì)平均求得表征系統(tǒng)宏觀性質(zhì)的物理量.因此分子動(dòng)力學(xué)計(jì)算過(guò)程可分為下面三個(gè)步驟:首先將牛頓方程進(jìn)行離散化,求解各個(gè)粒子在t+?t的坐標(biāo);接下來(lái)是計(jì)算不同時(shí)刻(t+ ?t,t+2?t,...,t+l?t)的某個(gè)物理量A(r1(t+?t),r2(t+?t),...,rN(t+?t));最后對(duì)各個(gè)物理量求平均值:

經(jīng)典分子動(dòng)力學(xué)采用的是經(jīng)驗(yàn)勢(shì)或者半經(jīng)驗(yàn)勢(shì),相互作用勢(shì)不需要另外求解,因此計(jì)算速度較快,能夠計(jì)算的粒子數(shù)量也較大,可以達(dá)到千萬(wàn)個(gè)原子量級(jí).常用的經(jīng)驗(yàn)勢(shì)函數(shù)有Lennard-Jones勢(shì)函數(shù)[55],Morse勢(shì)函數(shù)[56],Born和Mayer[57]提出的描述離子晶體中離子之間相互作用的Born-Mayer勢(shì).至今為止,已經(jīng)有各種各樣的勢(shì)函數(shù)形式被提出,這些勢(shì)函數(shù)只考慮兩個(gè)原子之間的相互作用,僅和兩個(gè)原子位置坐標(biāo)相關(guān),在許多無(wú)機(jī)化合物中可以有效使用,但對(duì)于金屬以及共價(jià)鍵結(jié)合的有機(jī)分子,不能給出精確的結(jié)果.在20世紀(jì)80年代初,各種考慮了多體相互作用的新勢(shì)函數(shù)被相繼提出.Daw和Baskes[58]基于密度泛函理論提出了嵌入原子模型(EAM);Finnis和Sinclair[59]根據(jù)密度函數(shù)二次矩理論提出了形式上與EAM基本相同的經(jīng)驗(yàn)F-S模型.Abell[60]在1985年提出了一種Morse形式的多原子間相互作用勢(shì)模型,比傳統(tǒng)的二體勢(shì)和三體勢(shì)模型更加精確.

傳統(tǒng)勢(shì)函數(shù)對(duì)處于常溫常壓下的物質(zhì)可以給出很好的描述,而對(duì)處于高溫等離子體狀態(tài)時(shí)的物質(zhì)則存在較大誤差.一些描述電子結(jié)構(gòu)的半經(jīng)驗(yàn)方法:比如TF模型、AA模型、OF模型等能夠在高溫下較精確地描述離子運(yùn)動(dòng)的勢(shì)場(chǎng).結(jié)合這些半經(jīng)驗(yàn)方法和分子動(dòng)力學(xué)方法便可得到描述高溫稠密物質(zhì)性質(zhì)的模型:比如TFMD模型[38,41]、AAMD模型[46,47]以及OFMD模型[44,45].但由于TF模型、AA模型以及OF模型的假設(shè)都是基于高溫的情況,因此這些模型只能在高溫下描述電子結(jié)構(gòu),在低溫情況下會(huì)帶來(lái)一定的誤差.

5.2 第一性原理分子動(dòng)力學(xué)方法

第一性原理分子動(dòng)力學(xué)方法是基于密度泛函理論計(jì)算原子間相互作用勢(shì)的一種量子分子動(dòng)力學(xué)方法.1985年,Car和Parrinello[14]把密度泛函理論應(yīng)用到分子動(dòng)力學(xué)中,提出了第一原理分子動(dòng)力學(xué)方法(ab initiomolecular dynamics).該方法極大地推動(dòng)了計(jì)算機(jī)模擬低溫(T≤10 eV)下物態(tài)性質(zhì)的發(fā)展,已成為計(jì)算機(jī)模擬最先進(jìn)和最重要的方法之一.由于不需要對(duì)離子間勢(shì)函數(shù)做任何經(jīng)驗(yàn)假設(shè),而是直接采用密度泛函理論求解離子間的相互作用勢(shì),因此該方法具有非常強(qiáng)的普適性,并得到了非常廣泛的應(yīng)用[61?67].

在第一性原理分子動(dòng)力學(xué)中,原子核間的相互作用勢(shì)U(r)由密度泛函理論進(jìn)行計(jì)算,整個(gè)分子動(dòng)力學(xué)過(guò)程可分為三步:第一步,對(duì)于給定的原子核坐標(biāo)ri,求解Kohn-Sham方程[68]得到粒子數(shù)密度分布,進(jìn)而得到相互作用勢(shì);第二步,計(jì)算每個(gè)原子核受到的力;第三步,求解牛頓運(yùn)動(dòng)方程

密度泛函理論基于絕熱近似,在求解薛定諤方程時(shí),電子和原子核的坐標(biāo)同時(shí)出現(xiàn)在電子和原子核相互作用項(xiàng)中,直接求解難度非常大.由于原子核質(zhì)量遠(yuǎn)大于電子質(zhì)量,原子核的運(yùn)動(dòng)速度相對(duì)電子的運(yùn)動(dòng)速度非常小.在考慮電子的運(yùn)動(dòng)時(shí),可以認(rèn)為電子絕熱于原子核的運(yùn)動(dòng);在考慮原子核運(yùn)動(dòng)時(shí),電子可以及時(shí)調(diào)整空間分布來(lái)適應(yīng)原子核的運(yùn)動(dòng),這就是絕熱近似或稱Born-Oppenheimer近似[69].通過(guò)絕熱近似,求解多粒子體系薛定諤方程可以簡(jiǎn)化為分別求解電子運(yùn)動(dòng)和原子核運(yùn)動(dòng)的問(wèn)題.密度泛函理論的基本思想是將原子、分子和固體的基態(tài)物理性質(zhì)用粒子密度函數(shù)來(lái)描述.密度泛函理論的基礎(chǔ)是Hohenberg和Kohn關(guān)于非均勻電子氣體理論的Hohenberg-Kohn定理[70].粒子數(shù)密度函數(shù)是確定多粒子系統(tǒng)基態(tài)物理性質(zhì)的基本變量,要確定系統(tǒng)基態(tài)物理性質(zhì)需要三個(gè)確定量:粒子數(shù)密度函數(shù)、動(dòng)能泛函、交換關(guān)聯(lián)能泛函.粒子數(shù)密度函數(shù)和動(dòng)能泛函可以采用Kohn-Sham方程進(jìn)行求解;而交換關(guān)聯(lián)能泛函一般采用局域密度近似[71]和廣義梯度近似[72]得到.

原則上,量子分子動(dòng)力學(xué)方法可以描述任意條件下物質(zhì)的物態(tài)性質(zhì),但第一性原理分子動(dòng)力學(xué)的應(yīng)用范圍一直局限于溫度較低的區(qū)域(10eV以下),主要有以下原因:第一是在高溫條件下需要非常大的計(jì)算量;第二是大部分第一性原理程序?yàn)榱藴p少計(jì)算量都是基于贗勢(shì)實(shí)現(xiàn)的,而在高溫高壓下對(duì)贗勢(shì)的標(biāo)準(zhǔn)和構(gòu)造要求較高;第三是高溫狀態(tài)難以達(dá)到熱平衡,計(jì)算過(guò)程中難以保證收斂性[73].當(dāng)溫度較高時(shí),許多電子成為非局域電子或自由電子,相對(duì)于離子運(yùn)動(dòng)的時(shí)間尺度,電子運(yùn)動(dòng)的時(shí)間尺度非常小,因此在每一分子動(dòng)力學(xué)時(shí)間步中電子與離子發(fā)生了多次碰撞.在溫度升高時(shí),自由電子增多,電子的運(yùn)動(dòng)速度增快,與離子碰撞次數(shù)也增多,這時(shí)碰撞引起離子位置的漲落是不可忽略的,Born-Oppenheimer近似也不再成立.傳統(tǒng)的第一性原理分子動(dòng)力學(xué)方法基于Born-Oppenheimer近似,并沒(méi)有考慮離子位置的漲落,因此在高溫下的計(jì)算不可避免地帶來(lái)一些問(wèn)題.Dai等[48?50]對(duì)電子與離子在高溫下的碰撞作用采用類似布朗運(yùn)動(dòng)的思想來(lái)處理,引入離子-電子碰撞項(xiàng),利用Langevin方程描述離子運(yùn)動(dòng),取得了很好的效果,并將第一性原理分子動(dòng)力學(xué)方法拓展到了理想等離子體、高能量密度物理區(qū)域,填補(bǔ)了低溫凝聚態(tài)與高溫稠密狀態(tài)之間的空白.

第一性原理分子動(dòng)力學(xué)方法和其改進(jìn)方法雖然能夠在很大密度以及溫度范圍內(nèi)精確地得到物質(zhì)的狀態(tài)參數(shù),但由于其計(jì)算的復(fù)雜性以及所耗計(jì)算機(jī)時(shí)巨大,在精確度要求不是特別高的流體模擬中不如一些經(jīng)驗(yàn)或者半經(jīng)驗(yàn)的模型簡(jiǎn)單方便.一些全局狀態(tài)方程就是在不同溫度密度區(qū)域結(jié)合各類實(shí)驗(yàn)結(jié)果、理論結(jié)果再采用插值方法獲得的,比如More等[74]所提出的QEOS(quotidian equation of state),以及美國(guó)洛斯阿拉莫斯國(guó)家實(shí)驗(yàn)室的SESAME數(shù)據(jù)庫(kù)[75].

5.3 全局狀態(tài)方程

高能量密度條件下的流體動(dòng)力學(xué)過(guò)程涉及非常寬泛的物態(tài)范圍,比如慣性約束聚變中所涉及的溫度范圍為10?3—104eV、密度范圍為10?3—102g/cm3,因此在流體動(dòng)力學(xué)模擬中需要全局狀態(tài)方程.全局狀態(tài)方程就是在非常大的密度和溫度范圍內(nèi)描述壓強(qiáng)P(ρ,T)和內(nèi)能E(ρ,T).盡管量子分子動(dòng)力學(xué)方法原則上可以對(duì)任意狀態(tài)的物質(zhì)提供統(tǒng)一的描述,但實(shí)際上,目前還沒(méi)有哪個(gè)模型能夠覆蓋所有區(qū)域和所有物質(zhì),從而不能滿足流體動(dòng)力學(xué)模擬的需求.在不同區(qū)域采用不同的近似模型,并在這些區(qū)域間進(jìn)行光滑的插值,從而構(gòu)筑全局狀態(tài)方程模型是目前所采用的主要方法.目前已有一些大范圍的狀態(tài)方程(EOS)模型,比如Bushman和Fortov[76],Godval和Sikka[77],Basko[78]所提出的各種模型以及More等[74]所提出的QEOS模型.還有一些是以表格形式給出的,比如SESAME數(shù)據(jù)庫(kù)[75].SESAME數(shù)據(jù)庫(kù)的狀態(tài)方程數(shù)據(jù)是在不同溫度密度區(qū)域結(jié)合各類實(shí)驗(yàn)、理論結(jié)果采用插值方法獲得的.SESAME數(shù)據(jù)庫(kù)和QEOS模型使用比較廣泛,它們也是很多文獻(xiàn)進(jìn)行狀態(tài)方程數(shù)據(jù)分析和比對(duì)的重要參考.

QEOS把自由能分為電子項(xiàng)、離子項(xiàng)以及半經(jīng)驗(yàn)的束縛修正項(xiàng),表示如下:

對(duì)于電子項(xiàng)的貢獻(xiàn)Fe,可采用TF模型進(jìn)行計(jì)算;對(duì)于離子項(xiàng)的貢獻(xiàn)Fi,在固體中把熱離子當(dāng)作聲子氣體處理,采用玻色-愛因斯坦統(tǒng)計(jì)進(jìn)行計(jì)算,流體中則采用Cowan模型處理.半經(jīng)驗(yàn)的束縛項(xiàng)自由能不依賴于溫度,可表示為

式中E0和b是兩個(gè)經(jīng)驗(yàn)參數(shù),它們滿足以下兩個(gè)關(guān)系式:第一個(gè)是在常溫和固體密度條件下,總壓強(qiáng)為0;第二個(gè)是需要知道物質(zhì)的體積模量以及物質(zhì)的初始密度,通過(guò)壓強(qiáng)和自由能的熱力學(xué)關(guān)系式再聯(lián)立(51)式可以得到壓強(qiáng)和兩個(gè)經(jīng)驗(yàn)參數(shù)的關(guān)系,再通過(guò)B=ρ0(?P/?ρ)ρ0在已知體積模量B以及初始密度ρ0的條件下可以得到另一個(gè)關(guān)于E0和b的關(guān)系式.利用這兩個(gè)關(guān)系式可求得兩個(gè)經(jīng)驗(yàn)參數(shù)E0和b的值.

QEOS模型在國(guó)際上得到了比較廣泛的使用,很多輻射流體力學(xué)程序使用QEOS模型的計(jì)算結(jié)果作為狀態(tài)方程參數(shù),QEOS有可以免費(fèi)獲得的計(jì)算程序[79].Atzeni和Meyer-ter-Vehn[80]將鋁的Hugoniot曲線的理論計(jì)算結(jié)果與實(shí)驗(yàn)結(jié)果進(jìn)行了比較,如圖6所示.可以看出,QEOS結(jié)果和SESAME結(jié)果在104GPa以下與實(shí)驗(yàn)所得到的數(shù)據(jù)[81?84]比較接近.在104GPa以上時(shí),兩者之間有一定的差別,但總體上和實(shí)驗(yàn)數(shù)據(jù)符合較好,這說(shuō)明QEOS和SESAME都能給出較精確的狀態(tài)方程參數(shù).

圖6 鋁的Hugoniot曲線[80],圖中給出了QEOS結(jié)果[74],SESAME結(jié)果[75]以及一些實(shí)驗(yàn)結(jié)果[81?84]Fig.6.Hugoniot curves of aluminum[80],the b lack solid line and red dashed line are the resu lts of QEOS and SESAME,respectively.The dots in figure are experimental resu lts[81?84].

圖7[85]給出了溫度為5和30 eV時(shí)不同密度下Al等離子體的狀態(tài)方程.圖中比較了TF模型[86],QEOS模型[74],SESAME數(shù)據(jù)庫(kù)[75],NPA模型[87]以及AAMD模型[88]結(jié)果.NPA(neutral pseudoatom)模型[87]采用密度泛函理論處理電子結(jié)構(gòu),采用HNC近似描述離子間相互作用.從圖7可以看出,30 eV時(shí)所有模型結(jié)果的符合度比5 eV時(shí)好.溫度較高時(shí),一些基于高溫條件的半經(jīng)驗(yàn)?zāi)Ｐ?比如TF模型以及AAMD模型,也可以得到比較精確的結(jié)果.從圖7還可以看出,在5 eV時(shí)TF模型以及AAMD模型所得到的結(jié)果在較高密度條件下明顯要比NPA模型結(jié)果、QMD結(jié)果以及SESAME數(shù)據(jù)庫(kù)結(jié)果高.

圖7 溫度為5和30 eV時(shí)不同密度下鋁的狀態(tài)方程[85],圖中分別是TF模型[86]、QEOS[74]、SESAME[75]、NPA模型[87]以及AAMD模型[88]的計(jì)算結(jié)果比較Fig.7.The EOS of aluminiumunder the temperatu re of 5 and 30 eV[85],the solid lineswith d iff erentmarks are the resu lts of TF[86],QEOS[74],SESAME[75],NPA[87]and AAMD[88]models.

圖8[89]給出了TFMD模型[90]、 OFMD模型[73]、AAMD模型[88]、QMD模型、QLMD模型[50]以及SESAME數(shù)據(jù)庫(kù)[75]在不同溫度下鐵的Hugoniot曲線.由于高溫下QMD模型不容易收斂,因此沒(méi)有給出高溫下QMD的結(jié)果.從圖8可以看出,在高溫段(100 eV以上),TFMD,AAMD,OFMD以及QLMD四個(gè)模型和SESAME數(shù)據(jù)庫(kù)符合得非常好,這說(shuō)明四個(gè)MD模型在高溫下都能較精確地給出鐵的狀態(tài)方程,但在低溫下,TFMD模型和OFMD模型結(jié)果明顯要比其他幾個(gè)模型要高.相對(duì)前兩個(gè)模型,AAMD模型能在更低的溫度下給出精確結(jié)果.這三種半經(jīng)驗(yàn)?zāi)Ｐ徒Y(jié)果在低溫下比量子分子動(dòng)力學(xué)結(jié)果要高,這說(shuō)明了低溫下使用半經(jīng)驗(yàn)?zāi)Ｐ偷木窒扌?溫度較低時(shí),電子交換、關(guān)聯(lián)作用需要考慮[91],而TFMD方法在計(jì)算電子壓強(qiáng)時(shí)沒(méi)有考慮交換、關(guān)聯(lián)作用所產(chǎn)生的負(fù)壓強(qiáng)對(duì)總壓強(qiáng)的貢獻(xiàn),所以TFMD在低溫區(qū)給出的壓強(qiáng)偏高.AA模型雖然考慮了電子的殼層結(jié)構(gòu),但也只是簡(jiǎn)單的結(jié)構(gòu)近似,因此低溫時(shí)AAMD的精確度也會(huì)下降,這時(shí)只有基于密度泛函理論的量子分子動(dòng)力學(xué)模型才能給出較精確的結(jié)果.在0.1 eV的低溫下,電子的自旋效應(yīng)不能被忽視,考慮了量子自旋效應(yīng)的QLMD模型結(jié)果比未考慮量子自旋效應(yīng)的結(jié)果要稍微高一些,與SESAME數(shù)據(jù)庫(kù)結(jié)果更接近.

圖8 不同溫度下鐵的Hugoniot曲線[89],圖中分別給出了考慮了交換作用的TFMD模型[90]、AAMD模型[88]、OFMD模型[73]、QMD模型以及QLMD[50]結(jié)果和SESAME數(shù)據(jù)庫(kù)[75]結(jié)果Fig.8.The Hugoniot cu rve of iron at diff erent temperatu res[89],the resu lts of TFMD[90],AAMD[88],OFMD[73],QMD,QLMD[50]and SESAME[75]models are presented,respectively.

圖9[89]針對(duì)鐵的Hugoniot曲線給出了TFD模型[90]、VAAQP模型[92](variational-averageatom-in-quantum-plasmas,基于平均原子模型近似)、QEOS模型、QLMD模型[50]以及SESAME數(shù)據(jù)庫(kù)和各種實(shí)驗(yàn)結(jié)果[81,93,94]的比對(duì).由于沖擊壓縮實(shí)驗(yàn)中一些物理量的測(cè)量比較困難,實(shí)驗(yàn)結(jié)果有較大的不確定性,因此各類實(shí)驗(yàn)點(diǎn)的數(shù)據(jù)比較分散.從圖9可以看出,不同的理論模型給出的結(jié)果也存在一定差異,尤其壓強(qiáng)在10 Mbar以下、固體密度在2.5倍以下的溫稠密物質(zhì)區(qū)域,差別非常明顯.在此區(qū)域,TFD模型以及VAAQP模型結(jié)果偏大.QEOS模型、QLMD模型以及SESAME數(shù)據(jù)庫(kù)結(jié)果與俄羅斯的Rusbank數(shù)據(jù)庫(kù)[93]的結(jié)果十分接近.在高壓高密度情況下,各種模型雖然有一定的差異,但總體而言與實(shí)驗(yàn)結(jié)果符合較好.這說(shuō)明在高溫情況下,各種半經(jīng)驗(yàn)的統(tǒng)計(jì)模型也能給出比較精確的結(jié)果.

圖9 鐵的Hugoniot線的實(shí)驗(yàn)和理論結(jié)果的對(duì)比[89]Fig.9.The Hugoniot cu rve of iron[89],the lines are the resu lts of various theoreticalmodels and the dots are the experimental results.

6 結(jié) 語(yǔ)

本文根據(jù)離子耦合系數(shù)以及電子簡(jiǎn)并度參數(shù)在溫度密度平面內(nèi)把等離子體狀態(tài)劃分為四個(gè)區(qū)域,對(duì)四個(gè)區(qū)域內(nèi)等離子體的狀態(tài)方程及其熱力學(xué)性質(zhì)進(jìn)行了簡(jiǎn)要總結(jié).在完全電離等離子體區(qū)域,溫度非常高,原子完全電離,粒子之間相互作用非常小,電子和離子都可以采用理想氣體狀態(tài)方程描述.在部分電離弱耦合等離子體區(qū)域,部分電子電離,粒子之間相互作用相對(duì)較小,狀態(tài)方程可以采用Saha方程求解,考慮粒子之間長(zhǎng)程Coulomb相互作用后可采用Debye-Hückel方程得到.在強(qiáng)簡(jiǎn)并等離子體區(qū)域,原子核之間的相互作用由于周圍電子的屏蔽可以忽略,原子核可采用理想氣體狀態(tài)方程描述,電子采用TF模型可以得到較精確的結(jié)果,一些改進(jìn)型的TF模型可以把TF模型拓寬到更低溫度和更低密度.在強(qiáng)耦合等離子體區(qū)域,溫度不是很高,密度在物質(zhì)固體密度附近,此時(shí)物質(zhì)處于溫稠或者熱稠等離子體區(qū)域.當(dāng)?shù)入x子體處于前面三個(gè)區(qū)域時(shí),都有較準(zhǔn)確的方法描述等離子體的狀態(tài)方程,但等離子體處于溫稠以及熱稠區(qū)域時(shí),還沒(méi)有較完善的理論進(jìn)行統(tǒng)一的描述.量子分子動(dòng)力學(xué)雖然原則上能夠統(tǒng)一描述此區(qū)域物質(zhì)的性質(zhì),但由于計(jì)算量大、高溫不收斂等缺點(diǎn),其使用也受到一定的限制.一些改進(jìn)的量子分子動(dòng)力學(xué)方法比如QLMD,雖然明顯拓寬了量子分子動(dòng)力學(xué)方法的使用范圍,但其使用條件以及準(zhǔn)確性仍需要進(jìn)一步驗(yàn)證確認(rèn),而且其計(jì)算量過(guò)大,需要經(jīng)過(guò)大量的計(jì)算和累計(jì)才能得到較密集的狀態(tài)方程插值表,以滿足流體動(dòng)力學(xué)過(guò)程模擬的需要.一些半經(jīng)驗(yàn)分子動(dòng)力學(xué)方法,比如TFMD,AAMD,OFMD等,計(jì)算量雖小,但在低溫下的計(jì)算結(jié)果誤差較大.可見,發(fā)展描述溫稠和熱稠區(qū)域的先進(jìn)理論模型十分必要.另一方面，流體動(dòng)力學(xué)模擬所使用的狀態(tài)方程需要兼顧計(jì)算量和精度.狀態(tài)方程計(jì)算方法過(guò)于復(fù)雜,計(jì)算量太大,滿足不了流體動(dòng)力學(xué)模擬的需要;而狀態(tài)方程參數(shù)精確度太低,也會(huì)對(duì)流體動(dòng)力學(xué)模擬結(jié)果的可靠性產(chǎn)生較大影響.因此,采用不同的簡(jiǎn)單模型在不同密度和溫度區(qū)域進(jìn)行插值,不失為一種快速獲得等離子體狀態(tài)方程的好方法.

[1]Fortov V E 2016Extreme States of Matter(Berlin:Springer)p9

[2]Xu J A,MaoHK,Bell P M1986Science232 1404

[3]Ming L C,BassettW A1974Rev.Sci.Instrum.45 1115

[4]NellisW J 2006Rep.Prog.Phys.69 1479

[5]Q ian X S 2007Lectures on Physics(Shanghai:Shanghai Jiaotong University Press)(in Chinese)[錢學(xué)森2007物理力學(xué)講義(上海:上海交通大學(xué)出版社)]

[6]Fortov V E 2007Phys.Uspek.50 333

[7]Lind l J D,Amend t P,Berger R L,G lendinning S G,G lenzer S H,Haan SW,Kau ffman R L,Landen OL,Su ter L J 2004Phys.P lasmas11 339

[8]Saha MN,Srivastava BN 1965ATreatise on Heat(Allahabad:The India Press)

[9]Debye P,Hückel E 1923Physik Z24 185

[10]Thomas L H1927Proc.Cambridge Philos.Soc.23 542

[11]Fermi E 1928Z.Phys.48 73

[12]Xu X S,Zhang W X 1986Theory Guide of Practica l Equations of State(Beijing:Science Press)(in Chinese)[徐錫申,張萬(wàn)箱 1986實(shí)用物態(tài)方程理論導(dǎo)引(北京:科學(xué)出版社)]

[13]Feynman R P,Metropolis N,Teller E 1949Phys.Rev.75 1561

[14]Car R,ParrinelloM1985Phys.Rev.Lett.55 2471

[15]Letter R 1955Phys.Rev.99 1854

[16]D irac P AM1930Proc.Camb.Phil.Soc.26 376

[17]Cowan R D,Ashkin J 1957Phys.Rev.105 144

[18]Bloch F 1929Zeits.F Phys.57 545

[19]McCarthy S L 1965Lawrence Livermore Laboratory ReportUCRL-14364

[20]Kirzhnits D A1957Sov.Phys.JETP5 64

[21]Tang W H,Zhang R Q 2008In troduction tothe Theory and Compution of Equations of State(Beijing:Higher Education Press)(in Chinese)[湯文輝,張若棋2008物態(tài)方程理論及計(jì)算概論(北京:高等教育出版社)]

[22]Bell AR 1980Rutherford and Appleton Laboratories ReportRL-80-091

[23]Zink JW 1968Phys.Rev.176 279

[24]Lee C M,Thorsos BI1978Phys.Rev.A17 2073

[25]Rozsnyai BF 1972Phys.Rev.A5 1137

[26]Rozsnyai BF 1982J.Quan t.Spectrosc.Radiat.Transf.27 211

[27]Rozsnyai BF,Lamou reux M1990J.Quan t.Spectrosc.Radiat.Transf.43 381

[28]Bauche-Arnou lt C,Bauche J,Klapisch 1978M.J.Opt.Soc.Am.68 1136

[29]Bauche-Arnou lt C,Bauche J,Klapisch M1979Phys.Rew.A20 2424

[30]Bar-ShalomA,Oreg J,Goldstein W H,Shvart D,Zigler A1989Phys.Rev.A40 3183

[31]Oreg J,Goldstein W H,Bar-ShalomA,Klapisch M1990J.Comput.Phys.91 460

[32]Iglesias C A,Rogers F J,W ilson BG 1987Astrophys.J.322 45

[33]Kurten KE,Ristig ML,Clark JW 1977Lett.Al NuovoCimen to20 313

[34]Kang HS,Ree F H1998Phys.Rev.E57 5988

[35]Rose D V,GenoniTC,W elch D R,C lark R E,Campbell R B,Meh lhorn TA,Flicker D G 2009Phys.P lasams16 102105

[36]D ror O,Nard i E,Rosen feld Y 1988Phys.Rev.A38 5801

[37]Fu rukawa H,N ishihara K1992Phys.Rev.A46 6596

[38]Zérah G,C lérouin J G,Pollock E L 1992Phys.Rev.Lett.69 446

[39]C lérouin J G,Pollock E L,Zérah G 1992Phys.Rev.A46 5130

[40]Flavien L,Gilles Z 2006Phys.Rev.E73 016403

[41]Danel J F,Kazand jian L,Zérah G 2006Phys.P lasmas13 092701

[42]Dharma-wardana MW C,Francois P 1982Phys.Rev.A26 2096

[43]Dharma-wardana MW C,MurilloMS 2008Phys.Rev.E77 026401

[44]Danel J F,Kazand jian L,Zérah G 2009Phys.Rev.E79 066408

[45]Horner D A,Lambert F,Kress J D,Collins L A2009Phys.Rev.B80 024305

[46]Hou Y,Jin F T,Yuan J M2006Phys.P lasmas13 093301

[47]Hou Y,Jin F T,Yuan J M2007J Phys.:Condens.Matter19 425204

[48]Dai J,Yuan J 2009Europhys.Lett.88 20001

[49]Dai J,Hou Y,Yuan J 2010Phys.Rev.Lett.104 245001

[50]Dai J,Hou Y,Yuan J 2010Astrophys.J.721 1158

[51]Militzer B2009Phys.Rev.B79 155105

[52]W unsch K,Vorberger J,Gericke D O2009Phys.Rev.E79 010201

[53]Gericke D O,W unsch K,G rinenkoA,Vorberger J 2010J.Phys.:Conf.Ser.220 012001

[54]Pey russe O,Mazevet S,Recou les V,Dorchies F,Harmand M,Levy A,Fuchs J,Mancic A,NakatsutsumiM,Renaudin P,Audebert P 2009CP,Atomic Processes in P lasmas1161 200

[55]Frenkel D,Smit B1996Understanding Molecular Simu lation(San D iego:Acedemic Press)

[56]Morse P M1929Phys.Rev.34 57

[57]Born M,Mayer J E 1931Z.Phys.75 1

[58]DawMS,Baskes MI1983Phys.Rev.Lett.50 1285

[59]Finnis MW,Sinclair J E 1984Philosophic Magazine A50 45

[60]Abell G C 1985Phys.Rev.B31 6184

[61]Collins L,Kwon I,Kress J,Trou llier N,Lynch D 1995Phys.Rev.E52 6202

[62]Desjarlais MP,Kress J D,Collins L A2002Phys.Rev.E66 025401

[63]Mazevet S,C lérouin J,Recou les V,Anglade P M,Zérah G 2005Phys.Rev.Lett.95 085002

[64]Mazevet S,Desjarlais MP,Collins L A,Kress D,Magee N H2005Phys.Rev.E71 016409

[65]C lérouin J,Noiret P 2008Phys.Rev.B78 224203

[66]Holst B,Redmer R 2008Phys.Rev.B77 184201

[67]Mazevet S,Zérah G 2008Phys.Rev.Lett.101 155001

[68]Kohn K,ShamL J 1965Phys.Rev.140 A1133

[69]Born M,Huang K1954Dynamical Theory of Crystal Lattice(Ox ford:Ox ford University Press)

[70]Hohenberg P,Kohn W 1964Phys.Rev.B136 864

[71]Hansen J P,McDonald IR 1976Theory of Simple Liquids(Lodon:Academic Press)

[72]PerdewJ P,Kieron B,Matthias E 1996Phys.Rev.Lett.77 3865

[73]Lambert F,Clérouin J,Zérah G 2006Phys.Rev.E73 016403

[74]More R M,W arren KH,Young D A,Zimmerman G B1988Phys.F luids31 3059

[75]Lyon S P,Johnson J D 1992SESAME:The Los Alamos National Laboratory Equation of State DatabaseReport No.LA-UR-92-3407

[76]Bushman AV,Fortov V E 1983Soviet Phys.Uspek.26 465

[77]Godval BK,Sikka S K1983Phys.Rep.102 121

[78]BaskoMM.Metallic 1985Soviet High Temperat.Phys.23 388

[79]KempAJ,Meyer-ter-Vehn J 1998Nucl.Instrum.Meth.Phys.Res.A415 674

[80]Atzeni S,Meyer-ter-Vehn J 2004The Physics of Inertia l Fusion(NewYork:Oxford University Press)

[81]Altshu ler L 1965Soviet Phys.Uspek.8 52

[82]Ragan C E 1982Phys.Rev.A25 3360

[83]V lad imirov AS,Voloshin N P,Nogin V N,Petrov tsev AV,SimonenkoV A1984JETP Lett.39 82

[84]Mitchell AC,Nellis W J 1981J.Appl.Phys.52 3363

[85]Hou Y 2009Ph.D.D issertation(Changsha:National University of Defence Technology)(in Chinese)[侯永2009博士學(xué)位論文(長(zhǎng)沙:國(guó)防科學(xué)技術(shù)大學(xué))]

[86]Letter R 1955Phys.Rev.99 1854

[87]Perrot F,Dharma-wardana MW C,Benage J 2002Phys.Rew.E65 046414

[88]Hou Y,Yuan J 2009Phys.Rev.E79 016402

[89]Dai J Y 2009Ph.D.D issertation(Changsha:National University of Defence Technology)(in Chinese)[戴佳鈺2010博士學(xué)位論文(長(zhǎng)沙:國(guó)防科學(xué)技術(shù)大學(xué))]

[90]Lambert F,C lérouin J,Mazevet S 2006Europhys.Lett.75 681

[91]Fromy P,Deu tsh C,Maynard G 1996Phys.P lasmas3 714

[92]Piron R,Blenski T2011Phys.Rev.E83 026403

[93]Shock W ave Database http://www.ihed.ras.ru/rusbank/gassim/[2003-07-13]

[94]Batani D,Morelli A,Tomasini M,Benuzzi-Mounaix A,Philippe F,Koenig M,Marchet B,Masclet I,Rabec M,Reverd in C,Caub le R,Celliers P,Collins G,Silva L D,Hall T,Moret M,Sacchi B,Baclet P,Cathala B2002Phys.Rev.Lett.88 235502

PACS:05.70.Ce,52.25.KnDOI:10.7498/aps.66.030505

Equations of state and thermodynamic properties of hot plasma

Tang Wen-Hui?Xu Bin-Bin Ran Xian-Wen Xu Zhi-Hong

(College of Science,National University of Defense Technology,Changsha 410073,China)(Received 17 October 2016;revised manuscript received 17 December 2016)

The equations of state(EOS)and the thermodynamics properties of plasma under high temperature are widely applied tothe fields of astrophysics,controllable fusion,weapon design and damage.In this paper wemain ly reviewthe theoreticalmodel and computing method of the EOS of hot plasma on diff erent density scales and temperature scales.For an ideal plasma,the interaction between ions can be ignored,the EOS is simple and the theories turn matured.Under the condition of extremely high temperature,ions are ionized completely and the EOSs of ions and electrons can be approximated by the EOS of ideal gas.W hen the temperature is not very high and ions are just partly ionized,the EOS can be obtained by Saha model or its modified model.W hen atoms are strongly compressed,the EOS can be calculated by Thomas-Fermimodel or its modified model.For the non-ideal plasma,there is a strong coupling between ions.Nounified theoreticalmodel can completely describe the interaction between ions at arbitrary density and arbitrary temperature.In principle,the quantummolecular dynamics(QMD)can accurately describe the EOS of plasma in large density range and large temperature range.However,due tothe enormous computation and the diffi culty in converging,it is diffi cult toapply QMD tothe plasma under high temperature.W ith simple computing method and small computation,classicalmolecu lar dynamics using semi-empirical potential can calcu late the EOS accurately at high temperature.However,it will produce great error at lower temperature.It is a simple and eff ective way toobtain a global EOS by using diff erent theoreticalmodels in diff erent density range and diff erent temperature range and by interpolating in the vacant density range and vacant temperature range.

hot plasma,equations of state,thermodynamics properties

10.7498/aps.66.030505

?通信作者.E-mail:wenhuitang@163.com

?Corresponding author.E-mail:wenhuitang@163.com