Double Chooz发表首个残余反中微子通量定量测量结果

停堆的核反应堆并不沉寂。堆芯和燃料池中积累的裂变产物会在数月乃至数年间继续衰变，其贝塔衰变链释放反中微子。Double Chooz合作组发表了这一残余通量的首个定量测量结果，测量利用了Chooz核电站两座反应堆同时停堆的时段。

来源：journals.aps.org

通俗解读

木火熄灭后，余烬仍会长时间泛红。核反应堆也有类似现象，只是它发出的微光肉眼无法看见：那是一股幽灵般的粒子流——反中微子，由留在堆芯内的辐照燃料以及储存在水池中的乏燃料组件内裂变产物的贝塔衰变链产生。

在法国阿登地区Chooz两座反应堆停堆期间的17.2天有效测量时间内，近探测器——距堆芯约400米，并有115米水当量的屏蔽——在扣除本底后分离出106 ± 18个残余候选事例，约合每天六个。详细模型预测，44%的信号来自冷却乏燃料的水池；探测器并未直接分离这一比例。

这项结果设想的实际用途是核材料核查：反中微子穿过物质时发生相互作用的概率极低，因此可以在不进入堆芯的情况下观察反应堆。不过必须注意：这一结果依赖数十立方米闪烁体和115米水当量的地下屏蔽。它验证了残余通量计算；并未证明便携式设备能在核电站地表运行。

技术解析

1. 如何捕捉反中微子，以及为何需要大量事例。 所用反应是逆贝塔衰变（IBD）：$\bar{\nu}_e + p \rightarrow e^+ + n$νˉe​+p→e++n。电子反中微子与自由质子——即闪烁体中的氢原子核——相互作用，产生一个中子和一个正电子。信号由两部分组成，正是这一点使得每天只有少数事例时仍能进行测量。正电子立即沉积动能，随后湮灭：这构成瞬发信号。中子散射后被俘获；Double Chooz的“全中子俘获”方法利用内部体积中钆、氢和碳上的俘获形成延迟信号。要求两种信号出现在0.5至800 µs的时间窗内，且相距小于1.2 m，可以排除大部分环境放射性。反中微子的能量由瞬发信号按$E_{\bar{\nu}_e} \simeq E_{\text{prompt}} + 0.78\ \text{MeV}$Eνˉe​​≃Eprompt​+0.78 MeV重建，而该反应的运动学阈值约为1.8 MeV：低于这一阈值时，该通道无法探测反中微子。

2. 残余微光从何而来：三条链，三只时钟。 反应堆运行时，其反中微子通量主要来自短寿命裂变产物；停堆后的数小时内，这一贡献会迅速下降。可探测的剩余信号主要来自三组母核—子核对：在辐照及随后的冷却过程中，长寿命母核持续产生寿命短得多的子核。子核每次产生后都会迅速消失，因此其活度主要跟随母核；停堆首先去除了原本遮蔽这些持久衰变链的主导通量。

   核素对	母核半衰期	子核：半衰期及最大能量$Q_\beta$Qβ​	作用
   $^{144}\mathrm{Ce} \rightarrow {}^{144}\mathrm{Pr}$144Ce→144Pr	285天	17.3 min · 3.00 MeV	约占54%的模型预期1–3 MeV信号
   $^{106}\mathrm{Ru} \rightarrow {}^{106}\mathrm{Rh}$106Ru→106Rh	372天	30.1 s · 3.54 MeV	约占模型预期信号的38%
   $^{90}\mathrm{Sr} \rightarrow {}^{90}\mathrm{Y}$90Sr→90Y	28.9年	3.19 h · 2.28 MeV	在模型的长期预测中成为主导（十年后> 90%）

   这些能量决定了其余一切：$Q_\beta$Qβ​为2至3.5 MeV，而IBD阈值为1.8 MeV，因此每个贝塔能谱中只有高能尾部可被探测。作者据此预计，可探测残余通量的98.7%位于3 MeV以下——这既说明为何要把分析限制在1–3 MeV的瞬发能量窗口，也解释了测量为何困难：这里的本底最为密集。

3. 1–3 MeV窗口究竟测量什么，以及验证它的对照。 把分析限制在预期出现信号的窗口内，会引出一个经典质疑：只在选择寻找的地方找到预期结果。作者以负对照回应。3 MeV以上的产额在扣除本底后与零相容。模型确实预测有一条小的残余尾部，但每天低于0.1个事例，因而低于当前灵敏度。该区域没有可测超额，这支持本底模型；却不能证明其中的物理通量严格为零。

