Trajectory optimization for dynamic couch rotation during volumetric modulated arc radiotherapy

Robust optimization and assessment of dynamic trajectory and mixed-beam arc radiotherapy: a preliminary study

Objective.Dynamic trajectory radiotherapy (DTRT) and dynamic mixed-beam arc therapy (DYMBARC) exploit non-coplanarity and, for DYMBARC, simultaneously optimized photon and electron beams. Margin concepts to account for set-up uncertainties during delivery are ill-defined for electron fields. We develop robust optimization for DTRT&DYMBARC and compare dosimetric plan quality and robustness for both techniques and both optimization strategies for four cases.Approach.Cases for different treatment sites and clinical target volume (CTV) to planning target volume (PTV) margins,m, were investigated. Dynamic gantry-table-collimator photon paths were optimized to minimize PTV/organ-at-risk (OAR) overlap in beam's-eye-view and minimize potential photon multileaf collimator (MLC) travel. For DYMBARC plans, non-isocentric partial electron arcs or static fields with shortened source-to-surface distance (80 cm) were added. Direct aperture optimization (DAO) was used to simultaneously optimize MLC-based intensity modulation for both photon and electron beams yielding deliverable PTV-based DTRT&DYMBARC plans. Robust-optimized plans used the same paths/arcs/fields. DAO with stochastic programming was used for set-up uncertainties with equal weights in all translational directions and magnitudeδsuch thatm= 0.7δ. Robust analysis considered random errors in all directions with or without an additional systematic error in the worst 3D direction for the adjacent OARs.Main results.Electron contribution was 7%-41% of target dose depending on the case and optimization strategy for DYMBARC. All techniques achieved similar CTV coverage in the nominal (no error) scenario. OAR sparing was overall better in the DYMBARC plans than in the DTRT plans and DYMBARC plans were generally more robust to the considered uncertainties. OAR sparing was better in the PTV-based than in robust-optimized plans for OARs abutting or overlapping with the target volume, but more affected by uncertainties.Significance.Better plan robustness can be achieved with robust optimization than with margins. Combining electron arcs/fields with non-coplanar photon trajectories further improves robustness and OAR sparing.

Organs-at-risk dose and normal tissue complication probability with dynamic trajectory radiotherapy (DTRT) for head and neck cancer

We compared dynamic trajectory radiotherapy (DTRT) to state-of-the-art volumetric modulated arc therapy (VMAT) for 46 head and neck cancer cases. DTRT had lower dose to salivary glands and swallowing structure, resulting in lower predicted xerostomia and dysphagia compared to VMAT. DTRT is deliverable on C-arm linacs with high dosimetric accuracy.

Noncoplanar Volumetric Modulated Arc Therapy for Hepatocellular Carcinoma Based on a Cage-Like Radiotherapy System: A Simulation Study

BACKGROUND

The incorporation of noncoplanar beam arrangements has been proposed in liver radiotherapy modalities, which can reduce the dose in normal tissues compared to coplanar techniques. Noncoplanar radiotherapy techniques for hepatocellular carcinoma treatment based on the Linac design have a limited effective arc angle to avoid collisions.

PURPOSE

To propose a novel noncoplanar volumetric modulated arc therapy technique based on a cage-like radiotherapy system and investigate its performance in hepatocellular carcinoma patients.

METHODS

The computed tomography was deflected 90° to meet the structure of a cage-like radiotherapy system and design the noncoplanar volumetric modulated arc therapy technique based on a cage-like radiotherapy system plan in the Pinnacle3 planning system. An noncoplanar volumetric modulated arc therapy technique based on a cage-like radiotherapy system plan was customized for each of 10 included hepatocellular carcinoma patients, with 6 dual arcs ranging from -30° to 30°. Six couch angles were set with an interval of 36° and distributed along with the longest diameter of planning target volume. The dosimetric parameters of noncoplanar volumetric modulated arc therapy technique based on a cage-like radiotherapy system plan were compared with the noncoplanar volumetric modulated arc therapy and volumetric modulated arc therapy plan.

RESULTS

The 3 radiotherapy techniques regarding planning target volume were statistically different for D98%, D2%, conformity index, and homogeneity index with χ² = 9.692, 14.600, 8.600, and 12.600, and P = .008, .001, .014, and .002, respectively. Further multiple comparisons revealed that noncoplanar volumetric modulated arc therapy technique based on a cage-like radiotherapy system significantly reduced the mean dose (P = .005) and V5 (P = .005) of the normal liver, the mean dose (P = .005) of the stomach, and V30 (P = .028) of the lung compared to noncoplanar volumetric modulated arc therapy. Noncoplanar volumetric modulated arc therapy technique based on a cage-like radiotherapy system significantly reduced the mean dose (P = .005) and V5 (P = .005) of the normal liver, the mean dose (P = .017) of the spinal cord, V50 (P = .043) of the duodenum, the maximum dose (P = .007) of the esophagus, and V30 (P = .047) of the whole lung compared to volumetric modulated arc therapy. The results indicate that noncoplanar volumetric modulated arc therapy technique based on a cage-like radiotherapy system protects the normal liver, stomach, and lung better than noncoplanar volumetric modulated arc therapy and protects the normal liver, spinal cord, duodenum, esophagus, and lung better than volumetric modulated arc therapy.

CONCLUSIONS

The noncoplanar volumetric modulated arc therapy technique based on a cage-like radiotherapy system technique with the arrangement of noncoplanar arcs provided optimal dosimetric gains compared with noncoplanar volumetric modulated arc therapy and volumetric modulated arc therapy, except for the heart. Noncoplanar volumetric modulated arc therapy technique based on a cage-like radiotherapy system should be considered in more clinically challenging cases.

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.

Sub-arc collimator angle optimization based on the conformity index heatmap for VMAT planning of multiple brain metastases SRS treatments

Objective

The aim of this study was to investigate the impact of collimator angle optimization in single-isocenter coplanar volume modulated arc therapy (VMAT) stereotactic radiosurgery (SRS) for multiple metastases with respect to dosimetric quality and treatment delivery efficiency. In particular, this is achieved by a novel algorithm of sub-arc collimator angle optimization (SACAO).

Methods

Twenty patients with multiple brain metastases were retrospectively included in this study. A multi-leaf collimator (MLC) conformity index (MCI) that is defined as the ratio of the area of target projection in the beam’s eye view (BEV) to the related area fitted by MLC was applied. Accordingly, for each control point, 180 MCI values were calculated with a collimator angle interval of 1°. A two-dimensional heatmap of MCI as a function of control point and collimator angle for each full arc was generated. The optimal segmentation of sub-arcs was achieved by avoiding the worst MCI at each control point. Then, the optimal collimator angle for each sub-arc would be determined by maximizing the summation of MCI. Each patient was scheduled to undergo single-center coplanar VMAT SRS based on either the novel SACAO algorithm or the conventional VMAT with static collimator angle (ST-VMAT). The dosimetric parameters, field sizes, and the monitoring units (Mus) were evaluated.

Results

The mean dose-volumetric parameters for the target volume of SACAO were comparable to ST-VMAT, while the conformity index (CI), homogeneity index (HI), and gradient index (GI) were reduced by SACAO. Improved sparing of organs at risk (OARs) was also obtained by SACAO. In particular, the SACAO method significantly (p < 0.01) reduced the field size (76.59 ± 32.55 vs. 131.95 ± 56.71 cm2) and MUs (655.35 ± 71.99 vs. 729.85 ± 73.52) by 41.11%.

Conclusions

The SACAO method could be superior in improving the CI, HI, and GI of the targets as well as normal tissue sparing for multiple brain metastases SRS. In particular, SACAO has the potential of increasing treatment efficiency in terms of field size and MU.

Comparison of anatomically informed class solution template trajectories with patient-specific trajectories for stereotactic radiosurgery and radiotherapy

