近些年，研究弱束缚核^{6, 7}Li和⁹Be等作为弹核参与的熔合反应一直是实验和理论核物理研究的热点课题^[1]。经过大量研究发现，在这些弱束缚核所参与的熔合反应中，与理论计算以及紧束缚核相比，完全熔合在库仑势垒以上能区有着明显的压低，在库仑势垒以下区域则有着可观的增强^[2-12]。这些弱束缚核对熔合产生的影响吸引着人们更多地去探究背后的原因与物理机制。当弱束缚核作为弹核参与核反应时，由于较低的破裂阈值，很容易在与靶核发生熔合反应之前便破裂成多个碎片；此外，弹核与靶核之间还会发生转移反应，生成类弹核和类靶核，比如，²⁰⁸Pb(⁷Li, ⁶Li)²⁰⁹Pb、²⁰⁸Pb(⁶Li, ⁵Li)²⁰⁹Pb^[13]等，这些类弹核同样也会发生破裂反应。因此在弱束缚核熔合反应过程中存在直接破裂和转移破裂两种破裂反应道。由于破裂反应道的存在，导致弹核和类弹核被靶核或者类靶核完全俘获的概率下降，从而压低了完全熔合截面(Complete Fusion, CF)，增大了非完全熔合截面(Incomplete Fusion, ICF)。但是对于不同破裂反应道，例如中子转移破裂、质子转移破裂、直接破裂等，发生的概率以及对熔合截面的影响是不相同的，因此进一步研究弱束缚核熔合反应过程中不同破裂反应道机制是十分必要的^[14-17]。

目前国内外已开展了相关弱束缚核破裂反应机制研究的实验，主要测量方法是通过在线带电粒子测量，鉴别出不同破裂反应道，并通过理论模型^[18-19]将各破裂反应道的破裂概率从垒下能区外推到垒上能区。本文将围绕以下部分展开论述：第1节：主要的实验测量方法，第2节：近年来国内外研究进展，第3节：本课题组开展的相关工作，第4节：总结与展望。

目前涉及到的^{6, 7}Li和⁹Be等弱束缚核破裂反应机制研究，主要是通过在线带电粒子测量方法，利用位置灵敏探测器(一般为硅条探测器)测量反应过程中发射的带电粒子，并给出粒子准确的散射角度θ和能量信息。根据实验需求，可以进一步开展符合测量，利用两体动力学计算挑选出真实符合事件，并结合相关物理量开展后续分析。例如通过两个破裂碎片的能量，就可以计算出破裂反应道的Q值。具体计算公式为

其中：$ E_1 $、$ E_2 $为两个破裂碎片的动能；$ E_{\rm{recoil}} $为靶核(可能发生转移反应产生类靶核)的反冲能；$ E_{\rm{lab}} $为入射粒子动能；其中$ E_{\rm{recoil}} $可以由动量守恒定理可知。至此我们可以计算得出不同破裂反应道的Q值，并很好地再现破裂反应的Q值谱。当反应前后的粒子都处于基态时，这时反应的Q值一般被写作$ Q_{\rm{gg}} $，$ Q_{\rm{gg}} $是一个已知量。当反应后的粒子处于激发态时，拥有激发能$ E_{\rm{x}} $，此时的Q值为

在实验中，探测到的碎片都是处于基态，反冲的靶核(类靶核)可能处于激发态，因此可以通过重建不同破裂反应Q值谱，进一步确定反冲的靶核(类靶核)能级信息，从而开展下一步研究。

在破裂反应中，除了靶核(类靶核)可能处于激发态，破裂前的弹核(类弹核)也可能处于激发态，由于弹核(类弹核)激发态寿命很短，我们很难直接探测到，一般采用破裂碎片的相对能量来进行鉴定。在弹核(类弹核)破裂过程中如果没有外力作用，破裂碎片的相对能量可以通过下式计算：

