Full Monte Carlo and measurement-based overall performance assessment of improved clinical implementation of eMC algorithm with emphasis on lower energy range

New version 13.6.23 of the electron Monte Carlo (eMC) algorithm in Varian Eclipse™ treatment planning system has a model for 4 MeV electron beam and some general improvements for dose calculation. This study provides the first overall accuracy assessment of this algorithm against full Monte Carlo (MC) simulations for electron beams from 4 MeV to 16 MeV with most emphasis on the lower energy range. Beams in a homogeneous water phantom and clinical treatment plans were investigated including measurements in the water phantom. Two different material sets were used with full MC: (1) the one applied in the eMC algorithm and (2) the one included in the Eclipse™ for other algorithms. The results of clinical treatment plans were also compared to those of the older eMC version 11.0.31. In the water phantom the dose differences against the full MC were mostly less than 3% with distance-to-agreement (DTA) values within 2 mm. Larger discrepancies were obtained in build-up regions, at depths near the maximum electron ranges and with small apertures. For the clinical treatment plans the overall dose differences were mostly within 3% or 2 mm with the first material set. Larger differences were observed for a large 4 MeV beam entering curved patient surface with extended SSD and also in regions of large dose gradients. Still the DTA values were within 3 mm. The discrepancies between the eMC and the full MC were generally larger for the second material set. The version 11.0.31 performed always inferiorly, when compared to the 13.6.23.     

High resolution ion chamber array delivery quality assurance for robotic radiosurgery: Commissioning and validation

PurposeHigh precision radiosurgery demands comprehensive delivery-quality-assurance techniques. The use of a liquid-filled ion-chamber-array for robotic-radiosurgery delivery-quality-assurance was investigated and validated using several test scenarios and routine patient plans.Methods and materialPreliminary evaluation consisted of beam profile validation and analysis of source–detector-distance and beam-incidence-angle response dependence. The delivery-quality-assurance analysis is performed in four steps: (1) Array-to-plan registration, (2) Evaluation with standard Gamma-Index criteria (local-dose-difference ⩽ 2%, distance-to-agreement ⩽ 2 mm, pass-rate ⩾ 90%), (3) Dose profile alignment and dose distribution shift until maximum pass-rate is found, and (4) Final evaluation with 1 mm distance-to-agreement criterion. Test scenarios consisted of intended phantom misalignments, dose miscalibrations, and undelivered Monitor Units. Preliminary method validation was performed on 55 clinical plans in five institutions.ResultsThe 1000SRS profile measurements showed sufficient agreement compared with a microDiamond detector for all collimator sizes. The relative response changes can be up to 2.2% per 10 cm source–detector-distance change, but remains within 1% for the clinically relevant source–detector-distance range. Planned and measured dose under different beam-incidence-angles showed deviations below 1% for angles between 0° and 80°. Small-intended errors were detected by 1 mm distance-to-agreement criterion while 2 mm criteria failed to reveal some of these deviations. All analyzed delivery-quality-assurance clinical patient plans were within our tight tolerance criteria.ConclusionWe demonstrated that a high-resolution liquid-filled ion-chamber-array can be suitable for robotic radiosurgery delivery-quality-assurance and that small errors can be detected with tight distance-to-agreement criterion. Further improvement may come from beam specific correction for incidence angle and source–detector-distance response.     

A novel method for dose distribution registration using fiducial marks made by a megavoltage beam in film dosimetry for intensity-modulated radiation therapy quality assurance

