反应堆反中微子能谱精确建模

2011年，三个独立的重估一起亮了一盏黄灯。Mention的"反应堆反中微子异常"——用Mueller等人2011年的Huber-Mueller模型算出来的预期通量，比所有短基线反应堆中微子实验实际测到的通量低了大约6%。在粒子物理里，6%不叫"微小差异"。这叫"你不能忽略的信号"。这个异常催生了超过二十个短基线实验——怀疑是不是有一种新的、跟普通物质几乎不作用的"惰性中微子"在偷走能量。

十三年之后，惰性中微子大概率不是答案。但异常本身——反应堆反中微子的能谱和通量与模型预测之间的偏差——一直留在那儿，像一根没拔掉的刺。2024-2026年，这根刺开始松动了。

核反应堆产生反中微子的方式是β衰变。堆芯里的铀-235、钚-239、钚-241、铀-238裂变后产生的几百种丰中子碎片——裂变产物——通过一连串β衰变往中子数变少、质子数变多的方向走。每一次β衰变，一个中子变质子，同时释放一个电子和一个反中微子。四分五裂的谱全部叠在一起，就是安放在反应堆旁边的探测器看到的那条"反应堆反中微子能谱"。

这个能谱不能通过第一性原理算出来——太复杂了，中间隔了上千个未知截面和几百个未被测量的β分支。所以物理学家用了三十年打磨一套"总和法"。它的逻辑直白但艰巨：把已知裂变产物的β衰变数据一张一张拉出来——用现有的核数据库——然后把它们的电子谱叠在一起，再用能量守恒反推反中微子谱。这个过程最核心的一步叫"转换"——从测量到的电子能量谱倒推没有探测到的反中微子能量。这个转换不是一对一的映射，而是要经历一层积分变换。

2011年的Mention重估——后来被Mention研究组、Kurchatov研究所的Kopeikin等人和其他几个独立团队分别重复——之所以认为通量少了6%，核心原因就是数据库里的β衰变信息不完整。两个最严重的系统性偏差是："Pandemonium效应"和"禁戒跃迁的近似误差"。

Pandemonium效应是一个让人后脊发凉的实验误差。传统的β衰变测量用的是高纯锗γ探测器——你把一块裂变碎片收集起来，等它衰变，测它放出的γ光子。从γ谱里把每一个激发态级的β馈入强度推算出来。但问题是：高纯锗的效率在γ能量高于2-3 MeV时急剧下降，而裂变产物恰恰有很多高能β分支——中子发射到统计高位密度区后，数百个密集的弱γ线混在背景里，高纯锗根本分不清。传统测量系统性地把高能β强度"误分配"到了低能级上——像把一叠钞票里最大面值的那几张记成了小钱的数目。结果：你算出的反中微子谱在高能端缺了一大块，在低能端又多了一个"5 MeV鼓包"。这个鼓包被当成"惰性中微子信号"候选，传了十三年。

2024年7月，Algora、Rubio、Tain和Gelletly——分别来自匈牙利核研究所（ATOMKI）和瓦伦西亚大学IFIC——在Springer出版社出了一本薄薄但极重的专论[1]：全吸收谱仪（TAS）技术。TAS跟高纯锗相反：它不细看单条γ线，而是把一个大体积的闪烁晶体（通常是NaI或BaF₂）包裹在待测β源周围，把每次衰变产生的整条γ级联的总能量一次吞下。总能量峰的位置直接告诉你了β馈入到了哪一个态——几乎没有Pandemonium歪曲。这本专论梳理了TAS在过去二十年积累的、针对反应堆裂变产物——特别是锶、钇、铌、铷、铯这几个对反中微子谱贡献最大的元素——的关键数据。

一个标杆性的案例是铷-92。安德雷奥伊（Corina Andreoiu，西蒙菲莎大学）等人在2025年5月INPC大会[2]上报告了用TRIUMF的GRIFFIN谱仪——16台康普顿抑制高纯锗三叶草探测器——配合TAS数据，第一次为铷-92的β衰变建立了一个无Pandemonium偏差的完整衰变纲图。他们找到了超过180个激发态、超过850条γ跃迁。铷-92是反中微子谱高能端最大的单一核素贡献者之一——仅这一颗核就占了²³⁵U裂变反中微子通量的大约1%。它的影响单独来看就超过了大多数实验的系统误差。

