Two new DOSXYZnrc sources for 4D Monte Carlo simulations of continuously variable beam configurations, with applications to RapidArc, VMAT, TomoTherapy and CyberKnife

Monte Carlo Modeling of Dynamic Tumor Tracking on a Gimbaled Linear Accelerator

Purpose and Aim

The Vero4DRT (Brainlab AG) linear accelerator is capable of dynamic tumor tracking (DTT) by panning/tilting the radiation beam to follow respiratory-induced tumor motion in real time. In this study, the panning/tilting motion is modeled in Monte Carlo (MC) for quality assurance (QA) of four-dimensional (4D) dose distributions created within the treatment planning system (TPS).

Materials and Methods

Step-and-shoot intensity-modulated radiation therapy plans were optimized for 10 previously treated liver patients. These plans were recalculated on multiple phases of a 4D computed tomography (4DCT) scan using MC while modeling panning/tilting. The dose distributions on each phase were accumulated to create a respiratory-weighted 4D dose distribution. Differences between the TPS and MC modeled doses were examined.

Results

On average, 4D dose calculations in MC showed the maximum dose of an organ at risk (OAR) to be 10% greater than the TPS' three-dimensional dose calculation (collapsed cone [CC] convolution algorithm) predicted. MC's 4D dose calculations showed that 6 out of 24 OARs could exceed their specified dose limits, and calculated their maximum dose to be 4% higher on average (up to 13%) than the TPS' 4D dose calculations. Dose differences between MC and the TPS were greatest in the beam penumbra region.

Conclusion

Modeling panning/tilting for DTT has been successfully modeled with MC and is a useful tool to QA respiratory-correlated 4D dose distributions. The dose differences between the TPS and MC calculations highlight the importance of using 4D MC to confirm the safety of OAR doses before DTT treatments.

A new DOSXYZnrc method for Monte Carlo simulations of 4D dose distributions

The purpose of this study is to present a novel method for generating Monte Carlo 4D dose distributions in a single DOSXYZnrc simulation. During a standard simulation, individual energy deposition events are summed up to generate a 3D dose distribution and their associated temporal information is discarded. This means that in order to determine dose distributions as a function of time, separate simulations would have to be run for each interval of interest. Consequently, it has not been clinically feasible until now to routinely perform Monte Carlo simulations of dose rate, time-resolved dose accumulation, or electronic portal imaging devices (EPID) cine-mode images for volumetric modulated arc therapy (VMAT) plans. To overcome this limitation, we modified DOSXYZnrc and defined new input and output variables that allow a time-like parameter associated with each particle history to be binned in a user-defined manner. Under the new code version, computation times are the same as for a standard simulation, and the time-integrated 4D dose is identical to the standard 3D dose. We present a comparison of scintillator measurements and Monte Carlo simulations for dose rate during a VMAT beam delivery, a study of dose rate in a VMAT total body irradiation plan, and simulations of transit (through-patient) EPID cine-mode images.

Monte Carlo methods for device simulations in radiation therapy

Monte Carlo (MC) simulations play an important role in radiotherapy, especially as a method to evaluate physical properties that are either impossible or difficult to measure. For example, MC simulations (MCSs) are used to aid in the design of radiotherapy devices or to understand their properties. The aim of this article is to review the MC method for device simulations in radiation therapy. After a brief history of the MC method and popular codes in medical physics, we review applications of the MC method to model treatment heads for neutral and charged particle radiation therapy as well as specific in-room devices for imaging and therapy purposes. We conclude by discussing the impact that MCSs had in this field and the role of MC in future device design.

Ion chamber and film-based quality assurance of mixed electron-photon radiation therapy

PURPOSE

In previous work, we demonstrated that mixed electron-photon radiation therapy (MBRT) produces treatment plans with improved normal tissue sparing and similar target coverage, when compared to photon-only plans. The purpose of this work was to validate the MBRT delivery process on a Varian TrueBeam accelerator and laying the groundwork for a patient-specific quality assurance (QA) protocol based on ion chamber point measurements and 2D film measurements.

METHODS

MC beam models used to calculate the MBRT dose distributions of each modality (photons/electrons) were validated with a single-angle beam MBRT treatment plan delivered on a slab of Solid Water phantom with a film positioned at a depth of 2 cm. The measured film absorbed dose was compared to the calculated dose. To validate clinical deliveries, a polymethyl methacrylate (PMMA) cylinder was machined and holes were made to fit an ionization chamber. A complex MBRT plan involving a photon arc and three electron delivery angles was created with the aim of reproducing a clinically realistic dose distribution in typical soft tissue sarcoma tumours of the extremities. The treatment plan was delivered on the PMMA cylinder. Point measurements were taken with an Exradin A1SL chamber at two nominal depths: 1.4 cm and 2.1 cm. The plan was also delivered on a second identical phantom with an insert at 2 cm depth, where a film was placed. An existing EGSnrc user-code, SPRRZnrc, was modified to calculate the stopping power ratios between any materials in the same voxelized geometry used for dose calculation purposes. This modified code, called SPRXYZnrc, was used to calculate a correction factor, k MBRT , accounting for the differences in electron fluence spectrum at the measurement point compared to that at reference conditions. The uncertainty associated with neglecting potential ionization chamber fluence perturbation correction factors using this approach was estimated.

RESULTS

The film measurement from the Solid Water phantom treatment plan was in good agreement with the simulated dose distribution, with a gamma pass rate of 96.1% for a 3%/2 mm criteria. For the PMMA phantom delivery, for the same gamma criteria, the pass rate was 97.3%. The ion chamber measurements of the total delivered dose agreed with the MC-simulated dose within 2.1%. The beam quality correction factors amounted to, at most, a 4% correction on the ion chamber measurement. However, individual contribution of low electron energies proved difficult to precisely measure due to their steep dose gradients, with disagreements of up to 28% ± 15% at 2.1 cm depth (6 MeV). Ion chamber measurement procedure of electron beams was achieved in less than 5 min, and the entire validation process including phantom setup was performed in less than 30 min.

CONCLUSION

The agreement between measured and simulated MBRT doses indicates that the dose distributions obtained from the MBRT treatment planning algorithm are realistically achievable. The SPRXYZnrc MC code allowed for convenient calculations of k MBRT simultaneously with the dose distributions, laying the groundwork for patient-specific QA protocol practical for clinical use. Further investigation is needed to establish the accuracy of our ionization chamber correction factors k MBRT calculations at low electron energies.

Monte Carlo as quality control tool of stereotactic body radiation therapy treatment plans

PURPOSE/OBJECTIVE

The objective of this study was to verify the accuracy of treatment plans of stereotactic body radiation therapy (SBRT) and to verify the feasibility of the use of Monte Carlo (MC) as quality control (QC) on a daily basis.

