A brief look at model-based dose calculation principles, practicalities, and promise

Being certain about uncertainties: a robust evaluation method for high-dose-rate prostate brachytherapy treatment plans including the combination of uncertainties

In high-dose-rate (HDR) prostate brachytherapy the combined effect of uncertainties cause a range of possible dose distributions deviating from the nominal plan, and which are not considered during treatment plan evaluation. This could lead to dosimetric misses for critical structures and overdosing of organs at risk. A robust evaluation method to assess the combination of uncertainties during plan evaluation is presented and demonstrated on one HDR prostate ultrasound treatment plan retrospectively. A range of uncertainty scenarios are simulated by changing six parameters in the nominal plan and calculating the corresponding dose distribution. Two methods are employed to change the parameters, a probabilistic approach using random number sampling to evaluate the likelihood of variation in dose distributions, and a combination of the most extreme possible values to access the worst-case dosimetric outcomes. One thousand probabilistic scenarios were run on the single treatment plan with 43.2% of scenarios passing seven of the eight clinical objectives. The prostate D90 had a standard deviation of 4.4%, with the worst case decreasing the dose by up to 27.2%. The urethra D10 was up to 29.3% higher than planned in the worst case. All DVH metrics in the probabilistic scenarios were found to be within acceptable clinical constraints for the plan under statistical tests for significance. The clinical significance of the results from the robust evaluation method presented on any individual treatment plan needs to be compared in the context of a historical data set that contains patient outcomes with robustness analysis data to ascertain a baseline acceptance.

AAPM WGDCAB Report 372: A joint AAPM, ESTRO, ABG, and ABS report on commissioning of model-based dose calculation algorithms in brachytherapy

The introduction of model-based dose calculation algorithms (MBDCAs) in brachytherapy provides an opportunity for a more accurate dose calculation and opens the possibility for novel, innovative treatment modalities. The joint AAPM, ESTRO, and ABG Task Group 186 (TG-186) report provided guidance to early adopters. However, the commissioning aspect of these algorithms was described only in general terms with no quantitative goals. This report, from the Working Group on Model-Based Dose Calculation Algorithms in Brachytherapy, introduced a field-tested approach to MBDCA commissioning. It is based on a set of well-characterized test cases for which reference Monte Carlo (MC) and vendor-specific MBDCA dose distributions are available in a Digital Imaging and Communications in Medicine-Radiotherapy (DICOM-RT) format to the clinical users. The key elements of the TG-186 commissioning workflow are now described in detail, and quantitative goals are provided. This approach leverages the well-known Brachytherapy Source Registry jointly managed by the AAPM and the Imaging and Radiation Oncology Core (IROC) Houston Quality Assurance Center (with associated links at ESTRO) to provide open access to test cases as well as step-by-step user guides. While the current report is limited to the two most widely commercially available MBDCAs and only for 192 Ir-based afterloading brachytherapy at this time, this report establishes a general framework that can easily be extended to other brachytherapy MBDCAs and brachytherapy sources. The AAPM, ESTRO, ABG, and ABS recommend that clinical medical physicists implement the workflow presented in this report to validate both the basic and the advanced dose calculation features of their commercial MBDCAs. Recommendations are also given to vendors to integrate advanced analysis tools into their brachytherapy treatment planning system to facilitate extensive dose comparisons. The use of the test cases for research and educational purposes is further encouraged.

Accuracy of Acuros[Formula: see text] BV as determined from GATE monte-carlo simulation

