MCNPX simulation of proton dose distribution in homogeneous and CT phantoms

A dose simulation system was constructed based on the MCNPX Monte Carlo package to simulate proton dose distribution in homogeneous and CT phantoms. Conversion from Hounsfield unit of a patient CT image set to material information necessary for Monte Carlo simulation is based on Schneider's approach. In order to validate this simulation system, inter-comparison of depth dose distributions among those obtained from the MCNPX, GEANT4 and FLUKA codes for a 160MeV monoenergetic proton beam incident normally on the surface of a homogeneous water phantom was performed. For dose validation within the CT phantom, direct comparison with measurement is infeasible. Instead, this study took the approach to indirectly compare the 50% ranges (R_50%) along the central axis by our system to the NIST CSDA ranges for beams with 160 and 115MeV energies. Comparison result within the homogeneous phantom shows good agreement. Differences of simulated R_50% among the three codes are less than 1mm. For results within the CT phantom, the MCNPX simulated water equivalent R_eq,50% are compatible with the CSDA water equivalent ranges from the NIST database with differences of 0.7 and 4.1mm for 160 and 115MeV beams, respectively.