高级检索

留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码
常昶, 杨磊, 林承键, 杨彦云, 温培威, 骆天鹏, 马军兵, 许世伟, 王康, 段芳芳, 马南茹, 贾会明, 杨峰, 黄大湖, 张明昊, 杨过, 杨赟, 莫腾欢. 近垒能区7Be+120Sn的准弹性散射研究[J]. 原子核物理评论, 2023, 40(3): 356-361. doi: 10.11804/NuclPhysRev.40.2023005
引用本文: 常昶, 杨磊, 林承键, 杨彦云, 温培威, 骆天鹏, 马军兵, 许世伟, 王康, 段芳芳, 马南茹, 贾会明, 杨峰, 黄大湖, 张明昊, 杨过, 杨赟, 莫腾欢. 近垒能区7Be+120Sn的准弹性散射研究[J]. 原子核物理评论, 2023, 40(3): 356-361. doi: 10.11804/NuclPhysRev.40.2023005
Chang CHANG, Lei YANG, Chengjian LIN, Yanyun YANG, Peiwei WEN, Tianpeng LUO, Junbing MA, Shiwei XU, Kang WANG, Fangfang DUAN, Nanru MA, Huiming JIA, Feng YANG, Dahu HUANG, Minghao ZHANG, Guo YANG, Yun YANG, Tenghuan MO. Quasielastic Scattering Study for the 7Be+120Sn System at the Energy Near the Coulomb Barrier[J]. Nuclear Physics Review, 2023, 40(3): 356-361. doi: 10.11804/NuclPhysRev.40.2023005
Citation: Chang CHANG, Lei YANG, Chengjian LIN, Yanyun YANG, Peiwei WEN, Tianpeng LUO, Junbing MA, Shiwei XU, Kang WANG, Fangfang DUAN, Nanru MA, Huiming JIA, Feng YANG, Dahu HUANG, Minghao ZHANG, Guo YANG, Yun YANG, Tenghuan MO. Quasielastic Scattering Study for the 7Be+120Sn System at the Energy Near the Coulomb Barrier[J]. Nuclear Physics Review, 2023, 40(3): 356-361. doi: 10.11804/NuclPhysRev.40.2023005

近垒能区7Be+120Sn的准弹性散射研究

doi: 10.11804/NuclPhysRev.40.2023005
基金项目: 国家重点研发计划项目(2018YFA0404404);国家自然科学基金资助项目(12235020, U2167204, 12175313, 12175314, U1732145, U1867212和119611012);稳定支持基础科研计划;中核集团领创项目(LC192209000701, LC202309000201);青年人才发展基金项目(YC2122120001);国防基础科研计划项目(JCKY2020201C002)
详细信息

Quasielastic Scattering Study for the 7Be+120Sn System at the Energy Near the Coulomb Barrier

