十三种锕系核素三能区辐射中子俘获截面的经验与半经验参数化：从EXFOR原始实验记录到数据集级不确定性评估

《Symmetry》：Empirical and Semi-Empirical Parameterization of Radiative Neutron Capture Cross Sections for Thirteen Actinides Across Three Energy Regimes

　　

研究人员提出一种仅基于实验数据的系统性参数化，用于描述对反应堆物理、燃料循环建模和嬗变重要的十三种锕系核素（²³²Th、²³³U、²³⁵U、²³⁸U、²³⁷Np、²³⁹Pu、²⁴⁰Pu、²⁴¹Pu、²⁴²Pu、²⁴¹Am、²⁴³Am、²⁴⁴Cm和²⁴⁶Cm）的辐射中子俘获截面σ_n,γ(E)。研究人员将90,299个EXFOR数据点约化为逐核素表格形式，并按照标准核数据约定划分的三个能区进行分析：热区（E < 0.5 eV）、共振区（分为RRR，0.5 eV–1 keV和URR，1–100 keV）以及快区（0.1–1.0 MeV）。在每个能区中，经验描述和半经验描述均在共同点不确定度模型下拟合同一数据；数据集整体归一化偏移的影响则通过引入逐数据集对数归一化参数后的重新拟合单独评估。模型通过约化χ²、AIC/BIC以及在URR中的数据集级留出验证进行比较，每个导出量均同时给出delta法不确定度和数据集级重采样95%区间。在热区，束缚能级加第一共振形式仅凭截面数据即可将²³⁹Pu、²⁴¹Pu、²³⁷Np、²⁴¹Am和²⁴³Am的第一共振能量恢复到优于3%；在URR，经Porter–Thomas修正的统计模型在AIC上对12个核素中的7个更优，对四个易裂变靶包括裂变竞争，但在数据集级留出验证中仅对5个中的3个更优；在快区，此处采用的闭式黑核统计模型对²³²Th、²³⁸U和²³⁷Np失效，而默认TALYS计算在一次归一化后重现²³²Th和²³⁸U的能量依赖，表明失效源于闭式模型近似而非统计模型理论。数据集级重采样区间最宽可达delta法不确定度的80倍，说明实验间差异而非计数统计主导了合并锕系俘获记录的不确定度。

研究背景与问题：锕系核素的辐射中子俘获决定辐照燃料中高锕系核素累积、热中子与快中子系统中子经济性、次锕系核素嬗变方案设计以及保障与燃耗测量解释。评价库虽给出受实验约束的推荐截面，但其实验记录来自不同实验室、技术、归一化、分辨率和年代，且核素间能区覆盖极不均匀。EXFOR（实验核反应数据交换格式）保留每个原始数据集标识，使合并记录既可按大样本分析，也可按相关实验活动集合分析；后者在物理上更正确，因为同一实验活动的数据点共享通量归一化、样品质量和探测效率，并非独立。研究人员试图回答一个刻意狭窄的问题：仅用透明、低参数公式，能从原始实验记录中提取什么，简单物理在何处失效。评价库输入被排除在拟合之外，仅作为后验基准；经验式提供模型无关参照，半经验式保留少量自由核参数，若模型正确，其拟合参数应重现独立列表核数据，这是纯经验拟合无法提供的可证伪检验。研究保留十三种核素，并因数据可用性排除²⁴⁵Cm，其EXFOR (n,γ)记录仅两个热区数据点。研究意义在于建立三能区一致参数化，量化数据集级系统差异对不确定度的支配作用，并为核数据评价与测量优先级提供诊断。

关键技术方法（不超过250字）：研究人员从IAEA EXFOR检索并编译90,299个(n,γ)点，其中48,276个带报告不确定度；保留EXFOR数据集标识（DSID）以支持数据集级留出和重采样。按热区E < 0.5 eV、RRR（resolved resonance region，可分辨共振区）0.5 eV–1 keV、URR（unresolved resonance region，不可分辨共振区）1–100 keV、快区0.1–1.0 MeV分治；在共同点不确定度模型下拟合经验与半经验式，缺失不确定度采用10%替代并做5%、20%和能区中位数敏感性检查；以每数据集对数归一化参数评估整体归一化稳健性；用约化χ²、AIC（Akaike information criterion，赤池信息准则）/BIC（Bayesian information criterion，贝叶斯信息准则）和URR数据集级留出验证比较模型；对1/E加权等lethargy分组平均；采用Porter–Thomas宽度涨落修正统计模型并对易裂变靶加入裂变竞争；以delta法和以EXFOR数据集为簇的cluster bootstrap（B = 400）给出不确定度；快区与默认TALYS-2.24计算对比。样本队列来源为EXFOR数据库。