Class solution template trajectories are used clinically for efficiency, safety, and reproducibility. The aim was to develop class solutions for single cranial metastases radiotherapy/radiosurgery based on intracranial target positioning and compare to patient-specific trajectories in the context of 4π optimization. Template trajectories were constructed based on the open-source Montreal Neurological Institute (MNI) average brain. The MNI brain was populated with evenly spaced spherical target volumes (2 cm diameter, N = 243) and organs-at-risk (OARs) were identified. Template trajectories were generated for six anatomical regions (frontal, medial, and posterior, each with laterality dependence) based on previously published 4π optimization methods. Volumetric modulated arc therapy (VMAT) treatment plans generated using anatomically informed template 4π trajectories and patientspecific 4π trajectories were compared against VMAT plans from a standard four-arc template. Four-arc optimization techniques were compared to the standard VMAT template by placing three spherical targets in each of six anatomical regions of a test patient. This yielded 54 plans to compare various plan quality metrics. Increasing plan technique complexity, the total number of OAR maximum dose reductions compared to the standard arc template for the 6 anatomical classes was 4+/-2 (OFIXEDc) and 7+/-2 (OFIXEDi). In 65.6% of all cases, optimized fixed-couch positions outperformed the standard-arc template. Of the three comparisons, the most complex (OFIXEDi) showed the greatest statistical significance compared to the least complex (VMATi) across 12 plan quality metrics of maximum dose to each OAR, V12Gy, total plan Monitor Units, conformity index, and gradient index (p < 0.00417). In approximately 70% of all cases, 4π optimization methods outperformed the standard-arc template in terms of maximum dose reduction to OAR, by exclusively changing the arc geometry. We conclude that a tradeoff exists between complexity of a class solution methodology compared to patient-specific methods for arc selection, in the context of plan quality improvement.

Organ-at-risk sparing with dynamic trajectory radiotherapy for head and neck cancer: comparison with volumetric arc therapy on a publicly available library of cases

Background

Dynamic trajectory radiotherapy (DTRT) extends volumetric modulated arc therapy (VMAT) with dynamic table and collimator rotation during beam-on. The aim of the study is to establish DTRT path-finding strategies, demonstrate deliverability and dosimetric accuracy and compare DTRT to state-of-the-art VMAT for common head and neck (HN) cancer cases.

Methods

A publicly available library of seven HN cases was created on an anthropomorphic phantom with all relevant organs-at-risk (OARs) delineated. DTRT plans were generated with beam incidences minimizing fractional target/OAR volume overlap and compared to VMAT. Deliverability and dosimetric validation was carried out on the phantom.

Results

DTRT and VMAT had similar target coverage. For three locoregionally advanced oropharyngeal carcinomas and one adenoid cystic carcinoma, mean dose to the contralateral salivary glands, pharynx and oral cavity was reduced by 2.5, 1.7 and 3.1 Gy respectively on average with DTRT compared to VMAT. For a locally recurrent nasopharyngeal carcinoma, D_0.03 cc to the ipsilateral optic nerve was above tolerance (54.0 Gy) for VMAT (54.8 Gy) but within tolerance for DTRT (53.3 Gy). For a laryngeal carcinoma, DTRT resulted in higher dose than VMAT to the pharynx and brachial plexus but lower dose to the upper oesophagus, thyroid gland and contralateral carotid artery. For a single vocal cord irradiation case, DTRT spared most OARs better than VMAT. All plans were delivered successfully on the phantom and dosimetric validation resulted in gamma passing rates of 93.9% and 95.8% (2%/2 mm criteria, 10% dose threshold).

Conclusions

This study provides a proof of principle of DTRT for common HN cases with plans that were deliverable on a C-arm linac with high accuracy. The comparison with VMAT indicates substantial OAR sparing could be achieved.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13014-022-02092-5.

Patient-Specific Collision Zones for 4π Trajectory Optimized Radiation Therapy

PURPOSE

The 4π methodology determines optimized non-coplanar sub arcs for stereotactic radiation therapy which minimize dose to organs-at-risk. Every combination of treatment angle is examined, but some angles are not appropriate as a collision would occur between the gantry and the couch or the gantry and the patient. Those combinations of couch and gantry angles are referred to as collision zones. A major barrier to applying 4π to stereotactic body radiation therapy (SBRT) is the unknown shape of the collision zones, which are significant as patients take up a large volume within the 4π sphere. This study presents a system which determines patient-specific collision zones, without additional clinical steps, to enable safe and deliverable non-coplanar treatment trajectories for SBRT patients.

METHODS

To augment patient's computed tomography (CT) scan, full body scans of patients in treatment position were acquired using an optical scanner. A library of a priori scans (N = 25) was created. Based on the patients treatment position and their body dimensions, a library scan is selected and registered to the CT scan of the patient. Next, a model of the couch and immobilization equipment is added to the patient model. This results in a patient model that is then aligned with a model of the treatment linac in a "virtual treatment room", where both components can be rotated to test for collisions. To test the collision detection algorithm, an end-to-end test was performed using a cranial phantom. The registration algorithm was tested by comparing the registered patient collision zones to those generated by using the patient's matching scan.

RESULTS

The collision detection algorithm was found to have a 97.80% accuracy, a 99.99% sensitivity and a 99.99% negative predictive value (NPV). Analysis of the registration algorithm determined that a 6 cm buffer was required to achieve a 99.65% mean sensitivity, where a sensitivity of unity is considered to be a requirement for safe treatment delivery. With a 6 cm buffer the mean accuracy was 86.70% and the mean NPV was 99.33%.

CONCLUSIONS

Our method of determining patient-specific collision zones can be accomplished with minimal user intervention based on an a priori library of body surface scans, thus enabling the safe application of 4π SBRT.

SBRT of ventricular tachycardia using 4pi optimized trajectories

PURPOSE

To investigate the possible advantages of using 4pi-optimized arc trajectories in stereotactic body radiation therapy of ventricular tachycardia (VT-SBRT) to minimize exposure of healthy tissues.

METHODS AND MATERIALS

Thorax computed tomography (CT) data for 15 patients were used for contouring organs at risk (OARs) and defining realistic planning target volumes (PTVs). A conventional trajectory plan, defined as two full coplanar arcs was compared to an optimized-trajectory plan provided by a 4pi algorithm that penalizes geometric overlap of PTV and OARs in the beam's-eye-view. A single fraction of 25 Gy was prescribed to the PTV in both plans and a comparison of dose sparing to OARs was performed based on comparisons of maximum, mean, and median dose.

RESULTS

A significant average reduction in maximum dose was observed for esophagus (18%), spinal cord (26%), and trachea (22%) when using 4pi-optimized trajectories. Mean doses were also found to decrease for esophagus (19%), spinal cord (33%), skin (18%), liver (59%), lungs (19%), trachea (43%), aorta (11%), inferior vena cava (25%), superior vena cava (33%), and pulmonary trunk (26%). A median dose reduction was observed for esophagus (40%), spinal cord (48%), skin (36%), liver (72%), lungs (41%), stomach (45%), trachea (53%), aorta (45%), superior vena cava (38%), pulmonary veins (32%), and pulmonary trunk (39%). No significant difference was observed for maximum dose (p = 0.650) and homogeneity index (p = 0.156) for the PTV. Average values of conformity number were 0.86 ± 0.05 and 0.77 ± 0.09 for the conventional and 4pi optimized plans respectively.

CONCLUSIONS

4pi optimized trajectories provided significant reduction to mean and median doses to cardiac structures close to the target but did not decrease maximum dose. Significant improvement in maximum, mean and median doses for noncardiac OARs makes 4pi optimized trajectories a suitable delivery technique for treating VT.

Comparison of non-coplanar optimization of static beams and arc trajectories for intensity-modulated treatments of meningioma cases