Funds: National Key R&D Program of China(2018YFA0404404); National Natural Science Foundation of China(12235020, U2167204, 12175313, 12175314, U1732145, U1867212, 11961131012); Continuous-Support Basic Scientific Research Proiect; Leading Innovation Project(LC192209000701, LC202309000201); Young Talent Development Foundation(YC212212000101); Basic Scientific Research Program under grant(JCKY2020201C002)
More Information
  • 摘要: 近垒能区弱束缚核的反应动力学是目前核物理研究热点之一。本工作使用大立体角硅探测器阵列测量了$^{7}{\rm{Be}}$+$^{120}{\rm{Sn}}$体系在48.05 MeV的准弹性散射,结合蒙特卡罗模拟得到其微分截面。基于光学模型,分别用频率方法和贝叶斯方法对角分布进行拟合。在前角区,两种方法给出了一致的结果;在后角区,频率方法的结果呈现明显的振荡结构,而贝叶斯方法的结果较为平滑,在靠近180°有振荡上升的趋势。
  • 图  1  具有大立体角覆盖范围的硅探测器阵列示意图(红线表示束流入射方向)(在线彩图)

    图  2  粒子鉴别图(在线彩图)

    (a) ${\text{TOF-}} E_{\text{r}}$可以将初级束 $ ^{7}{\rm{Li}} $、次级束$ ^{7}{\rm{Be}} $和束流中的$ ^{4}{\rm{He}} $区分; (b) $\varDelta E \text{-}E_{\rm{r}}$二维谱,$ ^{7}{\rm{Be}} $、$ ^{7}{\rm{Li}} $、$ ^{4}{\rm{He}} $和$ ^{3}{\rm{He}} $可以清晰的鉴别(为了方便区分$ ^{4}{\rm{He}} $和$ ^{3}{\rm{He}} $,用紫线示意两者理论分布)。

    图  3  (a) 0°~20°望远镜探测束流分布;(b) 将束流投影到x方向;(c) 将束流投影到y方向;红线为高斯拟合结果(在线彩图)

    图  4  蒙特卡罗模拟$ \theta $-$ \phi $二维图(在线彩图)

    (a) 各向同性粒子发射,模拟阵列立体角为36%4π; (b) 模拟卢瑟福散射,并考虑束流偏移。

    图  5  $ ^{7}{\rm{Be}} $+$ ^{120}{\rm{Sn}} $ 体系在48.05 MeV能量下的准弹性散射角分布(在线彩图)

    虚线为贝叶斯方法拟合的结果,蓝线是频率方法拟合的实验结果(使用的MINUIT程序 [28])。阴影区域为贝叶斯方法给出的1$ \sigma $置信区间。

    表  1  频率方法与贝叶斯方法拟合的光学势参数

    V0/MeV r0V/fm aV/fm W0/MeV r0W/MeV aW/MeV
    频率方法 39.57 1.14 0.77 49.00 1.06 0.83
    贝叶斯方法 27.72±9.41 0.26±0.05 0.67±0.38 28.95±9.04 1.23±0.09 0.69±0.10
    下载: 导出CSV
  • [1] KEELEY N, RAABE R, ALAMANOS N, et al. Prog Part Nucl Phys, 2007, 59(2): 579. doi:  10.1016/j.ppnp.2007.02.002
    [2] CANTO L, GOMES P, DONANGELO R, et al. Phys Rep, 2006, 424(1-2): 1. doi:  10.1016/j.physrep.2005.10.006
    [3] MAZZOCCO M. Int J Mod Phys E, 2010, 19(05n06): 977. doi:  10.1142/S0218301310015424
    [4] RAABE R, SIDA J, CHARVET J, et al. Nature, 2004, 431(7010): 823. doi:  10.1038/nature02984
    [5] KOLATA J J, AMRO H, BECCHETTI F D, et al. Phys Rev C, 2007, 75: 031302. doi:  10.1103/PhysRevC.75.031302
    [6] DI PIETRO A, RANDISI G, SCUDERI V, et al. Phys Rev Lett, 2010, 105: 022701. doi:  10.1103/PhysRevLett.105.022701
    [7] STANDYŁO L, ACOSTA L, ANGULO C, et al. Phys Rev C, 2013, 87: 064603. doi:  10.1103/PhysRevC.87.064603
    [8] JACKSON D. Rep Prog Phys, 1974, 37(1): 55. doi:  10.1088/0034-4885/37/1/002
    [9] AGUILERA E F, MARTINEZ-QUIROZ E, LIZCANO D, et al. Phys Rev C, 2009, 79: 021601. doi:  10.1103/PhysRevC.79.021601
    [10] JONSON B. Phys Rep, 2004, 389(1): 1. doi:  10.1016/j.physrep.2003.07.004
    [11] VERMA S, DAS J, JHINGAN A, et al. Eur Phys J A, 2010, 44(3): 385. doi:  10.1140/epja/i2010-10966-2
    [12] BARIONI A, ZAMORA J C, GUIMARÃES V, et al. Phys Rev C, 2011, 84: 014603. doi:  10.1103/PhysRevC.84.014603
    [13] ZAMORA J C, GUIMARÃES V, BARIONI A, et al. Phys Rev C, 2011, 84: 034611. doi:  10.1103/PhysRevC.84.034611
    [14] TABACARU G, AZHARI A, BRINKLEY J, et al. Phys Rev C, 2006, 73: 025808. doi:  10.1103/PhysRevC.73.025808
    [15] MORCELLE V, LICHTENTHÄLER R, LINARES R, et al. Phys Rev C, 2014, 89: 044611. doi:  10.1103/PhysRevC.89.044611
    [16] KALITA K, VERMA S, SINGH R, et al. Phys Rev C, 2006, 73: 024609. doi:  10.1103/PhysRevC.73.024609
    [17] LICHTENTHÄLER R, DE FARIA P, LÉPINE-SZILY A, et al. Eur Phys J Spec Top, 2007, 150(1): 27. doi:  10.1140/epjst/e2007-00257-9
    [18] MAZZOCCO M, KEELEY N, BOIANO A, et al. Phys Rev C, 2019, 100: 024602. doi:  10.1103/PhysRevC.100.024602
    [19] BASAK A, BILLAH M, KOBRA M, et al. Eur Phys Lett, 2011, 94(6): 62002. doi:  10.1209/0295-5075/94/62002
    [20] DE FARIA P N, LICHTENTHÄLER R, PIRES K C C, et al. Phys Rev C, 2010, 81: 044605. doi:  10.1103/PhysRevC.81.044605
    [21] YANG L, LIN C, YAMAGUCHI H, et al. Nat Commun, 2022, 13: 7193. doi:  10.1038/s41467-022-34767-8
    [22] ZHAN W, GUO Z, LIU G, et al. Sci China Math, 1999, 42(5): 528. doi:  10.1007/BF02882249
    [23] SUN Z, ZHAN W L, GUO Z Y, et al. Nucl Instr and Meth A, 2003, 503(3): 496. doi:  10.1016/S0168-9002(03)01005-2
    [24] HE J, XU S, MA P, et al. Nucl Instr and Meth A, 2012, 680: 43. doi:  10.1016/j.nima.2012.03.040
    [25] WANG D X, LIN C J, YANG L, et al. Compact 16-channel Integrated Charge-sensitive Preamplifier Module for Silicon Strip Detectors[Z].
    [26] WOODS R D, SAXON D S. Phys Rev, 1954, 95: 577. doi:  10.1103/PhysRev.95.577
    [27] THOMPSON I J. Comput Phys Commun, 1988, 7(4): 167.
    [28] JAMES F, ROOS M. Comput Phys Commun, 1975, 10(20): 343.
    [29] METROPOLIS N, ROSENBLUTH A W, ROSENBLUTH M N, et al. The Journal of Chemical Physics, 1953, 21(6): 1087. doi:  10.1063/1.1699114
    [30] YANG L, LIN C, ZHANG Y, et al. Phys Lett B, 2020, 807: 135540. doi:  10.1016/j.physletb.2020.135540
    [31] NIU Z, LIANG H. Phys Lett B, 2018, 778: 48. doi:  10.1016/j.physletb.2018.01.002
    [32] NIU Z, LIANG H, SUN B, et al. Phys Rev C, 2019, 99(6): 064307. doi:  10.1103/PhysRevC.99.064307
    [33] WANG Z A, PEI J, LIU Y, et al. Phys Rev Lett, 2019, 123(12): 122501. doi:  10.1103/PhysRevLett.123.122501
    [34] SIVIA D, CARLILE C, HOWELLS W, et al. Phys B, 1992, 182: 341. doi:  10.1016/0921-4526(92)90036-R
    [35] RONG C, RANGEL J, WU Y, et al. The European Physical Journal A, 2021, 57(4): 1.
    [36] KING G B, LOVELL A E, NEUFCOURT L, et al. Phys Rev Lett, 2019, 122: 232502. doi:  10.1103/PhysRevLett.122.232502
    [37] CANTO L, GOMES P, DONANGELO R, et al. Phys Rep, 2015, 596: 1. doi:  10.1016/j.physrep.2015.08.001
    [38] KOLATA J, GUIMARÃES V, AGUILERA E. Eur Phys J A, 2016, 52(5): 1.
  • 加载中
图(5) / 表 (1)
计量
  • 文章访问数:  176
  • HTML全文浏览量:  39
  • PDF下载量:  32
  • 被引次数: 0
出版历程
  • 收稿日期:  2023-01-12
  • 修回日期:  2023-01-29
  • 刊出日期:  2023-09-20

近垒能区7Be+120Sn的准弹性散射研究

doi: 10.11804/NuclPhysRev.40.2023005
    基金项目:  国家重点研发计划项目(2018YFA0404404);国家自然科学基金资助项目(12235020, U2167204, 12175313, 12175314, U1732145, U1867212和119611012);稳定支持基础科研计划;中核集团领创项目(LC192209000701, LC202309000201);青年人才发展基金项目(YC2122120001);国防基础科研计划项目(JCKY2020201C002)
    作者简介:

    常昶(1998−),男,山西大同人,硕士研究生,从事原子核反应研究;E-mail: changchangccc@outlook.com

    通讯作者: 杨磊,E-mail: yang_lei@ciae.ac.cn林承键,E-mail: cjlin@ciae.ac.cn
  • 中图分类号: O571.53

摘要: 近垒能区弱束缚核的反应动力学是目前核物理研究热点之一。本工作使用大立体角硅探测器阵列测量了$^{7}{\rm{Be}}$+$^{120}{\rm{Sn}}$体系在48.05 MeV的准弹性散射,结合蒙特卡罗模拟得到其微分截面。基于光学模型,分别用频率方法和贝叶斯方法对角分布进行拟合。在前角区,两种方法给出了一致的结果;在后角区,频率方法的结果呈现明显的振荡结构,而贝叶斯方法的结果较为平滑,在靠近180°有振荡上升的趋势。

