Tolerance limits and methodologies for IMRT measurement-based verification QA: Recommendations of AAPM Task Group No. 218

Evaluation of complexity and deliverability of prostate cancer treatment plans designed with a knowledge-based VMAT planning technique

PURPOSE

Knowledge-based planning (KBP) techniques have been reported to improve plan quality, efficiency, and consistency in radiation therapy. However, plan complexity and deliverability have not been addressed previously for treatment plans guided by an established in-house KBP system. The purpose of this work was to assess dosimetric, mechanical, and delivery properties of plans designed with a common KBP method for prostate cases treated via volumetric modulated arc therapy (VMAT).

METHODS

Thirty-one prostate patients previously treated with VMAT were replanned with an in-house KBP method based on the overlap volume histogram. VMAT plan complexities of the KBP plans and the reference clinical plans were quantified via monitor units, modulation complexity scores, the edge metric, and average leaf motion per degree of gantry rotation. Each set of plans was delivered to the same diode array and agreement between computed and measured dose distributions was evaluated using the gamma index. Varying percent dose-difference (1-3%) and distance-to-agreement (1 mm to 3 mm) thresholds were assessed for gamma analyses.

RESULTS

Knowledge-based planning (KBP) plans achieved average reductions of 6.4 Gy (P < 0.001) and 8.2 Gy (P < 0.001) in mean bladder and rectum dose compared to reference plans, while maintaining clinically acceptable target dose. However, KBP plans were significantly more complex than reference plans in each evaluated metric (P < 0.001). KBP plans also showed significant reductions (P < 0.05) in gamma passing rates at each evaluated criterion compared to reference plans.

CONCLUSIONS

While KBP plans had significantly reduced bladder and rectum dose, they were significantly more complex and had significantly worse quality assurance outcomes than reference plans. These results suggest caution should be taken when implementing an in-house KBP technique.

Technical Note: Dosimetric characterization of the dynamic beam flattening MLC sequence on a ring shaped, Jawless Linear Accelerator with double stacked MLC

PURPOSE

To characterize the dosimetric features and limitations of the dynamic beam flattening (DBF) on the Halcyon 2.0 linear accelerator (Varian Medical Systems).

METHODS

A predefined multi-leaf collimator (MLC) sequence was introduced and used to flatten the 6 MV flattening filter free (FFF) beam on the Halcyon 2.0. Dosimetric characterizations of the flattened beams, including beam flatness, symmetry, percent depth dose (PDD), output factor and MU linearity, were investigated. Flatness and symmetry were obtained from profile measurements with both radiographic films (EDR2) and a two dimensional ion-chamber array (IC Profiler, Sun Nuclear Corporation). MU linearity, output factors, and PDDs were measured in a water tank with a CC13 ion chamber (Scanditronix Wellhöfer, Nuremburg, Germany). In addition, the effect of the DBF sequence on 3D plan quality was evaluated by creating DBF plans for a 4-field box rectum and an AP/PA spine plan. Patient specific QA was performed on these plans.

RESULTS

At 100 cm SSD and 10 cm depth, a flatness of <3% was observed on both transversal and radial profiles for all square field sizes ≥10 cm with DBF. For both larger and smaller field sizes the flatness showed a tendency to increase as the fields got bigger or smaller, respectively. Similar trends in flatness were observed at all depths measured. All measured output factors for square field sizes ≥5 cm were within 1% of the TPS prediction. Linearity was ≤2.02% for all measurements. For both treatment sites, the MD judged the plans created for the Halcyon without the use of DBF not to be clinically acceptable, however considered both the TrueBeam plan and the Halcyon plan with the DBF sequence to be clinically acceptable.

CONCLUSIONS

The DBF sequence on the Halcyon and its characteristics were investigated. The analysis indicates that the DBF sequence can be used on the Halcyon to generate clinically acceptable 3D treatment plans.

Technical note: A modified gamma evaluation method for dose distribution comparisons

Purpose

In this work we have developed a novel method of dose distribution comparison, the inverse gamma (IG) evaluation, by modifying the commonly used gamma evaluation method.

Methods

The IG evaluation calculates the gamma criteria (dose difference criterion, ΔD, or distance‐to‐agreement criterion, Δd) that are needed to achieve a predefined pass rate or gamma agreement index (GAI). In‐house code for evaluating IG with a fixed ΔD of 3% was developed using Python (v3.5.2) and investigated using treatment plans and measurement data from 25 retrospective patient specific quality assurance tests (53 individual arcs).

Results

It was found that when the desired GAI was set to 95%, approximately three quarters of the arcs tested were able to achieve Δd within 1 mm (mean Δd: 0.7 ± 0.5 mm). The mean Δd required in order for all points to pass the gamma evaluation (i.e., GAI = 100%) was 4.5 ± 3.1 mm. The possibility of evaluating IG by fixing the Δd or ΔD/Δd, instead of fixing the ΔD at 3%, was also investigated.

Conclusion

The IG method and its indices have the potential to be implemented clinically to quantify the minimum dose and distance criteria based on a specified GAI. This method provides additional information to augment standard gamma evaluation results during patient specific quality assurance testing of individual treatment plans. The IG method also has the potential to be used in retrospective audits to determine an appropriate set of local gamma criteria and action levels based on a cohort of patient specific quality assurance plans.

Hybrid 3D analytical linear energy transfer calculation algorithm based on precalculated data from Monte Carlo simulations

PURPOSE

The dose-averaged linear energy transfer (LET_d ) for intensity-modulated proton therapy (IMPT) calculated by one-dimensional (1D) analytical models deviates from more accurate but time-consuming Monte Carlo (MC) simulations. We developed a fast hybrid three-dimensional (3D) analytical LET_d calculation that is more accurate than 1D analytical model.

METHODS

We used the Geant4 MC code to generate 3D LET_d distributions of monoenergetic proton beams in water for all energies and used a customized error function to fit the LET_d lateral profiles at various depths to the MC simulation. The 3D LET_d calculation kernel was a lookup table of these fitted coefficients, and LET_d was determined directly from spot energies and voxel coordinates during analytical dose calculations. We validated our new method by comparing the calculated LET_d distributions to MC results using 3D Gamma index analysis with 3%/2 mm criteria in 12 patient geometries. The significance of the improvement in Gamma index analysis passing rates over the 1D analytical model was determined using the Wilcoxon rank-sum test.

RESULTS

The passing rate of 3D Gamma analysis comparing LET_d distributions from the hybrid 3D method and the 1D method to MC simulations was significantly improved from 94.0% ± 2.5% to 98.0% ± 1.0% (P = 0.0003). The typical time to calculate dose and LET_d simultaneously using an Intel Xeon E5-2680 2.50 GHz workstation was approximately 2.5 min.

CONCLUSIONS

Our new method significantly improved the LET_d calculation accuracy compared to the 1D method while maintaining significantly shorter calculation time even comparing with the GPU-based fast MC code.

Transitioning from measurement-based to combined patient-specific quality assurance for intensity-modulated proton therapy

OBJECTIVE

This study is part of ongoing efforts aiming to transit from measurement-based to combined patient-specific quality assurance (PSQA) in intensity-modulated proton therapy (IMPT). A Monte Carlo (MC) dose-calculation algorithm is used to improve the independent dose calculation and to reveal the beam modeling deficiency of the analytical pencil beam (PB) algorithm.

METHODS

A set of representative clinical IMPT plans with suboptimal PSQA results were reviewed. Verification plans were recalculated using an MC algorithm developed in-house. Agreements of PB and MC calculations with measurements that quantified by the γ passing rate were compared.

RESULTS

The percentage of dose planes that met the clinical criteria for PSQA (>90% γ passing rate using 3%/3 mm criteria) increased from 71.40% in the original PB calculation to 95.14% in the MC recalculation. For fields without beam modifiers, nearly 100% of the dose planes exceeded the 95% γ passing rate threshold using the MC algorithm. The model deficiencies of the PB algorithm were found in the proximal and distal regions of the SOBP, where MC recalculation improved the γ passing rate by 11.27% (p < 0.001) and 16.80% (p < 0.001), respectively.

CONCLUSIONS

The MC algorithm substantially improved the γ passing rate for IMPT PSQA. Improved modeling of beam modifiers would enable the use of the MC algorithm for independent dose calculation, completely replacing additional depth measurements in IMPT PSQA program. For current users of the PB algorithm, further improving the long-tail modeling or using MC simulation to generate the dose correction factor is necessary.

ADVANCES IN KNOWLEDGE

We justified a change in clinical practice to achieve efficient combined PSQA in IMPT by using the MC algorithm that was experimentally validated in almost all the clinical scenarios in our center. Deficiencies in beam modeling of the current PB algorithm were identified and solutions to improve its dose-calculation accuracy were provided.

MLC parameters from static fields to VMAT plans: an evaluation in a RT-dedicated MC environment (PRIMO)

Background

PRIMO is a graphical environment based on PENELOPE Monte Carlo (MC) simulation of radiotherapy beams able to compute dose distribution in patients, from plans with different techniques. The dosimetric characteristics of an HD-120 MLC (Varian), simulated using PRIMO, were here compared with measurements, and also with Acuros calculations (in the Eclipse treatment planning system, Varian).

Materials and methods

A 10 MV FFF beam from a Varian EDGE linac equipped with the HD-120 MLC was used for this work. Initially, the linac head was simulated inside PRIMO, and validated against measurements in a water phantom. Then, a series of different MLC patterns were established to assess the MLC dosimetric characteristics. Those tests included: i) static fields: output factors from MLC shaped fields (2 × 2 to 10 × 10 cm²), alternate open and closed leaf pattern, MLC transmitted dose; ii) dynamic fields: dosimetric leaf gap (DLG) evaluated with sweeping gaps, tongue and groove (TG) effect assessed with profiles across alternate open and closed leaves moving across the field. The doses in the different tests were simulated in PRIMO and then compared with EBT3 film measurements in solid water phantom, as well as with Acuros calculations. Finally, MC in PRIMO and Acuros were compared in some clinical cases, summarizing the clinical complexity in view of a possible use of PRIMO as an independent dose calculation check.

Results

Static output factor MLC tests showed an agreement between MC calculated and measured OF of 0.5%. The dynamic tests presented DLG values of 0.033 ± 0.003 cm and 0.032 ± 0.006 cm for MC and measurements, respectively. Regarding the TG tests, a general agreement between the dose distributions of 1–2% was achieved, except for the extreme patterns (very small gaps/field sizes and high TG effect) were the agreement was about 4–5%. The analysis of the clinical cases, the Gamma agreement between MC in PRIMO and Acuros dose calculation in Eclipse was of 99.5 ± 0.2% for 3%/2 mm criteria of dose difference/distance to agreement.

Conclusions

MC simulations in the PRIMO environment were in agreement with measurements for the HD-120 MLC in a 10 MV FFF beam from a Varian EDGE linac. This result allowed to consistently compare clinical cases, showing the possible use of PRIMO as an independent dose calculation check tool.

Testing the methodology for a dosimetric end-to-end audit of IMRT/VMAT: results of IAEA multicentre and national studies

Introduction: Within an International Atomic Energy Agency (IAEA) co-ordinated research project (CRP), a remote end-to-end dosimetric quality audit for intensity modulated radiation therapy (IMRT)/ volumetric arc therapy (VMAT) was developed to verify the radiotherapy chain including imaging, treatment planning and dose delivery. The methodology as well as the results obtained in a multicentre pilot study and national trial runs conducted in close cooperation with dosimetry audit networks (DANs) of IAEA Member States are presented.Material and methods: A solid polystyrene phantom containing a dosimetry insert with an irregular solid water planning target volume (PTV) and organ at risk (OAR) was designed for this audit. The insert can be preloaded with radiochromic film and four thermoluminescent dosimeters (TLDs). For the audit, radiotherapy centres were asked to scan the phantom, contour the structures, create an IMRT/VMAT treatment plan and irradiate the phantom. The dose prescription was to deliver 4 Gy to the PTV in two fractions and to limit the OAR dose to a maximum of 2.8 Gy. The TLD measured doses and film measured dose distributions were compared with the TPS calculations.Results: Sixteen hospitals from 13 countries and 64 hospitals from 6 countries participated in the multicenter pilot study and in the national runs, respectively. The TLD results for the PTV were all within ±5% acceptance limit for the multicentre pilot study, whereas for national runs, 17 participants failed to meet this criterion. All measured doses in the OAR were below the treatment planning constraint. The film analysis identified seven plans in national runs below the 90% passing rate gamma criteria.Conclusion: The results proved that the methodology of the IMRT/VMAT dosimetric end-to-end audit was feasible for its intended purpose, i.e., the phantom design and materials were suitable; the phantom was easy to use and it was robust enough for shipment. Most importantly the audit methodology was capable of identifying suboptimal IMRT/VMAT delivery.

Commissioning and performance evaluation of RadCalc for the Elekta unity MRI-linac

Recent availability of MRI‐guided linear accelerators has introduced a number of clinical challenges, particularly in the context of online plan adaptation. Paramount among these is verification of plan quality prior to patient treatment. Currently, there are no commercial products available for monitor unit verification that fully support the newly FDA cleared Elekta Unity 1.5 T MRI‐linac. In this work, we investigate the accuracy and precision of RadCalc for this purpose, which is a software package that uses a Clarkson integration algorithm for point dose calculation. To this end, 18 IMRT patient plans (186 individual beams) were created and used for RadCalc point dose calculations. In comparison with the primary treatment planning system (Monaco), mean point dose deviations of 0.0 ± 1.0% (n = 18) and 1.7 ± 12.4% (n = 186) were obtained on a per‐plan and per‐beam basis, respectively. The dose plane comparison functionality within RadCalc was found to be highly inaccurate, however, modest improvements could be made by artificially shifting jaws and multi leaf collimator positions to account for the dosimetric shift due to the magnetic field (67.3% vs 96.5% mean 5%/5 mm gamma pass rate).

Treatment plan verification: A review on the comparison of dose distributions

PURPOSE

The aim of this review article is to provide a useful reference for dose comparison techniques within the frame of treatment plan verification. Each technique is presented with a general description given along with advantages and disadvantage and the rationale for its development.

METHODS

The review was conducted in PubMed from 1993 to 2019 including articles referring to the methodology of dose comparison for treatment plan verification.

RESULTS

The search identified thirty-one dose comparison methods that were categorized according to the number of physical parameters that take into account for dose comparison.

CONCLUSIONS

Among the available methods for the comparison of two dose distributions, the γ-analysis (gamma analysis) has been widely adopted as the gold standard in verification procedures. However, due to various intrinsic limitations of gamma index, the development of a better metric taking into account both statistical and in clinical parameters is required.

Two-dimensional solid-state array detectors: A technique for in vivo dose verification in a variable effective area

PURPOSE

We introduce a technique that employs a 2D detector in transmission mode (TM) to verify dose maps at a depth of dmax in Solid Water. TM measurements, when taken at a different surface-to-detector distance (SDD), allow for the area at dmax (in which the dose map is calculated) to be adjusted.

METHODS

We considered the detector prototype "MP512" (an array of 512 diode-sensitive volumes, 2 mm spatial resolution). Measurements in transmission mode were taken at SDDs in the range from 0.3 to 24 cm. Dose mode (DM) measurements were made at dmax in Solid Water. We considered radiation fields in the range from 2 × 2 cm2 to 10 × 10 cm2 , produced by 6 MV flattened photon beams; we derived a relationship between DM and TM measurements as a function of SDD and field size. The relationship was used to calculate, from TM measurements at 4 and 24 cm SDD, dose maps at dmax in fields of 1 × 1 cm2 and 4 × 4 cm2 , and in IMRT fields. Calculations were cross-checked (gamma analysis) with the treatment planning system and with measurements (MP512, films, ionization chamber).

RESULTS

In the square fields, calculations agreed with measurements to within ±2.36%. In the IMRT fields, using acceptance criteria of 3%/3 mm, 2%/2 mm, 1%/1 mm, calculations had respective gamma passing rates greater than 96.89%, 90.50%, 62.20% (for a 4 cm SSD); and greater than 97.22%, 93.80%, 59.00% (for a 24 cm SSD). Lower rates (1%/1 mm criterion) can be explained by submillimeter misalignments, dose averaging in calculations, noise artifacts in film dosimetry.

CONCLUSIONS

