中等质量核形状共存与壳效应研究

上一篇

下一篇

刘冬, 郭建友. 中等质量核形状共存与壳效应研究[J]. 物理学报, 2025, 74(11): 112102-1. doi: 10.7498/aps.74.20250095
引用本文: 刘冬, 郭建友. 中等质量核形状共存与壳效应研究[J]. 物理学报, 2025, 74(11): 112102-1. doi: 10.7498/aps.74.20250095
Dong LIU, Jianyou GUO. Shape coexistence and shell effect of medium mass nuclei[J]. Acta Physica Sinica, 2025, 74(11): 112102-1. doi: 10.7498/aps.74.20250095
Citation: Dong LIU, Jianyou GUO. Shape coexistence and shell effect of medium mass nuclei[J]. Acta Physica Sinica, 2025, 74(11): 112102-1. doi: 10.7498/aps.74.20250095

中等质量核形状共存与壳效应研究

    作者简介: 刘冬. E-mail: 1204097178@qq.com .
    通讯作者: E-mail: jianyou@ahu.edu.cn.
  • 中图分类号: 21.10.Ft, 21.10.Ky, 21.10.Hw

Shape coexistence and shell effect of medium mass nuclei

    Corresponding author: E-mail: jianyou@ahu.edu.cn.
  • MSC: 21.10.Ft, 21.10.Ky, 21.10.Hw

  • 摘要: 对原子核形状共存和壳效应的研究有助于人们深入理解原子核内部结构. 物理学家们在Zn, Ge, Se, Kr的同位素研究中, 发现了显著的形状共存现象与刚性三轴性特征. 为了深入探究形状共存现象及其对原子核基态性质的影响, 我们采用相对论Hartree-Bogoliubov理论中密度依赖的介子交换模型, 对N = 32—42的偶偶核Zn, Ge, Se, Kr同位素的基态性质进行了系统研究, 获得的势能面清晰地展现了这些同位素存在形状共存和三轴性特征. 计算获得了原子核的基态能量、形变参数、双中子分离能、中子半径、质子半径和电荷半径, 结果都支持N = 40为新幻数, 部分结果也支持N = 32, 34为新幻数. 尤其, 三轴形变在其中扮演着重要角色. 进一步, 我们探讨了壳效应与形状共存现象之间可能存在的关联及其对原子核基态性质的影响, 并分析了这些变化的物理机制.
  • 加载中
  • 图 1  采用DD-ME2有效相互作用, 通过约束四极变形的三轴RHB计算生成中子数N = 32, N = 34, N = 40的Zn, Ge, Se, Kr同位素的势能曲面. 所有势能面经过归一化处理, 最小值代表基态能量. 等高线由能量相同点连成, 相邻等高线之间的能差为0.6 MeV

    Figure 1.  The potential energy surfaces of Zn, Ge, Se, Kr isotopes, with neutron numbers N = 32, N = 34, and N = 40. These surfaces are generated through triaxial RHB calculations with constrained quadrupole deformation, employing the DD-ME2 effective interactions. The scale of the potential energy is consistent across all surfaces, and the lowest minimum represents the ground state. The contours join points on the surface with the same energy, and the separation between neighboring contours is 0.6 MeV.

    图 2  采用DD-ME2有效相互作用, 通过约束四极变形的三轴RHB计算生成中子数N = 32—42的Zn, Ge, Se, Kr偶偶核同位素的结合能与比结合能, 并与参考文献[38,39]中的实验值进行对比

    Figure 2.  The binding energy per nucleon and the total binding energy of N = 32–42 Zn, Ge, Se, Kr isotopes in the triaxial RHB calculations using both the DD-ME2 interactions, with comparisons made to experimental data from Ref. [38,39].

    图 3  采用DD-ME2有效相互作用, 通过约束四极变形的三轴RHB计算生成中子数N = 32—42的Zn, Ge, Se, Kr偶偶核同位素的双中子分离能, 并与参考文献[38,39]中的实验值进行对比

    Figure 3.  The two-neutron separation energy for N = 32–42 Zn, Ge, Se, Kr isotopes, obtained through triaxial RHB calculations using the DD-ME2 interactions, with comparisons made to experimental data from Ref. [38,39].

    图 4  采用DD-ME2有效相互作用, 通过约束四极变形的三轴RHB计算生成中子数N = 32—42的Zn, Ge, Se, Kr偶偶核同位素的中子、质子半径与电荷半径, 并与参考文献[38,39]中电荷半径的实验值进行对比

    Figure 4.  The neutron, proton, and charge radii for N = 32–42 Zn, Ge, Se, Kr isotopes in the triaxial RHB calculations using the DD-ME2 interactions, with comparisons made to charge radii experimental data from Ref. [38,39].

    表 1  势能曲面极小值点, 坐标(β, γ)为极小值点, 第一极小值能量最低, 对应基态位置

    Table 1.  The minima of the potential energy surfaces. The coordinates (β, γ) specify the locations of these minima. The primary minimum, which is the deepest, corresponds to the ground state of the nucleus.

    原子核 第一极小值 第二极小值 实验四极形变值β
    62Zn (0.25, 20) 0.216
    64Zn (0.23, 0) (0.24, 60) 0.236
    70Zn (0, 0) 0.216
    64Ge (0.26, 27) 0.259
    66Ge (0.25, 60) (0.24, 0) 0.172
    72Ge (0, 0) (0.21, 60) 0.240
    66Se (0.26, 60) (0.25, 0)
    68Se (0.27, 60) (0.25, 0) 0.242
    74Se (0.22, 60) (0, 0) 0.302
    68Kr (0.27, 38)
    70Kr (0.30, 60) (0.25, 0)
    76Kr (0, 0) (0.19, 60) 0.290
    下载: 导出CSV
  • [1] Singh P, Korten W, Hagen T W, Gorgen A, Grente L, Salsac M D, Farget F, Clément E, de France G, Braunroth T, Bruyneel B, Celikovic I, Delaune O, Dewald A, Dijon A 2018 Phys. Rev. Lett. 121 192501 doi: 10.1103/PhysRevLett.121.192501
    [2] Abusara H, Ahmad S 2017 Phys. Rev. C 96 064303
    [3] Cejnar P, Jolie J, Casten R F 2010 Rev. Mod. Phys. 82 2155 doi: 10.1103/RevModPhys.82.2155
    [4] Norman E B, Drobizhev A, Gharibyan N, Gregorich K E, Kolomensky Yu G, Sammis B N, Scielzo N D, Shusterman J A, Thomas K J 2024 Phys. Rev. C 109 055501 doi: 10.1103/PhysRevC.109.055501
    [5] Majola S N T, Shi Z, Song B Y, Li Z P, Zhang S Q, Bark R A, Sharpey-Schafer J F, Aschman D G, Bvumbi S P, Bucher T D, Cullen D M, Dinoko T S, Easton J E, Erasmus N, Greenlees P T 2019 Phys. Rev. C 100 044324 doi: 10.1103/PhysRevC.100.044324
    [6] Yang Y L, Zhao P W, Li Z P 2023 Phys. Rev. C 107 024308 doi: 10.1103/PhysRevC.107.024308
    [7] Hua H, Wu C Y, Cline D, Hayes A B, Teng R, Clark R M, Fallon P, Goergen A, Macchiavelli A O, Vetter K 2004 Phys. Rev. C 69 014317 doi: 10.1103/PhysRevC.69.014317
    [8] Cwiok S, Heenen P H, Nazarewicz W 2005 Nature 433 705 doi: 10.1038/nature03336
    [9] Ayangeakaa A D, Janssens R V F, Wu C Y, Allmond J M, Wood J L, Zhu S, Albers M, Almaraz-Calderon S, Bucher B, Carpenter M P, Chiara C J, Cline D, Crawford H L, Harker J, Hayes A B, Hoffman C R, Kay B P, Kolos K, Korichi A 2016 Phys. Lett. B 754 254 doi: 10.1016/j.physletb.2016.01.036
    [10] 圣宗强, 郭建友 2008 物理学报 57 1557 doi: 10.7498/aps.57.1557 Sheng Z Q, Guo J Y 2008 Acta Phys. Sin. 57 1557 doi: 10.7498/aps.57.1557
    [11] 焦朋, 郭建友, 方向正 2010 物理学报 59 2369 doi: 10.7498/aps.59.2369 Jiao P, Guo J Y, Fang X Z 2010 Acta Phys. Sin. 59 2369 doi: 10.7498/aps.59.2369
    [12] 王刚, 方向正, 郭建友 2012 物理学报 61 102101 doi: 10.7498/aps.61.102101 Wang G, Fang X Z, Guo J Y 2012 Acta Phys. Sin. 61 102101 doi: 10.7498/aps.61.102101
    [13] Nomura K, Rodriguez-Guzman R, Robledo L M 2016 Phys. Rev. C 94 044314 doi: 10.1103/PhysRevC.94.044314
    [14] Karim A, Siddiqui T A, Ahmad S 2022 Phys. At. Nucl. 85 588 doi: 10.1134/S1063778823010246
    [15] Zhang X Y, Niu Z M, Sun W, Xia X W 2023 Phys. Rev. C 108 024310
    [16] Garcia-Ramos J E, Arias J M, Dukelsky J 2014 Phys. Lett. B 736 333 doi: 10.1016/j.physletb.2014.07.038
    [17] 童红, 张春梅, 石筑一, 汪红, 倪绍勇 2010 物理学报 59 3136 doi: 10.7498/aps.59.3136 Tong H, Zhang C, Shi Z Y, Wang H, Ni S Y 2010 Acta Phys. Sin. 59 3136 doi: 10.7498/aps.59.3136
    [18] Bonatsos D, Assimakis I E, Minkov N, Martinou A 2017 Phys. Rev. C 95 064326 doi: 10.1103/PhysRevC.95.064326
    [19] 支启军 2011 物理学报 60 052101 doi: 10.7498/aps.60.052101 Zhi Q J 2011 Acta Phys. Sin. 60 052101 doi: 10.7498/aps.60.052101
    [20] Wu X H, Ren Z X, Zhao P W 2022 Phys. Rev. C 105 L031303
    [21] Gupta S, Bakshi R, Gupta S, Singh S, Bharti A, Bhat G H, Sheikh J A 2023 Eur. Phys. J. A 59 258 doi: 10.1140/epja/s10050-023-01166-6
    [22] 卢希庭, 江栋兴, 叶沿林 2000 原子核物理(北京: 原子能出版社)第192页 Lu X T, Jiang D X, Ye Y L 2000 Nuclear Physics (Beijing: Atomic Energy Press) p192
    [23] Wienholtz F, Beck D, Blaum K, Borgmann C, Breitenfeldt M, Cakirli R B, George S, Herfurth F, Holt J D, Kowalska M, Kreim S, Lunney D, Manea V, Menendez J, Neidherr D, Rosenbusch M, Schweikhard L, Schwenk A, Simonis J, Stanja J, Wolf R N, Zuber K 2013 Nature 498 346 doi: 10.1038/nature12226
    [24] Steppenbeck D, Takeuchi S, Aoi N, Doornenbal P, Matsushita M, Wang H, Baba H, Fukuda N, Go S, Honma M, Lee J, Matsui K, Michimasa S, Motobayashi T, Nishimura D, Otsuka T, Sakurai H, Shiga Y, Soderstrom P A, Sumikama T, Suzuki H, Taniuchi R, Utsuno Y, Valiente-Dobon J J, Yoneda K 2013 Nature 502 207 doi: 10.1038/nature12522
    [25] Michimasa S, Kobayashi M, Kiyokawa Y, Ota S, Ahn D S, Baba H, Berg G P A, Dozono M, Fukuda N, Furuno T, Ideguchi E, Inabe N, Kawabata T, Kawase S, Kisamori K, Kobayashi K, Kubo T, Kubota Y, Lee C S, Matsushita M, Miya H, Mizukami A, Nagakura H, Nishimura D, Oikawa H, Sakai H, Shimizu Y, Stolz A, Suzuki H, Takaki M, Takeda H, Takeuchi S, Tokieda H, Uesaka T, Yako K, Yamaguchi Y, Yanagisawa Y, Yokoyama R, Yoshida K, Shimoura S 2018 Phys. Rev. Lett. 121 022506 doi: 10.1103/PhysRevLett.121.022506
    [26] Liu J, Niu Y F, Long W H 2020 Phys. Lett. B 806 135524 doi: 10.1016/j.physletb.2020.135524
    [27] Zhang W, Huang J K, Sun T T, Peng J, Zhang S Q 2024 Chin. Phys. C 48 104105 doi: 10.1088/1674-1137/ad62dd
    [28] Ring P 1996 Prog. Part. Nucl. Phys. 37 193 doi: 10.1016/0146-6410(96)00054-3
    [29] Vretenar D, Afanasjev A V, Lalazissis G A, Ring P 2005 Phys. Rep. 409 101 doi: 10.1016/j.physrep.2004.10.001
    [30] Meng J, Toki H, Zhou S G, Zhang S Q, Long W H, Geng L S 2006 Prog. Part. Nucl. Phys. 57 470 doi: 10.1016/j.ppnp.2005.06.001
    [31] Liang H Z, Meng J, Zhou S G 2015 Phys. Rep. 570 1 doi: 10.1016/j.physrep.2014.12.005
    [32] Niksic T, Paar N, Vretenar D, Ring P 2014 Comput. Phys. Commun. 185 1808 doi: 10.1016/j.cpc.2014.02.027
    [33] Ring P, Schuck P 1981 Phys. Today 36 70
    [34] Staszack A, Stoitsov M, Baran A, Nazarewicz W 2010 Eur. Phys. J. A 46 85 doi: 10.1140/epja/i2010-11018-9
    [35] Tian Y, Ma Z Y, Ring P 2009 Phys. Lett. B 676 44 doi: 10.1016/j.physletb.2009.04.067
    [36] Niksic T, Ring P, Vretenar D, Tian Y, Ma Z Y 2010 Phys. Rev. C 81 054318 doi: 10.1103/PhysRevC.81.054318
    [37] 沈水法, 王华磊, 孟海燕, 阎玉鹏, 沈洁洁, 王飞鹏, 蒋海滨, 包莉娜 2021 物理学报 70 192101 doi: 10.7498/aps.70.20210187 Shen S F, Wang H L, Meng H Y, Yan Y P, Shen J J, Wang F P, Jiang H B, Bao L N 2021 Acta Phys. Sin. 70 192101 doi: 10.7498/aps.70.20210187
    [38] Wang M, Huang W J, Kondev F G, Audi G, Naimi S 2021 Chin. Phys. C 45 030003 doi: 10.1088/1674-1137/abddaf
    [39] Wang S J, Kanellakopoulos A, Yang X F, Bai S W, Billowes J, Bissell M L, Blaum K, Cheal B, Devlin C S, Garcia Ruiz R F, Han J Z, Heylen H, Kaufmann S, König K, Koszorús Á, Lechner S, Malbrunot-Ettenauer S, Nazarewicz W, Neugart R, Neyens G, Nörtershäuser W, Ratajczyk T, Reinhard P G, Rodríguez L V, Sels S, Xie L, Xu Z Y, Yordanov D T, Yu Y M 2024 Phys. Lett. B 856 138867 doi: 10.1016/j.physletb.2024.138867
    [40] El Adri M, Oulne M 2020 Int. J. Mod. Phys. E 29 2050089
    [41] 陈翠红, 李占奎, 王秀华, 李荣华, 方芳, 王柱生, 李海霞 2023 物理学报 72 122902 doi: 10.7498/aps.72.20230213 Chen C H, Li Z K, Wang X H, Li R H, Fang F, Wang Z S, Li H X 2023 Acta Phys. Sin. 72 122902 doi: 10.7498/aps.72.20230213
    [42] Liu Y, Wang R, Mushtaq Z, Tian Y , He X H, Qiu H, Chen X R 2025 Chin. Phys. C 49 034103
    [43] Enciu M, Liu H N, Obertelli A, Doornenbal P, Nowacki F, Ogata K, Poves A, Yoshida K, Achouri N L 2022 Phys. Rev. Lett. 129 262501 doi: 10.1103/PhysRevLett.129.262501
  • 加载中
