TOPAS Monte Carlo simulation for a scanning proton therapy system in SPHIC

Purpose: This study aims to use TOPAS Monte Carlo (MC) code to model the beamline of a scanned proton therapy system in Shanghai Proton and Heavy Ion Center (SPHIC). This tool can be useful for treatment plan quality assurance. Methods and materials: In-air lateral profiles of single proton spots at five positions along the Z-axis were measured, which can then be used to calculate the angular variance, spatial angular covariance, and spatial variance at the beam exit position. The mean energy and energy spread of the primary incident beam for 15 energies were tuned by comparing with measured relative integrated depth dose at the central axis. Three spread-out Bragg peak plans in homogeneous water cubes and one plan with a half anthropomorphic head phantom were created to verify the accuracy of the model. In addition, the comparison was also performed to RayStation TPS with an MC dose calculation engine for the plans. Results: The average difference of spot sizes was 1.3% for all five energies. For all 15 energies, the average range (D90) difference was 0.04 mm, the mean distal-fall-off (DDF) difference was 0.04 mm, and the mean full width at half maximum (FWHM) difference was 0.10 mm. For the three homogeneous cube plans, a maximum difference of 0.2 mm for the D90 was observed, and the absolute differences in the lateral penumbras for each plan all agreed within 1 mm. For the plan on the anthropomorphic half head phantom, the mean dose deviation is less than 1.2% compared to measurements. Conclusion: We successfully developed an MC model for the proton scanning system used in SPHIC, which can be helpful for dose verification of the Treatment Planning System (TPS) and quality assurance in the future.