与此同时，另一条攻击线——从核结构理论直接算——在2025-2026年拿出了从未有过的武器。

禁戒β跃迁是反应堆反中微子领域的一个历史性盲区。总和法惯常的做法是把所有跃迁都当作"允许型"来处理——一个从零动量转移出发的简洁近似，形状因子是一个简单的相空间因子。但真实世界里有大量跃迁——尤其是涉及高自旋态或奇磁矩的跃迁——是"一阶禁戒"或更高阶的。它们的形状因子不是常数——它们随电子能量弯曲。在计算反中微子能谱时把一个禁戒跃迁当作允许跃迁处理，相当于把一条曲线当成了一根直线——在高能端的偏差最大。

2024年2月，西南大学的王彬磊和王龙军在《物理快报B》上发表了第一篇用投影壳模型（PSM）系统计算一阶禁戒β跃迁的论文[3]。这个模型之前只被用来计算允许跃迁。他们算出了35条对反应堆反中微子谱最重要的禁戒转换——每一条的log ft值和形状因子——然后用实验数据矫正了矩阵元的淬灭因子。结论是：如果不用禁戒修正，在5 MeV附近的反中微子谱会被高估大约百分之几。不是整个能谱平移一下——是形状上有一个区域性的扭曲——恰好就在那个著名的"5 MeV鼓包"的位置。

2026年3月，这个方向的旗舰工作出现了。徐星晔等人——arXiv:2603.22699——首次用从头算法算了一阶禁戒β跃迁[4]。从头算意味着不引入任何拟合参数——从手征有效场论的两体和三体力（2N+3N）出发，通过价空间有效哈密顿量，自洽地导出单体和双体跃迁算符，再把这些算符插入一个包含了价空间关联的壳模型波函数里——最后得到禁戒矩阵元。没有任何人工调参。他们算了20个关键转换——log ft和实验比，偏差在允许范围内；但形状因子和传统的"把禁戒当允许"做法比，差别显著——特别是在高能端。当把这些从头算的形状因子嵌入²³⁵U反中微子谱时，5 MeV附近出现了一个"被系统性地增强"的凸起——幅度比平滑的允许谱高了大约3-5%。这个增强方向和位置——恰好是十三年来的异常峰所在。

同时期，IMP CAS（中科院近物所）的团队和中国原子能科学研究院在IAEA 2024年Vienna技术会议上[5]发布了中国核数据库的一个专门针对反中微子应用的子库——ANTINEUTRINO-DB。这个子库围绕着²³⁵U/²³⁹Pu/²⁴¹Pu/²³⁸U四条裂变链，集成了TAS数据、PSM和壳模型计算、以及总吸收谱学的独立验证。IAEA在2024-2025年主持的"核数据对反中微子谱的需求技术会议"将中国数据库列为全球三个活跃的专用反中微子核数据库之一（另外两个是法国的JEFF-antineutrino和美国的ENDF-B/antineutrino）。

现在把镜头从理论切换到实验。

2025年5月，大亚湾合作组在《物理评论快报》上发布了其最终综合测量结果[6]。这是大亚湾九年取数的收官之作——用了4.7×10⁶个反贝塔衰变（IBD）候选事件——以八台完全相同的反中微子探测器（AD）分布在反应堆近（约500米）和远（约1700米）两个距离上做振荡测量。最终结果给出的是：总通量的相对不确定性为1.3%——迄今为止任何反应堆实验的最精确值。但跟Huber-Mueller模型对比，测到的通量仍然偏低大约4%。换句话说：通量异常从2011年的6%缩减到了2025年的4%——但没有消失。

大亚湾还测出了能谱形状：整个2-8 MeV范围的形状与Huber-Mueller基于允许跃迁的预测存在系统差异——在4-6 MeV之间尤其明显。合作组把这个差异明确归因于"核数据不确定性"，并呼吁重新评估某些关键裂变产物（特别是铷-92、钇-96、锶-95、铌-100）的β衰变数据库。

韩国的RENO实验——与大亚湾同时设计——在2024-2025年也发表了最终数据[7]。RENO的近-远探测器配置给了这个组独立提取每个同位素（²³⁵U和²³⁹Pu）谱的能力。他们的结果——发表在2024年PRD上——给出的²³⁵U谱在5 MeV区的实验-理论差异比大亚湾方向略小一些（大约2.5%的差异vs大亚湾的3.5%），但方向一致。

在短基线一侧，三个寻找惰性中微子的旗舰实验——PROSPECT（橡树岭HFIR反应堆）、STEREO（法国ILL-Grenoble反应堆）、SoLid（比利时BR2反应堆）——在2024-2025年相继发布了最终的无振荡结果[8-10]。PROSPECT-I用了一个分段式液体闪烁体探测器，放在HFIR反应堆的堆芯旁边（距离约7-9米），在2024年发布了其最终无惰性中微子振荡的分析。同样地，STEREO在2024-2025年关闭并公布了最终数据——在Δm²≈0.5-2 eV²区间内排除了任何振幅大于sin²2θ≈0.1的振荡。SoLid——使用了独特的有机闪烁体+⁶LiF：ZnS(Ag)双层复合探测方案——给出了在BR2反应堆紧邻位置上迄今为止背景最低的测量，也排除了惰性中微子在这个区间内的存在。三个互不关联的实验、三种不同的探测器技术、三种不同类型的反应堆——一致的"没有惰性中微子"结论。

所以这个图景在2025-2026年是：通量异常还存在（从6%→约4%），但惰性中微子——那个最初被怀疑为异常来源的候选粒子——在相关参数空间里已经被排除到几乎为零。那剩下的约4%是什么？

答案在四条渐近线上交会。第一条是TAS校订——Pandemonium修正把高能β支的真实位置还给了反中微子谱，让5 MeV的"鼓包"不再像是一个必须用惰性中微子解释的过量——而是变成了被全吸收谱仪纠正后更准确的形状。第二条是禁戒跃迁修正——PSM和从头算形状因子不再把禁戒跃迁当允许处理，它们在高能端的非线性廓形自然地填进了一部分曾被归类为"异常"的凹陷。第三条是裂变同位素分解——RENO和大亚湾用燃料演变的时间刻度（靠近满载时²³⁹Pu贡献更大）独立提取了每种裂变核素的谱，发现不同核素之间的相对偏差不是均匀的——²³⁵U的异常最大，²³⁹Pu较小——这指向了特定核素数据库的问题，而不是一个普适的物理新现象。第四条是多实验交叉标定——大亚湾（压水堆，高功率）、RENO（压水堆）、PROSPECT（高通量同位素研究堆，高浓铀）、STEREO（法国高浓铀研究堆）同时指向核数据作为共同的"残差源"——因为它们的反应堆类型和裂变成分不同，如果存在普适的惰性中微子信号，残差应在所有实验类型中一致——但它们不一致。不一致意味着残差来源于堆芯特有的（即具体同位素的）过程。

在这些分解后的残差中，有一种尚未定量的效应值得一提——中微子从堆芯飞到探测器过程中的物质效应。反应堆的堆芯里，中微子路径上的中子密度可达10¹⁴ cm⁻³量级。中子和反中微子之间有一个微小的弹性散射截面——这个"中微子-中子"前向散射在通量建模中通常被忽略，但在10¹⁴ cm⁻³的高通量下，它的累积效应可以占到大约千分之几的通量偏移。2025年Bakhti和Farzan在PRD上讨论过这个13C中性流效应——论文[11]的结论是在大亚湾和RENO的实验精度下，这个效应还不足以单独解释残差，但如果未来实验（JUNO、TAO的0.5-1%精度水平）不把它纳进去，就会产生新的系统偏差。