图( 4) 表( 1)
计量
  • 文章访问数:  23
  • HTML全文浏览数:  23
  • PDF下载数:  1
  • 施引文献:  0
出版历程
  • 收稿日期:  2025-01-21
  • 刊出日期:  2025-06-05

中等质量核形状共存与壳效应研究

    通讯作者: E-mail: jianyou@ahu.edu.cn.
    作者简介: 刘冬. E-mail: 1204097178@qq.com
  • 安徽大学物理与光电工程学院, 合肥 230601

摘要: 对原子核形状共存和壳效应的研究有助于人们深入理解原子核内部结构. 物理学家们在Zn, Ge, Se, Kr的同位素研究中, 发现了显著的形状共存现象与刚性三轴性特征. 为了深入探究形状共存现象及其对原子核基态性质的影响, 我们采用相对论Hartree-Bogoliubov理论中密度依赖的介子交换模型, 对N = 32—42的偶偶核Zn, Ge, Se, Kr同位素的基态性质进行了系统研究, 获得的势能面清晰地展现了这些同位素存在形状共存和三轴性特征. 计算获得了原子核的基态能量、形变参数、双中子分离能、中子半径、质子半径和电荷半径, 结果都支持N = 40为新幻数, 部分结果也支持N = 32, 34为新幻数. 尤其, 三轴形变在其中扮演着重要角色. 进一步, 我们探讨了壳效应与形状共存现象之间可能存在的关联及其对原子核基态性质的影响, 并分析了这些变化的物理机制.