PurposePhotographic film is widely used for the dose distribution verification of intensity-modulated radiation therapy (IMRT). However, analysis for verification of the results is subjective. We present a novel method for marking the isocenter using irradiation from a megavoltage (MV) beam transmitted through slits in a multi-leaf collimator (MLC).MethodsWe evaluated the effect of the marking irradiation at 500 monitor units (MU) on the total transmission through the MLC using an ionization chamber and Radiochromic Film. Film dosimetry was performed for quality assurance (QA) of IMRT plans. Three methods of registration were used for each film: marking by irradiating with an MV beam through slits in the MLC (MLC-IC); marking with a fabricated phantom (Phantom-IC); and a subjective method based on isodose lines (Manual). Each method was subjected to local γ-analysis.ResultsThe effect of the marking irradiation on the total transmission was 0.16%, as measured by a ionization chamber at a 10-cm depth in a solid phantom, while the inter-leaf transmission was 0.3%, determined from the film. The mean pass rates for each registration method agreed within ±1% when the criteria used were a distance-to-agreement (DTA) of 3 mm and a dose difference (DD) of 3%. For DTA/DD criteria of 2 mm/3%, the pass rates in the sagittal plane were 96.09 ± 0.631% (MLC-IC), 96.27 ± 0.399% (Phantom-IC), and 95.62 ± 0.988% (Manual).ConclusionThe present method is a versatile and useful method of improving the objectivity of film dosimetry for IMRT QA.     

Variable angle stereo imaging for rapid patient position correction in an in-house real-time position monitoring system

PurposeTo develop and validate a variable angle stereo image based position correction methodology in an X-ray based in-house online position monitoring system.Materials and methodsA stereo imaging module that enables 3D position determination and couch correction of the patient based on images acquired at any arbitrary angle and arbitrary angular separation was developed and incorporated to the in-house SeedTracker real-time position monitoring system. The accuracy of the developed system was studied by imaging an anthropomorphic phantom implanted with radiopaque markers set to known offset positions from its reference position in an Elekta linear accelerator (LA) and associated XVI imaging system. The accuracy of the system was further validated using CBCT data set from 10 prostate SBRT patients. The time gains achieved with the stereo image based position correction was compared with the manual matching of seed positions in Digitally Reconstructed Radiographs (DRRs) and kV images in the Mosaiq record and verify system.ResultsBased on phantom and patient CBCT dataset study stereo imaging module implemented in the SeedTracker shown to have an accuracy of 0.1(σ = 0.5) mm in detecting the 3D position offset. The time comparison study showed that stereo image based methodology implemented in SeedTracker was a minimum of 80(4) s faster than the manual method implemented in Mosaiq R&V system with a maximum time saving of 146(6) s.ConclusionThe variable angle stereo image based position correction method was shown to be accurate and faster than the standard manual DRR–kV image based correction approach, leading to more efficient treatment.     

Evaluation criteria for film based intensity modulated radiation therapy quality assurance

The aim of this study was to use different gamma histogram criteria for the comparison of planned dose with irradiated dose distribution and find that what percent of pixels passing a certain criteria imitate a good quality plan. The dose was calculated for 156 patients by inverse planning optimization using the Corvus treatment planning system. Gafchromic films in combination with 2571 0.6 cm³ Farmer type ionization chamber and Farmer 2570/1 electrometer from NE Technology were used to measure the delivered dose in solid water phantom. All the measurements were performed on Varian CL21EX linear accelerator (Varian Medical Systems, Palo Alto, CA) fitted with a Millennium 120 leaf collimator. In this study the mean value of the percent of passing pixels within the region of interest under the criterion of 3% DD and 3 mm DTA is 90.2 ± 7.1% for head and neck cases and 92.2 ± 5.8% for non-head and neck cases. If we choose the criteria of 3% DD and 3 mm DTA then 96.3% head and neck plans have the percent of passing pixels ≥ 75% and 95.1% non-head and neck plans have the percent of passing pixels ≥ 80%. It is evident from the results of this study that the criterion of 5% DD and 3 mm DTA with the percent of passing pixels ≥ 90 for non-head and neck cases while the percent of passing pixels ≥ 85 for head and neck cases endorse that a plan is good. The results of this study may be useful for other institutions which use verification software and EBT films for patient specific IMRT QA.     

“Edge-on” MOSkin detector for stereotactic beam measurement and verification