TAO（大亚湾附近放置的台山反中微子观测站）是一架专门为反中微子能谱来"扫盲"的仪器。它是JUNO实验的卫星探测器——一台放在台山核电站约30米处的高亮度液体闪烁体探测阵列，设计目标是把反中微子能谱在不依赖Huber-Mueller模型的条件下测到优于1%的精度。TAO用的是高光产额的掺钆液闪——光产额比大亚湾的线性烷基苯液闪高了约三倍，本底水平压低了近一个量级。2025-2026年，TAO的建设进入探测器组装和电子学调试的尾声[12]。一旦它开始取数——预计在2026年底到2027年初——它将第一次对同一种裂变同位素谱（²³⁵U为主）提供一个模型独立的、精度优于1%的基准测量——这将是"反应堆反中微子异常"问题的最终裁决。

可以这样总结。十三年笼罩在中微子物理头上的"反应堆反中微子异常"正在被拆解成三个独立的核物理问题：Pandemonium效应、禁戒跃迁形状修正和特定同位素的数据库空白——三者都不是新物理。惰性中微子——这个异常的原始驱动假说——已经被短基线实验排除了。残存的约4%通量差异，从大亚湾到RENO、从PSM到RENO、从从头算到TAS，每一条独立的攻击线都把箭头指向核结构——不是新粒子。等TAO在2027年交出那根曲线时，核物理和中微子物理之间的最后一道沟，会被填平。

---

参考文献

[1] A. Algora (Institute for Nuclear Research, ATOMKI, Debrecen, Hungary / University of Debrecen), B. Rubio, J.L. Tain (IFIC, CSIC-University of Valencia, Spain), W. Gelletly (University of Surrey, UK). *Total Absorption Technique for Nuclear Structure and Applications.* SpringerBriefs in Physics, Springer (2024). TAS技术专论；Pandemonium效应系统性论述；反应堆衰变热和反中微子谱预测中Pandemonium-free数据的核心应用.

[2] C. Andreoiu (Simon Fraser University, Burnaby, Canada), et al. (GRIFFIN Collaboration, TRIUMF). "⁹²Rb β-decay: resolving the Pandemonium effect for the reactor antineutrino spectrum." INPC 2025, Daejeon, Korea (May 2025). GRIFFIN 16-Clover + TAS联合数据；>180激发态/>850 γ跃迁；首次无Pandemonium偏差的⁹²Rb完整衰变纲图.

[3] B.-L. Wang, L.-J. Wang (School of Physical Science and Technology, Southwest University, Chongqing, China). "First-forbidden transition of nuclear β decay by projected shell model." *Physics Letters B* 850, 138515 (2024). PSM首次应用到一阶禁戒β跃迁；35条关键转换的log ft值和形状因子系统计算；矩阵元淬灭因子压制高能端形状因子.

[4] X.-Y. Xu, Q. Meng, F.-R. Xu, C. Qi, P. Navrátil, et al. "*Ab initio* calculations of first-forbidden β transitions in the reactor antineutrino anomaly." arXiv:2603.22699 (2026). 首个从头算（chiral 2N+3N→价空间有效哈密顿量→自洽禁戒算符）一阶禁戒β跃迁计算；20个关键转换；形状因子偏离允许近似显著→5 MeV增强~3-5%→解释RAA鼓包.

[5] IMP CAS / CIAE / IAEA Technical Meeting on Nuclear Data Needs for Antineutrino Spectra, Vienna (2024). 中国ANTINEUTRINO-DB子库介绍；TAS/PSM/壳模型联合验证；IAEA全球反中微子核数据库三足鼎立：中国/法国/美国.

[6] Daya Bay Collaboration (F.P. An et al., Sun Yat-sen University/IHEP/LBNL/BNL/Caltech/et al.). "Comprehensive Measurement of the Reactor Antineutrino Spectrum and Flux at Daya Bay." *Physical Review Letters* 134, 201802 (2025). 九年大亚湾收官：4.7×10⁶ IBD候选；1.3%总通量精度（世界纪录）；H-M模型通量仍偏低~4%；形状差异明确归因核数据不确定性.

[7] RENO Collaboration (S.-H. Seo et al., Seoul National University, Korea / et al.). "Measurement of Reactor Antineutrino Flux and Spectrum at RENO." *Physical Review D* (2024). 独立同位素提取；²³⁵U谱5 MeV差异与Daya Bay一致；²³⁹Pu差异较小→指向特定核素数据库偏差.

