Benchmark measurements and simulations of dose perturbations due to metallic spheres in proton beams

•We compared measurements and Monte Carlo predictions of dose perturbations caused by the metal objects in proton beams.•Different Monte Carlo codes were used, including MCNPX, GEANT4 and Fast Dose Calculator.•Good agreement was found between measurements and Monte Carlo simulations.•The modification of multiple Coulomb scattering model in MCNPX code yielded improved accuracy.•Our results confirmed that Monte Carlo codes are well suited for predicting multiple Coulomb scattering in proton therapy.

Monte Carlo simulations are increasingly used for dose calculations in proton therapy due to its inherent accuracy. However, dosimetric deviations have been found using Monte Carlo code when high density materials are present in the proton beamline. The purpose of this work was to quantify the magnitude of dose perturbation caused by metal objects. We did this by comparing measurements and Monte Carlo predictions of dose perturbations caused by the presence of small metal spheres in several clinical proton therapy beams as functions of proton beam range and drift space. Monte Carlo codes MCNPX, GEANT4 and Fast Dose Calculator (FDC) were used. Generally good agreement was found between measurements and Monte Carlo predictions, with the average difference within 5% and maximum difference within 17%. The modification of multiple Coulomb scattering model in MCNPX code yielded improvement in accuracy and provided the best overall agreement with measurements. Our results confirmed that Monte Carlo codes are well suited for predicting multiple Coulomb scattering in proton therapy beams when short drift spaces are involved.