王再軍

（天津職業(yè)技術師范大學理學院，天津 300222）

豐中子奇A核37-47P自旋和宇稱的相對論平均場模型研究

王再軍

（天津職業(yè)技術師范大學理學院，天津 300222）

用相對論平均場模型結合NL3、TM2和NL-SH參數(shù)對P的奇A核同位素的2s1/2和1d3/2單質子能級變化趨勢以及自旋和宇稱進行了理論計算。發(fā)現(xiàn)豐中子同位素37-47P，尤其是47P中可能存在2s1/2和1d3/2單質子能級反轉。對于37P和39P，TM2和NL-SH參數(shù)給出的自旋和宇稱是（3/2）+，NL3給出的是（1/2）+；對于41-47P，NL3、TM2和NL-SH給出的自旋和宇稱都是（3/2）+。自旋和宇稱的結果只有部分和NUBASE2012中的估計值是一致的。

相對論平均場模型；自旋；宇稱；能級反轉；P同位素

自旋和宇稱是原子核的兩個基本性質。在大多數(shù)情況下，研究原子核的其他性質時，或是研究有原子核參與的核過程時，所涉及的原子核的自旋和宇稱都需要是已知的。對于穩(wěn)定原子核，自旋和宇稱的實驗測量早已不是問題?，F(xiàn)在，所有已知穩(wěn)定的原子核和大多數(shù)長壽命原子核的自旋和宇稱都已確定。而對于很多不穩(wěn)定的原子核，由于其壽命短，目前實驗上準確測量它們的自旋和宇稱還很困難。在最新發(fā)布的NUBASE2012[1]核數(shù)據(jù)表中，大多數(shù)短壽命的原子核的自旋和宇稱都是根據(jù)相同質子數(shù)或相同中子數(shù)的臨近原子核的自旋和宇稱趨勢估算得來的。本文將用相對論平均場理論對P原子核的奇核子同位素進行計算，對其豐中子同位素37-47P的自旋和宇稱進行預測。P原子核的奇核子同位素由偶數(shù)個中子和15個質子組成。根據(jù)原子核的殼模型理論，這些同位素的自旋和宇稱應由其最外面的一個價質子決定。但是，由于P原子核屬于s-d殼區(qū)，其某些豐中子同位素中可能會出現(xiàn)2s1/2能級和1d3/2能級反轉，這就會導致這些同位素的自旋和宇稱發(fā)生變化。本文將用相對論平均場理論結合NL3[2-3]、TM2[4]和NL-SH[5]參數(shù)計算這些變化，以期使這些計算結果能對豐中子P原子核的研究提供有用的參考，并且能對相對論平均場理論及其不同的參數(shù)對不穩(wěn)定原子核的有效性檢驗提供新的結果。

1 相對論平均場理論簡述

原子核相對論平均場模型現(xiàn)在已經(jīng)發(fā)展成為一個標準的原子核結構理論，被廣泛用來描述原子核的結構。大量的計算表明[6-9]，該理論可以很好地描述穩(wěn)定原子核甚至是不穩(wěn)定原子核的基態(tài)和低激發(fā)態(tài)的一些性質。該理論的詳細內容可以參見文獻[6-10]，這里只對其主要內容做一簡述。在相對論平均場模型中，包含核子場Ψ、σ介子場、ω介子場、ρ介子和電磁場Aμ作用項的原子核體系的拉格朗日密度為：

其中：

這里：m、mσ、mω和mρ分別為核子、σ介子、ω介子和ρ介子的質量；α=1/137為電磁作用常數(shù)。核子場和介子場之間的有效作用常數(shù)分別為gσ、gω和gρ。g2和g3為σ介子場的非線性耦合常數(shù)，c3表示ω介子場的自作用常數(shù)；τa和τ3分別為泡利矩陣及其第三分量。

在不考慮費米海的情況下，利用變分原理可以得到一組耦合方程為：

這里：M*（r）=m+gσσ（r）為σ介子的有效質量；V（r）是洛侖茲矢量的類時分量；ρs、ρv和ρp分別為標量粒子，重子和質子的密度，ρ3為中子和質子的密度差。上面這組耦合方程可以用自洽疊代方法進行數(shù)值求解得到核子的分布，然后可以進一步計算原子核的結合能、質子和中子分布的均方根半徑、單粒子能級等。

2 數(shù)值結果和討論

這一部分給出相對論平均場理論對奇A核25-47P的計算結果，并對結果進行討論。表1列出了NL3、TM2和NL-SH參數(shù)計算的平均結合能（E/A）和電荷分布的均方根半徑（Rrms），同時也給出了最新發(fā)布的NUBASE2012核數(shù)據(jù)表中的結果。

表1 25-47P的平均結合能和電荷分布的均方根半徑，以及NUBASE2012核數(shù)據(jù)表中的結果和31P的電荷分布均方根半徑的實驗數(shù)據(jù)