It is possible to perform TM measurements at the SSD which produces the best fit between the area at dmax in which the dose map is calculated and the size of the monitored target.

Early clinical experience with varian halcyon V2 linear accelerator: Dual-isocenter IMRT planning and delivery with portal dosimetry for gynecological cancer treatments

PURPOSE

Varian Halcyon linear accelerator version 2 (The Halcyon 2.0) was recently released with new upgraded features. The aim of this study was to report our clinical experience with Halcyon 2.0 for a dual-isocenter intensity-modulated radiation therapy (IMRT) planning and delivery for gynecological cancer patients and examine the feasibility of in vivo portal dosimetry.

METHODS

Twelve gynecological cancer patients were treated with extended-field IMRT technique using two isocenters on Halcyon 2.0 to treat pelvis and pelvic/or para-aortic nodes region. The prescription dose was 45 Gy in 25 fractions (fxs) with simultaneous integrated boost (SIB) dose of 55 or 57.5 Gy in 25 fxs to involved nodes. All treatment plans, pretreatment patient-specific QA and treatment delivery records including daily in vivo portal dosimetry were retrospectively reviewed. For in vivo daily portal dosimetry analysis, each fraction was compared to the reference baseline (1st fraction) using gamma analysis criteria of 4 %/4 mm with 90% of total pixels in the portal image planar dose.

RESULTS

All 12 extended-field IMRT plans met the planning criteria and delivered as planned (a total of 300 fractions). Conformity Index (CI) for the primary target was achieved with the range of 0.99-1.14. For organs at risks, most were well within the dose volume criteria. Treatment delivery time was from 5.0 to 6.5 min. Interfractional in vivo dose variation exceeded gamma analysis threshold for 8 fractions out of total 300 (2.7%). These eight fractions were found to have a relatively large difference in small bowel filling and SSD change at the isocenter compared to the baseline.

CONCLUSION

Halcyon 2.0 is effective to create complex extended-field IMRT plans using two isocenters with efficient delivery. Also Halcyon in vivo dosimetry is feasible for daily treatment monitoring for organ motion, internal or external anatomy, and body weight which could further lead to adaptive radiation therapy.

Field-in-field breast planning for a jawless, double-stack MLC LINAC using flattening-filter-free beams

BACKGROUND

This study intends to develop an efficient field-in-field (FiF) planning technique with the Eclipse treatment planning system (TPS) to determine the feasibility of using the Halcyon treatment delivery system for 3D treatment of breast cancer.

METHODS

Ten treatment plans were prepared on the Halcyon treatment planning system and compared to the same patients' clinically delivered TrueBeam plans which used flattened 6 MV and 10 MV beams. Patients selected for this study were treated via simple, tangential breast irradiation and did not receive radiotherapy of the supraclavicular or internal mammary lymph nodes. Planning target volumes (PTV) volumes ranged from 519 cc to 1211 cc with a mean target volume of 877 cc. Several planning techniques involving collimator, gantry rotation, and number of FiF segments were investigated as well as the use of the dynamically flattened beam (DFB) - a predefined MLC pattern that is designed to provide a flattened beam profile at 10 cm depth on a standard water phantom. For comparison, the clinically delivered TrueBeam plans remained unaltered except for normalization of the target coverage to more readily compare the two treatment delivery techniques.

RESULTS

Using the physician defined PTV, normalized such that 98% of the volume was covered by 95% of the prescribed dose, the Halcyon plans were deemed clinically acceptable and comparable to the TrueBeam plans by the radiation oncologist. Resulting average global maximum doses in the test patients were identical between the TrueBeam and Halcyon plans (108% of Rx) and a mean PTV dose of 102.5% vs 101.6%, respectively.

CONCLUSIONS

From this study a practical and efficient planning method for delivering 3D conformal breast radiotherapy using the Halcyon linear accelerator has been developed. When normalized to the clinically desired coverage, hot spots were maintained to acceptable levels and overall plan quality was comparable to plans delivered on conventional C-arm LINACs.

Comparison of patient-specific intensity modulated radiation therapy quality assurance for the prostate across multiple institutions

Aim

To evaluate the success of a patient-specific intensity modulated radiation therapy (IMRT) quality assurance (QA) practice for prostate cancer patients across multiple institutions using a questionnaire survey.

Background

The IMRT QA practice involves different methods of dose distribution verification and analysis at different institutions.

Materials and Methods

Two full-arc volumetric modulated arc therapy (VMAT) plan and 7 fixed-gantry IMRT plan with DMLC were used for patient specific QA across 22 institutions. The same computed tomography image and structure set were used for all plans. Each institution recalculated the dose distribution with fixed monitor units and without any modification. Single-point dose measurement with a cylindrical ionization chamber and dose distribution verification with a multi-detector or radiochromic film were performed, according to the QA process at each institution.

Results

Twenty-two institutions performed the patient-specific IMRT QA verifications. With a single-point dose measurement at the isocenter, the average difference between the calculated and measured doses was 0.5 ± 1.9%. For the comparison of dose distributions, 18 institutions used a two or three-dimensional array detector, while the others used Gafchromic film. In the γ test with dose difference/distance-to-agreement criteria of 3%-3 mm and 2%-2 mm with a 30% dose threshold, the median gamma pass rates were 99.3% (range: 41.7%-100.0%) and 96.4% (range: 29.4%-100.0%), respectively.

Conclusion

This survey was an informative trial to understand the verification status of patient-specific IMRT QA measurements for prostate cancer. In most institutions, the point dose measurement and dose distribution differences met the desired criteria.

Impact of delivery characteristics on dose delivery accuracy of volumetric modulated arc therapy for different treatment sites

This study aimed to investigate the impact of delivery characteristics on the dose delivery accuracy of volumetric modulated arc therapy (VMAT) for different treatment sites. The pretreatment quality assurance (QA) results of 344 VMAT patients diagnosed with gynecological (GYN), head and neck (H&N), rectal or prostate cancer were randomly chosen in this study. Ten metrics reflecting VMAT delivery characteristics were extracted from the QA plans. Compared with GYN and rectal plans, H&N and prostate plans had higher aperture complexity and monitor units (MU), and smaller aperture area. Prostate plans had the smallest aperture area and lowest leaf speed compared with other plans (P < 0.001). No differences in gantry speed were found among the four sites. The gamma passing rates (GPRs) of GYN, rectal and H&N plans were inversely associated with union aperture area (UAA) and leaf speed (Pearson's r: -0.39 to -0.68). GPRs of prostate plans were inversely correlated with aperture complexity, MU and small aperture score (SAS) (absolute Pearson's r: 0.34 to 0.49). Significant differences in GPR between high SAS and low SAS subgroups were found only when leaf speed was <0.42 cm s-1 (P < 0.001). No association of GPR with gantry speed was found in four sites. Leaf speed was more strongly associated with UAA. Aperture complexity and MU were more strongly associated with SAS. VMAT plans from different sites have distinct delivery characteristics. Affecting dose delivery accuracy, leaf speed is the key factor for GYN, rectal and H&N plans, while aperture complexity, MU and small apertures have a higher influence on prostate plans.

Analysis of dose comparison techniques for patient-specific quality assurance in radiation therapy

Purpose

Gamma evaluation is the most commonly used technique for comparison of dose distributions for patient‐specific pretreatment quality assurance in radiation therapy. Alternative dose comparison techniques have been developed but not widely implemented. This study aimed to compare and evaluate the performance of several previously published alternatives to the gamma evaluation technique, by systematically evaluating a large number of patient‐specific quality assurance results.

Methods

The agreement indices (or pass rates) for global and local gamma evaluation, maximum allowed dose difference (MADD) and divide and conquer (D&C) techniques were calculated using a selection of acceptance criteria for 429 patient‐specific pretreatment quality assurance measurements. Regression analysis was used to quantify the similarity of behavior of each technique, to determine whether possible variations in sensitivity might be present.

Results

The results demonstrated that the behavior of D&C gamma analysis and MADD box analysis differs from any other dose comparison techniques, whereas MADD gamma analysis exhibits similar performance to the standard global gamma analysis. Local gamma analysis had the least variation in behavior with criteria selection. Agreement indices calculated for 2%/2 mm and 2%/3 mm, and 3%/2 mm and 3%/3 mm were correlated for most comparison techniques.

Conclusion

