趙先鋒

(滁州學(xué)院電子信息工程系，安徽 滁州 239000)

在包括σ,ω,ρ介子及σ介子自相互作用的相對(duì)論平均場(chǎng)理論(RMF)中，需要確定核子與介子的耦合常數(shù)。 近些年來(lái)，在相對(duì)論平均場(chǎng)理論的研究中，人們得到了多組介子與核子的耦合參數(shù)，如NL3、TM1、TM2、NLSH、GL85和GL97等[1-5]。這五個(gè)耦合常數(shù)中，其中的gρ是描述無(wú)限核物質(zhì)中質(zhì)子和中子的不對(duì)稱(chēng)性質(zhì)的，它聯(lián)系于飽和核物質(zhì)的密度ρ0、有效質(zhì)量m*，尤其是對(duì)稱(chēng)能量系數(shù)asym[5]。由于實(shí)驗(yàn)核物理的不斷發(fā)展，對(duì)于飽和核物質(zhì)對(duì)稱(chēng)能的研究乃是一個(gè)非?；钴S的領(lǐng)域。

對(duì)稱(chēng)能的研究可以通過(guò)多種實(shí)驗(yàn)方法進(jìn)行，其中最重要的一種是通過(guò)重離子碰撞獲得相關(guān)數(shù)據(jù)[6-15]。2004年，Tsang等[16]在美國(guó)密歇根州立大學(xué)的國(guó)家超導(dǎo)回旋加速器實(shí)驗(yàn)室(NSCL/MSU)進(jìn)行了112Sn和114Sn的重離子碰撞試驗(yàn)，給出了較新的對(duì)稱(chēng)能實(shí)驗(yàn)數(shù)據(jù)。

本文從相對(duì)論平均場(chǎng)理論出發(fā)，導(dǎo)出ρ介子與核子的耦合參數(shù)gρ的代數(shù)關(guān)系式； 并對(duì)相對(duì)論平均場(chǎng)理論計(jì)算中曾用到的幾組ρ介子耦合參數(shù)進(jìn)行分析；采用NSCL/MSU給出的核子最新實(shí)驗(yàn)數(shù)據(jù)，對(duì)ρ介子與核子的耦合參數(shù)gρ進(jìn)行計(jì)算并對(duì)結(jié)果進(jìn)行討論。

1 耦合常數(shù)gρ代數(shù)表達(dá)式的推導(dǎo)

核物質(zhì)的拉格朗日密度為[5]

(1)

式中，Ψ是核子的Dirac旋量，對(duì)應(yīng)的質(zhì)量為m。σ,ω,ρ分別是σ介子、ω介子和ρ介子的場(chǎng)算符，mσ,mω,mρ是這些介子的質(zhì)量.gσ,gω,gρ,g2,g3為核子耦合常數(shù)，它們可以利用飽和核物質(zhì)的性質(zhì)，如密度ρ0、束縛能B/A、壓縮系數(shù)K、有效質(zhì)量m*和對(duì)稱(chēng)能量系數(shù)asym定出。

對(duì)于核物質(zhì)，應(yīng)用相對(duì)論平均場(chǎng)近似后，得到由中子和質(zhì)子非對(duì)稱(chēng)性所貢獻(xiàn)的能量密度為[5]

(2)

式中k為核子的費(fèi)米動(dòng)量。

令核物質(zhì)核子數(shù)密度為ρ，中子數(shù)密度為ρn，質(zhì)子數(shù)密度為ρp，則有

ρ=ρn+ρp

(3)

再令

t=(ρn-ρp)

(4)

利用ρ介子的平均值方程

(5)

則有

(6)

故由中子和質(zhì)子非對(duì)稱(chēng)性所貢獻(xiàn)的單個(gè)核子能量密度為

(7)

利用積分公式

并令核子有效質(zhì)量為

m*=m-gσσ

(8)

則(7)式積分后化為

(9)

對(duì)稱(chēng)能系數(shù)為單核子能量對(duì)t的二階偏導(dǎo)數(shù)

(10)

易求

(11)

為了計(jì)算(9)式后面幾項(xiàng)對(duì)t的二階偏導(dǎo)，令

(12)

聯(lián)立(3)、(4)式得

(13)

費(fèi)米動(dòng)量和飽和核密度的關(guān)系是

(14)

把(14)式代入(13)式，得

(15)

所以，

(16)

由此，得

(17)

由(15)式可知，t=0即意味著kn=kp=kF，因此

(18)

同理可得，

(19)

又因?yàn)?/p>

(20)

所以

(21)