The American Association of Physicists in Medicine's Task Group No.43 has provided a standardised dose calculation methodology that is now the international benchmark for all brachytherapy dosimetry publications and treatment planning systems. However, limitations in this methodology has seen the development of Model-Based Dose Calculation Algorithms (MBDCA). In 2009, Varian (Varian Medical Systems, Palo Alto, CA, USA) released Acuros[Formula: see text] BrachyVision (ABV) which calculates dose by explicitly solving the Linear Boltzmann Transport Equation. In this study we have assessed the accuracy of ABV dose calculations within a range of materials relevant to high dose rate brachytherapy with an iridium-192 ([Formula: see text]Ir) source. Accuracy assessment has been achieved by implementing a modelled GamaMed Plus [Formula: see text]Ir source within a series of phantoms using the GEANT4 Application for Emission Tomography (GATE) to calculate dose for comparison with dose as determined by ABV. Comparisons between GATE and ABV were made using point-to-point profile comparisons and 1D gamma analysis. Source validation results yielded good agreement with published data. Spectrum and TG43U1 comparisons showed no major differences, with TG43U1 comparisons agreeing within ± 1%. Point-to-point comparisons showed large differences between GATE and ABV near the source and in low density materials. 1D gamma analysis pass criteria of 2%/1 mm and 2%/2 mm yielded passing rates ranging between 51.72-100% and 62.07-100% respectively. A critical analysis of this study's results suggest that ABV is unable to accurately calculate doses in low density materials. Furthermore, spatial accuracy of dose near the source is within 2 mm.

A review of dosimetric impact of implementation of model-based dose calculation algorithms (MBDCAs) for HDR brachytherapy

To obtain dose distributions more physically representative to the patient anatomy in brachytherapy, calculation algorithms that can account for heterogeneity are required. The current standard AAPM Task Group No 43 (TG-43) dose calculation formalism has some clinically relevant dosimetric limitations. Lack of tissue heterogeneity and scattered dose corrections are the major weaknesses of the TG-43 formalism and could lead to systematic dose errors in target volumes and organs at risk. Over the last decade, model-based dose calculation algorithms (MBDCAs) have been clinically offered as complementary algorithms beyond the TG43 formalism for high dose rate (HDR) brachytherapy treatment planning. These algorithms provide enhanced dose calculation accuracy by using the information in the patient's computed tomography images, which allows modeling the patient's geometry, material compositions, and the treatment applicator. Several researchers have investigated the implementation of MBDCAs in HDR brachytherapy for dose optimization, but moving toward using them as primary algorithms for dose calculations is still lagging. Therefore, an overview of up-to-date research is needed to familiarize clinicians with the current status of the MBDCAs for different cancers in HDR brachytherapy. In this paper, we review the MBDCAs for HDR brachytherapy from a dosimetric perspective. Treatment sites covered include breast, gynecological, lung, head and neck, esophagus, liver, prostate, and skin cancers. Moreover, we discuss the current status of implementation of MBDCAs and the challenges.

Comprehensive methodology for commissioning modern 3D-image-based treatment planning systems for high dose rate gynaecological brachytherapy: A review

PURPOSE

Correct commissioning of treatment planning systems (TPSs) is important for reducing treatment failure events. There is currently no comprehensive and robust methodology available for TPS commissioning in modern brachytherapy. This review aimed to develop a comprehensive template for commissioning modern 3D-image-based brachytherapy TPSs for high dose rate (HDR) gynaecological applications.

METHODS

The literature relevant to TPS commissioning, including both external beam radiation therapy (EBRT) and brachytherapy, as well as guidelines by the International Atomic Energy Agency (IAEA), the American Association of Physicists in Medicine (AAPM), and the European Society for Radiotherapy and Oncology (ESTRO) were searched, studied and appraised. The applied relevant EBRT TPS commissioning tests were applied to brachytherapy. The developed template aimed to cover all dosimetric and non-dosimetric issues.

RESULTS

The essential commissioning items could be categorized into six parts: geometry, dose calculation, plan evaluation tools, plan optimization, TPS output, and end-to-end verification. The final template consists of 43 items. This paper presents the purpose and role of each test, as well as tolerance limits, to facilitate the use of the template.

CONCLUSION

The information and recommendations available in a collection of publications over many years have been reviewed in order to develop a comprehensive template for commissioning complex modern 3D-image-based brachytherapy TPSs for HDR gynaecological applications. The up-to-date and concise information contained in the template can aid brachytherapy physicists during TPS commissioning as well as devising a regular quality assurance program and allocation of time and resources.