Dosimetry in small radiation field is challenging and complicated because of dose volume averaging and beam perturbations in a detector. We evaluated the suitability of the “Edge-on” MOSkin (MOSFET) detector in small radiation field measurement. We also tested the feasibility for dosimetric verification in stereotactic radiosurgery (SRS) and stereotactic radiotherapy (SRT). “Edge-on” MOSkin detector was calibrated and the reproducibility and linearity were determined. Lateral dose profiles and output factors were measured using the “Edge-on” MOSkin detector, ionization chamber, SRS diode and EBT2 film. Dosimetric verification was carried out on two SRS and five SRT plans. In dose profile measurements, the “Edge-on” MOSkin measurements concurred with EBT2 film measurements. It showed full width at half maximum of the dose profile with average difference of 0.11 mm and penumbral width with difference of ±0.2 mm for all SRS cones as compared to EBT2 film measurement. For output factor measurements, a 1.1% difference was observed between the “Edge-on” MOSkin detector and EBT2 film for 4 mm SRS cone. The “Edge-on” MOSkin detector provided reproducible measurements for dose verification in real-time. The measured doses concurred with the calculated dose for SRS (within 1%) and SRT (within 3%). A set of output correction factors for the “Edge-on” MOSkin detector for small radiation fields were derived from EBT2 film measurement and presented. This study showed that the “Edge-on” MOSkin detector is a suitable tool for dose verification in small radiation field.     

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.     

In phantom assessment of superficial doses under TomoTherapy irradiation

PurposeAim of this work is the assessment of build-up and superficial doses of different clinical Head&Neck plans delivered with Helical TomoTherapy (HT) (Accuray, Sunnyvale, CA). Depth dose profiles and superficial dose points were measured in order to evaluate the Treatment Planning System (TPS) capability of an accurate dose modeling in regions of disequilibrium. Geometries and scattering conditions were investigated, similar to the ones generally encountered in clinical treatments.MethodsMeasurements were performed with two dosimeters: Gafchromic® EBT3 films (Ashland Inc., Wayne, NJ) and a synthetic single crystal diamond detector (PTW-Frieburg microDiamond, MD). A modified version of the Alderson RANDO phantom was employed to house the detectors. A comparison with TPS data was carried out in terms of dose difference (DD) and distance-to-agreement (DTA).ResultsDD between calculated data and MD measurements are within 4% even in points with high spatial dose variation. For depth profiles, EBT3 data show a DDmax of 3.3% and DTAmax of 2.2 mm, in low and high gradient regions, respectively, and compare well with MD data. EBT3 superficial points always results in measured doses lower than TPS evaluated ones, with a maximum DTA value of 1.5 mm.ConclusionsDoses measured with the two devices are in good agreement and compare well with calculated data. The deviations found in the present work are within the reference tolerance level, suggesting that the HT TPS is capable of a precise dose estimation both in superficial regions and in correspondence with interfaces between air and PMMA.     

Exit fluence analysis using portal dosimetry in volumetric modulated arc therapy

AimIn measuring exit fluences, there are several sources of deviations which include the changes in the entrance fluence, changes in the detector response and patient orientation or geometry. The purpose of this work is to quantify these sources of errors.BackgroundThe use of the volumetric modulated arc therapy treatment with the help of image guidance in radiotherapy results in high accuracy of delivering complex dose distributions while sparing critical organs. The transit dosimetry has the potential of Verifying dose delivery by the linac, Multileaf collimator positional accuracy and the calculation of dose to a patient or phantom.Materials and methodsThe quantification of errors caused by a machine delivery is done by comparing static and arc picket fence test for 30 days. A RapidArc plan, created for the pelvis site was delivered without and with Rando phantom and exit portal images were acquired. The day to day dose variation were analysed by comparing the daily exit dose images during the course of treatment. The gamma criterion used for analysis is 3% dose difference and 3 mm distance to agreement with a threshold of 10% of maximum dose.ResultsThe maximum standard deviation for the static and arc picket fence test fields were 0.19 CU and 1.3 CU, respectively. The delivery of the RapidArc plans without a phantom shows the maximum standard deviation of 1.85 CU and the maximum gamma value of 0.59. The maximum gamma value for the RapidArc plan delivered with the phantom was found to be 1.2. The largest observed fluence deviation during the delivery to patient was 5.7% and the maximum standard deviation was 4.1 CU.ConclusionIt is found from this study that the variation due to patient anatomy and interfraction organ motion is significant.     

