Evaluation of clinically applied treatment beams with respect to bunker shielding parameters for a Cyberknife M6.

Compared to a conventional linear accelerator, the Cyberknife (CK) is a unique system with respect to radiation protection shielding and the variety and number of non-coplanar beams are two key components regarding this aspect. In this work, a framework to assess the direction distribution and modulation factor (MF) of clinically applied treatment beams of a CyberKnife M6 is developed. Database filtering options allow studying the influence of different parameters such as collimator types, treatment sites or different bunker sizes. A distribution of monitor units (MU) is generated by projecting treatment beams onto the walls, floor and ceiling of the CyberKnife bunker. This distribution is found to be highly heterogeneous and depending, among other parameters, on the bunker size. For our bunker design, 10%-13% of the MUs are delivered to the right and left wall, each. The floor receives more than 64% of the applied MUs, while the wall behind the patient's head is not hit by primary treatment beams. Between 0% and 5% of the total MUs are delivered to the wall at the patient's feet. This number highly depends on the treatment site, e.g., for extracranial patients no beams hit that wall. Collimator choice was found to have minor influence on the distribution of MUs. On the other hand, the MF depends on the collimator type as well as on the treatment site. The MFs (delivered MU/prescribed dose) for all treatments, all MLC treatments, cranial and extracranial treatments are 8.3, 6.4, 7.7, and 9.9 MU/cGy, respectively. The developed framework allows assessing and monitoring important parameters regarding radiation protection of a CK-M6 using the actually applied treatment beams. Furthermore, it enables evaluating different clinical and constructional situations using the filtering options.

Combining Monte Carlo methods with coherent wave optics for the simulation of phase-sensitive X-ray imaging.

Phase-sensitive X-ray imaging shows a high sensitivity towards electron density variations, making it well suited for imaging of soft tissue matter. However, there are still open questions about the details of the image formation process. Here, a framework for numerical simulations of phase-sensitive X-ray imaging is presented, which takes both particle- and wave-like properties of X-rays into consideration. A split approach is presented where we combine a Monte Carlo method (MC) based sample part with a wave optics simulation based propagation part, leading to a framework that takes both particle- and wave-like properties into account. The framework can be adapted to different phase-sensitive imaging methods and has been validated through comparisons with experiments for grating interferometry and propagation-based imaging. The validation of the framework shows that the combination of wave optics and MC has been successfully implemented and yields good agreement between measurements and simulations. This demonstrates that the physical processes relevant for developing a deeper understanding of scattering in the context of phase-sensitive imaging are modelled in a sufficiently accurate manner. The framework can be used for the simulation of phase-sensitive X-ray imaging, for instance for the simulation of grating interferometry or propagation-based imaging.

Monte Carlo dose calculation for photon and electron radiotherapy on dynamically deforming anatomy.

