Energy-differential measurement of the $^{\mathrm{nat}}$C(n,p) and $^{\mathrm{nat}}$C(n,d) reactions at the n_TOF facility at CERN

本文在 n_TOF 设施上测量了天然碳的(n,p) 和(n,d) 反应在高达 25 MeV 能量范围内的微分截面，结果显示与主要评价库存在显著差异，但与 TALYS-2.0 计算结果表现出意外的一致性，尤其是对于(n,p) 反应。

要点

想象一下，你正在试图理解一种特定类型的“子弹”（中子）如何与一个非常常见的靶材发生相互作用：一块碳块（就像铅笔中的石墨）。当这些子弹击中碳时，它们有时会敲出更小的碎片，比如微小的弹珠（质子）或稍重的弹珠（氘核）。

欧洲核子研究中心（CERN）的 n_TOF 设施（一台向靶材发射中子的巨型机器）的科学家们决定精确测量这种情况发生的频率以及涉及的能量。他们专注于两种特定的反应：

1. (n,p) 反应：中子击中碳，一个质子飞出。
2. (n,d) 反应：中子击中碳，一个氘核（一个质子和一个中子粘在一起）飞出。

以下是他们做了什么、如何做的以及发现了什么的简单故事。

设置：高速摄像机与碳铅笔

科学家们没有使用普通的相机，而是使用了“飞行时间”技术。想象一条长 182.5 米的跑道。

• 他们向一个铅靶发射了一束质子，产生了一股中子流。
• 这些中子在长长的跑道上奔跑。
• 因为它们速度很快，到达终点所需的时间告诉科学家们它们确切拥有多少能量。中子越快 = 能量越高。

在这条跑道的中间，他们放置了一片非常薄的天然碳片（大约和人类头发一样厚）。这片碳片周围环绕着两组硅望远镜。把这些望远镜想象成高科技的三明治探测器。

• 第一层（薄切片）：一层非常薄的硅，测量粒子穿过时损失了多少能量（就像减速带）。
• 第二层（厚切片）：一层较厚的硅，捕获粒子并测量其剩余的总能量。

通过比较“减速带”能量与“总能量”，科学家们可以区分质子和氘核，尽管它们看起来非常相似。这就像通过它们撞击墙壁后的反弹方式来区分乒乓球和高尔夫球。

挑战：整理混乱的数据

他们收集的数据是一种混乱的混合体。当中子击中碳时，它并不会只产生一个干净的结果。它可能会让剩余的碳核处于一种“兴奋”状态（激发态），就像钟被敲击后会以特定的音调鸣响一样。

• 原子核可能处于其“平静”状态（基态）或各种“激发”状态。
• 每种状态产生的粒子具有略微不同的能量和方向。

为了理清这些，科学家们必须使用计算机模型（TALYS-2.0）。把这个模型想象成一本预测碳核行为的高级食谱书。他们不只使用一种食谱；他们尝试了480 种不同的食谱变体，以观察结果变化了多少。这至关重要，因为如果食谱错了，他们的测量结果也会错。

他们还使用了人工智能（神经网络）。由于数据中的粒子彼此非常接近，人眼无法轻易将质子与氘核分开。他们训练计算机识别每种粒子类型的独特“指纹”，就像俱乐部里一位非常聪明的保镖，确切知道谁应该排在哪个队伍里。

重大发现：“缺失”的能量

当科学家们最终计算出结果时，他们发现了一些令人惊讶的事情。

“图书馆”与“现实世界”
科学家们通常依赖“数据图书馆”（就像物理书籍图书馆）来告诉他们当中子击中碳时会发生什么。这些图书馆被用于设计核反应堆、医疗设备和太空防护罩。

• 预期：图书馆表示反应应该以某种频率发生（特定的“截面”）。
• 现实：n_TOF 团队发现，反应发生的频率显著高于图书馆的预测，尤其是对于质子反应。

这就像天气预报说下雨的概率是 10%，但当你走到外面时，却下起了倾盆大雨。现有的“预报”（数据图书馆）低估了这场风暴。

