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

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

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

Efficiency enhancements of a Monte Carlo beamlet based treatment planning process: implementation and parameter study.

OBJECTIVE The computational effort to perform beamlet calculation, plan optimization and final dose calculation of a treatment planning process (TPP) generating intensity modulated treatment plans is enormous, especially if Monte Carlo (MC) simulations are used for dose calculation. The goal of this work is to improve the computational efficiency of a fully MC based TPP for static and dynamic photon, electron and mixed photon-electron treatment techniques by implementing multiple methods and studying the influence of their parameters. APPROACH A framework is implemented calculating MC beamlets efficiently in parallel on each available CPU core. The user can specify the desired statistical uncertainty of the beamlets, a fractional sparse dose threshold to save beamlets in a sparse format and minimal distances to the PTV surface from which 2x2x2=8 (medium) or even 4x4x4=64 (large) voxels are merged. The compromise between final plan quality and computational efficiency of beamlet calculation and optimization is studied for several parameter values to find a reasonable trade-off. For this purpose, four clinical and one academic case are considered with different treatment techniques. MAIN RESULTS Setting the statistical uncertainty to 5% (photon beamlets) and 15% (electron beamlets), the fractional sparse dose threshold relative to the maximal beamlet dose to 0.1% and minimal distances for medium and large voxels to the PTV to 1 cm and 2 cm, respectively, does not lead to substantial degradation in final plan quality. Only OAR sparing is slightly degraded. Furthermore, computation times are reduced by about 58% (photon beamlets), 88% (electron beamlets) and 96% (optimization) compared to using 2.5% (photon beamlets) and 5% (electron beamlets) statistical uncertainty and no sparse format nor voxel merging. SIGNIFICANCE Several methods are implemented improving computational efficiency of beamlet calculation and plan optimization of a fully MC based TPP without substantial degradation in final plan quality

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

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

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 ( ) between two adjacent control points ( ) and maximum absolute change of G ( ). Four DTRT plans are created with different maximum G&∆G: &  = 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 in the lateral, longitudinal, and vertical directions, and machine uncertainties between 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 to investigate the impact of the GT rotation gradient on plan quality is successfully demonstrated. With increasing , slight (case C, : up to‍-1‍Gy) and substantial (case A, : up to -9.3 Gy, case‍B, : 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 96% (2%‍global/2 mm Gamma passing rate) with the dose calculation. With increasing , 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

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

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

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

Assessing dose rate distributions in VMAT plans

Dose rate is an essential factor in radiobiology. As modern radiotherapy delivery techniques such as volumetric modulated arc therapy (VMAT) introduce dynamic modulation of the dose rate, it is important to assess the changes in dose rate. Both the rate of monitor units per minute (MU rate) and collimation are varied over the course of a fraction, leading to different dose rates in every voxel of the calculation volume at any point in time during dose delivery. Given the radiotherapy plan and machine specific limitations, a VMAT treatment plan can be split into arc sectors between Digital Imaging and Communications in Medicine control points (CPs) of constant and known MU rate. By calculating dose distributions in each of these arc sectors independently and multiplying them with the MU rate, the dose rate in every single voxel at every time point during the fraction can be calculated. Independently calculated and then summed dose distributions per arc sector were compared to the whole arc dose calculation for validation. Dose measurements and video analysis were performed to validate the calculated datasets. A clinical head and neck, cranial and liver case were analyzed using the tool developed. Measurement validation of synthetic test cases showed linac agreement to precalculated arc sector times within ±0.4 s and doses ±0.1 MU (one standard deviation). Two methods for the visualization of dose rate datasets were developed: the first method plots a two-dimensional (2D) histogram of the number of voxels receiving a given dose rate over the course of the arc treatment delivery. In similarity to treatment planning system display of dose, the second method displays the dose rate as color wash on top of the corresponding computed tomography image, allowing the user to scroll through the variation over time. Examining clinical cases showed dose rates spread over a continuous spectrum, with mean dose rates hardly exceeding 100 cGy min(-1) for conventional fractionation. A tool to analyze dose rate distributions in VMAT plans with sub-second accuracy was successfully developed and validated. Dose rates encountered in clinical VMAT test cases show a continuous spectrum with a mean less than or near 100 cGy min(-1) for conventional fractionation