Study of efficiency in five-field and field-by-field intensity modulated radiation therapy (IMRT) plan using DOSXYZnrc Monte Carlo code

Implementation of a modern treatment technique, such as IMRT, has been improved. In line with that, Monte Carlo (MC) simulations of this technique require the ability of complex beam configurations modelling with respect to the patient. The source 20 DOSXYZnrc with the dynamic and step and shoot technique can be used to simulate the modality. However, they have a different process to obtain the dose distribution in a certain phantom. This study aimed to compare the simulation efficiency and isodose dose distribution in a water phantom from various beam angles and multileaf collimator (MLC) positions in an IMRT plan using source 20. The 30 × 30 × 30 cm³ phantom was irradiated by Varian Clinac iX10MV photon beam with various field sizes from 2 × 2 to 6 × 6 cm² using some beam angles 5°, 30°, 90°, 180°, and 300° and maintaining the source to surface distance (SSD) of 100 cm. The field-by-field and five-field methods were used to obtain the 3-dimensional (3D) dose distribution. The dose distribution of these methods was compared using the gamma index, DVH analysis, and simulation efficiency. Higher efficiency is better because it implies that it takes less time to reach a given uncertainty. The implementation of source 20 has been validated, with similar results, with validated source in DOSXYZnrc. The identical 3D-dimensions dose distributions using source 20 for dynamic and step and shoot were observed. Two simulations used the same number of histories with the statistical uncertainty of less than 3%. The step and shoot technique was more efficient than the dynamic simulation.     

Comparative analysis of several detectors for the measurement of radiation transmission and leakage from a multileaf collimator

The multileaf collimator (MLC) is the standard device used to shape radiation beams for 3-d conformal and intensity-modulated radiation therapy (IMRT). Due to the inherent properties of MLC, there is a small amount of radiation transmitted through the leaves, called radiation transmission (RT). Accurate measurements of this radiation are required to commission and validate IMRT-capable treatment planning systems because this radiation may impact the dosimetry of IMRT-calculated dose distributions. This work compares several detectors in the measurement of RT for a micro-multileaf collimation system. The results show that there are statistically significant differences in the measured RT values between detectors from 3.5 to 12.5% for the same MLC model and less than 0.2% relative to the isocentre dose for an open reference field. However, although small in magnitude, these differences may impact the dosimetry of IMRT treatment planning by up to 1.78 Gy to the healthy tissue surrounding the target for a treatment of 60 Gy in 30 fractions. By the later, these differences must be included as a source of uncertainty in IMRT dose delivery. Also, it must be established which detector offers the most reliable results in the measurement of the RT by using Monte Carlo simulation methods.     

Monte Carlo dose verification of VMAT treatment plans using Elekta Agility 160-leaf MLC

In this study, we verified volumetric modulated arc therapy (VMAT) plans in an Elekta Synergy system with an integrated Agility 160-leaf multileaf collimator (MLC) by comparing them with Monte Carlo (MC)-calculated dose distributions using the AAPM TG-119 structure sets. The head configuration of the linear accelerator with the integrated MLC was simulated with the EGSnrc/BEAMnrc code. Firstly, the dosimetric properties of the MLC were evaluated with the MC technique and film measurements. Next, VMAT plans were created with the Pinnacle³ treatment planning system (TPS) for four regions in the AAPM TG-119 structures. They were then verified by comparing them with MC-calculated dose distributions using dose volume histograms (DVHs) and three-dimensional (3D) gamma analysis. The MC simulations for the Agility MLC dosimetric properties were in acceptable agreement with measurements. TPS-VMAT plans using TG-119 structure sets agreed with MC dose distributions within 2% in the comparison of D₉₅ in planning target volumes (PTVs) evaluated from DVHs. In contrast, higher dose regions such as D₂₀, D₁₀, and D₅ in PTVs for TPS tended to be smaller than MC values. This tendency was particularly noticeable for mock head and neck with complicated structures. In 3D gamma analysis, the passing rates with 3%/3mm criteria in PTVs were ≥99%, except for mock head and neck (89.5%). All passing rates for organs at risk (OARs) were in acceptable agreement of >96%. It is useful to verify dose distributions of PTVs and OARs in TPS-VMAT plans by using MC dose calculations and 3D gamma analysis.     

Incorrect dosimetric leaf separation in IMRT and VMAT treatment planning: Clinical impact and correlation with pretreatment quality assurance

PurposeDynamic treatment planning algorithms use a dosimetric leaf separation (DLS) parameter to model the multi-leaf collimator (MLC) characteristics. Here, we quantify the dosimetric impact of an incorrect DLS parameter and investigate whether common pretreatment quality assurance (QA) methods can detect this effect.Methods16 treatment plans with intensity modulated radiation therapy (IMRT) or volumetric modulated arc therapy (VMAT) technique for multiple treatment sites were calculated with a correct and incorrect setting of the DLS, corresponding to a MLC gap difference of 0.5 mm. Pretreatment verification QA was performed with a bi-planar diode array phantom and the electronic portal imaging device (EPID). Measurements were compared to the correct and incorrect planned doses using gamma evaluation with both global (G) and local (L) normalization. Correlation, specificity and sensitivity between the dose volume histogram (DVH) points for the planning target volume (PTV) and the gamma passing rates were calculated.ResultsThe change in PTV and organs at risk DVH parameters were 0.4–4.1%. Good correlation (>0.83) between the PTV_mean dose deviation and measured gamma passing rates was observed. Optimal gamma settings with 3%L/3 mm (per beam and composite plan) and 3%G/2 mm (composite plan) for the diode array phantom and 2%G/2 mm (composite plan) for the EPID system were found. Global normalization and per beam ROC analysis of the diode array phantom showed an area under the curve <0.6.ConclusionsA DLS error can worsen pretreatment QA using gamma analysis with reasonable credibility for the composite plan. A low detectability was demonstrated for a 3%G/3 mm per beam gamma setting.     

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

Purpose/objectiveThe 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/methodsUsing 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.ResultsValidation 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%.ConclusionSBRT treatment plans of Monaco agreed with MC results. It possible to use MC for treatment plans verifications as independent QC tool.     

Dosimetric and Monte Carlo verification of jaws-only IMRT plans calculated by the Collapsed Cone Convolution algorithm for head and neck cancers

AimThe aim of this study is to verify the Prowess Panther jaws-only intensity modulated radiation therapy (JO-IMRT) treatment planning (TP) by comparing the TP dose distributions for head-and-neck (H&N) cancer with the ones simulated by Monte Carlo (MC).BackgroundTo date, dose distributions planned using JO-IMRT for H&N patients were found superior to the corresponding three-dimensional conformal radiotherapy (3D-CRT) plans. Dosimetry of the JO-IMRT plans were also experimentally verified using an ionization chamber, MapCHECK 2, and Octavius 4D and good agreements were shown.Materials and methodsDose distributions of 15 JO-IMRT plans of nasopharyngeal patients were recalculated using the EGSnrc Monte Carlo code. The clinical photon beams were simulated using the BEAMnrc. The absorbed dose to patients treated by fixed-field IMRT was computed using the DOSXYZnrc. The simulated dose distributions were then compared with the ones calculated by the Collapsed Cone Convolution (CCC) algorithm on the TPS, using the relative dose error comparison and the gamma index using global methods implemented in PTW-VeriSoft with 3%/3 mm, 2%/2 mm, 1%/1 mm criteria.ResultsThere is a good agreement between the MC and TPS dose. The average gamma passing rates were 93.3 ± 3.1%, 92.8 ± 3.2%, 92.4 ± 3.4% based on the 3%/3 mm, 2%/2 mm, 1%/1 mm criteria, respectively.ConclusionsAccording to the results, it is concluded that the CCC algorithm was adequate for most of the IMRT H&N cases where the target was not immediately adjacent to the critical structures.     

Intelligence-guided beam angle optimization in treatment planning of intensity-modulated radiation therapy

An intelligence guided approach based on fuzzy inference system (FIS) was proposed to automate beam angle optimization in treatment planning of intensity-modulated radiation therapy (IMRT). The model of FIS is built on inference rules in describing the relationship between dose quality of IMRT plan and irradiated region of anatomical structure. Dose quality of IMRT plan is quantified by the difference between calculated and constraint doses of the anatomical structures in an IMRT plan. Irradiated region of anatomical structure is characterized by the metric, covered region of interest, which is the region of an anatomical structure under radiation field while beam’s eye-view is conform to target volume. Initially, an IMRT plan is created with a single beam. The dose difference is calculated for the input of FIS and the output of FIS is obtained with processing of fuzzy inference. Later, a set of candidate beams is generated for replacing the current beam. This process continues until no candidate beams is found. Then the next beam is added to the IMRT plan and optimized in the same way as the previous beam. The new beam keeps adding to the IMRT plan until the allowed beam number is reached. Two spinal cases were investigated in this study. The preliminary results show that dose quality of IMRT plans achieved by this approach is better than those achieved by the default approach with equally spaced beam setting. It is effective to find the optimal beam combination of IMRT plan with the intelligence-guided approach.     

Correlation between gamma passing rate and complexity of IMRT plan due to MLC position errors

PurposeThis study evaluates the correlation between the susceptibility of the γ passing rate of IMRT plans to the multi-leaf collimator (MLC) position errors and a quantitative plan complexity metric.MethodsTwenty patients were selected for this study. For each patient, two IMRT plans were generated using sliding window and step-&-shoot techniques, respectively. Modulation complexity score (MCS) was calculated for all IMRT plans, and symmetric MLC leaf bank errors, ranging from 0.3 mm to 1 mm, were introduced. Original and modified plans were delivered using Varian’s Clinac iX. The obtained dose distribution using ArcCHECK was then compared with the TPS calculated dose distribution of the original plans. 3D gamma analysis was performed for each verification with passing criteria of 2%/2 mm. The γ passing rate decreasing gradient were calculated to evaluate relationship between variation of γ passing rate due to MLC errors and complexity.ResultsA linear regression analysis was applied between γ gradient and complexity, and the results showed a linear correlation (R² = 0.81 and 0.82 for open and closed MLC error types, respectively) indicating the more complex plans are more susceptible to MLC leaf bank errors. Meanwhile, correlation of re-normalized γ passing rate and complexity for all errors scenarios also presented a strong correlation (r > 0.75).ConclusionThe statistics results revealed variation relationship of dosimetry robust of plans with various complexities to MLC errors. Our results also suggested that the observed susceptibility is independent of the delivery techniques.     

Multiple-source models for electron beams of a medical linear accelerator using BEAMDP computer code

AimThe aim of this work was to develop multiple-source models for electron beams of the NEPTUN 10PC medical linear accelerator using the BEAMDP computer code.BackgroundOne of the most accurate techniques of radiotherapy dose calculation is the Monte Carlo (MC) simulation of radiation transport, which requires detailed information of the beam in the form of a phase-space file. The computing time required to simulate the beam data and obtain phase-space files from a clinical accelerator is significant. Calculation of dose distributions using multiple-source models is an alternative method to phase-space data as direct input to the dose calculation system.Materials and methodsMonte Carlo simulation of accelerator head was done in which a record was kept of the particle phase-space regarding the details of the particle history. Multiple-source models were built from the phase-space files of Monte Carlo simulations. These simplified beam models were used to generate Monte Carlo dose calculations and to compare those calculations with phase-space data for electron beams.ResultsComparison of the measured and calculated dose distributions using the phase-space files and multiple-source models for three electron beam energies showed that the measured and calculated values match well each other throughout the curves.ConclusionIt was found that dose distributions calculated using both the multiple-source models and the phase-space data agree within 1.3%, demonstrating that the models can be used for dosimetry research purposes and dose calculations in radiotherapy.     

Comparative Planning of Flattening-Filter-Free and Flat Beam IMRT for Hypopharynx Cancer as a Function of Beam and Segment Number

Although highly conformal dose distributions can be achieved by IMRT planning, this often requires a large number of segments or beams, resulting in increased treatment times. While flattening-filter-free beams offer a higher dose rate, even more segments may be required to create homogeneous target coverage. Therefore, it is worthwhile to systematically investigate the dependence of plan quality on gantry angles and number of segments for flat vs. FFF beams in IMRT planning. For the practical example of hypopharynx cancer, we present a planning study of flat vs. FFF beams using three different configurations of gantry angles and different segment numbers. The two beams are very similar in physical properties, and are hence well-suited for comparative planning. Starting with a set of plans of equal quality for flat and FFF beams, we assess how far the number of segments can be reduced before the plan quality is markedly compromised, and compare monitor units and treatment times for the resulting plans. As long as a sufficiently large number of segments is permitted, all planning scenarios give good results, independently of gantry angles and flat or FFF beams. For smaller numbers of segments, plan quality decreases both for flat and FFF energies; this effect is stronger for fewer gantry angles and for FFF beams. For low segment numbers, FFF plans are generally worse than the corresponding flat beam plans, but they are less sensitive to a decrease in segment number if many gantry angles are used (18 beams); in this case the quality of flat and FFF plans remains comparable even for few segments.     

Efficiency of a novel non-monotonic segmented leaf sequence delivery of Varian MLC for non-split IMRT fields

AimDevelopment of bidirectional non-monotonic segmented leaf sequence (NSLS) MLC delivery technique compatible with Varian MLC for non-split IMRT fields reducing total monitor units (TotalMU) and the number of segments (NS) simultaneously and assessment of its efficiency using a plan scoring index (PSI).Materials and methodsThe optimal fluence of IMRT plans of ten patients of lung carcinoma, calculated using Eclipse TPS version 11.0 (Varian Medical Systems, Palo Alto, CA, USA), was used to generate the segmented MLC fields using our newly developed equally spaced (ES) reducing level and NSLS algorithms in MATLAB® version 2011b for 6–10 intensity levels. These MLC fields were imported into the plans with the same field setup and the final dose was recalculated. The results were compared with those of commercially available multiple static segments (MSS) leaf motion calculation (LMC) algorithm and few previously published algorithms. Plan scoring index (PSI) and degree of modulation (DoM) was calculated to compare the quality of different plans for the same patient.ResultsThe average differences in TotalMU and NS with respect to MSS algorithm are −3.80% and −14.28% for the NSLS algorithm, respectively. The calculated average PSI and DoM is 0.75, 2.51 and 0.91, 2.41 for the MSS and NSLS algorithms, respectively.ConclusionsIMRT plans generated using the NSLS algorithm resulted in the best PSI, DoM values among all the leaf sequencing algorithms. Our proposed NSLS algorithm allows bidirectional delivery in Varian medical linear accelerator which is not commercially available. NSLS algorithm is efficient in reducing the TotalMU and NS with equivalent plan quality as that of MSS.     

Determination of initial electron parameters by means of Monte Carlo simulations for the Siemens Artiste Linac 6 MV photon beam

AimIn this study, we investigated initial electron parameters of Siemens Artiste Linac with 6 MV photon beam using the Monte Carlo method.BackgroundIt is essential to define all the characteristics of initial electrons hitting the target, i.e. mean energy and full width of half maximum (FWHM) of the spatial distribution intensity, which is needed to run Monte Carlo simulations. The Monte Carlo is the most accurate method for simulation of radiotherapy treatments.Materials and methodsLinac head geometry was modeled using the BEAMnrc code. The phase space files were used as input file to DOSXYZnrc simulation to determine the dose distribution in a water phantom. We obtained percent depth dose curves and the lateral dose profile. All the results were obtained at 100 cm of SSD and for a 10 × 10 cm² field.ResultsWe concluded that there existed a good conformity between Monte Carlo simulation and measurement data when we used electron mean energy of 6.3 MeV and 0.30 cm FWHM value as initial parameters. We observed that FWHM values had very little effect on PDD and we found that the electron mean energy and FWHM values affected the lateral dose profile. However, these effects are between tolerance values.ConclusionsThe initial parameters especially depend on components of a linac head. The phase space file which was obtained from Monte Carlo Simulation for a linac can be used as calculation of scattering, MLC leakage, to compare dose distribution on patients and in various studies.     

Objective function based ranking method for selection of optimal beam angles in IMRT

We propose a novel method for the selection of optimal beam angles in Intensity Modulated Radiation Therapy (IMRT). The proposed approach uses an objective function based metric called “target-to-critical organ objective function ratio” to find out the optimal gantry angles. The beams are ranked based on this metric and are accordingly chosen for IMRT optimization. We have used the Pinnacle TPS (Philips Medical System V 16.2) for performing the IMRT optimization. In order to validate our approach, we have applied it in four clinical cases: Head and Neck, Lung, Abdomen and Prostate. Basically, for all clinical cases, two set of plans were created with same clinical objectives, namely Equal angle plan (EA Plan) and Suitable angle Plan (SA Plan). In the EA plans, the beam angles were placed in an equiangular manner starting from the gantry angle of 0°. In the corresponding SA plans, the beam angles were decided using the guidance provided by the algorithm. The reduction in OAR mean dose and max dose obtained in SA plans is about 3 to 16% and 3 to 15% respectively depending upon the treatment site while obtaining equal target coverage as compared to their EA counterparts. It takes approximately 15–25 min to find the optimal beam angles. The results obtained from the clinical cases indicate that the plan quality is considerably improved when the beam angles are optimized using the proposed method.     

Influence of the Presence of Tissue Expanders on Energy Deposition for Post-Mastectomy Radiotherapy

An increasing number of studies have shown that post-mastectomy radiotherapy presents benefits associated with the patients survival and a significant fraction of the treated patients makes use of tissue expanders for breast reconstruction. Some models of tissue expanders have a magnetic disk on their surface that constitutes heterogeneity in the radiation field, which can affect the dose distribution during the radiotherapy treatment. In this study, the influence of a metallic heterogeneity positioned in a breast tissue expander was evaluated by means of Monte Carlo simulations using the MCNPX code and using Eclipse treatment planning system. Deposited energy values were calculated in structures which have clinical importance for the treatment. Additionally, the effect in the absorbed energy due to backscattering and attenuation of the incident beam caused by the heterogeneity, as well as due to the expansion of the prosthesis, was evaluated in target structures for a 6 MV photon beam by simulations. The dose distributions for a breast treatment were calculated using a convolution/superposition algorithm from the Eclipse treatment planning system. When compared with the smallest breast expander volume, underdosage of 7% was found for the largest volume of breast implant, in the case of frontal irradiation of the chest wall, by Monte Carlo simulations. No significant changes were found in dose distributions for the presence of the heterogeneity during the treatment planning of irradiation with an opposed pair of beams. Even considering the limitation of the treatment planning system, the results obtained with its use confirm those ones found by Monte Carlo simulations for a tangent beam irradiation. The presence of a heterogeneity didńt alters the dose distributions on treatment structures. The underdosage of 7% observed with Monte Carlo simulations were found for irradiation at 0°, not used frequently in a clinical routine.     

Measurement of depth-dose of linear accelerator and simulation by use of Geant4 computer code

Radiation therapy is an established method of cancer treatment. New technologies in cancer radiotherapy need a more accurate computation of the dose delivered in the radiotherapy treatment plan. This study presents some results of a Geant4-based application for simulation of the absorbed dose distribution given by a medical linear accelerator (LINAC). The LINAC geometry is accurately described in the Monte Carlo code with use of the accelerator manufacturer''s specifications. The capability of the software for evaluating the dose distribution has been verified by comparisons with measurements in a water phantom; the comparisons were performed for percentage depth dose (PDD) and profiles for various field sizes and depths, for a 6-MV electron beam. Experimental and calculated dose values were in good agreement both in PDD and in transverse sections of the water phantom.     

Monte carlo model and output factors of elekta infinity™ 6 and 10 MV photon beam