MATERIAL/METHODS

Using EGSnrc, a MC model of Agility™ linear accelerator was created. Various measurements (Percentage depth dose (PDD), Profiles and Output factors) were done for different fields sizes from 1x1 up to 40x40 (cm²). An iterative model optimization was performed to achieve adequate parameters of MC simulation. 40 SBRT patient's dosimetry plans were calculated by Monaco™ 3.1.1. CT images, RT-STRUCT and RT-PLAN files from Monaco™ being used as input for Moderato MC code. Finally, dose volume histogram (DVH) and paired t-tests for each contour were used for dosimetry comparison of the Monaco™ and MC.

RESULTS

Validation of MC model was successful, as <2% difference comparing to measurements for all field's sizes. The main energy of electron source incident on the target was 5.8 MeV, and the full width at half maximum (FWHM) of Gaussian electron source were 0.09 and 0.2 (cm) in X and Y directions, respectively. For 40 treatment plan comparisons, the minimum absolute difference of mean dose of planning treatment planning (PTV) was 0.1% while the maximum was 6.3%. The minimum absolute difference of Max dose of PTV was 0.2% while the maximum was 8.1%.

CONCLUSION

SBRT treatment plans of Monaco agreed with MC results. It possible to use MC for treatment plans verifications as independent QC tool.

GPU-accelerated Monte Carlo simulation of MV-CBCT

Monte Carlo simulation (MCS) is one of the most accurate computation methods for dose calculation and image formation in radiation therapy. However, the high computational complexity and long execution time of MCS limits its broad use. In this paper, we present a novel strategy to accelerate MCS using a graphic processing unit (GPU), and we demonstrate the application in mega-voltage (MV) cone-beam computed tomography (CBCT) simulation. A new framework that generates a series of MV projections from a single simulation run is designed specifically for MV-CBCT acquisition. A Geant4-based GPU code for photon simulation is incorporated into the framework for the simulation of photon transport through a phantom volume. The FastEPID method, which accelerates the simulation of MV images, is modified and integrated into the framework. The proposed GPU-based simulation strategy was tested for its accuracy and efficiency in a Catphan 604 phantom and an anthropomorphic pelvis phantom with beam energies at 2.5 MV, 6 MV, and 6 MV FFF. In all cases, the proposed GPU-based simulation demonstrated great simulation accuracy and excellent agreement with measurement and CPU-based simulation in terms of reconstructed image qualities. The MV-CBCT simulation was accelerated by factors of roughly 900-2300 using an NVIDIA Tesla V100 GPU card against a 2.5 GHz AMD Opteron™ Processor 6380.

Monte Carlo simulation of 6-MV dynamic wave VMAT deliveries by Vero4DRT linear accelerator using EGSnrc moving sources

The commissioning and benchmark of a Monte Carlo (MC) model of the 6‐MV Brainlab‐Mitsubishi Vero4DRT linear accelerator for the purpose of quality assurance of clinical dynamic wave arc (DWA) treatment plans is reported. Open‐source MC applications based on EGSnrc particle transport codes are used to simulate the medical linear accelerator head components. Complex radiotherapy irradiations can be simulated in a single MC run using a shared library format combined with BEAMnrc “source20.” Electron energy tuning is achieved by comparing measured vs simulated percentage depth doses (PDDs) for MLC‐defined field sizes in a water phantom. Electron spot size tuning is achieved by comparing measured and simulated inplane and crossplane beam profiles. DWA treatment plans generated from RayStation (RaySearch) treatment planning system (TPS) are simulated on voxelized (2.5 mm³) patient CT datasets. Planning target volume (PTV) and organs at risk (OAR) dose–volume histograms (DVHs) are compared to TPS‐calculated doses for clinically deliverable dynamic volumetric modulated arc therapy (VMAT) trajectories. MC simulations with an electron beam energy of 5.9 MeV and spot size FWHM of 1.9 mm had the closest agreement with measurement. DWA beam deliveries simulated on patient CT datasets results in DVH agreement with TPS‐calculated doses. PTV coverage agreed within 0.1% and OAR max doses (to 0.035 cc volume) agreed within 1 Gy. This MC model can be used as an independent dose calculation from the TPS and as a quality assurance tool for complex, dynamic radiotherapy treatment deliveries. Full patient CT treatment simulations are performed in a single Monte Carlo run in 23 min. Simulations are run in parallel using the Condor High‐Throughput Computing software¹ on a cluster of eight servers. Each server has two physical processors (Intel Xeon CPU E5‐2650 0 @2.00 GHz), with 8 cores per CPU and two threads per core for 256 calculation nodes.

Monte Carlo calculations of radiotherapy dose in "homogeneous" anatomy

Given the substantial literature on the use of Monte Carlo (MC) simulations to verify treatment planning system (TPS) calculations of radiotherapy dose in heterogeneous regions, such as head and neck and lung, this study investigated the potential value of running MC simulations of radiotherapy treatments of nominally homogeneous pelvic anatomy. A pre-existing in-house MC job submission and analysis system, built around BEAMnrc and DOSXYZnrc, was used to evaluate the dosimetric accuracy of a sample of 12 pelvic volumetric arc therapy (VMAT) treatments, planned using the Varian Eclipse TPS, where dose was calculated with both the Analytical Anisotropic Algorithm (AAA) and the Acuros (AXB) algorithm. In-house TADA (Treatment And Dose Assessor) software was used to evaluate treatment plan complexity, in terms of the small aperture score (SAS), modulation index (MI) and a novel exposed leaf score (ELS/ELA). Results showed that the TPS generally achieved closer agreement with the MC dose distribution when treatments were planned for smaller (single-organ) targets rather than larger targets that included nodes or metastases. Analysis of these MC results with reference to the complexity metrics indicated that while AXB was useful for reducing dosimetric uncertainties associated with density heterogeneity, the residual TPS dose calculation uncertainties resulted from treatment plan complexity and TPS model simplicity. The results of this study demonstrate the value of using MC methods to recalculate and check the dose calculations provided by commercial radiotherapy TPSs, even when the treated anatomy is assumed to be comparatively homogeneous, such as in the pelvic region.

Clinical implementation of a Monte Carlo based independent TPS dose checking system