同理，

(22)

把(9)、(11)、(18)、(19)、(21)和(22)式代入(10)式，得對(duì)稱(chēng)能系數(shù)為

(23)

由此得ρ介子的耦合參數(shù)為

(24)

由(14)式，上式亦可寫(xiě)為

(25)

式(24)、(25)即為計(jì)算ρ介子耦合參數(shù)的代數(shù)表達(dá)式，據(jù)此，ρ介子的耦合參數(shù)gρ可以利用飽和核物質(zhì)的密度ρ0、對(duì)稱(chēng)能系數(shù)asym和有效質(zhì)量m*給出。

按照上述公式求出的ρ介子耦合參數(shù)gρ，適用于飽和核物質(zhì)；至于是否可以外推到如中子星這樣的高密度物質(zhì)中去，要視實(shí)驗(yàn)或觀測(cè)數(shù)據(jù)而作相應(yīng)調(diào)整。當(dāng)然，把計(jì)算出的核子耦合參數(shù)應(yīng)用于中子星物質(zhì)時(shí)，還需要考慮到中子星物質(zhì)的化學(xué)平衡條件、電中性條件和粒子數(shù)守恒條件。

2 ρ介子耦合參數(shù)的分析及新的耦合參數(shù)的計(jì)算

在相對(duì)論平均場(chǎng)理論的計(jì)算中，人們?cè)褂眠^(guò)多組耦合參數(shù)(見(jiàn)表1)。這幾組數(shù)據(jù)，有的是直接由上述公式計(jì)算得出，在核物理及核天體物理的理論計(jì)算中與實(shí)驗(yàn)或觀測(cè)結(jié)果符合得較好；有的是由上述公式算出后，其結(jié)果與實(shí)驗(yàn)或觀測(cè)相去較遠(yuǎn)，因此，又根據(jù)有限核的單粒子性質(zhì)作了相應(yīng)調(diào)整[5]。

表1 相對(duì)論平均場(chǎng)理論計(jì)算中用過(guò)的ρ介子的耦合參數(shù)

對(duì)于耦合參數(shù)組GL97，文獻(xiàn)[5]給出的ρ介子的耦合參數(shù)gρ數(shù)值是4.791，與我們的計(jì)算值4.794 2相符，因此，文獻(xiàn)給出的該值是直接由計(jì)算得來(lái)的。對(duì)于耦合參數(shù)組GL85，文獻(xiàn)[4]給出的gρ數(shù)值是6.2，與我們的計(jì)算值6.35相差2.4%，文獻(xiàn)給出的值是直接由計(jì)算得來(lái)并作了很小調(diào)整(見(jiàn)圖1)。

圖1 耦合參數(shù)組GL97和GL85給出的gρ與本文計(jì)算值的比較

圖2為耦合參數(shù)組NL3，TM1，TM2和NLSH給出的gρ值與本文計(jì)算值的比較. 由圖2及表1可見(jiàn)，文獻(xiàn)中給出的gρ值分別為：NL3-1.339 1，TM1-1.409 2，TM2-1.437 4，NLSH-1.284 9；而本文根據(jù)相應(yīng)文獻(xiàn)提供的核子參數(shù)計(jì)算的gρ值分別為：NL3-5.399 8，TM1-5.649 4，TM2-5.712，NLSH-5.153 4，二者相差太遠(yuǎn). 原因是NL3和TM2適合描述輕核，TM1適合描述Z>20的重核，而NLSH仔細(xì)考慮了核的同位旋性質(zhì)，適合計(jì)算遠(yuǎn)離β穩(wěn)定線的核[5]，直接計(jì)算值必須要利用有限核的單粒子性質(zhì)作調(diào)整，才能與無(wú)限核物質(zhì)和核天體物理的實(shí)驗(yàn)和觀測(cè)結(jié)果相符。

圖2 耦合參數(shù)組NL3，TM1，TM2和NLSH給出的gρ與本文計(jì)算值的比較

2004年，在NSCL/MSU進(jìn)行的112Sn和114Sn重離子碰撞試驗(yàn)給出的對(duì)稱(chēng)能等有關(guān)核子實(shí)驗(yàn)數(shù)據(jù)見(jiàn)表2[16]。核子有效質(zhì)量，當(dāng)采用MDI物態(tài)方程(the MDI EOS)時(shí)取為m*/m=0.67，當(dāng)采用Skyrme擬合(the Skyrme fit)時(shí)取為m*/m=0.77。 利用NSCL/MSU實(shí)驗(yàn)數(shù)據(jù)計(jì)算的ρ介子的耦合參數(shù)如表2和圖3，它們介于GL85、GL97與NL3、TM1、TM2和NLSH之間，此計(jì)算值是否可以外推到中子星等高密度核物質(zhì)中去，是否與核物理及核天體物理的實(shí)驗(yàn)或觀測(cè)結(jié)果符合，都需具體的理論計(jì)算與實(shí)驗(yàn)或觀測(cè)結(jié)果來(lái)驗(yàn)證。如不符合，尚需對(duì)其進(jìn)行適當(dāng)調(diào)整。