表中31P的電荷分布均方根半徑的實驗值來自文獻[11]。另外，31P電荷分布的均方根半徑還有2個實驗結果，分別是3.193 fm和3.187 fm[12]。對比3套參數(shù)給出的平均結合能和電荷分布的均方根半徑以及NUBASE2012給出的數(shù)據(jù)，可以得到如下結果：①對于平均結合能，NL3、TM2和NL-SH3套參數(shù)給出的結果的絕對差值不超過0.18 MeV，相對差值不超過2.19%；給出的平均結合能與NUBASE2012給出的數(shù)據(jù)的絕對差值不超過0.16 MeV，相對差值不超過2.33%。②對于電荷分布的均方根半徑，3套參數(shù)給出的結果的絕對差值不超過0.079 fm，相對差值不超過2.38%。③對于31P，3套參數(shù)給出的平均結合能的理論值與實驗值的絕對偏離只有0.156 MeV，相對偏離為1.87%；給出的電荷分布的均方根半徑與實驗值的絕對偏離只有0.047 fm，相對差值為1.47%。為了能看出3套參數(shù)給出的理論值與實驗值的總體變化關系和趨勢，在圖1中給出了平均結合能的理論值和實驗值隨中子數(shù)的變化曲線，圖2給出了電荷分布的均方根半徑隨中子數(shù)的變化曲線，圖中的空心圓點為31P的實驗值。

從圖1中看到，3套參數(shù)給出平均結合能和NUBASE2012給出的結果都相差很小，趨勢一致；從圖2中看到，3套參數(shù)給出的電荷分布的均方根半徑變化趨勢也一致，NL-SH和NL3給出的結果比較接近，但TM2給出的結果相對偏大，對于31P，實驗值更接近于NL-SH和NL3的預言值。從上述分析得出，不論是穩(wěn)定同位素還是不穩(wěn)定同位素，相對論平均場理論對P原子核的計算結果總體是可靠的。但進一步對單核子能級的計算發(fā)現(xiàn)，3套參數(shù)對37-47P的計算結果有一定的差別。

圖1 平均結合能的理論值和實驗值隨中子數(shù)的變化曲線

圖2 3套參數(shù)給出的電荷分布的均方根半徑變化趨勢

表2列出了用3套參數(shù)所計算的25-47P的2s1/2和1d3/2單質子能級，表中ε（2s1/2）和ε（1d3/2）分別代表2s1/2和1d3/2能級。這里著重分析能級的變化趨勢。從表中看到，對于25-35P，3套參數(shù)均預言2s1/2能級低于1d3/2能級。而對于37P和39P，3套參數(shù)的結果就有了差別。其中NL3預言2s1/2能級低于1d3/2能級，而NL-SH、TM2則預言2s1/2能級高于1d3/2能級，也就是在37P和39P中可能存在2s1/2和1d3/2能級次序的反轉。對于41-47P，3套參數(shù)再一次給出相同的結果，但都是2s1/2能級高于1d3/2能級。也就是說，3套參數(shù)都預言41-47P中存在2s1/2和1d3/2能級次序的反轉。由此可知，不同的參數(shù)對37-47P的預言還是有很大差別的，這個差別表現(xiàn)在能級的具體數(shù)值上盡管很小，但卻對另一個物理量產(chǎn)生了本質的影響，那就是原子核的自旋和宇稱。同時這也說明，原子核中子的過剩對費米面附近的質子的運動狀態(tài)會產(chǎn)生影響。

表2 25-47P的2s1/2和1d3/2單質子能級

表3 NL3、TM2和NL-SH預言的自旋和宇稱值，以及NBUBASE2012中給出的結果

粒子的運動總是傾向于使自己處于能量最低的狀態(tài)。因此，2s1/2和1d3/2能級次序的改變，會對37-47P中質子的運動狀態(tài)產(chǎn)生重要影響。根據(jù)原子核的殼模型理論和表2的結果，不難給出由NL3、TM2和NL-SH參數(shù)所預言的25-47P的自旋和宇稱，結果見表3。同時，在表3中也列出了NUBASE2012的結果，其中#號表示估計結果。對表中的結果對比分析，發(fā)現(xiàn)一些很有意義的現(xiàn)象。從表中可以看到，對于25-35P，3套參數(shù)的預言結果與NUBASE2012給出的結果一致，這包括所有有實驗結果的同位素27-35P。這再一次說明相對論平均場對穩(wěn)定核是可靠的。而對于其余的同位素，3套參數(shù)和NUBASE2012給出的結果各有不同。對于37P和39P，NL3給出的結果是（1/2）+，與NUBASE2012給出的結果一致；而TM2和NL-SH參數(shù)給出的結果為（3/2）+。這是一個很有意義的結果。通常情況下，相對論平均場的不同參數(shù)在計算同一原子核時所給出的結合能、平均結合能、電荷均方根半徑等結果不會表現(xiàn)出本質的差別，因此很難從這些結果去判斷哪一套參數(shù)更合適。而自旋和宇稱的變化確是本質的。因此，這個結果就給判斷哪一套參數(shù)更適合于描述P原子核提供了參考。對于41-45P，3套參數(shù)所預言的結果均為（3/2）+，與NUBASE2012給出的結果不一致。對于47P，NL3、TM2和NL-SH3套參數(shù)所預言的結果均為（3/2）+，同時也和NBUBASE2012給出的結果一致。3套參數(shù)的理論結果均與實驗結果相同這種一致性，強烈預示在47P中可能存在2s1/2和1d3/2單質子能級次序反轉。