Two methods for non-coplanar beam direction optimization, one for static beams and another for arc trajectories, were proposed for intracranial tumours. The results of the beam angle optimizations were compared with the beam directions used in the clinical plans. Ten meningioma cases already treated were selected for this retrospective planning study. Algorithms for non-coplanar beam angle optimization (BAO) and arc trajectory optimization (ATO) were used to generate the corresponding plans. A plan quality score, calculated by a graphical method for plan assessment and comparison, was used to guide the beam angle optimization process. For each patient, the clinical plans (CLIN), created with the static beam orientations used for treatment, and coplanar VMAT approximated plans (VMAT) were also generated. To make fair plan comparisons, all plan optimizations were performed in an automated multicriteria calculation engine and the dosimetric plan quality was assessed. BAO and ATO plans presented, on average, moderate global plan score improvements over VMAT and CLIN plans. Nevertheless, while BAO and CLIN plans assured a more efficient OARs sparing, the ATO and VMAT plans presented a higher coverage and conformity of the PTV. Globally, all plans presented high-quality dose distributions. No statistically significant quality differences were found, on average, between BAO, ATO and CLIN plans. However, automated plan solution optimizations (BAO or ATO) may improve plan generation efficiency and standardization. In some individual patients, plan quality improvements were achieved with ATO plans, demonstrating the possible benefits of this automated optimized delivery technique.

Fully automated noncoplanar radiation therapy treatment planning

PURPOSE

To perform fully automated noncoplanar (NC) treatment planning, we propose a method called NC-POPS to produce NC plans using the Pareto optimal projection search (POPS) algorithm.

METHODS

NC radiation therapy treatment planning has the potential to improve dosimetric quality as compared to traditional coplanar techniques. Likewise, automated treatment planning algorithms can reduce a planner's active treatment planning time and remove inter-planner variability. Our NC-POPS algorithm extends the original POPS algorithm to the NC setting with potential applications to both intensity-modulated radiation therapy (IMRT) and volumetric modulated arc therapy (VMAT). The proposed algorithm consists of two main parts: (1) NC beam angle optimization (BAO) and (2) fully automated inverse planning using the POPS algorithm.

RESULTS

We evaluate the performance of NC-POPS by comparing between various NC and coplanar configurations. To evaluate plan quality, we compute the homogeneity index (HI), conformity index (CI), and dose-volume histogram statistics for various organs-at-risk (OARs). As compared to the evaluated coplanar baseline methods, the proposed NC-POPS method achieves significantly better OAR sparing, comparable or better dose conformity, and similar dose homogeneity.

CONCLUSIONS

Our proposed NC-POPS algorithm provides a modular approach for fully automated treatment planning of NC IMRT cases with the potential to substantially improve treatment planning workflow and plan quality.

Verification of the delivered patient radiation dose for non-coplanar beam therapy

PURPOSE

There is an increased interest in using non-coplanar beams for radiotherapy, including SBRT and SRS. This approach can significantly reduce doses to organs-at-risk, however, it requires stringent quality assurance, especially when a dynamic treatment couch is used. In this work, new functionality that allows using non-coplanar beam arrangements in addition to conventional coplanar beams was added and validated to the previously developed in vivo dose verification system.

METHODS

The existing program code was modified to manage the additional treatment couch parameters: angle and positions. Ten non-coplanar test plans that use a static couch were created in the treatment planning system. Also, two plans that use a dynamic treatment couch were created and delivered using Varian Developer mode, since the treatment planning system does not support a dynamic couch. All non-coplanar test trajectories were delivered on a simple geometric phantom, using an Edge linear accelerator (Varian Medical Systems) with the megavoltage imager deployed and acquiring megavoltage transmission images that were used to calculate the delivered 3D dose distributions in the phantom with the updated dose calculation algorithm. The reconstructed dose distributions were compared using the 3D chi-comparison test with 2%/2mm tolerances to the corresponding reference dose distributions obtained from the treatment planning system.

RESULTS

The chi-comparison test resulted in at least a 97.0% pass rate over the entire 3D volume for all tested trajectories. For static gantry, static couch non-coplanar fields, and non-coplanar arcs using dynamic couch the pass rates observed were at least 98%, while for the static couch, non-transverse coplanar arc fields, pass rates were at least 97%.

CONCLUSIONS

A model-based 3D dose calculation algorithm has been extended and validated for a variety of non-coplanar beam trajectories of different complexities. This system can potentially be applied for quality assurance of treatment delivery systems that use complex, non-coplanar beam arrangements.

Dosimetric comparison of coplanar and noncoplanar beam arrangements for radiotherapy of patients with lung cancer: A meta-analysis

Purpose

Radiotherapy plays an important role in the treatment of lung cancer, and both coplanar beam arrangements (CBA) and noncoplanar beam arrangements (NCBA) are adopted in clinic practice. The aim of this study is to answer the question whether NCBA are dosimetrically superior to CBA.

Methods

Search of publications were performed in PubMed, Web of Science, and the Cochran Library till March 2020. The searching terms were as following: ((noncoplanar) or ("non coplanar") or ("4pi") or ("4π")) AND (("lung cancer") or ("lung tumor") or ("lung carcinoma")) AND ((radiotherapy) or ("radiation therapy")). The included studies and extracted data were manually screened. All forest and funnel plots were carried out with RevMan software, and the Egger’s regression asymmetry tests were conducted with STATA software.

Results

Nine studies were included and evaluated in the meta‐analysis and treatment plans were designed with both CBA and NCBA. For the planning target volumes (PTV), D98%, D2%, the conformity index (CI), and the gradient index (GI) had no statistically significant difference. For organs‐at‐risk (OAR), V20 of the whole lung and the maximum dose of the spinal cord were significantly reduced in NCBA plans compared with CBA ones. But V10, V5, and mean dose of the whole lung, the maximum dose of the heart, and the maximum dose of the esophagus exhibited no significant difference when the two types of beam arrangements were compared.

Conclusion

After combining multicenter results, NCBA plans have significant advantages in reducing V20 of the whole lung and max dose of spinal cord.

Radiobiological effects of the interruption time with Monte Carlo Simulation on multiple fields in photon beams

PURPOSE

The interruption time is the irradiation interruption that occurs at sites and operations such as the gantry, collimator, couch rotation, and patient setup within the field in radiotherapy. However, the radiobiological effect of prolonging the treatment time by the interruption time for tumor cells is little evaluated. We investigated the effect of the interruption time on the radiobiological effectiveness with photon beams based on a modified microdosimetric kinetic (mMK) model.

METHODS

The dose-mean lineal energy y_D (keV/µm) of 6-MV photon beams was calculated by the particle and heavy ion transport system (PHITS). We set the absorbed dose to 2 or 8 Gy, and the interruption time (τ) was set to 1, 3, 5, 10, 30, and 60 min. The biological parameters such as α_0, β_0, and DNA repair constant rate (a + c) values were acquired from a human non-small-cell lung cancer cell line (NCI-H460) for the mMK model. We used two-field and four-field irradiation with a constant dose rate (3 Gy/min); the photon beams were paused for interruption time τ. We calculated the relative biological effectiveness (RBE) to evaluate the interruption time's effect compared with no interrupted as a reference.

RESULTS

The y_D of 6-MV photon beams was 2.32 (keV/µm), and there was little effect by changing the water depth (standard deviation was 0.01). The RBE with four-field irradiation for 8 Gy was decreased to 0.997, 0.975, 0.900, and 0.836 τ = 1, 10, 30, 60 min, respectively. In addition, the RBE was affected by the repair constant rate (a + c) value, the greater the decrease in RBE with the longer the interruption time when the (a + c) value was large.

CONCLUSION

The ~10-min interruption of 6-MV photon beams did not significantly impact the radiobiological effectiveness, since the RBE decrease was <3%. Nevertheless, the RBE's effect on tumor cells was decreased about 30% by increasing the 60 min interruption time at 8 Gy with four-field irradiation. It is thus necessary to make the interruption time as short as possible.

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.

Many-isocenter optimization for robotic radiotherapy

