[论文解读] The Geant4 Hadronic Verification Suite for the Cascade Energy Range

V. Ivanchenko, G. Folger+3 more|ArXiv.org|Jun 2, 2003

Nuclear reactor physics and engineering参考文献 8被引用 7

一句话总结

本论文提出了一套基于Geant4的验证套件，用于级联能区（高达几GeV）的强子模型验证，利用实验截面数据验证质子引发反应中中子和π介子的产生。该套件可对模型预测（Bertini模型与二体碰撞模型）与测量结果进行系统比较，结果显示对微分截面和π⁺/π⁻产生比值的再现性良好，但在较高能量下归一化值被低估了2–3倍。

ABSTRACT

A Geant4 hadronic process verification suite has been designed to test and optimize Geant4 hadronic models in the cascade energy range. It focuses on quantities relevant to the LHC radiation environment and spallation source targets. The general structure of the suite is presented, including the user interface, stages of verification, management of experimental data, event generation, and comparison of results to data. Verification results for the newly released Binary cascade and Bertini cascade models are presented.

研究动机与目标

• 开发Geant4强子模型在级联能区（阈值至几GeV）的系统性验证框架，对LHC辐射和散裂靶应用至关重要。
• 实现对Geant4强子物理模型（特别是新发布的Bertini和二体碰撞模型）的独立测试与优化。
• 提供一个模块化、可扩展的接口，用于比较模型生成的事件数据与高精度实验测量的微分和双微分截面数据。
• 确保对次级粒子能谱（中子、π介子）模拟的准确性，以满足加速器和探测器设计的需求。
• 将验证范围扩展至更高入射能量及更多粒子种类（如π介子、中子），超越初始的质子束实验。

提出的方法

• 该套件使用通用接口连接到Geant4的Level 2强子物理框架，实现对末态生成器的独立测试。
• 通过定义入射束流能量和靶核，生成单个事件，利用分档相空间区间追踪次级粒子的运动学参数（能量、动量、角度、重子数）。
• 截面通过公式 Δσᵢ = σₜₒₜ × Nᵢ / N 计算，双微分截面估算为 d²σ/dEdΩ = Δσᵢ / (ΔE × ΔΩ)。
• 实验数据来自EXFOR数据库，经重格式化以供PAW分析，对低能谱进行重采样，当仅有双微分数据时采用数值积分。
• 用户界面由宏文件和kumac脚本驱动，支持线性与对数分档，输出以HBOOK和ASCII格式存储以供分析。
• 在多种测试案例中执行验证，包括不同束流能量、靶核和相互作用模型，结果与已发表测量数据进行对比。

实验结果

研究问题

• RQ1 Geant4中的二体碰撞模型和Bertini碰撞模型在100–800 MeV质子轰击下，对中子双微分截面的再现精度如何？
• RQ2 在597–730 MeV质子轰击下，这些模型对π⁺/π⁻产生比值的再现程度如何？
• RQ3 这些模型在不同靶材料（Be、C、Al、Fe、Ni、Pb）中对次级粒子能谱的能量和角度依赖性的预测能力如何？
• RQ4 随着入射能量增加，模型在归一化和形状上与实验数据的一致性表现如何？
• RQ5 该验证套件能否作为可扩展、模块化的框架，用于未来Geant4强子模型的验证？

主要发现

• 二体碰撞模型在Be、Al、Fe和Pb上对质子散射中中子产生微分截面的再现具有良好的形状和峰值位置一致性，但在较高能量下归一化值被低估了2–3倍。
• Bertini碰撞模型在中子产额和角分布方面与实验数据表现出相似的一致性，尤其在256 MeV质子轰击铝靶时，对峰值结构的再现效果良好。
• 二体碰撞模型在多个靶核（C、Al、Ni、Cu）上对π⁺/π⁻产生比值的再现效果良好，正确捕捉了前向角π⁺占主导的现象。
• π介子产生截面的能量依赖性被定性地再现，但整体归一化值被低估了2–3倍，且在较高角度处偏差随能量增加而加剧。
• 验证套件成功验证了新的Geant4强子模型，并已扩展至更高入射能量（最高达15 GeV）和更多粒子种类（π介子、中子）。
• 该套件作为模块化、可重用的验证框架效果显著，目前已用于Geant4强子物理的进一步开发。