Development of a Monte Carlo based robustness calculation and evaluation tool.

Evaluating plan robustness is a key step in radiotherapy. To develop a flexible Monte Carlo (MC)-based robustness calculation and evaluation tool to assess and quantify dosimetric robustness of intensity-modulated radiotherapy (IMRT) treatment plans by exploring the impact of systematic and random uncertainties resulting from patient setup, patient anatomy changes, and mechanical limitations of machine components. The robustness tool consists of two parts: the first part includes automated MC dose calculation of multiple user-defined uncertainty scenarios to populate a robustness space. An uncertainty scenario is defined by a certain combination of uncertainties in patient setup, rigid intrafraction motion and in mechanical steering of the following machine components: angles of gantry, collimator, table-yaw, table-pitch, table-roll, translational positions of jaws, multileaf-collimator (MLC) banks, and single MLC leaves. The Swiss Monte Carlo Plan (SMCP) is integrated in this tool to serve as the backbone for the MC dose calculations incorporating the uncertainties. The calculated dose distributions serve as input for the second part of the tool, handling the quantitative evaluation of the dosimetric impact of the uncertainties. A graphical user interface (GUI) is developed to simultaneously evaluate the uncertainty scenarios according to user-specified conditions based on dose-volume histogram (DVH) parameters, fast and exact gamma analysis, and dose differences. Additionally, a robustness index (RI) is introduced with the aim to simultaneously evaluate and condense dosimetric robustness against multiple uncertainties into one number. The RI is defined as the ratio of scenarios passing the conditions on the dose distributions. Weighting of the scenarios in the robustness space is possible to consider their likelihood of occurrence. The robustness tool is applied on IMRT, a volumetric modulated arc therapy (VMAT), a dynamic trajectory radiotherapy (DTRT), and a dynamic mixed beam radiotherapy (DYMBER) plan for a brain case to evaluate the robustness to uncertainties of gantry-, table-, collimator angle, MLC, and intrafraction motion. Additionally, the robustness of the IMRT, VMAT, and DTRT plan against patient setup uncertainties are compared. The robustness tool is validated by Delta4 measurements for scenarios including all uncertainty types available. The robustness tool performs simultaneous calculation of uncertainty scenarios, and the GUI enables their fast evaluation. For all evaluated plans and uncertainties, the planning target volume (PTV) margin prevented major clinical target volume (CTV) coverage deterioration (maximum observed standard deviation of D 98 % CTV $D98{\% _{{\rm{CTV}}}}$ was 1.3 Gy). OARs close to the PTV experienced larger dosimetric deviations (maximum observed standard deviation of D 2 % chiasma $D2{\% _{{\rm{chiasma}}}}$ was 14.5 Gy). Robustness comparison by RI evaluation against patient setup uncertainties revealed better dosimetric robustness of the VMAT and DTRT plans as compared to the IMRT plan. Delta4 validation measurements agreed with calculations by >96% gamma-passing rate (3% global/2 mm). The robustness tool was successfully implemented. Calculation and evaluation of uncertainty scenarios with the robustness tool were demonstrated on a brain case. Effects of patient and machine-specific uncertainties and the combination thereof on the dose distribution are evaluated in a user-friendly GUI to quantitatively assess and compare treatment plans and their robustness.

Loebner HA ，Volken W ，Mueller S ，Bertholet J ，Mackeprang PH ，Guyer G ，Aebersold DM ，Stampanoni MFM ，Manser P ，Fix MK ... -  《-》

Impact of the gradient in gantry-table rotation on dynamic trajectory radiotherapy plan quality.