Despite significant dosimetric gains, clinical implementation of the 4π non-coplanar radiotherapy on the widely available C-arm gantry system is hindered by limited clearance, and the need to perform complex coordinated gantry and couch motion. A robotic radiotherapy platform would be conducive to such treatment but a new conflict between field size and MLC modulation resolution needs to be managed for versatile applications. This study investigates the dosimetry and delivery efficiency of purposefully creating many isocenters to achieve simultaneously high MLC modulation resolution and large tumor coverage. An integrated optimization framework was proposed for simultaneous beam orientation optimization (BOO), isocenter selection, and fluence map optimization (FMO). The framework includes a least-square dose fidelity objective, a total variation term for regularizing the fluence smoothness, and a group sparsity term for beam selection. A minimal number of isocenters were identified for efficient target coverage. Colliding beams excluded, high-resolution small-field 4π intensity-modulated radiotherapy (IMRT) treatment plans with 50 cm source-to-isocenter distance (SID-50) on 10 Head and Neck (H&N) cancer patients were compared with low-resolution large-field plans with 100 cm SID (SID-100). With the same or better target coverage, the average reduction of [Dmean, Dmax] of 20-beam SID-50 plans from 20-beam SID-100 plans were [2.09 Gy, 1.19 Gy] for organs at risk (OARs) overall, [3.05 Gy, 0.04 Gy] for parotid gland, [3.62 Gy, 5.19 Gy] for larynx, and [3.27 Gy, 1.10 Gy] for mandible. R50 and integral dose were reduced by 5.3% and 9.6%, respectively. Wilcoxon signed-rank test showed significant difference (p   <  0.05) in planning target volume (PTV) homogeneity, PTV Dmax, R50, Integral dose, and OAR Dmean and Dmax. The estimated delivery time of 20-beam [SID-50, SID-100] plans were [19, 18] min and [14, 9] min, assuming 5 fractions and 30 fractions, respectively. With clinically acceptable delivery efficiency, many-isocenter optimization is dosimetrically desirable for treating large targets with high modulation resolution on the robotic platform.

Toward the combined optimization of dynamic axes (CODA) for stereotactic radiotherapy and radiosurgery using fixed couch trajectories

PURPOSE

To develop a novel system for patient-specific combined optimization of couch, collimator, and gantry angles for use in volumetric modulated arc therapy (VMAT) treatment planning. The system was designed to produce highly compact dose distributions by extensively sampling the 4π space. Automated fixed couch trajectory planning was used to reduce normal tissue doses by avoiding beams-eye-view (BEV) overlap with organs-at-risk (OARs) and improve monitor unit (MU) efficiency through collimator angle optimization.

METHODS

By merging distinct BEV objective functions used to optimize the couch rotation angle and collimator angle, a three-dimensional (3D) cost space (the CODA cube) was constructed with axes of gantry, couch, and collimator rotation angles. At each voxel in this CODA cube, the cost of implementing this combination of axes positions in fixed couch trajectories was quantified. The CODA cube was sampled and explored using a modified constrained Bellman-Ford algorithm to suggest low-cost fixed candidate arcs on each plane of the space, from which 10-arcs are chosen throughout the 3D space using a k-means clustering algorithm. These fixed couch trajectories were then imported into the Eclipse treatment planning system (v.11) and inverse-optimized according to clinical standards. Eight artificial cranial targets were contoured in a test-patient anatomy, and seven treatment plans were generated from combinations of three and four targets. The CODA cube optimized plans were compared to standard 4-arc VMAT plans for cranial stereotactic radiotherapy/surgery that were optimized for the same sets of targets; maximum dose to each OAR, V12Gy to normal brain, conformity, and total MUs were compared. Both planning methods were inverse-optimized with identical dosimetric objectives.

RESULTS

CODA plans resulted in a reduction in maximum dose to OARs of 20.6% (P < 0.01), with maximum brainstem dose decreased by 2.63 Gy (P = 0.031) on average when compared to the standard arc arrangement. The mean reduction in total MU was 8.6% (P = 0.156), the mean increase in the inverse of the van't Riet conformation number was 0.1%, (P = 0.67) and the mean decrease in normal brain tissue receiving 12 Gy or higher was 3.9% (P = 0.16), when compared to the standard VMAT arc configuration (n = 7).

CONCLUSIONS

The optimization of couch, collimator, and gantry angles simultaneously using a 3D optimization space achieved improvement on multiple clinical metrics when compared to the standard VMAT arc configuration. A statistically significant sparing to OAR maximum doses was seen. Combining these optimizations may yield superior results to independent optimization.

Treatment planning optimization with beam motion modeling for dynamic arc delivery of SBRT using Cyberknife with multileaf collimation

PURPOSE

The use of dynamic arcs for delivery of stereotactic body radiation therapy (SBRT) on Cyberknife is investigated, with a view to improving treatment times. This study investigates the required modeling of robot and multileaf collimator (MLC) motion between control points in the trajectory and then uses this to develop an optimization method for treatment planning of a dynamic arc with Cyberknife. The resulting plans are compared in terms of dose-volume histograms and estimated treatment times with those produced by a conventional beam arrangement.

METHODS

Five SBRT patient cases (prostate A - conventional, prostate B - brachytherapy-type, lung, liver, and partial left breast) were retrospectively studied. A suitable arc trajectory with control points spaced at 5° was proposed and treatment plans were produced for typical clinical protocols. The optimization consisted of a fluence optimization, segmentation, and direct aperture optimization using a gradient descent method. Dose delivered by the moving MLC was either taken to be the dose delivered discretely at the control points or modeled using effective fluence delivered between control points. The accuracy of calculated dose was assessed by recalculating after optimization using five interpolated beams and 100 interpolated apertures between each optimization control point. The resulting plans were compared using dose-volume histograms and estimated treatment times with those for a conventional Cyberknife beam arrangement.

RESULTS

If optimization is performed based on discrete doses delivered at the arc control points, large differences of up to 40% of the prescribed dose are seen when recalculating with interpolation. When the effective fluence between control points is taken into account during optimization, dosimetric differences are <2% for most structures when the plans are recalculated using intermediate nodes, but there are differences of up to 15% peripherally. Treatment plan quality is comparable between the arc trajectory and conventional body path. All plans meet the relevant clinical goals, with the exception of specific structures which overlap with the planning target volume. Median estimated treatment time is 355 s (range 235-672 s) for arc delivery and 675 s (range 554-1025 s) for conventional delivery.

CONCLUSIONS

The method of using effective fluence to model MLC motion between control points is sufficiently accurate to provide for accurate inverse planning of dynamic arcs with Cyberknife. The proposed arcing method produces treatment plans with comparable quality to the body path, with reduced estimated treatment delivery time.

Parallel beamlet dose calculation via beamlet contexts in a distributed multi-GPU framework

PURPOSE

Dose calculation is one of the most computationally intensive, yet essential tasks in the treatment planning process. With the recent interest in automatic beam orientation and arc trajectory optimization techniques, there is a great need for more efficient model-based dose calculation algorithms that can accommodate hundreds to thousands of beam candidates at once. Foundational work has shown the translation of dose calculation algorithms to graphical processing units (GPUs), lending to remarkable gains in processing efficiency. But these methods provide parallelization of dose for only a single beamlet, serializing the calculation of multiple beamlets and under-utilizing the potential of modern GPUs. In this paper, the authors propose a framework enabling parallel computation of many beamlet doses using a novel beamlet context transformation and further embed this approach in a scalable network of multi-GPU computational nodes.

METHODS

The proposed context-based transformation separates beamlet-local density and TERMA into distinct beamlet contexts that independently provide sufficient data for beamlet dose calculation. Beamlet contexts are arranged in a composite context array with dosimetric isolation, and the context array is subjected to a GPU collapsed-cone convolution superposition procedure, producing the set of beamlet-specific dose distributions in a single pass. Dose from each context is converted to a sparse representation for efficient storage and retrieval during treatment plan optimization. The context radius is a new parameter permitting flexibility between the speed and fidelity of the dose calculation process. A distributed manager-worker architecture is constructed around the context-based GPU dose calculation approach supporting an arbitrary number of worker nodes and resident GPUs. Phantom experiments were executed to verify the accuracy of the context-based approach compared to Monte Carlo and a reference CPU-CCCS implementation for single beamlets and broad beams composed by addition of beamlets. Dose for representative 4π beam sets was calculated in lung and prostate cases to compare its efficiency with that of an existing beamlet-sequential GPU-CCCS implementation. Code profiling was also performed to evaluate the scalability of the framework across many networked GPUs.