Radiation oncology treatment centers looking to compare between different dose comparison techniques, criteria or lower dose thresholds may apply the results of this study to estimate the expected change in calculated agreement indices and possible variation in sensitivity to delivery dose errors.

DVH analysis using a transmission detector and model-based dose verification system as a comprehensive pretreatment QA tool for VMAT plans: Clinical experience and results

Purpose

Dose volume histogram (DVH)‐based analysis is utilized as a pretreatment quality assurance tool to determine clinical relevance from measured dose which is difficult in conventional gamma‐based analysis. In this study, we report our clinical experience with an ionization‐based transmission detector and model‐based verification system, using DVH analysis, as a comprehensive pretreatment QA tool for complex volumetric modulated arc therapy plans.

Methods and Materials

Seventy‐three subsequent treatment plans categorized into four clinical sites (Head and Neck, Thorax, Abdomen, and Pelvis) were evaluated. The average dose (D_mean) and dose received by 1% (D₁) of the planning target volumes (PTVs) and organs at risks (OARs) calculated using the treatment planning system (TPS) were compared to a computed (model‐based) and reconstructed dose, from the measured fluence, using DVH analysis. The correlation between gamma (3% 3 mm) and DVH‐based analysis for targets was evaluated. Furthermore, confidence and action limits for detector and verification systems were established.

Results

Linear regression confirmed an excellent correlation between TPS planned and computed dose using a model‐based verification system (r² = 1). The average percentage difference between TPS calculated and reconstructed dose for PTVs achieved using DVH analysis for each site is as follows: Head and Neck — 0.57 ± 2.8% (D_mean) and 2.6 ± 2.7% (D₁), Abdomen — 0.19 ± 2.8% and 1.64 ± 2.2%, Thorax — 0.24 ± 2.1% and 3.12 ± 2.8%, Pelvis 0.37 ± 2.4% and 1.16 ± 2.3%, respectively. The average percentage of passed gamma values achieved was above 95% for all cases. However, no correlation was observed between gamma passing rates and DVH difference (%) for PTVs (r² = 0.11). The results demonstrate a confidence limit of 5% (D_mean and D₁) for PTVs using DVH analysis for both computed and reconstructed dose distribution.

Conclusion

DVH analysis of treatment plan using a model‐based verification system and transmission detector provided useful information on clinical relevance for all cases and could be used as a comprehensive pretreatment patient‐specific QA tool.

Evaluating the sensitivity of Halcyon's automatic transit image acquisition for treatment error detection: A phantom study using static IMRT

PURPOSE

The Varian Halcyon™ electronic portal imaging detector is always in-line with the beam and automatically acquires transit images for every patient with full-field coverage. These images could be used for "every patient, every monitor unit" quality assurance (QA) and eventually adaptive radiotherapy. This study evaluated the imager's sensitivity to potential clinical errors and day-to-day variations from clinical exit images.

METHODS

Open and modulated fields were delivered for each potential error. To evaluate output changes, monitor units were scaled by 2%-10% and delivered to solid water slabs and a homogeneous CIRS phantom. To mimic weight changes, 0.5-5.0 cm of buildup was added to the solid water. To evaluate positioning changes, a homogeneous and heterogeneous CIRS phantom were shifted 2-10 cm and 0.2-1.5 cm, respectively. For each test, mean relative differences (MRDs) and standard deviations in the pixel-difference histograms (σ_RD ) between test and baseline images were calculated. Lateral shift magnitudes were calculated using cross-correlation and edge-detection filtration. To assess patient variations, MRD and σ_RD were calculated from six prostate patients' daily exit images and compared between fractions with and without gas present.

RESULTS

MRDs responded linearly to output and buildup changes with a standard deviation of 0.3%, implying a 1% output change and 0.2 cm changes in buildup could be detected with 2.5σ confidence. Shifting the homogenous phantom laterally resulted in detectable MRD and σ_RD changes, and the cross-correlation function calculated the shift to within 0.5 mm for the heterogeneous phantom. MRD and σ_RD values were significantly associated with the presence of gas for five of the six patients.

CONCLUSIONS

Rapid analyses of automatically acquired Halcyon™ exit images could detect mid-treatment changes with high sensitivity, though appropriate thresholds will need to be set. This study presents the first steps toward developing effortless image evaluation for all aspects of every patient's treatment.

Complexity metrics for IMRT and VMAT plans: a review of current literature and applications

Modulated radiotherapy with multileaf collimators is widely used to improve target conformity and normal tissue sparing. This introduced an additional degree of complexity, studied by multiple teams through different properties. Three categories of complexity metrics were considered in this review: fluence, deliverability and accuracy metrics. The first part of this review is dedicated to the inventory of these complexity metrics. Different applications of these metrics emerged. Influencing the optimizer by integrating complexity metrics into the cost function has been little explored and requires more investigations. In modern treatment planning system, it remains confined to MUs or treatment time limitation. A large majority of studies calculated metrics only for analysis, without plan modification. The main application was to streamline the patient specific quality assurance workload, investigating the capability of complexity metrics to predict patient specific quality assurance results. Additionally complexity metrics were used to analyze behaviour of TPS optimizer, compare TPS, operators and plan properties, and perform multicentre audit. Their potential was also explored in the context of adaptive radiotherapy and automation planning. The second part of the review gives an overview of these studies based on the complexity metrics.

Feasibility of a GATE Monte Carlo platform in a clinical pretreatment QA system for VMAT treatment plans using TrueBeam with an HD120 multileaf collimator

Purpose

To evaluate the quality of patient‐specific complicated treatment plans, including commercialized treatment planning systems (TPS) and commissioned beam data, we developed a process of quality assurance (QA) using a Monte Carlo (MC) platform. Specifically, we constructed an interface system that automatically converts treatment plan and dose matrix data in digital imaging and communications in medicine to an MC dose‐calculation engine. The clinical feasibility of the system was evaluated.

Materials and Methods

A dose‐calculation engine based on GATE v8.1 was embedded in our QA system and in a parallel computing system to significantly reduce the computation time. The QA system automatically converts parameters in volumetric‐modulated arc therapy (VMAT) plans to files for dose calculation using GATE. The system then calculates dose maps. Energies of 6 MV, 10 MV, 6 MV flattening filter free (FFF), and 10 MV FFF from a TrueBeam with HD120 were modeled and commissioned. To evaluate the beam models, percentage depth dose (PDD) values, MC calculation profiles, and measured beam data were compared at various depths (D_max, 5 cm, 10 cm, and 20 cm), field sizes, and energies. To evaluate the feasibility of the QA system for clinical use, doses measured for clinical VMAT plans using films were compared to dose maps calculated using our MC‐based QA system.

Results

A LINAC QA system was analyzed by PDD and profile according to the secondary collimator and multileaf collimator (MLC). Values for MC calculations and TPS beam data obtained using CC13 ion chamber (IBA Dosimetry, Germany) were consistent within 1.0%. Clinical validation using a gamma index was performed for VMAT treatment plans using a solid water phantom and arbitrary patient data. The gamma evaluation results (with criteria of 3%/3 mm) were 98.1%, 99.1%, 99.2%, and 97.1% for energies of 6 MV, 10 MV, 6 MV FFF, and 10 MV FFF, respectively.

Conclusions

We constructed an MC‐based QA system for evaluating patient treatment plans and evaluated its feasibility in clinical practice. We observed robust agreement between dose calculations from our QA system and measurements for VMAT plans. Our QA system could be useful in other clinical settings, such as small‐field SRS procedures or analyses of secondary cancer risk, for which dose calculations using TPS are difficult to verify.

Evaluation of interfraction setup variations for postmastectomy radiation therapy using EPID-based in vivo dosimetry