研究结果：3.1 热区（E < 0.5 eV）。3.1.1 经验与半经验描述：经验式为幂律σ_γ(E) = A·E^?b，b = 1/2时给出理想1/v行为，并解析导出Westcott因子g_W；半经验式为束缚能级加第一正能共振的Breit–Wigner形式，五个参数在不提供列表共振参数下拟合。3.1.2 结果与共振参数恢复：七个核素可拟合，BIC与AIC选择一致；²³²Th和²³³U因最低共振远高于热区，幂律更优；五个低共振核素半经验式强优，χ²_ν降低2.4倍（²⁴³Am）至450倍（²³⁹Pu）；第一共振能量恢复全程优于3%，五个中四个优于1%，拟合宽度44–112 meV量级正确；²³²Th和²³³U多起点拟合正确塌缩至束缚能级项。3.1.3 热截面、Westcott因子与共振积分：0.0253 eV截面与IAEA参考比较，12个可比值中9个在3%内、5个在1%内；²⁴¹Am偏高19%但重采样不显著，²⁴³Am偏高16%且重采样后仍真实，²³³U偏高8%仅因七个热区点。3.1.4 稳健性与不确定度量化：激活逐数据集归一化后σ₀移动小于2.1%、共振能量移动小于0.1%；重采样区间可比delta法宽至80倍；²⁴¹Am重采样区间526–702 b包含587 ± 12 b，²⁴³Am区间84.3–88.6 b排除75.1 ± 1.8 b；²³⁵U和²³⁸U热拟合在窄能 span 不可识别；²³⁹Pu和²⁴¹Pu重采样退化为两个数据集。3.2 共振区。3.2.1 分组平均与RRR/URR划分：0.5 eV–100 keV按1/E权重、25 bins/decade等lethargy分组，部分共振积分由组平均求和。3.2.2 RRR：经验幂律拟合组平均，十二个中十个R²为负，相对RMSE在十一个中为107–677%，说明RRR不宜用单一平滑幂律，brrr仅作粗斜率指标。3.2.3 URR：统计模型保留s、p波并作Porter–Thomas宽度涨落修正，W为0.689–0.965；易裂变靶加入裂变竞争后，AIC偏好仅²³³U变为经验式（ΔAIC = 4.6），²³⁵U经验式仍优但裕度减小（ΔAIC = 0.7），²³⁹Pu和²⁴¹Pu统计模型仍优（ΔAIC = 10.9和5.6）；统计模型在12个中7个获AIC/BIC支持。3.2.4 物理有效性与RIPL-3基准：宽度涨落修正解除²³⁹Pu的S₀边界限制（无裂变S₀ = (1.40 ± 0.28)×10^?4，含裂变S₀ = (3.70 ± 1.17)×10^?4，约为RIPL-3的(1.20 ± 0.10)×10^?4三倍）；²³⁷Np由二倍亏缺改善至5%一致但样本内拟合变差；²³²Th和²³⁸U的S₀仅升13–14%，不足以补四倍差距；R_γ呈系统性奇偶特征，偶偶靶拟合/参考比约0.98–2.12，奇A靶含裂变拟合为0.042–0.382。3.2.5 数据集级留出验证：仅五个核素，3/5不具推断权重，但信息准则排序未在样本外重现；²³⁸U统计模型虽被AIC/BIC强优，其留出RMSLE（root-mean-square logarithmic error，均方根对数误差）反比经验式差21%。3.2.6 ²³²Th和²³⁸U强度函数与数据集归一化：其S₀低于RIPL-3三至四倍，远超5–15%可信归一化误差，表明s + p统计核存在真实局限，需完整Hauser–Feshbach处理。3.3 快区（0.1–1.0 MeV）。3.3.1 范围、筛查与经验描述：仅²³²Th、²³⁸U、²³⁷Np、²³⁹Pu、²⁴¹Am有可用数据，剔除三点；拟合对数二次和指数式，R²仅0.31–0.72，方差分解的组内相关ρ_I ≈ 0.5，约一半未解释方差为实验活动间差异。3.3.2 闭式Hauser–Feshbach模型失败：黑核形成截面、背移Fermi气能级密度、Lorentzian E1 γ强度函数和自由a、C的闭式模型对²³²Th、²³⁸U、²³⁷Np给出负R²（?1.14、?1.04、?0.41），a/A达下界0.034而既定锕系值约0.11 MeV^?1；默认TALYS-2.24无拟合参数对²³²Th R² = ?0.25，单次归一化后对²³²Th、²³⁸U、²³⁹Pu的R²为0.32、0.32、0.57，幅度偏低1.19、1.64、2.08倍，²³⁷Np仍不被重现（R² = ?0.05）。3.3.3 归一化处理及诊断价值：拟合逐数据集归一化因子在快区跨因子二，²³²Th 18/24和²³⁸U 24/34数据集尺度偏移超过5%，许多达0.5–0.6，远超可信实验归一化误差，χ²_ν仍为98和144，说明闭式模型失败被归一化自由度强化。3.4 锕系链系统学：36项检验预期约1.8个假阳性，无一项通过Bonferroni校正（α = 0.0014）；候选关系包括log₁₀σ₀与S_n、σ₀与奇偶性、brrr与A及(N ? Z)/A、部分RI与A，均仅作筛选。中子–质子不对称I = (N ? Z)/A与R_γ基准比较为描述性诊断，不证明对称破缺。3.5 三能区拟合质量：热区半经验式对²³⁹Pu、²⁴¹Pu、²³⁷Np、²⁴¹Am使RMSE降低5–53倍，R²达0.885–0.999，²⁴³Am因归一化不一致仍不满意；RRR为核心负结果，十一个中相对RMSE 107–677%、十个中R²为负；URR相对RMSE 6–41%、十个中R² > 0.90，宽度涨落修正降低八个核素RMSE；快区两种经验式几乎不可区分，闭式模型对²³²Th、²³⁸U、²³⁷Np差25–102%且R²为负，相对RMSE 7–41%为三区中最弱约束。

