基于 Gogny相互作用的非對稱核物質(zhì)同位旋的相關(guān)性質(zhì)

陳 融, 楊艷芳, 陳列文

(1.上海大學(xué)理學(xué)院,上海 200444;2.上海交通大學(xué) 物理系,上海 200240)

基于 Gogny相互作用的非對稱核物質(zhì)同位旋的相關(guān)性質(zhì)

陳 融1,2, 楊艷芳1, 陳列文2

(1.上海大學(xué)理學(xué)院,上海 200444;2.上海交通大學(xué) 物理系,上海 200240)

基于有限程的 Gogny相互作用,系統(tǒng)研究同位旋非對稱核物質(zhì)的狀態(tài)方程 (equation of state,EOS)及其單粒子勢和它們的同位旋相關(guān)性.通過對比各類模型以及實驗值,說明現(xiàn)有的 D1,D1S和 D1N這 3組 Gogny參數(shù)的不足之處,指出探尋新的 Gogny核力參數(shù)的必要性.

Gogny相互作用;非對稱核物質(zhì);狀態(tài)方程;同位旋

本工作研究的是零溫均勻無窮大非對稱核物質(zhì),該系統(tǒng)僅由中子和質(zhì)子組成,且只考慮強相互作用.雖然在實際情況下,人們研究的往往是有限大核物質(zhì),但本研究依然可以根據(jù)均勻無窮大零溫系統(tǒng)這個再考慮模型推導(dǎo)出的一系列物理量,結(jié)合表面能、庫侖的相互作用等進(jìn)行修正,進(jìn)而給出和實驗相吻合的結(jié)果.同時,很多真實的物理系統(tǒng)可以很好地看成是無窮大核物質(zhì),如中子星的內(nèi)部、重原子核(如208Pb)內(nèi)部等.

1 Gogny相互作用

Gogny兩體相互作用的一般形式為[2]

式中,r=r1-r2,Pσ為自旋交換算符,Pτ為同位旋交換算符.從式 (1)可以看到,Gogny相互作用由 11個含量綱參數(shù) (W1,W2,B1,B2,H1,H2,M1,M2,μ1,μ2,t0)和 2個不含量綱的參數(shù) (x0,α)共 13個參數(shù)構(gòu)成 (此式省略了自旋軌道耦合項,因為該項不適用于均勻無窮大核物質(zhì)).

文獻(xiàn)[2-9]均未明確給出基于 Gogny相互作用的單粒子勢和每核子能量的具體表達(dá)式.本研究給出了基于 Gogny相互作用的零溫均勻無窮大非對稱核物質(zhì)的單粒子勢和每核子能量的表達(dá)式,這對研究相關(guān)的物理量具有重要意義.同時,本工作也將系統(tǒng)地研究同位旋非對稱核物質(zhì)的狀態(tài)方程和單粒子勢及其同位旋相關(guān)性.

2 計算與分析

本研究根據(jù)已有的 3組 Gogny相互作用參數(shù)D1[2],D1S[10],D1N[11],并結(jié)合其他核力模型和實驗數(shù)據(jù),分析了非對稱核物質(zhì)中每核子能量隨密度的變化關(guān)系、單粒子勢隨單核子能量的變化關(guān)系、對稱勢隨核子動能的變化關(guān)系、等壓壓縮系數(shù)隨同位旋不對稱度的變化關(guān)系、有效質(zhì)量的同位旋劈裂以及對稱能隨密度的變化關(guān)系等.此外,還計算了飽和密度處對稱核物質(zhì)性質(zhì)以及對稱能的高階特征參數(shù),討論了等壓不可壓縮系數(shù)同位旋依賴的拋物線近似.本研究中所有的物理量均采用自然單位制.

2.1 每核子能量與對稱能

每核子能量 E(ρ,δ)是在數(shù)密度為ρ=ρn+ρp,同位旋不對稱度為δ=(ρn-ρp)/ρ的核物質(zhì)中,平均每個核子的總能量.非對稱核物質(zhì)的狀態(tài)方程(equation of state,EOS)就是每核子能量 E(ρ,δ)的表達(dá)式.基于 Gogny相互作用的 EOS為

式中 ,fτ=2h-3Θ [pf(τ) -p]為零溫時質(zhì)子 (或中子 )的 Fermi分布函數(shù) (Θ [pf(τ) -p]為階躍函數(shù) ),RT=t0αρα+1[3-(2x0+1)δ2]/8為核物質(zhì)的重排項,pFn,pFp分別為中子和質(zhì)子的 Fermi動量.