Postmastectomy radiation therapy is technically difficult and can be considered one of the most complex techniques concerning patient setup reproducibility. Slight patient setup variations — particularly when high‐conformal treatment techniques are used — can adversely affect the accuracy of the delivered dose and the patient outcome. This research aims to investigate the inter‐fraction setup variations occurring in two different scenarios of clinical practice: at the reference and at the current patient setups, when an image‐guided system is used or not used, respectively. The results were used with the secondary aim of assessing the robustness of the patient setup procedure in use. Forty eight patients treated with volumetric modulated arc and intensity modulated therapies were included in this study. EPID‐based in vivo dosimetry (IVD) was performed at the reference setup concomitantly with the weekly cone beam computed tomography acquisition and during the daily current setup. Three indices were analyzed: the ratio R between the reconstructed and planned isocenter doses, γ% and the mean value of γ from a transit dosimetry based on a two‐dimensional γ‐analysis of the electronic portal images using 5% and 5 mm as dose difference and distance to agreement gamma criteria; they were considered in tolerance if R was within 5%, γ% > 90% and γmean < 0.4. One thousand and sixteen EPID‐based IVD were analyzed and 6.3% resulted out of the tolerance level. Setup errors represented the main cause of this off tolerance with an occurrence rate of 72.2%. The percentage of results out of tolerance obtained at the current setup was three times greater (9.5% vs 3.1%) than the one obtained at the reference setup, indicating weaknesses in the setup procedure. This study highlights an EPID‐based IVD system's utility in the radiotherapy routine as part of the patient’s treatment quality controls and to optimize (or confirm) the performed setup procedures’ accuracy.

Comprehensive evaluation of the high-resolution diode array for SRS dosimetry

A high-resolution diode array has been comprehensively evaluated. It consists of 1013 point diode detectors arranged on the two 7.7 × 7.7 cm² printed circuit boards (PCBs). The PCBs are aligned face to face in such a way that the active volumes of all diodes are in the same plane. All individual correction factors required for accurate dosimetry have been validated for conventional and flattening filter free (FFF) 6MV beams. That included diode response equalization, linearity, repetition rate dependence, field size dependence, angular dependence at the central axis and off-axis in the transverse, sagittal, and multiple arbitrary planes. In the end-to-end tests the array and radiochromic film dose distributions for SRS-type multiple-target plans were compared. In the equalization test (180° rotation), the average percent dose error between the normal and rotated positions for all diodes was 0.01% ± 0.1% (range -0.3 to 0.4%) and -0.01% ± 0.2% (range -0.9 to 0.9%) for 6 MV and 6MV FFF beams, respectively. For the axial angular response, corrected dose stayed within 2% from the ion chamber for all gantry angles, until the beam direction approached the detector plane. In azimuthal direction, the device agreed with the scintillator within 1% for both energies. For multiple combinations of couch and gantry angles, the average percent errors were -0.00% ± 0.6% (range: -2.1% to 1.6%) and -0.1% ± 0.5% (range -1.6% to 2.1%) for the 6MV and 6MV FFF beams, respectively. The measured output factors were largely within 2% of the scintillator, except for the 5 mm 6MV beam showing a 3.2% deviation. The 2%/1 mm gamma analysis of composite SRS measurements produced the 97.2 ± 1.3% (range 95.8-98.5%) average passing rate against film. Submillimeter (≤0.5 mm) dose profile alignment with film was demonstrated in all cases.

Modeling Elekta VersaHD using the Varian Eclipse treatment planning system for photon beams: A single-institution experience

The aim of this study was to report a single‐institution experience and commissioning data for Elekta VersaHD linear accelerators (LINACs) for photon beams in the Eclipse treatment planning system (TPS). Two VersaHD LINACs equipped with 160‐leaf collimators were commissioned. For each energy, the percent‐depth‐dose (PDD) curves, beam profiles, output factors, leaf transmission factors and dosimetric leaf gaps (DLGs) were acquired in accordance with the AAPM task group reports No. 45 and No. 106 and the vendor‐supplied documents. The measured data were imported into Eclipse TPS to build a VersaHD beam model. The model was validated by creating treatment plans spanning over the full‐spectrum of treatment sites and techniques used in our clinic. The quality assurance measurements were performed using MatriXX, ionization chamber, and radiochromic film. The DLG values were iteratively adjusted to optimize the agreement between planned and measured doses. Mobius, an independent LINAC logfile‐based quality assurance tool, was also commissioned both for routine intensity‐modulated radiation therapy (IMRT) QA and as a secondary check for the Eclipse VersaHD model. The Eclipse‐generated VersaHD model was in excellent agreement with the measured PDD curves and beam profiles. The measured leaf transmission factors were less than 0.5% for all energies. The model validation study yielded absolute point dose agreement between ionization chamber measurements and Eclipse within ±4% for all cases. The comparison between Mobius and Eclipse, and between Mobius and ionization chamber measurements lead to absolute point dose agreement within ±5%. The corresponding 3D dose distributions evaluated with 3%global/2mm gamma criteria resulted in larger than 90% passing rates for all plans. The Eclipse TPS can model VersaHD LINACs with clinically acceptable accuracy. The model validation study and comparisons with Mobius demonstrated that the modeling of VersaHD in Eclipse necessitates further improvement to provide dosimetric accuracy on par with Varian LINACs.

Departmental action limits for TQA energy variations defined by means of statistical process control methods

The purpose of this study is to define departmental action limits for energy percentage variation measured by means of step-wedge helical Tomotherapy quality assurance module. Individual charts using the Statistical Process Control techniques have been used to identify retrospectively out-of-control situations ascribable to documented actions performed on the Tomotherapy system. Using the in-control data of our analysis process capability indices (cp, cpk, cpm and cpmk) are calculated in order to document the real working condition of the Tomotherapy system. Our findings indicate use of an action limit of 1.0% for energy percentage variation difference between the measured and reference output is a good working condition of a Tomotherapy system. cp and cpk indices are suggested as good indices that correctly report the system capability. A method for calculating and reporting Tomotherapy action limits for the integrated self-checking TQA energy check was shown in this study. SPC technique has proven to be efficient in defining departmental action limits from retrospective data for TQA energy measurements, hence optimally enabling corrective improvements in the process of quality assurance.

Machine Learning for Patient-Specific Quality Assurance of VMAT: Prediction and Classification Accuracy

PURPOSE

To assess the accuracy of machine learning to predict and classify quality assurance (QA) results for volumetric modulated arc therapy (VMAT) plans.

METHODS AND MATERIALS

Three hundred three VMAT plans, including 176 gynecologic cancer and 127 head and neck cancer plans, were chosen in this study. Fifty-four complexity metrics were extracted from the QA plans and considered as inputs. Patient-specific QA was performed, and gamma passing rates (GPRs) were used as outputs. One Poisson lasso (PL) regression model was developed, aiming to predict individual GPR, and 1 random forest (RF) classification model was developed to classify QA results as "pass" or "fail." Both technical validation (TV) and clinical validation (CV) were used to evaluate the model reliability. GPR prediction accuracy of PL and classification performance of PL and RF were evaluated.

RESULTS

In TV, the mean prediction error of PL was 1.81%, 2.39%, and 4.18% at 3%/3 mm, 3%/2 mm, and 2%/2 mm, respectively. No significant differences in prediction errors between TV and CV were observed. In QA results classification, PL had a higher specificity (accurately identifying plans that can pass QA), whereas RF had a higher sensitivity (accurately identifying plans that may fail QA). By using 90% as the action limit at a 3%/2 mm criterion, the specificity of PL and RF was 97.5% and 87.7% in TV and 100% and 71.4% in CV, respectively. The sensitivity of PL and RF was 31.6% and 100% in TV and 33.3% and 100% in CV, respectively. With 100% sensitivity, the QA workload of 81.2% of plans in TV and 62.5% of plans in CV could be reduced by RF.

CONCLUSIONS

The PL model could accurately predict GPR for most VMAT plans. The RF model with 100% sensitivity was preferred for QA results classification. Machine learning can be a useful tool to assist VMAT QA and reduce QA workload.

Dosimetric impact and detectability of multi-leaf collimator positioning errors on Varian Halcyon

The purpose of this study is to investigate the dosimetric impact of multi‐leaf collimator (MLC) positioning errors on a Varian Halcyon for both random and systematic errors, and to evaluate the effectiveness of portal dosimetry quality assurance in catching clinically significant changes caused by these errors. Both random and systematic errors were purposely added to 11 physician‐approved head and neck volumetric modulated arc therapy (VMAT) treatment plans, yielding a total of 99 unique plans. Plans were then delivered on a preclinical Varian Halcyon linear accelerator and the fluence was captured by an opposed portal dosimeter. When comparing dose–volume histogram (DVH) values of plans with introduced MLC errors to known good plans, clinically significant changes to target structures quickly emerged for plans with systematic errors, while random errors caused less change. For both error types, the magnitude of clinically significant changes increased as error size increased. Portal dosimetry was able to detect all systematic errors, while random errors of ±5 mm or less were unlikely to be detected. Best detection of clinically significant errors, while minimizing false positives, was achieved by following the recommendations of AAPM TG‐218. Furthermore, high‐ to moderate correlation was found between dose DVH metrics for normal tissues surrounding the target and portal dosimetry pass rates. Therefore, it may be concluded that portal dosimetry on the Halcyon is robust enough to detect errors in MLC positioning before they introduce clinically significant changes to VMAT treatment plans.