式中$ Q_{\rm{BU}} $、$E_{\rm{P,\,x}}$分别为破裂反应的反应能和弹核(类弹核)激发能；Δ$E_{\rm{T,\,C}}$为相对能修正项，与弹核(类弹核)破裂位置有关。当弹核(类弹核)远离靶核(类靶核)破裂时，不会受到库仑场的作用，此时Δ$E_{\rm{T,\,C}}$为0，对应的$ E_{\rm{rel}}$可以呈现出一个峰结构。而当其靠近靶核(类靶核)破裂时，会受到库仑场的影响，改变碎片移动方向，此时Δ$E_{\rm{T,\,C}}$不为0，是一个连续分布的值，对应的$ E_{\rm{rel}}$会呈现一个很宽的分布。因此破裂碎片的相对能谱可以很好地描述弹核(类弹核)的末态，当然在这个分析过程中，我们需要考虑末态寿命的影响。

如图1所示，A为弹核，B为靶核，C为类弹核，D为类靶核，e、f分别为两个破裂碎片。破裂碎片的相对能量可以表达为$ E_{\rm{rel}}$ = $E_{\rm{1,\,c}}$+$E_{\rm{2,\,c}}$，其中$E_{\rm{1,\,c}}$、$E_{\rm{2,\,c}}$是两个碎片在以类弹核为质心系下的能量，可以通过已知量和测量的物理量进一步得到碎片相对能的计算公式：

图  1  类弹核破裂成碎片示意图

其中：$ m_{\rm{1}} $、$ m_{\rm{2}} $分别为两个破裂碎片质量；$ \theta_{12} $为两个破裂碎片在实验室系下的相对夹角，可以通过破裂碎片的散射角和方位角计算得到。通过两体动力学计算，我们还可以得到弹核(类弹核)在实验室下发射角度、两个破裂碎片在以类弹核为质心的相对动量方向夹角等物理信息。通过对在线带电粒子的动力学测量，我们可以得到很多有用的物理信息，并完成带电粒子来源分析和破裂反应道鉴别，从而进一步帮助我们更加清楚地重现破裂反应过程^[19-21]。

本节主要介绍近二十年来国内外关于弱束缚核作为炮弹的破裂反应研究。一方面利用符合测量方法，从实验上鉴别不同的破裂反应道，探究不同破裂反应道对于熔合压低的影响；一方面利用两体动力学计算得到的物理信息，确定弹核(类弹核)、靶核(类靶核)的末态信息，提取近靶破裂和共振破裂截面，进一步探究破裂反应机制；另外，破裂反应引起的非完全熔合过程与转移反应过程都可能会对完全熔合截面造成影响^[22]。目前围绕这个问题利用其他符合测量方法，开展了一些实验研究，并得出了相关结论。