Model-Based Dose Calculation Algorithms for Brachytherapy Dosimetry

The purpose of this study was to review the limitations of dose calculation formalisms for photon-emitting brachytherapy sources based on the American Association of Physicists in Medicine (AAPM) Task Group No. 43 (TG-43) report and to provide recommendations to transition to model-based dose calculation algorithms. Additionally, an overview of these algorithms and approaches is presented. The influence of tissue and seed/applicator heterogeneities on brachytherapy dose distributions for breast, gynecologic, head and neck, rectum, and prostate cancers as well as eye plaques and electronic brachytherapy treatments were investigated by comparing dose calculations based on the TG-43 formalism and model-based dose calculation algorithms.

Deployment and performance of model-based dose calculation algorithm in 192Ir shielded cylinder brachytherapy

PURPOSE

The purpose of the study was to integrate model-based dose calculation algorithm (MBDCA) and 3-D planning into our institutional capabilities for clinical cases with inherent heterogeneities, namely shielded cylinder cases, which were previously performed using templates.

METHODS AND MATERIALS

AcurosBV (Varian Medical Systems) was selected as MBDCA, and we compared results against the American Association of Physicists in Medicine working group (WG) reference Test Case #4, which contains a 36-mm-diameter 180-degree shielded cylinder. The last five clinically used template plans, as generated with ABACUS (Varian Medical Systems), were compared against MBDCA results. Clinical plans used 20-, 23-, and 26-mm-diameter cylinders, prescribed from 5 to 7 Gy, 50- to 110-mm active length, 7 to 21 dwell positions, with 5- or 10-mm spacing.

RESULTS

Compared with the WG reference plan, AcurosBV produced excellent agreement, within 1% at reference points. Larger deviations arose only within the applicator itself. Historical plans generated using ABACUS had higher point dose than AcurosBV by 3-4% or 2-3% using transport within medium at prescription points, with dose to medium or water, respectively.

CONCLUSIONS

To verify the accuracy of our MBDCA algorithm, we benchmarked against the WG data set available for shielded cylinders. We discovered a 3-4% difference in dose from our established historical templates, which is easily outweighed by daily positioning uncertainties. Dose distributions from MBDCA were used to assess and validate existing historical templates.

The effect of density overrides on magnetic resonance-guided radiation therapy planning for lung cancer

Background and Purpose

Inverse treatment planning for lung cancer can be challenging since density heterogeneities may appear inside the planning target volume (PTV). One method to improve the quality of intensity modulation is the override of low density tissues inside the PTV during plan optimization. For magnetic resonance-guided radiation therapy (MRgRT), where the influence of the magnetic field on secondary electrons is sensitive to the tissue density, the reliability of density overrides has not yet been proven. This work, therefore, gains a first insight into density override strategies for MRgRT.

Material and methods

Monte Carlo-based treatment plans for five lung cancer patients were generated based on free-breathing CTs and two density override strategies. Different magnetic field configurations were considered with their effect being accounted for during optimization. Optimized plans were forward calculated to 4D-CTs and accumulated for the comparison of planned and expected delivered dose.

Results

For MRgRT, density overrides led to a discrepancy between the delivered and planned dose. The tumor volume coverage deteriorated for perpendicular magnetic fields of 1.5 T to 93.6% (D_98%). For inline fields a maximal increase of 2.2% was found for the mean dose. In terms of organs at risk, a maximal sparing of 0.6 Gy and 0.9 Gy was observed for lung and heart, respectively.

Conclusions

In this work, first results on the effect of density overrides on treatment planning for MRgRT are presented. It was observed that the underestimation of magnetic field effects in overridden densities during treatment planning resulted in an altered delivered dose, depending on the field strength and orientation.

A comparative assessment of inhomogeneity and finite patient dimension effects in 60Co and 192Ir high-dose-rate brachytherapy