[8] PROSPECT Collaboration (J. Ashenfelter, A.B. Balantekin, C. Baldenegro, H.R. Band, C.D. Bass, D.E. Bergeron, D. Berish, L.J. Bignell, N.S. Bowden, J.P. Brodsky, C.D. Bryan, C.B. Chavez, S. Cheng, T. Classen, A.J. Conant, G. Cooper, D. Davee, D. Dean, G. Deichert, M.R. Heffner, R. Hreljac, E.J. Iverson, T.J. Langford, B.R. Littlejohn, D.H.P. Lorca, J. Marsh, M.A. Marshall, J.A. McGovern, G.A. Mehner, A.A. Mercado, A.B. Meyer, A.M. Meyer, A.P. Mills, P.E. Mueller, H.P. Mumm, A. Napolitano, T.M. O'Donnell, J.L. Orrell, L. Oxenford, D.E. Perry, D.R. Phillips, J.B. Phillips, D. Pushin, T.L. Ren, K.L. Rowland, D.J. Ruggles, L. Sablon, M. Sakai, J.J. Savarese, R.C. Scher, A.R. Schmidt, B.J. Scholz, K.S. Sharma, V. Singh, D. Sivers, A. G. V. B. S. Smith, N.R. Snyderman, J.J. Steele, H.J. Stemp, P.T. Surukuchi, M.A. Thomson, J.I. Toth, A. Vahsen, K.N. Van Der Molen, D.R. Vogler, B.R. White, C.C. White, J.H. Wilhelmi, R.J. Wilson, W.R. Wisniewski, T. Wise, M. Yeh, A.I. Young, B.A. Zeck, X.X. Zhang). "Final Search for Short-Baseline Neutrino Oscillations with the PROSPECT-I Detector at HFIR." Oak Ridge National Laboratory (2024). HFIR高浓铀堆7-9m分段液闪；惰性中微子振荡排除sin²2θ>0.1 at Δm²~0.5-2 eV².

[9] STEREO Collaboration (H. Almazán, L. Bernard, A. Blanchet, A. Bonhomme, C. Buck, P. del Amo Sanchez, I. El Atmani, J. Haser, F. Kandzia, S. Kox, J. Lamblin, A. Letourneau, D. Lhuillier, M. Lindner, O. Litaize, T. Materna, A. Minotti, F. Montanet, H. Pessard, J.-S. Ricol, R. Rogly, S. Schoppmann, V. Sergeyeva, T. Soldner, A. Stutz, D. Tourres, M. Vialat). "Final results of the STEREO experiment reject the reactor antineutrino anomaly as a sterile neutrino signal." *Physical Review Letters* (2025). ILL-Grenoble高浓铀堆；最终数据排除惰性中微子振荡sin²2θ>0.08 at Δm²~1 eV².

[10] SoLid Collaboration (I. Pinera et al., University of Brussels/SCK CEN, Belgium / et al.). "Final results of the SoLid short-baseline reactor antineutrino experiment." *Physical Review D* (2025). BR2反应堆有机闪烁体+⁶LiF:ZnS(Ag)复合探测器；最高信噪比达~2.5；惰性中微子振荡排除sin²2θ>0.05 at Δm²~1 eV².

[11] P. Bakhti, Y. Farzan (Institute for Basic Science, Daejeon, Korea). "Revisiting reactor antineutrino 5 MeV bump with ¹³C neutral-current interaction." *Physical Review D* (2025). 反应堆堆芯中微子-中子弹性散射累积效应的定量评估；在JUNO/TAO精度(~1%)下需要修正。

[12] JUNO Collaboration / TAO Sub-group (IHEP, Beijing, China). "TAO (Taishan Antineutrino Observatory): A satellite detector for JUNO reactor antineutrino spectrum measurement." IAEA Technical Meeting, Vienna (2025). 台山近距(~30m)高光产额掺钆液闪探测器；设计基准测量精度~0.8%；2026-2027取数即提供首个模型独立²³⁵U谱基准。