圖 1給出了 3組參數(shù) (D1,D1S和 D1N)下,不同δ時的每核子能量隨密度的變化關(guān)系.由圖 1可見,對于 3組 Gogny相互作用參數(shù),E(ρ,δ)隨同位旋不對稱度δ的提高而增加,這說明 E(ρ,δ)是同位旋相關(guān)的,當(dāng)δ=0時,即為對稱核物質(zhì).研究發(fā)現(xiàn),對任意一組參數(shù),δ=0時的 E(ρ,δ)均在 ρ0=0.16 fm-3左右達(dá)到極小值,約為 -16 MeV(即ρ0=0.16 fm-3稱為飽和密度),和實驗值相符.當(dāng)δ較大時,E(ρ,δ)在密度很小情況下會出現(xiàn)先增后減再增的現(xiàn)象,甚至當(dāng)δ大于某個臨界值δc后,E(ρ,δ)會隨密度單調(diào)上升,此時就不存在飽和密度.計算表明,D1參數(shù)的δc約為 0.79,D1S的δc約為 0.78,D1N的δc約為 0.83.

圖 1 不同δ下 3組參數(shù)每核子能量隨密度的變化Fig.1 Relationship between density and energy per nucleon w ith d ifferentδ

在核天文學(xué)中,對稱能是與中子星中質(zhì)子的比例、中子星的演化、中子星的半徑和密度等密切相關(guān)的重要物理量,因此,有關(guān)對稱能的研究一直是人們十分感興趣的問題[1].

將每核子能量 E(ρ,δ)泰勒展開[12],其中的就是非對稱核物質(zhì)的對稱能Esym(ρ).可以看出,對稱能為核子數(shù)密度ρ的函數(shù),與δ無關(guān).圖 2給出了 3組參數(shù)的對稱能隨密度的變化關(guān)系.可以看出,3組參數(shù)的對稱能隨密度都呈先增后減的趨勢.在ρ0=0.16 fm-3附近,3組參數(shù)的對稱能和實驗值 30 MeV吻合得很好.目前,關(guān)于對稱能在密度較大區(qū)域究竟是增還是減,并不十分清楚.

圖 2 3組參數(shù)下對稱能隨密度的變化Fig.2 Relationship between density and symmetry energy under three param eter sets of Gogny interaction

2.2 單粒子勢(光學(xué)勢)隨核子能量的變化

單粒子勢 U(ρ,δ,p,τ)=Uτ表示動量為 p的核子在核子數(shù)密度為ρ,同位旋不對稱度為δ的核物質(zhì)中感受到的所有來自其他核子的作用.基于Gogny相互作用的單粒子勢可表示為

式中,τ為同位旋第三分量的大小,中子τn=-1/2,質(zhì)子 τp=-1/2,kfτ=(3π2ρτ)1/3為質(zhì)子 (或中子 )的Fermi波矢 ,pfτ= ?kfτ為質(zhì)子 (或中子 )的 Fermi動量.

圖 3給出了單粒子勢 (光學(xué)勢[13-15])隨單核子能量 E-M的變化關(guān)系,其中單核子能量 E-M=可見,3組Gogny參數(shù)在能量較大時,單粒子勢始終小于 0 MeV,這與實驗值 30 MeV不符 (見圖 3中 Hama[16]的 Fit 1和 Fit 2).這說明現(xiàn)有的 3組 Gogny相互作用參數(shù)還有待改進(jìn).

圖 3 單粒子勢(光學(xué)勢)隨核子能量的變化Fig.3 Relationship between single nucleon energy and single par ticle poten tial(optical poten tial)

2.3 對稱勢

研究非對稱核物質(zhì)中的對稱勢,對理解核物理和天體物理中的許多重要物理問題具有重要意義[1].零溫時的對稱勢 Usym定義為

研究表明,Usym(ρ,p)的大小對δ不敏感.為方便起見,取δ=0.5,則對稱勢簡化為

