Clinical validation of a GPU-based Monte Carlo dose engine of a commercial treatment planning system for pencil beam scanning proton therapy.

PURPOSE: To perform the validation of the GPU-based (Graphical Processing Unit based) proton Monte Carlo (MC) dose engine implemented in a commercial TPS (RayStation 10B) and to report final dose calculation times for clinical cases. MATERIALS AND METHODS: 440 patients treated at the Proton Therapy Center of Trento, Italy, between 2018 and 2019 were selected for this study. 636 approved plans with 3361 beams computed with the clinically implemented CPU-MC dose engine (version 4.2 and 4.5), were used for the validation of the new algorithm. For each beam, the dose was recalculated using the new GPU-MC dose engine with the initial CPU computation settings and compared to the original CPU-MC dose. Beam dose difference distributions were studied to ensure that the two dose distributions were equal within the expected fluctuations of the MC statistical uncertainty (s) of each computation. Plan dose distributions were compared with respect to the dosimetric indices D98, D50 and D1 of all ROIs defined as targets. A complete assessment of the computation time as a function of s and dose grid voxel size was done. RESULTS: The median over all mean beam dose differences between CPU- and GPU-MC was -0.01% and the median of the corresponding standard deviations was close to (√2s) both for simulations with an s of 0.5% and 1.0% per beam. This shows that the two dose distributions can be considered equal. All the DVH indices showed an average difference below 0.04%. About half of the plans were computed with 1.0% statistical uncertainty on a 2 mm dose calculation grid, for which the median computation time was 5.2 s. The median computational speed for all plans in the study was 8.4 million protons/second. CONCLUSION: A validation of a clinical MC algorithm running on GPU was performed on a large pool of patients treated with pencil beam scanning proton therapy. We demonstrated that the differences with the previous CPU-based MC were only due to the intrinsic statistical fluctuations of the MC method, which translated to insignificant differences on plan dose level. The significant increase in dose calculation speed is expected to facilitate new clinical workflows.

Clinical implementation of pencil beam scanning proton therapy for liver cancer with forced deep expiration breath hold.

PURPOSE: To present our technique for liver cancer treatments with proton therapy in pencil beam scanning mode and to evaluate the impact of uncertainties on plan quality. MATERIALS AND METHODS: Seventeen patients affected by liver cancer were included in this study. Patients were imaged and treated in forced breath-hold using the Active Breathing Coordinator system and monitored with an optical tracking system. Three simulation CTs were acquired to estimate the anatomical variability between breath-holds and generate an internal target volume (ITV). The treatment plans were optimized with a Single Field Optimization technique aimed at minimizing the use of range shifter. Plan robustness was tested simulating systematic range and setup uncertainties, as well as the interplay effect between breath-holds. The appropriateness of margin was further verified based on the actual positioning data acquired during treatment. RESULTS: The dose distributions of the nominal plans achieved a satisfactory target coverage in 11 out of 17 patients, while in the remaining 6 D95 to the PTV was affected by the constraint on mean liver dose. The constraints for all other organs at risk were always within tolerances. The interplay effect had a limited impact on the dose distributions: the worst case scenario showed a D95 reduction in the ITV < 3.9 GyRBE and no OAR with D1 > 105% of the prescription dose. The robustness analysis showed that for 13 out of 17 patients the ITV coverage in terms of D95 was better than D95 of the PTV in the nominal plan. For the remaining 4 patients, the maximum difference between ITV D95 and PTV D95 was ≤0.7% even for the largest simulated setup error and it was deemed clinically acceptable. Hot spots in the OARs were always lower than 105% of the prescription dose. Positioning images confirmed that the breath hold technique and the PTV margin were adequate to compensate for inter- and intra-breath-hold variations in liver position. CONCLUSION: We designed and clinically applied a technique for the treatment of liver cancer with proton pencil beam scanning in forced deep expiration breath-hold. The initial data on plan robustness and patient positioning suggest that the choices in terms of planning technique and treatment margins are able to reach the desired balance between target coverage and organ at risk sparing.