0 引 言
相对论能区和中高能区的核电荷拾取反应在很多基础研究和应用领域中都具有重大的意义, 例如天体物理学、 宇宙辐射和太空飞行, 以及放射性束流的产生。在BNL、 SPS、 GSI和LBL的不同能区的重离子作用于各种靶的核电荷拾取反应都有相关研究
[1‐7]。理论上讲, 电荷拾取反应应该是边缘以及表面两种作用过程的贡献, 对于边缘作用, 电荷拾取反应截面
应正比于
; 对于表面作用电荷拾取反应截面
应正比于
。综合两种效应的贡献,
应正比于
。Ren
[1]和Nilsen
[2]等用大量的实验数据得到, 1 A GeV能区附近的电荷拾取反应截面对射弹核质量
的幂函数关系的指数大体上在2倍因子范围内, 可以用表达式
(以mb为单位)给出, 其中
表示外围碰撞。Sümmerer等
[3]通过测量电荷拾取反应产物中同位素的分布, 说明了电荷拾取反应截面对射弹核质量的二次依赖关系很大程度上是由于初级碎片带电粒子的蒸发, 即射弹核的一个中子转换成了一个质子, 涉及到Δ共振。同样地, Bachelier
[4], Roy-Stephan
[5]等利用磁谱仪测量了0.95 A GeV的
20Ne作用于Al靶得到
20Na产物的电荷转移反应中的能量转移, 检测到两个峰, 一个是由np→pn引起的低激发能峰, 一个是由靶核激发到Δ共振引起的300 MeV的峰值。
在中能和低能区且低于费米能的情况下, 电荷拾取反应的机制是转移反应
[6], 即质子拾取和中子剥离过程。在此能量下, 射弹核速度小于原子核内核子的费米速度, 接触时费米球重叠, 使得核子可以从一个费米海跳到另一个费米海。但在相对论能量下, 射弹核和靶核的费米范围完全不重叠, 这就阻止了靶核质子向射弹核的转移, 则Δ共振的形成和衰退在高能区的电荷拾取反应过程中是最可能的基本过程。在这个过程中, 射弹核的一个中子转化成一个质子, 通过
反应, 质子被射弹核吸收, 释放出
[7]。综上所述, 在中能和高能区的核-核碰撞过程中, 质子转移反应和Δ共振的形成与衰变这两种机制应该是同时发挥作用。
目前的数据集中在高能区的电荷拾取反应截面的研究, 对于中低能区的数据研究较少。为了研究电荷拾取反应与射弹核质量之间的关系, 需要排除能量和靶核质量的影响。其中一种方法是在有限的能量区域内作用单一的靶。本文研究的是一定能量范围(400 A MeV能区附近) 的12C、 40Ar、 56Fe、 84Kr射弹分别作用于Al、 C、 CH2、 H靶的核电荷拾取反应截面与射弹核质量的关系。
1 实验部分
1) 利用日本放射线医学综合研究所重离子医用加速器(NIRS HIMAC)产生的射弹核束流分批照射Al、 C、 CH
2与CR-39固体核径迹探测器组成的复合靶(垂直照射)。每种复合靶由4个靶和10片CR-39放置排列而成, 相邻两个靶中间放置两片CR-39固体核径迹探测器。
图 1 为靶与CR-39固体核径迹探测器
[8]放置的示意图。
实验所用Al、 C、 CH2与CR-39组成的复合靶的面积为5.0×5.0 cm2, 束流通量约为3 000离子/cm2(具体根据射弹核束流的电荷数而定, 以减少束流径迹的相互重叠)。这里的束流能量为出加速器时的能量, 束流到达各个靶上表面时的能量可以利用各个靶的厚度及CR-39探测器的厚度通过SRIM计算软件分别计算得到。每个靶上表面的束流能量的不确定度根据束流线长度的不确定度和CR-39探测器的不确定度计算得来, 结果小于1 A MeV。
2) 带电粒子进入CR-39探测器后, 会在它的径迹区域形成一定的辐射损伤, 形成潜径迹, 将CR-39探测器放在温度为70 ℃, 浓度为7 mol/L的NaOH溶液中蚀刻30 h。蚀刻完成后, 把CR-39放置在去离子的纯净水中清洗10 h左右。重离子及其碎片在CR-39探测器中形成的径迹经过化学蚀刻后在CR-39上下表面的截面为圆形或接近圆形(垂直照射)。
3) 利用HSP-1000高影像图像分析显微镜扫描CR-39蚀刻径迹探测器中记录的核径迹。
4) 利用图像分析软件PitFit测量CR-39蚀刻径迹探测器记录的各种核碎片蚀刻径迹的坐标、 面积等信息。由于PitFit分析软件对径迹进行自动拟合的过程中会出现径迹识别错误的现象, 所以在自动拟合后本实验采用手动再进行一遍修正。手动拟合过程是一个非常耗时的过程, 但是它可以使检测率达到100%。
5) 先对分析软件测量出来的CR-39固体蚀刻核径迹探测器中的蚀刻径迹截面面积进行修正
[9],然后用截面面积大小与弹核碎片所带电荷数成正比的方法来判定弹核碎片的电荷 (
)。对于薄靶中弹核的电荷拾取反应, 弹核碎片的速度本质上与弹核的速度是一样的, 则有如下推论: 不论是弹核还是碎片核的蚀刻坑的大小或者短轴与它的电荷之间是单一线性函数关系。由于粒子在探测器上蚀刻径迹截面的面积与粒子的电荷数
呈线性关系, 所以本实验用蚀刻径迹在CR-39表面的截面面积分布来确定带电粒子的电荷, 从而对带电粒子进行鉴别。
以下以462 A MeV的
56Fe作用于CH
2靶为例对具体实验过程加以说明。
图 2 是462 A MeV的
56Fe作用于CH
2靶后的
56Fe束流和碎片的蚀刻径迹面积分布图。
56Fe束流和一直到电荷数
的碎片可以清楚地看到峰值, 而对于电荷数
的碎片则无法从图中识别。本实验采用了7个高斯分布叠加拟合, 得到了
56Fe束流和电荷数
的碎片的平均蚀刻径迹面积和误差。
图 3 是
图 2 中的
56Fe束流和电荷数
的碎片的电荷与高斯拟合蚀刻径迹面积的关系图。可以看到, 蚀刻径迹面积随电荷数增加而线性增加。利用这个线性拟合, 得到
以及
。根据这个关系可以确定其他碎片的电荷。
按照同样的方法, 确定了每一种束流在不同靶上的束流和碎片的电荷与蚀刻径迹面积的依赖关系, 并识别出了相应的弹核碎片的电荷。
6) 根据图像分析软件PitFit测到的CR-39蚀刻径迹探测器记录的各种弹核碎片蚀刻径迹的坐标, 利用径迹重建技术
[10](包括蚀刻径迹截面面积修正、 标准点选取、 利用最小二乘法进行径迹重建以及径迹匹配四个主要步骤)确定发生核碎裂以及电荷拾取的弹核和碎片的数量, 并计算电荷拾取反应截面。在径迹匹配这个关键步骤中, 由于径迹匹配数据量巨大, 本实验采用Fortran语言程序来完成径迹匹配。由于测量值与理论值之间存在误差, 故要求在合理的误差范围内完成径迹匹配,
图 4 为462 A MeV的
56Fe作用于CH
2靶的坐标差值
和
的分布图, 在4倍标准偏差范围内
搜索匹配径迹, 这样, 库仑散射效应的影响便可以忽略不计, 使得匹配精度可以达到99.99%以上。
图 4 中, 每个图像中的光滑曲线是高斯分布拟合。当在4倍标准偏差范围内匹配径迹时, 会出现以下5种情况:
1) 在4倍标准偏差范围内有1个匹配径迹, 而且匹配径迹的面积在束流的面积范围内, 则考虑没有发生碎裂。
2) 在4倍标准偏差范围内有1个匹配径迹, 而且匹配径迹的面积小于束流的面积范围, 则考虑发生了碎裂。
3) 在4倍标准偏差范围内有2~3个匹配径迹, 而且匹配径迹的电荷之和与束流电荷相同或略小, 也考虑发生了碎裂。
4) 在4倍标准偏差范围内有1个匹配径迹, 而且匹配径迹的面积大于束流的面积范围, 则考虑发生了电荷拾取反应, 如
图 5 中的21212号编码的核径迹, 也就是本文所研究的对象。
5) 在4倍标准偏差范围内没有匹配径迹, 则认为是一个完全碎裂事件, 只是产生的碎片电荷数低于CR-39探测器的阈值。
2 计算与分析
为了测量电荷拾取反应截面, 最主要的实验要求就是获得一个足够的电荷分辨率, 以便于从大量的穿过靶后未发生电荷变化的射弹核中提取出相对稀少的电荷拾取的碎片。为此, 本文选择了上节实验部分中径迹匹配中的情况(1)和情况(4)来分析。
图 6(a) 是
56Fe束流作用在CH
2靶, 过靶后未发生反应的
56Fe束流粒子和电荷增大的碎片的蚀刻径迹面积分布图。
图 6(b) 是
图 6(a) 从中间断开的放大图, 可以更清楚地显示蚀刻径迹面积分布。对过靶后未发生反应的束流粒子的蚀刻径迹面积分布图进行高斯拟合, 然后确定电荷拾取事件的数量。
式中: 为靶的原子序数; 为靶的厚度; 为靶的密度; 为阿伏伽德罗常数; 为过靶后没有发生电荷变化的束流粒子的数量; 为电荷拾取事件的数量。
表 1 是本实验组以前实验测量和分析计算得到的400 A MeV能区附近不同能量
12C、
40Ar、
56Fe、
84Kr射弹作用于Al、 C、 CH
2靶的核电荷拾取反应截面的相关数据, 这里的误差为统计误差。
表 1 中列出了过靶后未发生电荷变化的束流粒子数量和发生电荷拾取反应的粒子数量, 电荷拾取反应截面由
式(1)计算得出。有关实验过程、 轨迹重建以及数据处理的详细信息可参见之前的研究工作
[12‐14]。
对于同一种靶, 考虑到在一定束流能区内, 核电荷拾取反应截面在实验误差内一致, 所以本文对表中所包含的同种靶的电荷拾取反应截面取平均值。
表 1~
表 4 中的mean数值为同种靶的电荷拾取反应截面平均值。同时, 利用C和CH
2靶上的实验结果, 用
.5
关系式分别计算出不同射弹作用于H靶的电荷拾取反应截面, 如4个表中所示。H靶上的束流能量为C靶和CH
2靶前表面束流能量的平均值。
整理
表 1 ~
表 4 的实验数据中的电荷拾取反应平均值mean, 可以得到一定能量范围(400 A MeV能区附近)内, 同一种靶不同射弹的电荷拾取反应截面数据如
表 5。
对
表 5 中的数据分别进行幂函数和线性函数拟合,
图 7 和
图 8 分别为本文所研究能量区域的不同靶的电荷拾取反应截面对射弹核质量的幂函数和线性函数关系拟合图。
幂函数拟合关系式为
式中:
和
为拟合参数, 具体见
表 6。相应的最小值
也见
表 6。
线性函数拟合关系式为
式中:
和
为拟合参数, 具体见
表 7。相应的最小值
也见
表 7。
由
图 7 及
表 6 可以看出, 本实验结果利用幂函数
式(2)来拟合的幂拟合参数
在-0.01~0.71之间。这种对于射弹核质量
P 的依赖关系, 与在更高束流能量更重射弹核观察到的结果
[1‐2,5]不同。由
图 8 及
表 7 可以看出, 本实验结果利用线性函数
式(3)拟合的拟合参数
在-0.004~0.24之间, 此结果在一定程度上证实电荷拾取反应过程应为边缘及表面反应过程的集体贡献。
3 结 论
对于一种特定类型的相互作用, 核截面对质量
的依赖强度长期以来都被用作相关机制的线索。在核物质中, 反应截面对
的线性依赖关系通常表示射弹核中长的平均自由程; 对
的依赖关系表示非常短的平均自由程, 因此是表面相互作用; 对
的依赖关系表示周长处的相互作用。本实验结果在一定程度上证实电荷拾取反应过程应为边缘及表面反应过程的集体贡献, 但这两种过程又不能完全解释结果。在中高能重离子电荷拾取反应中, 因为是高激发能且是边缘碰撞, 所以会出现一种反应机制就是射弹核中的一个中子被激发到Δ共振态, 随后马上衰变为一个质子, 然后释放出一个π介子, 它的反应截面应与射弹核的中子密度和分布有很大的关系。Δ共振态(
=3/2,
=3/2)是核子的一种激发态, 这个激发态是通过旋转一个夸克的自旋和同位旋得到的。在中高能重离子碰撞的电荷拾取反应中, 主要是Δ共振引起了自旋-同位旋反应
[19]。质子转移反应和Δ共振的形成与衰变, 在中能区的核-核碰撞过程中应是同时发挥作用的。
山西省青年科技研究基金资助项目(202203021212396)