一线希望
有趣的是，他们新的、更详细的测量结果与 TALYS-2.0 计算机模型 的预测非常吻合。这表明计算机模型一直以来都是正确的，但“图书馆”（科学家们使用的书籍）包含的是过时或不正确的信息。

为什么这很重要？

这篇论文解释说，这不仅仅是一场理论游戏。碳无处不在：

• 在我们的身体里：它是组织的主要组成部分。
• 在医学中：它被用于癌症治疗（强子疗法）。
• 在太空中：它被用于卫星的防护罩。

当高能中子在这些环境中击中碳时，它们会产生次级粒子。如果我们不知道这种情况发生的确切频率，我们就无法准确计算患者接受的辐射剂量，或者太空飞船防护罩的效果如何。

结论

该团队以高精度成功测量了这些反应，从反应开始的时刻（约 14-15 MeV）一直到 25 MeV。

• 他们证明了反应发生的频率高于当前标准数据所表明的。
• 他们证实，他们的结果与特定的计算机模型（TALYS-2.0）一致，但与当今工程师和医生使用的主要数据图书馆不一致。

简而言之，他们取了一片非常薄的碳片，用高速中子轰击它，利用人工智能和超级计算机对碎片进行分类，并发现碳对中子反应的“规则手册”需要进行重大更新。

技术摘要

以下是论文《在 CERN n_TOF 设施上对 natC(n,p) 和 natC(n,d) 反应进行能量微分测量》的详细技术总结。

1. 问题陈述与动机

轻核（特别是碳）的中子诱发反应对于基础核物理、医疗应用（强子治疗）、辐射屏蔽以及中子敏感金刚石探测器的校准至关重要。

• 差距： 尽管 $^{nat}$C(n,p) 和 $^{nat}$C(n,d) 反应此前已得到研究，但在反应阈值上方（10–25 MeV）的能量范围内的截面数据仍缺乏有效约束。现有的实验数据存在显著的不一致性，主要核数据评估库（如 ENDF/B-VIII.1、TENDL-2023）彼此之间以及与积分测量结果之间往往存在分歧。
• 触发因素： 此前在 n_TOF 进行的一次积分截面测量得出的数值显著高于任何现有评估值，这表明现有库中可能存在系统性误差，或者需要更精确的微分数据。
• 目标： 对 $^{nat}$C(n,p) 和 $^{nat}$C(n,d) 反应进行高达 25 MeV 的高精度能量微分测量，以解决差异并提供基准数据。

2. 方法论

实验装置

• 设施： 实验在 CERN 的 n_TOF（中子飞行时间） 设施进行，具体位于 EAR1（实验区 1），飞行路径为 182.5 m。这条长路径确保了在高能中子下的高能量分辨率。
• 靶材： 使用天然碳（$^{nat}$C）石墨样品（5 cm $\times$ 5 cm，厚度 0.25 mm），倾斜 45° 放置，有效厚度为 0.35 mm。天然碳由约 98.9% 的 $^{12}$C 和约 1.1% 的 $^{13}$C 组成。
• 探测器： 两个位置灵敏的 $\Delta E$-$E$ 硅望远镜 被放置在束流之外（覆盖 20° 至 140° 的角度）。
  • 每个望远镜由一层 20 $\mu$m 厚的 $\Delta E$ 层和一层 300 $\mu$m 厚的 $E$ 层组成。
  • 探测器被分割为总共 64 条带（每个望远镜 32 条），实现了颗粒化探测。
• 束流： 来自 CERN 质子同步加速器的脉冲 20 GeV 质子束轰击铅靶，通过散裂产生宽谱中子。

数据分析与重建

分析涉及四个主要步骤：