讨论部分总结：讨论指出，热区是半经验模型明确改进之处，低共振存在时仅凭截面数据恢复共振能量优于3%，无低共振时正确塌缩至束缚能级项；URR部分有效，AIC对12个中7个支持而数据集级留出仅5个中3个支持，二者差异本身说明AIC优势中有多少来自物理、多少来自吸收数据集系统学；RRR和快区不适用，RRR组平均后仍无幂律描述，快区闭式黑核近似即使允许24–34个自由归一化也不能描述数据，而默认TALYS经单次归一化可重现部分核素能量依赖，故缺陷在近似而非统计模型理论；²⁴⁵Cm因仅两个热区点被排除，是重要测量缺口；R_γ奇偶特征指向两个J通道塌缩为一个有效?D?是剩余主要模型缺陷，RIPL-3的?D₀?可约束R_γ并打破S₀–R_γ退化；delta法与重采样之间80倍差距表明合并EXFOR记录不确定度由独立实验活动系统差异主导，导出量应附数据集级重采样区间，少于三个独立数据集的核素应标记为无有意义的经验不确定度。局限包括未按实验来源加权、单?近似使URR偏差可达4–5倍、d波项近似、裂变竞争仅对易裂变靶以固定?Γ_f?/?Γ_γ?纳入、省略Moldauer弹性增强项、归一化参数可吸收模型误差、重采样需至少三个独立数据集、十三种核素不足以确立通过多重比较校正的系统学。方法扩展包括联合估计缺失不确定度尺度、按年代和技术做交叉方差分解、重积分评价库以进行加权比较、用RIPL-3约束的J相关能级间距做显式自旋加权、d波扩展、快区TALYS/EMPIRE完整反应码计算以及更多独立实验活动。

研究结论翻译：该分析基于90,299个EXFOR数据点并仅拟合实验数据，为十三种锕系核素建立了三能区一致的经验与半经验(n,γ)截面参数化。作为外部验证，导出的热截面与IAEA参考汇编在十二个有参考值核素中的九个内一致于3%。在热区，半经验束缚能级加共振模型仅凭截面数据将²³⁹Pu、²⁴¹Pu、²³⁷Np、²⁴¹Am和²⁴³Am的第一共振能量恢复到优于3%，相对幂律基线使χ²_ν降低2.4至450倍；闭式Westcott因子g_W = (2/√π)·Γ(2 ? b)由经验幂律解析得到。可分辨共振区不存在此类幂律描述：即使作1/E组平均，十二个核素中十个的R²仍为负。在URR，经宽度涨落修正的双分波统计模型对易裂变靶包括裂变竞争，在AIC上对十二个中七个更优，但数据集级留出仅五个中三个更优；这一差距表明信息准则在异质实验数据上可能高估物理模型。Porter–Thomas修正解除了²³⁹Pu受边界限制的强度函数（AIC 120.9→65.5），并使²³⁷Np与RIPL-3一致至5%，但仅使²³²Th和²³⁸U的S₀升高13–14%，不能解释其四倍亏缺。拟合R_γ呈系统性奇偶模式，且在对易裂变靶纳入裂变竞争后仍存在，表明单J近似是剩余主要模型缺陷；按表8对易裂变靶纳入裂变竞争后，其S₀为RIPL-3的1.4–3.1倍，因此²³⁹Pu的接近一致仅在省略裂变时获得。在快区，此处闭式黑核模型不足：黑核近似使能级密度参数偏离二至三倍，²³²Th和²³⁸U的χ²_ν仍为98和144；默认TALYS计算重现²³²Th、²³⁸U和²³⁹Pu能量依赖但幅度低估1.2–2.1倍，故失败归因于闭式模型近似，而绝对预测仍需调整参数；约一半快区未解释方差来自实验活动间差异。所报锕系链系统学未通过多重比较校正，仅具假设生成意义。核心定量结果来自不确定度比较：数据集级重采样区间最宽可达delta法不确定度的80倍，表明锕系俘获语料不确定度压倒性由独立测量活动间系统差异而非计数统计主导。因此主要实践建议是：任何从合并EXFOR数据提取的锕系俘获导出量，若未附数据集级重采样区间，不应发表。