Purpose

To perform a comparative study of heterogeneities and finite patient dimension effects in ⁶⁰Co and ¹⁹²Ir high-dose-rate (HDR) brachytherapy.

Material and methods

Clinically equivalent plans were prepared for 19 cases (8 breast, 5 esophagus, 6 gynecologic) using the Ir2.A85-2 and the Co0.A86 HDR sources, with a TG-43 based treatment planning system (TPS). Phase space files were obtained for the two source designs using MCNP6, and validated through comparison to a single source dosimetry results in the literature. Dose to water, taking into account the patient specific anatomy and materials (D_w,m), was calculated for all plans using MCNP6, with input files prepared using the BrachyGuide software tool to analyze information from DICOM RT plan exports.

Results

A general TG-43 dose overestimation was observed, except for the lungs, with a greater magnitude for ¹⁹²Ir. The distribution of percentage differences between TG-43 and Monte Carlo (MC) in dose volume histogram (DVH) indices for the planning target volume (PTV) presented small median values (about 2%) for both ⁶⁰Co and ¹⁹²Ir, with a greater dispersion for ¹⁹²Ir. Regarding the organs at risk (OARs), median percentage differences for breast V_50% were 3% (5%) for ⁶⁰Co (¹⁹²Ir). Differences in median skin D_2cc were found comparable, with a larger dispersion for ¹⁹²Ir, and the same applied to the lung D_10cc and the aorta D_2cc. TG-43 overestimates D_2cc for the rectum and the sigmoid, with median differences from MC within 2% and a greater dispersion for ¹⁹²Ir. For the bladder, the median of the difference is greater for ⁶⁰Co (~2%) than for ¹⁹²Ir (~0.75%), demonstrating however a greater dispersion again for ¹⁹²Ir.

Conclusions

The magnitude of differences observed between TG-43 based and MC dosimetry and their smaller dispersion relative to ¹⁹²Ir, suggest that ⁶⁰Co HDR sources are more amenable to the TG-43 assumptions in clinical treatment planning dosimetry.

Initial evaluation of Advanced Collapsed cone Engine dose calculations in water medium for I-125 seeds and COMS eye plaques

PURPOSE

To investigate the dose calculation accuracy in water medium of the Advanced Collapsed cone Engine (ACE) for three sizes of COMS eye plaques loaded with low-energy I-125 seeds.

METHODS

A model of the Oncura 6711 I-125 seed was created for use with ACE in Oncentra^® Brachy (OcB) using primary-scatter separated (PSS) point dose kernel and Task Group (TG) 43 datasets. COMS eye plaque models of diameters 12, 16, and 20 mm were introduced into the OcB applicator library based on 3D CAD drawings of the plaques and Silastic inserts. To perform TG-186 level 1 commissioning, treatment plans were created in OcB for a single source in water and for each COMS plaque in water for two scenarios: with only one centrally loaded seed, or with all seed positions loaded. ACE dose calculations were performed in high accuracy mode with a 0.5 × 0.5 × 0.5 mm³ calculation grid. The resulting dose data were evaluated against Monte Carlo (MC) calculated doses obtained with MCNP6, using both local and global percent differences.

RESULTS

ACE doses around the source for the single seed in water agreed with MC doses on average within < 5% inside a 6 × 6 × 6 cm³ region, and within < 1.5% inside a 2 × 2 × 2 cm³ region. The PSS data were generated at a higher resolution within 2 cm from the source, resulting in this improved agreement closer to the source due to fewer approximations in the ACE dose calculation. Average differences in both investigated plaque loading patterns in front of the plaques and on the plaque central axes were ≤ 2.5%, though larger differences (up to 12%) were found near the plaque lip.

CONCLUSIONS

Overall, good agreement was found between ACE and MC dose calculations for a single I-125 seed and in front of the COMS plaques in water. More complex scenarios need to be investigated to determine how well ACE handles heterogeneous patient materials.