SU-E-T-474: Monte Carlo Phase Space Production to Model Magnetically Scanned Proton Beams for IMPT

PURPOSE

Accurate dose predictions in proton beam therapy using magnetically scanned beams are highly dependent on the accurate modeling of the lateral dose profiles. This study was performed to provide proton phase spaces for Monte Carlo simulations, used to accurately simulate doses at distances up to 12 cm from the central axis of the beam.

METHODS

Measured lateral dose profiles at various depths in water were compared to Monte Carlo simulations of doses for 90 discreet initial proton energies. Phase spaces were produced using a one dimensional energy distribution, and a combination of several two dimensional spatial and directional distributions. Simulations were performed iteratively using variations in the initial phase space distributions to achieve acceptable agreement between measured and simulated lateral dose profiles, i.e. differences in FWHM < 0.5 mm and dose differences less that 0.1% at distances up to 12.5 cm.

RESULTS

90 phase spaces of proton sources for different initial beam energies were created for use in Monte Carlo simulations of scanned proton beam therapy patient plans. At a depth of 2 cm in water, the simulated and measured FWHM of the lateral dose profiles differed in in-plane direction by an average of 0.05 mm, in cross-plane direction by 0.13 mm. All simulated profiles were within 0.1% of the measured doses at distances between 2cm and 12.5 cm from the central beam axis.

CONCLUSIONS

A library of 90 phase space files has been created to accurately simulate magnetically scanned proton beams for IMPT, providing accurate dose distributions up to 12 cm distance from the central beam axis. This project is supported in part by P01CA021239 from the National Cancer Institute. The content is solely the responsibility of the authors and does not necessarily represent the official views of the National Cancer Institute or the National Institutes of Health.

SU-E-T-624: Comparison of PTV+PRV-Based Optimization and Robust Optimization in Intensity-Modulated Proton Therapy

PURPOSE

Robust optimization leads to IMPT plans that are more robust than and superior in optimality compared to PTV-based optimized plans. Robust optimization incorporates setup and range uncertainties, which implicitly adds margins to targets and organs-at-risk (OARs); whereas PTV-based optimization only considers setup uncertainties and adds margins only to targets in practice. The purpose of this work is to determine if the superiority of robustly optimized plans is due to not assigning margins to OARs during PTV-based optimization.

METHODS

Plan robustness and optimality of the PTV plus Planning organs-at-Risk Volume (PRVs)-based plans and robustly optimized plans were compared for 5 head and neck cancer cases and one rhabdomyosarcoma case. The PRVs were generated by expansion from OARs by 3 mm. 9 different dose distributions were computed - one each for ± setup uncertainties along three spatial directions, for ± range uncertainty, and the nominal dose distribution. The worst-case dose distribution was obtained by assigning the lowest dose among the 9 doses to each voxel in the target and the highest dose to each voxel outside the target. The DVHs from the worst-case dose were used to assess the plan optimality and robustness. D1cc doses for spinal cord and brainstem, mean doses Dmean for oral cavity and parotids, and D1% doses for other organs were also used to assess plan optimality. D5% and D95% doses are used to assess target dose coverage and homogeneity.

RESULTS

For H&N cases, PTV+PRV-based optimization was inferior to robust optimization. However, PTV+PRV-based optimization yielded plans that spared OARs better than PTV-based optimization, although the target dose robustness and homogeneity were comparable to the PTV-based optimization. The same conclusions are also valid in the rhabdomyosarcoma case.

CONCLUSIONS

We find that the PTV+PRV method can partly improve plan optimality, but it is still inferior to robust optimization method. This research is supported by National Cancer Institute (NCI) grant P01CA021239, the University Cancer Foundation via the Institutional Research Grant program at the University of Texas MD Anderson Cancer Center, and MD Anderson’s cancer center support grant CA016672.

SU-E-T-625: Robustness Evaluation and Robust Optimization of IMPT Plans Based on Per-Voxel Standard Deviation of Dose Distributions

PURPOSE

Proton dose distributions, IMPT in particular, are highly sensitive to setup and range uncertainties. We report a novel method, based on per-voxel standard deviation (SD) of dose distributions, to evaluate the robustness of proton plans and to robustly optimize IMPT plans to render them less sensitive to uncertainties.

METHODS

For each optimization iteration, nine dose distributions are computed - the nominal one, and one each for ± setup uncertainties along x, y and z axes and for ± range uncertainty. SD of dose in each voxel is used to create SD-volume histogram (SVH) for each structure. SVH may be considered a quantitative representation of the robustness of the dose distribution. For optimization, the desired robustness may be specified in terms of an SD-volume (SV) constraint on the CTV and incorporated as a term in the objective function. Results of optimization with and without this constraint were compared in terms of plan optimality and robustness using the so called'worst case' dose distributions; which are obtained by assigning the lowest among the nine doses to each voxel in the clinical target volume (CTV) and the highest to normal tissue voxels outside the CTV. The SVH curve and the area under it for each structure were used as quantitative measures of robustness. Penalty parameter of SV constraint may be varied to control the tradeoff between robustness and plan optimality. We applied these methods to one case each of H&N and lung.

RESULTS

In both cases, we found that imposing SV constraint improved plan robustness but at the cost of normal tissue sparing.

CONCLUSIONS

SVH-based optimization and evaluation is an effective tool for robustness evaluation and robust optimization of IMPT plans. Studies need to be conducted to test the methods for larger cohorts of patients and for other sites. This research is supported by National Cancer Institute (NCI) grant P01CA021239, the University Cancer Foundation via the Institutional Research Grant program at the University of Texas MD Anderson Cancer Center, and MD Anderson’s cancer center support grant CA016672.

WE-G-BRCD-02: Speed Up Standard Optimization of IMPT Treatment Planning Through the Application of GPUs

PURPOSE

While intensity-modulated proton therapy (IMPT) has great potential to deliver highly conformal tumoricidal dose to targets whilst minimizing dose to nearby organs-at-risk, IMPT optimization is very time consuming and memory extensive due to finer dose grids and a large number of energy layers used compared to intensity-modulated radiation therapy (IMRT). In this presentation, for the first time, a new approach is introduced to speed up the IMPT treatment planning through application of parallel computing with Graphic Processor Units (GPUs).

METHODS

Parallel computation with GPUs, which are affordable and can be plugged in a workstation easily, is potentially a good way to improve the computation efficiency. In our approach, we used the standard quadratic objective function to optimize the intensity map of beamlets. The objective function and gradient equations, which are the most time consuming parts of the optimization, were calculated with GPUs. We compared the computation time of optimization done by an Intel ® Core™ i7 CPU and that by the same CPU accelerated by GPUs (TESLA C1060). The influence matrix was pre- calculated before optimization with an in-house proton pencil beam dose calculation engine. Two clinical cases were studied: one base-of-skull (BOS) case (clivus chordoma) and one prostate case (adenocarcinoma). The dose volume histogram (DVH) data for the tumor and critical organs were derived for comparison of optimization results generated by CPU and GPUs.

RESULTS

For the BOS case, application of GPUs for the optimization and overall gained 54 and 36.5 times speedup. For the prostate case, application of GPUs for the optimization and overall gained 69 and 28.5 times speedup.

CONCLUSIONS

The application of GPUs for the parallel computing of IMPT treatment plan optimization can dramatically improve the computation efficiency. The optimization time can be reduced from typically half to one hour to only several minutes.

SU-E-T-394: Comparison of Planned Dose Distribution Vs. Delivered Dose Distribution for Both IMRT and Proton Therapy Using Weekly Repeat 4DCT Data Sets

PURPOSE

To evaluate treatment deviations and the impact of treatment modality in the presence of breathing motion and anatomical changes during the course of lung cancer radiotherapy.

METHODS

Two non-small cell lung cancer patients were enrolled in a randomized clinical trial to compare IMRT and proton therapy. To rigorously evaluate the impact of motion and anatomical changes, we used a '5D' dose accumulation approach to sum dose distributions from phase-to-phase and week-to-week to the reference (end-expiration) phase of the original planning 4DCT data set. Six to eight weekly 4DCT data sets that consisted of 10 breathing phases were acquired during the treatment course. The original plan was re-calculated for each phase and deformably mapped to the reference phase to compare the 'delivered' dose distribution with the planned dose distribution of both the IMRT plan and the proton plan for each patient. DVHs derived from delivered dose distribution were compared to that from the planned dose distribution.

RESULTS

The delivered dose showed 3% and 2% increase in the dose to the CTV for IMRT and proton plan respectively. Target coverage remained acceptable despite tumor shrinkage from 29% to 49%. The doses to normal structures, such as lung and heart, increased more in the proton plan than in the IMRT plan. The V20 of the total lung volume increased by 4% and 6% from the delivered dose compared to the planned dose for IMRT and proton plan, respectively.

CONCLUSIONS

The results showed sufficient target coverage was maintained for both modalities. Increases in lung dose were observed in both modalities, but more in the proton arm, perhaps due to weight loss and tumor shrinkage. Adaptive proton therapy strategy is recommended to minimize normal tissue doses. Supported in part by NCI P01 CA021239-29A1.

SU-E-T-396: Beam Angle Optimization for IMPT Comparing Gantries and Fixed Beam Lines

PURPOSE

To explore the potential of beam angle optimization (BAO) for IMPT and compare fixed beamlines with gantries.

METHODS

For three patients with challenging intracranial lesions, we generate reference IMPT treatment plans applying three manually selected beam orientations and treatment plans applying three optimized beam orientations considering five scenarios: (1) patients are in supine position and the treatment room features (1.a) a horizontal beamline, (1.b) a horizontal, 45°, and vertical beamline, (1.c) a gantry, (2) patients are in supine or seated position and the treatment room features (2.a) a horizontal beamline, or (2.b) a horizontal, 45°, and vertical beamline. We use a genetic algorithm that considers up to 1,400 non-coplanar candidate beams and evaluates 10,000 beam ensembles for one BAO. Beam orientations that may compromise the robustness of treatment plans are excluded before the optimization based on an objective measure of existing tissue heterogeneities.

RESULTS

The optimized beam ensembles exhibit certain similarities even though the sets of candidate beams differ significantly for the five scenarios. Compared to manually selected beam orientations, they provide improved OAR sparing and equivalent target coverage. Compared to one another, they yield comparable target conformity (deviations of the conformity number <1%), target homogeneity (standard deviations of the target dose <0.8 Gy), and sparing of OARs (deviations of average mean and maximum doses in OARs +/- 1 Gy). Using a gantry, however, the integral dose can be reduced by 5-15% compared to a horizontal beamline with patients in supine position. For the investigated cases comparable reductions can be achieved by also irradiating in seated position with a horizontal, 45°, and vertical beamline.

CONCLUSIONS

BAO has the potential to provide beneficial IMPT treatment plans. Compared to fixed beamlines, gantries yield only modest effects regarding OAR sparing but may enable a significant reduction of integral dose for individual patients.

Dosimetric evaluation of 120-leaf multileaf collimator in a Varian linear accelerator with 6-MV and 18-MV photon beams

In this study the dosimetric characteristics of 120-leaf multileaf collimators (MLCs) were evaluated for 6-MV and 18-MV photon beams. The dose rate, percentage depth dose, surface dose, dose in the build-up region, beam profile, flatness, symmetry, and penumbra width were measured using three field-defining methods: (i) 'Jaw only', (ii) 'MLC only', and (iii) 'MLC+Jaw'. Analysis of dose rate shows that the dose rate for 'MLC only' field was higher than that for 'Jaw only" and 'MLC+Jaw' fields in both the energies. The 'percentage of difference' of dose rates between 'MLC only' and 'MLC+Jaw' was (0.9% to 4.4%) and (1.14% to 7%) for 6 MV and 18 MV respectively. The surface dose and dose in the build-up region were more pronounced for 'MLC only' fields for both energies, and no significant difference was found in percentage depth dose beyond dmax for both energies. Beam profiles show that flatness and symmetry for both the energies were less than the 3%. The penumbra width for 'MLC only' field was more than that for the other two field-defining methods by (1 to 2 mm) and (0.8 to 1.3 mm) for 6-MV and 18-MV photon beams respectively. Analysis of 'width of 50% dose level' of the beam profiles at dmax to reflect the field size shows 1 to 2 mm more for 6-MV photons and 2.2 to 2.4 mm morefor 18-MV photons for 'MLC only' fields. The results of this study suggest that the characteristics of 120-leaf MLC system with 6 MV and 18 MV are same in all aspects except the surface dose, penumbra, dose in the build-up region, and width of 50% dose levels.

Respiratory gated simultaneous integrated boost-intensity modulated radiotherapy (SIB-IMRT) after breast conservative surgery for carcinoma of the breast: The Salmaniya Medical complex experience

PURPOSE

To present our clinical experience using SIB-IMRT Technique for Intact Breast cancer.

MATERIALS AND METHODS

A retrospective review of 45 cases of Stage I-IV breast cancer patients treated with SIB-IMRT with respiratory gating after Conservative treatments from 25th November 2008 to 16th February 2010. The most common fractionation was 1.8 Gy to Ipsilateral Breast tissue and 2.2 Gy to the lumpectomy cavity giving whole breast dose as 50.4 Gy and Lumpectomy cavity dose as 61.6 Gy over 28 fractions concomitantly. Respiratory gating was done and CT-images were taken in inspiratory breath hold position.

RESULTS

A total of 45 patients with breast cancer - stage I (17.7%), II (71%), III (8.9%), IV (2.2%) were treated with SIB- IMRT with respiratory gated radiotherapy. Out of 45 patients, 24 are of left sided breast cancer and 21 are of right sided breast cancer patients. The median, Dose maximum (D-max) in SIB-IMRT is 106.2% of prescribed lumpectomy site dose. The median isodose line prescribed to PTV-2 is 100%. The Conformity index (CI) is 0.9688 (median value) and Homogeneity index (HI) 1.06 (median). The median ipsilateral lung, mean dose is 21.66 Gy and V-20 is 37.4%. For left sided cases the median value of mean heart dose, V-30 and V-40 are 22.98 Gy, 23.45% and 9.45 % respectively. Acute skin toxicity was of Grade-I in 2.2 %, Grade-II in 64.4 %, Grade-III in 31.1 %, and Grade-IV in 2.2 %. The global Breast cosmoses were seen excellent in majority (93%) of case at median follow up of 8 months duration.

CONCLUSIONS

Breast SIB-IMRT Technique is feasible and comparable with other treatment techniques with reduced treatment duration by six fractions. At median follow up of 8 months the skin toxicity and cosmoses are excellent in high percentage of cases.

Does kV-MV dual-energy computed tomography have an advantage in determining proton stopping power ratios in patients?

Conventional kilovoltage (kV) x-ray-based dual-energy CT (DECT) imaging using two different x-ray energy spectra is sensitive to image noise and beam hardening effects. The purpose of this study was to evaluate the theoretical advantage of the DECT method for determining proton stopping power ratios (SPRs) using a combination of kV and megavoltage (MV) x-ray energies. We investigated three representative x-ray energy pairs: 100 and 140 kVp comprised the kV-kV pair, 100 kVp and 1 MV comprised the kV-MV pair, and two 1 MV x-ray beams-one with and one without external filtration-comprised the MV-MV pair. The SPRs of 34 human tissues were determined using the DECT method with these three x-ray energy pairs. Small perturbations were introduced into the CT numbers and x-ray spectra used for the DECT calculation to simulate the effects of random noise and beam hardening. An error propagation analysis was performed on the DECT calculation algorithm to investigate the propagation of CT number uncertainty to final SPR estimation and to suggest the best x-ray energy combination. We found that the DECT method using each of the three beam pairs achieved similar accuracy in determining the SPRs of human tissues in ideal conditions. However, when CT number uncertainties and artifacts such as imaging noise and beam hardening effects were considered, the kV-MV DECT improved the accuracy of SPR estimation substantially over the kV-kV or MV-MV DECT methods. Furthermore, our error propagation analysis showed that the combination of 100 kVp and 1 MV beams was close to the optimal selection when using the DECT method to determine SPRs. Overall, the kV-MV combination makes the DECT method more robust in resolving the effective atomic numbers for biological tissues than the traditional kV-kV DECT method.