Quantification of residual dose estimation error on log file-based patient dose calculation

PurposeThe log file-based patient dose estimation includes a residual dose estimation error caused by leaf miscalibration, which cannot be reflected on the estimated dose. The purpose of this study is to determine this residual dose estimation error.Methods and materialsModified log files for seven head-and-neck and prostate volumetric modulated arc therapy (VMAT) plans simulating leaf miscalibration were generated by shifting both leaf banks (systematic leaf gap errors: ±2.0, ±1.0, and ±0.5 mm in opposite directions and systematic leaf shifts: ±1.0 mm in the same direction) using MATLAB-based (MathWorks, Natick, MA) in-house software. The generated modified and non-modified log files were imported back into the treatment planning system and recalculated. Subsequently, the generalized equivalent uniform dose (gEUD) was quantified for the definition of the planning target volume (PTV) and organs at risks.ResultsFor MLC leaves calibrated within ±0.5 mm, the quantified residual dose estimation errors that obtained from the slope of the linear regression of gEUD changes between non- and modified log file doses per leaf gap are in head-and-neck plans 1.32 ± 0.27% and 0.82 ± 0.17 Gy for PTV and spinal cord, respectively, and in prostate plans 1.22 ± 0.36%, 0.95 ± 0.14 Gy, and 0.45 ± 0.08 Gy for PTV, rectum, and bladder, respectively.ConclusionsIn this work, we determine the residual dose estimation errors for VMAT delivery using the log file-based patient dose calculation according to the MLC calibration accuracy.     

Hippocampal dose during Linac-based stereotactic radiotherapy for brain metastases: An observational study

IntroductionAim of the present study is to evaluate homolateral and contralateral hippocampus (H-H, C-H, respectively) dose during Fractionated Stereotactic Radiotherapy (FSRT) or Radiosurgery (SRS) for brain metastases (BM).Materials & methodsPatients with BM < 5, size ≤ 30 mm, KPS ≥ 80 and a life expectancy > 3 months, were considered for SRS/FSRT (total dose 15–30 Gy, 1–5 fractions). For each BM, a Flattening Filter Free (FFF) Volumetric Modulated Arc Therapy (VMAT) plan was generated with one or two arcs. Hippocampi were not considered during optimizations phase and were contoured and evaluated retrospectively in terms of dose: the Dmedian, Dmean, D0.1cc and the V1Gy, V2Gy, V5Gy and V10Gy were analyzed.ResultsFrom April 2014 to December 2015, 81 BM were treated with FFF-FSRT/SRS. For the H-H, the average values of Dmedian, Dmean and D0.1cc were 1.5Gy, 1.54Gy and 2.2Gy, respectively, while the V1Gy, V2Gy, V5Gy and V10Gy values were 25%, 8.9%, 8.9% and 2.1%, respectively. For the C–H, the average Dmedian, Dmean and D0.1 cc were 0.7Gy, 0.7Gy, 0.9Gy, respectively, while the average values of V1Gy, V2Gy, V5Gy and V10Gy were 18%, 10.2%, 2.8% and 1.4%, respectively. Tumor dimension, tumor cranial-caudal length and the distance between BM and H-H were correlated to Dmedian, Dmean and D0.1cc. For C-H, only the distance from PTV was correlated with a dose reduction.ConclusionDuring FFF-FSRT/SRS, hippocampus received a negligible dose. Despite its clinical significance is still under evaluation, in patients with a long life expectancy, H-H should be considered during Linac-based FSRT/SRS.     

Utilization of cone-beam CT for offline evaluation of target volume coverage during prostate image-guided radiotherapy based on bony anatomy alignment