BACKGROUND: Dose calculation in radiotherapy aims to accurately estimate and assess the dose distribution of a treatment plan. Monte Carlo (MC) dose calculation is considered the gold standard owing to its ability to accurately simulate particle transport in inhomogeneous media. However, uncertainties such as the patient's dynamically deforming anatomy can still lead to differences between the delivered and planned dose distribution. PURPOSE: Development and validation of a deformable voxel geometry for MC dose calculations (DefVoxMC) to account for dynamic deformation in the dose calculation process of photon- and electron-based radiotherapy treatment plans for clinically motivated cases. METHODS: DefVoxMC relies on the subdivision of a regular voxel geometry into dodecahedrons. It allows shifting the dodecahedrons' corner points according to the deformation in the patient's anatomy using deformation vector fields (DVF). DefVoxMC is integrated into the Swiss Monte Carlo Plan (SMCP) to allow the MC dose calculation of photon- and electron-based treatment plans on the deformable voxel geometry. DefVoxMC is validated in two steps. A compression test and a Fano test are performed in silico. Delta4 (for photon beams) and EBT4 film measurements in a cubic PMMA phantom (for electron beams) are performed on a TrueBeam in Developer Mode for clinically motivated treatment plans. During these measurements, table motion is used to mimic rigid dynamic patient motion. The measured and calculated dose distributions are compared using gamma passing rate (GPR) (3% / 2 mm (global), 10% threshold). DefVoxMC is used to study the impact of patient-recorded breathing motion on the dose distribution for clinically motivated lung and breast cases, each prescribed 50 Gy to 50% of the target volume. A volumetric modulated arc therapy (VMAT) and an arc mixed-beam radiotherapy (Arc-MBRT) plan are created for the lung and breast case, respectively. For the dose calculation, the dynamic deformation of the patient's anatomy is described by DVFs obtained from deformable image registration of the different phases of 4DCTs. The resulting dose distributions are compared to the ones of the static situation using dose-volume histograms and dose differences. RESULTS: DefVoxMC is successfully integrated into the SMCP to enable the MC dose calculation of photon- and electron-based treatments on a dynamically deforming patient anatomy. The compression and the Fano test agree within 1.0% and 0.1% with the expected result, respectively. Delta4 and EBT4 film measurements agree with the calculated dose by a GPR >95%. For the clinically motivated cases, breathing motion resulted in areas with a dose increase of up to 26.9 Gy (lung) and up to 7.6 Gy (breast) compared to the static situation. The largest dose differences are observed in high-dose-gradient regions perpendicular to the beam plane, consequently decreasing the planning target volume coverage (V95%) by 4.2% for the lung case and 2.0% for the breast case. CONCLUSIONS: A novel method for MC dose calculation for photon- and electron-based treatments on dynamically deforming anatomy is successfully developed and validated. Applying DefVoxMC to clinically motivated cases, we found that breathing motion has non-negligible impact on the dosimetric plan quality.

A dosimetrically motivated pathfinding approach for non-isocentric dynamic trajectory radiotherapy.

Objective.Non-isocentric dynamic trajectory radiotherapy (DTRT) involves dynamic table translations in synchrony with intensity modulation and dynamic gantry, table, and/or collimator rotation. This work aims to develop and evaluate a novel dosimetrically motivated path determination technique for non-isocentric DTRT.Approach.The path determination considers all available beam directions, given on a user-specified grid of gantry angle, table angle, and longitudinal, vertical, and lateral table position. Additionally, the source-to-target distance of all beam directions can be extended by moving the table away from the gantry along the central beam axis to increase the collision-free space. The path determination uses a column generation algorithm to iteratively add beam directions to paths until a user-defined total path length is reached. A subsequent direct aperture optimization of the intensity modulation along the paths creates deliverable plans. Non-isocentric DTRT plans using the path determination and using a manual path setup were created for a craniospinal and a spinal irradiation case. Furthermore, VMAT, isocentric DTRT, and non-isocentric DTRT plans are created for a breast, head and neck (H&N), and esophagus case. Additionally, a HyperArc plan is created for the H&N case. The plans are compared in terms of the dosimetric treatment plan quality and estimated delivery time.Main results.For the craniospinal and spinal irradiation case, using path determination results in dose distributions with improved conformity but a slightly worse target homogeneity compared to manual path setup. The non-isocentric DTRT plans maintained target coverage while reducing the mean dose to organs-at-risk on average by 1.7 Gy (breast), 1.0 Gy (H&N), and 1.6 Gy (esophagus) compared to the VMAT plans and by 0.8 Gy (breast), 0.6 Gy (H&N), and 0.8 Gy (esophagus) compared to the isocentric DTRT plans.Significance.A general dosimetrically motivated path determination applicable to non-isocentric DTRT plans is successfully developed, further advancing the treatment planning for non-isocentric DTRT.

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.

Dosimetric optimization for dynamic mixed beam arc therapy (DYMBARC).