The increase in complexity of treatment plans over time through modalities such as intensity modulated radiotherapy (IMRT) and volumetric modulated arc therapy (VMAT) has often not been met with an increase in capability of the secondary dose calculation checking systems typically used to verify the treatment planning system. Monte Carlo (MC) codes such as EGSnrc have become easily available and are capable of performing calculations of highly complex radiotherapy treatments. This educational note demonstrates a method for implementing and using a fully automated system for performing and analysing full MC calculations of conformal, IMRT and VMAT radiotherapy plans. Example calculations were based on BEAMnrc/DOSXYZnrc and are performed automatically after either uploading exported plan DICOM data through a Python-based web interface, or exporting DICOM data to a monitored network location. This note demonstrates how completed MC calculations can then be analysed using an automatically generated dose point comparison report, or easily re-imported back into the treatment planning system. Agreement between the TPS and MC calculation was an improvement on agreement between RadCalc and the TPS, with differences ranging from 1.2 to 5.5% between RadCalc and the treatment planning system (TPS), and 0.1-1.7% between MC and TPS. Comparison of the dose-volume histogram (DVH) parameters [Formula: see text], [Formula: see text], [Formula: see text], and [Formula: see text] for the example VMAT plans showed agreement for the mean planning target volume dose within [Formula: see text], [Formula: see text] and [Formula: see text] generally within [Formula: see text] with the exception of a brain case, and [Formula: see text] within [Formula: see text]. Overall, this note provides a demonstration of a system that has been integrated well into existing clinical workflow, and has been shown to be a valuable additional tool in the secondary checking of treatment plan calculations.

Validation of 4D Monte Carlo dose calculations using a programmable deformable lung phantom

PURPOSE

To validate the accuracy of 4D Monte Carlo (4DMC) simulations to calculate dose deliveries to a deforming anatomy in the presence of realistic respiratory motion traces. A previously developed deformable lung phantom comprising an elastic tumor was modified to enable programming of arbitrary motion profiles. 4D simulations of the dose delivered to the phantom were compared with the measurements.

METHODS

The deformable lung phantom moving with irregular breathing patterns was irradiated using static and VMAT beam deliveries. Using the RADPOS 4D dosimetry system, point doses were measured inside and outside the tumor. Dose profiles were acquired using films along the motion path of the tumor (S-I). In addition to dose measurements, RADPOS was used to record the motion of the tumor during dose deliveries. Dose measurements were then compared against 4DMC simulations with EGSnrc/4DdefDOSXYZnrc using the recorded tumor motion.

RESULTS

The agreements between dose profiles from measurements and simulations were determined to be within 2%/2 mm. Point dose agreements were within 2σ of experimental and/or positional/dose reading uncertainties. 4DMC simulations were shown to accurately predict the sensitivity of delivered dose to the starting phase of breathing motions. We have demonstrated that our 4DMC method, combined with RADPOS, can accurately simulate realistic dose deliveries to a deforming anatomy moving with realistic breathing traces. This 4DMC tool has the potential to be used as a quality assurance tool to verify treatments involving respiratory motion. Adaptive treatment delivery is another area that may benefit from the potential of this 4DMC tool.

Large-scale dosimetric assessment of Monte Carlo recalculated doses for lung robotic stereotactic body radiation therapy

Owing to its short computation time and simplicity, the Ray-Tracing algorithm (RAT) has long been used to calculate dose distributions for the CyberKnife system. However, it is known that RAT fails to fully account for tissue heterogeneity and is therefore inaccurate in the lung. The aim of this study is to make a dosimetric assessment of 219 non-small cell lung cancer CyberKnife plans by recalculating their dose distributions using an independent Monte Carlo (MC) method. For plans initially calculated by RAT without heterogeneity corrections, target coverage was found to be significantly compromised when considering MC doses. Only 35.4% of plans were found to comply to their prescription doses. If the normal tissue dose limits were respected in the treatment planning dose, the MC recalculated dose did not exceed these limits in over 97% of the plans. Comparison of RAT and recalculated-MC doses confirmed the overestimation of RAT doses observed in previous studies. An inverse correlation between the RAT/MC dose ratio and the target size was also found to be statistically significant (p<10^-4), consistent with other studies. In addition, the inaccuracy and variability in target coverage incurred from dose calculations using RAT without heterogeneity corrections was demonstrated. On average, no clinically relevant differences were observed between MC-calculated dose-to-water and dose-to-medium for all tissues investigated (⩽1%). Patients receiving a dose D_95% larger than 119 Gy in EQD2₁₀ (or ≈52 Gy in 3 fractions) as recalculated by MC were observed to have significantly superior loco-regional progression-free survival rates (p=0.02) with a hazard ratio of 3.45 (95%CI: 1.14-10.5).

Monte Carlo dose verification for a single-isocenter VMAT plan in multiple brain metastases

The purpose of this study was to verify the accuracy of dose calculation algorithms of a treatment planning system for a single-isocenter volumetric modulated arc therapy (VMAT) plan in multiple brain metastases, by comparing the dose distributions of treatment planning system with those of Monte Carlo (MC) simulations. We used a multitarget phantom containing 9 acrylic balls with a diameter of 15.9 mm inside a Lucy phantom measuring 17 × 17 × 17 cm³. Seven VMAT plans were created using the multitarget phantom: 1 multitarget plan (MTP) and 6 single target plans (STP). Three of the STP plans had a large jaw field setting, almost equivalent to that of the MTP, while the other plans had a jaw field setting fitted to each planning target volume. The isocenter for all VMAT plans was set to the center of the phantom. The VMAT dose distributions were calculated using the analytical anisotropic algorithm (AAA) and were also recalculated through Acuros XB (AXB) and MC simulations under the same irradiation conditions. The AAA and AXB methods tended to overestimate dosage compared with the MC method in the MTP and in STPs with large jaw field settings. The dose distribution in single-isocenter VMAT plans for multiple brain metastases was influenced by jaw field settings. Finally, we concluded that MC-VMAT dose calculations are useful for 3D dose verification of single-isocenter VMAT plans for multiple brain metastases.

EGSnrc application for IMRT planning

The aim of this study was to describe a detailed instruction of intensity modulated radiotherapy (IMRT) planning simulation using BEAMnrc-DOSXYZnrc code system (EGSnrc package) and present a new graphical user interface based on MATLAB code (The MathWorks) to combine more than one. 3ddose file which were obtained from the IMRT plan. This study was performed in four phases: the commissioning of Varian Clinac iX6 MV, the simulation of IMRT planning in EGSnrc, the creation of in-house VDOSE GUI, and the analysis of the isodose contour and dose volume histogram (DVH) curve from several beam angles. The plan paramaters in sequence and control point files were extracted from the planning data in Tan Tock Seng Hospital Singapore (multileaf collimator (MLC) leaf positions - bank A and bank B, gantry angles, coordinate of isocenters, and MU indexes). VDOSE GUI which was created in this study can display the distribution dose curve in each slice and beam angle. Dose distributions from various MLC settings and beam angles yield different dose distributions even though they used the same number of simulated particles. This was due to the differences in the MLC leaf openings in every field. The value of the relative dose error between the two dose ditributions for "body" was 51.23 %. The Monte Carlo (MC) data was normalized with the maximum dose but the analytical anisotropic algorithm (AAA) data was normalized by the dose in the isocenter. In this study, we have presented a Monte Carlo simulation framework for IMRT dose calculation using DOSXYZnrc source 21. Further studies are needed in conducting IMRT simulations using EGSnrc to minimize the different dose error and dose volume histogram deviation.