RESULTS

The dosimetric accuracy of the context-based method displays <1.35% and 2.35% average error from the existing serialized CPU-CCCS algorithm and Monte Carlo simulation for beamlet-specific PDDs in water and slab phantoms, respectively. The context-based method demonstrates substantial speedup of up to two orders of magnitude over the beamlet-sequential GPU-CCCS method in the tested configurations. The context-based framework demonstrates near linear scaling in the number of distributed compute nodes and GPUs employed, indicating that it is flexible enough to meet the performance requirements of most users by simply increasing the hardware utilization.

CONCLUSIONS

The context-based approach demonstrates a new expectation of performance for beamlet-based dose calculation methods. This approach has been successful in accelerating the dose calculation process for very large-scale treatment planning problems - such as automatic 4π IMRT beam orientation and VMAT arc trajectory selection, with hundreds of thousands of beamlets - in clinically feasible timeframes. The flexibility of this framework makes it as a strong candidate for use in a variety of other very large-scale treatment planning tasks and clinical workflows.

Recent developments in non-coplanar radiotherapy

This paper gives an overview of recent developments in non-coplanar intensity modulated radiotherapy (IMRT) and volumetric modulated arc therapy (VMAT). Modern linear accelerators are capable of automating motion around multiple axes, allowing efficient delivery of highly non-coplanar radiotherapy techniques. Novel techniques developed for C-arm and non-standard linac geometries, methods of optimization, and clinical applications are reviewed. The additional degrees of freedom are shown to increase the therapeutic ratio, either through dose escalation to the target or dose reduction to functionally important organs at risk, by multiple research groups. Although significant work is still needed to translate these new non-coplanar radiotherapy techniques into the clinic, clinical implementation should be prioritized. Recent developments in non-coplanar radiotherapy demonstrate that it continues to have a place in modern cancer treatment.

Dosimetric accuracy of dynamic couch rotation during volumetric modulated arc therapy (DCR-VMAT) for primary brain tumours

Radiotherapy treatment plans using dynamic couch rotation during volumetric modulated arc therapy (DCR-VMAT) reduce the dose to organs at risk (OARs) compared to coplanar VMAT, while maintaining the dose to the planning target volume (PTV). This paper seeks to validate this finding with measurements. DCR-VMAT treatment plans were produced for five patients with primary brain tumours and delivered using a commercial linear accelerator (linac). Dosimetric accuracy was assessed using point dose and radiochromic film measurements. Linac-recorded mechanical errors were assessed by extracting deviations from log files for multi-leaf collimator (MLC), couch, and gantry positions every 20 ms. Dose distributions, reconstructed from every fifth log file sample, were calculated and used to determine deviations from the treatment plans. Median (range) treatment delivery times were 125 s (123-133 s) for DCR-VMAT, compared to 78 s (64-130 s) for coplanar VMAT. Absolute point doses were 0.8% (0.6%-1.7%) higher than prediction. For coronal and sagittal films, respectively, 99.2% (96.7%-100%) and 98.1% (92.9%-99.0%) of pixels above a 20% low dose threshold reported gamma  <1 for 3% and 3 mm criteria. Log file analysis showed similar gantry rotation root-mean-square error (RMSE) for VMAT and DCR-VMAT. Couch rotation RMSE for DCR-VMAT was 0.091° (0.086-0.102°). For delivered dose reconstructions, 100% of pixels above a 5% low dose threshold reported gamma  <1 for 2% and 2 mm criteria in all cases. DCR-VMAT, for the primary brain tumour cases studied, can be delivered accurately using a commercial linac.

Beam selection for stereotactic ablative radiotherapy using Cyberknife with multileaf collimation

The Cyberknife system (Accuray Inc., Sunnyvale, CA) enables radiotherapy using stereotactic ablative body radiotherapy (SABR) with a large number of non-coplanar beam orientations. Recently, a multileaf collimator has also been available to allow flexibility in field shaping. This work aims to evaluate the quality of treatment plans obtainable with the multileaf collimator. Specifically, the aim is to find a subset of beam orientations from a predetermined set of candidate directions, such that the treatment quality is maintained but the treatment time is reduced. An evolutionary algorithm is used to successively refine a randomly selected starting set of beam orientations. By using an efficient computational framework, clinically useful solutions can be found in several hours. It is found that 15 beam orientations are able to provide treatment quality which approaches that of the candidate beam set of 110 beam orientations, but with approximately half of the estimated treatment time. Choice of an efficient subset of beam orientations offers the possibility to improve the patient experience and maximise the number of patients treated.

Collision-Free Path Planning and Delivery Sequence Optimization in Noncoplanar Radiation Therapy

Radiation therapy is among the top three cancer treatments in current medical services. The novel noncoplanar radiation therapy which claimed the best characteristics in almost all dosimetric properties encountered the challenges of the potential collision and the long time delivering. In this paper, we proposed a brand new scheme which uses a combined method of the collision avoidance path planning based on an improved probability roadmap method (PRM) and the delivery sequence optimization based on a modified genetic algorithm (GA) to solve the problems in noncoplanar radiation therapy. A uniform sampling strategy, an improved connection strategy, and an efficient local planner are introduced to optimize the roadmap result and accelerate the roadmap construction. The GA is improved by the elitist selection, the local search strategy, and the similar substitution strategy to achieve a better performance both in convergence rate and optimal solution. Experiments are carried out on the simulation platform with typical therapy system models. The results show that our proposed methods work well with the radiation therapy system in a compact working area. Collision is avoided and time consumption is reduced. We believe that our proposed algorithms could solve the problems in current radiation therapy and promote their clinic applications.

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.

Part 1: Optimization and evaluation of dynamic trajectory radiotherapy

PURPOSE

Although volumetric modulated arc therapy (VMAT) is a well-accepted treatment technique in radiotherapy using a coplanar delivery approach, VMAT might be further improved by including dynamic table and collimator rotations leading to dynamic trajectory radiotherapy (DTRT). In this work, an optimization procedure for DTRT was developed and the potential benefit of DTRT was investigated for different treatment sites.

METHODS

For this purpose, a dedicated optimization framework for DTRT was developed using the Eclipse Scripting Research Application Programming Interface (ESRAPI). The contours of the target and organs at risk (OARs) structures were exported by applying the ESRAPI and were used to determine the fractional volume-overlap of the OARs with the target from several potential beam directions. Thereby, an additional weighting was applied taking into account the relative position of the OAR with respect to the target and radiation beam, that is, penalizing directions where the OAR is proximal to the target. The resulting two-dimensional gantry-table map was used as input for an A* path finding algorithm returning an optimized gantry-table path. Thereby, the process is also taking into account CT scan length and collision restrictions. The A* algorithm was used again to determine the dynamic collimator angle path by optimizing the area between the MLC leaves and the target contour for each gantry-table path leading to gantry-collimator paths. The resulting gantry-table and gantry-collimator paths are combined and serve as input for the intensity modulation optimization using a research VMAT optimizer and the ESRAPI resulting in dynamic trajectories. This procedure was evaluated for five clinically motivated cases: two head and neck, one lung, one esophagus, and one prostate. Final dose calculations were performed using the Swiss Monte Carlo Plan (SMCP). Resulting dose distributions for the DTRT treatment plans and for the standard VMAT plans were compared based on dose distributions and dose volume histogram (DVH) parameters. For this comparison, the dose distribution for the VMAT plans were recalculated using the SMCP. In addition, the suitability of the delivery of a DTRT treatment plan was demonstrated by means of gafchromic film measurements on a TrueBeam linear accelerator.