BACKGROUND: Non-coplanarity and mixed beam modality could be combined to further enhance dosimetric treatment plan quality. We introduce dynamic mixed beam arc therapy (DYMBARC) as an innovative technique that combines non-coplanar photon and electron arcs, dynamic gantry and collimator rotations, and intensity modulation with photon multileaf collimator (MLC). However, finding favorable beam directions for DYMBARC is challenging due to the large solution space, machine component constraints, and optimization parameters, posing a highly non-convex optimization problem. PURPOSE: To establish DYMBARC and solve the pathfinding challenge by employing direct aperture optimization (DAO) to determine the table angles and gantry angle ranges of photon and electron arcs for different clinically motivated cases. METHODS: The method starts by generating a grid of beam directions based on user-defined resolutions along the gantry and table angle axes for each beam quality considered. Beam directions causing collisions or entering through the end of CT are excluded. For electrons, a fixed source-to-surface distance of 80 cm is used to reduce in-air scatter. Electron beam energies with insufficient range to reach the target or beam directions impinging on the table before reaching the patient are excluded. The remaining beam directions form the pathfinding solution space. Promising photon and electron MLC-defined apertures, with associated monitor unit (MU) weights, are iteratively added using a hybrid-DAO algorithm. This algorithm combines column generation to add apertures and simulated annealing to further refine aperture shapes and weights. Apertures are added until the requested number of paths are formed and the user-defined maximum total gantry angle range is reached. Paths are resampled to a finer gantry angle resolution and subject to DAO for simultaneous optimization of beam intensities along the photon/electron arcs. Subsequent final dose calculation and MU weight reoptimization result in a deliverable DYMBARC plan. DYMBARC plans are created for three clinically motivated cases (brain, breast, and pelvis) and compared to DYMBARC variants: colli-DTRT (dynamic collimator trajectory radiotherapy) using non-coplanar photon arcs; and Arc-MBRT (mixed beam radiotherapy) using photons and electrons but restricted to coplanar setup. Additionally, a manually defined volumetric modulated arc therapy (VMAT) setup serves as a reference clinical technique. Dose distributions, dose-volume histograms, and dosimetric endpoints are evaluated. Dosimetric validation with radiochromic film measurements (gamma evaluation, 3% / 2 mm (global), 10% dose threshold) is performed on a TrueBeam system in developer mode for one case. RESULTS: While maintaining similar target coverage and homogeneity, DYMBARC reduced mean doses to organs-at-risk compared to VMAT by an average of 3.2, 0.5, and 2.9 Gy for the brain, breast, and pelvis cases, respectively. Similar or smaller mean dose reductions were observed for Arc-MBRT or colli-DTRT, compared to VMAT. Electron contributions to the mean planning target volume dose ranged from 2% to 34% for DYMBARC and from 11% to 40% for Arc-MBRT. Measurement validation showed >99.7% gamma passing rate. CONCLUSIONS: DYMBARC was successfully established using a dosimetrically optimized pathfinding approach, combining non-coplanarity with mixed beam modality. DYMBARC facilitated the determination of photon and electron contributions on a case-by-case basis, enhancing more personalized treatment modalities.

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

Background and purpose: 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. Materials and methods: 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 %). Results: 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). Conclusions: 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.

Dosimetrically motivated beam-angle optimization for non-coplanar arc radiotherapy with and without dynamic collimator rotation.

BACKGROUND: Non-coplanar techniques have shown to improve the achievable dose distribution compared to standard coplanar techniques for multiple treatment sites but finding optimal beam directions is challenging. Dynamic collimator trajectory radiotherapy (colli-DTRT) is a new intensity modulated radiotherapy technique that uses non-coplanar partial arcs and dynamic collimator rotation. PURPOSE: To solve the beam angle optimization (BAO) problem for colli-DTRT and non-coplanar VMAT (NC-VMAT) by determining the table-angle and the gantry-angle ranges of the partial arcs through iterative 4π fluence map optimization (FMO) and beam direction elimination. METHODS: BAO considers all available beam directions sampled on a gantry-table map with the collimator angle aligned to the superior-inferior axis (colli-DTRT) or static (NC-VMAT). First, FMO is performed, and beam directions are scored based on their contributions to the objective function. The map is thresholded to remove the least contributing beam directions, and arc candidates are formed by adjacent beam directions with the same table angle. Next, FMO and arc candidate trimming, based on objective function penalty score, is performed iteratively until a desired total gantry angle range is reached. Direct aperture optimization on the final set of colli-DTRT or NC-VMAT arcs generates deliverable plans. colli-DTRT and NC-VMAT plans were created for seven clinically-motivated cases with targets in the head and neck (two cases), brain, esophagus, lung, breast, and prostate. colli-DTRT and NC-VMAT were compared to coplanar VMAT plans as well as to class-solution non-coplanar VMAT plans for the brain and head and neck cases. Dosimetric validation was performed for one colli-DTRT (head and neck) and one NC-VMAT (breast) plan using film measurements. RESULTS: Target coverage and conformity was similar for all techniques. colli-DTRT and NC-VMAT plans had improved dosimetric performance compared to coplanar VMAT for all treatment sites except prostate where all techniques were equivalent. For the head and neck and brain cases, mean dose reduction-in percentage of the prescription dose-to parallel organs was on average 0.7% (colli-DTRT), 0.8% (NC-VMAT) and 0.4% (class-solution) compared to VMAT. The reduction in D2% for the serial organs was on average 1.7% (colli-DTRT), 2.0% (NC-VMAT) and 0.9% (class-solution). For the esophagus, lung, and breast cases, mean dose reduction to parallel organs was on average 0.2% (colli-DTRT) and 0.3% (NC-VMAT) compared to VMAT. The reduction in D2% for the serial organs was on average 1.3% (colli-DTRT) and 0.9% (NC-VMAT). Estimated delivery times for colli-DTRT and NC-VMAT were below 4 min for a full gantry angle range of 720°, including transitions between arcs, except for the brain case where multiple arcs covered the whole table angle range. These times are in the same order as the class-solution for the head and neck and brain cases. Total optimization times were 25%-107% longer for colli-DTRT, including BAO, compared to VMAT. CONCLUSIONS: We successfully developed dosimetrically motivated BAO for colli-DTRT and NC-VMAT treatment planning. colli-DTRT and NC-VMAT are applicable to multiple treatment sites, including body sites, with beneficial or equivalent dosimetric performances compared to coplanar VMAT and reasonable delivery times.

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

BACKGROUND: 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. PURPOSE: 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. METHODS: 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. RESULTS: 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. CONCLUSION: 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.

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

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

Semi-classical Monte Carlo algorithm for the simulation of X-ray grating interferometry.

Traditional simulation techniques such as wave optics methods and Monte Carlo (MC) particle transport cannot model both interference and inelastic scattering phenomena within one framework. Based on the rules of quantum mechanics to calculate probabilities, we propose a new semi-classical MC algorithm for efficient and simultaneous modeling of scattering and interference processes. The similarities to MC particle transport allow the implementation as a flexible c++ object oriented extension of EGSnrc-a well-established MC toolkit. In addition to previously proposed Huygens principle based transport through optics components, new variance reduction techniques for the transport through gratings are presented as transport options to achieve the required improvement in speed and memory costs necessary for an efficient exploration (system design-dose estimations) of the medical implementation of X-ray grating interferometry (GI), an emerging imaging technique currently subject of tremendous efforts towards clinical translation. The feasibility of simulation of interference effects is confirmed in four academic cases and an experimental table-top GI setup. Comparison with conventional MC transport show that deposited energy features of EGSnrc are conserved.

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

BACKGROUND: Evaluating plan robustness is a key step in radiotherapy. PURPOSE: 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. METHODS: 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. RESULTS: 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). CONCLUSIONS: 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.

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, D0.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.

Monte Carlo implementation, validation, and characterization of a 120 leaf MLC.