AimTo assess target volume coverage during prostate image-guided radiotherapy based on bony anatomy alignment and to assess possibility of safety margin reduction.BackgroundImplementation of IGRT should influence safety margins. Utilization of cone-beam CT provides current 3D anatomic information directly in irradiation position. Such information enables reconstruction of the actual dose distribution.Materials and methodsSeventeen prostate patients were treated with daily bony anatomy image-guidance. Cone-beam CT (CBCT) scans were acquired once a week immediately after bony anatomy alignment. After the prostate, seminal vesicles, rectum and bladder were contoured, the delivered dose distribution was reconstructed. Target dose coverage was evaluated by the proportion of the CTV encompassed by the 95% isodose. Original plans employed a 1 cm safety margin. Alternative plans assuming a smaller 7 mm margin between CTV and PTV were evaluated in the same way. Rectal and bladder volumes were compared with the initial ones. Rectal and bladder volumes irradiated with doses higher than 75 Gy, 70 Gy, 60 Gy, 50 Gy and 40 Gy were analyzed.ResultsIn 12% of reconstructed plans the prostate coverage was not sufficient. The prostate underdosage was observed in 5 patients. Coverage of seminal vesicles was not satisfactory in 3% of plans. Most of the target underdosage corresponded to excessive rectal or bladder filling. Evaluation of alternative plans assuming a smaller 7 mm margin revealed 22% and 11% of plans where prostate and seminal vesicles coverage, respectively, was compromised. These were distributed over 8 and 7 patients, respectively.ConclusionSufficient dose coverage of target volumes was not achieved for all patients. Reducing of safety margin is not acceptable. Initial rectal and bladder volumes cannot be considered representative for subsequent treatment.     

Measurement of the influence of titanium hip prosthesis on therapeutic electron beam dose distributions in a novel pelvic phantom

PurposeTo investigate the degree of 18 and 22 MeV electron beam dose perturbations caused by unilateral hip titanium (Ti) prosthesis.MethodsMeasurements were acquired using Gafchromic EBT2 film in a novel pelvic phantom made out of Nylon-12 slices in which a Ti-prosthesis is embedded. Dose perturbations were measured and compared using depth doses for 8 × 8, 10 × 10 and 11 × 11 cm² applicator-defined field sizes at 95 cm source-surface-distance (SSD). Comparisons were also made between film data at 100 cm SSD for a 10 × 10 cm² field and dose calculations made on CMS XiO treatment planning system utilizing the pencil beam algorithm. The extent of dose deviations caused by the Ti prosthesis based on film data was quantified through the dose enhancement factor (DEF), defined as the ratio of the dose influenced by the prosthesis and the unchanged beam.ResultsAt the interface between Nylon-12 and the Ti implant on the prosthesis entrance side, the dose increased to values of 21 ± 1% and 23 ± 1% for 18 and 22 MeV electron beams, respectively. DEFs increased with increasing electron energy and field size, and were found to fall off quickly with distance from the nylon-prosthesis interface. A comparison of film and XiO depth dose data for 18 and 22 MeV gave relative errors of 20% and 25%, respectively.ConclusionThis study outlines the lack of accuracy of the XiO TPS for electron planning in highly heterogeneous media. So a dosimetric error of 20–25% could influence clinical outcome.     

DEMAT: A multi-institutional dosimetry audit of rotational and static intensity-modulated radiotherapy

PurposeStatic beam intensity-modulated-radiation-therapy (IMRT) and/or Volumetric-Modulated-Arc-Therapy (VMAT) are now available in many regional radiotherapy departments. The aim of this multi-institutional audit was to design a new methodology based on radiochromic films to perform an independent quality control.MethodsA set of data were sent to all participating centres for two clinical localizations: prostate and Head and Neck (H&N) cancers. The agreement between calculations and measurements was verified in the Octavius phantom (PTW) by point measurements using ionization chambers and by 2D measurements using EBT3 radiochromic films. Due to uncertainties in the whole procedure, criteria were set to 5% and 3% in local dose and 3 mm in distance excluding doses lower than 10% of the maximum doses. No normalization point or area was used for the quantitative analysis.Results13 radiotherapy centres participated in this audit involving 28 plans (12 IMRT, 16 VMAT). For point measurements, mean errors were −0.18 ± 1.54% and 0.00 ± 1.58% for prostate and H&N cases respectively. For 2D measurements with 5%/3 mm criteria, gamma map analysis showed a pixel pass rate higher than 95% for prostate and H&N. Mean gamma index was lower than 0.4 for prostate and 0.5 for H&N. Both techniques yielded similar results.ConclusionThis study showed the feasibility of an independent quality control by peers for conventional IMRT and VMAT. Results from all participating centres were found to be in good agreement. This regional study demonstrated the feasibility of our new methodology based on radiochromic films without dose normalization on a specific point.     