由于中子和質(zhì)子的質(zhì)量十分接近,因此,可以看成是質(zhì)量 m0=939 MeV/c2的核子.圖 4給出了對稱勢隨核子動能 Ekin的變化關(guān)系.由圖可知,D1給出的結(jié)果與唯象的動量相關(guān)相互作用 (momentumdependent interaction, MD I) 模 型[13-14,17-18]、DBHF(Dirac-Brueckner-Hartree-Fock)理論、非相對論BHF(Brueckner-Hartree-Fock)理論以及相對論沖擊近似 (relativistic impulse app roximation,RIA)等微觀理論給出的結(jié)果類似[15],這與目前公認(rèn)的結(jié)果比較接近:即隨著核子動能的增大,中子的單粒子勢會逐漸變小,最終會小于質(zhì)子的單粒子勢,使得對稱勢為負(fù)[15].相比較而言,D1S和 D1N這兩組參數(shù)雖然在核子動能很低的情況下給出的結(jié)果與微觀理論的結(jié)果比較吻合,但當(dāng)動能繼續(xù)增大時,特別是動能大于 100 MeV后,對稱勢將大大偏離實驗值.

圖 4 對稱勢隨核子動能的變化Fig.4 Relationship between nucleon k inetic energy and symmetry potential

2.4 有效質(zhì)量

在核物理中,另一表征單粒子勢性質(zhì)的物理量為核子的有效質(zhì)量m＊τ.有效質(zhì)量反應(yīng)了核子單粒子勢的動量相關(guān)性[1],其定義式為

在 Gogny相互作用中,非對稱核物質(zhì)中質(zhì)子和中子的有效質(zhì)量隨動量 p的變化是不同的,這被稱為有效質(zhì)量的同位旋劈裂.

圖 5給出了δ=0.5時 3組參數(shù)飽和密度時的有效質(zhì)量隨動量的同位旋劈裂.圖 6給出了δ=0.5時,在各密度的 Fermi動量處的有效質(zhì)量隨密度的同位旋劈裂.由圖 5和圖 6可見,D1N給出的最新結(jié)果與 D1,D1S區(qū)別較大.在圖 5中,當(dāng)動量很大時,3組參數(shù)所給出的質(zhì)子和中子的有效質(zhì)量都趨于它們的靜質(zhì)量,而最新的D1N參數(shù)表現(xiàn)出的中子有效質(zhì)量始終小于質(zhì)子的有效質(zhì)量.在圖 6中,D1N給出的結(jié)果與 D1,D1S區(qū)別更為明顯.隨密度增大,D1和D1S參數(shù)的同位旋劈裂較緩慢.而 D1N參數(shù)的同位旋劈裂則隨密度的增大變化得很快,并且中子的有效質(zhì)量幾乎總是小于質(zhì)子的有效質(zhì)量.但是在豐中子核物質(zhì)中,已有的研究認(rèn)為中子的有效質(zhì)量要比質(zhì)子大[12],這樣看來,D1N的結(jié)果不夠理想.

圖 5 有效質(zhì)量隨動量的同位旋劈裂Fig.5 Isosp in splitting of effectivemass w ith momentum

圖 6 Ferm i動量處有效質(zhì)量隨密度的同位旋劈裂Fig.6 Isosp in splitting of effectivemass w ith density at Ferm imomentum

2.5 飽和密度處的性質(zhì)

飽和密度 (每核子能量取極小值時的密度)處的核物質(zhì)的性質(zhì)一直是本研究比較關(guān)心的問題.這里主要研究不可壓縮系數(shù) (等壓壓縮系數(shù))在飽和密度處的性質(zhì).最近的研究表明,δ的高階項對于確定極其豐中子的中子星外殼內(nèi)緣分離出液態(tài)核時的轉(zhuǎn)變密度起著非常重要的作用[12].

不可壓縮系數(shù) K為描述核物質(zhì)性質(zhì)的一個重要參量[12],可以用來描述系統(tǒng)的力學(xué)穩(wěn)定性等性質(zhì),其定義式為式中 ,為壓強.K(ρ,δ)在飽和密度 ρ0處可寫為 Ksat(δ),稱為等壓壓縮系數(shù) (isobaric incompressibility)[19].

為了研究等壓壓縮系數(shù)的性質(zhì),將對稱核物質(zhì)的每核子能量 E0(ρ)=E(ρ,δ=0)在 ρ0處展開到第四階,即式中,反映了在具有不對稱度δ的非對稱核物質(zhì)中,飽和密度處的不可壓縮系數(shù)相對于對稱核物質(zhì) (δ=0)飽和密度處的不可壓縮系數(shù)的偏移,或者說,反映了等壓壓縮系數(shù)Ksat(δ)對同位旋的依賴性.圖 7和圖 8分別給出了 3組參數(shù)下的等壓壓縮系數(shù)隨δ以及δ2的變化關(guān)系.

圖 7 3組參數(shù)等壓壓縮系數(shù)隨δ變化Fig.7 Isobar ic incompressib ility vs.δunder three param eter sets

由圖 7可見 ,對于 3組 Gogny參數(shù) ,Ksat(δ)隨δ的增大總是減小的,最終在臨界值δc處,Ksat(δ)仍然大于 0,這說明在飽和密度處,系統(tǒng)總是力學(xué)穩(wěn)定的.由圖 8可見,(Ksat-K0)隨δ2的變化曲線幾乎可以用 Ksat,2δ2這條直線來代替 ,這說明 Ksat(δ)很好地滿足了拋物線近似,也說明當(dāng)斜率參量L比較小時,Ksat(δ)的高階量對 Ksat(δ)影響不大.

圖 8 3組參數(shù)的等壓壓縮系數(shù)隨δ2的拋物線近似Fig.8 Parabolic approximation betweenδ2and isobar ic incompressibility under three param eter sets

最后,根據(jù)計算結(jié)果,本研究在表1中給出一些高階量和關(guān)鍵參數(shù)的具體值,其中ρ0為δ=0時的對稱核物質(zhì)的飽和密度,E0/A為對稱核物質(zhì)飽和密度處的每核子能量的值.可以看出,3組參數(shù)下的ρ0,E0/A,Esym(ρ0)的值與實驗值比較吻合.但是需要指出的是,這 3組參數(shù)下斜率參量 L的最大值均不超過 33.6 MeV,這和實驗值 L=60.0 MeV相比,還是偏小的.這也說明現(xiàn)有的 3組 Gogny相互作用參數(shù)還有待改進(jìn).

表1 3組參數(shù)下一些高階量和關(guān)鍵參數(shù)Table 1 Values under three sets of Gogny parameter s

3 結(jié) 束 語

本工作基于 Gogny相互作用,研究了零溫均勻無窮大非對稱核物質(zhì)的同位旋相關(guān)性.結(jié)果發(fā)現(xiàn),現(xiàn)有的 D1,D1S和 D1N 3組參數(shù)都與實驗值 (趨于30 MeV)差距比較大,不能很好地描述非對稱核物質(zhì)的單粒子勢;同時,D1S和 D1N也不能很好地描述對稱勢.另外,從研究的有效質(zhì)量同位旋劈裂的結(jié)果來看,D1N參數(shù)給出的結(jié)果和預(yù)期值相差較大.此外,還研究了等壓不可壓縮系數(shù) Ksat(δ)的同位旋隨不對稱度δ的變化關(guān)系,討論了 Ksat(δ)的同位旋依賴的拋物線近似.結(jié)果發(fā)現(xiàn),Ksat(δ)能很好地滿足拋物線近似,這說明 Ksat(δ)的高階量對于這 3組Gogny相互作用參數(shù)的影響不大.本研究還計算了飽和密度處的對稱核物質(zhì)的性質(zhì)及對稱能的高階特征參數(shù),列出了一些關(guān)鍵值.過小的L值也說明了尋找新的 Gogny相互作用參數(shù)的重要性.

從對比各類模型以及實驗值的綜合結(jié)果來看,現(xiàn)有的D1,D1S和D1N這 3組 Gogny相互作用參數(shù)還存在不足,所以給出新的更好的 Gogny相互作用參數(shù)是十分必要的.

[1] DECHECHARGé J,GOGNY D.Hartree-Fock-Bogolyubov calculationswith the D1 effective interaction on spherical nuclei[J].Physical Review C,1980,21(4):1568-1593.

[2] SKYRME T.The effective nuclear potential[J].Nuclear Physics,1959,9(4):615-634.

[3] DAS C B,GUPTA S D,GALE C,et al.Momentum dependence of symmetry potential in asymmetric nuclear matter for transport model calculations[J]. Physical Review C,2003,67(3):034611.

[4] ZHANG Y J, SU R K, SONG H Q, et al.Thermodynamical p roperties and Coulomb instabilities in hot nuclear systems with the Gogny interaction[J].Physical Review C,1996,54(3):1137-1143.

[5] áNGELES M,GARCíA P.Magnetization of a neutron p lasmawith skyrme and Gogny forces in the presence of a strongmagnetic field[J].Physical Review C,2008,77(6):065806.

[6] KUCHAREK H,RING P,SCHUCK P.Pairing properties of nuclear matter at finite temperature[J].Z Phys A:Atomic Nuclei,1989,334:119-124.

[7] 陳融.非對稱核物質(zhì)的狀態(tài)方程 [D].上海:上海大學(xué),2009.