Electron modulated arc therapy (EMAT) using photon MLC for postmastectomy chest wall treatment I: Monte Carlo-based dosimetric characterizations

PURPOSE

To study the dosimetric properties of electron arc beams delivered by photon-beam multi-leaf collimators (pMLC) in electron modulated arc therapy (EMAT) for postmastectomy chest wall treatments.

METHODS

Using the Monte Carlo method, we simulated a 2100EX Varian linear accelerator and verified the beam models in a water tank. Dosimetric characterizations were performed on cylindrical water phantoms of elliptical bases with various field sizes, arc ranges and source-to-surface distances (SSDs) for 6, 9 and 12 MeV beam energy.

RESULTS

The arc beam has a higher bremsstrahlung dose than the static beam at the isocenter due to crossfire, but choosing a field size greater than 5 cm effectively reduces the bremsstrahlung dose. The depths of the 90% maximum dose located at 1.7, 2.8 and 4.1 cm for 6, 9 and 12 MeV, respectively, are similar to those of the static beams and independent of the field size and arc range.

CONCLUSION

Based on the study, we recommend using the 5 cm field width for electron arc beams considering both bremsstrahlung dose at the isocenter and the arc profile penumbra. To ensure sufficient PTV edge coverage, we recommend a field length extension of at least 4 cm from PTV's edge for all beam energies and an arc extension of around 7°, 5°, and 5° for beam energies 6, 9, and 12 MeV, respectively. These dosimetric characterizations are the basis of pMLC-delivered EMAT treatment planning for postmastectomy chest wall patients.

Radiological tissue equivalence of deformable silicone-based chemical radiation dosimeters (FlexyDos3D)

FlexyDos3D, a silicone‐based chemical radiation dosimeter, has great potential to serve as a three‐dimensional (3D) deformable dosimetric tool to verify complex dose distributions delivered by modern radiotherapy techniques. To facilitate its clinical application, its radiological tissue needs to be clarified. In this study we investigated its tissue‐equivalence in comparison with water and Solid Water (RMI457). We found that its effective and mean atomic numbers were 40% and 20% higher and the total interaction probabilities for kV x‐ray photons were larger than those of water respectively. To assess the influence of its over‐response to kV photons, its HU value was measured by kV computed tomography (CT) and was found higher than all the soft‐tissue substitutes. When applied for dose calculation without correction, this effect led to an 8% overestimation in electron density via HU‐value mapping and 0.65% underestimation in target dose. Furthermore, depth dose curves (PDDs) and off‐axis ratios (profiles) at various beam conditions as well as the dose distribution of a full‐arc VMAT plan in FlexyDos3D and reference materials were simulated by Monte Carlo, where the results showed great agreement. As indicated, FlexyDos3D exhibits excellent radiological water‐equivalence for clinical MV x‐ray dosimetry, while its nonwater‐equivalent effect for low energy x‐ray dosimetry requires necessary correction. The key findings of this study provide pertinent reference for further FlexyDos3D characterization research.

Evaluating small field dosimetry with the Acuros XB (AXB) and analytical anisotropic algorithm (AAA) dose calculation algorithms in the eclipse treatment planning system

Background:

An increasing number of external beam treatment modalities including intensity modulated radiation therapy, volumetric modulated arc therapy (VMAT) and stereotactic radiosurgery uses very small fields for treatment planning and delivery. However, there are major challenges in small photon field dosimetry, due to the partial occlusion of the direct photon beam source’s view from the measurement point, lack of lateral charged particle equilibrium, steep dose-rate gradient and volume averaging effect of the detector response and variation of the energy fluence in the lateral direction of the beam. Therefore, experimental measurements of dosimetric parameters such as percent depth doses (PDDs), beam profiles and relative output factors (ROFs) for small fields continue to be a challenge.

Materials and Methods:

In this study, we used a homogeneous water phantom and the heterogeneous anthropomorphic stereotactic end-to-end verification (STEEV) head phantom for all dose measurements and calculations. PDDs, lateral dose profiles and ROFs were calculated in the Eclipse Treatment Planning System version 13·6 using the Acuros XB (AXB) and the analytical anisotropic algorithms (AAAs) in a homogenous water phantom. Monte Carlo (MC) simulations and measurements using the Exradin W1 Scintillator were also accomplished for four photon energies: 6 MV, 6FFF, 10 MV and 10FFF. Two VMAT treatment plans were generated for two different targets: one located in the brain and the other in the neck (close to the trachea) in the head phantom (CIRS, Norfolk, VA, USA). A Varian Truebeam linear accelerator (Varian, Palo Alto, CA, USA) was used for all treatment deliveries. Calculated results with AXB and AAA were compared with MC simulations and measurements.

Results:

The average difference of PDDs between W1 Exradin Scintillator measurements and MC simulations, AAA and AXB algorithm calculations were 1·2, 2·4 and 3·2%, respectively, for all field sizes and energies. AXB and AAA showed differences in ROF of about 0·3 and 2·9%, respectively, compared with W1 Exradin Scintillator measured values. For the target located in the brain in the head phantom, the average dose difference between W1 Exradin Scintillator and the MC simulations, AAA and AXB were 0·2, 3·2 and 2·7%, respectively, for all field sizes. Similarly, for the target located in the neck, the respective dose differences were 3·8, 5·7 and 3·5%.

Conclusion:

In this study, we compared dosimetric parameters such as PDD, beam profile and ROFs in water phantom and isocenter point dose measurements in an anthropomorphic head phantom representing a patient. We observed that measurements using the W1 Exradin scintillator agreed well with MC simulations and can be used efficiently for dosimetric parameters such as PDDs and dose profiles and patient-specific quality assurance measurements for small fields. In both homogenous and heterogeneous media, the AXB algorithm dose prediction agrees well with MC and measurements and was found to be superior to the AAA algorithm.

Dependence of volume dose indices on dose calculation algorithms for VMAT-SBRT plans for peripheral lung tumor