To improve organ at risk (OAR) sparing, dynamic trajectory radiotherapy (DTRT) extends VMAT by dynamic table and collimator rotation during beam-on. However, comprehensive investigations regarding the impact of the gantry-table (GT) rotation gradient on the DTRT plan quality have not been conducted. To investigate the impact of a user-defined GT rotation gradient on plan quality of DTRT plans in terms of dosimetric plan quality, dosimetric robustness, deliverability, and delivery time. The dynamic trajectories of DTRT are described by GT and gantry-collimator paths. The GT path is determined by minimizing the overlap of OARs with planning target volume (PTV). This approach is extended to consider a GT rotation gradient by means of a maximum gradient of the path ( G m a x ${G}_{max}$ ) between two adjacent control points ( G = | Δ table angle / Δ gantry angle | $G = | \Delta {{\mathrm{table\ angle}}/\Delta {\mathrm{gantry\ angle}}} |$ ) and maximum absolute change of G ( Δ G m a x ${{\Delta}}{G}_{max}$ ). Four DTRT plans are created with different maximum G&∆G: G m a x ${G}_{max}$ & Δ G m a x ${{\Delta}}{G}_{max}$  = 0.5&0.125 (DTRT-1), 1&0.125 (DTRT-2), 3&0.125 (DTRT-3) and 3&1‍(DTRT-4), including 3-4 dynamic trajectories, for three clinically motivated cases in the head and neck and brain region (A, B, and C). A reference VMAT plan for each case is created. For all plans, plan quality is assessed and compared. Dosimetric plan quality is evaluated by target coverage, conformity, and OAR sparing. Dosimetric robustness is evaluated against systematic and random patient-setup uncertainties between ± 3 mm $ \pm 3\ {\mathrm{mm}}$ in the lateral, longitudinal, and vertical directions, and machine uncertainties between ± 4 ∘ $ \pm 4^\circ \ $ in the dynamically rotating machine components (gantry, table, collimator rotation). Delivery time is recorded. Deliverability and delivery accuracy on a TrueBeam are assessed by logfile analysis for all plans and additionally verified by film measurements for one case. All dose calculations are Monte Carlo based. The extension of the DTRT planning process with user-defined G m a x & Δ G m a x ${G}_{max}\& {{\Delta}}{G}_{max}$ to investigate the impact of the GT rotation gradient on plan quality is successfully demonstrated. With increasing G m a x & Δ G m a x ${G}_{max}\& {{\Delta}}{G}_{max}$ , slight (case C, D m e a n , p a r o t i d l . ${D}_{mean,\ parotid\ l.}$ : up to‍-1‍Gy) and substantial (case A, D 0.03 c m 3 , o p t i c n e r v e r . ${D}_{0.03c{m}^3,\ optic\ nerve\ r.}$ : up to -9.3 Gy, case‍B, D m e a n , b r a i n $\ {D}_{mean,\ brain}$ : up to -4.7‍Gy) improvements in OAR sparing are observed compared to VMAT, while maintaining similar target coverage. All plans are delivered on the TrueBeam. Expected and actual machine position values recorded in the logfiles deviated by <0.2° for gantry, table and collimator rotation. The film measurements agreed by >96% (2%‍global/2 mm Gamma passing rate) with the dose calculation. With increasing G m a x & Δ G m a x ${G}_{max}\& {{\Delta}}{G}_{max}$ , delivery time is prolonged by <2 min/trajectory (DTRT-4) compared to VMAT and DTRT-1. The DTRT plans for case A and B and the VMAT plan for case C plan reveal the best dosimetric robustness for the considered uncertainties. The impact of the GT rotation gradient on DTRT plan quality is comprehensively investigated for three cases in the head and neck and brain region. Increasing freedom in this gradient improves dosimetric plan quality at the cost of increased delivery time for the investigated cases. No clear dependency of GT rotation gradient on dosimetric robustness is observed.

Loebner HA ，Mueller S ，Volken W ，Wallimann P ，Aebersold DM ，Stampanoni MFM ，Fix MK ，Manser P ... -  《-》

Part 1: Optimization and evaluation of dynamic trajectory radiotherapy.

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. 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. 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%. 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.

Fix MK ，Frei D ，Volken W ，Terribilini D ，Mueller S ，Elicin O ，Hemmatazad H ，Aebersold DM ，Manser P ... -  《-》

Robustness analysis of dynamic trajectory radiotherapy and volumetric modulated arc therapy plans for head and neck cancer.

Dynamic trajectory radiotherapy (DTRT) has been shown to improve healthy tissue sparing compared to volumetric arc therapy (VMAT). This study aimed to assess and compare the robustness of DTRT and VMAT treatment-plans for head and neck (H&N) cancer to patient-setup (PS) and machine-positioning uncertainties. The robustness of DTRT and VMAT plans previously created for 46 H&N cases, prescribed 50-70 Gy to 95 % of the planning-target-volume, was assessed. For this purpose, dose distributions were recalculated using Monte Carlo, including uncertainties in PS (translation and rotation) and machine-positioning (gantry-, table-, collimator-rotation and multi-leaf collimator (MLC)). Plan robustness was evaluated by the uncertainties' impact on normal tissue complication probabilities (NTCP) for xerostomia and dysphagia and on dose-volume endpoints. Differences between DTRT and VMAT plan robustness were compared using Wilcoxon matched-pair signed-rank test (α = 5 %). Average NTCP for moderate-to-severe xerostomia and grade ≥ II dysphagia was lower for DTRT than VMAT in the nominal scenario (0.5 %, p = 0.01; 2.1 %, p < 0.01) and for all investigated uncertainties, except MLC positioning, where the difference was not significant. Average differences compared to the nominal scenario were ≤ 3.5 Gy for rotational PS (≤ 3°) and machine-positioning (≤ 2°) uncertainties, <7 Gy for translational PS uncertainties (≤ 5 mm) and < 20 Gy for MLC-positioning uncertainties (≤ 5 mm). DTRT and VMAT plan robustness to the investigated uncertainties depended on uncertainty direction and location of the structure-of-interest to the target. NTCP remained on average lower for DTRT than VMAT even when considering uncertainties.

Loebner HA ，Bertholet J ，Mackeprang PH ，Volken W ，Elicin O ，Mueller S ，Guyer G ，Aebersold DM ，Stampanoni MFM ，Fix MK ，Manser P ... -  《-》

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.

Bertholet J ，Guyer G ，Mueller S ，Loebner HA ，Volken W ，Aebersold DM ，Manser P ，Fix MK ... -  《-》