-
摘要: 通过重离子熔合蒸发反应82Se(18O,p4n)95Nb,布居了95Nb的高自旋态。利用多探头中国联合伽马探测阵列进行了在束γ测量,研究了95Nb核的高自旋态能级结构。基于γ-γ符合关系扩展了95Nb的能级纲图,发现了25条新的γ跃迁和21条新的能级,并调整了部分能级的摆放位置。将新的能级结构与壳模型计算进行了比较和讨论。结果表明,要正确的描述95Nb核的高自旋态(自旋值大于39/2)能级结构,应考虑质子跨越38子壳 f5/2→g9/2核心激发。Abstract: High-spin level structure of 95Nb has been investigated using the multi-detector array of the conjoint gamma array in China via the 82Se(18O, p4n)95Nb. Based on γ-γ coincidence relationships, the level scheme of 95Nb has been modified and extended with 25 new γ rays and 15 new levels. The new level structure of 95Nb has been compared with the shell model calculations. It is suggested that the proton core excitation(f5/2→g9/2) across $Z = 38 $ sub-shell should be taken into account in order to adequately describe the high-spin level structure in 95Nb above spin of 39/2.
-
Keywords:
- high-spin states /
- level scheme /
- linear polarization /
- shell model
-
0. 引言
由于Z=38 亚壳与N = 50幻数壳的存在,A~90,N~50核区的核素大多具有近球形核结构,低激发态主要表现出单粒子运动的特性[1]。对于N = 50和56之间的核结构,通常理解为这些壳层之上核子的单粒子激发[2-5]。随着中子数的增加,A≈100的原子核一般会经历从球形向形变核的转变,当中子数超过56时,这个质量区域的原子核形状会发生快速的形状转变,价核子数目的增多和激发能的增加更可能会使原子核的高自旋态出现集体运动和核芯激发现象。
通过本组之前对A~90质量区Nb同位素的研究发现,对于N>50的稳定线以及稳定线附近的92, 93Nb核[4-5],中等自旋能级可由Z = 38的核芯激发或N = 56的亚壳层激发形成,更高自旋的能级会涉及N = 50的核芯激发。然而,由于实验的困难,稳定线右侧的丰中子Nb同位素高自旋态信息非常缺乏。Bucurescu等[1]利用82Se(16O, 1p2n)95Nb反应,16O束流能量为48 MeV,布局了95Nb的高自旋态,由于入射粒子能量较低,所得纲图信息较少,并未对高自旋能级进行指定。
因此,进一步扩展丰中子侧Nb同位素能级纲图,揭示其单粒子激发和集体运动的共存和竞争,是很有意义的。本文将详细介绍95Nb高自旋态的实验研究结果,在扩展其能级纲图的基础上,通过把实验结果和壳模型计算结果相比较,探讨95Nb的激发机制。
1. 实验和数据分析
在HI-13串列加速器上,利用中国联合伽马探测阵列,最多可放置40个Ge探测器,BGO和HPGe探测器分别与法兰相连,径向距离可通过螺杆调节,如图1所示。采用82Se(18O, 1p4n)95Nb熔合蒸发反应,18O束流能量为82和88 MeV,对95Nb的能级结构进行研究,平均束流强度为90 nA,靶由厚度为0.85 mg/cm2的82Se和4.45 mg/cm2厚的金衬组成。探测系统是由28台带反康的高纯锗(或Clover)探测器和1 台小平面探测器组成。其中6 台探测器放置在与束流夹角为150°的方向,7台探测器放置在60°方向,90°和120°方向各8台。
用60Co,133Ba和152Eu标准放射源对探测器进行了能量和效率刻度。将实验数据离线反演成γ-γ二维对称矩阵,用于分析γ-γ符合关系,提取γ跃迁的强度和建立能级纲图。处理实验数据后得到的95Nb纲图,如图2所示,图中标红射线是在此次实验中新发现的射线,标蓝射线调整了摆放顺序。
为进一步确定新发现能级的自旋,利用90°和150°方向上的探测器构建二维γ-γ DCO矩阵, DCO比例系数为
$$ {{R}}_{{\rm{DCO}}} = \frac{{I_\theta ^{{\gamma _1}}({\text{gated by }} {\gamma _2} {\text{ at }} {{90}^ \circ })}}{{I_{{{90}^ \circ }}^{{\gamma _1}}({\text{gated by }} {\gamma _2} {\text{ at }} \theta )}}\frac{{\varepsilon _\theta ^{{\gamma _1}}\varepsilon _{{{90}^ \circ }}^{{\gamma _2}}}}{{\varepsilon _{{{90}^ \circ }}^{{\gamma _1}}\varepsilon _\theta ^{{\gamma _2}}}} \text{,} $$ (1) 其中:γ2是用来开窗的,已知极性的电四极(ΔI=2)跃迁,与该跃迁存在符合关系的纯四极跃迁的RDCO值在1.0附近,而纯偶极跃迁的值在0.5附近;公式中$I_\theta ^{{\gamma _1}}$和$I_{90^\circ }^{{\gamma _1}}$分别是在150°和90°方向上开窗射线γ2得到的射线γ1的强度;θ为小角度方向150°;ε为90°或150°方向不同能量跃迁的探测器效率。
此外,本工作还利用线极化测量方法确定γ跃迁的电磁特性。利用放置在90°的分辨率较好的Clover探测器测量平行和垂直方向的康普顿散射事件来确定入射光子的极化度(非对称度),把符合事件处理为平行散射事件矩阵和垂直散射事件矩阵。非对称度公式[6]为
$$ A(\theta ) = \frac{{a(E_\gamma ){N_ \bot } - {N_\parallel }}}{{a(E_\gamma ){N_ \bot } + {N_\parallel }}} \text{,} $$ (2) 其中:a($E_\gamma $)为极化效率刻度曲线,是Clover探测器对平行和垂直散射的γ跃迁探测效率的比值;N⊥和$N_\parallel$分别为实验测量的与反应平面垂直和平行的康普顿散射光子的计数。每个Clover探测器可以分为4块晶体,利用90°方向上的Clover探测器(此方向A(θ)有最大值)。实验时伽马射线入射到任意一晶体上,与之相邻的两个晶体探测到的康普顿散射光子被认定分别对应于垂直和平行的康普顿事件。此外,四块晶体探测效率不同,需要对探测器的效率进行修正。作极化探测效率刻度曲线时,使用没有被极化的放射源133Ba和152Eu进行刻度,选取水平和垂直事件进行计算,公式定义为:a($E_\gamma $) = N||/N⊥,将得到的不同能量下的a($E_\gamma $)进行直线拟合。拟合公式为:$a({E_\gamma }) = 0.995\,17 - 2.662\,62{{\rm e}^{ - 5}} \times {E_\gamma }$。A(θ)的计算值若为正,该跃迁为电跃迁,若为负,则为磁跃迁。图3为利用824 keV伽马跃迁开门得到的824.6,844.3和872.2 keV跃迁在垂直(V)和平行方向(P)的能谱,从图中可以看到,824.6,844.3和872.2 keV的水平方向计数明显小于垂直方向计数,为电跃迁。结合之前的DCO测量结果,γ跃迁DCO比值如图4所示,可判断它们三者均为电四极跃迁,以此类推可得到其他γ跃迁的电磁特性。
2. 结果与讨论
2.1 能级纲图
本工作之前,Bucurescu等[1]利用82Se (16O, 1p2n) 95Nb反应,将95Nb自旋和能量推高至33/2 ћ态和8.7 MeV。如图2所示,我们用已知的γ射线开窗寻找新的γ射线,在其基础上发现了约25条新的γ射线、21个新的能级,建立了2个新的序列(序列2和3),自旋推高到了51/2 ћ,激发能约11 MeV;此外,根据符合和强度关系调整了部分能级位置的摆放。图5为此次基态跃迁能量824 keV开窗谱和678.5 keV开窗谱。
在本次实验的678.5 keV开窗谱中,如图5(b)所示,发现了675.8 keV,并且872.2 keV与自身没有级联关系,利用708.1 keV和795.0 keV的级联加和关系和能量强度,修改了678.5 keV的位置在两条824.2 keV之间,并且根据能量强度关系,确定675.8 keV在能级8 692.6 keV之上。在序列1不同γ跃迁的开窗谱中,发现了更高自旋态49/2+之上的1 914.1 keV跃迁。在能级4 069.3 keV之上,分别发现了557.6,511.0,971.9,1 019.7 keV的退激序列。
在本次实验824.2 keV的开窗谱中,如图5(a)所示,可以观测到1 358.9,716.6 keV的正宇称退激序列,通过级联关系、强度关系和级联加和关系确定了716.6 keV往上的两条退激路径。与此同时,在708.1 keV开窗谱中,发现了1 219.7, 1 171.8, 1 183.0和659.3 keV的级联顺序,且它们的级联加和关系和序列2相呼应。这增加了序列3的合理性。对于基态跃迁,708.1和1 358.9 keV只与824.6 keV有级联关系,这能够从实验上确定824.6 keV为基态跃迁。
2.2 壳模型计算
为了解释95Nb的激发能级的实验观测值,我们利用壳模型程序Nushellx[7]对95Nb的激发能级进行计算。选取的壳模型空间为GWB,相互作用势为GWBXG[8]。壳模型空间总共由10 条轨道组成:π(f5/2,p3/2,p1/2,g9/2),ν(p1/2,g9/2,g7/2,d5/2,d3/2,s1/2)。每条轨道所对应的单粒子能量分别为:−5.322,−6.144,−3.941,−1.250,−0.696,−2.597,5.159,1.830,4.261和1.741 MeV。考虑到壳模型计算的复杂性和耗时性,对计算的组态空间进行截断,只将各态的主要组态成分考虑进去,而其他次要的组态成分则忽略。实验结果和壳模型计算结果如图6所示。
95Nb费米能级为πg9/2和νd5/2,相对于Z = 38和N = 50的核心外面有三个价质子和四个价中子。由于质子高jg9/2轨道单粒子激发能较低,而中子高jg7/2轨道的激发能较高,所以95Nb低位能级主要由π(p1/2 g9/2)轨道激发。正如壳模型的计算结果所给出的那样,基态9/2+由质子轨道1g9/2上未成对的质子表示。随着自旋的增加,23/2+及以上开始涉及质子轨道2p1/2质子对向轨道1g9/2激发,例如,31/2+和35/2+态,它们的质子组态都为π(1g9/2)3。当自旋更高时,将涉及中子轨道d5/2上的中子跨56亚壳向更高能量的3s1/2轨道激发,例如,35/2+态,中子组态为ν(1g7/2)1(2d5/2)2(3s1/2)1。在高自旋态时,如39/2+和45/2+态时,涉及质子轨道1f5/2上的质子跨Z = 38壳的核芯激发,质子组态为π(1f5/2)5(2p1/2)1(1g9/2)3。
从图6中可以看出,壳模型计算整体符合较好,对于高自旋能级,实验能级和计算结果相差较大。通常情况下,特定能级的计算值与实验值的偏差程度取决于被忽略的组态成分对这一能级波函数的贡献程度。忽略了特定能级的重要组态成分将会导致这条能级的计算值比实验值大很多。对于本工作中的壳模型计算结果在高激发态出现的明显差异,最主要原因可能是在GWB相互作用空间没有考虑h11/2轨道。因此,进一步进行壳模型计算的意义在于,能够揭示其中子激发与质子激发的竞争和相互作用。
3. 总结
利用熔合蒸发反应82Se (18O, p4n) 95Nb布居了95Nb的高自旋态,扩展了95Nb的能级方案,对95Nb核的高自旋态能级纲图做了重要修改,新发现了25条γ跃迁和21条能级。利用DCO和线极化测量方法,确定了部分能级的自旋和宇称。为了更好地理解95Nb的能级结构,使用Nushellx程序进行了壳模型计算,结果表明,95Nb较低自旋能级组态的主要成分由若干个πg9/2和νd5/2价核子构成,随着自选增加出现π(p1/2→g9/2)激发,在较高自旋能级中,将涉及ν(d5/2→g7/2,s1/2)激发,当自旋大于39/2+和45/2+出现了质子跨Z = 38的核心激发。
致谢 本工作是在中国原子能科学研究院串列加速器上完成的,感谢加速器运维组为实验提供了束流,同时感谢各位老师同学对实验各方面的帮助。
-
[1] BUCURESCU D, PODOLYÁK Zs, RUSU C, et al. Phys Rev C, 2005, 71: 034315. DOI: 10.1103/PhysRevC.71.034315
[2] CUI Xingzhu, ZHU Lihua, WU Xiaoguang, et al. Phys Rev C, 2005, 72: 044322. DOI: 10.1103/PhysRevC.72.044322
[3] LUO Pengwei, WU Xiaoguang, SUN Huibin, et al. Phys Rev C, 2014, 89: 034318. DOI: 10.1103/PhysRevC.89.034318
[4] WU Yiheng, LU Jingbin, LUO Pengwei, et al. Chin Phys Lett, 2014, 34: 042102. DOI: 10.1088/0256-307X/31/4/042102
[5] WU Y H, LU J B, ZHEN R, et al. Phys Rev C, 2022, 105: 034344. DOI: 10.1103/PhysRevC.105.034344
[6] 方永得, 张玉虎, 周小红, 等. 高能物理与核物理, 2007, 31(10): 938. FANG Y D, ZHANG Y D, ZHOU X H, et al. High Energy Phys and Nucl Phys, 2007, 31(10): 938. (in Chinese)
[7] BROWN B A, RAE W D M. Nuclear Data Sheets, 2014, 120: 115. DOI: 10.1016
[8] ZHANG Changhua, WANG Shunjin, GU Jinnan, et al. Phys Rev C, 1999, 60: 054316. DOI: 10.1103/PhysRevC.60.054316