Prediction of dosimetric accuracy for VMAT plans using plan complexity parameters via machine learning

PURPOSE

The dosimetric accuracies of volumetric modulated arc therapy (VMAT) plans were predicted using plan complexity parameters via machine learning.

METHODS

The dataset consisted of 600 cases of clinical VMAT plans from a single institution. The predictor variables (n = 28) for each plan included complexity parameters, machine type, and photon beam energy. Dosimetric measurements were performed using a helical diode array (ArcCHECK), and the dosimetric accuracy of the passing rates for a 5% dose difference (DD5%) and gamma index of 3%/3 mm (γ3%/3 mm) were predicted using three machine learning models: regression tree analysis (RTA), multiple regression analysis (MRA), and neural networks (NNs). First, the prediction models were applied to 500 cases of the VMAT plans. Then, the dosimetric accuracy was predicted using each model for the remaining 100 cases (evaluation dataset). The error between the predicted and measured passing rates was evaluated.

RESULTS

For the 600 cases, the mean ± standard deviation of the measured passing rates was 92.3% ± 9.1% and 96.8% ± 3.1% for DD5% and γ3%/3 mm, respectively. For the evaluation dataset, the mean ± standard deviation of the prediction errors for DD5% and γ3%/3 mm was 0.5% ± 3.0% and 0.6% ± 2.4% for RTA, 0.0% ± 2.9% and 0.5% ± 2.4% for MRA, and -0.2% ± 2.7% and -0.2% ± 2.1% for NN, respectively.

CONCLUSIONS

NNs performed slightly better than RTA and MRA in terms of prediction error. These findings may contribute to increasing the efficiency of patient-specific quality-assurance procedures.

Use of metrics to quantify IMRT and VMAT treatment plan complexity: A systematic review and perspectives

PURPOSE

Fixed-field intensity modulated radiation therapy (FF-IMRT) or volumetric modulated arc therapy (VMAT) beams complexity is due to fluence fluctuation. Pre-treatment Quality Assurance (PTQA) failure could be linked to it. Several plan complexity metrics (PCM) have been published to quantify this complexity but in a heterogeneous formalism. This review proposes to gather different PCM and to discuss their eventual PTQA failure identifier abilities.

METHODS AND MATERIALS

A systematic literature search and outcome extraction from MEDLINE/PubMed (National Center for Biotechnology Information, NCBI) was performed. First, a list and a synthesis of available PCM is made in a homogeneous formalism. Second, main results relying on the link between PCM and PTQA results but also on other uses are listed.

RESULTS

A total of 163 studies were identified and n = 19 were selected after inclusion and exclusion criteria application. Difference is made between fluence and degree of freedom (DOF)-based PCM. Results about the PCM potential as PTQA failure identifier are described and synthesized. Others uses are also found in quality, big data, machine learning and audit procedure.

CONCLUSIONS

A state of the art is made thanks to this homogeneous PCM classification. For now, PCM should be seen as a planning procedure quality indicator although PTQA failure identifier results are mitigated. However limited clinical use seems possible for some cases. Yet, addressing the general PTQA failure prediction case could be possible with the big data or machine learning help.

On the need for tuning the dosimetric leaf gap for stereotactic treatment plans in the Eclipse treatment planning system

The dosimetric leaf gap (DLG) and tongue-and-groove (T&G) effects are critical aspects in the modeling of multileaf collimators (MLC) in the treatment planning system (TPS). In this study, we investigated the dosimetric impact of limitations associated with the T&G modeling in stereotactic plans and its relationship with the need for tuning the DLG in the Eclipse TPS. Measurements were carried out using Varian TrueBeam STx systems from two different institutions. Test fields presenting MLC patterns with several MLC gap sizes (meanGap) and different amounts of T&G effect (TGi) were first evaluated. Secondly, dynamic conformal arc (DCA) and volumetric modulated arc therapy (VMAT) deliveries of stereotactic cases were analyzed in terms of meanGap and TGi. Two DLG values were used in the TPS: the measured DLG (DLGmeas ) and an optimal DLG (DLGopt ). Measured and calculated doses were compared according to dose differences and gamma passing rates (GPR) with strict local gamma criteria of 2%/2 mm. The discrepancies were analyzed for DLGmeas and DLGopt , and their relationships with both TGi and meanGap were investigated. DCA arcs involved significantly lower TGi and larger meanGap than VMAT arcs (P < 0.0001). By using DLGmeas in the TPS, the dose discrepancies increased as TGi increased and meanGap decreased for both test fields and clinical plans. Dose discrepancies dramatically increased with the ratio TGi/meanGap. Adjusting the DLG value was then required to achieve acceptable calculations and configuring the TPS with DLGopt led to an excellent agreement with median GPRs (2%/2 mm) > 99% for both institutions. We also showed that DLGopt could be obtained from the results of the test fields. We demonstrated that the need for tuning the DLG is due to the limitations of the T&G modeling in the Eclipse TPS. A set of sweeping gap tests modified to incorporate T&G effects can be used to determine the optimal DLG value.

Transit and non-transit 3D EPID dosimetry versus detector arrays for patient specific QA

PURPOSE

Despite their availability and simplicity of use, Electronic Portal Imaging Devices (EPIDs) have not yet replaced detector arrays for patient specific QA in 3D. The purpose of this study is to perform a large scale dosimetric evaluation of transit and non-transit EPID dosimetry against absolute dose measurements in 3D.

METHODS

After evaluating basic dosimetric characteristics of the EPID and two detector arrays (Octavius 1500 and Octavius 1000SRS ), 3D dose distributions for 68 VMAT arcs, and 10 IMRT plans were reconstructed within the same phantom geometry using transit EPID dosimetry, non-transit EPID dosimetry, and the Octavius 4D system. The reconstructed 3D dose distributions were directly compared by γ-analysis (2L2 = 2% local/2 mm and 3G2 = 3% global/2 mm, 50% isodose) and by the percentage difference in median dose to the high dose volume (%∆HDVD 50 ).

RESULTS

Regarding dose rate dependency, dose linearity, and field size dependence, the agreement between EPID dosimetry and the two detector arrays was found to be within 1.0%. In the 2L2 γ-comparison with Octavius 4D dose distributions, the average γ-pass rate value was 92.2 ± 5.2%(1SD) and 94.1 ± 4.3%(1SD) for transit and non-transit EPID dosimetry, respectively. 3G2 γ-pass rate values were higher than 95% in 150/156 cases. %∆HDVD 50 values were within 2% in 134/156 cases and within 3% in 155/156 cases. With regard to the clinical classification of alerts, 97.5% of the treatments were equally classified by EPID dosimetry and Octavius 4D.

CONCLUSION

Transit and non-transit EPID dosimetry are equivalent in dosimetric terms to conventional detector arrays for patient specific QA. Non-transit 3D EPID dosimetry can be readily used for pre-treatment patient specific QA of IMRT and VMAT, eliminating the need of phantom positioning.

Impact of plan parameters and modulation indices on patient-specific QA results for standard and stereotactic VMAT

PURPOSE

To demonstrate the impact of modulation indices and plan parameters on the gamma passing rates (GPR) of patient-specific quality assurance of standard and stereotactic volumetric modulated arc therapy (VMAT) plans.

METHODS

A total of 758 patients' QA plans were utilized, including standard VMAT plans with Trilogy (n = 87, group A) and TreuBeam STx (n = 332, group B), and 339 stereotactic VMAT plans with TrueBeam STx (group C). Modulation indices were obtained considering the speed and acceleration of the multileaf collimator (MLC) (MI_s, MI_a), and MLC, gantry speed, and dose rate changes (MI_t). The mean aperture size (MA), monitor unit (MU), and amount of jaw tracking (%JT) were acquired. Gamma analysis was performed with 2 mm/2% and 1 mm/2% for the standard and stereotactic VMAT plans, respectively. Statistical analyses were performed to investigate the correlation between modulation index/plan parameters and GPR.