表2 NSCL/MSU給出的核子參數(shù)及本文據(jù)此計(jì)算的ρ介子耦合參數(shù)

圖3 利用NSCL/MSU重離子碰撞數(shù)據(jù)計(jì)算的ρ介子耦合參數(shù)

3 結(jié) 論

本文從相對(duì)論平均場(chǎng)理論出發(fā)，導(dǎo)出了ρ介子與核子的耦合參數(shù)gρ與飽和核物質(zhì)密度ρ0、對(duì)稱(chēng)能量系數(shù)asym及核子有效質(zhì)量m*之間的代數(shù)關(guān)系式；由得到的公式對(duì)相對(duì)論平均場(chǎng)理論計(jì)算中曾用到的幾組ρ介子耦合參數(shù)進(jìn)行了分析。 采用NSCL/MSU給出的核子最新實(shí)驗(yàn)數(shù)據(jù)，對(duì)ρ介子與核子的耦合參數(shù)gρ進(jìn)行了計(jì)算并作了初步討論。

參考文獻(xiàn)：

[1]LALAZISSIS G A, KONIG J, RING P. A new parameterization for the Lagrangian density of relativistic mean field theory[J]. Phys Rev C, 1997, 55(1):540-543.

[2]SUGAHARA Y, TOKI H. Relativistic mean field theory for lambda hypernuclei and neutron stars[J]. Prog Theor Phys, 1994, 92(4):803-813.

[3]SHARMA M M，NAGARAJAN M A , RING P. Rho meson coupling in the relativistic mean field theory and description of exotic nuclei[J]. Phys Lett B, 1993, 312(4):377-381.

[4]GLENDENNING N K. Neutron stars are giant hypernuclei?[J]. Ap J, 1985,293: 470.

[5]GLENDENNING N K. Compact stars: neuclear physics,particle physics, and general relativety[M]. New York:Springer-Verlag, 1997.

[6]TSANG M B , FRIEDMAN W A , GELBKE C K, et al. Isotopic scaling in nuclear reactions[J]. Phys Rev Lett, 2001, 86(22):5023-5026.

[7]ONO A , DANIELEWICZ P, FRIEDMAN W A , et al. Isospin fractionation and isoscaling in dynamical simulations of nuclear collisions[J]. Phys Rev C, 2003, 68(5):051601.

[8]LI Bao an. Neutron-proton differential flow as a probe of isospin-dependence of the nuclear equation of state[J]. Phys Rev Lett, 2000, 85(20):4221-4224.

[9]LI Bao an. Probing the high density behaviour of nuclear symmetry[J]. Phys Rev Lett, 2002, 88:192701.

[10]CHEN Lie wen, GRECO V, KO C M, et al. Effects of symmetry energy on two-nucleon correlation functions in heavy-ion collisions induced by neutron-rich nuclei[J]. Phys Rev Lett, 2003, 90(16): 162701.

[11]TAN W P, LI Bao an , DONANGELO R, et al. Fragment isotope distributions and the isospin dependent equation of state[J]. Phys Rev C, 2001, 64(5): 051901(R).

[12]MULLER H, SEROT B D. Phase transitions in warm, asymmetric nuclear matter[J]. Phys Rev C , 1995, 52(4):2072-2091.

[13]LI Bao an, KO C M. Chemical and mechanical instability in hot isospin-asymmetric nuclear matter[J]. Nucl Phys A， 1997(618): 498-508.

[14]BARAN V, COLONNA M, TORO M D, et al. Spinodal decomposition of low-density asymmetric nuclear matter[J]. Nucl Phys A, 1998, 632(2):287-303.

[15]XU H S, TSANG M B, LIU T X, et al. Isospin fractionation in nuclear multifragmentation[J]. Phys Rev Lett, 2000, 85(4):716-719.

[16]TSANG M B, LIU T X, SHI L, et al. Isospin diffusion and the nuclear symmetry energy in heavy ion reactions[J]. Phys Rev Lett， 2004, 92(6):062701.