2002年，Hinde等^[23]给出了⁹Be+^natPb实验中各破裂反应道的概率^[23]，并认为⁸Be的破裂可以分为两种：一种是瞬时破裂(Prompt Breakup)，一种是长寿命的基态⁸Be($ T_{1/2} $≈10⁻¹⁶ s)的延时破裂(Delayed Breakup)，并得出结论：弱束缚核熔合反应中垒上完全熔合压低现象是由反应中的瞬时破裂导致的，此外推测⁹Be垒上完全熔合压低因子与靶核的电荷息息相关。2010年，Rafiei等^[24]报道了同样的实验测量结果，束流为⁹Be，靶核分别为：²⁰⁹Bi、²⁰⁸PbS、¹⁹⁶Pt、¹⁸⁶W、¹⁶⁸Er、¹⁴⁴Sm，并通过两体动力学计算，对符合的α粒子进行Q值重构，并计算得到了破裂碎片相对能谱。通过将符合α粒子的Q值谱与相对能量$ E_{\rm{rel}}$谱进行关联，鉴别出了不同破裂反应道。图2展示了⁹Be+²⁰⁹Bi在$E_{\rm{cm}}=35.38\;\rm{MeV}$ 和26.73 MeV能量下Q-$ E_{\rm{rel}}$的关联谱，并认为图中的5个区域分别对应不同的破裂反应道。区域1事件来源于基态⁸Be的破裂；区域2事件来源于激发态⁸Be的破裂；区域3事件来源于激发态⁹Be的破裂，破裂产物除了探测到的α粒子以外还有一个中子，破裂过程可以认为是激发态⁹Be先破裂成基态⁸Be加一个中子，后续基态⁸Be再破裂成两个α粒子；区域4的事件来源于5/2⁻共振态⁹Be的破裂；区域5的事件来源于激发态⁹Be的破裂。并认为区域2和区域5的事件为瞬时破裂，可以对垒上完全熔合压低起作用，得出了与Hinde等^[23]相同的结论。弱束缚核熔合反应中的垒上完全熔合压低现象主要是由弱束缚核转移核子形成的类弹核发生瞬时破裂导致的。Luong等^[13]在2011年发表的文章中，利用带电粒子符合测量的方法，分别得到⁶Li+²⁰⁸Pb、⁷Li+²⁰⁸Pb系统下不同破裂反应道的Q值谱和破裂碎片相对能谱，进而提取了靶核(类靶核)、弹核(类弹核)的激发态信息，给出了各个破裂反应道中瞬时破裂的占比，并指出转移反应产生类弹核的破裂是占主要的。在分析过程中引入破裂时标的概念，如图3所示，进而鉴别了瞬时破裂和延时破裂过程，为进一步探究破裂反应机制提供了新思路。

图  2  ⁹Be+²⁰⁹Bi中符合α粒子Q-$ E_{\rm{rel}}$关联谱^[24]

图  3  碎片相对能$ E_{\rm{rel}}$与破裂距离、破裂反应时间关联图^[13]

Luong等^[19]基于2011年的工作，2013年系统整理了^{6, 7}Li在¹⁴⁴Sm,、^{207, 208}Pb、²⁰⁹Bi靶核上的实验结果，通过符合测量，重构了不同体系下的破裂反应Q值谱和相对能$ E_{\rm{rel}}$谱，最后得出了各个反应体系的瞬时破裂中每个破裂反应道占有的比例，从而得到每个体系中对熔合影响最大的破裂反应道。Kalkal等^[20]选用了轻质量的靶⁵⁸Ni和⁶⁴Zn，束流为^{6, 7}Li，通过对破裂碎片角度和能量进行符合测量，并将得到的破裂碎片相对能谱与其他重靶体系下的结果进行对比，发现⁷Li在撞击轻质量靶不会破裂成α和t，而⁶Li撞击不同质量的靶破裂成α和d的相对能量谱也有着很大差别。为了进一步区分不同破裂成分，Kalkal进行了破裂碎片的角度关联，明显观察到了^{6, 7}Li延时破裂带子结构，分别对应⁶Li的3⁺态和⁷Li的7/2⁻态，这种处于类弹核共振态下的延时破裂又称之为共振破裂，而靠近类靶核的瞬时破裂又称之为近靶破裂，两种成分在角度关联谱上可以很好地被鉴别开，如图4所示。因此从实验角度上得到了不同体系下共振破裂和近靶破裂的比例，并进行了模拟验证。最后得出相关结论：⁶Li在与不同质量的靶发生反应时，直接破裂的模式是极为不同的，在与轻质量靶发生反应时共振破裂起到了更大的作用；而且⁷Li和⁶Li与重质量靶发生反应时的直接破裂的模式也不一样。Kalkal分析认为这可能是和⁶Li的激发态相比，⁷Li的激发态的激发能远大于直接破裂的阈值。另外，在⁶Li与轻质量靶发生反应时，直接破裂中的近靶破裂的比例极小，这表明直接破裂反应道不会对熔合反应产生影响。

图  4  ^{6, 7}Li+²⁰⁹Bi在29 MeV下直接破裂角关联结果^[20]