RESULTS

Spearman's rank correlation to GPRs with MI_s, MI_a, and MI_t, were -0.44, -0.45, and -0.46 for group A; -0.39, -0.37, and -0.38 for group B; and -0.04, -0.11, and -0.10 for group C, respectively. While MU and MA showed significant correlations in all groups, %JT showed a significant correlation only with stereotactic VMAT plans. The most influential parameter combinations were MU-MA (r_s = 0.50), MI_s-%JT (r_s = 0.43), and MU-%JT (r_s = 0.38) for groups A, B, and C, respectively.

CONCLUSIONS

MLC modulation mostly affected the GPR in the delivery of standard VMAT plans, while MU and %JT showed more importance in stereotactic VMAT plans.

Evaluation of the Impact of the Linac MLC and Gantry Sag in volumetric modulated arc therapy

PURPOSE

Mechanical sag in the radiotherapy linear accelerator gantry and multi-leaf collimator (MLC) carriage effectively causes systematic deviations in the isocenter with respect to gantry angle. To minimize the impact of this error on treatment, a tolerance value of a 1-mm mechanical isocenter shift is commonly accepted for intensity-modulated radiation therapy quality assurance (QA). However, this tolerance value has not been firmly established for volumetric modulated arc therapy (VMAT) treatments. The purpose of this study is therefore to evaluate the impact of gantry and MLC carriage sag on VMAT clinical performance.

METHODS

A published dataset of Elekta and Varian sag measurements served as a starting point for the investigation. Typical sag profiles were chosen and modeled as continuous isocenter deviations in three dimensions. The data were then incorporated into existing Digital Imaging and Communications in Medicine protocol, extended for radiotherapy plans via a "beam-splitting" algorithm. Three treatment sites were investigated in parallel: head and neck, prostate, and prostate with surrounding lymph nodes. Monte Carlo-simulated dose distributions were obtained for varying magnifications of the modeled sag. The resulting dose distributions, including that for no error, were compared qualitatively and quantitatively, against multiple metrics.

RESULTS

The dose-volume histograms (DVHs) for all plans exhibited a decrease in planning target volume (PTV) dose uniformity with increasing sag magnification, whereas dose to organs at risk exhibited no coherent trend. The prostate plan was shown to be the most vulnerable to mechanical sag across all considered metrics. However, all plans with peak isocenter deviation less than 1 mm were well within typical cutoff points for each metric.

CONCLUSIONS

All avenues of investigation presented substantiate the commonly accepted tolerance value of a 1-mm peak isocenter shift in annual linac QA.

Modulation indices and plan delivery accuracy of volumetric modulated arc therapy

PURPOSE

We evaluated the performance of various modulation indices (MI) for volumetric modulated arc therapy (VMAT) to predict plan delivery accuracy.

METHODS

The specific indices evaluated were MI quantifying the mechanical uncertainty (MI_t ), MI quantifying the mechanical and dose calculation uncertainties (MI_c ), MI for station parameter optimized radiation therapy (MI_SPORT ), modulation complexity score for VMAT (MCS_v ), leaf travel modulation complexity score (LTMCS), plan averaged beam area (PA), plan averaged beam irregularity (PI), plan averaged beam modulation (PM), and plan normalized monitor unit (PMU) to predict VMAT delivery accuracy. By utilizing 240 VMAT plans generated with the Trilogy and TrueBeam STx, Spearman's rank correlation coefficients (r) were calculated between the MIs and measures of conventional methods.

RESULTS

For the Trilogy system, MI_c showed the highest r values with gamma passing rates (GPRs) (r = -0.624 with P < 0.001 for MapCHECK2 and r = -0.655 with P < 0.001 for ArcCHECK). For TrueBeam STx, MI_c also showed the highest r values with GPRs (r = -0.625 with P < 0.001 for the MapCHECK2 and r = -0.561 with P < 0.001 for the ArcCHECK). The MI_t and MI_c showed the highest r values to the MLC position errors for the Trilogy and TrueBeam STx systems (r = 0.770 with P < 0.001 and r = 0.712 with P < 0.001, respectively). The PA showed the highest percent of r values (P < 0.05) to differences in the dose-volume parameters between original VMAT plans and actual deliveries for the Trilogy systems (30.9%). Both the MI_t and MI_c showed the highest percent of r values (P < 0.05) to differences in the dose-volume parameters between original VMAT plans and actual deliveries for the TrueBeam STx systems (31.8%).

CONCLUSION

To comprehensively review the results, the MI_c showed the best performance to predict the VMAT delivery accuracy.

Predicting VMAT patient-specific QA results using a support vector classifier trained on treatment plan characteristics and linac QC metrics

The use of treatment plan characteristics to predict patient-specific quality assurance (QA) measurement results has recently been reported as a strategy to help facilitate automated pre-treatment verification workflows or to provide a virtual assessment of delivery quality. The goal of this work is to investigate the potential of using treatment plan characteristics and linac performance metrics (i.e. quality control test results) in combination with machine learning techniques to predict the results of VMAT patient-specific QA measurements. Using features that describe treatment plan complexity and linac performance metrics, we trained a linear support vector classifier (SVC) to classify the results of VMAT patient-specific QA measurements. The 'targets' in this model were simple classes representing median dose difference between measured and expected dose distributions-'hot' if the median dose deviation was  >1%, 'cold' if it was  <-1%, and 'normal' if it was within  ±1%. A total of 1620 unique patient-specific QA measurements were available for model development and testing. 75% of the data were used to develop and cross-validate the model, and the remaining 25% were used for an independent assessment of model performance. For the model development phase, a recursive feature elimination (RFE) cross-validation technique was used to eliminate unimportant features. Model performance was assessed using receiver operator characteristic (ROC) curve metrics. Of the ten features found to be most predictive of patient-specific QA measurement results, half were derived from treatment plan characteristics and half from quality control (QC) metrics characterizing linac performance. The model achieved a micro-averaged area under the ROC curve of 0.93, and a macro-averaged area under the ROC curve of 0.88. This work demonstrates the potential of using both treatment plan characteristics and routine linac QC results in the development of machine learning models for VMAT patient-specific QA measurements.

National survey of patient specific IMRT quality assurance in China

BACKGROUND

To analyze and present the China's national survey on patient-specific IMRT quality assurance (QA).

METHODS

A national survey was conducted in all radiotherapy centers in China to collect comprehensive information on status of IMRT QA practice, including machine, technique, equipment, issues and suggestions.

RESULTS

Four hundred and three centers responded to this survey, accounting for 56.92% of all the centers implementing IMRT in China. The total number of medical physicists and the total number of patients treated with IMRT annually in these centers was 1599 and 305,000 respectively. All centers implemented measurement-based verification. Point dose verification and 2D dose verification was implemented in 331 and 399 centers, respectively. Three hundred forty-eight centers had 2D arrays, and 52 centers had detector devices designed to measure VMAT beams. EPID and film were used in 78 and 70 centers, respectively. Seventeen and 20 centers used log file and 3D DVH analysis, respectively. One hundred sixty-eight centers performed measurement-based verification not for each patient based on different selection criteria. The techniques and methods varied significantly in both point dose and dose distribution verification, from evaluation metrics, criteria, tolerance limit, and steps to check failed IMRT QA plans. Major issues identified in this survey were the limited resources of physicists, QA devices, and linacs.

CONCLUSIONS

IMRT QA was implemented in all the surveyed centers. The practice of IMRT QA varied significantly between centers. An increase in personnel, QA devices and linacs is highly desired. National standard, guideline, regulation and training programs are urgently needed in China for consistent and effective implementation of IMRT QA.

Assessing the feasibility of single target radiosurgery quality assurance with portal dosimetry

Purpose

To assess the feasibility of using portal dosimetry (PD) for pre‐treatment quality assurance of single target, flattening filter free (FFF), volumetric arc therapy intracranial radiosurgery plans.

Methods