PURPOSE: Recently, the new high definition multileaf collimator (HD120 MLC) was commercialized by Varian Medical Systems providing high resolution in the center section of the treatment field. The aim of this work is to investigate the characteristics of the HD120 MLC using Monte Carlo (MC) methods. METHODS: Based on the information of the manufacturer, the HD120 MLC was implemented into the already existing Swiss MC Plan (SMCP). The implementation has been configured by adjusting the physical density and the air gap between adjacent leaves in order to match transmission profile measurements for 6 and 15 MV beams of a Novalis TX. These measurements have been performed in water using gafchromic films and an ionization chamber at an SSD of 95 cm and a depth of 5 cm. The implementation was validated by comparing diamond measured and calculated penumbra values (80%-20%) for different field sizes and water depths. Additionally, measured and calculated dose distributions for a head and neck IMRT case using the DELTA(4) phantom have been compared. The validated HD120 MLC implementation has been used for its physical characterization. For this purpose, phase space (PS) files have been generated below the fully closed multileaf collimator (MLC) of a 40 × 22 cm(2) field size for 6 and 15 MV. The PS files have been analyzed in terms of energy spectra, mean energy, fluence, and energy fluence in the direction perpendicular to the MLC leaves and have been compared with the corresponding data using the well established Varian 80 leaf (MLC80) and Millennium M120 (M120 MLC) MLCs. Additionally, the impact of the tongue and groove design of the MLCs on dose has been characterized. RESULTS: Calculated transmission values for the HD120 MLC are 1.25% and 1.34% in the central part of the field for the 6 and 15 MV beam, respectively. The corresponding ionization chamber measurements result in a transmission of 1.20% and 1.35%. Good agreement has been found for the comparison between transmission profiles resulting from MC simulations and film measurements. The simulated and measured values for the penumbra agreed within <0.5 mm for all field sizes, depths, and beam energies, and a good agreement has been found between the measured and the calculated dose distributions for the IMRT case. The total energy spectra are almost identical for the three MLCs. However, the mean energy, fluence and energy fluence are significantly different. Due to the different leaf widths of the MLCs, the shape of these distributions is different, each representing its leave structure. Due to the increase in width from the inner to the outer HD120 MLC leaves, the fluence and energy fluence clearly decrease below the outer leaves. The MLC80 and the M120 MLC resulted in an increase of the fluence and energy fluence compared with those resulted for the HD120 MLC. The dose reduction can exceed 20% compared with the dose of the open field due to the tongue and groove design of the HD120 MLC. CONCLUSIONS: The HD120 MLC has been successfully implemented into the SMCP. Comparisons between MC calculations and measurements show very good agreement. The SMCP is now able to calculate accurate dose distributions for treatment plans using the HD120 MLC.

VMC++ validation for photon beams in the energy range of 20-1000 keV.

PURPOSE: In high energy teletherapy, VMC++ is known to be a very accurate and efficient Monte Carlo (MC) code. In principle, the MC method is also a powerful dose calculation tool in other areas in radiation oncology, e.g., brachytherapy or orthovoltage radiotherapy. However, VMC++ is not validated for the low-energy range of such applications. This work aims in the validation of the VMC++ MC code for photon beams in the energy range between 20 and 1000 keV. METHODS: Dose calculations were performed in different 40 x 40 x 40 cm3 phantoms of different materials. Dose distributions of monoenergetic (ranging from 20 to 1000 keV) 10 x 10 and 2 x 2 cm2 parallel beams were calculated. Voxel sizes of 4 x 4 x 4 and 1 x 1 x 1 mm3 were used for the dose calculations. The resulting dose distributions were compared to those calculated using EGSnrc, which is used as a golden standard in this work. RESULTS: At energies between 100 and 1000 keV, EGSnrc and VMC++ calculated dose distributions agree within the statistical uncertainty of about 1% (1sigma). At energies < or = 50 keV, dose differences of up to 1.6% (in % of D(max)) occur when VMC++ and EGSnrc are compared. Turning off Rayleigh scattering, binding effects for Compton scattering, and the atomic relaxation after photoelectric absorption in EGSnrc (all not implemented in VMC++) leads to an agreement between both MC codes within statistical uncertainty. Further, using the KERMA approximation feature implemented in VMC++ leads to very efficient simulations in the energy range between 20 and 1000 keV. CONCLUSIONS: Further improvements for very low energies in accuracy of VMC++ could be achieved by implementing Rayleigh scattering, binding effects for Compton scattering, and the atomic relaxation after photoelectric absorption. Implementation into VMC++ of KERMA approximation has been validated.