The purpose of this study was to investigate the dependence of volume dose indices on dose calculation algorithms for volumetric modulated arc therapy (VMAT) for stereotactic body radiotherapy (SBRT) plans to treat peripheral lung tumors by comparing them with those of Monte Carlo (MC) calculations. VMAT-SBRT plans for peripheral lung tumors were created using the Eclipse treatment planning system (TPS) for 24 patients with nonsmall cell lung cancer. VMAT dose distributions for gross tumor volume (GTV), internal target volume (ITV), and planning target volume (PTV) were calculated using the analytical anisotropic algorithm (AAA), the Acuros XB (AXB) algorithm, and a MC algorithm. VMAT dose distributions of the 3 algorithms were compared using their volume dose indices from dose volume histograms (DVHs), a dose difference map, and 3-dimensional gamma analysis. The DVHs for GTV and ITV from AAA, AXB, and MC were in good agreement. The difference between the ITV and PTV volume dose indices from AAA and MC increased as D₉₈, D₉₅, D₈₀, D₅₀, and D₂. In particular, the difference between D₉₈ for PTV from AAA and MC was up to 48%. A >5% difference between D₉₅ for PTV from AAA and MC was 11 patients, but only 2 patients for ITV. The volume dose indices for AXB were near those of MC. AAA tended to overestimate the PTV volume dose indices compared to AXB and MC. Thus, we propose that the volume dose indices for the ITV be used because they are independent of dose calculation algorithms.

Estimating the uncertainty of calculated out-of-field organ dose from a commercial treatment planning system

Therapeutic radiation to cancer patients is accompanied by unintended radiation to organs outside the treatment field. It is known that the model-based dose algorithm has limitation in calculating the out-of-field doses. This study evaluated the out-of-field dose calculated by the Varian Eclipse treatment planning system (v.11 with AAA algorithm) in realistic treatment plans with the goal of estimating the uncertainties of calculated organ doses. Photon beam phase-space files for TrueBeam linear accelerator were provided by Varian. These were used as incident sources in EGSnrc Monte Carlo simulations of radiation transport through the downstream jaws and MLC. Dynamic movements of the MLC leaves were fully modeled based on treatment plans using IMRT or VMAT techniques. The Monte Carlo calculated out-of-field doses were then compared with those calculated by Eclipse. The dose comparisons were performed for different beam energies and treatment sites, including head-and-neck, lung, and pelvis. For 6 MV (FF/FFF), 10 MV (FF/FFF), and 15 MV (FF) beams, Eclipse underestimated out-of-field local doses by 30%-50% compared with Monte Carlo calculations when the local dose was <1% of prescribed dose. The accuracy of out-of-field dose calculations using Eclipse is improved when collimator jaws were set at the smallest possible aperture for MLC openings. The Eclipse system consistently underestimates out-of-field dose by a factor of 2 for all beam energies studied at the local dose level of less than 1% of prescribed dose. These findings are useful in providing information on the uncertainties of out-of-field organ doses calculated by Eclipse treatment planning system.

A fast jaw-tracking model for VMAT and IMRT Monte Carlo simulations

Modern radiotherapy techniques involve routine use of volumetric arc therapy (VMAT) and intensity modulated radiotherapy (IMRT) with jaw‐tracking – dynamic motion of the secondary collimators (jaws) in tandem with multi‐leaf collimators (MLCs). These modalities require accurate dose calculations for the purposes of treatment planning and dose verification. Monte Carlo (MC) methods for radiotherapy dose calculation are widely accepted as capable of achieving high accuracy. This paper presents an efficiency‐enhancement method for secondary collimator modeling, presented in the context of a tool for MC‐based dose second checks. The model constitutes an accuracy trade‐off in the source model for the sake of efficiency enhancement, but maintains the advantages of MC transport in patient heterogeneities. The secondary collimator model is called Flat‐Absorbing‐Jaw‐Tracking (FAJT). Transmission through and scatter from the secondary collimators is neglected, and jaws are modeled as perfectly absorbing planes. To couple the motion of secondary collimators with MLCs for jaw‐tracking, the FAJT model was built into the VCU‐MLC model. Gamma‐index analysis of the dose distributions from FAJT against the full BEAMnrc MC simulations showed over 99% pass rate for a range of open fields, two clinical IMRT, and one VMAT treatment plan, for 2%/2 mm criteria above 10%. Using FAJT, the simulation speed of the secondary collimators for open fields increased by a factor of 237, 1489, and 1395 for 4 × 4, 10 × 10, and 30 × 30 cm², respectively. In general, clinically oriented simulation times are reduced from “hours” to “minutes” on identical hardware. Results for nine representative clinical cases (seven with jaw‐tracking) are presented. The average 2%/2 mm γ‐test success rate above the 80% isodose was 96.8% when tested against the EPIDose electronic portal image‐based dose reconstruction method and 97.3% against the Eclipse analytical anisotropic algorithm.

An automated Monte Carlo QC system for volumetric modulated arc therapy: Possibilities and challenges

PURPOSE

To develop and implement an automated Monte Carlo (MC) system for patient specific VMAT quality control in a patient geometry that generates treatment planning system (TPS) compliant DICOM objects and includes a module for 3D analysis of dose deviations. Also, the aims were to recommend diagnose specific tolerance criteria and an evaluation procedure.

METHODS

The EGSnrc code package formed the basis for development of the MC system. The workflow consists of a number of modules connected to a TPS by means of manual DICOM exports and imports which were executed sequentially without user interaction. DVH comparison was performed in the TPS. In addition, MC- and TPS dose distributions were analysed by applying the normalized dose difference (NDD) formalism. NDD failure maps and a pass rate for a certain threshold were obtained. 170 clinical plans (prostate, thorax, head-and-neck and gynecological) were selected for analysis.

RESULTS

Agreement within 1.5% was found between clinical- and MC data for the mean dose to the target volumes and within 3% for parameters more sensitive to the shape of the DVH e.g. D_98% PTV. Regarding the NDD analysis, tolerance criteria 2%/3 mm were established for prostate plans and 3%/3 mm for the rest of the cases.

CONCLUSIONS

An automated MC system was developed and implemented. Evaluation procedure is recommended with NDD-analysis as a first step. For pass rate < 95%, the evaluation continues with comparison of DVH parameters. For deviations larger than 2%, a visual inspection of the clinical- and MC dose distributions is performed.

A simple technique to improve calculated skin dose accuracy in a commercial treatment planning system

Radiation dermatitis during radiotherapy is correlated with skin dose and is a common clinical problem for head and neck and thoracic cancer patients. Therefore, accurate prediction of skin dose during treatment planning is clinically important. The objective of this study is to evaluate the accuracy of skin dose calculated by a commercial treatment planning system (TPS). We evaluated the accuracy of skin dose calculations by the anisotropic analytical algorithm (AAA) implemented in Varian Eclipse (V.11) system. Skin dose is calculated as mean dose to a contoured structure of 0.5 cm thickness from the surface. The EGSnrc Monte Carlo (MC) simulations are utilized for the evaluation. The 6, 10 and 15 MV photon beams investigated are from a Varian TrueBeam linear accelerator. The accuracy of the MC dose calculations was validated by phantom measurements with optically stimulated luminescence detectors. The calculation accuracy of patient skin doses is studied by using CT based radiotherapy treatment plans including 3D conformal, static gantry IMRT, and VMAT treatment techniques. Results show the Varian Eclipse system underestimates skin doses by up to 14% of prescription dose for the patients studied when external body contour starts at the patient's skin. The external body contour is used in a treatment planning system to calculate dose distributions. The calculation accuracy of skin dose with Eclipse can be considerably improved to within 4% of target dose by extending the external body contour by 1 to 2 cm from the patient's skin. Dose delivered to deeper target volumes or organs at risk are not affected. Although Eclipse treatment planning system has its limitations in predicting patient skin dose, this study shows the calculation accuracy can be considerably improved to an acceptable level by extending the external body contour without affecting the dose calculation accuracy to the treatment target and internal organs at risk. This is achieved by moving the calculation entry point away from the skin.

3D VMAT Verification Based on Monte Carlo Log File Simulation with Experimental Feedback from Film Dosimetry

A model based on a specific phantom, called QuAArC, has been designed for the evaluation of planning and verification systems of complex radiotherapy treatments, such as volumetric modulated arc therapy (VMAT). This model uses the high accuracy provided by the Monte Carlo (MC) simulation of log files and allows the experimental feedback from the high spatial resolution of films hosted in QuAArC. This cylindrical phantom was specifically designed to host films rolled at different radial distances able to take into account the entrance fluence and the 3D dose distribution. Ionization chamber measurements are also included in the feedback process for absolute dose considerations. In this way, automated MC simulation of treatment log files is implemented to calculate the actual delivery geometries, while the monitor units are experimentally adjusted to reconstruct the dose-volume histogram (DVH) on the patient CT. Prostate and head and neck clinical cases, previously planned with Monaco and Pinnacle treatment planning systems and verified with two different commercial systems (Delta4 and COMPASS), were selected in order to test operational feasibility of the proposed model. The proper operation of the feedback procedure was proved through the achieved high agreement between reconstructed dose distributions and the film measurements (global gamma passing rates > 90% for the 2%/2 mm criteria). The necessary discretization level of the log file for dose calculation and the potential mismatching between calculated control points and detection grid in the verification process were discussed. Besides the effect of dose calculation accuracy of the analytic algorithm implemented in treatment planning systems for a dynamic technique, it was discussed the importance of the detection density level and its location in VMAT specific phantom to obtain a more reliable DVH in the patient CT. The proposed model also showed enough robustness and efficiency to be considered as a pre-treatment VMAT verification system.

Quantification of dose differences between two versions of Acuros XB algorithm compared to Monte Carlo simulations--the effect on clinical patient treatment planning

A commercialized implementation of linear Boltzmann transport equation solver, the Acuros XB algorithm (AXB), represents a class of most advanced type 'c' photon radiotherapy dose calculation algorithms. The purpose of the study was to quantify the effects of the modifications implemented in the more recent version 11 of the AXB (AXB11) compared to the first commercial implementation, version 10 of the AXB (AXB10), in various anatomical regions in clinical treatment planning. Both versions of the AXB were part of Varian's Eclipse clinical treatment planning system and treatment plans for 10 patients were created using intensity-modulated radiotherapy (IMRT) and volumetric-modulated arc radiotherapy (VMAT). The plans were first created with the AXB10 and then recalculated with the AXB11 and full Monte Carlo (MC) simulations. Considering the full MC simulations as reference, a DVH analysis for gross tumor and planning target volumes (GTV and PTV) and organs at risk was performed, and also 3D gamma agreement index (GAI) values within a 15% isodose region and for the PTV were determined. Although differences up to 12% in DVH analysis were seen between the MC simulations and the AXB, based on the results of this study no general conclusion can be drawn that the modifications made in the AXB11 compared to the AXB10 would imply that the dose calculation accuracy of the AXB10 would be inferior to the AXB11 in the clinical patient treatment planning. The only clear improvement with the AXB11 over the AXB10 is the dose calculation accuracy in air cavities. In general, no large deviations are present in the DVH analysis results between the two versions of the algorithm, and the results of 3D gamma analysis do not favor one or the other. Thus it may be concluded that the results of the comprehensive studies assessing the accuracy of the AXB10 may be extended to the AXB11.

The accuracy of Acuros XB algorithm for radiation beams traversing a metallic hip implant - comparison with measurements and Monte Carlo calculations

In this study, the clinical benefit of the improved accuracy of the Acuros XB (AXB) algorithm, implemented in a commercial radiotherapy treatment planning system (TPS), Varian Eclipse, was demonstrated with beams traversing a high‐Z material. This is also the first study assessing the accuracy of the AXB algorithm applying volumetric modulated arc therapy (VMAT) technique compared to full Monte Carlo (MC) simulations. In the first phase the AXB algorithm was benchmarked against point dosimetry, film dosimetry, and full MC calculation in a water‐filled anthropometric phantom with a unilateral hip implant. Also the validity of the full MC calculation used as reference method was demonstrated. The dose calculations were performed both in original computed tomography (CT) dataset, which included artifacts, and in corrected CT dataset, where constant Hounsfield unit (HU) value assignment for all the materials was made. In the second phase, a clinical treatment plan was prepared for a prostate cancer patient with a unilateral hip implant. The plan applied a hybrid VMAT technique that included partial arcs that avoided passing through the implant and static beams traversing the implant. Ultimately, the AXB‐calculated dose distribution was compared to the recalculation by the full MC simulation to assess the accuracy of the AXB algorithm in clinical setting. A recalculation with the anisotropic analytical algorithm (AAA) was also performed to quantify the benefit of the improved dose calculation accuracy of type ‘c’ algorithm (AXB) over type ‘b’ algorithm (AAA). The agreement between the AXB algorithm and the full MC model was very good inside and in the vicinity of the implant and elsewhere, which verifies the accuracy of the AXB algorithm for patient plans with beams traversing through high‐Z material, whereas the AAA produced larger discrepancies.

PACS numbers: 87.55.‐x, 87.55.D‐, 87.55.K‐, 87.55.kd, 87.55.Qr

MCTP system model based on linear programming optimization of apertures obtained from sequencing patient image data maps

PURPOSE

The authors present a hybrid direct multileaf collimator (MLC) aperture optimization model exclusively based on sequencing of patient imaging data to be implemented on a Monte Carlo treatment planning system (MC-TPS) to allow the explicit radiation transport simulation of advanced radiotherapy treatments with optimal results in efficient times for clinical practice.