4. 本底，以及地下屏蔽为何重要。 每天只有六个事例时，任何因素都不可忽略。三类本底分别处理。偶然符合——两次互不相关的能量沉积，一次被当作瞬发信号，另一次被当作延迟信号——由抑制能力超过400的神经网络去除，残余部分再根据孤立能量沉积率进行原位测量。相关本底来自宇宙线缪子：一种是岩石内散裂产生的快中子，它先模拟瞬发信号再被俘获；另一种是$^{9}\mathrm{Li}$9Li之类的β-n同位素，由缪子在闪烁体的$^{12}\mathrm{C}$12C上产生散裂，其衰变本身就能给出电子—中子对。后者通过与母缪子的时空相关性加以标记。近探测器的115 m水当量屏蔽降低了缪子通量，从而减少这些同位素的产生；任何否决系统都无法将其彻底消除。

5. 预测：逐个燃料组件进行的盘存计算。 预期数目不是数量级估算，而是详细模拟。预测的IBD能谱写作
   $\frac{dN_{\bar{\nu}_e}}{dE_{\text{vis}}} = \sum_i \frac{N_p}{4\pi L_i^2}\iint \varphi_i(E,t)\,\sigma_{\text{IBD}}(E)\,P_{ee}(E,L_i)\,\mathcal{D}(E_{\text{vis}}|E)\,dE\,dt,$dEvis​dNνˉe​​​=i∑​4πLi2​Np​​∬φi​(E,t)σIBD​(E)Pee​(E,Li​)D(Evis​∣E)dEdt,
   其中$N_p$Np​是靶质子数，$L_i$Li​是到源$i$i的距离，$\varphi_i$φi​是随时间变化的发射能谱，$\sigma_{\text{IBD}}$σIBD​是截面，$P_{ee}$Pee​是振荡后的存活概率，$\mathcal{D}$D是探测器响应。裂变产物活度由APOLLO-2.8.4和DARWIN-3代码结合JEFF-3.1.1核数据库计算；贝塔能谱来自BESTIOLE库，该库采用费米理论，并纳入电磁修正、有限核尺寸、原子屏蔽以及允许跃迁和禁戒跃迁的形状因子。每个燃料组件的辐照历史均单独模拟，其整个寿命期的热功率一直跟踪到每段停堆期开始。

6. 模型为何将44%归因于燃料池。 两座储存池靠近堆芯，因此从约400 m外的探测器看，其几何因子$1/L^2$1/L2仍与堆芯相近。每个组件经过多个辐照循环后进入水池，并在那里冷却多年。停堆堆芯内是快速衰变的较新盘存；燃料池则累积了不同冷却年龄的组件，其中一些同位素可长期存在。停堆时间延长时，模型预测的燃料池贡献总和会接近堆芯贡献。积分计算给出堆芯56%、燃料池44%；现有几何布局无法分别测量这两个比例。

7. 预测受何限制，以及不能由此得出什么结论。 预测事例数的归一化不确定度达到7.4%，其中谱建模为6.0%，基线几何为2.9%，裂变产物盘存为2.1%，燃料池内组件盘存为2.0%，探测器贡献低于0.8%。主导项来自$^{144}\mathrm{Pr}$144Pr贝塔跃迁的形状，其不确定度取详细核结构计算与所谓$\xi$ξ近似之差。这7.4%既不是数据—模型符合程度的精度，也不是盘存推断的性能指标：本研究只是把计数与预测相比较，并未反向估计燃料质量或组成。

图示符号：$N_{\rm obs}$Nobs​为测得的残余计数，$N_{\rm pred}$Npred​为模型预测，$H_0(B)$H0​(B)为仅有本底的假设。$5.9\,\sigma$5.9σ比较的是$N_{\rm obs}$Nobs​与$H_0(B)$H0​(B)，而不是$N_{\rm pred}$Npred​。

为何能够实现

测量成功依赖一个时机和一种几何条件。时机方面：2017年，Chooz两座反应堆曾四次同时停堆，其中一次持续20.8天——这是不可缺少的条件，因为满功率时近探测器每天记录约900个IBD事例，而这里每天只有6.2 ± 1.1个残余事例。几何方面：探测器已经建成、校准并得到充分理解，安装在400 m外；多年的振荡测量已经确定了它的系统效应。

数值比较很直接。此前对2011–2012年数据的分析只收集到约二十个候选事例——不足以提取残余能谱。在这里，扣除本底后的106 ± 18个残余候选事例与预期的88 ± 7个相比较：差值为$\Delta = 18 \pm 19$Δ=18±19个事例，即两者相符。远探测器测得27 ± 13个，而预期为14 ± 1个；误差条较大，结论也宽松得多。

局限十分明确。信号分析覆盖1–3 MeV窗口。3 MeV以上，扣除本底后的产额与零相容；模型则保留一条延伸至约4.5 MeV的小残余尾部，其每天低于0.1个事例，本样本无法分辨。7.4%的归一化不确定度描述的是预测；实验并未证明盘存推断的精度。