[8] FELDMEIER H,L INDNER J.Field-dependent coup ling strength for scalar fields[J].Zeitschrift für Physik A Hadrons and Nuclei,1991,341(1):83-88.

[9] DANIELEW ICZ P. Determination of the mean-field momentum-dependence using ellip tic flow[J].Nuclear Physics A,2000,673:375-410.

[10] CHAPPERT F,GIROD M,HILARIRE S.Towards a new Gogny force parameterization: impact of the neutron matter equation of state[J].Physics Letters B,2008,668:420-424.

[11] SONG H Q,ZHENG GD,SU R K.Critical phenomena in nuclear matterwith Gogny interaction[J].J Phys G:Nucl Part Phys,1990,16:1861-1871.

[12] RING P,SCHUCK P.The nuclear many-body p roblem[M].New York:Sp ringer-Verlag,1980:147-177.

[13] LáSZLó P,CSERNA I G F,CHARLES G,et al.Nuclear equation of state with momentum-dependent interactions[J].Physical Review C,1992,46(2):736-747.

[14] CHEN L W,KO C M,L I B A.Determination of the stiffness of the nuclear symmetry energy from isospin diffusion[J].Phys Rev Lett,2005,94:032701.

[15] L IB A,CHEN L W,KO CM.Recentprogress and new challenges in isospin physics with heavy-ion reactions[J].Physics Reports,2008,464:113-281.

[16] COOPER E D,HAMA S,CLARK B C,et al.Global Dirac phenomenology for proton-nucleuselastic scattering[J].Physical Review C,1993,47(1):297-311.

[17] 徐駿.非對稱核物質(zhì)熱力學(xué)性質(zhì)、中子星結(jié)構(gòu)和蒙特卡羅算法相關(guān)問題的研究[D].上海:上海交通大學(xué),2008.

[18] WELKE G M,PRAKASH M,KUO T T S,et al.Azimuthal distributions in heavy ion collisions and the nuclear equation of state[J]. Physical Review C,1988,38(5):2101-2107.

[19] CHEN L W,CA I B J,KO C M,et al.Higher-order effects on the incompressibility of isospin asymmetric nuclear matter[J].Physical Review C,2009,80(1):014322.

Isosp in Dependent Proper ties of Asymmetr ic Nuclear M atter Using Gogny Interaction

CHEN Rong1,2, YANG Yan-fang1, CHEN Lie-wen2
(1.College of Sciences,ShanghaiUniversity,Shanghai200444,China;2.Department of Physics,Shanghai Jiao Tong University,Shanghai200240,China)

Using the finite-range Gogny interaction,the equation of state(EOS)and the single-particle potential as well as their isosp in dependence in isosp in asymmetric nuclear matter is studied.By comparison with othermodelsand empirical values,we show the limitation of the Gogny interactionsD1,D1S and D1N,and necessity of a new parameterization of Gogny interaction.

Gogny interaction;asymmetric nuclear matters;equation of state(EOS);isospin

O 571

A

1007-2861(2011)02-0164-06

10.3969/j.issn.1007-2861.2011.02.011

2009-09-23

國家自然科學(xué)基金資助項目(10575071,10675082);上海市青年科技啟明星計劃資助項目(06QA 14024),國家重點基礎(chǔ)研究發(fā)展計劃 (973計劃)資助項目 (2007CB815004);上海市重點學(xué)科建設(shè)資助項目(S30105)

楊艷芳 (1973～),女,副教授,博士,研究方向為物理學(xué).E-mail:yangyangzhaoyb@yahoo.com.cn

(編輯:劉志強)

自量子色動力學(xué) (quantum chromodynamics,QCD)被提出以來,人們對強相互作用有了進(jìn)一步的認(rèn)識.雖然QCD可以描述強相互作用,理論上也應(yīng)該可以推導(dǎo)出核力,但是在使用 QCD解決具體問題時,其復(fù)雜性令其難以直接求解.而核物理研究的一個極為重要的部分,就是找到真正的“核力”.為此,人們構(gòu)造了含經(jīng)驗參量的簡化模型來處理與核力相關(guān)的具體問題,如基于零程的 Skyrme相互作用的Hartree-Fock理論和基于介子交換模型的相對論平均場理論等.

1980年,Gogny等[1]提出了 Gogny相互作用,修正了 Skyrme作用的一些不足之處,并對 Gogny相互作用進(jìn)行了計算,給出了與實驗吻合得更好的結(jié)果.因此,本研究使用 Gogny相互作用.