METHODS

The planning system (called CARMEN) is a full MC-TPS, controlled through aMATLAB interface, which is based on the sequencing of a novel map, called "biophysical" map, which is generated from enhanced image data of patients to achieve a set of segments actually deliverable. In order to reduce the required computation time, the conventional fluence map has been replaced by the biophysical map which is sequenced to provide direct apertures that will later be weighted by means of an optimization algorithm based on linear programming. A ray-casting algorithm throughout the patient CT assembles information about the found structures, the mass thickness crossed, as well as PET values. Data are recorded to generate a biophysical map for each gantry angle. These maps are the input files for a home-made sequencer developed to take into account the interactions of photons and electrons with the MLC. For each linac (Axesse of Elekta and Primus of Siemens) and energy beam studied (6, 9, 12, 15 MeV and 6 MV), phase space files were simulated with the EGSnrc/BEAMnrc code. The dose calculation in patient was carried out with the BEAMDOSE code. This code is a modified version of EGSnrc/DOSXYZnrc able to calculate the beamlet dose in order to combine them with different weights during the optimization process.

RESULTS

Three complex radiotherapy treatments were selected to check the reliability of CARMEN in situations where the MC calculation can offer an added value: A head-and-neck case (Case I) with three targets delineated on PET/CT images and a demanding dose-escalation; a partial breast irradiation case (Case II) solved with photon and electron modulated beams (IMRT + MERT); and a prostatic bed case (Case III) with a pronounced concave-shaped PTV by using volumetric modulated arc therapy. In the three cases, the required target prescription doses and constraints on organs at risk were fulfilled in a short enough time to allow routine clinical implementation. The quality assurance protocol followed to check CARMEN system showed a high agreement with the experimental measurements.

CONCLUSIONS

A Monte Carlo treatment planning model exclusively based on maps performed from patient imaging data has been presented. The sequencing of these maps allows obtaining deliverable apertures which are weighted for modulation under a linear programming formulation. The model is able to solve complex radiotherapy treatments with high accuracy in an efficient computation time.

The accuracy of the out-of-field dose calculations using a model based algorithm in a commercial treatment planning system

The out-of-field dose can be clinically important as it relates to the dose of the organ-at-risk, although the accuracy of its calculation in commercial radiotherapy treatment planning systems (TPSs) receives less attention. This study evaluates the uncertainties of out-of-field dose calculated with a model based dose calculation algorithm, anisotropic analytical algorithm (AAA), implemented in a commercial radiotherapy TPS, Varian Eclipse V10, by using Monte Carlo (MC) simulations, in which the entire accelerator head is modeled including the multi-leaf collimators. The MC calculated out-of-field doses were validated by experimental measurements. The dose calculations were performed in a water phantom as well as CT based patient geometries and both static and highly modulated intensity-modulated radiation therapy (IMRT) fields were evaluated. We compared the calculated out-of-field doses, defined as lower than 5% of the prescription dose, in four H&N cancer patients and two lung cancer patients treated with volumetric modulated arc therapy (VMAT) and IMRT techniques. The results show that the discrepancy of calculated out-of-field dose profiles between AAA and the MC depends on the depth and is generally less than 1% for in water phantom comparisons and in CT based patient dose calculations for static field and IMRT. In cases of VMAT plans, the difference between AAA and MC is <0.5%. The clinical impact resulting from the error on the calculated organ doses were analyzed by using dose-volume histograms. Although the AAA algorithm significantly underestimated the out-of-field doses, the clinical impact on the calculated organ doses in out-of-field regions may not be significant in practice due to very low out-of-field doses relative to the target dose.

Monte Carlo modeling of HD120 multileaf collimator on Varian TrueBeam linear accelerator for verification of 6X and 6X FFF VMAT SABR treatment plans

A Monte Carlo (MC) validation of the vendor-supplied Varian TrueBeam 6 MV flattened (6X) phase-space file and the first implementation of the Siebers-Keall MC MLC model as applied to the HD120 MLC (for 6X flat and 6X flattening filter-free (6X FFF) beams) are described. The MC model is validated in the context of VMAT patient-specific quality assurance. The Monte Carlo commissioning process involves: 1) validating the calculated open-field percentage depth doses (PDDs), profiles, and output factors (OF), 2) adapting the Siebers-Keall MLC model to match the new HD120-MLC geometry and material composition, 3) determining the absolute dose conversion factor for the MC calculation, and 4) validating this entire linac/MLC in the context of dose calculation verification for clinical VMAT plans. MC PDDs for the 6X beams agree with the measured data to within 2.0% for field sizes ranging from 2 × 2 to 40 × 40 cm2. Measured and MC profiles show agreement in the 50% field width and the 80%-20% penumbra region to within 1.3 mm for all square field sizes. MC OFs for the 2 to 40 cm2 square fields agree with measurement to within 1.6%. Verification of VMAT SABR lung, liver, and vertebra plans demonstrate that measured and MC ion chamber doses agree within 0.6% for the 6X beam and within 2.0% for the 6X FFF beam. A 3D gamma factor analysis demonstrates that for the 6X beam, &gt; 99% of voxels meet the pass criteria (3%/3 mm). For the 6X FFF beam, &gt; 94% of voxels meet this criteria. The TrueBeam accelerator delivering 6X and 6X FFF beams with the HD120 MLC can be modeled in Monte Carlo to provide an independent 3D dose calculation for clinical VMAT plans. This quality assurance tool has been used clinically to verify over 140 6X and 16 6X FFF TrueBeam treatment plans.

Monte Carlo comparison of superficial dose between flattening filter free and flattened beams

This study investigates the superficial dose from FFF beams in comparison with the conventional flattened ones using a Monte Carlo (MC) method. Published phase-space files which incorporated real geometry of a TrueBeam accelerator were used for the dose calculation in phantom and clinical cases. The photon fluence on the central axis is 3 times that of a flattened beam for a 6 MV FFF beam and 5 times for a 10 MV beam. The mean energy across the field in air at the phantom surface is 0.92-0.95 MeV for the 6 MV FFF beam and 1.18-1.30 MeV for the corresponding flattened beam. At 10 MV, the values are 1.52-1.72 and 2.15-2.87 MeV for the FFF and flattened beams, respectively. The phantom dose at the depth of 1 mm in the 6 MV FFF beam is 6% ± 2.5% (of the maximum dose) higher compared to the flattened beam for a 25 × 25 cm(2) field and 14.6% ± 1.9% for the 2 × 2 cm(2) field. For the 10 MV beam, the corresponding differences are 3.4% ± 1.5% and 10.7% ± 0.6%. The skin dose difference at selected points on the patient's surface between the plans using FFF and flattened beams in the head-and-neck case was 6.5% ± 2.3% (1SD), and for the breast case it was 6.4% ± 2.3%. The Monte Carlo simulations showed that due to the lower mean energy in the FFF beam, the clinical superficial dose is higher without the flattening filter compared to the flattened beam.

