Encrypted login | home

Program Information

Monte Carlo Calculation of Source Terms and Attenuation Lengths for Neutrons Produced by 50-200 MeV Protons On Brass

no image available
J Ramos-Mendez

J Ramos-Mendez1*, H Paganetti2 , B Faddegon1 , (1) University of California San Francisco, San Francisco, CA, (2) Massachusetts General Hospital, Boston, MA


SU-E-T-554 (Sunday, July 12, 2015) 3:00 PM - 6:00 PM Room: Exhibit Hall

Purpose: We used TOPAS (TOPAS wraps and extends Geant4 for medical physicists) to compare Geant4 physics models with published data for neutron shielding calculations. Subsequently, we calculated the source terms and attenuation lengths (shielding data) of the total ambient dose equivalent (TADE) in concrete for neutrons produced by protons in brass.

Methods: Stage1: The Bertini and Binary nuclear models available in Geant4 were compared with published attenuation at depth of the TADE in concrete and iron. Stage2: Shielding data of the TADE in concrete was calculated for 50-200 MeV proton beams on brass. Stage3: Shielding data from Stage2 was extrapolated for 235 MeV proton beams. This data was used in a point-line-source analytical model to calculate the ambient dose per unit therapeutic dose at two locations inside one treatment room at the Francis H Burr Proton Therapy Center. Finally, we compared these results with experimental data and full TOPAS simulations.

Results: At larger angles (~130o) the TADE in concrete calculated with the Bertini model was about 9 times larger than that calculated with the Binary model. The attenuation length in concrete calculated with the Binary model agreed with published data within 7%±0.4% (statistical uncertainty) for the deepest regions and 5%±0.1% for shallower regions. For iron the agreement was within 3%±0.1%. The ambient dose per therapeutic dose calculated with the Binary model, relative to the experimental data, was a ratio of 0.93±0.16 and 1.23±0.24 for two locations. The analytical model overestimated the dose by four orders of magnitude. These differences are attributed to the complexity of the geometry.

Conclusion: The Binary and Bertini models gave comparable results, with the Binary model giving the best agreement with published data at large angle. Shielding data we calculated using the Binary model is useful for fast shielding calculations with other analytical models.

Funding Support, Disclosures, and Conflict of Interest: his work was supported by National Cancer Institute Grant R01CA140735

Contact Email: