On mixed electron–photon radiation therapy optimization using the column generation approach

Quality assurance for mixed electron-photon beam radiation therapy using treatment log files and MapCHECK

BACKGROUND

Mixed photon-electron beam radiotherapy (MBRT) is a technique that combines the use of both photons and electrons in one single treatment plan to exploit their advantageous and complimentary characteristics. Compared to other photon treatment modalities, it has been shown that the MBRT technique contributes to better target coverage and organ-at-risk (OAR) sparing. However, the use of combined photons and electrons in one delivery makes the technique more complex and a well-established quality assurance (QA) protocol for MBRT is essential.

PURPOSE

To investigate the feasibility of using MapCHECK and log file-dose reconstruction for MBRT plan verification and to recommend a patient-specific quality assurance (PSQA) protocol for MBRT.

METHODS

MBRT plans were robustly optimized for five soft-tissue sarcoma (STS) patients. Each plan comprised step-and-shoot deliveries of a six MV photon beam and a combination of five electron beam energies at an SAD of 100 cm. The plans were delivered to the MapCHECK device with collapsed gantry angle and the 2D dose distributions at the detector depth were measured. To simulate the expected dose distribution delivered to the MapCHECK, a MapCHECK computational phantom was modeled in EGSnrc based on vendor-supplied blueprint information. The dose to the detectors in the model was scored using the DOSXYZnrc user code. The agreement between the measured and the simulated dose distribution was evaluated using 2D gamma analysis with a gamma criterion of 3%/2 mm and a low dose threshold of 10%. One of the plans was selected and delivered with a rotating gantry angle for trajectory log file collection. To evaluate the potential interlinac and intralinac differences, the plan was delivered repeatedly on three linacs. From the collected log files, delivery parameters were retrieved to recalculate the 3D dose distributions in the patient's anatomy with DOSXYZnrc. The recalculated mean dose to the clinical target volume (CTV) and OARs from all deliveries were computed and compared with the planned dose in terms of percentage difference. To validate the accuracy of log file-based QA, the log file-recalculated dose was also compared with film measurement.

RESULTS

The agreement of the total dose distribution between the MapCHECK measurement and simulation showed gamma passing rates of above 97% for all five MBRT plans. In the log file-dose recalculation, the difference between the recalculated and the planned dose to the CTV and OARs was below 1% for all deliveries. No significant inter- or intralinac differences were observed. The log file-dose had a gamma passing rate of 98.6% compared to film measurement.

CONCLUSION

Both the MapCHECK measurements and log file-dose recalculations showed excellent agreement with the expected dose distribution. This study demonstrates the potential of using MapCHECK and log files as MBRT QA tools.

Delivery time reduction for mixed photon-electron radiotherapy by using photon MLC collimated electron arcs

Objective. Electron arcs in mixed-beam radiotherapy (Arc-MBRT) consisting of intensity-modulated electron arcs with dynamic gantry rotation potentially reduce the delivery time compared to mixed-beam radiotherapy containing electron beams with static gantry angle (Static-MBRT). This study aims to develop and investigate a treatment planning process (TPP) for photon multileaf collimator (pMLC) based Arc-MBRT.Approach. An existing TPP for Static-MBRT plans is extended to integrate electron arcs with a dynamic gantry rotation and intensity modulation using a sliding window technique. The TPP consists of a manual setup of electron arcs, and either static photon beams or photon arcs, shortening of the source-to-surface distance for the electron arcs, initial intensity modulation optimization, selection of a user-defined number of electron beam energies based on dose contribution to the target volume and finally, simultaneous photon and electron intensity modulation optimization followed by full Monte Carlo dose calculation. Arc-MBRT plans, Static-MBRT plans, and photon-only plans were created and compared for four breast cases. Dosimetric validation of two Arc-MBRT plans was performed using film measurements.Main results. The generated Arc-MBRT plans are dosimetrically similar to the Static-MBRT plans while outperforming the photon-only plans. The mean heart dose is reduced by 32% on average in the MBRT plans compared to the photon-only plans. The estimated delivery times of the Arc-MBRT plans are similar to the photon-only plans but less than half the time of the Static-MBRT plans. Measured and calculated dose distributions agree with a gamma passing rate of over 98% (3% global, 2 mm) for both delivered Arc-MBRT plans.Significance. A TPP for Arc-MBRT is successfully developed and Arc-MBRT plans showed the potential to improve the dosimetric plan quality similar as Static-MBRT while maintaining short delivery times of photon-only treatments. This further facilitates integration of pMLC-based MBRT into clinical practice.

Robust mixed electron-photon radiation therapy planning for soft tissue sarcoma

BACKGROUND

Mixed electron-photon beam radiation therapy (MBRT) is an emerging technique in which external electron and photon beams are simultaneously optimized into a single treatment plan. MBRT exploits the steep dose falloff and high surface dose of electrons while maintaining target conformity by leveraging the sharp penumbra of photons.

PURPOSE

This study investigates the dosimetric benefits of MBRT for soft tissue sarcoma (STS) patients.

MATERIAL AND METHODS

A retrospective cohort of 22 STS of the lower extremity treated with conventional photon-based Volumetric Modulated Arc Therapy (VMAT) were replanned with MBRT. Both VMAT and MBRT treatments were planned on the Varian TrueBeam linac using the Millenium multi-leaf collimator. No electron applicator, cutout or additional collimating devices were used for electron beams of MBRT plans. MBRT plans were optimized to use a combination of 6 MV photons and five electron energies (6, 9, 12, 16, 20 MeV) by a robust column generation algorithm. Electron beams in this study were planned at standard 100 cm source-axis distance (SAD). The dose to the clinical target volume (CTV), bone, normal tissue strip and other organs-at-risk (OARs) were compared using a Wilcoxon signed-rank test.

RESULTS

As part of the original VMAT treatment, tissue-equivalent bolus was required in 10 of the 22 patients. MBRT plans did not require bolus by virtue of the higher electron entrance dose. CTV coverage by the prescription dose was found to be clinically equivalent between plans of either modality:V 50Gy $V_{\text{50Gy}}$ (MBRT) = 97.9 ± 0.2% versusV 50Gy $V_{\text{50Gy}}$ (VMAT) = 98.1 ± 0.6% (p=0.34). Evaluating the absolute paired difference between doses to OARs in MBRT and VMAT plans, we observed lowerV 20Gy $V_{\text{20Gy}}$ to normal tissue in MBRT plans by 14.9 ± 3.2% (p < 10 - 6 $p<10^{-6}$ ). Similarly,V 50Gy $V_{\text{50Gy}}$ to bone was found to be decreased by 8.2 ± 4.0% (p < 10 - 3 $p<10^{-3}$ ) of the bone volume.

CONCLUSION

For STS with subcutaneous involvement, MBRT offers statistically significant sparing of OARs without sacrificing target coverage when compared to VMAT. MBRT plans are deliverable on conventional linacs without the use of electron applicators, shortened source-to-surface distance (SSD) or bolus. This study shows that MBRT is a logistically feasible technique with clear dosimetric benefits.

Applying the column generation method to the intensity modulated high dose rate brachytherapy inverse planning problem

Objective.Intensity modulated high dose rate brachytherapy (IMBT) is a rapidly developing application of brachytherapy where anisotropic dose distributions can be produced at each source dwell position. This technique is made possible by placing rotating metallic shields inside brachytherapy needles or catheters. By dynamically directing the radiation towards the tumours and away from the healthy tissues, a more conformal dose distribution can be obtained. The resulting treatment planning involves optimizing dwell position and shield angle (DPSA). The aim of this study was to investigate the column generation method for IMBT treatment plan optimization.Approach.A column generation optimization algorithm was developed to optimize the dwell times and shield angles. A retrospective study was performed on 10 prostate cases using RapidBrachyMCTPS. At every iteration, the plan was optimized with the chosen DPSA which would best improve the cost function that was added to the plan. The optimization process was stopped when the remaining DPSAs would not add value to the plan to limit the plan complexity.Main results.The average number of DPSAs and voxels were 2270 and 7997, respectively. The column generation approach yielded near-optimal treatment plans by using only 11% of available DPSAs on average in ten prostate cases. The coverage and organs at risk constraints passed in all ten cases.Significance.The column generation method produced high-quality deliverable prostate IMBT plans. The treatment plan quality reached a plateau, where adding more DPSAs had a minimal effect on dose volume histogram parameters. The iterative nature of the column generation method allows early termination of the treatment plan creation process as soon as the dosimetric indices from dose volume histogram satisfy the clinical requirements or if their values stabilize.

Auto-commissioning of a Monte Carlo electron beam model with application to photon MLC shaped electron fields

Objective.Presently electron beam treatments are delivered using dedicated applicators. An alternative is the usage of the already installed photon multileaf collimator (pMLC) enabling efficient electron treatments. Currently, the commissioning of beam models is a manual and time-consuming process. In this work an auto-commissioning procedure for the Monte Carlo (MC) beam model part representing the beam above the pMLC is developed for TrueBeam systems with electron energies from 6 to 22 MeV.Approach.The analytical part of the electron beam model includes a main source representing the primary beam and a jaw source representing the head scatter contribution each consisting of an electron and a photon component, while MC radiation transport is performed for the pMLC. The auto-commissioning of this analytical part relies on information pre-determined from MC simulations, in-air dose profiles and absolute dose measurements in water for different field sizes and source to surface distances (SSDs). For validation calculated and measured dose distributions in water were compared for different field sizes, SSDs and beam energies for eight TrueBeam systems. Furthermore, a sternum case in an anthropomorphic phantom was considered and calculated and measured dose distributions were compared at different SSDs.Main results.Instead of the manual commissioning taking up to several days of calculation time and several hours of user time, the auto-commissioning is carried out in a few minutes. Measured and calculated dose distributions agree generally within 3% of maximum dose or 2 mm. The gamma passing rates for the sternum case ranged from 96% to 99% (3% (global)/2 mm criteria, 10% threshold).Significance.The auto-commissioning procedure was successfully implemented and applied to eight TrueBeam systems. The newly developed user-friendly auto-commissioning procedure allows an efficient commissioning of an MC electron beam model and eases the usage of advanced electron radiotherapy utilizing the pMLC for beam shaping.

A novel stochastic optimization method for handling misalignments of proton and photon doses in combined treatments

Objective. Combined proton–photon treatments, where most fractions are delivered with photons and only a few are delivered with protons, may represent a practical approach to optimally use limited proton resources. It has been shown that, when organs at risk (OARs) are located within or near the tumor, the optimal multi-modality treatment uses protons to hypofractionate parts of the target volume and photons to achieve near-uniform fractionation in dose-limiting healthy tissues, thus exploiting the fractionation effect. These plans may be sensitive to range and setup errors, especially misalignments between proton and photon doses. Thus, we developed a novel stochastic optimization method to directly incorporate these uncertainties into the biologically effective dose (BED)-based simultaneous optimization of proton and photon plans. Approach. The method considers the expected value E b and standard deviation σ b of the cumulative BED b in every voxel of a structure. For the target, a piecewise quadratic penalty function of the form b min − E b − 2 σ b + 2 is minimized, aiming for plans in which the expected BED minus two times the standard deviation exceeds the prescribed BED b min . Analogously, E b + 2 σ b − b max + 2 is considered for OARs. Main results. Using a spinal metastasis case and a liver cancer patient, it is demonstrated that the novel stochastic optimization method yields robust combined treatment plans. Tumor coverage and a good sparing of the main OARs are maintained despite range and setup errors, and especially misalignments between proton and photon doses. This is achieved without explicitly considering all combinations of proton and photon error scenarios. Significance. Concerns about range and setup errors for safe clinical implementation of optimized proton–photon radiotherapy can be addressed through an appropriate stochastic planning method.

Enabling non-isocentric dynamic trajectory radiotherapy by integration of dynamic table translations

Objective. The purpose of this study is to develop a treatment planning process (TPP) for non-isocentric dynamic trajectory radiotherapy (DTRT) using dynamic gantry rotation, collimator rotation, table rotation, longitudinal, vertical and lateral table translations and intensity modulation and to validate the dosimetric accuracy. Approach. The TPP consists of two steps. First, a path describing the dynamic gantry rotation, collimator rotation and dynamic table rotation and translations is determined. Second, an optimization of the intensity modulation along the path is performed. We demonstrate the TPP for three use cases. First, a non-isocentric DTRT plan for a brain case is compared to an isocentric DTRT plan in terms of dosimetric plan quality and delivery time. Second, a non-isocentric DTRT plan for a craniospinal irradiation (CSI) case is compared to a multi-isocentric intensity modulated radiotherapy (IMRT) plan. Third, a non-isocentric DTRT plan for a bilateral breast case is compared to a multi-isocentric volumetric modulated arc therapy (VMAT) plan. The non-isocentric DTRT plans are delivered on a TrueBeam in developer mode and their dosimetric accuracy is validated using radiochromic films. Main results. The non-isocentric DTRT plan for the brain case is similar in dosimetric plan quality and delivery time to the isocentric DTRT plan but is expected to reduce the risk of collisions. The DTRT plan for the CSI case shows similar dosimetric plan quality while reducing the delivery time by 45% in comparison with the IMRT plan. The DTRT plan for the breast case showed better treatment plan quality in comparison with the VMAT plan. The gamma passing rates between the measured and calculated dose distributions are higher than 95% for all three plans. Significance. The versatile benefits of non-isocentric DTRT are demonstrated with three use cases, namely reduction of collision risk, reduced setup and delivery time and improved dosimetric plan quality.

FLASH radiotherapy treatment planning and models for electron beams

The FLASH effect designates normal tissue sparing at ultra-high dose rate (UHDR, >40 Gy/s) compared to conventional dose rate (∼0.1 Gy/s) irradiation while maintaining tumour control and has the potential to improve the therapeutic ratio of radiotherapy (RT). UHDR high-energy electron (HEE, 4-20 MeV) beams are currently a mainstay for investigating the clinical potential of FLASH RT for superficial tumours. In the future very-high energy electron (VHEE, 50-250 MeV) UHDR beams may be used to treat deep-seated tumours. UHDR HEE treatment planning focused at its initial stage on accurate dosimetric modelling of converted and dedicated UHDR electron RT devices for the clinical transfer of FLASH RT. VHEE treatment planning demonstrated promising dosimetric performance compared to clinical photon RT techniques in silico and was used to evaluate and optimise the design of novel VHEE RT devices. Multiple metrics and models have been proposed for a quantitative description of the FLASH effect in treatment planning, but an improved experimental characterization and understanding of the FLASH effect is needed to allow for an accurate and validated modelling of the effect in treatment planning. The importance of treatment planning for electron FLASH RT will augment as the field moves forward to treat more complex clinical indications and target sites. In this review, TPS developments in HEE and VHEE are presented considering beam models, characteristics, and future FLASH applications.

A hybrid column generation and simulated annealing algorithm for direct aperture optimization

The purpose of this work was to develop a hybrid column generation (CG) and simulated annealing (SA) algorithm for direct aperture optimization (H-DAO) and to show its effectiveness in generating high quality treatment plans for intensity modulated radiation therapy (IMRT) and mixed photon-electron beam radiotherapy (MBRT). The H-DAO overcomes limitations of the CG-DAO with two features improving aperture selection (branch-feature) and enabling aperture shape changes during optimization (SA-feature). The H-DAO algorithm iteratively adds apertures to the plan. At each iteration, a branch is created for each field provided. First, each branch determines the most promising aperture of its assigned field and adds it to a copy of the current apertures. Afterwards, the apertures of each branch undergo an MU-weight optimization followed by an SA-based simultaneous shape and MU-weight optimization and a second MU-weight optimization. The next H-DAO iteration continues the branch with the lowest objective function value. IMRT and MBRT treatment plans for an academic, a brain and a head and neck case generated using the CG-DAO and H-DAO were compared. For every investigated case and both IMRT and MBRT, the H-DAO leads to a faster convergence of the objective function value with number of apertures compared to the CG-DAO. In particular, the H-DAO needs about half the apertures to reach the same objective function value as the CG-DAO. The average aperture areas are 27% smaller for H-DAO than for CG-DAO leading to a slightly larger discrepancy between optimized and final dose. However, a dosimetric benefit remains. The H-DAO was successfully developed and applied to IMRT and MBRT. The faster convergence with number of apertures of the H-DAO compared to the CG-DAO allows to select a better compromise between plan quality and number of apertures.

Graph neural networks and deep reinforcement learning for simultaneous beam orientation and trajectory optimization of Cyberknife

Objective. Despite the high-quality treatment, the long treatment time of the Cyberknife system is believed to be a drawback. The high flexibility of its robotic arm requires meticulous path-finding algorithms to deliver the prescribed dose in the shortest time.Approach. We proposed a Deep Q-learning based on Graph Neural Networks to find the subset of the beams and the order to traverse them. A complex reward function is defined to minimize the distance covered by the robotic arm while avoiding the selection of close beams. Individual beam scores are also generated based on their effect on the beam intensity and are incorporated in the reward function. Main results. The performance of the presented method is evaluated on three clinical cases suffering from lung cancer. Applying this approach leads to an average of 35% reduction in the treatment time while delivering the prescribed dose provided by the physicians.Significance. Shorter treatment times result in a better treatment experience for individual patients, reduces discomfort and the sides effects of inadvertent movements for them. Additionally, it creates the opportunity to treat a higher number of patients in a given time period at the radiation therapy centers.

Implementation and experimental validation of a robust hybrid direct aperture optimization approach for mixed-beam radiotherapy

Purpose

The objectives of the work presented in this paper were to (1) implement a robust‐optimization method for deliverable mixed‐beam radiotherapy (MBRT) plans within a previously developed MBRT planning framework; (2) perform an experimental validation of the delivery of robust‐optimized MBRT plans; and (3) compare PTV‐based and robust‐optimized MBRT plans in terms of target dose robustness and organs at risk (OAR) sparing for clinical head and neck and brain patient cases.

Methods

A robust‐optimization method, which accounts for translational setup errors, was implemented within a previously developed treatment planning framework for MBRT. The framework uses a hybrid direct aperture optimization method combining column generation and simulated annealing. A robust plan was developed and then delivered to an anthropomorphic head phantom using the Developer Mode of a TrueBeam linac. Planar dose distributions were measured and compared to the planned dose. Robust‐optimized and PTV‐based plans were developed for three clinical patient cases consisting of two head and neck cases and one brain case. The plans were compared in terms of the robustness to 5 mm shifts of the target volume dose as well as in terms of OAR sparing.

Results

Using a gamma criterion of 3%/2 mm and a dose threshold of 10%, the agreement between film measurements and dose calculations was better than 97.7% for the total plan and better than 95.5% for the electron component of the plan. For the two head and neck patient cases, the average clinical target volume (CTV) dose homogeneity index (V95%–V107%) over all the considered setup error scenarios was on average 19% lower for the PTV‐based plans and it had a larger standard deviation. The robust‐optimized plans achieved, on average, a 20% reduction in the OAR doses compared to the PTV‐based plans. For the brain patient case, the CTV dose homogeneity index was similar for the two plans, while the OAR doses were 22% lower, on average, for the robust‐optimized plan. No clear trend in terms of electron contributions was found across the three patient cases, although robust‐optimized plans tended toward higher electron beam energies.

Conclusions

A framework for robust optimization of deliverable MBRT plans has been developed and validated. PTV‐based MBRT were found to not be robust to setup errors, while the dose delivered by the robust‐optimized plans were clinically acceptable for all considered error scenarios and had better OAR sparing. This study shows that the robust optimization is a promising alternative to conventional PTV margins for MBRT.

Ion chamber and film-based quality assurance of mixed electron-photon radiation therapy

PURPOSE

In previous work, we demonstrated that mixed electron-photon radiation therapy (MBRT) produces treatment plans with improved normal tissue sparing and similar target coverage, when compared to photon-only plans. The purpose of this work was to validate the MBRT delivery process on a Varian TrueBeam accelerator and laying the groundwork for a patient-specific quality assurance (QA) protocol based on ion chamber point measurements and 2D film measurements.

METHODS

MC beam models used to calculate the MBRT dose distributions of each modality (photons/electrons) were validated with a single-angle beam MBRT treatment plan delivered on a slab of Solid Water phantom with a film positioned at a depth of 2 cm. The measured film absorbed dose was compared to the calculated dose. To validate clinical deliveries, a polymethyl methacrylate (PMMA) cylinder was machined and holes were made to fit an ionization chamber. A complex MBRT plan involving a photon arc and three electron delivery angles was created with the aim of reproducing a clinically realistic dose distribution in typical soft tissue sarcoma tumours of the extremities. The treatment plan was delivered on the PMMA cylinder. Point measurements were taken with an Exradin A1SL chamber at two nominal depths: 1.4 cm and 2.1 cm. The plan was also delivered on a second identical phantom with an insert at 2 cm depth, where a film was placed. An existing EGSnrc user-code, SPRRZnrc, was modified to calculate the stopping power ratios between any materials in the same voxelized geometry used for dose calculation purposes. This modified code, called SPRXYZnrc, was used to calculate a correction factor, k MBRT , accounting for the differences in electron fluence spectrum at the measurement point compared to that at reference conditions. The uncertainty associated with neglecting potential ionization chamber fluence perturbation correction factors using this approach was estimated.

RESULTS

The film measurement from the Solid Water phantom treatment plan was in good agreement with the simulated dose distribution, with a gamma pass rate of 96.1% for a 3%/2 mm criteria. For the PMMA phantom delivery, for the same gamma criteria, the pass rate was 97.3%. The ion chamber measurements of the total delivered dose agreed with the MC-simulated dose within 2.1%. The beam quality correction factors amounted to, at most, a 4% correction on the ion chamber measurement. However, individual contribution of low electron energies proved difficult to precisely measure due to their steep dose gradients, with disagreements of up to 28% ± 15% at 2.1 cm depth (6 MeV). Ion chamber measurement procedure of electron beams was achieved in less than 5 min, and the entire validation process including phantom setup was performed in less than 30 min.

CONCLUSION

The agreement between measured and simulated MBRT doses indicates that the dose distributions obtained from the MBRT treatment planning algorithm are realistically achievable. The SPRXYZnrc MC code allowed for convenient calculations of k MBRT simultaneously with the dose distributions, laying the groundwork for patient-specific QA protocol practical for clinical use. Further investigation is needed to establish the accuracy of our ionization chamber correction factors k MBRT calculations at low electron energies.

Photon beam energy dependent single-arc volumetric modulated arc optimization

PURPOSE

The purpose of this work was to present a new single-arc mixed photon (6&18MV) VMAT (SAMP) optimization framework that concurrently optimizes for two photon energies with corresponding partial arc lengths.

METHODS AND MATERIALS

Owing to simultaneous optimization of energy dependent intensity maps and corresponding arc locations, the proposed model poses nonlinearity. Unique relaxation constraints based on McCormick approximations were introduced for linearization. Energy dependent intensity maps were then decomposed to generate apertures. Feasibility of the proposed framework was tested on a sample of ten prostate cancer cases with lateral separation ranging from 34 cm (case no.1) to 52 cm (case no.6). The SAMP plans were compared against single energy (6MV) VMAT (SE) plans through dose volume histograms (DVHs) and radiobiological parameters including normal tissue complication probability (NTCP) and equivalent uniform dose (EUD).

RESULTS

The contribution of higher energy photon beam optimized by the algorithm demonstrated an increase for cases with a lateral separation >40 cm. SAMP-VMAT notably improved bladder and rectum sparing in large size cases. Compared to single energy, SAMP-VMAT plans reduced bladder and rectum NTCP in cases with large lateral separation. With the exception of one case, SAMP-VMAT either improved or maintained femoral heads compared to SE-VMAT. SAMP-VMAT reduced the nontarget tissue integral dose in all ten cases.

CONCLUSIONS

A single-arc VMAT optimization framework comprising mixed photon energy partial arcs was presented. Overall results underline the feasibility and potential of the proposed approach for improving OAR sparing in large size patients without compromising the target homogeneity and coverage.

TriB-RT: Simultaneous optimization of photon, electron and proton beams

PURPOSE

To develop a novel treatment planning process (TPP) with simultaneous optimization of modulated photon, electron and proton beams for improved treatment plan quality in radiotherapy.

METHODS

A framework for fluence map optimization of Monte Carlo (MC) calculated beamlet dose distributions is developed to generate treatment plans consisting of photon, electron and spot scanning proton fields. Initially, in-house intensity modulated proton therapy (IMPT) plans are compared to proton plans created by a commercial treatment planning system (TPS). A triple beam radiotherapy (TriB-RT) plan is generated for an exemplary academic case and the dose contributions of the three particle types are investigated. To investigate the dosimetric potential, a TriB-RT plan is compared to an in-house IMPT plan for two clinically motivated cases. Benefits of TriB-RT for a fixed proton beam line with a single proton field are investigated.

RESULTS

In-house optimized IMPT are of at least equal or better quality than TPS-generated proton plans, and MC-based optimization shows dosimetric advantages for inhomogeneous situations. Concerning TriB-RT, for the academic case, the resulting plan shows substantial contribution of all particle types. For the clinically motivated case, improved sparing of organs at risk close to the target volume is achieved compared to IMPT (e.g. myelon and brainstem [Formula: see text] -37%) at cost of an increased low dose bath (healthy tissue V _10% +22%). In the scenario of a fixed proton beam line, TriB-RT plans are able to compensate the loss in degrees of freedom to substantially improve plan quality compared to a single field proton plan.

CONCLUSION

A novel TPP which simultaneously optimizes photon, electron and proton beams was successfully developed. TriB-RT shows the potential for improved treatment plan quality and is especially promising for cost-effective single-room proton solutions with a fixed beamline in combination with a conventional linac delivering photon and electron fields.

Simultaneous trajectory generation and volumetric modulated arc therapy optimization

PURPOSE

Trajectory-based treatment planning involves the combination of a gantry-couch trajectory with volumetric modulated arc therapy (VMAT) treatment plan optimization. This work presents the implementation of an optimization methodology that generates a trajectory simultaneous with treatment plan optimization (simTr-VMAT).

METHODS

The optimization algorithm is based on the column generation approach, in which a treatment plan is iteratively constructed through the solution of a subproblem called the "pricing problem." The property of the pricing problem to rank candidate apertures based on their associated price is leveraged to select an optimal aperture while simultaneously determining the trajectory path. A progressively increasing gantry-couch grid resolution is used to provide an initial coarse sampling of the angular solution space while maintaining fine control point spacing with the final treatment plan. The trajectory optimization was applied and compared to coplanar VMAT treatment plans for a lung patient, a glioblastoma patient, and a prostate patient. Algorithm validation was performed through the generation of 5000 random trajectories and optimization using column generation VMAT for each patient case, representing the solution space for the trajectory optimization problem. The simTr-VMAT trajectories were compared against these random trajectories based on a quality metric that prefers trajectories with few control points and low objective function value over long, inefficient trajectories.

RESULTS

For the lung patient, the simTr-VMAT plan resulted in a decrease of the mean dose of 1.5 and 1.0 Gy to the heart and ipsilateral lung, respectively. For the glioblastoma patient, the simTr-VMAT plan resulted in improved planning target volume coverage with a decrease in mean dose to the eyes, lens, nose, and contralateral temporal lobe between 2 and 7 Gy. The prostate patient showed no clinically relevant dosimetric improvement. The simTr-VMAT treatment plans ranked at the 99.6, 96.3, and 99.4 percentiles compared to the distribution of randomly generated trajectories for the lung, glioblastoma, and prostate patients, respectively.

CONCLUSION

The simTr-VMAT optimization methodology resulted in treatment plans with equivalent or improved dosimetric outcomes compared to coplanar VMAT treatment plans, with the trajectories resulting from the optimization ranking among the optimal trajectories for each patient case.

Trajectory-based VMAT for cranial targets with delivery at shortened SAD

INTRODUCTION

Trajectory-based volumetric modulated arc therapy (tr-VMAT) treatment plans enable the option for noncoplanar delivery yielding steeper dose gradients and increased sparing of critical structures compared to conventional treatment plans. The addition of translational couch motion to shorten the effective source-to-axis distance (SAD) may result in improved delivery precision and an increased effective dose rate. In this work, tr-VMAT treatment plans using a noncoplanar "baseball stitch" trajectory were implemented, applied to patients presented with cranial targets, and compared to the clinical treatment plans.

METHODS

A treatment planning workflow was implemented: (a) beamlet doses were calculated for control points defined along a baseball stitch trajectory using a collapsed-cone convolution-superposition algorithm; (b) VMAT treatment plans were optimized using the column generation approach; (c) a final dose distribution was calculated in Varian Eclipse using the analytical anisotropic algorithm by importing the optimized treatment plan parameters. Tr-VMAT plans were optimized for ten patients presented with cranial targets at both standard and shortened SAD, and compared to the clinical treatment plans through isodose distributions, dose-volume histograms, and dosimetric indices. The control point specifications of the optimized tr-VMAT plans were used to estimate the delivery time.

RESULTS

The optimized tr-VMAT plans with both shortened and standard SAD delivery yielded a comparable plan quality to the clinical treatment plans. A statistically significant benefit was observed for dose gradient index and monitor unit efficiency for shortened SAD tr-VMAT plans, while improved target volume conformity was observed for the clinical treatment plan (P ≤ 0.05). A clear dosimetric benefit was not demonstrated between tr-VMAT delivery at shortened SAD compared to standard SAD, but shortened SAD delivery yielded a fraction size-dependent reduction in the estimated delivery time.

CONCLUSION

The implementation of "baseball stitch" tr-VMAT treatment plans to patients presented with cranial targets demonstrated comparable plan quality to clinical treatment plans. The delivery at shortened SAD produced a fraction size-dependent decrease in estimated delivery time.

Simultaneous optimization of mixed photon energy beams in volumetric modulated arc therapy

PURPOSE

Despite the availability of multiple energy photon beams on clinical linear accelerators, volumetric modulated arc therapy (VMAT) optimization is currently limited to a single photon beam. The purpose of this work was to present a proof-of-principle study on an algorithm for simultaneous optimization of mixed photon beams for VMAT (MP - VMAT), utilizing an additional photon energy as an additional degree of freedom.

METHODS

The MP - VMAT optimization algorithm is presented as a two-step heuristic approach. First, a convex linear programming problem is solved for simultaneous optimization of nonuniform dual energy intensity maps (DEIMs) for an angular resolution of 36 equi-spaced beam segments. Subsequently, for a given gantry speed schedule, the second step aims to best replicate each DEIM by dispersing MP - VMAT apertures along with their corresponding intensities over their respective beam segment. This constitutes a nonlinear problem, which is linearized using McCormick relaxation. The final large-scale mixed integer linear programming (MILP) dispersion model ensures a contiguous and smooth transition of multileaf collimators (MLCs) from one beam segment to the next. To demonstrate the proof-of-principle, we first compared the quality of dose volume histograms (DVHs) of MP - VMAT to the ones calculated from 36 DEIMs following the step 1 of MP - VMAT model. Additionally, the MLCs motion violations were evaluated for the complete 360° gantry rotation for gantry speeds ranging from 1 to 6° per second. The quality of MP - VMAT plans were also compared to conventional single energy VMAT plans via DVH, homogeneity index (HI), and conformity number (CN) for two prostate cases.

RESULTS

The MP - VMAT model resulted in a successful convergence of DVHs relative to the ones from DEIMs with HI and CN of 0.05 and 0.9, respectively, for 1 and 2° per second gantry speed schedules. In replicating the DEIMs, the MILP dispersion model was able to achieve optimality for almost all segments at 1° per second and for majority of segments at 2° per second. Although, DVHs quality was slightly inferior for 3° per second gantry speed, the target conformity of 0.9 and heterogeneity of 0.08 were achievable even for the suboptimal solutions. No violations of the MLC constraints were observed throughout the complete 360 degree arc rotation for any gantry speed schedule, thereby confirming MILP dispersion model. For the two prostate cases, the results showed MP - VMAT's ability to achieve substantial dose reduction in rectum and bladder while yielding similar target coverage compared to single energy VMAT. Bladder volume was mostly spared in low-to-intermediate dose region. Rectal volume sparing (3 % to 12 %) was observed in the intermediate (from 25 to 50 Gy) dose region.

CONCLUSION

We demonstrate the first formalism of a large-scale simultaneous optimization of mixed photon energy beams for VMAT. Dosimetric comparison of MP - VMAT to single energy VMAT demonstrated potential advantages of using mixed photon energy beams for prostate plans, thus providing an impetus for further testing on a large clinical cohort.

Integrating DVH criteria into a column generation algorithm for VMAT treatment planning

Volumetric-modulated arc therapy (VMAT) treatment planning is an efficient treatment technique with a high degree of flexibility in terms of dose rate, gantry speed, and aperture shapes during rotation around the patient. However, the dynamic nature of VMAT results in a large-scale nonconvex optimization problem. Determining the priority of the tissues and voxels to obtain clinically acceptable treatment plans poses additional challenges for VMAT optimization. The main purpose of this paper is to develop an automatic planning approach integrating dose-volume histogram (DVH) criteria in direct aperture optimization for VMAT, by adjusting the model parameters during the algorithm. The proposed algorithm is based on column generation, an optimization technique that sequentially generates the apertures and optimizes the corresponding intensities. We take the advantage of iterative procedure in this method to modify the weight vector of the penalty function based on the DVH criteria and decrease the use of trial-and-error in the search for clinically acceptable plans. We evaluate the efficiency of the algorithm and treatment quality using a clinical prostate case and a challenging head-and-neck case. In both cases, we generate 15 random initial weight vectors to assess the robustness of the algorithm. In the prostate case, our methodology obtained clinically acceptable plans in all instances with only a 10% increase in the computational time, while simple VMAT optimization found just three acceptable plans. To have an idea with respect to the existing software, we compared the obtained DVH to a commercial software. The quality of the diagrams of the proposed method, especially for the healthy tissues, is significantly better while the computational time is less. In the head-and-neck case, 93.3% of the clinically acceptable plans are obtained while no plan was acceptable in simple VMAT. In sum, the results demonstrate the ability of the proposed optimization algorithm to obtain clinically acceptable plans without human intervention and also its robustness to weight parameters. Moreover, our proposed weight adjustment procedure proves to reduce the symmetry in the solution space and the time required for the post-optimization phase.

Inverse optimization of low-cost kilovoltage x-ray arc therapy plans

PURPOSE

The objective of this work was to investigate the benefits of using inverse optimization treatment planning for kilovoltage arc therapy (KVAT) and to assess the dosimetric limitations of KVAT.

METHODS

Monte Carlo (MC) calculated, inversely optimized KVAT plans of spherical, idealized breast, lung, and prostate lesions were calculated using the EGSnrc/BEAMnrc and DOSXYZnrc MC codes. The dose delivered with the KVAT system, which generates 200-225 kV photon beamlets, was calculated and inversely optimized using an optimization framework developed at McGill University. KVAT dose distributions were compared with inversely optimized and MC generated megavoltage (MV) volumetric modulated arc therapy (VMAT) plans as a reference. Prescription doses delivered to 95% of the planning target volume (PTV) were 38.5 (10 fractions), 60 (30 fractions) and 73.8 (41 fractions) Gy for the breast, lung and prostate patients, respectively. Dose distributions, dose volume histograms, and PTV homogeneity indices were used to evaluate KVAT and VMAT plans based on RTOG protocols.

RESULTS

All organ-at-risk (OAR) doses were within prescribed dose limits for KVAT and VMAT plans. Generally, KVAT plans delivered higher doses to OARs. For example, due to the lower energy of KVAT, 50% of the rib volume received 12.9 Gy from KVAT while only receiving 2.5 Gy from VMAT. OAR doses were especially high for the KVAT prostate plan due to the presence of large volumes of bony anatomy, which illustrates a limitation of the KVAT system. The KVAT treatment times per fraction for the breast, lung and prostate patients were 2.8, 2.6 and 5.5 min, respectively.

CONCLUSIONS

The inversely optimized KVAT plans presented in this work have demonstrated the ability of our novel low-cost, kilovoltage x-ray therapy system to safely treat deep-seated spherical lesions in breast and lung patients while meeting RTOG dose constraints on OARs.

Dosimetric evaluation of a novel electron–photon mixed beam, produced by a medical linear accelerator

Aim

This study deals with the characteristics of simultaneous photon and electron beams in homogenous and inhomogeneous phantoms by experimental and Monte Carlo dosimetry, for therapeutic purposes. Materials and methods: Both 16 and 20 MeV high-energy electron beams were used as the original beam to strike perforated lead sheets to produce the mixed beam. The dosimetry results were achieved by measurement in an ion chamber in a water phantom and film dosimetry in a Perspex nasal phantom, and then compared with those calculated through a simulation approach. To evaluate two-dimensional dose distribution in the inhomogeneous medium, the dose–area histogram was obtained.

Results

The highest percentage of photon contribution in mixed beam was found to be 36% for 2-mm thickness of lead layer with holes diameter of 0·2 cm for a 20 MeV primary electron energy. For small fields, the percentage depth dose parameters variations were found to be similar to pure electron beam within ±2%. The most feasible flatness in beam profile was 11% for pure electron and 7% for the mixed beam. Penumbra changes as function of depth was about ten times better than in pure electron field.

Conclusions

The results present some dosimetric advantages that can make this study a platform for the production of simultaneous mixed beams in future linear accelerators (LINACs), which through redesign of the LINAC head, which could lead to setup error reduction and a decrease of intra-fractional tumour cells repair.

Part 2: Dynamic mixed beam radiotherapy (DYMBER): Photon dynamic trajectories combined with modulated electron beams

PURPOSE

The purpose of this study was to develop a treatment technique for dynamic mixed beam radiotherapy (DYMBER) utilizing increased degrees of freedom (DoF) of a conventional treatment unit including different particle types (photons and electrons), intensity and energy modulation and dynamic gantry, table, and collimator rotations.

METHODS

A treatment planning process has been developed to create DYMBER plans combining photon dynamic trajectories (DTs) and step and shoot electron apertures collimated with the photon multileaf collimator (pMLC). A gantry-table path is determined for the photon DTs with minimized overlap of the organs at risk (OARs) with the target. In addition, an associated dynamic collimator rotation is established with minimized area between the pMLC leaves and the target contour. pMLC sequences of photon DTs and electron pMLC apertures are then simultaneously optimized using direct aperture optimization (DAO). Subsequently, the final dose distribution of the electron pMLC apertures is calculated using the Swiss Monte Carlo Plan (SMCP). The pMLC sequences of the photon DTs are then re-optimized with a finer control point resolution and with the final electron dose distribution taken into account. Afterwards, the final photon dose distribution is calculated also using the SMCP and summed together with the one of the electrons. This process is applied for a brain and two head and neck cases. The resulting DYMBER dose distributions are compared to those of dynamic trajectory radiotherapy (DTRT) plans consisting only of photon DTs and clinically applied VMAT plans. Furthermore, the deliverability of the DYMBER plans is verified in terms of dosimetric accuracy, delivery time and collision avoidance. For this purpose, The DYMBER plans are delivered to Gafchromic EBT3 films placed in an anthropomorphic head phantom on a Varian TrueBeam linear accelerator.

RESULTS

For each case, the dose homogeneity in the target is similar or better for DYMBER compared to DTRT and VMAT. Averaged over all three cases, the mean dose to the parallel OARs is 16% and 28% lower, D2% to the serial OARs is 17% and 37% lower and V10% to normal tissue is 12% and 4% lower for the DYMBER plans compared to the DTRT and VMAT plans, respectively. The DYMBER plans are delivered without collision and with a 4-5 min longer delivery time than the VMAT plans. The absolute dose measurements are compared to calculation by gamma analysis using 2% (global)/2 mm criteria with passing rates of at least 99%.

CONCLUSIONS

A treatment technique for DYMBER has been successfully developed and verified for its deliverability. The dosimetric superiority of DYMBER over DTRT and VMAT indicates utilizing increased DoF to be the key to improve brain and head and neck radiation treatments in future.