English Abstract

    • 近年来, 核物理学界对原子核基态的形状共存与壳效应展现出了浓厚的兴趣, 这些奇特现象推动了对原子核复杂结构的实验研究. 实验中观测到的原子核展现出多样化的形状特征[1], 不仅丰富了我们对原子核多样性的认识, 同时也对当前各种理论框架构成了重大挑战, 要求理论模型能够更为精确和全面地描述这些现象[25]. 原子核是由核子组成的极其复杂的量子多体系统, 在内禀系中, 由于对称性的自发破缺, 原子核可以有各种形状, 大多数原子核的基态具有轴对称的四极形变, 表现为长椭球形或扁椭球形. 然而, 也存在一些原子核的形态更为复杂, 并不遵循这种轴对称的规律. 物理学家们为了解决这一问题, 引入了三轴形变的概念[6,7]. 同时, 原子核轴对称性被打破, 也是导致形状共存现象的产生的原因之一[8,9].

      为了深入探究这些现象的本质, 核物理学家已发展出许多理论方法. 2008年, 圣宗强等[10]、焦朋等[11]和王刚等[12]在相对论平均场理论下用NL3参数组系统地研究了多个同位素链中的偶偶核, 并在后续对原子核形状演变做了更进一步的探讨. 2016年, Nomura等[13]基于Gogny-D1 M能量密度泛函的自洽平均场近似所提供的微观数据, 运用相互作用玻色子模型, 对质量数A = 100附近的原子核的形状相变进行了详细研究, 发现了Zr和Sr原子核中存在着显著的长椭球形和扁椭球形共存的现象. 2022年, Karim等[14]在协变密度泛函数理论(CDFT)的框架下, 利用密度依赖的介子交换和点耦合模型, 系统地研究了I同位素(Z = 53, N = 70—80)的基态性质, 结果清楚地展现了原子核的形状相变和形状共存现象. 2023年, Zhang等[15]采用连续谱中的形变相对论Hartree-Bogoliubov (RHB)理论研究了Kr和Sr同位素的基态性质. 同时也预测了在这些同位素链中可能存在形状共存现象. 此外, 集体哈密顿量以及基于机器学习的密度泛函理论等多种方法也被用于研究原子核的形状共存[1621].

      原子核的形状主要由其核子数以及核子之间的相互作用决定. 原子核中的质子和中子两两配对, 当质子或中子壳层被完全占满时, 即达到所谓的“满壳”状态, 此时原子核具有相对稳定的能量状态, 核子之间的相互作用力达到了一种平衡, 原子核的基态通常呈现为球形[22]. 当添加或减少核子时, 这种平衡状态可能会被打破, 从而导致原子核发生形变, 形成其他非球形形状, 例如长椭球形或扁椭球形等. 在研究不同原子核沿着同位素链的形状相变时, 我们通过理论计算发现当中子数为传统幻数时, 例如N = 50, 原子核的形状呈现为标准的球形. 当中子数为新幻数时, 例如N = 32, 34, 40[2326], 原子核通常会出现形状共存与三轴性特征. 值得注意的是, 幻数的稳定性可能显著依赖于核区: 近期研究表明, 在超镄核区(Z = 102—110, N≈258), 由于极端质子-中子比与强库仑排斥效应, 可能存在新的壳层闭合现象[27]. 这种幻数随核区变化的特性进一步凸显了核结构对形变、连续谱效应等微观机制的敏感性. 这让我们对新幻数的研究产生了兴趣. 在本文中, 我们基于协变密度泛函理论, 深入探究了中子数包含新幻数的Zn, Ge, Se, Kr同位素的基态性质, 此外还探究了这些同位素中可能存在的形状共存与三轴性. 通过采用具有密度依赖的介子交换的RHB理论, 我们计算了这些同位素的势能面. 在此基础上我们还获得了这些同位素的基态形变值、结合能、双中子分离能、质子半径、中子半径以及电荷半径, 计算的结果都支持N = 40为新幻数, 部分结果也支持N = 32, 34为新幻数. 该方法实现了对原子核的基态性质的自洽描述, 为理解原子核复杂结构提供了新的视角与见解.

    • CDFT在描述稳定核和奇特核的基态性质方面非常有效. 它对原子核形状演化及相变和多种形状共存等现象提供了特殊的视角. 在我们目前的研究中, 采用了基于介子交换相互作用的CDFT, 详细内容可见对此方面有开创性贡献的文献[2831]. 在基于介子交换相互作用的CDFT的这一理论框架下, 我们选用了DD-ME2参数集进行RHB计算. 在下面, 我们将对理论框架进行概述. 在采用DD-ME2参数组进行RHB计算时, 该体系的拉格朗日密度表示为

      其中$ \psi $表示狄拉克旋量, $ M $表示核子质量. 符号$ m_\sigma ^{} $, $ m_\omega ^{} $$ m_\rho ^{} $分别表示$ \sigma $, $ \omega $$ \rho $介子的质量. 相应的耦合常数分别为$ {g_\sigma } $, $ {g_\omega } $$ {g_\rho } $. $ e $与质子的电荷有关. 场张量$ {\varOmega _{\mu \nu }} $, $ {\vec R^{\mu \nu }} $$ {F_{\mu \nu }} $分别表示为

      对于满足时间反演对称性的静态原子核而言, 其内部不产生流. 因此, 空间分量$ \omega $, $ \vec \rho $$ A $消失了. 剩下的是类时间分量$ {\omega _0} $, $ {\vec \rho _0} $$ {A_0} $. 由同位旋守恒可得只有$ {\rho _3} $分量存在. 经典变分原理可以给出核子运动的狄拉克方程:

      以及介子和光子运动的克莱因-高登方程:

      其中相互作用势分为标量势与矢量势:

      其中

      $ {\rho _{tv}} $代表同位旋矢量密度, 即质子和中子矢量密度之差. 更多细节可以参考文献[31].

    • 本研究基于上述理论框架开展计算, 采用了DD-ME2参数组与可分离对力. 计算使用代码可在Niksic等[32]论文中下载. 在讨论结果之前, 需要先解释数值方面的问题. 本文利用三轴相对论Hartree-Bogoliubov(RHB)模型, 对总能量进行约束计算, 这个过程使我们能够准确定位全局基态最小值并跟踪基态形状的演变. 使用三维谐振子基求解RHB方程. 对于玻色子, 壳层的数量固定为NB = 20, 费米子的壳层数则固定为Nf = 12. 约束计算是通过对轴向和三轴质量四极矩施加约束来实现的, 使用二次约束的方法[33]使用函数的变化不受限制.

      式中, $ \langle \hat H\rangle $为总能量; $ \langle {{{\hat Q}_{2\mu }}} \rangle $为四极算符的期望值.

      $ {q_{2\mu }} $是多极矩的约束值, $ {C_{2\mu }} $是相应的刚度常数. 二次约束为系统增加了一个额外力项$\displaystyle \sum\nolimits_{\mu = 0, 2} {{\lambda _\mu }} {Q_{2\mu }} $, 其中$ {\lambda _\mu } = 2{C_{2\mu }}( {\langle {{{\hat Q}_{2\mu }}} \rangle - {q_{2\mu }}} ). $此外, 四极矩$ {\hat Q_{2\mu }} $与约束值$ {q_{2\mu }} $之间的差取决于刚度常数的值, 即$ {C_{2\mu }} $的值越小, 四极矩与相应的约束值的偏差越大. 然而, 增大刚度常数常常会导致自洽计算无法收敛. 我们使用增广拉格朗日方法[34]确保了自洽计算过程的收敛性. 在数值计算的过程中, β范围从0到0.6, 步长为0.05, γ从0°到60°, 步长为5. 我们使用了Tian等[35]提出的在坐标空间中可分离的有限范围配对相互作用, 由下式给出:

      式中, $ {\boldsymbol{R }}= {1}/{2}\left( {{{\boldsymbol{r}}_1} + {{\boldsymbol{r}}_2}} \right) $为质量中心; $ {\boldsymbol{r}} = \left( {{{\boldsymbol{r}}_1} - {{\boldsymbol{r}}_2}} \right) $为相对坐标; $ P(r) $为形状因子, 由下式给出:

      式中, $ G $为配对强度; $ a $为配对宽度. 具体数值$ G $ = 728 MeV·fm3, $ a $ = 0.644 fm. 此数据通过Gogny力D1S参数确定[35], 更加具体的推导过程可在Niksic等[36]论文中查看. 在此基础上, 我们计算了基态形变下的几种物理量, 其中包括比结合能, 双中子分离能, 中子半径, 质子半径以及电荷半径.

      图1中可发现, N = 40时, 70Zn, 72Ge, 74Se, 76Kr都算出了球形极小值, 这与我们的预期相同. 当壳层闭合时, 原子核的能量最低, 也最稳定. 核子在空间中的分布趋向于均匀, 从而形成球形核. 但是除了70Zn外, 所有核都出现了形状共存现象或者三轴性. N = 32时, 62Zn, 64Ge, 68Kr都出现了明显的三轴形变, 这与其他理论模型预测的结果一致(例如, 沈水法等[37]通过自洽推转壳模型计算表明, 64Ge在基态和集体转动态下可能存在γ = –30°的三轴形变), 66Se也出现了形状共存现象. N = 34时, 64Zn, 66Ge, 68Se, 70Kr均出现了形状共存现象. 具体形变值可以查看表1. 我们计算得出的形变值与传统的幻数核之间有明显的差异. 例如N = 50时, 80Zn, 82Ge, 84Se, 8 6Kr均为标准球形核, 其形变值均为0. 当然, 原子核的形变受到核子数, 核子间相互作用等多种因素的影响, 实验测量中也没有发现完美的球形核, 理论计算出原子核的形变值不为零, 并不能直接推断出壳层一定不闭合. 因此我们可以尝试通过其他物理量, 例如比结合能, 双中子分离能等来间接分析这个问题.

      图2展示了Zn, Ge, Se, Kr四条同位素链的结合能与比结合能. 图2中清晰地展示了原子核的结合能随着中子数的增加而增加, 实验数据与理论计算符合得很好. 而在N = 40时, 70Zn, 72Ge的比结合能均达到最大值, 而74Se, 76Kr只是接近最大值, 这可能与库仑斥力有关, 质子之间存在库仑斥力, 这种斥力会随着质子数的增加而增大. 当质子数过多时, 库仑斥力会削弱核子之间的相互作用, 从而影响比结合能. 也可能与原子核的形变有关, 当N = 40时, 70Zn, 72Ge计算所得到的第一极小值都是球形, 而74Se的第一极小值是一个明显的椭球形, 这种形变会影响核子之间的相互作用与原子核的稳定性, 从而导致比结合能的变化. 但是我们认为影响最大的原因是质子数的区别. 除了中子数外, 质子数对原子核的稳定性也有着很大的影响. 当质子数和中子数都满足特定条件时, 原子核才会达到最大的比结合能, 这点在70Zn上体现得尤为明显.

      双中子分离能被认为是研究原子核壳层结构的重要物理量. 它指的是从原子核内部分离出两个中子所需要的能量, 其计算公式为

      其中EB(Z, N)表示具有Z个质子和N个中子的原子核的结合能. 我们知道在传统幻数中, 如N = 50时, 双中子分离能会出现明显的下降. 这是因为幻数核对应的核子填充了特定的壳层, 形成了满壳层或接近满壳层的结构. 再增加中子会导致这个额外的中子无法填充到已有的满壳层中, 而是需要进入下一个更高的能级. 这个额外的中子与原子核的相互作用相对较弱, 因此移除这两个中子所需的能量会突然减小, 即双中子分离能快速下降. 在El Adri和Oulne[40]的研究中, 就展示了Ge, Se, Kr, Sr同位素链在N = 50与N = 70时, 双中子分离能所产生的明显的下降现象. 同样, 图3中展示了Zn, Ge, Se, Kr四条同位素链的双中子分离能随着中子数N的变化而减小. 在N = 32时, 62Zn与64Ge理论预测的双中子分离能都出现了明显的下降; N = 34时, 68Se与70Kr的双中子分离能都出现了明显的下降; N = 40时, 70Zn, 72Ge, 76Kr的双中子分离能都出现了明显的下降. 值得注意的是, 导致双中子分离能突变的原因不只是新幻数, 也可能与形变有关, 比如出现形状共存现象时, 长椭球形和扁椭球形的竞争和混合也可能会导致分离能减小, 这一点在74Se体现得尤为明显. 在74Se中存在较强的球形核与长椭球形的竞争, 如果考虑第二极小值, 74Se的双中子分离能也会发生突变. 通过与实验数据进行对比, 可以观察到在部分接近质子滴线的核素中, 实验数据的趋势与理论预测存在一定偏离. 这种差异可能源于两方面原因, 一是接近质子滴线的核素测量难度大, 实验存在明显的误差; 二是理论与实验所预测的形变不同. 同时, 大部分同位素的理论计算与实验数据表现出良好的一致性.

      同样地, 原子核中子、质子半径以及电荷半径也同样常被用来揭示原子核的壳层结构. 具体来说, 当核内的质子和中子数接近幻数时, 质子和中子的分布更加紧凑, 半径较小, 反之, 未闭合的壳层会导致较大的半径和较不规则的形状. 电荷半径能够反映出原子核的整体大小, 通常也与壳层结构紧密相关. 由于质子带正电, 所以电荷半径是通过质子半径计算得来: $ {R}_{\text{c}} = \sqrt{{R}_{\text{p}}^{2}+0.64}, $ 其中0.64为体积修正, 通过质子体积得出. 在图4(a)中, 可以观察到N = 34时, 70Kr的中子半径出现了明显的突变, 在图4(b)中与图4(c)中, 观察到N = 34时, 64Zn, 68Se的电荷半径均发生了突变, N = 40时, 70Zn, 72Ge, 74Se, 76Kr的电荷半径都处于一个转折的位置, 这个结果与前面计算出的形变值、比结合能与双中子分离能都有很好的对应. 在图4(c)中加入了电荷半径的实验值, 可以观察到Zn, Ge, Se同位素链的理论预测与实验数据整体符合良好, 但Kr同位素链的电荷半径变化趋势与实验观测存在显著差异, 导致这一现象的原因可能是当前RHB理论高估了N = 40中子闭壳效应对质子分布的刚性约束. 尽管计算显示该核素维持球形基态(β≈0), 但本文引用的实验数据在此处仍有明显的β形变值. 从而使实验测量的电荷半径(Rc) 大于理论值. 当然, 现有实验数据的稀疏性进一步限制了趋势分析的可靠性: 在N = 34—40的区间内, 已发表的Kr同位素电荷半径仅有3个数据点, 并且靠近质子滴线导致电荷半径测量依赖间接方法, 这会引入约0.05 fm的系统误差. 可能要等待新一代放射性束流装置(如FRIB, HIAF)的高精度数据来进一步分析该问题[41,42].

    • 本文采用具有密度依赖的介子交换相互作用的RHB理论系统地研究了Zn, Ge, Se, Kr(32<N<42)偶偶核同位素的基态性质. 计算获得的势能面清楚地揭示了原子核的三轴性与形状共存现象. 通过精确定位势能面上最小值的位置, 得到了原子核基态的结合能、形变参数、双中子分离能, 以及中子、质子、电荷半径等物理性质, 讨论了形状共存与三轴性对这些原子核基态性质的影响, 以及这些原子核基态性质与壳层闭合之间存在的关联. 结果表明, 新幻数N = 32, 34, 40在一些原子核基态性质上展现出较强的稳定性. 值得注意的是, N = 32, 34的潜在幻数性质呈现出显著的核素依赖性, 例如: 在N = 32时, 62Zn和64Ge的双中子分离能出现了明显的下降趋势, 但相同中子数的64Se未表现出类似突变; 在N = 34时, 68Se和70Kr的双中子分离能显著下降, 而70Kr的电荷半径同时出现反常增大. 这种幻数效应随核区变化的现象与近期实验观测一致[23,24,43], 表明幻数的稳定性可能受形变、连续谱效应等核结构因素的调制[26]. 在N = 40时, 计算得到的形变参数、比结合能、双中子分离能、电荷半径所表现出的特征都支持N = 40为新幻数这一结论. 然而, 我们发现一些原子核在中子数N = 32, 34, 40处具有三轴性并且存在形状共存现象, 这与传统幻数核的球形闭壳特征存在明显差异. 未来工作可结合张量力修正、大尺度壳模型计算及新一代放射性束流装置(如FRIB, HIAF)的高精度数据, 进一步厘清核力参数化、质子-中子平衡与演生对称性的关联, 为极端条件下核结构的研究提供更完备的理论框架.

    参考文献 (43)

目录

/

返回文章
返回