Dose calculation of dynamic trajectory radiotherapy using Monte Carlo.

PURPOSE: Using volumetric modulated arc therapy (VMAT) delivery technique gantry position, multi-leaf collimator (MLC) as well as dose rate change dynamically during the application. However, additional components can be dynamically altered throughout the dose delivery such as the collimator or the couch. Thus, the degrees of freedom increase allowing almost arbitrary dynamic trajectories for the beam. While the dose delivery of such dynamic trajectories for linear accelerators is technically possible, there is currently no dose calculation and validation tool available. Thus, the aim of this work is to develop a dose calculation and verification tool for dynamic trajectories using Monte Carlo (MC) methods. METHODS: The dose calculation for dynamic trajectories is implemented in the previously developed Swiss Monte Carlo Plan (SMCP). SMCP interfaces the treatment planning system Eclipse with a MC dose calculation algorithm and is already able to handle dynamic MLC and gantry rotations. Hence, the additional dynamic components, namely the collimator and the couch, are described similarly to the dynamic MLC by defining data pairs of positions of the dynamic component and the corresponding MU-fractions. For validation purposes, measurements are performed with the Delta4 phantom and film measurements using the developer mode on a TrueBeam linear accelerator. These measured dose distributions are then compared with the corresponding calculations using SMCP. First, simple academic cases applying one-dimensional movements are investigated and second, more complex dynamic trajectories with several simultaneously moving components are compared considering academic cases as well as a clinically motivated prostate case. RESULTS: The dose calculation for dynamic trajectories is successfully implemented into SMCP. The comparisons between the measured and calculated dose distributions for the simple as well as for the more complex situations show an agreement which is generally within 3% of the maximum dose or 3mm. The required computation time for the dose calculation remains the same when the additional dynamic moving components are included. CONCLUSION: The results obtained for the dose comparisons for simple and complex situations suggest that the extended SMCP is an accurate dose calculation and efficient verification tool for dynamic trajectory radiotherapy. This work was supported by Varian Medical Systems.

Adaptive step size algorithm to increase efficiency of proton macro Monte Carlo dose calculation.

PURPOSE: To provide fast and accurate dose calculation in voxelized geometries for proton radiation therapy by implementing an adaptive step size algorithm in the proton macro Monte Carlo (pMMC) method. METHODS: The in-house developed local-to-global MMC method for proton dose calculation is extended with an adaptive step size algorithm for efficient proton transport through a voxelized geometry by sampling transport parameters from a pre-simulated database. Adaptive choice of an adequate slab size in dependence of material interfaces in the proton's longitudinal and lateral vicinity is investigated. The dose calculation algorithm is validated against the non-adaptive pMMC and full MC simulation for pencil and broad beams with various energies impinging on academic phantoms as well as a head and neck patient CT. RESULTS: For material interfaces perpendicular to a proton's direction, choice of nearest neighbor slab thickness shows best trade-off between dosimetric accuracy and calculation efficiency. Adaptive reduction of chosen slab size is shown to be required for material interfaces closer than 0.5 mm in lateral direction. For the academic phantoms, dose differences of within 1% or 1 mm compared to full Geant4 MC simulation are found, while achieving an efficiency gain of up to a factor of 5.6 compared to the non-adaptive algorithm and 284 compared to Geant4. For the head and neck patient CT, dose differences are within 1% or 1 mm with an efficiency gain factor of up to 3.4 compared to the non-adaptive algorithm and 145 compared to Geant4. CONCLUSION: An adaptive step size algorithm for proton macro Monte Carlo was implemented and evaluated. The dose calculation provides the accuracy of full MC simulations, while achieving an efficiency gain factor of three compared to the non-adaptive algorithm and two orders of magnitude compared to full MC for a complex patient CT.