1. 中子能量校准： 中子能量（$E_n$）通过相对于质子脉冲（通过壁电流监测器同步）的飞行时间（TOF）确定。对于 >4 MeV 的能量，束流的分辨率函数采用高斯分布建模。
2. 粒子鉴别： 采用 机器学习（神经网络） 方法，根据 $\Delta E$-$E$ 特征区分质子、氘核、氚核和背景粒子。这是必要的，因为质子和氘核在参数空间中的模式非常接近。
3. 截面重建：
   • 该问题被表述为一个逆问题（第一类 Fredholm 积分方程），将观测计数与底层截面联系起来。
   • 模型输入： 由于子核（$^{11}$B、$^{12}$B、$^{13}$B）激发态的角分布和分支比未直接测量，分析依赖于 TALYS-2.0 核反应代码的计算结果。
   • 不确定度量化： 为了解决模型依赖性，作者进行了 480 组不同的 TALYS 模型集 计算（改变能级密度模型、光学模型和伽马强度函数等参数）。他们还测试了“人工”分布（各向同性角分布和简化的分支比），以建立保守的系统不确定度。
   • 同位素分离： 由于 $^{13}$C 丰度较低且无法按同位素区分反应，最终结果针对 天然碳（$^{nat}$C）进行重建，并按天然丰度加权。

3. 主要贡献

• 高精度微分数据： 提供了从阈值到 25 MeV 的 $^{nat}$C(n,p) 和 $^{nat}$C(n,d) 的首个高分辨率、能量微分截面数据。
• 先进分析技术： 成功将神经网络应用于具有挑战性的 $\Delta E$-$E$ 区域中的粒子识别，并利用严格的统计框架（带拉格朗日乘子的最小二乘最小化）在考虑中子束分辨率函数的情况下重建截面。
• 系统不确定度评估： 使用 384 组 TALYS 模型集和人工分布进行了详尽的敏感性分析，表明模型相关的不确定度对于 (n,p) 反应低于 10%，对于 (n,d) 反应低于 20%。
• TALYS 验证： 证明虽然标准评估库失效，但只要选择正确的模型参数，特定的 TALYS-2.0 配置可以准确重现实验数据。

4. 主要结果

• 截面值：
  • $^{nat}$C(n,p)： 测量的截面 显著高于（幅度很大）主要评估库（TENDL-2023、ENDF/B-VIII.1、JEFF）的预测值。数值在 22 MeV 处达到峰值，约为 75 mb。
  • $^{nat}$C(n,d)： 数据也显示出与库的差异，尽管与某些数据集（如 TENDL-2015）的吻合度优于 (n,p) 反应。
• 与库的比较：
  • n_TOF 的结果与 TENDL-2023、ENDF/B-VIII.1 和 JEFF-3.3 ** largely 不一致**。
  • 结果显示与特定的 TALYS-2.0 计算（特别是使用后移费米气体模型计算能级密度的那些）存在 意想不到的吻合。
• 与过去数据的比较：
  • 新数据与 Kreger & Kern (1959) 和 Bobyr 等人 (1972) 的旧测量值吻合良好。
  • 它与当前库的基础（如 Rimmer & Fisher 数据）相矛盾。
  • 它证实了 n_TOF 早期的 积分测量，验证了高截面是一个真实的物理现象，而非积分方法的假象。

5. 意义与影响

• 核数据评估： 研究结果需要对主要核数据库中的 $^{nat}$C(n,p) 和 (n,d) 评估进行修订。当前对这些截面的低估可能导致辐射输运模拟中的误差。
• 医疗应用（强子治疗）： 准确的截面对于计算质子和离子治疗期间健康组织中的次级中子剂量至关重要。较高的测量截面意味着目前对次级中子的剂量估计可能被低估。
• 探测器开发： 结果对于金刚石中子探测器的设计和校准至关重要，这些探测器依赖于对快中子范围内反应阈值和截面的精确了解。
• 天体物理与聚变： 改进的数据有助于聚变设施（如 IFMIF-DONES）和空间剂量计的屏蔽计算。

总之，这项工作提供了一个决定性的实验基准，挑战了现有的核数据库，并强调了将高精度微分测量与先进理论建模（TALYS）相结合，以解决轻核反应物理中长期存在的差异的重要性。