TOPAS Monte Carlo simulation for double scattering proton therapy and dosimetric evaluation.

Purpose: To construct and commission a double scattering (DS) proton beam model in TOPAS Monte Carlo (MC) code. Dose comparisons of MC calculations to the measured and treatment planning system (TPS) calculated dose were performed.

Methods: The TOPAS nozzle model was based on the manufacturer blueprints. Nozzle set-up and beam current modulations were calculated using room-specific calibration data. This model was implemented to reproduce pristine peaks, spread-out Bragg peaks (SOBP) and lateral profiles. A stair-shaped target plan in water phantom was calculated and compared to measured data to verify range compensator (RC) modeling.

Results: TOPAS calculated pristine peaks agreed well with measurements, with accuracies of 0.03 cm for range R and 0.05 cm for distal dose fall-off (DDF). The calculated SOBP range, modulation width and DDF differences between MC calculations and measurements were within 0.05 cm, 0.5 cm and 0.03 cm respectively. MC calculated lateral penumbra agreed well with measured data, with difference less than 0.05 cm. For RC calculation, TPS underestimated the additional depth dose tail due to the nuclear halo effect. Lateral doses by TPS were 10% lower than measurement outside the target, while maximum difference of MC calculation was within 2%. At deeper depths inside the target volume, TPS overestimated doses by up to 25% while TOPAS predicted the dose to within 5% of measurements.

Conclusion: We have successfully developed and commissioned a MC based DS nozzle model. The performance of dose accuracy by TOPAS was superior to TPS, especially for highly inhomogeneous compensator.