English Abstract

常昶, 杨磊, 林承键, 杨彦云, 温培威, 骆天鹏, 马军兵, 许世伟, 王康, 段芳芳, 马南茹, 贾会明, 杨峰, 黄大湖, 张明昊, 杨过, 杨赟, 莫腾欢. 近垒能区7Be+120Sn的准弹性散射研究[J]. 原子核物理评论, 2023, 40(3): 356-361. doi: 10.11804/NuclPhysRev.40.2023005
引用本文: 常昶, 杨磊, 林承键, 杨彦云, 温培威, 骆天鹏, 马军兵, 许世伟, 王康, 段芳芳, 马南茹, 贾会明, 杨峰, 黄大湖, 张明昊, 杨过, 杨赟, 莫腾欢. 近垒能区7Be+120Sn的准弹性散射研究[J]. 原子核物理评论, 2023, 40(3): 356-361. doi: 10.11804/NuclPhysRev.40.2023005
Chang CHANG, Lei YANG, Chengjian LIN, Yanyun YANG, Peiwei WEN, Tianpeng LUO, Junbing MA, Shiwei XU, Kang WANG, Fangfang DUAN, Nanru MA, Huiming JIA, Feng YANG, Dahu HUANG, Minghao ZHANG, Guo YANG, Yun YANG, Tenghuan MO. Quasielastic Scattering Study for the 7Be+120Sn System at the Energy Near the Coulomb Barrier[J]. Nuclear Physics Review, 2023, 40(3): 356-361. doi: 10.11804/NuclPhysRev.40.2023005
Citation: Chang CHANG, Lei YANG, Chengjian LIN, Yanyun YANG, Peiwei WEN, Tianpeng LUO, Junbing MA, Shiwei XU, Kang WANG, Fangfang DUAN, Nanru MA, Huiming JIA, Feng YANG, Dahu HUANG, Minghao ZHANG, Guo YANG, Yun YANG, Tenghuan MO. Quasielastic Scattering Study for the 7Be+120Sn System at the Energy Near the Coulomb Barrier[J]. Nuclear Physics Review, 2023, 40(3): 356-361. doi: 10.11804/NuclPhysRev.40.2023005
    • 在过去的二十年里,随着加速器的发展,研究人员具备了产生远离稳定线核素的能力,开展了大量放射性束物理方面的研究工作,大大扩展了核物理研究的领域。库仑势垒附近能区弱束缚核的反应动力学研究,引起了低能核物理学界的极大兴趣[1-3]。弱束缚核结合能较低,通常具有奇特的结构,如团簇结构、中子晕和质子晕等[4-5]。弱束缚核在与靶核相互作用时较容易破裂,破裂之后可能会影响其他反应道,比如弹性散射过程发现有库仑虹部分消失或全部消失、四分之一角较小等[6-7]。弹性散射是核碰撞的一个基本的过程,主要发生在两核密度无交叠或者少量交叠的区域。虽然弹性散射过程比较简单,但它是研究核-核相互作用以及原子核的基本性质,如核物质的分布(尺寸、表面弥散等)的基本实验方法[8]。同时,通过测量弹性散射微分截面,可以提取总反应截面,这对于揭示弱束缚核反应机制至关重要。此外,通过测量弹性散射角分布,还可以进一步研究其它反应过程,如非弹、破裂、转移等,对弹性散射的耦合效应。在广义上,准弹性散射包括了弹性散射、非弹性散射、少数核子转移等接近弹性散射的周边反应过程。

      为了更清楚地了解近垒能区弱束缚核的反应动力学,我们选择$ ^{7}{\rm{Be}} $作为研究对象。$ ^{7}{\rm{Be}} $是位于质子滴线的放射性核,半衰期为53.2 d,容易破裂为$ ^{3}{\rm{He}} $$ ^{4}{\rm{He}} $集团,分离能为1.586 MeV。$ ^{7}{\rm{Be}} $因为是质子晕$ ^{8}{\rm{B}} $的核心,所以更加令人关注。对$ ^{7}{\rm{Be}} $的研究可以帮助我们更好地了解$ ^{8}{\rm{B}}$的反应机制[9-10]$ ^{8}{\rm{B}} $ 的质子的分离能只有0.138 MeV,具有质子晕结构,另外,$ ^{8}{\rm{B}} $在天体物理上也有着重要的意义。从实验角度来看,$ ^{7}{\rm{Be}} $直接破裂的产物$ ^{3}{\rm{He}} $$ ^{4}{\rm{He}} $都是带电粒子,很容易被探测到。因此,$ ^{7}{\rm{Be}} $也是研究近垒能区反应机制的理想对象。

      目前研究人员已经测量了$ ^{7}{\rm{Be}} $$ ^{7}{\rm{Li}} $, $ ^{12}{\rm{C}} $, $ ^{27}{\rm{Al}} $, $ ^{51}{\rm{V}} $, $ ^{58}{\rm{Ni}} $, 和 $ ^{208}{\rm{Pb}} $靶核的准弹性散射[10-18]。总体来说,除了$ ^{7}{\rm{Be}} $+$ ^{208}{\rm{Pb}} $的散射体系外,$ ^{7}{\rm{Be}} $的弹散数据均表明连续态对弹散的耦合效应并不显著。本工作中,我们测量了$ ^{7}{\rm{Be}} $+$ ^{120}{\rm{Sn}} $体系在48.05 MeV能量下的弹性散射角分布。选择$ ^{120}{\rm{Sn}} $作为靶核,是因为目前仍然缺少关于$ ^{7}{\rm{Be}} $与中等质量靶核的反应数据;同时该体系库仑相互作用和核相互作用都起着重要作用,使我们能够研究两者之间的相干效应。此外,人们发表了一系列弱束缚核在$ ^{120}{\rm{Sn}} $靶上的弹性散射数据,如$ ^{6, 7}{\rm{Li}} $[19], $ ^{6}{\rm{He}} $ [20], $ ^{8}{\rm{B}} $[21]。因此可以通过与同一靶核上的这些数据进行比较来进一步研究$ ^{7}{\rm{Be}} $的反应动力学。

    • 本实验是在中国科学院近代物理研究所的放射性离子束流线(RIBLL1)[22-23]上完成的。使用60.2 MeV 的$ ^{7}{\rm{Li}}^{3+} $初级束,流强为150~280 pnA,通过$ ^{1}{\rm{H}} $($ ^{7}{\rm{Li}} $, $ ^{7}{\rm{Be}} $)n反应($ Q $=−1.64 MeV)产生$ ^{7}{\rm{Be}} $次级束。氢气靶 [24]长80 mm,窗为2.5 μm的Havar膜,气压为800×102 Pa,液氦制冷,温度约为80 K。两个厚度为10 μm的塑料闪烁体(${\rm{C}} {_9}{\rm{H}}{_{10}}$)分别安装在RIBLL1的${\rm{T}} _1 $${\rm{T}} _2 $靶室,用来测量粒子的飞行时间(TOF)。经过RIBLL1的纯化和传输,$ ^{7}{\rm{Be}} $次级束进而轰击2.78 mg/cm2厚的自支撑同位素$ ^{120}{\rm{Sn}} $靶,丰度为98%,靶中心能量为48.05 MeV。次级束流强为$1\times{10^5}$ pps,纯度约为90%,主要污染物为初级束 $ ^{7}{\rm{Li}} $

      利用由9组硅望远镜组成的大立体角探测器阵列测量反应产物,如图1所示。阵列可分为前角和后角两部分。前角的每组望远镜包括三层硅探测器:第一层和第二层是单面硅条探测器(SSSD),厚度分别为20和300 μm,然后是厚度为1 500 μm的四分硅探测器(QSD)。后角的望远镜有着不同的结构:第一层是40 μm双面硅条探测器(DSSD),其次是两个厚的QSD,厚度分别为300和1 500 μm。SSSD和DSSD上每根硅条的宽度为3 mm,间隙是0.1 mm,灵敏探测面积是48 mm×48 mm,QSD的灵敏面积为50 mm×50 mm。前置放大器采用高性能电荷灵敏前置放大器 [25],与探测器密切配合,处于真空室中,以降低电子学噪声。本套探测系统共有502路信号,具有模块化、设计紧凑、几何效率高的优点。数据获取系统的触发条件是所有硅探测器的或逻辑(在0°处的硅探测器的触发做了100分除,避免计数率太高使得获取堵死)。TOF-$ E_{\text{r}} $(由0°望远镜测得的能量)如图2(a)所示,可以将次级束$ ^{7}{\rm{Be}} $与初级束$ { }^{7} {\rm{Li}}$$ ^{4}{\rm{He}} $分开。图2(b)显示了由一个前角望远镜测到的$\varDelta E{\text{-}}E_{\text{r}}$能谱。可以看出,$ ^{7}{\rm{Be}} $$ ^{7}{\rm{Li}} $以及反应产物$ ^{4}{\rm{He}} $$ ^{3}{\rm{He}} $可以清晰地鉴别。为了确定探测器的几何效率,我们利用蒙特卡罗方法模拟粒子各向同性发射。结果表明,阵列可以覆盖0°~21.2°, 29.8°~68.8° 和101.7°~164.2°的极角范围,立体角覆盖为36% 4π。望远镜总能量分辨约为3.1%,不足以分辨$ ^{7}{\rm{Be}} $的第一激发态($ E_{\rm ex} $=0.429 MeV, $ J^{\pi} $=1/2)。因此,得到的散射结果被认为是准弹性散射角分布。

      图  1  具有大立体角覆盖范围的硅探测器阵列示意图(红线表示束流入射方向)(在线彩图)

      图  2  粒子鉴别图(在线彩图)

      弹性散射微分截面的角分布由下式得到:

      $$ \frac{{\rm{d}} \sigma(\theta)}{{\rm{d}} \sigma_{\text {Ruth }}(\theta)}=C \frac{N(\theta)_{\text {expt }}}{N(\theta)_{\text {Ruth }}}, $$ (1)

      式中:$ {N(\theta)_{\text {expt }}} $是实验室角度$\theta $$ ^{7}{\rm{Be}} $弹性事件的产率;$ {N(\theta)_{\text {Ruth }}} $是蒙特卡罗模拟的卢瑟福散射产率。归一化常数C是一个全局归一化因子,它通过假设15°处的弹性散射截面为卢瑟福散射截面来确定。卢瑟福散射公式为

      $$ \frac{{\rm{d}} \sigma^{}_{\text {Ruth }}(\theta_{{\rm{c.m.}}})}{{\rm{d}} \varOmega}=\left(\frac{1}{4 {\pi} \varepsilon_0} \frac{Z_1 Z_2 e^2}{4 E}\right)^2 \frac{1}{\sin ^4 \frac{\theta_{{\rm{c.m.}}}}{2}}, $$ (2)

      其中:$ Z_1 $$ Z_2 $分别为弹核和靶核的质子数;$ e $为电子电荷,$ E $为弹核入射动能;$ \theta_{{\rm{c.m.}}} $为质心系角度,得到的截面为质心系截面。在数据处理过程中需要从实验室系转换到质心系,公式为

      $$ \frac{\sigma^{}_{{\rm{L}}}\left(\theta_{{\rm{L}}}\right)}{\sigma_{{\rm{c.m.}}}\left(\theta_{{\rm{c.m.}}}\right)} =\frac{\left(1+\gamma^2+2 \gamma \cos \theta_{{\rm{c.m.}}}\right)^{3 / 2}}{1+\gamma \cos \theta_{{\rm{c.m.}}}}, $$ (3)

      其中$ \gamma $为弹核和靶核质量数的比值。从卢瑟福散射公式中看到,如果$ \theta $在趋近于0的时候,截面会趋近于无限大,这显然和事实不符。因为卢瑟福散射公式在推导的过程中认为只发生单次散射并忽略核外电子的作用。而实际上在散射角很小时,核外电子与弹核的相互作用已经不能忽略,发生多次小角散射的概率也很大。因此,需要在合适的角度截断。在本工作的模拟过程中,截断角度设为8°。本阵列在0°放置有望远镜,可以方便地确定束流飘移位置[图3(a)],将其投影到xy方向见图3(b)、(c),可以确定束流中心x方向向右偏离3 mm,半高宽约为12 mm,y方向向上偏离3 mm,半高宽约为15 mm。根据束流偏移可以对散射角度做修正,图4(a)、(b)分别是模拟粒子各向同性发射和卢瑟福散射并考虑束流偏移的结果。

      图  3  (a) 0°~20°望远镜探测束流分布;(b) 将束流投影到x方向;(c) 将束流投影到y方向;红线为高斯拟合结果(在线彩图)

      图  4  蒙特卡罗模拟$ \theta $-$ \phi $二维图(在线彩图)

    • 光学模型(Optical Model)是核反应理论中最基础的理论模型,主要用来描述核反应的独立粒子阶段,这个阶段也是所有核反应都要经历的阶段。光学模型认为核碰撞过程可以描述为在弹靶系统组成的平均场中的运动,平均的核势场可以用一个复数的有效势来表示:

      $$ \begin{array}{l} U_{{\rm{N}}}(r)=V(r)+{\rm{i}} W(r), \end{array} $$ (4)

      $ V(r) $表示透射或者反射,虚部$ W(r) $表示吸收, 这个过程与光在介质中的传播类似,所以类似于光学中复折射率的方法。光学模型可以解决入射粒子在靶核平均场中的散射和吸收问题,利用它可以计算核反应的弹性散射截面和吸收截面(反应截面)。

      对于唯象光学势,实部和虚部可以写成:

      $$ \begin{array}{l} {\boldsymbol{U}}_{{\rm{N}}}(r)=V_0 f_V(r)+{\rm{i}} W_0 f_W(r), \end{array} $$ (5)

      其中:$ V_{0} $$ W_{0} $分别为实部、虚部深度;$ f(r) $是形状因子,通常取[26]

      $$ f(r)=\left[1+\exp \left(r-r_0\left(A_{{\rm{p}}}+A_{{\rm{t}}}\right)^{1 / 3} / a_0\right)\right]^{-1}, $$ (6)

      $ r_{0} $$ a_{0} $分别为作用势的半径和弥散参数,为可调参数;$ A_{\rm p} $$ A_{\rm t} $分别为弹核和靶核的质量数。

      基于OM,我们用FRESCO程序[27]对实验角分布做了拟合。FRESCO是一个耦合道计算程序,基于直接核反应理论框架,求解被认为在反应中起重要作用的原子核某种组分以及它们的相互作用势的薛定谔方程。拟合采用了2种方法:传统的频率方法和贝叶斯方法,见图5。频率方法是指在参数空间中逐步减小$ \chi^{2} $值来找到最优参数的方法,本文采用了MIGRAD算法[28]。贝叶斯方法是一种基于贝叶斯定理的统计推断方法。它的基本概念是利用先验概率和观测数据来计算后验概率,从而对未知参数或假设进行推断。在贝叶斯方法中,我们首先对未知参数或假设的先验分布进行设定,然后利用观测数据来更新这些分布,得到后验分布。这样,我们可以通过后验分布来对未知参数或假设进行推断和决策。根据贝叶斯定理,模型参数的后验概率分布函数可以写为

      图  5  $ ^{7}{\rm{Be}} $+$ ^{120}{\rm{Sn}} $ 体系在48.05 MeV能量下的准弹性散射角分布(在线彩图)

      $$ P ( x | D ) = \frac { P ( D | x ) P ( x ) } { P ( D ) } , $$ (7)

      P(x)为在没有给定数据集D的情况下,模型参数$ x $的先验分布,表示对参数集的先验知识。$ P(D|x) $是似然函数,通过比较实验数据与参数集为$ x $的模型值给出,采用了高斯分布的形式。$ P(D) $为归一化常数,保证后验分布是一个有效的积分为$ 1 $的概率密度。后验分布$ P(x|D) $代表了考虑数据集D后参数$ x $的最可能的分布。尽管贝叶斯定理的概念如此简单,但是在实际计算中,后验分布通常都是复杂高维的分布,很难进行抽样,这使得早期贝叶斯方法的使用受到了极大的限制。马尔可夫链蒙特卡罗方法(MCMC)的提出解决了这个问题[29]

      频率方法适合拟合少量参数,参数多时由于极小点众多,结果往往都在初值附近。这个方法的一个重要结论是,光学势的一些参数之间存在着很强的关联性。相比之下,贝叶斯方法通常依赖于马尔可夫链蒙特卡罗(MCMC)方法来对参数空间进行采样,并从似然函数和先验分布的乘积中获得后验预测,可以对参数的分布有一个整体展现,最后给出的结果是整体的期望值,但受先验分布的影响较大[30]。贝叶斯方法已广泛应用于核物理,并在弹性散射研究中与传统的频率方法进行了对比[30-36],该方法可以基于实验数据对光学势参数及其1$ \sigma $置信区间进行理论预测(图5中的阴影区域),这也是其优点之一,可以给出好的误差范围。两种方法拟合的光学势,见表1。对比发现,两种方法的光学势参数有一定的差别,尤其是贝叶斯方法给出的$r^{}_{0V}$太小、两者$ V_0 $$ W_0 $相差较大,贝叶斯的光学势参数普遍较小,频率方法的结果和文献基本一致。

      表 1  频率方法与贝叶斯方法拟合的光学势参数

      V0/MeV r0V/fm aV/fm W0/MeV r0W/MeV aW/MeV
      频率方法 39.57 1.14 0.77 49.00 1.06 0.83
      贝叶斯方法 27.72±9.41 0.26±0.05 0.67±0.38 28.95±9.04 1.23±0.09 0.69±0.10

      在前角区,两种方法给出了一致的结果,表面核力和库仑力相互干涉导致的库仑虹抑制不明显;在后角区,频率方法的结果呈现明显的振荡结构,可能是内部核力和库仑力相互干涉,也可能是近边散射和远边散射相干叠加,而贝叶斯方法的结果较为平滑,是对多个振荡的平均,但在近180°有振荡上升的趋势。OM拟合结果在大角度下与实验值有很大差异,在之前的研究中也出现了这样的问题[37-38],可能是因为非弹性过程和多步过程等因素的影响,需要进一步研究。

    • 本工作测量了弱束缚$ ^{7}{\rm{Be}} $+$ ^{120}{\rm{Sn}} $体系准弹性散射角分布。我们设计了立体角覆盖为36% 4π的硅探测器阵列。本阵列在0°放置了一组望远镜,在蒙特卡罗模拟卢瑟福散射时,可以利用这组望远镜进行束流漂移修正。对于实验结果,我们基于OM使用传统的频率方法和贝叶斯方法拟合,前角区两者结果相似,库仑虹抑制不明显,说明$ ^{7}{\rm{Be}} $的破裂道耦合效应较弱;后角区两者差异明显,且都与实验数据存在一定偏差。本实验补充了$ ^{7}{\rm{Be}} $与中重核的实验数据,将来可以比较$ ^{7}{\rm{Be}} $与不同质量靶核的系统学。在$ ^{7}{\rm{Be}} $的反应动力学方面,尤其是破裂碎片之间的关联还需要在未来进一步研究。

参考文献 (38)

目录

    /

    返回文章
    返回