A Monte Carlo approach to validation of FFF VMAT treatment plans for the TrueBeam linac

PURPOSE

To commission and benchmark a vendor-supplied (Varian Medical Systems) Monte Carlo phase-space data for the 6 MV flattening filter free (FFF) energy mode on a TrueBeam linear accelerator for the purpose of quality assurance of clinical volumetric modulated arc therapy (VMAT) treatment plans. A method for rendering the phase-space data compatible with BEAMnrc/DOSXYZnrc simulation software package is presented.

METHODS

Monte Carlo (MC) simulations were performed to benchmark the TrueBeam 6 MV FFF phase space data that have been released by the Varian MC Research team. The simulations to benchmark the phase space data were done in three steps. First, the original phase space which was created on a cylindrical surface was converted into a format that was compatible with BEAMnrc. Second, BEAMnrc was used to create field size specific phase spaces located underneath the jaws. Third, doses were calculated with DOSXYZnrc in a water phantom for fields ranging from 1 × 1 to 40 × 40 cm(2). Calculated percent depth doses (PDD), transverse profiles, and output factors were compared with measurements for all the fields simulated. After completing the benchmarking study, three stereotactic body radiotherapy (SBRT) VMAT plans created with the Eclipse treatment planning system (TPS) were calculated with Monte Carlo. Ion chamber and film measurements were also performed on these plans. 3D gamma analysis was used to compare Monte Carlo calculation with TPS calculations and with film measurement.

RESULTS

For the benchmarking study, MC calculated and measured values agreed within 1% and 1.5% for PDDs and in-field transverse profiles, respectively, for field sizes >1 × 1 cm(2). Agreements in the 80%-20% penumbra widths were better than 2 mm for all the fields that were compared. With the exception of the 1 × 1 cm(2) field, the agreement between measured and calculated output factors was within 1%. It is of note that excellent agreement in output factors for all field sizes including highly asymmetric fields was achieved without accounting for backscatter into the beam monitor chamber. For the SBRT VMAT plans, the agreement between Monte Carlo and ion chamber point dose measurements was within 1%. Excellent agreement between Monte Carlo, treatment planning system and Gafchromic film dose distribution was observed with over 99% of the points in the high dose volume passing the 3%, 3 mm gamma test.

CONCLUSIONS

The authors have presented a method for making the Varian IAEA compliant 6 MV FFF phase space file of the TrueBeam linac compatible with BEAMnrc/DOSXYZnrc. After benchmarking the modified phase space against measurement, they have demonstrated its potential for use in MC based quality assurance of complex delivery techniques.

Forcing lateral electron disequilibrium to spare lung tissue: a novel technique for stereotactic body radiation therapy of lung cancer

Stereotactic body radiation therapy (SBRT) has quickly become a preferred treatment option for early-stage lung cancer patients who are ineligible for surgery. This technique uses tightly conformed megavoltage (MV) x-ray beams to irradiate a tumour with ablative doses in only a few treatment fractions. Small high energy x-ray fields can cause lateral electron disequilibrium (LED) to occur within low density media, which can reduce tumour dose. These dose effects may be challenging to predict using analytic dose calculation algorithms, especially at higher beam energies. As a result, previous authors have suggested using low energy photons (<10 MV) and larger fields (>5 × 5 cm(2)) for lung cancer patients to avoid the negative dosimetric effects of LED. In this work, we propose a new form of SBRT, described as LED-optimized SBRT (LED-SBRT), which utilizes radiotherapy (RT) parameters designed to cause LED to advantage. It will be shown that LED-SBRT creates enhanced dose gradients at the tumour/lung interface, which can be used to manipulate tumour dose, and/or normal lung dose. To demonstrate the potential benefits of LED-SBRT, the DOSXYZnrc (National Research Council of Canada, Ottawa, ON) Monte Carlo (MC) software was used to calculate dose within a cylindrical phantom and a typical lung patient. 6 MV or 18 MV x-ray fields were focused onto a small tumour volume (diameter ∼1 cm). For the phantom, square fields of 1 × 1 cm(2), 3 × 3 cm(2), or 5 × 5 cm(2) were applied. However, in the patient, 3 × 1 cm(2), 3 × 2 cm(2), 3 × 2.5 cm(2), or 3 × 3 cm(2) field sizes were used in simulations to assure target coverage in the superior-inferior direction. To mimic a 180° SBRT arc in the (symmetric) phantom, a single beam profile was calculated, rotated, and beams were summed at 1° segments to accumulate an arc dose distribution. For the patient, a 360° arc was modelled with 36 equally weighted (and spaced) fields focused on the tumour centre. A planning target volume (PTV) was generated by considering the extent of tumour motion over the patient's breathing cycle and set-up uncertainties. All patient dose results were normalized such that at least 95% of the PTV received at least 54 Gy (i.e. D95 = 54 Gy). Further, we introduce 'LED maps' as a novel clinical tool to compare the magnitude of LED resulting from the various SBRT arc plans. Results from the phantom simulation suggest that the best lung sparing occurred for RT parameters that cause severe LED. For equal tumour dose coverage, normal lung dose (2 cm outside the target region) was reduced from 92% to 23%, comparing results between the 18 MV (5 × 5 cm(2)) and 18 MV (1 × 1 cm(2)) arc simulations. In addition to reduced lung dose for the 18 MV (1 × 1 cm(2)) arc, maximal tumour dose increased beyond 125%. Thus, LED can create steep dose gradients to spare normal lung, while increasing tumour dose levels (if desired). In the patient simulation, a LED-optimized arc plan was designed using either 18 MV (3 × 1 cm(2)) or 6 MV (3 × 3cm(2)) beams. Both plans met the D95 dose coverage requirement for the target. However, the LED-optimized plan increased the maximum, mean, and minimum dose within the PTV by as much as 80 Gy, 11 Gy, and 3 Gy, respectively. Despite increased tumour dose levels, the 18 MV (3 × 1 cm(2)) arc plan improved or maintained the V20, V5, and mean lung dose metrics compared to the 6 MV (3 × 3 cm(2)) simulation. We conclude that LED-SBRT has the potential to increase dose gradients, and dose levels within a small lung tumour. The magnitude of tumour dose increase or lung sparing can be optimized through manipulation of RT parameters (e.g. beam energy and field size).