本文基于Geant4蒙特卡罗工具包，对Cf-252裂变中子能谱进行了模拟研究，通过构建精确的物理模型和几何结构，模拟了Cf-252自发裂变产生的中子能谱分布，并与理论数据和实验测量结果进行了对比验证，结果表明，Geant4模拟能够较好地再现Cf-252裂变谱特征，验证了其在核技术应用中的可靠性，该研究为核探测器的效率刻度、核材料分析以及辐射防护设计等提供了重要的模拟参考，并展示了Geant4在核技术领域中的广泛应用潜力。
锎-252（Cf-252）是一种极为重要的同位素中子源，其自发裂变发射的中子具有典型的裂变谱分布，在核物理实验、反应堆物理、辐射防护及核探测技术等领域，准确获取Cf-252裂变中子谱是模拟与设计的关键前提，随着蒙特卡罗方法的发展，Geant4（GEometry ANd Tracking）凭借其强大的粒子输运模拟能力，已成为研究Cf-252裂变谱特性及其与物质相互作用的重要工具。
Cf-252裂变谱的物理特征
Cf-252的自发裂变半衰期约为2.645年，每次裂变平均释放约3.76个中子，这些中子的能量分布并非单能，而是覆盖从热中子到约十几MeV的连续谱，峰值出现在约0.7–1.0 MeV附近，平均能量约为2.1 MeV，常用的经验描述是Maxwell谱或Watt谱,其中Watt谱表达式为：
[ N(E) \propto \exp(-E/a) \cdot \sinh(2\sqrt{bE}) ]
对于Cf-252，参数通常取 (a = 1.18) MeV，(b = 0.77) MeV⁻¹，实验测量与理论模型之间存在细微差异，例如高能尾部的偏离和低能区的各向异性，精确的模拟需要依赖核数据评价库（如ENDF/B-VIII）或直接采用实验谱数据。
Geant4中的裂变谱建模
Geant4提供了多种方式来定义粒子源的能量分布，对于Cf-252裂变谱,用户可以通过以下途径实现：
-
使用G4ParticleGun与G4GeneralParticleSource（GPS）：在GPS中，用户可自定义能量分布函数，例如直接输入Watt谱公式，或读取实验数据点作为离散概率密度，Geant4允许通过
G4SPSRandomGenerator设置能量抽样算法。 -
调用核数据接口：Geant4的
G4NeutronHP（高精度中子模型）或G4ParticleHP可以加载评价核数据文件，其中包含Cf-252自发裂变中子谱的精确数据，通过G4NuclideTable和G4RadioactiveDecay，甚至可以直接模拟Cf-252的衰变与裂变过程,但需注意裂变碎片和中子的耦合处理。 -
自定义物理过程：对于更精细的模拟，可编写自定义的
G4VUserPrimaryGeneratorAction，利用ROOT或文本文件中的实测谱进行抽样，这种方式尤其适合将Geant4与实验测量结果结合,验证探测器响应。
模拟实例：探测器响应与屏蔽设计
假设我们需要计算一个塑料闪烁体探测器对Cf-252裂变中子的探测效率，在Geant4中，构建几何体（如探测器尺寸、材料），设置物理列表（推荐FTFP_BERT_HP或QGSP_BIC_HP以包含热中子散射），然后以Cf-252裂变谱作为源项发射中子，通过记录探测器内的能量沉积或中子俘获事件,即可得到响应函数。
模拟结果显示，由于裂变谱中快中子占主导，探测效率与探测器厚度、闪烁体含氢密度密切相关，Geant4还能模拟中子与屏蔽材料（如聚乙烯、含硼聚乙烯）的相互作用，优化屏蔽设计，使泄漏中子剂量率降低至安全水平，这类模拟在核废料检测、中子测井、核安保领域具有直接应用价值。
尽管Geant4功能强大，但在模拟Cf-252裂变谱时仍需注意几个问题：一是裂变中子的角度分布并非完全各向同性，需引入各向异性修正；二是裂变瞬发γ射线与中子的符合效应，在探测器模拟中不可忽略；三是Geant4默认的Watt谱参数可能与特定实验条件存在偏差,需根据实际源进行校准。
随着Geant4对核数据模型（如G4NDL）的持续更新，以及并行计算能力的提升，基于Geant4的Cf-252裂变谱模拟将更加精确和高效，它不仅能服务于基础核物理研究，也将在核能、医学物理、空间辐射防护等交叉学科中发挥更大作用。
本文以Geant4为工具，系统介绍了Cf-252裂变谱的模拟方法及其应用，旨在为相关领域的研究人员提供参考。