RESULTS

DVHs for the target volumes showed similar or improved coverage and dose homogeneity for DTRT compared with VMAT using equal or less number of dynamic trajectories for DTRT than arcs for VMAT for all cases studied. Depending on the case, improvements in mean and maximum dose for the DTRT plans were achieved for almost all OARs compared with the VMAT plans. Improvements in DTRT treatment plans for mean and maximum dose compared to VMAT plans were up to 16% and 38% relative to the prescribed dose, respectively. The measured and calculated dose values resulted in a passing rate of more than 99.5% for the two-dimensional gamma analysis using 2% and 2 mm criteria and a threshold of 10%.

CONCLUSIONS

DTRT plans for different treatment sites were generated and compared with VMAT plans. The delivery is suitable and dose comparisons demonstrate a high potential of DTRT to reduce dose to OARs using less dynamic trajectories than arcs, while target coverage is preserved.

A novel optimization framework for VMAT with dynamic gantry couch rotation

Existing volumetric modulated arc therapy (VMAT) optimization using coplanar arcs is highly efficient but usually dosimetrically inferior to intensity modulated radiation therapy (IMRT) with optimized non-coplanar beams. To achieve both dosimetric quality and delivery efficiency, we proposed in this study, a novel integrated optimization method for non-coplanar VMAT (4πVMAT). 4πVMAT with direct aperture optimization (DAO) was achieved by utilizing a least square dose fidelity objective, along with an anisotropic total variation term for regularizing the fluence smoothness, a single segment term for imposing simple apertures, and a group sparsity term for selecting beam angles. Continuous gantry/couch angle trajectories were selected using the Dijkstra's algorithm, where the edge and node costs were determined based on the maximal gantry rotation speed and the estimated fluence map at the current iteration, respectively. The couch-gantry-patient collision space was calculated based on actual machine geometry and a human subject 3D surface. Beams leading to collision are excluded from the DAO and beam trajectory selection (BTS). An alternating optimization strategy was implemented to solve the integrated DAO and BTS problem. The feasibility of 4πVMAT using one full-arc or two full-arcs was tested on nine patients with brain, lung, or prostate cancer. The plan was compared against a coplanar VMAT (2πVMAT) plan using one additional arc and collimator rotation. Compared to 2πVMAT, 4πVMAT reduced the average maximum and mean organs-at-risk dose by 9.63% and 3.08% of the prescription dose with the same target coverage. R50 was reduced by 23.0%. Maximum doses to the dose limiting organs, such as the brainstem, the major vessels, and the proximal bronchus, were reduced by 8.1 Gy (64.8%), 16.3 Gy (41.5%), and 19.83 Gy (55.5%), respectively. The novel 4πVMAT approach affords efficient delivery of non-coplanar arc trajectories that lead to dosimetric improvements compared with coplanar VMAT using more arcs.

Machine-specific quality assurance procedure for stereotactic treatments with dynamic couch rotations

OBJECTIVE

We present a method in which the treatment couch's accuracy is measured using the electronic portal imaging device (EPID) and a phantom of our own construction. Using this phantom, we were able to quantify the treatment couch walkout, and the rotation angle accuracy for both static and dynamic couch treatments. These measurements were used to provide an accurate measure of the treatment couch isocenter as well as to verify the couch rotation angle recorded in the trajectory log.

METHODS

The phantom was constructed using a polystyrene slab in which five ball bearings of 4 mm diameter are placed on the same plane at varying radii (0, 2.8, 4.4, 5.6, and 6.7 cm). The couch was rotated through its full extent (-90, 90 degrees) while MV images were acquired continuously. The couch rotational accuracy was calculated using a least squares minimization which fit the locations of the BBs to their expected locations relative to reference setup conditions. Using this approach, rotation angle and isocenter walkout was calculated in three dimensions. These measurements were used to quantify the accuracy of the couch as well as to validate the Varian TrueBeam trajectory logs. Additionally, a method for an EPID-based couch star-shot measurement was developed and compared with the traditional film-based method.

RESULTS

The measured couch center of rotation consisted of a cloud of points clustered around the room isocenter within 0.7 mm distance. The trajectory log couch angle values agreed with those recorded in the DICOM header of the EPID images to the third significant digit and the couch rotation angles recorded in the trajectory log and DICOM header agreed with the calculated values to 0.08 degrees. Comparison of couch star-shot measurement developed in this study with film-based star-shot measurements gave an agreement to within 0.2 mm.

CONCLUSION

We have developed a quality assurance method for the treatment couch which is simple, accurate, and enables the user to access a multitude of consistent data with a single measurement. Using this method, we have shown that the treatment couch is accurate for both static and dynamic stereotactic deliveries.

Trajectory optimization in radiotherapy using sectioning (TORUS)

PURPOSE

A challenging problem in trajectory optimization for radiotherapy is properly handling the synchronization of the medical accelerators dynamic delivery. The initial coarse sampling of control points implemented in a Progressive Resolution Optimization type approach (VMAT) routinely results in MLC aperture forming contention issues as the sampling resolution increases. This work presents an approach to optimize continuous, beam-on radiation trajectories through exploration of the anatomical topology present in the patient and formation of a novel dual-metric graph optimization problem.

METHODS

This work presents a new perspective on trajectory optimization in radiotherapy using the concept of sectioning (TORUS). TORUS avoids degradation of 3D dose optimization quality by mapping the connectedness of target regions from the BEV perspective throughout the space of deliverable coordinates. This connectedness information is then incorporated into a graph optimization problem to define ideal trajectories. The unique usage of two distance functions in this graph optimization permits the TORUS algorithm to generate efficient dynamic trajectories for delivery while maximizing the angular flux through all PTV voxels. 3D dose optimization is performed for trajectories using a commercial TPS progressive resolution optimizer.

RESULTS

The TORUS algorithm is applied to three example treatments: chest wall, scalp, and the TG-119 C-shape phantom. When static collimator coplanar trajectories are generated for the chest wall and scalp cases, the TORUS trajectories are found to outperform both 7 field IMRT and 2 arc VMAT plans in delivery time, organ at risk sparing, conformality, and homogeneity. For the TG-119 phantom, when static couch and collimator non-coplanar trajectories are optimized, TORUS trajectories have superior sparing of the central core avoidance with shorter delivery times, with similar dose conformality and homogeneity.

CONCLUSIONS

The TORUS algorithm is able to automatically generate trajectories having improved plan quality and delivery time over standard IMRT and VMAT treatments. TORUS offers an exciting and promising avenue forward toward increasing dynamic capabilities in radiation delivery.

Trajectory modulated arc therapy using quasi-continuous couch motion layered on top of volumetric modulated arc therapy in left breast and chest wall irradiation: a feasibility study

Aim

To investigate the dosimetric advantage of quasi-continuous couch motion-enabled trajectory modulated arc radiotherapy therapy (TMAT) over the coplanar tangential partial arcs volumetric modulated arc radiotherapy (VMAT) for treating left breast and chest wall patients.

Method

Treatment plans of 43 patients who received radiotherapy for left breast (17) or for left chest wall (26) using coplanar partial tangential arcs VMAT (reference plan) were considered for this study. For each patient, in addition to the treatment plan, a TMAT plan was also generated using quasi-continuous couch rotation. The TMAT plan consisted of original two 30° tangential arc beams and two supplementary beams having a couch rotation of ±10°, ±20° and ±30°, respectively. The difference in PTV volume coverage (PTV V95%) between TMAT plan and VMAT plan was calculated for all the cases and normalised to the plan’s prescription dose. Similarly, differences in PTV_V105% and several dose-volume parameters related to organs at risk (OAR) were also computed and tabulated.

Result

