-
本实验是在中国科学院近代物理研究所的放射性离子束流线(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−)。因此,得到的散射结果被认为是准弹性散射角分布。弹性散射微分截面的角分布由下式得到:
$$ \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)],将其投影到x和y方向见图3(b)、(c),可以确定束流中心x方向向右偏离3 mm,半高宽约为12 mm,y方向向上偏离3 mm,半高宽约为15 mm。根据束流偏移可以对散射角度做修正,图4(a)、(b)分别是模拟粒子各向同性发射和卢瑟福散射并考虑束流偏移的结果。 -
光学模型(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]。贝叶斯方法是一种基于贝叶斯定理的统计推断方法。它的基本概念是利用先验概率和观测数据来计算后验概率,从而对未知参数或假设进行推断。在贝叶斯方法中,我们首先对未知参数或假设的先验分布进行设定,然后利用观测数据来更新这些分布,得到后验分布。这样,我们可以通过后验分布来对未知参数或假设进行推断和决策。根据贝叶斯定理,模型参数的后验概率分布函数可以写为$$ 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],可能是因为非弹性过程和多步过程等因素的影响,需要进一步研究。
Quasielastic Scattering Study for the 7Be+120Sn System at the Energy Near the Coulomb Barrier
-
摘要: 近垒能区弱束缚核的反应动力学是目前核物理研究热点之一。本工作使用大立体角硅探测器阵列测量了
$^{7}{\rm{Be}}$ +$^{120}{\rm{Sn}}$ 体系在48.05 MeV的准弹性散射,结合蒙特卡罗模拟得到其微分截面。基于光学模型,分别用频率方法和贝叶斯方法对角分布进行拟合。在前角区,两种方法给出了一致的结果;在后角区,频率方法的结果呈现明显的振荡结构,而贝叶斯方法的结果较为平滑,在靠近180°有振荡上升的趋势。Abstract: The reaction kinetics of weakly bound nuclei in the nearbarrier energy region is currently one of the hotspots in nuclear physics research.The quasielastic scattering of the$^{7}{\rm{Be}}$ +$^{120}{\rm{Sn}}$ system at 48.05 MeV was measured using a large solid angle covered silicon detector array, and its differential cross section was obtained in combination with Monte Carlo simulations. Based on the optical model, the angular distributions were fitted by the frequentist method and the Bayesian method, respectively. At the forward angles, the two methods give consistent results; at the backward angles, the results of the frequentist method show an obvious oscillatory structure, while the results of the Bayesian method are smooth, with an oscillatory upward trend near 180°. -
图 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 -
[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.