论文标题：Open-Source Optimization of Hybrid Monte Carlo Methods for Fast Response Modeling of NaI (Tl) and HPGe Gamma Detectors

论文链接：https://www.mdpi.com/2673-4362/5/3/19

期刊名： Journal of Nuclear Engineering

期刊主页：https://www.mdpi.com/journal/jne

在核工程领域，伽马探测器的响应建模长期面临高成本、出口管制及计算资源消耗大等挑战。传统的MCNP和GADRAS等行业标准软件虽能提供高保真结果，但申请周期长达数月、对非专业用户不友好，且传统蒙特卡罗方法需数小时甚至数天完成复杂模拟。针对上述瓶颈，美国普渡大学核工程学院的Matthew Niichel与Stylianos Chatzidakis团队在Journal of Nuclear Engineering发表最新成果，开发了一款基于MATLAB的开放源应用程序，通过混合蒙特卡罗方法实现了NaI(Tl)和HPGe伽马探测器的快速响应建模，并配备图形用户界面，大幅降低了使用门槛与计算成本。

研究过程与结果

研究团队采用混合蒙特卡罗方法，将解析关系与传统蒙特卡罗随机采样相结合。该方法在初始化阶段预分配散射角分布（基于Klein Nishina公式），随后通过条件判断与随机数比较，依次模拟光電效应、康普顿散射和电子对效应三种伽马-物质相互作用。下图展示了算法中决定伽马事件类型的判断逻辑：

伽马-物质相互作用事件算法比较步骤框图

研究基于实测数据对模型进行了全面验证。实验采用约10微居里的Co 60、Cs 137和Na 22源，在NaI和HPGe探测器上分别采集100 s、300 s和600 s的能谱，并与模拟结果对比。结果表明：模拟谱与实验谱在康普顿连续区和光电峰位处高度吻合，相对误差在关键特征区域约45%–80%，考虑到随机噪声属可接受范围。此外，模型允许用户自定义探测器分辨率（0–100%）、同位素年龄（天/年）、计数时间（1–3000 s）及最多三种同位素组合，并支持以CSV格式导出能谱数据。

计算性能方面，在16核Dell i7处理器上模拟单同位素仅需0.13–0.48秒，在4核HP i7上约0.8–2.1秒，相比传统蒙特卡罗方法快2–3个数量级。与GADRAS对比，两者能谱架构高度相似，但普渡代码完全开放源且无需申请。

研究总结

该研究明确指出：开放源混合蒙特卡罗方法可在保证与实验及行业标准代码可比的响应精度前提下，将计算时间缩短至秒级甚至亚秒级，并彻底规避出口管制与高昂成本问题。所开发的MATLAB应用程序“GammaGenius”已发布于MATLAB File Exchange，可广泛应用于核工程教学、快速能谱预测及探测器响应数据集生成。

研究团队未来计划：(1)扩展同位素库至与GADRAS相当；(2)增加能量校准或偏压模拟功能；(3)引入中子探测器（如BF?）响应模块。虽无意替代MCNP等专业软件，但该工作为核工程教育与小规模研究提供了便捷、低成本、可复现的探测器建模新工具。

期刊介绍

主编：Dan Gabriel Cacuci, University of South Carolina, USA

期刊创刊于2020年, 最新Impact Factor 1.2，CiteScore 2.6，是一个国际开放获取期刊，发表同行评审的论文。期刊涵盖了与核和辐射过程的科学和应用相关的原创研究、想法和进展。发文类型包括original research papers, reviews, communications, brief reports, opinions, technical notes, editorials等。

截至目前，Journal of Nuclear Engineering期刊已被Scopus, ESCI (Web of Science), EBSCO等数据库收录。