TMAT shows an increment in the PTV dose coverage V95% with respect to reference plan by 4·7±2·5% when averaged overall prescription dose levels. Mean PTV dose (averaged overall prescription levels) for reference and TMAT plan was 4638·6±423·8 and 4793·5±447·2 cGy, respectively, and statistically insignificant (p=0·06). However mean PTV_V105% values for TMAT and for reference plans were 6·7±4·8 and 7·2±5·2%, respectively, and were not statistically different (p=0·85). Mean heart dose in TMAT was less than in VMAT plans, but not significantly. As regarding D1% to heart, TMAT plan was again found to be better with a mean difference of 137·1 cGy over VMAT plan. Other parameters evaluated were: mean dose and D1% to contralateral breast, and V20 Gy and V5 Gy for lung.

Conclusion

TMAT plans were found to be better than VMAT plans in terms of PTV coverage and D1% for heart. For evaluated dose parameters apart from PTV coverage and D1% to the heart, no significant differences were observed. Thus, TMAT plans yielded better dose distribution in terms of PTV dose coverage, hot spots and OAR doses.

A simple and robust trajectory-based stereotactic radiosurgery treatment

INTRODUCTION

We present a Trajectory-based Volumetric Modulated Arc Therapy (TVMAT) technique for Stereotactic Radiosurgery (SRS) that takes advantage of a modern linacs ability to modulate dose rate and move the couch dynamically. In addition, we investigate the quality of the developed TVMAT method and the dosimetric accuracy of the technique.

METHODS

The main feature of the TVMAT technique is a standard beam trajectory formed by dynamic motion of the treatment couch and the linac gantry. The couch rotates slowly through 180 degrees while the gantry delivers radiation through continuous sweeps of the gantry. The number of partial arcs that constitute the trajectory can be varied between two and eight and as the number of partial arcs increases, the trajectory more finely samples 4π geometry. Along these trajectories, the multi-leaf collimator (MLC) and dose rate are optimized through an inverse planning framework. The TVMAT method was tested on ten cranial SRS patients who were previously treated with the Dynamic Conformal Arc (DCA) technique. The plans were compared with the DCA and a four- arc VMAT technique with regards to dose to the OAR, dose falloff, V12Gy, and V4Gy. Validation measurements were performed using ion-chamber and Gafchromic film. In addition, the trajectory-log files were analyzed and compared with the treatment plan beam data.

RESULTS

The TVMAT treatment plans were successfully delivered with a treatment time between 3-8 min which mostly depended on total cumulated dose. Ion chamber measurements had an average measured error of 1.1 ± 0.6% and a maximum value of 2.2% of the delivered dose. The 2%, 2 mm gamma pass rates for the film measurements were 96% or greater. In a preliminary comparison of ten patients who underwent SRS treatments with the DCA technique, the TVMAT and VMAT techniques were able to produce plans with comparable dose falloff and OAR doses, while achieving better dose conformality, V4Gy and V12Gy when compared to the original DCA plans. The improvement of the TVMAT plans were as follows (mean % improvement ± standard err): Conformity (10 ± 2%), V4 (20 ± 20%), V12 (27 ± 10%), volume weighted mean dose to organs at risk (13 ± 13%), homogeneity index (2 ± 2%) and falloff (4 ± 2%).

CONCLUSION

We have developed and validated a trajectory-based dose delivery method which has dose distribution improvements while having a treatment time of 3-8 min. In addition, it has the potential for a simpler planning experience while maintaining an accurate delivery on the Varian Truebeam Linac.

Non-coplanar trajectories to improve organ at risk sparing in volumetric modulated arc therapy for primary brain tumors

BACKGROUND AND PURPOSE

To evaluate non-coplanar volumetric modulated arc radiotherapy (VMAT) trajectories for organ at risk (OAR) sparing in primary brain tumor radiotherapy.

MATERIALS AND METHODS

Fifteen patients were planned using coplanar VMAT and compared against non-coplanar VMAT plans for three trajectory optimization techniques. A geometric heuristic technique (GH) combined beam scoring and Dijkstra's algorithm to minimize the importance-weighted sum of OAR volumes irradiated. Fluence optimization was used to perform a local search around coplanar and GH trajectories, producing fluence-based local search (FBLS) and FBLS+GH trajectories respectively.

RESULTS

GH, FBLS, and FBLS+GH trajectories reduced doses to the contralateral globe, optic nerve, hippocampus, temporal lobe, and cochlea. However, FBLS increased dose to the ipsilateral lens, optic nerve and globe. Compared to GH, FBLS+GH increased dose to the ipsilateral temporal lobe and hippocampus, contralateral optics, and the brainstem and body. GH and FBLS+GH trajectories reduced bilateral hippocampi normal tissue complication probability (p=0.028 and p=0.043, respectively). All techniques reduced PTV conformity; GH and FBLS+GH trajectories reduced homogeneity but less so for FBLS+GH.

CONCLUSIONS

The geometric heuristic technique best spared OARs and reduced normal tissue complication probability, however incorporating fluence information into non-coplanar trajectory optimization maintained PTV homogeneity.

The development and verification of a highly accurate collision prediction model for automated noncoplanar plan delivery

PURPOSE

Significant dosimetric benefits had been previously demonstrated in highly noncoplanar treatment plans. In this study, the authors developed and verified an individualized collision model for the purpose of delivering highly noncoplanar radiotherapy and tested the feasibility of total delivery automation with Varian TrueBeam developer mode.

METHODS

A hand-held 3D scanner was used to capture the surfaces of an anthropomorphic phantom and a human subject, which were positioned with a computer-aided design model of a TrueBeam machine to create a detailed virtual geometrical collision model. The collision model included gantry, collimator, and couch motion degrees of freedom. The accuracy of the 3D scanner was validated by scanning a rigid cubical phantom with known dimensions. The collision model was then validated by generating 300 linear accelerator orientations corresponding to 300 gantry-to-couch and gantry-to-phantom distances, and comparing the corresponding distance measurements to their corresponding models. The linear accelerator orientations reflected uniformly sampled noncoplanar beam angles to the head, lung, and prostate. The distance discrepancies between measurements on the physical and virtual systems were used to estimate treatment-site-specific safety buffer distances with 0.1%, 0.01%, and 0.001% probability of collision between the gantry and couch or phantom. Plans containing 20 noncoplanar beams to the brain, lung, and prostate optimized via an in-house noncoplanar radiotherapy platform were converted into XML script for automated delivery and the entire delivery was recorded and timed to demonstrate the feasibility of automated delivery.

RESULTS

The 3D scanner measured the dimension of the 14 cm cubic phantom within 0.5 mm. The maximal absolute discrepancy between machine and model measurements for gantry-to-couch and gantry-to-phantom was 0.95 and 2.97 cm, respectively. The reduced accuracy of gantry-to-phantom measurements was attributed to phantom setup errors due to the slightly deformable and flexible phantom extremities. The estimated site-specific safety buffer distance with 0.001% probability of collision for (gantry-to-couch, gantry-to-phantom) was (1.23 cm, 3.35 cm), (1.01 cm, 3.99 cm), and (2.19 cm, 5.73 cm) for treatment to the head, lung, and prostate, respectively. Automated delivery to all three treatment sites was completed in 15 min and collision free using a digital Linac.

CONCLUSIONS

An individualized collision prediction model for the purpose of noncoplanar beam delivery was developed and verified. With the model, the study has demonstrated the feasibility of predicting deliverable beams for an individual patient and then guiding fully automated noncoplanar treatment delivery. This work motivates development of clinical workflows and quality assurance procedures to allow more extensive use and automation of noncoplanar beam geometries.

Initial characterization, dosimetric benchmark and performance validation of Dynamic Wave Arc

Background

Dynamic Wave Arc (DWA) is a clinical approach designed to maximize the versatility of Vero SBRT system by synchronizing the gantry-ring noncoplanar movement with D-MLC optimization. The purpose of this study was to verify the delivery accuracy of DWA approach and to evaluate the potential dosimetric benefits.

Methods