Whole-ventricular irradiation for intracranial germ cell tumors: Dosimetric comparison of pencil beam scanned protons, intensity-modulated radiotherapy and volumetric-modulated arc therapy.

BACKGROUND: Whole-ventricular radiotherapy (WV-RT) followed by a boost to the tumor bed (WV-RT/TB) is recommended for intracranial germ cell tumors (IGCT). As the critical brain areas are mainly in the target volume vicinity, it is unclear if protons indeed substantially spare neurofunctional organs at risk (NOAR). Therefore, a dosimetric comparison study of WV-RT/TB was conducted to assess whether proton or photon radiotherapy achieves better NOAR sparing. METHODS: Eleven children with GCT received 24 Gy(RBE) WV-RT and a boost up to 40 Gy(RBE) in 25 fractions of 1.6 Gy(RBE) with pencil beam scanning proton therapy (PBS-PT). Intensity-modulated radiotherapy (IMRT) and volumetric-modulated arc therapy (VMAT) plans were generated for these patients. NOAR were delineated and treatment plans were compared for target volume coverage (TVC), homogeneity index (HI), inhomogeneity coefficient (IC) and (N)OAR sparing. RESULTS: TVC was comparable for all three modalities. Compared to IMRT and VMAT, PBS-PT showed statistically significant optimized IC, as well as dose reduction, among others, in mean and integral dose to the: normal brain (-35.2%, -32.7%; -35.2%, -33.0%, respectively), cerebellum (-53.7%, -33.1%; -53.6%, -32.7%) and right temporal lobe (-14.5%, -31.9%; -14.7%, -29.9%). The Willis' circle was better protected with PBS-PT than IMRT (-7.1%; -7.8%). The left hippocampus sparing was higher with IMRT. Compared to VMAT, the dose to the hippocampi, amygdalae and temporal lobes was significantly decreased in the IMRT plans. CONCLUSIONS: Dosimetric comparison of WV-RT/TB in IGCT suggests PBS-PT's advantage over photons in conformality and NOAR sparing, whereas IMRT's superiority over VMAT, thus potentially minimizing long-term sequelae.

Evaluation of a new software prototype for frameless radiosurgery of arteriovenous malformations.

BACKGROUND: In order to locate an arteriovenous malformation, typically, a digital subtraction angiography (DSA) is carried out. To use the DSA for target definition an accurate image registration between CT and DSA is required. Carrying out a non-invasive, frameless procedure, registration of the 2D-DSA images with the CT is critical. A new software prototype is enabling this frameless procedure. The aim of this work was to evaluate the prototype in terms of targeting accuracy and reliability based on phantom measurements as well as with the aid of patient data. In addition, the user's ability to recognize registration mismatches and quality was assessed. METHODS: Targeting accuracy was measured with a simple cubic, as well as with an anthropomorphic head phantom. Clearly defined academic targets within the phantoms were contoured on the CT. These reference structures were compared with the structures generated within the prototype. A similar approach was used with patient data, where the clinically contoured target served as the reference structure. An important error source decreasing the target accuracy comes from registration errors between CT and 2D-DSA. For that reason, the tools in BC provided to the user to check these registrations are very important. In order to check if the user is able to recognize registration errors, a set of different registration errors was introduced to the correctly registered CT and 2D-DSA image data sets of three different patients. Each of six different users rated the whole set of registrations within the prototype. RESULTS: The target accuracy of the prototype was found to be below 0.04 cm for the cubic phantom and below 0.05 cm for the anthropomorphic head phantom. The mean target accuracy for the 15 patient cases was found to be below 0.3 cm. In the registration verification part, almost all introduced registration errors above 1° or 0.1 cm were detected by the six users. Nevertheless, in order to quantify and categorize the possibility to detect mismatches in the registration process more data needs to be evaluated. CONCLUSION: Our study shows, that the prototype is a useful tool that has the potential to fill the gap towards a frameless procedure when treating AVMs with the aid of 2D-DSA images in radiosurgery. The target accuracy of the prototype is similar to other systems already established in clinical routine.