3 結束語

用相對論平均場理論結合NL3、TM2和NL-SH參數(shù)對P原子核的奇核子同位素的2s1/2和1d3/2單質子能級隨中子數(shù)的變化趨勢以及自旋和宇稱進行了理論計算，并與最新發(fā)布的NBUBASE2012給出的結果進行了對比討論。研究發(fā)現(xiàn)，在豐中子奇核子同位素37-47P中，尤其是47P中，可能存在著2s1/2和1d3/2單質子能級反轉。對于25-35P的自旋和宇稱，3套參數(shù)和NBUBASE2012給出的結果都是（1/2）+，與實驗結果一致；對于37-47P的自旋和宇稱，3套參數(shù)和NBUBASE2012給出的結果卻不盡相同。理論預言37P和39P的自旋和宇稱可能是（1/2）+，也可能是（3/2）+，與NBUBASE2012給出的結果有差別；理論預言41-45P的自旋和宇稱是（3/2）+，與NBUBASE2012給出的結果完全不同；47P的自旋和宇稱為（3/2）+，與NBUBASE2012給出的結果完全一致。該研究結果可以為豐中子P原子核以及鄰近原子核性質的研究提供參考，同時也可以為相對論平均場理論及其不同的參數(shù)對不穩(wěn)定原子核的有效性檢驗提供新的結果。

［1］ AUDI G，KONDEV F G，WANG M，et al.The Nubase2012 evaluation of nuclear properties［J］．Chinese Physics C，2012，36：1157-1286.

［2］ LALAZISSIS G A，KONIG J，RING P.New parametrization for the Lagrangian density of relativistic mean field theory［J］．Phys Rev C，1997，55：540.

［3］ PATYK Z，BARAN A，BERGER J F，et al.Masses and radii ofsphericalnucleicalculated in various microscopic approaches［J］．Phys Rev C，1999，59：704.

［4］ SUGAHARA Y，TOKI H.Relativistic mean-field theory for unstable nuclei with non-linear σ and ω terms［J］．Nucl Phys A，1994，579：557.

［5］ SHARMA M M，NAGARAJIAN M A，RING P.Rho meson coupling in the relativistic mean field theory and description of exotic nuclei［J］．Phys Lett B，1993，312：377.

［6］ WALECKA J D.A theory of highly condensed matter［J］．Annals of Physics（NY），1994，83：491-529.

［7］ SEROT B D，WALECKA J D.The relativistic nuclear manybody problem［J］.Advances in Nuclear Physics，1986，16：1-327.

［8］ HOROWITZ C J，SEROT B D.Selfconsistent hartree description of finite nuclei in a relativistic quantum field theory［J］.Nuclear Physics A，1981，368：503-528.

［9］ MA Z Y，SHI H L，CHEN B Q.Isovector meson contribution in the relativistic Hartree-Fock approach for finite nuclei［J］. Physical Review C，1994，50：3170-3173.

［10］REN Z Z，MITTIN W，SARAZIN F.Relativistic mean-field study on proton skins and proton halos in exotic nuclei［J］. Nuclear Physics A，1999，652：250-270.

［11］ WESSELING J，JAGE C W，LAPIKAS L，et al.2s1/2occupancies in30Si，31P，and32S［J］.Phys Rev C，1997，55：2773.

［12］VRIES H，JAGER C W.Nuclear charge-density-distribution parameters from elastic electron scattering［J］.At Data Nucl Data Tables，1987，36：495.

Spins and parities of neutron-rich odd-A37-47P and relativistic mean field model

WANG Zai-jun
（School of Science，Tianjin University of Technology and Education，Tianjin 300222，China）

The 2s1/2and 1d3/2single-proton energy levels and the spin and parity values of the odd-A P isotopes are studied by using the relativistic mean field model with the NL3，TM2 and NL-SH parameters.Results show that the 2s1/2and 1d3/2proton level inversion may exist in the neutron-rich isotopes37-47P，especially in47P.For the spins and parities，the predictions for37P and39P with TM2 and NL-SH parameters are （3/2）+，and those with NL3 are （1/2）+；the predictions for41-47P with TM2，NL-SH and NL3 are all （3/2）+.Only part of the results of the predicted spin and parity values agree with the estimates given by the NUBASE2012.

relativistic mean field model；spin；parity；level inversion；P isotopes

O571.21

A

2095－0926（2014）01－0001－04

2014-02-25

國家自然科學基金項目（11275138，10675090）；天津職業(yè)技術師范大學科研發(fā)展基金項目（KJYB11-3）.

王再軍（1964—），男，教授，研究方向為理論核物理.