Respiratory motion is traditionally assessed using tumor motion magnitude. In proton therapy, respiratory motion causes density variations along the beam path that result in uncertainties of proton range. This work has investigated the use of water‐equivalent thickness (WET) to quantitatively assess the effects of respiratory motion on calculated dose in passively scattered proton therapy (PSPT). A cohort of 29 locally advanced non‐small cell lung cancer patients treated with 87 PSPT treatment fields were selected for analysis. The variation in WET (ΔWET) along each field was calculated between exhale and inhale phases of the simulation four‐dimensional computed tomography. The change in calculated dose (ΔDose) between full‐inhale and full‐exhale phase was quantified for each field using dose differences, 3D gamma analysis, and differential area under the curve (ΔAUC) analysis. Pearson correlation coefficients were calculated between ΔDose and ΔWET. Three PSPT plans were redesigned using field angles to minimize variations in ΔWET and compared to the original plans. The median ΔWET over 87 treatment fields ranged from 1‐9 mm, while the ΔWET 95th percentile value ranged up to 42 mm. The ΔWET was significantly correlated to the ΔDose for all metrics analyzed. The patient plans that were redesigned using ΔWET analysis to select field angles were more robust to the effects of respiratory motion, as ΔAUC values were reduced by more than 60% in all three cases. The tumor motion magnitude alone does not capture the potential dosimetric error due to respiratory motion because the proton range is sensitive to the motion of all patient anatomy. The use of ΔWET has been demonstrated to identify situations where respiratory motion can impact the calculated dose. Angular analysis of ΔWET may be capable of designing radiotherapy plans that are more robust to the effects of respiratory motion.PACS number(s): 87.55.‐x
Evaluation of the systematic error in using 3D dose calculation in scanning beam proton therapy for lung cancer
The objective of this study was to evaluate and understand the systematic error between the planned three‐dimensional (3D) dose and the delivered dose to patient in scanning beam proton therapy for lung tumors. Single‐field and multifield optimized scanning beam proton therapy plans were generated for ten patients with stage II‐III lung cancer with a mix of tumor motion and size. 3D doses in CT datasets for different respiratory phases and the time‐weighted average CT, as well as the four‐dimensional (4D) doses were computed for both plans. The 3D and 4D dose differences for the targets and different organs at risk were compared using dose‐volume histogram (DVH) and voxel‐based techniques, and correlated with the extent of tumor motion. The gross tumor volume (GTV) dose was maintained in all 3D and 4D doses, using the internal GTV override technique. The DVH and voxel‐based techniques are highly correlated. The mean dose error and the standard deviation of dose error for all target volumes were both less than 1.5% for all but one patient. However, the point dose difference between the 3D and 4D doses was up to 6% for the GTV and greater than 10% for the clinical and planning target volumes. Changes in the 4D and 3D doses were not correlated with tumor motion. The planning technique (single‐field or multifield optimized) did not affect the observed systematic error. In conclusion, the dose error in 3D dose calculation varies from patient to patient and does not correlate with lung tumor motion. Therefore, patient‐specific evaluation of the 4D dose is important for scanning beam proton therapy for lung tumors.PACS number: 87.55.D