Search

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 (Formula presented.) was 1.3 Gy). OARs close to the PTV experienced larger dosimetric deviations (maximum observed standard deviation of (Formula presented.) 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., Medical Physics, 49 (7), ISSN:0094-2405, ISSN:2473-4209, ISSN:1522-8541

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.

Purpose The objectives of the work presented in this paper were to (1) implement a robust-optimization method for deliverable mixed-beam radiotherapy (MBRT) plans within a previously developed MBRT planning framework; (2) perform an experimental validation of the delivery of robust-optimized MBRT plans; and (3) compare PTV-based and robust-optimized MBRT plans in terms of target dose robustness and organs at risk (OAR) sparing for clinical head and neck and brain patient cases. Methods A robust-optimization method, which accounts for translational setup errors, was implemented within a previously developed treatment planning framework for MBRT. The framework uses a hybrid direct aperture optimization method combining column generation and simulated annealing. A robust plan was developed and then delivered to an anthropomorphic head phantom using the Developer Mode of a TrueBeam linac. Planar dose distributions were measured and compared to the planned dose. Robust-optimized and PTV-based plans were developed for three clinical patient cases consisting of two head and neck cases and one brain case. The plans were compared in terms of the robustness to 5 mm shifts of the target volume dose as well as in terms of OAR sparing. Results Using a gamma criterion of 3%/2 mm and a dose threshold of 10%, the agreement between film measurements and dose calculations was better than 97.7% for the total plan and better than 95.5% for the electron component of the plan. For the two head and neck patient cases, the average clinical target volume (CTV) dose homogeneity index (V95%–V107%) over all the considered setup error scenarios was on average 19% lower for the PTV-based plans and it had a larger standard deviation. The robust-optimized plans achieved, on average, a 20% reduction in the OAR doses compared to the PTV-based plans. For the brain patient case, the CTV dose homogeneity index was similar for the two plans, while the OAR doses were 22% lower, on average, for the robust-optimized plan. No clear trend in terms of electron contributions was found across the three patient cases, although robust-optimized plans tended toward higher electron beam energies. Conclusions A framework for robust optimization of deliverable MBRT plans has been developed and validated. PTV-based MBRT were found to not be robust to setup errors, while the dose delivered by the robust-optimized plans were clinically acceptable for all considered error scenarios and had better OAR sparing. This study shows that the robust optimization is a promising alternative to conventional PTV margins for MBRT., Medical Physics, 48 (11), ISSN:0094-2405, ISSN:2473-4209, ISSN:1522-8541

Background Stereotactic radiosurgery (SRS) has been recognized as a first-line treatment option for small to moderate sized vestibular schwannoma (VS). Our aim is to evaluate the impact of SRS doses and other patient and disease characteristics on vestibular function in patients with VS. Methods Data on VS patients treated with single-fraction SRS to 12 Gy were retrospectively reviewed. No dose constraints were given to the vestibule during optimization in treatment planning. Patient and tumor characteristics, pre- and post-SRS vestibular examination results and patient-reported dizziness were assessed from patient records. Results Fifty-three patients were analyzed. Median follow-up was 32 months (range, 6–79). The median minimum, mean and maximum vestibular doses were 2.6 ± 1.6 Gy, 6.7 ± 2.8 Gy, and 11 ± 3.6 Gy, respectively. On univariate analysis, Koos grade (p = 0.04; OR: 3.45; 95% CI 1.01–11.81), tumor volume (median 6.1 cm3; range, 0.8–38; p = 0.01; OR: 4.85; 95% CI 1.43–16.49), presence of pre-SRS dizziness (p = 0.02; OR: 3.98; 95% CI 1.19–13.24) and minimum vestibular dose (p = 0.033; OR: 1.55; 95% CI 1.03–2.32) showed a significant association with patient-reported dizziness. On multivariate analysis, minimum vestibular dose remained significant (p = 0.02; OR: 1.75; 95% CI 1.05–2.89). Patients with improved caloric function had received significantly lower mean (1.5 ± 0.7 Gy, p = 0.01) and maximum doses (4 ± 1.5 Gy, p = 0.01) to the vestibule. Conclusions Our results reveal that 5 Gy and above minimum vestibular doses significantly worsened dizziness. Additionally, mean and maximum doses received by the vestibule were significantly lower in patients who had improved caloric function. Further investigations are needed to determine dose-volume parameters and their effects on vestibular toxicity.

Die DVT unterscheidet sich von allen anderen radiologischen Verfahren: Auf der einen Seite von den „einfacheren“ Methoden wie intraoralen Rontgenaufnahmen und Panoramaschichtaufnahmen. Beides zweidimensionale Techniken, welche stets mit Uberlagerungseffekten und geometrischen Verzerrungen belastet sind. Auf der anderen Seite stehen „aufwendigere“ Methoden wie Computertomographie, Kernspintomographie, Szintigraphie und auch Positronenemissionstomographie. Sie alle sind in Teilbereichen resp. fur entsprechende Spezialfragestellungen der DVT uberlegen. Bezuglich Verfugbarkeit, Kosten und teilweise auch Strahlenbelastung des Patienten stehen sie aber im Verhaltnis zur DVT ungunstig da. An eine Chair-side-Anwendung, wie sie in zahnarztlichen Praxen regelmasig gefragt ist, ist bei diesen Techniken nicht zu denken. Beim Thema Bildqualitat ist zu beachten, dass ein subjektiv „schoneres“ Bild dem weniger schonen nicht unbedingt objektiv diagnostisch uberlegen ist.

MOTIVATION Progress in the field of MR-guided radiotherapy has triggered the need for fast and accurate dose calculation in presence of magnetic fields. The aim of this work is to satisfy this need by extending the Macro Monte Carlo (MMC) method to enable dose calculation for photon, electron and proton beams in a magnetic field. METHODS The MMC method is based on the transport of particles in macroscopic steps through an absorber by sampling the relevant physical quantities from a pre-calculated database containing probability distribution functions. To enable MMC particle transport in a magnetic field, a transformation accounting for the Lorentz force is applied for each macro step by rotating the sampled position and direction around the magnetic field vector. The transformed position and direction distributions on local geometries are validated against full MC for electron and proton pencil beams. To enable photon dose calculation, an in-house MC algorithm is used for photon transport and interaction. Emerging secondary charged particles are passed to MMC for transport and energy deposition. The extended MMC dose calculation accuracy and efficiency is assessed by comparison with EGSnrc (photon and electron beams) and Geant4 (proton beam) calculated dose distributions of different energies and homogeneous magnetic fields for broad beams impinging on water phantoms with bone and lung inhomogeneities. RESULTS The geometric transformation on the local geometries is able to reproduce the results of full MC for all investigated settings (difference in mean value and standard deviation < 1%). MMC calculated dose distributions in a homogeneous magnetic field are in agreement with EGSnrc and Geant4, respectively, with gamma passing rates > 99.6% (global 2%, 2 mm and 10% threshold criteria) for all situations. MMC achieves a substantial efficiency gain of up to a factor of 21 (photon beam), 66 (electron beam) and 356 (proton beam) compared to EGSnrc or Geant4. CONCLUSION Efficient and accurate dose calculation in magnetic fields was successfully enabled by utilizing the developed extended MMC transport method for photon, electron and proton beams.

Purpose. To develop a novel treatment planning process (TPP) with simultaneous optimization of modulated photon, electron and proton beams for improved treatment plan quality in radiotherapy. Methods. A framework for fluence map optimization of Monte Carlo (MC) calculated beamlet dose distributions is developed to generate treatment plans consisting of photon, electron and spot scanning proton fields. Initially, in-house intensity modulated proton therapy (IMPT) plans are compared to proton plans created by a commercial treatment planning system (TPS). A triple beam radiotherapy (TriB-RT) plan is generated for an exemplary academic case and the dose contributions of the three particle types are investigated. To investigate the dosimetric potential, a TriB-RT plan is compared to an in-house IMPT plan for two clinically motivated cases. Benefits of TriB-RT for a fixed proton beam line with a single proton field are investigated. Results. In-house optimized IMPT are of at least equal or better quality than TPS-generated proton plans, and MC-based optimization shows dosimetric advantages for inhomogeneous situations. Concerning TriB-RT, for the academic case, the resulting plan shows substantial contribution of all particle types. For the clinically motivated case, improved sparing of organs at risk close to the target volume is achieved compared to IMPT (e.g. myelon and brainstem D m a x −37%) at cost of an increased low dose bath (healthy tissue V 10% +22%). In the scenario of a fixed proton beam line, TriB-RT plans are able to compensate the loss in degrees of freedom to substantially improve plan quality compared to a single field proton plan. Conclusion. A novel TPP which simultaneously optimizes photon, electron and proton beams was successfully developed. TriB-RT shows the potential for improved treatment plan quality and is especially promising for cost-effective single-room proton solutions with a fixed beamline in combination with a conventional linac delivering photon and electron fields.

Background Automated brain tumor segmentation methods are computational algorithms that yield tumor delineation from, in this case, multimodal magnetic resonance imaging (MRI). We present an automated segmentation method and its results for resection cavity (RC) in glioblastoma multiforme (GBM) patients using deep learning (DL) technologies. Methods Post-operative, T1w with and without contrast, T2w and fluid attenuated inversion recovery MRI studies of 30 GBM patients were included. Three radiation oncologists manually delineated the RC to obtain a reference segmentation. We developed a DL cavity segmentation method, which utilizes all four MRI sequences and the reference segmentation to learn to perform RC delineations. We evaluated the segmentation method in terms of Dice coefficient (DC) and estimated volume measurements. Results Median DC of the three radiation oncologist were 0.85 (interquartile range [IQR]: 0.08), 0.84 (IQR: 0.07), and 0.86 (IQR: 0.07). The results of the automatic segmentation compared to the three different raters were 0.83 (IQR: 0.14), 0.81 (IQR: 0.12), and 0.81 (IQR: 0.13) which was significantly lower compared to the DC among raters (chi-square = 11.63, p = 0.04). We did not detect a statistically significant difference of the measured RC volumes for the different raters and the automated method (Kruskal-Wallis test: chi-square = 1.46, p = 0.69). The main sources of error were due to signal inhomogeneity and similar intensity patterns between cavity and brain tissues. Conclusions The proposed DL approach yields promising results for automated RC segmentation in this proof of concept study. Compared to human experts, the DC are still subpar.

Purpose: In the current era of MRI-linac radiotherapy, dose optimization with arbitrary dose distributions is a reality. For the first time, we present new and targeted experiments and modeling to aid in evaluating the potential dose improvements offered with an electron beam mode during MRI-linac radiotherapy. Methods: Small collimated (1 cm diameter and 1.5 × 1.5 cm2 square) electron beams (6, 12 and 20 MeV) from a clinical linear accelerator (Varian Clinac 2100C) are incident perpendicular and parallel to the strong and localized magnetic fields (0–0.7 T) generated by a permanent magnet device. Gafchromic EBT3 film is placed inside a slab phantom to measure two-dimensional dose distributions. A benchmarked and comprehensive Monte Carlo model (Geant4) is established to directly compare with experiments. Results: With perpendicular fields a 5% narrowing of the beam FWHM and a 10 mm reduction in the 15% isodose penetration is seen for the 20 MeV beam. In the inline setup the penumbral width is reduced by up to 20%, and a local central dose enhancement of 100% is observed. Monte Carlo simulations are in agreement with the measured dose distributions (2% or 2 mm). Conclusion: A new range of experiments have been performed to offer insight into how an electron beam mode could offer additional choices in MRI-linac radiotherapy. The work extends on historic studies to bring a successful unified experimental and Monte Carlo modeling approach for studying small field electron beam dosimetry inside magnetic fields. The results suggest further work, particularly on the inline magnetic field scenario. © 2020 Associazione Italiana di Fisica Medica, Physica Medica, 78, ISSN:1120-1797, ISSN:1724-191X

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.

Micro computed tomography (µCT) scanners are used to create high-resolution images and to quantify properties of the scanned objects. While modern µCT scanners benefit from the cone beam geometry, they are compromised by scatter radiation. This work aims to develop a Monte Carlo (MC) model of a µCT scanner in order to characterize the scatter radiation in the detector plane. The EGS++ framework with the MC code EGSnrc was used to simulate the particle transport through the main components of the XtremeCT (SCANCO Medical AG, Switzerland). The developed MC model was based on specific information of the manufacturer and was validated against measurements. The primary and the scatter radiation were analyzed and by implementing a dedicated tracing method, the scatter radiation was subdivided into different scatter components. The comparisons of measured and simulated transmission values for different absorber and filter combinations result in a mean difference of 0.2% ± 1.4%, with a maximal local difference of 3.4%. The reconstructed image of the phantom based on measurements agrees well with the image reconstructed using the MC model. The local contribution of scattered radiation is up to 10% of the total radiation in the detector plane and most of the scattered particles result from interactions in the scanned object. The MC simulations show that scatter radiation contains information about the structure of the object. In conclusion, a MC model for a µCT scanner was successfully validated and applied to analyze the characteristics of the scatter radiation for a µCT scanner.

The aim of this work is to develop and investigate an inverse treatment planning process (TPP) for mixed beam radiotherapy (MBRT) capable of performing simultaneous optimization of photon and electron apertures. A simulated annealing based direct aperture optimization (DAO) is implemented to perform simultaneous optimization of photon and electron apertures, both shaped with the photon multileaf collimator (pMLC). Validated beam models are used as input for Monte Carlo dose calculations. Consideration of photon pMLC transmission during DAO and a weight re-optimization of the apertures after deliverable dose calculation are utilized to efficiently reduce the differences between optimized and deliverable dose distributions. The TPP for MBRT is evaluated for an academic situation with a superficial and an enlarged PTV in the depth, a left chest wall case including the internal mammary chain and a squamous cell carcinoma case. Deliverable dose distributions of MBRT plans are compared to those of modulated electron radiotherapy (MERT), photon IMRT and if available to those of clinical VMAT plans. The generated MBRT plans dosimetrically outperform the MERT, photon IMRT and VMAT plans for all investigated situations. For the clinical cases of the left chest wall and the squamous cell carcinoma, the MBRT plans cover the PTV similarly or more homogeneously than the VMAT plans, while OARs are spared considerably better with average reductions of the mean dose to parallel OARs and D 2% to serial OARs by 54% and 26%, respectively. Moreover, the low dose bath expressed as V 10% to normal tissue is substantially reduced by up to 45% compared to the VMAT plans. A TPP for MBRT including simultaneous optimization is successfully implemented and the dosimetric superiority of MBRT plans over MERT, photon IMRT and VMAT plans is demonstrated for academic and clinical situations including superficial targets with and without deep-seated part.

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., Radiation Oncology, 14 (1)

Background Vendor-independent Monte Carlo (MC) dose calculation (IDC) for patient-specific quality assurance of multi-leaf collimator (MLC) based CyberKnife treatments is used to benchmark and validate the commercial MC dose calculation engine for MLC based treatments built into the CyberKnife treatment planning system (Precision MC). Methods The benchmark included dose profiles in water in 15 mm depth and depth dose curves of rectangular MLC shaped fields ranging from 7.6 mm × 7.7 mm to 115.0 mm × 100.1 mm, which were compared between IDC, Precision MC and measurements in terms of dose difference and distance to agreement. Dose distributions of three phantom cases and seven clinical lung cases were calculated using both IDC and Precision MC. The lung PTVs ranged from 14 cm3 to 93 cm3. Quantitative comparison of these dose distributions was performed using dose-volume parameters and 3D gamma analysis with 2% global dose difference and 1 mm distance criteria and a global 10% dose threshold. Time to calculate dose distributions was recorded and efficiency was assessed. Results Absolute dose profiles in 15 mm depth in water showed agreement between Precision MC and IDC within 3.1% or 1 mm. Depth dose curves agreed within 2.3% / 1 mm. For the phantom and clinical lung cases, mean PTV doses differed from − 1.0 to + 2.3% between IDC and Precision MC and gamma passing rates were > =98.1% for all multiple beam treatment plans. For the lung cases, lung V20 agreed within ±1.5%. Calculation times ranged from 2.2 min (for 39 cm3 PTV at 1.0 × 1.0 × 2.5 mm3 native CT resolution) to 8.1 min (93 cm3 at 1.1 × 1.1 × 1.0 mm3), at 2% uncertainty for Precision MC for the 7 examined lung cases and 4–6 h for IDC, which, however, is not optimized for efficiency but used as a gold standard for accuracy. Conclusions Both accuracy and efficiency of Precision MC in the context of MLC based planning for the CyberKnife M6 system were benchmarked against MC based IDC framework. Precision MC is used in clinical practice at our institute.

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 3 mm. 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.

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.

OBJECTIVE In this study, the "Progressive Resolution Optimizer PRO3" (Varian Medical Systems) is compared to the previous version "PRO2" with respect to its potential to improve dose sparing to the organs at risk (OAR) and dose coverage of the PTV for head and neck cancer patients. MATERIALS AND METHODS For eight head and neck cancer patients, volumetric modulated arc therapy (VMAT) treatment plans were generated in this study. All cases have 2-3 phases and the total prescribed dose (PD) was 60-72Gy in the PTV. The study is mainly focused on the phase 1 plans, which all have an identical PD of 54Gy, and complex PTV structures with an overlap to the parotids. Optimization was performed based on planning objectives for the PTV according to ICRU83, and with minimal dose to spinal cord, and parotids outside PTV. In order to assess the quality of the optimization algorithms, an identical set of constraints was used for both, PRO2 and PRO3. The resulting treatment plans were investigated with respect to dose distribution based on the analysis of the dose volume histograms. RESULTS For the phase 1 plans (PD=54Gy) the near maximum dose D2% of the spinal cord, could be minimized to 22±5 Gy with PRO3, as compared to 32±12Gy with PRO2, averaged for all patients. The mean dose to the parotids was also lower in PRO3 plans compared to PRO2, but the differences were less pronounced. A PTV coverage of V95%=97±1% could be reached with PRO3, as compared to 86±5% with PRO2. In clinical routine, these PRO2 plans would require modifications to obtain better PTV coverage at the cost of higher OAR doses. CONCLUSION A comparison between PRO3 and PRO2 optimization algorithms was performed for eight head and neck cancer patients. In general, the quality of VMAT plans for head and neck patients are improved with PRO3 as compared to PRO2. The dose to OARs can be reduced significantly, especially for the spinal cord. These reductions are achieved with better PTV coverage as compared to PRO2. The improved spinal cord sparing offers new opportunities for all types of paraspinal tumors and for re-irradiation of recurrent tumors or second malignancies., Zielsetzung In dieser Studie wird der „Progressive Resolution Optimizer PRO3“ (Varian Medical Systems) mit der Vorgängerversion „PRO2“ gezielt auf sein Potential bei Kopf-Hals-Tumoren untersucht. Im Vordergrund stehen die Dosisabdeckung des PTV und die Schonung der Risikoorgane (OAR). Methoden Für acht Patienten wurden VMAT-Pläne mit beiden Algorithmen erstellt. Die Gesamtdosis im PTV von 60-72 Gy ist in allen Fällen auf 2-3 Phasen verteilt, die nacheinander abgestrahlt werden. Die Studie ist fokussiert auf die erste Bestrahlungsphase, die in allen Fällen eine komplexe PTV-Struktur mit Planungsdosis von 54 Gy hatte und teilweise mit den Parotiden überlappte. Grundlage für die Optimierung sind die ICRU83-Richtlinien für das PTV und das Planungsziel, das Rückenmark und die Parotiden außerhalb des PTV maximal zu schonen. In einem ersten Schritt wurde durch interaktive Planung mit PRO3 ein möglichst optimaler Plan generiert. Die dabei gewählten Grenzwerte („dose constraints“) für PTV und Risikoorgane wurden in einem zweiten Schritt für den Vergleich der Optimierungsalgorithmen verwendet. Die Optimierung wurde also für diese nun fest vorgegebenen Dosisgrenzwerte sowohl mit PRO3 wiederholt als auch mit PRO2 durchgeführt, ohne zusätzliche Interaktion des Planers. Die resultierenden beiden Bestrahlungspläne wurden anhand der Dosis-Volumen-Histogramme verglichen. Ergebnisse Ein markanter Unterschied in Planqualität ergibt sich aus dem Vergleich der Rückenmarkdosis in der ersten Bestrahlungsphase mit 54 Gy Dosis im PTV. Die maximum-nahe Dosis D2% im Rückenmark konnte mit PRO3 auf 22±5 Gy reduziert werden, während die Optimierung mit identischen Randbedingungen bei PRO2 eine Dosis D2% von 32±12 Gy ergab. Auch die mittlere Dosis für die Parotiden war kleiner mit PRO3, allerdings mit einem weniger großen Unterschied. Eine PTV-Abdeckung von V95% = 97±1% konnte mit PRO3 erreicht werden. Mit PRO2 lag diese Abdeckung dagegen nur bei V95% = 86±5%. Die Mehrzahl der so generierten PRO2-Pläne würde in klinischer Routine nicht akzeptiert und müsste durch die Anpassung der Randbedingungen modifiziert werden, so dass die Dosis auf die Risikoorgane entsprechend höher wird. Schlussfolgerung Ein Vergleich der Optimierungsalgorithmen PRO2 und PRO3 wurde bei acht Patienten mit Kopf-Hals-Tumoren durchgeführt. Hierbei wurde mit dem Algorithmus PRO3 eine deutlich bessere Planqualität erzielt als mit PRO2. Die Verbesserung führt zu einer tieferen Belastung der Risikoorgane, mit dem deutlichsten Unterschied beim Rückenmark. Diese Schonung geht nicht zu Lasten der PTV-Abdeckung. Vielmehr wurde bei allen PRO3-Plänen eine bessere PTV-Abdeckung erzielt als mit PRO2. Die verbesserte Möglichkeit der Rückenmarkschonung, auch für komplexe PTV, eröffnet neue Möglichkeiten der Bestrahlung paraspinaler Tumore. Durch die Schonung des Rückenmarks ist zudem eine erneute Bestrahlung bei wiederauftretenden Tumoren möglich.

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, D to the serial OARs is 17% and 37% lower and V 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.

For patients with inoperable liver metastases, intra-operative liver high dose-rate brachytherapy (HDR-BT) is a promising technology enabling delivery of a high radiation dose to the tumor, while sparing healthy tissue. Liver brachytherapy has been described in the literature as safe and effective for the treatment of primary or secondary hepatic malignancies. It is preferred over other ablative techniques for lesions that are either larger than 4 cm or located in close proximity to large vessels or the common bile duct. In contrast to external beam radiation techniques, organ movements do not affect the size of the irradiated volume in intra-operative HDR-BT and new technical solutions exist to support image guidance for intra-operative HDR-BT. We have retrospectively analyzed anonymized CT datasets of 5 patients who underwent open liver surgery (resection and/or ablation) in order to test whether the accuracy of a new image-guidance method specifically adapted for intra-operative HDR-BT is high enough to use it in similar situations and whether patients could potentially benefit from navigation-guided intra-operative needle placement for liver HDR-BT.

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., Journal of Applied Clinical Medical Physics, 19 (1), ISSN:1526-9914

The impact of several physical quantities on the spatial resolution of an X-ray scintillating pixel detector for a micro cone beam CT (µCBCT) is investigated and discussed. The XtremeCT from SCANCO Medical AG was simulated using the EGSnrc/EGS++ Monte Carlo (MC) framework and extensively benchmarked in a previous work. The resolution of the detector was determined by simulating a titanium knife-edge to obtain the edge spread function (ESF) and the modulation transfer function (MTF). Propagation of the scintillation light through the scintillator and its coupling into the fiber optics system was taken into account. The contribution of particles scattered in the main scanner components to the detector signal is very low and does not affect the spatial resolution of the detector. The resolution obtained from the energy deposition in the scintillator without any blurring due to the propagation of the scintillation light into the fiber optics array was 31 µm. By assuming isotropic light propagation in the scintillator, the resolution degraded to 360 µm. A simple light propagation model taking into account the impact of the scintillator's columnar microstructures was developed and compared with the MANTIS Monte Carlo simulation package. By reducing the width of the model's light propagation kernel by a factor of 2 compared to the isotropic case, the detector resolution can be improved to 83 µm, which corresponds well to the measured resolution of 86 µm. The resolution of the detector is limited mainly by the propagation of the scintillation light through the scintillator layer. It offers the greatest potential to improve the resolution of the µCBCT imaging system.

This work aims to develop, implement and validate a Monte Carlo (MC)-based independent dose calculation (IDC) framework to perform patient-specific quality assurance (QA) for multi-leaf collimator (MLC)-based CyberKnife® (Accuray Inc., Sunnyvale, CA) treatment plans. The IDC framework uses an XML-format treatment plan as exported from the treatment planning system (TPS) and DICOM format patient CT data, an MC beam model using phase spaces, CyberKnife MLC beam modifier transport using the EGS++ class library, a beam sampling and coordinate transformation engine and dose scoring using DOSXYZnrc. The framework is validated against dose profiles and depth dose curves of single beams with varying field sizes in a water tank in units of cGy/Monitor Unit and against a 2D dose distribution of a full prostate treatment plan measured with Gafchromic EBT3 (Ashland Advanced Materials, Bridgewater, NJ) film in a homogeneous water-equivalent slab phantom. The film measurement is compared to IDC results by gamma analysis using 2% (global)/2 mm criteria. Further, the dose distribution of the clinical treatment plan in the patient CT is compared to TPS calculation by gamma analysis using the same criteria. Dose profiles from IDC calculation in a homogeneous water phantom agree within 2.3% of the global max dose or 1 mm distance to agreement to measurements for all except the smallest field size. Comparing the film measurement to calculated dose, 99.9% of all voxels pass gamma analysis, comparing dose calculated by the IDC framework to TPS calculated dose for the clinical prostate plan shows 99.0% passing rate. IDC calculated dose is found to be up to 5.6% lower than dose calculated by the TPS in this case near metal fiducial markers. An MC-based modular IDC framework was successfully developed, implemented and validated against measurements and is now available to perform patient-specific QA by IDC.

Standard electron treatments are currently still performed using standard or molded patient-specific cut-outs placed in the electron applicator. Replacing cut-outs and electron applicators with a photon multileaf collimator (pMLC) for electron beam collimation would make standard electron treatments more efficient and would facilitate advanced treatment techniques like modulated electron radiotherapy (MERT) and mixed beam radiotherapy (MBRT). In this work, a multiple source Monte Carlo beam model for pMLC shaped electron beams commissioned at a source-to-surface distance (SSD) of 70 cm is extended for SSDs of up to 100 cm and validated for several Varian treatment units with field sizes typically used for standard electron treatments. Measurements and dose calculations agree generally within 3% of the maximal dose or 2 mm distance to agreement. To evaluate the dosimetric consequences of using pMLC collimated electron beams for standard electron treatments, pMLC-based and cut-out-based treatment plans are created for a left and a right breast boost, a sternum, a testis and a parotid gland case. The treatment plans consist of a single electron field, either alone (1E) or in combination with two 3D conformal tangential photon fields (1E2X). For each case, a pMLC plan with similar treatment plan quality in terms of dose homogeneity to the target and absolute mean dose values to the organs at risk (OARs) compared to a cut-out plan is found. The absolute mean dose to an OAR is slightly increased for pMLC-based compared to cut-out-based 1E plans if the OAR is located laterally close to the target with respect to beam direction, or if a 6 MeV electron beam is used at an extended SSD. In conclusion, treatment plans using cut-out collimation can be replaced by plans of similar treatment plan quality using pMLC collimation with accurately calculated dose distributions.

Purpose: Advanced radiation therapy requires highly sophisticated dose calculation algorithms such as finite element based Boltzmann solvers or Monte Carlo(MC) methods. MC is commonly accepted as the golden standard method for dose calculation in high energy treatments and thus it is used for benchmarking other algorithms. In this work the quality of dose distribution calculated using the Boltzmann solver based AcurosXB algorithm within Eclipse (Varian Medical Systems) is investigated for volumetric modulated arc treatment (VMAT) plans involving the high definition MLC (HD120 MLC) by comparing doses with the validated Swiss Monte Carlo Plan (SMCP).Materials & Methods: Within SMCP and Eclipse using AcurosXB, 10 VMAT H&N patient plans and corresponding verification plans were recalculated using fixed MUs. In SMCP, radiation transport and dose calculation were performed using VMC++ with a statistical uncertainty of 1%. The voxel size was 2.5 mm for SMCP and AcurosXB and the same material composition data was used for CT conversion. Dose volume histograms (DVH) were used in order to quantify the difference between the dose distributions of the patient plans. In addition, calculated verification plans were compared with measurements carried out with the Delta4 system (Scandidos) by using the gamma evaluation with 3%/3 mm criteria of points having a dose larger than 20% of isocenter dose.Results: DVHs for the patient plans showed good agreement between SMCP and AcurosXB calculations. Overall AcurosXB lead to an underestimation of the median dose values by about 1%. For measured total dose distributions of the verification plans on average 98.6% and 99.0% of the points fullfil the gamma criteria for the dose calculated using AcurosXB and SMCP, respectively. Conclusions: Resulting AcurosXB dose distributions for VMAT H&N plans involving a HD120 MLC are in good agreement with calculated SMCP dose distributions. Conflict of Interest: This work was supported by Varian Medical Systems