A PD algorithm was created for a 10X FFF beam on a Varian Edge linear accelerator (Varian Inc, Palo Alto, CA, USA). Treatment plans that were previously evaluated with Gafchromic EBT‐XD (Ashland, Bridgewater, NJ, USA) film were measured via PD and analyzed with the ARIA Portal Dosimetry workspace. Absolute dose evaluation for film and PD was done by computing the mean dose in the region receiving greater than or equal to 90% of the max dose and comparing to the mean dose in the same region calculated by the treatment planning system (TPS). Gamma analysis with 10% threshold and 3%/2 mm passing criteria was performed on film and portal images.

Results

Thirty‐six PD verification plans were delivered and analyzed. The average PD to TPS dose was 0.989 ± 0.01 while film to TPS dose was 1.026 ± 0.01. All PD plans passed the gamma analysis with 100% of points having gamma <1. Overall, PD to TPS dose agreement was found to be target size dependent. As target size decreases, PD to TPS dose ratio decreased from 1.004 for targets with diameters between 15–31 mm and 0.978 for targets with diameters less than 15 mm.

Conclusion

The agreement of PD to TPS mean dose in the high dose region was found to be dependent on target size. Film measurements did not exhibit size dependence. All PD plans passed the 3%/2 mm gamma analysis, but caution should be used when using PD to assess overall dosimetric accuracy of the treatment plan for small targets.

Distributive quality assurance and delivery of stereotactic ablative radiotherapy treatments amongst beam matched linear accelerators: A feasibility study

PURPOSE

Beam matching occurs on all linacs to some degree and when two are more are matched to each other, patients are able to be transferred between machines. Quality assurance of plans can also be performed "distributively" on any of the matched linacs. The degree to which machines are matched and how this translates to like delivery of plans has been the focus of a number of studies. This concept has not yet been explored for stereotactic techniques which require a higher degree of accuracy. This study proposes beam matching criteria which allows for the distributive delivery and quality assurance of stereotactic body radiotherapy (SBRT) plans.

METHOD

Two clinically relevant and complex volumetric modulated arc therapy (VMAT) SBRT spine and lung plans were chosen as benchmarking cases. These were delivered on nine previously beam matched linacs with quality assurance performed through ArcCheck and film exposure in the sagittal plane. Measured doses were compared to their treatment planning system predictions through gamma analysis at a range of criteria.

RESULTS

Despite differences in beam match parameters and variations in small fields, all nine linacs produced accurate deliveries with a tight deviation in the population sample. Pass rates were well above suggested tolerances at the recommended gamma criterion. Film was able to detect dose errors to a greater degree than ArcCheck.

CONCLUSION

Distributive quality assurance and delivery of stereotactic ablative radiotherapy treatments amongst beam matched linacs is certainly feasible provided the linacs are matched to a strict protocol like that suggested in this study and regular quality assurance is performed on the matched fleet. Distributive quality assurance and delivery of SBRT provides the possibility of efficiency gains for physicists as well as treatment staff.

Commissioning of the Mobius3D independent dose verification system for TomoTherapy

In radiation therapy, a secondary independent dose verification is an important component of a quality control system. Mobius3D calculates three‐dimensional (3D) patient dose using reference beam data and a collapsed cone convolution algorithm and analyzes dose‐volume histogram automatically. There are currently no published data on commissioning and determining tolerance levels of Mobius3D for TomoTherapy. To verify the calculation accuracy and adjust the parameters of this system, we compared the measured dose using an ion chamber and film in a phantom with the dose calculated using Mobius3D for nine helical intensity‐modulated radiation therapy plans, each with three nominal field widths. We also compared 126 treatment plans used in our institution to treat prostate, head‐and‐neck, and esophagus tumors based on dose calculations by treatment planning system for given dose indices and 3D gamma passing rates with those produced by Mobius3D. On the basis of these results, we showed that the action and tolerance levels at the average dose for the planning target volume (PTV) at each treatment site are at μ ± 2σ and μ ± 3σ, respectively. After adjusting parameters, the dose difference ratio on average was −0.2 ± 0.6% using ion chamber and gamma passing rate with the criteria of 3% and 3 mm on average was 98.8 ± 1.4% using film. We also established action and tolerance levels for the PTV at the prostate, head‐and‐neck, esophagus, and for the organ at risk at all treatment sites. Mobius3D calculations thus provide an accurate secondary dose verification system that can be commissioned easily and immediately after installation. Before clinical use, the Mobius3D system needs to be commissioned using the treatment plans for patients treated in each institution to determine the calculational accuracy and establish tolerances for each treatment site and dose index.

Comparison of gamma index based on dosimetric error and clinically relevant dose-volume index based on three-dimensional dose prediction in breast intensity-modulated radiation therapy

Background

Measurement-guided dose reconstruction has lately attracted significant attention because it can predict the delivered patient dose distribution. Although the treatment planning system (TPS) uses sophisticated algorithm to calculate the dose distribution, the calculation accuracy depends on the particular TPS used. This study aimed to investigate the relationship between the gamma passing rate (GPR) and the clinically relevant dose–volume index based on the predicted 3D patient dose distribution derived from two TPSs (XiO, RayStation).

Methods

Twenty-one breast intensity-modulated radiation therapy plans were inversely optimized using XiO. With the same plans, both TPSs calculated the planned dose distribution. We conducted per-beam measurements on the coronal plane using a 2D array detector and analyzed the difference in 2D GPRs between the measured and planned doses by commercial software. Using in-house software, we calculated the predicted 3D patient dose distribution and derived the predicted 3D GPR, the predicted per-organ 3D GPR, and the predicted clinically relevant dose–volume indices [dose–volume histogram metrics and the value of the tumor-control probability/normal tissue complication probability of the planning target volume and organs at risk]. The results derived from XiO were compared with those from RayStation.

Results

While the mean 2D GPRs derived from both TPSs were 98.1% (XiO) and 100% (RayStation), the mean predicted 3D GPRs of ipsilateral lung (73.3% [XiO] and 85.9% [RayStation]; p < 0.001) had no correlation with 2D GPRs under the 3% global/3 mm criterion. Besides, this significant difference in terms of referenced TPS between XiO and RayStation could be explained by the fact that the error of predicted V_5Gy of ipsilateral lung derived from XiO (29.6%) was significantly larger than that derived from RayStation (− 0.2%; p < 0.001).

Conclusions

GPR is useful as a patient quality assurance to detect dosimetric errors; however, it does not necessarily contain detailed information on errors. Using the predicted clinically relevant dose–volume indices, the clinical interpretation of dosimetric errors can be obtained. We conclude that a clinically relevant dose–volume index based on the predicted 3D patient dose distribution could add to the clinical and biological considerations in the GPR, if we can guarantee the dose calculation accuracy of referenced TPS.

Epid-based in vivo dose verification for lung stereotactic treatments delivered with multiple breath-hold segmented volumetric modulated arc therapy

We evaluated an EPID-based in-vivo dosimetry (IVD) method for the dose verification and the treatment reproducibility of lung SBRT-VMAT treatments in clinical routine. Ten patients with lung metastases treated with Elekta VMAT technique were enrolled. All patients were irradiated in five consecutive fractions, with total doses of 50 Gy. Set-up was carried out with the Elekta stereotactic body frame. Eight patients were simulated and treated using the Active Breath Control (ABC) system, a spirometer enabling patients to maintain a breath-hold at a predetermined lung volume. Two patients were simulated and treated in free-breathing using an abdominal compressor. IVD was performed using the SOFTDISO software. IVD tests were evaluated by means of (a) ratio R between daily in-vivo isocenter dose and planned dose and (b) γ-analysis between EPID integral portal images in terms of percentage of points with γ-value smaller than one (γ_% ) and mean γ-values (γ_mean ) using a 3%(global)/3 mm criteria. Alert criteria of ±5% for R ratio, γ_%  < 90%, and γ_mean  > 0.67 were chosen. 50 transit EPID images were acquired. For the patients treated with ABC spirometer, the results reported a high level of accuracy in dose delivery with 100% of tests within ±5%. The γ-analysis showed a mean value of γ_mean equal to 0.21 (range: 0.04-0.56) and a mean γ_% equal to 96.9 (range: 78-100). Relevant discrepancies were observed only for the two patients treated without ABC, mainly due to a blurring dose effect due to residual respiratory motion. Our method provided a fast and accurate procedure in clinical routine for verifying delivered dose as well as for detecting errors.