Fast protocol for radiochromic film dosimetry using a cloud computing web application

PurposeTo investigate the feasibility of a fast protocol for radiochromic film dosimetry to verify intensity-modulated radiotherapy (IMRT) plans.Method and materialsEBT3 film dosimetry was conducted in this study using the triple-channel method implemented in the cloud computing application (Radiochromic.com). We described a fast protocol for radiochromic film dosimetry to obtain measurement results within 1 h.Ten IMRT plans were delivered to evaluate the feasibility of the fast protocol. The dose distribution of the verification film was derived at 15, 30, 45 min using the fast protocol and also at 24 h after completing the irradiation. The four dose maps obtained per plan were compared using global and local gamma index (5%/3 mm) with the calculated one by the treatment planning system. Gamma passing rates obtained for 15, 30 and 45 min post-exposure were compared with those obtained after 24 h.ResultsSmall differences respect to the 24 h protocol were found in the gamma passing rates obtained for films digitized at 15 min (global: 99.6% ± 0.9% vs. 99.7% ± 0.5%; local: 96.3% ± 3.4% vs. 96.3% ± 3.8%), at 30 min (global: 99.5% ± 0.9% vs. 99.7% ± 0.5%; local: 96.5% ± 3.2% vs. 96.3 ± 3.8%) and at 45 min (global: 99.2% ± 1.5% vs. 99.7% ± 0.5%; local: 96.1% ± 3.8% vs. 96.3 ± 3.8%).ConclusionsThe fast protocol permits dosimetric results within 1 h when IMRT plans are verified, with similar results as those reported by the standard 24 h protocol.     

Dosimetric comparison of RapidPlan and manually optimized plans in volumetric modulated arc therapy for prostate cancer

PurposeThis study evaluated whether RapidPlan based plans (RP plans) created by a single optimization, are usable in volumetric modulated arc therapy (VMAT) for patients with prostate cancer.MethodsWe used 51 previously administered VMAT plans to train a RP model. Thirty RP plans were created by a single optimization without planner intervention during optimization. Differences between RP plans and clinical manual optimization (CMO) plans created by an experienced planner for the same patients were analyzed (Wilcoxon tests) in terms of homogeneity index (HI), conformation number (CN), D_95%, and D_2% to planning target volume (PTV), mean dose, V_50Gy, V_70Gy, V_75Gy, and V_78Gy to rectum and bladder, monitor unit (MU), and multi-leaf collimator (MLC) sequence complexity.ResultsRP and CMO values for PTV D_95%, PTV D_2%, HI, and CN were significantly similar (p < 0.05 for all). RP mean dose, V_50Gy, and V_70Gy to rectum were superior or comparable to CMO values; RP V_75Gy and V_78Gy were higher than in CMO plans (p < 0.05). RP bladder dose-volume parameter values (except V_78Gy) were lower than in CMO plans (p < 0.05). MU values were RP: 730 ± 55 MU and CMO: 580 ± 37 MU (p < 0.05); and MLC sequence complexity scores were RP: 0.25 ± 0.02 and CMO: 0.35 ± 0.03 (p < 0.05).ConclusionsRP plans created by a single optimization were clinically acceptable in VMAT for patient with prostate cancer. Our simple model could reduce optimization time, independently of planner’s skill and knowledge.