必须明确区分已经证明与尚在设想的内容。已经证明的是：残余通量以$5.9\,\sigma$5.9σ与本底相区分，且计数$106\pm18$106±18与预测$88\pm7$88±7相容，因为$\Delta=18\pm19$Δ=18±19。尚未证明的是：燃料的运行级监测。这项验证使用30 m³以上的闪烁体，置于115 m水当量屏蔽之下，在400 m处每天只获得约六个事例。小型地表仪器仍需控制宇宙线本底，并证明其对盘存变化的灵敏度。

因果链

堆芯内发生裂变 → 辐照期间积累长寿命母核（$^{144}\mathrm{Ce}$144Ce、$^{106}\mathrm{Ru}$106Ru、$^{90}\mathrm{Sr}$90Sr），并由其产生短寿命贝塔子核 → 反应堆停堆，主导通量迅速消失 → 已在堆芯和转移至燃料池的组件中形成的衰变链继续存在 → 四次同时停堆合计24.4天，发生在2017年 → 近探测器有效测量17.2天 → 得到$106\pm18$106±18个残余候选事例，位于1至3 MeV之间，与仅有本底的假设在$5.9\,\sigma$5.9σ水平不相容 → 3 MeV以上预测有一条小尾部，每天低于0.1个事例且未能分辨 → 与逐组件模型相容（$88\pm7$88±7，$\Delta=18\pm19$Δ=18±19） → 首个检验残余通量计算的定量基准 → 未来研究分离堆芯与燃料池贡献，并开发专用探测器。

一则轶事

Double Chooz曾利用2011–2012年的数据尝试过这项测量。当时合作组只收集到约二十个残余候选事例：近探测器尚未投入使用，两座反应堆同时停堆的时间也太短，无法提取能谱。2017年近探测器已经运行，四次同时停堆中有一次持续20.8天，共提供17.2天有效测量时间；这两项变化使定量测量成为可能。

影响与当前证据

该测量首次公开提供了残余通量求和计算与实际计数之间的定量比较。预测的归一化不确定度为7.4%；观测比较则为$\Delta=18\pm19$Δ=18±19个事例。两者不可混淆。

在应用层面，预测的来源分配是最具直接用途的新信息：模型将停堆场址信号的44%归因于燃料池，而非堆芯。因此，只模拟仍在堆芯中的燃料会使监测方法严重低估预期信号。不过，论文无法量化小型地表探测器能达到的灵敏度：文中没有此类数据，作者也没有从地下装置的活性体积作出外推。

研究者视角——开放问题

以下是可能的延伸方向，并非本研究的结果。

• 决定性的能谱检验。 若积累更多统计量，在数月内跟踪$^{144}\mathrm{Ce}$144Ce与$^{106}\mathrm{Ru}$106Ru的相对演化，可以检验模型的时间形式。要观察$^{90}\mathrm{Sr}$90Sr–$^{90}\mathrm{Y}$90Y核素对占主导则需要多年：预测只有在超过十年后才高于90%。
• 利用几何布局分离燃料池与堆芯。 此处56%/44%是模型预测，而不是测得的分离结果：在400 m处，38 m的距离差无法分辨。把探测器放得近得多，或在不同角度布置两个探测器，可以通过实验区分两种贡献，而不是依靠计算分配。
• 核结构瓶颈。 主导不确定度并非来自探测器，而是$^{144}\mathrm{Pr}$144Pr贝塔跃迁的形状。专门测量该贝塔能谱可以降低6.0%这一项，但研究没有量化届时模型最终能达到的精度。

来源

以下参考资料已在八月16日（2026年）的事实核查审计中验证。

• Double Chooz Collaboration，“First Measurement of Neutrino Emissions from Spent Nuclear Fuel by the Double Chooz Experiment”，Physical Review Letters 137，061803，八月4日（2026年）——同行评议论文。DOI：10.1103/dr26-j19g
• 同一工作的预印本版本：arXiv:2510.04869——已查阅全文（摘要、方法、系统不确定度表）。

背景参考

1. P. Vogel与J. F. Beacom——Angular distribution of neutron inverse beta decay——逆贝塔衰变的参考运动学和截面。

置信度声明

对从同行评议论文中复现的计数、显著性、能量窗口、几何参数和不确定度分量具有高置信度。对应用范围具有中等置信度：论文讨论核保障监督，但既未测试小型地表探测器，也未测试盘存推断。归因于燃料池的44%是模型输出，而非测得的分离结果。这项首个定量分析尚无独立重复。唯一性检查：登记表A001–A018中没有相同的DOI、arXiv标识符或等价的核心结果；结论为UNIQUE。DOI、日期、公式、链接及无重复均已核验。