第30卷增刊Ⅱ2006年12月
高能物理与核物理
HIGH
ENERGY
PHYSICS
AND
NUCLEAR
PHYSICS
Vol.30,Supp.Ⅱ
Dec.,2006
相对论平均场理论对重核中结团现象的研究
王士虎
1)
郭建友
合肥
230023)
(安徽大学物理与材料科学学院
摘要用相对论平均场理论结合结团模型对重核结团现象进行了研究.首先将原子核分为两个部分,然后用相对论平均场模型对原子核的核子按束缚能的高低来分布,结果获得了与实验相一致的结团.这表明相对论平均场模型也适用于重核结团现象的研究.关键词相对论平均场重核结团
1引言
重核结团的研究中占据很重要的地位.自从
[1]
2理论框架
相对论平均场理论(RMF)
[15,16]
是一个相当成功
Sandulescu等预言重的结团自发裂变以来,人们陆续发现了很多重的结团自发裂变现象(14C,
28
20
[2]
的理论,它解决了原子核领域的很多难题,我们在这里同样是用它来解决重核结团的.首先,我们给出相对论平均场理论(RMF)的拉格朗日密度:
L=iγµ∂µ−m−gσσ−gωγµ−gργµτρµ−
•
O,
24
Ne,
Mg,
32,34
Si···),导致了对重结团结构的关注.目
[4]
[5]
[3]
[15,17—20]
前人们建立了许多结团模型,如统一模型,超非对称裂变模型分析(ASAFM),有效液滴模型(ELDM),密度依赖的结团模型(DDCM)
[6—11]
等.
[12]
除此之外,人们建立了既包括结团的描述又包括平均场的描述的费密子的分子动力学方法(FMD)反对称分子动力学方法(AMD)等
[14]
[13]
1−τ11µ32eγµAµΨ+∂σ∂µσ−m2σσ−222111
g2σ3−g3σ4−ωµνωµν+344112µ
mσωωµν+c3(ωµωµ)2−241µν11µρρµν+m2ρµ−AµνAµν,ρρ424
•
,
等,也已经可以很
好地解决轻核结团结构的模型.最近,P.Arumugam
用相对论平均场模型(RMF)给出的分布密度也对轻核结团结构给出了一个清晰的物理图象,并且与实验结果很好的符合.但是,由于有些重核中存在较强的自旋–轨道耦合,他们无法用由相对论平均场模型(RMF)给出的分布密度得到这些重核中的结团结构
[14]
(1)
其中Ψ是核子的Dirac旋量,对应的质量为m.ω和ρ分别是isoscalar-scalar和isoscalar-vector介子.σ是isovector-vector介子.光子场A提供了原子核中的电磁属性.相关张量如下:
ωµν=∂ων−∂ωµ,Aµν=∂Aν−∂Aµ,
ρµν=∂ρν−∂ρµ−2gρρµ×ρν,
µ
ν
µ
ν
µ
ν
.
本文根据结团放射性的理论中的一些基本假设(原子核分成子核和结团两部分,这两部分与单个原子核具有相似的物质密度分布)
[6—11]
将原子核分成了内(2)(3)(4)
外两部分,利用相对论平均场模型(RMF)对这两个部分结合能进行了计算,很自然的导出了重核结团结构.通过对目前已知的大部分具有重核结团放射性的重核进行计算,发现与实验符合的比较好.
将拉格朗日密度代入Euler-Lagrange方程,我们很快得到核子的场方程.由于我们只涉及到原子核的基态
1)E-mail:[email protected]
87—89
88高能物理与核物理(HEP&NP)第30卷
性质,我们可以设介子场和光子场是静态的经典场,核子在经典场中做独立运动.
3
计算细节
对于参数的选择,我们对原子核采用了较好的
符合重核性质的一组参数NL3.值得注意的是,对关联对重核的定量计算是十分重要的,所以我们采用了简洁而有效的固定能隙的BCS(Barden-Cooper-Schrieffer)方法处理对关联:∆=12
.将对相互作用
作为微扰项,对不成对的核子采用泡利阻塞的办法处理.在这里,我们采用变形轴对称谐振子基展开的办法(主壳Nb=Nf=20),参数β2由核电四极矩Q给出.这里的Q为
Q=Qn+Qp=
9
5π
AR2β2,(5)
其中R=1.2A1
fm,输入的形变参数β的初值选取任意合理的值.
我们以如下方式考虑要解决的问题.先假设将某个原子核分成两部分,则这两部分应可视为两个新的独立的原子核(当然这时原子核不一定稳定)
[6—11]
.由
于原子核之间是弱关联而原子核内部是强关联,所以两个原子核之间的作用对于原子核内部的强关联而言只是微扰.而任意系统在无外界能量激发下基态的粒子总是向更高结合能处迁跃,所以两个新原子核的最外层核子的迁跃运动仅仅取决于两个原子核最外层结合能的高低.当任意一个原子核的最外层不再俘获或失去核子时,两个原子核就稳定了.此时,比较小的原子核就应该是分裂出来的结团.而原子核内部的结合能可以根据相对论平均场模型(RMF)计算,所以我们用相对论平均场模型(RMF)可以算出这2个原子核的核束缚能(大核用形变核计算,小核用球形核计算).由此可以确定可能出现的结团以及结团的结构.更复杂的问题也能依此解决.
为简化问题,实际计算中可以先假设分裂出一个质子.又由于对力的存在,2个核子在一起时的结合能比单个核子高,因此实际计算中考虑以2个核子为一个基本单位:
S2p=EB(Z−2,A−2)−EB(Z,A),(6)S2n=EB(Z,A−2)−EB(Z,A),
(7)
其中EB(Z,A)是质子数为Z,质量数为A的核束缚能.与此同时,我们还将实验所得出的结合能与之进行了
比较.
4
结果及讨论
我们知道结团要形成,必须它本身比较稳定.中
子数目的范围随质子数而定,而质子数是可以连续变化的.所以,为了减少计算量我们首先来看看中子数随质子数改变的情况.由于纵坐标是能级差,所以当某一点为负值,而它的右面为正值时该点代表的结团既不俘获新的核子,其最外层核子也不被别的核俘获(注:由于每个原子核里都有4He结团,所以有些原子核省略了,纵坐标的值是大核最外层中子或质子的结
合能减去小核最外层中子或质子的结合能).
图1
Numericalrelationsbetweenneutronsandprotonsinsmallparts——
224Ra.
图2Numericalrelationsbetweenneutronsandprotonsinsmallparts——224Ra.
从图1,2可知,224
Ra中有4He结团(这点我们在
其他原子核的分析中可以看到,因为每个原子核都一样,所以我们省略了)还有18O−4He=14C结团.其中,我们注意到16C对最外层的两个核子的束缚极极小,
增刊Ⅱ王士虎等:相对论平均场理论对重核中结团现象的研究89
在小结团受到大结团作用时(该作用相对结团内部作用来说极小)会失去两个中子.所以就不考虑20O时最外层质子的情况.
子核的结构获得更深入的了解.由于本理论只涉及到结合能的计算,在具体应用中几乎不存在限制条件.
但是,我们也注意到某些核子的分析结果与实验存在一定差距.这有可能是由于原子核衰变时使α结团发生了分裂,也有可能是由于(RMF)的误差.该理论的成功说明考虑了BCS的RMF可以对重核的结团结构作出较好的描述,这也为结团现象的研究开辟了一条新的途径.
5结论
我们通过以上分析可以看到大部分本文所讨论的
原子核都符合重核的结团放射性实验的数据.更进一步得到了重核内部的结团结构.这样,我们可以对原
参考文献(References)
1GreinerW.Z.Phys.,1994,A349:315
2SandulescuA,PoenaruDN,GreinerWetal.Sov.J.Part.Nucl.,1980,11:528
3RoseHJ,JonesGA.Nature(London),1984,307:2454PoenaruDNetal.Z.Phys.,1986,A325:435
5DuarteSBetal.At.DataNucl.DataTables,2002,80:235
6XUChang,RENZhong-Zhou.Phys.Rev.,2003,C68:034319
7XUChang,RENZhong-Zhou.Phys.Rev.,2004,C69:024614
8RENZhong-Zhou,XUChang,WANGYong-Qinetal.Nucl.Phys.,2004,A738:318
9XUChang,RENZhong-Zhou.Nucl.Phys.,2005,A753:174
10XUChang,RENZhong-Zhou.Nucl.Phys.,2005,A760:
303
11RENZhong-Zhou,XUChang,WANGZai-Jun.Phys.Rev.,
2004,C70:034304
12FeldmeierH,SchnackJ.Rev.Mod.Phys.,2000,72:65513Kanada-En’YoY,HoriuchiH.Prog.Theor.Phys.,2001,
142(Suppl.):205
14ArumugamP,SharmaBK,PatraSK.Phys.Rev.,2005,
C71:064308
15SerotBD,WaleckaJD.Adv.Nucl.Phys.,1986,16:116SerotBD.Rep.Prog.Phys.,1992,55:1855
17GambhirYK,RingP,ThimetA.Ann.Phys.,1990,198:
132
18SugaharaY,TokiH.Nucl.Phys.,1994,A579:55719PandaPK,PatraSK,ReinhardtJetal.Int.J.Mod.
Phys.,1997,E6:307
20GmucaS.Nucl.Phys.,1992,A547:447
RelativisticMeanFieldStudyofClusteringinHeavyNuclei
WANGShi-Hu1)
GUOJian-You
(DepartmentofPhysics,AnhuiUniversity,Hefei230039,China)
AbstractTheclusteringphenomenoninheavynucleiisdiscussedinthispaperwithintherelativisticmeanfield(RMF)approach.TheclusterscanbeobservedinexperimentsareobtainedthroughthecalculatednucleonbindingenergydistributionsaccordingtotheRMFtheoryaftertheatomicnucleusisdividedintotwoparts,insteadoftheusualmethodofanalysistothecalculatednucleondensitydistributionsbecauseofthestrongerspin-orbitcoupleinheavynuclei.TheresultsindicatethattheRMFtheoryiswellsuitedforthestudyofclusteringinheavynuclei.KeywordsRMF,clusteringinheavynuclei
1)E-mail:[email protected]
第30卷增刊Ⅱ2006年12月
高能物理与核物理
HIGH
ENERGY
PHYSICS
AND
NUCLEAR
PHYSICS
Vol.30,Supp.Ⅱ
Dec.,2006
相对论平均场理论对重核中结团现象的研究
王士虎
1)
郭建友
合肥
230023)
(安徽大学物理与材料科学学院
摘要用相对论平均场理论结合结团模型对重核结团现象进行了研究.首先将原子核分为两个部分,然后用相对论平均场模型对原子核的核子按束缚能的高低来分布,结果获得了与实验相一致的结团.这表明相对论平均场模型也适用于重核结团现象的研究.关键词相对论平均场重核结团
1引言
重核结团的研究中占据很重要的地位.自从
[1]
2理论框架
相对论平均场理论(RMF)
[15,16]
是一个相当成功
Sandulescu等预言重的结团自发裂变以来,人们陆续发现了很多重的结团自发裂变现象(14C,
28
20
[2]
的理论,它解决了原子核领域的很多难题,我们在这里同样是用它来解决重核结团的.首先,我们给出相对论平均场理论(RMF)的拉格朗日密度:
L=iγµ∂µ−m−gσσ−gωγµ−gργµτρµ−
•
O,
24
Ne,
Mg,
32,34
Si···),导致了对重结团结构的关注.目
[4]
[5]
[3]
[15,17—20]
前人们建立了许多结团模型,如统一模型,超非对称裂变模型分析(ASAFM),有效液滴模型(ELDM),密度依赖的结团模型(DDCM)
[6—11]
等.
[12]
除此之外,人们建立了既包括结团的描述又包括平均场的描述的费密子的分子动力学方法(FMD)反对称分子动力学方法(AMD)等
[14]
[13]
1−τ11µ32eγµAµΨ+∂σ∂µσ−m2σσ−222111
g2σ3−g3σ4−ωµνωµν+344112µ
mσωωµν+c3(ωµωµ)2−241µν11µρρµν+m2ρµ−AµνAµν,ρρ424
•
,
等,也已经可以很
好地解决轻核结团结构的模型.最近,P.Arumugam
用相对论平均场模型(RMF)给出的分布密度也对轻核结团结构给出了一个清晰的物理图象,并且与实验结果很好的符合.但是,由于有些重核中存在较强的自旋–轨道耦合,他们无法用由相对论平均场模型(RMF)给出的分布密度得到这些重核中的结团结构
[14]
(1)
其中Ψ是核子的Dirac旋量,对应的质量为m.ω和ρ分别是isoscalar-scalar和isoscalar-vector介子.σ是isovector-vector介子.光子场A提供了原子核中的电磁属性.相关张量如下:
ωµν=∂ων−∂ωµ,Aµν=∂Aν−∂Aµ,
ρµν=∂ρν−∂ρµ−2gρρµ×ρν,
µ
ν
µ
ν
µ
ν
.
本文根据结团放射性的理论中的一些基本假设(原子核分成子核和结团两部分,这两部分与单个原子核具有相似的物质密度分布)
[6—11]
将原子核分成了内(2)(3)(4)
外两部分,利用相对论平均场模型(RMF)对这两个部分结合能进行了计算,很自然的导出了重核结团结构.通过对目前已知的大部分具有重核结团放射性的重核进行计算,发现与实验符合的比较好.
将拉格朗日密度代入Euler-Lagrange方程,我们很快得到核子的场方程.由于我们只涉及到原子核的基态
1)E-mail:[email protected]
87—89
88高能物理与核物理(HEP&NP)第30卷
性质,我们可以设介子场和光子场是静态的经典场,核子在经典场中做独立运动.
3
计算细节
对于参数的选择,我们对原子核采用了较好的
符合重核性质的一组参数NL3.值得注意的是,对关联对重核的定量计算是十分重要的,所以我们采用了简洁而有效的固定能隙的BCS(Barden-Cooper-Schrieffer)方法处理对关联:∆=12
.将对相互作用
作为微扰项,对不成对的核子采用泡利阻塞的办法处理.在这里,我们采用变形轴对称谐振子基展开的办法(主壳Nb=Nf=20),参数β2由核电四极矩Q给出.这里的Q为
Q=Qn+Qp=
9
5π
AR2β2,(5)
其中R=1.2A1
fm,输入的形变参数β的初值选取任意合理的值.
我们以如下方式考虑要解决的问题.先假设将某个原子核分成两部分,则这两部分应可视为两个新的独立的原子核(当然这时原子核不一定稳定)
[6—11]
.由
于原子核之间是弱关联而原子核内部是强关联,所以两个原子核之间的作用对于原子核内部的强关联而言只是微扰.而任意系统在无外界能量激发下基态的粒子总是向更高结合能处迁跃,所以两个新原子核的最外层核子的迁跃运动仅仅取决于两个原子核最外层结合能的高低.当任意一个原子核的最外层不再俘获或失去核子时,两个原子核就稳定了.此时,比较小的原子核就应该是分裂出来的结团.而原子核内部的结合能可以根据相对论平均场模型(RMF)计算,所以我们用相对论平均场模型(RMF)可以算出这2个原子核的核束缚能(大核用形变核计算,小核用球形核计算).由此可以确定可能出现的结团以及结团的结构.更复杂的问题也能依此解决.
为简化问题,实际计算中可以先假设分裂出一个质子.又由于对力的存在,2个核子在一起时的结合能比单个核子高,因此实际计算中考虑以2个核子为一个基本单位:
S2p=EB(Z−2,A−2)−EB(Z,A),(6)S2n=EB(Z,A−2)−EB(Z,A),
(7)
其中EB(Z,A)是质子数为Z,质量数为A的核束缚能.与此同时,我们还将实验所得出的结合能与之进行了
比较.
4
结果及讨论
我们知道结团要形成,必须它本身比较稳定.中
子数目的范围随质子数而定,而质子数是可以连续变化的.所以,为了减少计算量我们首先来看看中子数随质子数改变的情况.由于纵坐标是能级差,所以当某一点为负值,而它的右面为正值时该点代表的结团既不俘获新的核子,其最外层核子也不被别的核俘获(注:由于每个原子核里都有4He结团,所以有些原子核省略了,纵坐标的值是大核最外层中子或质子的结
合能减去小核最外层中子或质子的结合能).
图1
Numericalrelationsbetweenneutronsandprotonsinsmallparts——
224Ra.
图2Numericalrelationsbetweenneutronsandprotonsinsmallparts——224Ra.
从图1,2可知,224
Ra中有4He结团(这点我们在
其他原子核的分析中可以看到,因为每个原子核都一样,所以我们省略了)还有18O−4He=14C结团.其中,我们注意到16C对最外层的两个核子的束缚极极小,
增刊Ⅱ王士虎等:相对论平均场理论对重核中结团现象的研究89
在小结团受到大结团作用时(该作用相对结团内部作用来说极小)会失去两个中子.所以就不考虑20O时最外层质子的情况.
子核的结构获得更深入的了解.由于本理论只涉及到结合能的计算,在具体应用中几乎不存在限制条件.
但是,我们也注意到某些核子的分析结果与实验存在一定差距.这有可能是由于原子核衰变时使α结团发生了分裂,也有可能是由于(RMF)的误差.该理论的成功说明考虑了BCS的RMF可以对重核的结团结构作出较好的描述,这也为结团现象的研究开辟了一条新的途径.
5结论
我们通过以上分析可以看到大部分本文所讨论的
原子核都符合重核的结团放射性实验的数据.更进一步得到了重核内部的结团结构.这样,我们可以对原
参考文献(References)
1GreinerW.Z.Phys.,1994,A349:315
2SandulescuA,PoenaruDN,GreinerWetal.Sov.J.Part.Nucl.,1980,11:528
3RoseHJ,JonesGA.Nature(London),1984,307:2454PoenaruDNetal.Z.Phys.,1986,A325:435
5DuarteSBetal.At.DataNucl.DataTables,2002,80:235
6XUChang,RENZhong-Zhou.Phys.Rev.,2003,C68:034319
7XUChang,RENZhong-Zhou.Phys.Rev.,2004,C69:024614
8RENZhong-Zhou,XUChang,WANGYong-Qinetal.Nucl.Phys.,2004,A738:318
9XUChang,RENZhong-Zhou.Nucl.Phys.,2005,A753:174
10XUChang,RENZhong-Zhou.Nucl.Phys.,2005,A760:
303
11RENZhong-Zhou,XUChang,WANGZai-Jun.Phys.Rev.,
2004,C70:034304
12FeldmeierH,SchnackJ.Rev.Mod.Phys.,2000,72:65513Kanada-En’YoY,HoriuchiH.Prog.Theor.Phys.,2001,
142(Suppl.):205
14ArumugamP,SharmaBK,PatraSK.Phys.Rev.,2005,
C71:064308
15SerotBD,WaleckaJD.Adv.Nucl.Phys.,1986,16:116SerotBD.Rep.Prog.Phys.,1992,55:1855
17GambhirYK,RingP,ThimetA.Ann.Phys.,1990,198:
132
18SugaharaY,TokiH.Nucl.Phys.,1994,A579:55719PandaPK,PatraSK,ReinhardtJetal.Int.J.Mod.
Phys.,1997,E6:307
20GmucaS.Nucl.Phys.,1992,A547:447
RelativisticMeanFieldStudyofClusteringinHeavyNuclei
WANGShi-Hu1)
GUOJian-You
(DepartmentofPhysics,AnhuiUniversity,Hefei230039,China)
AbstractTheclusteringphenomenoninheavynucleiisdiscussedinthispaperwithintherelativisticmeanfield(RMF)approach.TheclusterscanbeobservedinexperimentsareobtainedthroughthecalculatednucleonbindingenergydistributionsaccordingtotheRMFtheoryaftertheatomicnucleusisdividedintotwoparts,insteadoftheusualmethodofanalysistothecalculatednucleondensitydistributionsbecauseofthestrongerspin-orbitcoupleinheavynuclei.TheresultsindicatethattheRMFtheoryiswellsuitedforthestudyofclusteringinheavynuclei.KeywordsRMF,clusteringinheavynuclei
1)E-mail:[email protected]