AimThis study aimed to commission the Elekta Infinity™ working in 6 and 10 MV photon beam installed in Concord International Hospital, Singapore, and compare the OFs between MC simulation and measurement using PTW semiflex and microDiamond detector for small field sizes.Material and MethodsThere are two main steps in this study: modelling of Linac 6 and 10 MV photon beam and analysis of the output factors for field size 2 × 2–10 × 10 cm². The EGSnrc/BEAMnrc-DOSXYZnrc code was used to model and characterize the Linac and to calculate the dose distributions in a water phantom. The dose distribution and OFs were compared to the measurement data in the same condition.ResultsThe commissioning process was only conducted for a 10 × 10 cm² field size. The PDD obtained from MC simulation showed a good agreement with the measurement. The local dose difference of PDDs was less than 2% for 6 and 10 MV. The initial electron energy was 5.2 and 9.4 MeV for 6 and 10 MV photon beam, respectively. This Linac model can be used for dose calculation in other situations and different field sizes because this Linac has been commissioned and validated using Monte Carlo simulation. The 10 MV Linac produces higher electron contamination than that of 6 MV.ConclusionsThe Linac model in this study was acceptable. The most important result in this work comes from OFs resulted from MC calculation. This value was more significant than the OFs from measurement using semiflex and microDiamond for all beam energy and field sizes because of the CPE phenomenon.     

Visualisation of Respiratory Tumour Motion and Co-Moving Isodose Lines in the Context of Respiratory Gating,IMRT and Flattening-Filter-Free Beams

Respiratory motion during percutaneous radiotherapy can be considered based on respiration-correlated computed tomography (4DCT). However, most treatment planning systems perform the dose calculation based on a single primary CT data set, even though cine mode displays may allow for a visualisation of the complete breathing cycle. This might create the mistaken impression that the dose distribution were independent of tumour motion. We present a movie visualisation technique with the aim to direct attention to the fact that the dose distribution migrates to some degree with the tumour and discuss consequences for gated treatment, IMRT plans and flattening-filter-free beams. This is a feasibility test for a visualisation of tumour and isodose motion. Ten respiratory phases are distinguished on the CT, and the dose distribution from a stationary IMRT plan is calculated on each phase, to be integrated into a movie of tumour and dose motion during breathing. For one example patient out of the sample of five lesions, the plan is compared with a gated treatment plan with respect to tumour coverage and lung sparing. The interplay-effect for small segments in the IMRT plan is estimated. While the high dose rate, together with the cone-shaped beam profile, makes the use of flattening-filter-free beams more problematic for conformal and IMRT treatment, it can be the option of choice if gated treatment is preferred. The different effects of respiratory motion, dose build-up and beam properties (segments and flatness) for gated vs. un-gated treatment can best be considered if planning is performed on the full 4DCT data set, which may be an incentive for future developments of treatment planning systems.     

Variance-based sensitivity analysis of biological uncertainties in carbon ion therapy

PurposeBiological models to estimate the relative biological effectiveness (RBE) or the equivalent dose in 2 Gy fractions (EQD2) are needed for treatment planning and plan evaluation in carbon ion therapy. We present a model-independent, Monte Carlo based sensitivity analysis (SA) approach to quantify the impact of different uncertainties on the biological models.Methods and materialsThe Monte Carlo based SA is used for the evaluation of variations in biological parameters. The key property of this SA is the high number of simulation runs, each with randomized input parameters, allowing for a statistical variance-based ranking of the input variations. The potential of this SA is shown in a simplified one-dimensional treatment plan optimization. Physical properties of carbon ion beams (e.g. fragmentation) are simulated using the Monte Carlo code FLUKA. To estimate biological effects of ion beams compared to X-rays, we use the Local Effect Model (LEM) in the framework of the linear-quadratic (LQ) model. Currently, only uncertainties in the output of the biological models are taken into account.Results/conclusionsThe presented SA is suitable for evaluation of the impact of variations in biological parameters. Major advantages are the possibility to access and display the sensitivity of the evaluated quantity on several parameter variations at the same time. Main challenges for later use in three-dimensional treatment plan evaluation are computational time and memory usage. The presented SA can be performed with any analytical or numerical function and hence be applied to any biological model used in carbon ion therapy.