Santra等^[25]在⁶Li+²⁰⁹Bi实验中对发射出来的α粒子和d粒子进行符合测量，从实验上给出了⁶Li-α+d共振破裂和近靶破裂的截面角分布(其文章中描述为Sequential、Direct破裂)，如图5所示，并用理论模型CDCC(Continuum Discretized Coupled Channels)进行计算。认为在⁶Li-α+d破裂反应道中，3⁺态⁶Li的共振破裂是占主要的。将共振破裂截面和近靶破裂截面结果同单举α截面结果比较，认为在这个过程中可能存在其他破裂反应道的贡献，因此进一步研究其破裂反应机制是十分有意义的。Chattopadhyay^[26]在近期的⁶Li+¹¹²Sn工作中，观察到了破裂碎片α和³He的符合事件，发现了新的破裂反应道，并在相对能谱中首次观察到了⁷Be的7/2⁻和5/2⁻激发态，分别给出了近靶破裂和7/2⁻和5/2⁻共振破裂截面角分布，与理论模型计算结果吻合较好，为单举α粒子的来源提供了一种新的解释。

图  5  单举α截面角分布和相关破裂反应道截面角分布^[25]

为了对弱束缚核破裂机制进一步研究，Shrivastava等^[27]在⁷Li+¹⁹⁸Pt实验中，利用带电粒子-γ射线符合测量方法，通过γ射线和带电粒子开窗，确定了类靶核信息，并提取了相关反应道截面。同时观察到^{5, 6}He大集团转移的现象，并利用相关动力学理论模型对实验数据进行了验证和计算，阐明了⁷Li不同集团结构在近垒能区的动力学计算中的重要性，为弱束缚核破裂机制的实验研究提供了新的测量方法。Cook等^[21]在2016年的工作中利用动力学模型模拟了⁹Be在不同靶上的破裂过程，并考虑了破裂时标的影响，得出结论，考虑类弹核激发态寿命对于重现破裂过程是十分重要的，并通过计算指出，破裂反应对于完全熔合截面压低因子的影响只有9%，没有达到理论上的30%，推测主要的压低是因为一步转移过程引起的。为了进一步验证该结论，Cook等^[22]开展了⁷Li+²⁰⁹Bi实验研究，通过带电粒子符合测量方法得到了两个破裂粒子均未被俘获的所有破裂反应道截面，约为(36±1) mb，其反应道贡献的α粒子称为No-Capture Breakup (NCBU) α。在该反应道截面中扣除不会对完全熔合截面产生压低的延迟破裂成分后，该截面约为20 mb，相比于已有的相同体系下的非完全熔合截面(302±21) mb^[15]，仅占一小部分，不能完全解释非完全熔合截面来源，指出应该有其他反应道对于非完全熔合截面的贡献。得到实验中没有粒子与之符合的α粒子能量和角度分布二维谱，称为Unaccompanied α，并与NCBU α分布对比，发现Unaccompanied α呈现能量更高、角度更前倾现象。通过模型计算给出了当其中一个破裂碎片被俘获后，产生的俘获破裂α粒子角分布，称之为Breakup Capture(BUC) α，其分布在后角成峰，与现有的Unaccompanied α角分布趋势明显不同，如图6所示，并指出这个矛盾是由于α粒子来源于不同的反应机制造成的。经过进一步研究发现，非完全熔合过程生成的Po同位素，仅有一小部分来源于弹核破裂再熔合过程，主要是来源于t的直接转移过程，由此得出结论：完全熔合截面的压低主要是弱束缚核中集团转移的影响，而不是弹核破裂机制影响的。

图  6  ⁷Li+²⁰⁹Bi实验中α粒子角分布^[22]

