## 中文标题：解析核B样条解谱法用于反冲质子望远镜快中子谱恢复：蒙特卡罗响应矩阵的一种确定性替代方案

《Radiation Measurements》：Analytical-kernel B-spline unfolding of fast-neutron spectra from a proton recoil telescope: a deterministic alternative to Monte Carlo response matrices.

　　

## 反冲质子望远镜（Proton Recoil Telescope, PRT）的闭式解析响应核替代蒙特卡罗方法构建响应矩阵；采用带吉洪诺夫罚项与L曲线（L-curve）正则化参数自动选取的B样条解谱方法，并精确施加非负性约束；在2.45 MeV和14.1 MeV能量点基于Geant4的模型外（out-of-model）测试中，中子线位恢复精度优于4%；在匹配算子（matched-operator）合成测试中，分离良好的峰位定位精度优于0.15 MeV；基于泊松加权的确定性计算流程可在标准工作站上于数秒内完成全部解谱分析。

## 论文解读：解析核B样条解谱法实现快中子谱的快速确定性反演

### 一、研究背景与问题

快中子能谱测量在聚变等离子体诊断、辐射防护和核物理实验中具有核心地位。反冲质子望远镜（PRT）通过测量中子与含氢辐射体发生弹性散射所产生的反冲质子来反推中子能谱，其测量原理可追溯至Bernardini和Oppenheimer的早期工作。反冲质子能谱P(E_p)与待求中子能谱N(E_n)之间满足第一类Fredholm积分方程，其中响应核n(E_p, E_n)描述了入射中子能量与探测质子能量之间的映射关系。求解该方程属于典型的不适定逆问题（ill-posed inverse problem），需要同时提供两条关键要素：精确的响应核（response kernel）与能在噪声存在下稳定反演的解谱算法。

传统方法中，响应矩阵的构建高度依赖蒙特卡罗输运模拟（如Geant4、MCNP），其计算成本高昂（耗时数小时至数天）且结果带有统计涨落；而解谱算法方面，尽管已有最小二乘、最大熵、期望最大化（EM）及贝叶斯方法等多种方案，但往往面临正则化参数选择困难、计算开销大或可能产生负值非物理解等问题。Trunec等人曾提出基于B样条展开的解谱框架，但响应矩阵仍依赖蒙特卡罗模拟。Cordella等人则推导了等半径、共轴PRT响应的闭式解析表达式，采用椭圆积分形式并包含Bragg-Kleeman质子能损关系，与Geant4基准符合良好。将上述两项工作相结合，即可实现无需蒙特卡罗输入的解析响应矩阵构建，这正是本文研究的出发点。该研究发表于《Radiation Measurements》。

### 二、主要技术方法

研究人员采用以下关键技术方法：基于椭圆积分闭式表达式的确定性解析响应核构建（源自Cordella等人的等半径共轴PRT响应模型，包含Bragg-Kleeman能损关系）；三次B样条基函数展开中子谱；带二阶导数Tikhonov罚项和L曲线正则化参数自动选择的有约束最小二乘反演；有界变量最小二乘（Bounded-Variable Least Squares, BVLS）精确施加非负性约束；泊松加权白化处理计数统计特性；bootstrap重采样进行不确定度量化。模型外验证采用独立Geant4输运模拟生成反冲质子数据，Geant4版本为11.1.1（2.45 MeV）与11.3.2（14.1 MeV），物理列表均采用FTFP_BERT_HP。

### 三、研究结果

**合成基准测试：条件数与噪声传播。** 研究人员在匹配算子配置下测试了三峰参考谱（峰位4.0、8.0和13.5 MeV），L曲线自动选取正则化参数λ=1.8×10^-7。三个峰均被正确恢复，最大峰位误差仅0.13 MeV（约3%），约化卡方χ²_ν=0.99表明拟合在统计上自洽。响应矩阵条件数约1.45×10⁴，属于中等程度不适定。在覆盖分离良好、非对称、紧密间隔、混合宽度、五峰及重度重叠等六类合成谱的测试中，除重度重叠案例外，峰位误差均不超过0.13 MeV，χ²_ν接近1。

**基于Geant4数据的模型外验证。** 该测试打破了匹配算子配置——数据由独立Geant4输运生成而响应矩阵由解析核构建，因而可检验核本身的物理正确性。2.45 MeV能量点覆盖五种辐射体厚度（5-100 μm），解析核恢复的中子线位在2.35-2.40 MeV之间（偏差≤4%），且无虚假结构；14.1 MeV能量点恢复峰位13.9 MeV（偏差-1.4%）。联合双线测试成功分辨两条谱线，但D-T/D-D强度比约为真实值的2.4倍，定量揭示了核的各向同性质心系发射假设对14 MeV附近谱线间绝对强度的限制。稳健性测试表明，λ变化30倍、本底扣除变化±40%均不影响峰位恢复。

**不确定性量化与期望最大化比较。** 500次泊松bootstrap重采样给出逐点95%置信区间，真实谱在各能量点均落于置信带内。置信带在峰位处较宽（5-10%），体现了正则化反演放大不确定性的特征。与EM基线相比，B样条解的χ²_ν=0.99统计一致；而200次迭代的EM虽然原始残差更低（δ=0.0134对0.0486），但χ²_ν=0.32表明其已过度拟合计数噪声。B样条解约2秒完成，优于EM约8秒且无需人工选择迭代次数。

### 四、讨论与结论

讨论部分将方法局限性分为两类。第一类继承自解析核本身假设：要求共轴等半径几何、薄靶单散射近似、非相对论1-20 MeV能区，忽略γ本底、多次散射和碳反冲；且模型外验证目前仅覆盖单一准直几何、两种线能量和五种厚度。第二类源于正则化反演内在限制：能量分辨率受bin宽度约束，最小可分辨峰间距约1-1.5 MeV（峰宽0.5-0.6 MeV时）。计算性能方面，完整分析流程约5秒完成，内存占用小于1 MB，相比蒙特卡罗响应构建加速3-4个数量级。

研究结论为：该工作将Cordella等人的闭式椭圆积分PRT响应与Trunec等人的B样条解谱框架相结合，证明了响应矩阵可解析构建而无需蒙特卡罗输入。其新贡献包括：解析、无噪声、可复现的响应矩阵构建；BVLS方法精确强制非负性（完成了前人遗留问题）；L曲线自动正则化替代交叉验证；泊松加权反演公式；以及首次独立模型外验证——在2.45 MeV（五种辐射体厚度）和14.1 MeV下，解析核恢复中子线位的精度在4%以内。该方法提供了1-20 MeV快中子能谱测量的快速确定性框架，在聚变诊断、辐射防护和核物理领域具有重要应用前景。