摘要：本文基于核内级联物理过程,采用Monte Carlo方法发展了一款质子、中子以及π介子的粒子输运程序.基本物理模型基于适当简化和核内级联Bertini模型,同时借鉴了INCL模型质心系下的角微分分布以克服Bertini模型之不足,即采用Monte Carlo方法模拟核子与核子、核子与π介子间的弹性散射、非弹性散射等过程,粒子相互作用时,核子密度随半径变化且作用截面参考Bertini模型22类实验截面数据,出射粒子散射角在质心系下的抽样遵从INCL模型所确定的微分分布.可模拟45—3500 MeV的中子、质子或2500 MeV以下π介子引起的核内级联过程.入射粒子能量在60—378 MeV范围内反应截面理论计算值与已有实验数据、以及在65—3000 MeV较宽能区范围内反应截面、出射粒子增殖比、微分截面和剩余核等计算结果与MCNPX,GEANT 4和PHITS模拟结果符合较好.

关键词：核子 核内级联 monte carlo方法 

单位：中国科学院高能物理研究所; 北京100049