近期本课题组在中国原子能科学研究院的串列加速器上，使用30、40和47 MeV的^{6, 7}Li束流轰击210 μg/cm²的²⁰⁹Bi靶^[28]，通过覆盖大立体角的硅探测器阵列来实现带电粒子的符合测量。通过望远镜探测器系统清晰地实现了带电粒子的鉴别，并通过动力学符合测量计算，如图7所示，成功地在近垒和垒上能区鉴别出α+α，α+t，α+d，α+p四种破裂反应道，并给出了不同破裂反应道事件在总的破裂事件中的比例。值得注意的是，在30 MeV，⁶Li+²⁰⁹Bi实验数据中，我们新观察到了α+t破裂符合事件，说明在近垒能区，⁶Li的1n拾取过程是不可以被忽略的。目前课题组在⁶Li+²⁰⁹Bi实验数据中新发现了α+³He破裂反应道，进一步验证了反应过程中⁷Be的产生，为该体系下的1p拾取提供实验证据。我们利用重构的破裂反应Q值谱，可以进一步确定²¹⁰Bi、²⁰⁸Bi、²⁰⁸Pb等类靶核的末态信息，如图8所示，展示了三个能量点下⁷Li+²⁰⁹Bi中α+α符合事件重构的Q值谱，$Q_{\rm{gg}}=13.57\;\rm{MeV}$ ，我们可观察到²⁰⁸Pb的不同激发态结构，从而为我们进一步分析单核子转移机制提供了新的手段。

图  7  ^{6, 7}Li+²⁰⁹Bi不同入射能量下相对能谱对比^[28]

图  8  ⁷Li+²⁰⁹Bi不同入射能量下α+α符合事件Q值谱

利用两体动力学计算得到的角度信息，我们可以给出破裂碎片角度关联图像，如图9所示，图中$ \theta_{12} $为两个破裂碎片在实验室系下出射方向相对夹角，β为两个碎片在以类弹核为质心下相对动量方向夹角，具体可参见图1。我们可以明显观察到⁶Li+²⁰⁹Bi中⁸Be的基态和第一激发态结构，结合破裂相对能谱，我们可以清晰地分辨瞬时破裂和延时破裂结构，并可进一步提取截面信息。值得注意的是，与其他开展的工作相比，我们在30 MeV ⁶Li+²⁰⁹Bi数据中，首次观测到了⁸Be的第一激发态结构，发现随着入射能量的提升，⁸Be共振破裂成分逐渐增加。对比^{6, 7}Li+²⁰⁹Bi中α+t破裂道相对能谱，如图10所示，发现相比于库仑激发的⁷Li，来自⁶Li转移后的⁷Li可以布局到更高的激发态，目前截面和理论计算工作仍在进行中。

图  9  30，40 MeV下⁶Li+²⁰⁹Bi中α+α符合角关联谱及相对能谱

图  10  40 MeV下 ⁶Li+²⁰⁹Bi(a)与⁷Li+²⁰⁹Bi(b)数据α+t破裂道相对能谱

本文主要介绍了目前弱束缚核破裂反应机制的测量方法，并整理了目前国内外相关实验研究。由于弱束缚核的结合能较低，其作为弹核参与的重离子熔合反应过程存在众多机制。目前开展的实验研究中，主要是通过符合测量等多种手段，对熔合反应过程中存在的各种破裂过程进行了鉴别，并解释了单举α粒子的来源，一定程度上解释了各种破裂过程对于熔合压低的影响。但是，目前的实验数据主要集中在垒下或者近垒能区，对于垒上的破裂机制研究数据较少；从破裂反应Q值谱提供的物理图像，结合带电粒子-γ符合测量方法，很好地描述破裂反应后类靶核末态信息，从而为单核子转移机制研究等提供新的视角；我们可以从相对能谱中提取弹核(类弹核)的末态信息，结合末态寿命，鉴别瞬时破裂和延时破裂成分，进一步对熔合截面压低做出解释，同时可以进一步探究背后的破裂反应机制。

由于破裂反应引起的非完全熔合过程与转移反应过程都可能会对完全熔合截面造成影响^[22]，如果可以进一步鉴别和提取两个反应过程的截面，会极大地帮助我们理解这两个过程对于完全熔合截面的影响。但是，目前从实验的角度上实现非完全熔合与转移反应产物仍存在困难，我们需要在实验装置和测量方法等方面进行提升。在以后的研究工作中，本课题组也会进一步利用带电粒子-γ符合测量方法开展原子核反应机制方面的研究。