DWA is an extended form of VMAT with a continuous varying ring position. The main difference in the optimization modules of VMAT and DWA is during the angular spacing, where the DWA algorithm does not consider the gantry spacing, but only the Euclidian norm of the ring and gantry angle. A preclinical version of RayStation v4.6 (RaySearch Laboratories, Sweden) was used to create patient specific wave arc trajectories for 31 patients with various anatomical tumor regions (prostate, oligometatstatic cases, centrally-located non-small cell lung cancer (NSCLC) and locally advanced pancreatic cancer-LAPC). DWA was benchmarked against the current clinical approaches and coplanar VMAT. Each plan was evaluated with regards to dose distribution, modulation complexity (MCS), monitor units and treatment time efficiency. The delivery accuracy was evaluated using a 2D diode array that takes in consideration the multi-dimensionality of DWA during dose reconstruction.

Results

In centrally-located NSCLC cases, DWA improved the low dose spillage with 20 %, while the target coverage was increased with 17 % compared to 3D CRT. The structures that significantly benefited from using DWA were proximal bronchus and esophagus, with the maximal dose being reduced by 17 % and 24 %, respectively. For prostate and LAPC, neither technique seemed clearly superior to the other; however, DWA reduced with more than 65 % of the delivery time over IMRT. A steeper dose gradient outside the target was observed for all treatment sites (p < 0.01) with DWA. Except the oligometastatic cases, where the DWA-MCSs indicate a higher modulation, both DWA and VMAT modalities provide plans of similar complexity. The average ɣ (3 % /3 mm) passing rate for DWA plans was 99.2 ± 1 % (range from 96.8 to 100 %).

Conclusions

DWA proven to be a fully functional treatment technique, allowing additional flexibility in dose shaping, while preserving dosimetrically robust delivery and treatment times comparable with coplanar VMAT.

Electronic supplementary material

The online version of this article (doi:10.1186/s13014-016-0633-7) contains supplementary material, which is available to authorized users.

A modular approach to intensity-modulated arc therapy optimization with noncoplanar trajectories

Utilizing noncoplanar beam angles in volumetric modulated arc therapy (VMAT) has the potential to combine the benefits of arc therapy, such as short treatment times, with the benefits of noncoplanar intensity modulated radiotherapy (IMRT) plans, such as improved organ sparing. Recently, vendors introduced treatment machines that allow for simultaneous couch and gantry motion during beam delivery to make noncoplanar VMAT treatments possible. Our aim is to provide a reliable optimization method for noncoplanar isocentric arc therapy plan optimization. The proposed solution is modular in the sense that it can incorporate different existing beam angle selection and coplanar arc therapy optimization methods. Treatment planning is performed in three steps. First, a number of promising noncoplanar beam directions are selected using an iterative beam selection heuristic; these beams serve as anchor points of the arc therapy trajectory. In the second step, continuous gantry/couch angle trajectories are optimized using a simple combinatorial optimization model to define a beam trajectory that efficiently visits each of the anchor points. Treatment time is controlled by limiting the time the beam needs to trace the prescribed trajectory. In the third and final step, an optimal arc therapy plan is found along the prescribed beam trajectory. In principle any existing arc therapy optimization method could be incorporated into this step; for this work we use a sliding window VMAT algorithm. The approach is demonstrated using two particularly challenging cases. The first one is a lung SBRT patient whose planning goals could not be satisfied with fewer than nine noncoplanar IMRT fields when the patient was treated in the clinic. The second one is a brain tumor patient, where the target volume overlaps with the optic nerves and the chiasm and it is directly adjacent to the brainstem. Both cases illustrate that the large number of angles utilized by isocentric noncoplanar VMAT plans can help improve dose conformity, homogeneity, and organ sparing simultaneously using the same beam trajectory length and delivery time as a coplanar VMAT plan.

Noncoplanar VMAT for nasopharyngeal tumors: Plan quality versus treatment time

PURPOSE

The authors investigated the potential of optimized noncoplanar irradiation trajectories for volumetric modulated arc therapy (VMAT) treatments of nasopharyngeal patients and studied the trade-off between treatment plan quality and delivery time in radiation therapy.

METHODS

For three nasopharyngeal patients, the authors generated treatment plans for nine different delivery scenarios using dedicated optimization methods. They compared these scenarios according to dose characteristics, number of beam directions, and estimated delivery times. In particular, the authors generated the following treatment plans: (1) a 4π plan, which is a not sequenced, fluence optimized plan that uses beam directions from approximately 1400 noncoplanar directions and marks a theoretical upper limit of the treatment plan quality, (2) a coplanar 2π plan with 72 coplanar beam directions as pendant to the noncoplanar 4π plan, (3) a coplanar VMAT plan, (4) a coplanar step and shoot (SnS) plan, (5) a beam angle optimized (BAO) coplanar SnS IMRT plan, (6) a noncoplanar BAO SnS plan, (7) a VMAT plan with rotated treatment couch, (8) a noncoplanar VMAT plan with an optimized great circle around the patient, and (9) a noncoplanar BAO VMAT plan with an arbitrary trajectory around the patient.

RESULTS

VMAT using optimized noncoplanar irradiation trajectories reduced the mean and maximum doses in organs at risk compared to coplanar VMAT plans by 19% on average while the target coverage remains constant. A coplanar BAO SnS plan was superior to coplanar SnS or VMAT; however, noncoplanar plans like a noncoplanar BAO SnS plan or noncoplanar VMAT yielded a better plan quality than the best coplanar 2π plan. The treatment plan quality of VMAT plans depended on the length of the trajectory. The delivery times of noncoplanar VMAT plans were estimated to be 6.5 min in average; 1.6 min longer than a coplanar plan but on average 2.8 min faster than a noncoplanar SnS plan with comparable treatment plan quality.

CONCLUSIONS

The authors' study reconfirms the dosimetric benefits of noncoplanar irradiation of nasopharyngeal tumors. Both SnS using optimized noncoplanar beam ensembles and VMAT using an optimized, arbitrary, noncoplanar trajectory enabled dose reductions in organs at risk compared to coplanar SnS and VMAT. Using great circles or simple couch rotations to implement noncoplanar VMAT, however, was not sufficient to yield meaningful improvements in treatment plan quality. The authors estimate that noncoplanar VMAT using arbitrary optimized irradiation trajectories comes at an increased delivery time compared to coplanar VMAT yet at a decreased delivery time compared to noncoplanar SnS IMRT.

Optimization approaches to volumetric modulated arc therapy planning

Volumetric modulated arc therapy (VMAT) has found widespread clinical application in recent years. A large number of treatment planning studies have evaluated the potential for VMAT for different disease sites based on the currently available commercial implementations of VMAT planning. In contrast, literature on the underlying mathematical optimization methods used in treatment planning is scarce. VMAT planning represents a challenging large scale optimization problem. In contrast to fluence map optimization in intensity-modulated radiotherapy planning for static beams, VMAT planning represents a nonconvex optimization problem. In this paper, the authors review the state-of-the-art in VMAT planning from an algorithmic perspective. Different approaches to VMAT optimization, including arc sequencing methods, extensions of direct aperture optimization, and direct optimization of leaf trajectories are reviewed. Their advantages and limitations are outlined and recommendations for improvements are discussed.

Treatment techniques to reduce cardiac irradiation for breast cancer patients treated with breast-conserving surgery and radiation therapy: a review

Thousands of women diagnosed with breast cancer each year receive breast-conserving surgery followed by adjuvant radiation therapy. For women with left-sided breast cancer, there is risk of potential cardiotoxicity from the radiation therapy. As data have become available to quantify the risk of cardiotoxicity from radiation, strategies have also developed to reduce the dose of radiation to the heart without compromising radiation dose to the breast. Several broad categories of techniques to reduce cardiac radiation doses include breath hold techniques, prone positioning, intensity-modulated radiation therapy, and accelerated partial breast irradiation, as well as many small techniques to improve traditional three-dimensional conformal radiation therapy. This review summarizes the published scientific literature on the various techniques to decrease cardiac irradiation in women treated to the left breast for breast cancer after breast-conserving surgery.