The dosimetric accuracy of a commercial model-based dose calculation algorithm in modeling a six-groove direction modulated brachytherapy tandem applicator

Objective.With advancements in high-dose rate brachytherapy, the clinical translation of intensity modulated brachytherapy (IMBT) innovations necessitates utilization of model-based dose calculation algorithms (MBDCA) for accurate and rapid dose calculations. This study uniquely benchmarks a commercial MBDCA, BrachyVision ACUROSTM(BVA), against Monte Carlo (MC) simulations, evaluating dose distributions for a novel IMBT applicator, termed as thesix-grooveDirection Modulated Brachytherapy (DMBT) tandem, expanding beyond previous focus on partially shielded vaginal cylinder applicators, through a novel methodology.Approach.The DMBT tandem applicator, made of a tungsten alloy with six evenly spaced grooves, was simulated using the GEANT4 MC code. Subsequently, two main scenarios were created using the BVA and reproduced by the MC simulations: 'Source at the Center of the Water Phantom (SACWP)' and 'Source at the Middle of the Applicator (SAMA)' for three cubical virtual water phantoms (20 cm)3, (30 cm)3, and (40 cm)3. A track length estimator was utilized for dose calculation and 2D/3D scoring were performed. The difference in isodose surfaces/lines (i.e. coverage) at each voxel,ΔDIsodose Levels/Lines, was thus calculated for relevant normalization points (rref).Results.The coverage was comparable, based on 2D scoring, between the BVA and MC isodose surfaces/lines for the region of clinical relevance, (i.e. within 5 cm radius from the source) withΔDIsodose Lines(rref: 1 cm from the source) falling within 2% for the two scenarios for all phantom sizes. For the phantom (20 cm)3,ΔDIsodose Levels(3D scoring) recorded the range [-3.0% +6.5%] ([-7.4% +7.3%]) for 95% of the voxels of the same scoring volume for the SACWP (SAMA) scenario.Significance.The results indicated that the BVA could render comparable coverage as compared to the MC simulations in the region of clinical relevance for different phantom sizes.ΔDIsodose Linesmay offer an advantageous metric for evaluation of MBDCAs in clinical setting.

On the impact of improved dose calculation accuracy in clinical treatment planning for superficial high-dose-rate brachytherapy of extensive scalp lesions

TG-43-based dose calculations disregard tissue heterogeneities and finite scatter conditions, prompting the introduction of model-based dose calculation algorithms (MBDCAs) to improve accuracy in high-dose-rate (HDR) brachytherapy. This study evaluated the effectiveness of MBDCAs over TG-43 in HDR 192Ir brachytherapy of extended scalp lesions. Treatment planning dose calculations were compared with Monte Carlo (MC) data. TG-43 exhibited a dose overestimation ranging from 10% to 23% as the distance from the implant increased, while a better agreement from 2% to 6% was observed between the MBDCA and MC, supporting the adoption of MBDCAs for dose calculations in broad scalp lesions.

Investigations on the beam quality correction factor for ionization chambers in high-energy brachytherapy dosimetry

Objective. To enhance the investigations on MC calculated beam quality correction factors of thimble ionization chambers from high-energy brachytherapy sources and to develop reliable reference conditions in source and detector setups in water.Approach. The response of five different ionization chambers from PTW-Freiburg and Standard Imaging was investigated for irradiation by a high dose rate Ir-192 Flexisource in water. For a setup in a Beamscan water phantom, Monte Carlo simulations were performed to calculate correction factors for the chamber readings. After exact positioning of source and detector the absorbed dose rate at the TG-43 reference point at one centimeter nominal distance from the source was measured using these factors and compared to the specification of the calibration certificate. The Monte Carlo calculations were performed using the restricted cema formalism to gain further insight into the chamber response. Calculations were performed for the sensitive volume of the chambers, determined by the methods currently used in investigations of dosimetry in magnetic fields.Main results. Measured dose rates and values from the calibration certificate agreed within the combined uncertainty (k= 2) for all chambers except for one case in which the full air cavity was simulated. The chambers showed a distinct directional dependence. With the restricted cema formalism calculations it was possible to examine volume averaging and energy dependence of the perturbation factors contributing to the beam quality correction factor also differential in energy.Significance. This work determined beam quality correction factors to measure the absorbed dose rate from a brachytherapy source in terms of absorbed dose to water for a variety of ionization chambers. For the accurate dosimetry of brachytherapy sources with ionization chambers it is advisable to use correction factors based on the sensitive volume of the chambers and to take account for the directional dependence of chamber response.

The effect of vaginal cylinder inhomogeneity on the HDR brachytherapy dose calculations using Monte Carlo simulations

PURPOSE

To analytically assess the heterogeneity effect of vaginal cylinders (VC) made of high-density plastics on dose calculations, considering the prescription point (surface or 5 mm beyond the surface), and benchmark the accuracy of a commercial model-based dose calculation (MBDC) algorithm using Monte Carlo (MC) simulations.

METHODS AND MATERIALS

The GEANT4 MC code was used to simulate a commercial 192 Ir HDR source and VC, with diameters ranging from 20 to 35 mm, inside a virtual water phantom. Standard plans were generated from a commercial treatment planning system [TPS-BrachyVision ACUROS (BV)] optimized for a treatment length of 5 cm through two dose calculation approaches: (1) assuming all the environment as water (i.e., Dw,w-MC & Dw,w-TG43 ) and (2) accounting for the heterogeneity of VC applicators (i.e., Dw,w-App-MC & Dw,w-App-MBDC ). The compared isodose lines, and dose & energy difference maps were extracted for analysis. In addition, the dose difference on the peripheral surface, along the applicator and at middle of treatment length, as well as apical tip was evaluated.

RESULTS

The Dw,w-App-MC results indicated that the VC heterogeneity can cause a dose reduction of (up to) % 6.8 on average (for all sizes) on the peripheral surface, translating to 1 mm shrinkage of the isodose lines compared to Dw,w-MC . In addition, the results denoted that BV overestimates the dose on the peripheral surface and apical tip of about 3.7% and 17.9%, respectively, (i.e., Dw,w-App-MBDC vs Dw,w-App-MC ) when prescribing to the surface. However, the difference between the two were negligible at the prescription point when prescribing to 5 mm beyond the surface.

CONCLUSION

The VCs' heterogeneity could cause dose reduction when prescribing dose to the surface of the applicator, and hence increases the level of uncertainty. Thus, reviewing the TG43 results, in addition to ACUROS, becomes prudent, when evaluating the surface coverage at the apex.

Dosimetric Impact of Air Pockets in the Vaginal Cuff Brachytherapy Using Model-based Dose Calculation Algorithm

Background

Endometrial cancer is the most common disease of the female reproductive system. Vaginal cuff brachytherapy (VCB) has intrinsic advantages compared to external beam therapy when treated with radiation. A single-channel cylinder is a standard applicator in VCB. The present study aims to estimate a change in the dose to vaginal mucosa due to air pockets between the cylinder and vaginal mucosa by calculating with the Acuros BV algorithm and comparing it to the Task Group 43 (TG-43) algorithm.

Materials and Methods

Patients who presented with air packets were included retrospectively. For each patient, three plans were created: the first plan used TG-43, the second plan used dose recalculation with Acuros BV, and the third plan was generated by re-optimization by Acuros BV. On the same axial computed tomography image, the point doses at the cylinder's surface and the displaced mucosa were recorded and the ratios were then estimated.

Results

The average volume of air pockets was 0.08 cc (range of 0.01-0.3 cc), and 84% of air pockets displaced the vaginal mucosa by ≥0.2 cm. The average ratios of dose were 0.77 ± 0.09 (1 standard deviation [SD]) and 0.78 ± 0.09 (1 SD) for TG-43 and Acuros BV algorithms, respectively. Due to the air pocket, mucosa received a reduced dose by an average of 22.72% and an average of 23.29% for TG-43 and Acuros BV, respectively. The maximum displacement of mucosa and the ratio of doses were negatively correlated for both. In the Optimized Acuros BV plan, total dwell time increased by 1.8% but no considerable change in the dose ratios.

Conclusion

The calculated dose of mucous membrane forced out of the cylinder surface by air pockets by the Acuros BV algorithm was nonsignificantly different from TG-43. Therefore, even in the presence of air pockets, the TG-43 algorithm for calculating the VCB dose is appropriate.

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.

Dose perturbation effects by metal hip prosthesis in gynaecological 192Ir HDR brachytherapy

The present study was conducted to elucidate the effects of hip prostheses in ¹⁹²Ir HDR brachytherapy and determine dose uncertainties introduced by the treatment planning. A gynaecological phantom irradiated using Nucletron ¹⁹²Ir microSelectron HDR source was modeled using MCNP5 code. Three hip materials considered in this study were water, bone, and metal prosthesis. According to the obtained results, a dose perturbation was observed within the medium with a higher atomic number, which reduced the dose to the nearby region.

Accuracy of Acuros TM 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 AcurosTM 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 (192 Ir) source. Accuracy assessment has been achieved by implementing a modelled GamaMed Plus192 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.

Evaluating dosimetric accuracy of the 6 MV calibration on EBT3 film in the use of Ir-192 high dose rate brachytherapy

Purpose

To evaluate the dosimetric accuracy of EBT3 film calibrated with a 6 MV beam for high dose rate brachytherapy and propose a novel method for direct film calibration with an Ir‐192 source.

Methods

The 6 MV calibration was performed in water on a linear accelerator (linac). The Ir‐192 calibration was accomplished by irradiating the film wrapped around a cylinder applicator with an Ir‐192 source. All films were scanned 1‐day post‐irradiation to acquire calibration curves for all three (red, blue, and green) channels. The Ir‐192 calibration films were also used for single‐dose comparison. Moreover, an independent test film under a H.A.M. applicator was irradiated and the 2D dose distribution was obtained separately for each calibration using the red channel data. Gamma analysis and point‐by‐point profile comparison were performed to evaluate the performance of both calibrations. The uncertainty budget for each calibration system was analyzed.

Results

The red channel had the best performance for both calibration systems in the single‐dose comparison. We found a significant 4.89% difference from the reference for doses <250 cGy using the 6 MV calibration, while the difference was only 0.87% for doses >600 cGy. Gamma analysis of the 2D dose distribution showed the Ir‐192 calibration had a higher passing rate of 91.9% for the 1 mm/2% criterion, compared to 83.5% for the 6 MV calibration. Most failing points were in the low‐dose region (<200 cGy). The point‐by‐point profile comparison reported a discrepancy of 2%–3.6% between the Ir‐192 and 6 MV calibrations in this low‐dose region. The linac‐ and Ir‐192‐based dosimetry systems had an uncertainty of 4.1% (k = 2) and 5.66% (k = 2), respectively.

Conclusions

Direct calibration of EBT3 films with an Ir‐192 source is feasible and reliable, while the dosimetric accuracy of 6 MV calibration depends on the dose range. The Ir‐192 calibration should be used when the measurement dose range is below 250 cGy.

Personalized brachytherapy dose reconstruction using deep learning

BACKGROUND AND PURPOSE

Accurate calculation of the absorbed dose delivered to the tumor and normal tissues improves treatment gain factor, which is the major advantage of brachytherapy over external radiation therapy. To address the simplifications of TG-43 assumptions that ignore the dosimetric impact of medium heterogeneities, we proposed a deep learning (DL)-based approach, which improves the accuracy while requiring a reasonable computation time.

MATERIALS AND METHODS

We developed a Monte Carlo (MC)-based personalized brachytherapy dosimetry simulator (PBrDoseSim), deployed to generate patient-specific dose distributions. A deep neural network (DNN) was trained to predict personalized dose distributions derived from MC simulations, serving as ground truth. The paired channel input used for the training is composed of dose distribution kernel in water medium along with the full-volumetric density maps obtained from CT images reflecting medium heterogeneity.

RESULTS

The predicted single-dwell dose kernels were in good agreement with MC-based kernels serving as reference, achieving a mean relative absolute error (MRAE) and mean absolute error (MAE) of 1.16 ± 0.42% and 4.2 ± 2.7 × 10^-4 (Gy.sec^-1/voxel), respectively. The MRAE of the dose volume histograms (DVHs) between the DNN and MC calculations in the clinical target volume were 1.8 ± 0.86%, 0.56 ± 0.56%, and 1.48 ± 0.72% for D90, V150, and V100, respectively. For bladder, sigmoid, and rectum, the MRAE of D5cc between the DNN and MC calculations were 2.7 ± 1.7%, 1.9 ± 1.3%, and 2.1 ± 1.7%, respectively.

CONCLUSION

The proposed DNN-based personalized brachytherapy dosimetry approach exhibited comparable performance to the MC method while overcoming the computational burden of MC calculations and oversimplifications of TG-43.

Characterization of Dosimetric Differences in Strut-Adjusted Volume Implant Treatment Plans Calculated With TG-43 Formalism and a Model-Based Dose Calculation Algorithm

PURPOSE

To comprehensively characterize dosimetric differences between calculations with a commercial model-based dose calculation algorithm (MBDCA) and the TG-43 formalism in application to accelerated partial breast irradiation (APBI) with the strut-adjusted volume implant (SAVI) applicator.

METHODS

Dose for 100 patients treated with the SAVI applicator was recalculated with an MBDCA for comparison to dose calculated via TG-43. For every pair of dose calculations, dose-volume histogram (DVH) metrics including V90%, V95%, V100%, V150%, and V200% for the PTV_EVAL were compared. Features were defined for each case including (1) applicator size, (2) ratio between PTV_EVAL contour and 1-cm rind surrounding SAVI applicator, (3) ratio between dwell time in central catheter and total dwell time, and (4) mean computed tomography (CT) number within the lumpectomy cavity. Wilcoxon rank sum tests were performed to test whether treatment plans could be stratified according to feature values into groups with statistically significant dosimetry differences between MBDCA and TG-43.

RESULTS

For all DVH metrics, differences between TG-43 and MBDCA calculations were statistically significant (P < .05). Minimum (maximum) relative percent differences between the MBDCA and TG-43 for V90%, V95%, and V100% were -2.1% (0.1%), -3.1% (-0.1%), and -5.0% (-0.5%), respectively. The median relative percent difference in mean PTV_EVAL dose between the MBDCA and TG-43 was -3.9%, with minimum (maximum) difference of -6.5% (-1.8%). For V90%, V95%, and V100%, plan quality worsened beyond defined thresholds in 26, 23, and 31 cases with no instances of coverage improvement. Features 1, 2, and 4 were shown to be able to stratify treatment plans into groups with statistically significant differences in dosimetry metrics between MBDCA and TG-43.

CONCLUSIONS

Investigated dose metrics for SAVI treatments were found to be systematically lower with MBDCA calculation in comparison to TG-43. Plans could be stratified according to several features by the magnitude of dosimetric differences between these calculations.

Validation of the TOPAS Monte Carlo toolkit for HDR brachytherapy simulations

PURPOSE

The goal of this work is to validate the user-friendly Geant4-based Monte Carlo toolkit TOol for PArticle Simulation (TOPAS) for brachytherapy applications.

METHODS AND MATERIALS

Brachytherapy simulations performed with TOPAS were systematically compared with published TG-186 reference data. The photon emission energy spectrum, the air-kerma strength, and the dose-rate constant of the model-based dose calculation algorithm (MBDCA)-WG generic Ir-192 source were extracted. For dose calculations, a track-length estimator was implemented. The four Joint AAPM/ESTRO/ABG MBDCA-WG test cases were evaluated through histograms of the local and global dose difference volumes. A prostate, a palliative lung, and a breast case were simulated. For each case, the dose ratio map, the histogram of the global dose difference volume, and cumulative dose-volume histograms were calculated.

RESULTS

The air-kerma strength was (9.772 ± 0.001) × 10-8 U Bq-1 (within 0.3% of the reference value). The dose-rate constant was 1.1107 ± 0.0005 cGy h-1 U-1 (within 0.01% of the reference value). For all cases, at least 96.9% of voxels had a local dose difference within [-1%, 1%] and at least 99.9% of voxels had a global dose difference within [-0.1%, 0.1%]. The implemented track-length estimator scorer was more efficient than the default analog dose scorer by a factor of 237. For all clinical cases, at least 97.5% of voxels had a global dose difference within [-1%, 1%]. Dose-volume histograms were consistent with the reference data.

CONCLUSIONS

TOPAS was validated for high-dose-rate brachytherapy simulations following the TG-186 recommended approach for MBDCAs. Built on top of Geant4, TOPAS provides broad access to a state-of-the-art Monte Carlo code for brachytherapy simulations.

Impact of GBBS algorithm on post-mastectomy scar boost irradiation of breast using catheter flap

Purpose

Post-mastectomy radiation therapy significantly reduces locoregional recurrence rates, which can be achieved with external beam radiotherapy delivered to chest wall, followed by scar irradiation either by electron or high-dose-rate (HDR) mould brachytherapy. The present study evaluates dosimetric advantage of Acuros^® BV, a TG-186 MBDCA, over TG-43 formalism using ¹⁹²Ir source for HDR brachytherapy in chest wall scar boost using catheter flap.

Material and methods

A total of 25 patients, free of cardiac and pulmonary co-morbidities, who met the inclusion criteria were involved in the study. Catheter flap made of silicon with 20 channels was used to deliver a total dose of 7.5 Gy/3 fx by HDR surface mould brachytherapy to delineated scar volume. Plan was optimized with iterative method to obtain desired results with TG-43 formalism, followed by Acuros^® BV (GBBS algorithm) without altering dwell positions or time. The two algorithm plans were analyzed qualitatively and quantitatively with dose-volume histograms.

Results

The mean D_98% CTV-HDR_evl coverage decreased by 1.16% compared to TG-43, and near-maximum dose decreased by 8.18% (p = 0.000), mean D_max dose to CTV-HDR_evl, and mean D_mean dose was lesser by 6.25% (p = 0.000) and 10.82% (p = 0.000), respectively, compared to TG-43. Heart D_2% showed significant results, whereas D_median (cGy) revealed very significant difference. A 5 mm thick skin contour showed statistically significant results (p = 0.000) for V_150% and V_200%.

Conclusions

The presented data showed how Acuros^® BV, algorithm-based calculation in scar boost irradiation of breast, accounting for a mass density of the medium and scatter condition, considered actual dose prediction in a medium.

Clinically-implementable template plans for multidwell treatments using Leipzig-style applicators in 192Ir surface brachytherapy

PURPOSE

Multiple dwell positions ("multidwell") within a Leipzig-style applicator can be used to increase dose uniformity and treatment area. Model-based dose calculation algorithms (MBDCAs) are necessary for accurate calculations involving these applicators because of their nonwater equivalency and complex geometry. The purpose of this work was to create template plans from MBDCA calculations and present their dwell times and positions for users of these applicators without access to MBDCAs.

METHODS AND MATERIALS

The Leipzig-style solid applicator model within our treatment planning system was used to design template plans. Five template plans, normalized to 0.3 cm depth within a water phantom, were optimized using the treatment planning system MBDCA. Each template plan contained unique dwell positions, times, and active lengths (0.5-1.5 cm). A single-dwell distribution was optimized for comparison. The stem of this applicator stops within the shell; therefore, one template plan contained an intrafraction rotation to determine the largest dose distribution achievable. Effects of imperfect applicator rotation were quantified by inserting rotational offsets and comparing the V_100%, D_95%, and minimum dose coverage for planning target volumes created from 80%/90% isodose lines.

RESULTS

The 90% (80%) isodose line dimensions at 0.3 cm depth for single-dwell increased from 0.94 × 0.97 (1.53 × 1.57) cm² to 2.09 × 1.24 (2.75 × 1.88) cm² in the largest template plan. Manually inserted angular offsets up to ±10° for the template plan requiring rotation preserved V_100%, D_95%, and minimum dose within 2.0%, 1.9%, and 8.0%, respectively.

CONCLUSION

A set of template plans was created to provide accessibility to the multidwell methodology, even for users without access to MBDCAs. Each template plan should be reviewed before clinical implementation.

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.

Feasibility of using multiple-dwell positions in 192Ir Leipzig-style brachytherapy surface applicators to expand target coverage and clinical application

PURPOSE

Leipzig-style applicators for surface brachytherapy are traditionally used with a single-source dwell position. This study explores the feasibility of using multiple-source dwell positions ("multidwell") to improve the dose coverage and applicability of Leipzig-style applicators.

METHODS AND MATERIALS

A virtual model of the Leipzig-style applicator was commissioned for a model-based dose calculation algorithm (MBDCA) and compared against American Association of Physicists in Medicine working group 186 benchmarking data sets and ionization chamber point measurements. An absolute dosimetry technique based on radiochromic film was used to validate both single-dwell and multidwell plans.

RESULTS

Dose distributions generated from the MBDCA-based virtual model were consistent with working group data sets, ion chamber measurements, and radiochromic film analysis. In one multidwell configuration, at 3 mm prescription depth, the 80% isodose width was increased to 25 mm, compared with 15 mm in the same dimension for a single-dwell delivery. In the same multidwell configuration, the flatness, measured as >98% isodose line, was more than doubled to 8 mm, compared with 3 mm in the same dimension. For multidwell plans, 2-D planar agreement between radiochromic film and MBDCA exceeded 93% in gamma analysis (3%/1 mm criteria). Submillimeter positional agreement was found, with a total dosimetric uncertainty of 4.5% estimated for the entire system.

CONCLUSIONS

Leipzig-style surface applicators with multiple-source dwell positions have been benchmarked against radiochromic film dosimetry. Results show that the clinically viable coverage area can be increased significantly.

Dosimetric comparison of surface mould HDR brachytherapy with VMAT

INTRODUCTION

The aim of this study was to investigate the dosimetric differences between surface mould high-dose-rate (HDR) brachytherapy and external beam volumetric-modulated arc therapy (VMAT) for two treatment sites.

METHODS

Previously treated HDR brachytherapy surface mould scalp (n = 4) and lower leg (n = 3) treatments were retrospectively analysed. The VMAT plans were optimised using an additional 3-mm setup margin on the clinical target volume (CTV) of the previously treated HDR plans. The HDR plans were calculated and normalised using the TG-43 formalism and recalculated with Acuros BV (AC).

RESULTS

On average, the mean brain and normal tissue doses were reduced by 44.8% and 27.4% for scalp and lower leg VMAT cases, respectively, when compared to AC calculated HDR plans. For VMAT plans, the average dose to a 1-mm thick skin structure deep to the target volume was not any lower than that in AC HDR plans. On average, the CTV coverage was 13.8% and 9.6% lower for scalp cases with AC dose calculation than with TG-43 and 8.3% and 5.3% lower for lower leg cases if 0- or 1-cm backscatter material was applied above the catheters, respectively.

CONCLUSIONS

VMAT is a feasible treatment option in the case of extensive skin malignancies of the scalp and lower leg. Uncertainties related to delivered dose with HDR brachytherapy when using the TG-43 dose calculation model or possible air gaps between the mould and skin favour the use of VMAT. The potential soft tissue deformation needs to be considered if VMAT is used.

The validity of Acuros BV and TG-43 for high-dose-rate brachytherapy superficial mold treatments

PURPOSE

The purpose of this work is to validate the Acuros BV dose calculation algorithm for high-dose-rate (HDR) brachytherapy superficial mold treatments in the absence of full scatter conditions and compare this with TG-43 dose calculations. We also investigate the impact of additional back scatter material (bolus) applied above surface molds to the dose distributions under the mold.

METHODS AND MATERIALS

The absorbed dose at various depths was compared for simulations performed using either TG-43 or Acuros BV dose calculations. Parameter variations included treatment area, thickness of the bolus, and surface shape (flat or spherical). Film measurements were carried out in a flat phantom.

RESULTS

Acuros BV calculations and film measurements agreed within 1.5% but were up to 15% lower than TG-43 dose calculations when no bolus was applied above the treatment catheters. The difference in dose at the prescription depth (1 cm below the central catheter) increased with increasing treatment area: 3.3% difference for a 3 × 3.5 cm² source loading area, 7.4% for 8 × 9 cm², and 13.4% for 18 × 19 cm². The dose overestimation of the TG-43 model decreased when bolus was added above the treatment catheters.

CONCLUSIONS

The TG-43 dosimetry formalism cannot model surface mold treatments in the absence of full scatter conditions within 5% for loading areas larger than approximately 5 × 5 cm². The TG-43 model results in an overestimation of the delivered dose, which increases with treatment area. This confirms the need for model-based dose calculation algorithms as discussed in TG-186.

Dosimetric and radiobiological comparison of TG-43 and Monte Carlo calculations in 192Ir breast brachytherapy applications

PURPOSE

To investigate the clinical significance of introducing model based dose calculation algorithms (MBDCAs) as an alternative to TG-43 in ¹⁹²Ir interstitial breast brachytherapy.

MATERIALS AND METHODS

A 57 patient cohort was used in a retrospective comparison between TG-43 based dosimetry data exported from a treatment planning system and Monte Carlo (MC) dosimetry performed using MCNP v. 6.1 with plan and anatomy information in DICOM-RT format. Comparison was performed for the target, ipsilateral lung, heart, skin, breast and ribs, using dose distributions, dose-volume histograms (DVH) and plan quality indices clinically used for plan evaluation, as well as radiobiological parameters.

RESULTS

TG-43 overestimation of target DVH parameters is statistically significant but small (less than 2% for the target coverage indices and 4% for homogeneity indices, on average). Significant dose differences (>5%) were observed close to the skin and at relatively large distances from the implant leading to a TG-43 dose overestimation for the organs at risk. These differences correspond to low dose regions (<50% of the prescribed dose), being less than 2% of the prescribed dose. Detected dosimetric differences did not induce clinically significant differences in calculated tumor control probabilities (mean absolute difference <0.2%) and normal tissue complication probabilities.

CONCLUSION

While TG-43 shows a statistically significant overestimation of most indices used for plan evaluation, differences are small and therefore not clinically significant. Improved MBDCA dosimetry could be important for re-irradiation, technique inter-comparison and/or the assessment of secondary cancer induction risk, where accurate dosimetry in the whole patient anatomy is of the essence.

Comparison of image-based three-dimensional treatment planning using AcurosTM BV and AAPM TG-43 algorithm for intracavitary brachytherapy of carcinoma cervix

Aim

To compare the image-based three-dimensional treatment planning using AcurosTM BV and AAPM TG-43 algorithm for intracavitary brachytherapy of carcinoma cervix.

Materials and methods

Twenty-seven patients with cancer cervix, stage IIB or IIIB with vaginal involvement limited to the upper third of the vagina was included into the study. Intracavitary treatments with the patient in this study done with computed tomography and magnetic resonance imaging compatible ring applicator. Groupe European de Curietherapie and European Society for Therapeutic Radiology and Oncology recommended doses to target volumes and organs at risk compared using dose volume histogram.

Results

The mean value of Point ‘A’ dose was compared between AcurosTM BV and TG-43, which indicates 0·13% difference. The differences in the mean dose to gross tumour volume for various volumes are V100% 0·28%, V150% 1·22% and V200% 1·03%; all volumes showed small difference but statistical significant (p<0·05). The mean dose of high-risk clinical target volume (HRCTV) D90 using AcurosTM BV was 8·47 Gy, which was 1·63% less compared with TG-43. The mean point A dose using AcurosTM BV is 1·04 times the dose to D90 of mean HRCTV. The same difference was observed in comparison with TG43. D2cc and D0·1cc of the bladder, rectum and sigmoid showed a statistically significant difference (p<0·05) in comparison with TG-43.

Conclusion

The differences in dosimetric parameters between the AcurosTM BV and TG-43 proved to be statistically significant. The difference is very small, and they are clinically insignificant.

Impact of heterogeneity-corrected dose calculation using a grid-based Boltzmann solver on breast and cervix cancer brachytherapy

Purpose

To analyze the impact of heterogeneity-corrected dose calculation on dosimetric quality parameters in gynecological and breast brachytherapy using Acuros, a grid-based Boltzmann equation solver (GBBS), and to evaluate the shielding effects of different cervix brachytherapy applicators.

Material and methods

Calculations with TG-43 and Acuros were based on computed tomography (CT) retrospectively, for 10 cases of accelerated partial breast irradiation and 9 cervix cancer cases treated with tandem-ring applicators. Phantom CT-scans of different applicators (plastic and titanium) were acquired. For breast cases the V_20Gyαβ3 to lung, the D_0.1cm³, D_1cm³, D_2cm³ to rib, the D_0.1cm³, D_1cm³, D_10cm³ to skin, and D_max for all structures were reported. For cervix cases, the D_0.1cm³, D_2cm³ to bladder, rectum and sigmoid, and the D₅₀, D₉₀, D₉₈, V₁₀₀ for the CTV_HR were reported. For the phantom study, surrogates for target and organ at risk were created for a similar dose volume histogram (DVH) analysis. Absorbed dose and equivalent dose to 2 Gy fractionation (EQD2) were used for comparison.

Results

Calculations with TG-43 overestimated the dose for all dosimetric indices investigated. For breast, a decrease of ~8% was found for D_10cm³ to the skin and 5% for D_2cm³ to rib, resulting in a difference ~ –1.5 Gy EQD2 for overall treatment. Smaller effects were found for cervix cases with the plastic applicator, with up to –2% (–0.2 Gy EQD2) per fraction for organs at risk and –0.5% (–0.3 Gy EQD2) per fraction for CTV_HR. The shielding effect of the titanium applicator resulted in a decrease of 2% for D_2cm³ to the organ at risk versus 0.7% for plastic.

Conclusions

Lower doses were reported when calculating with Acuros compared to TG-43. Differences in dose parameters were larger in breast cases. A lower impact on clinical dose parameters was found for the cervix cases. Applicator material causes systematic shielding effects that can be taken into account.

A generic high-dose rate (192)Ir brachytherapy source for evaluation of model-based dose calculations beyond the TG-43 formalism

PURPOSE

In order to facilitate a smooth transition for brachytherapy dose calculations from the American Association of Physicists in Medicine (AAPM) Task Group No. 43 (TG-43) formalism to model-based dose calculation algorithms (MBDCAs), treatment planning systems (TPSs) using a MBDCA require a set of well-defined test case plans characterized by Monte Carlo (MC) methods. This also permits direct dose comparison to TG-43 reference data. Such test case plans should be made available for use in the software commissioning process performed by clinical end users. To this end, a hypothetical, generic high-dose rate (HDR) (192)Ir source and a virtual water phantom were designed, which can be imported into a TPS.

METHODS

A hypothetical, generic HDR (192)Ir source was designed based on commercially available sources as well as a virtual, cubic water phantom that can be imported into any TPS in DICOM format. The dose distribution of the generic (192)Ir source when placed at the center of the cubic phantom, and away from the center under altered scatter conditions, was evaluated using two commercial MBDCAs [Oncentra(®) Brachy with advanced collapsed-cone engine (ACE) and BrachyVision ACUROS™ ]. Dose comparisons were performed using state-of-the-art MC codes for radiation transport, including ALGEBRA, BrachyDose, GEANT4, MCNP5, MCNP6, and PENELOPE2008. The methodologies adhered to recommendations in the AAPM TG-229 report on high-energy brachytherapy source dosimetry. TG-43 dosimetry parameters, an along-away dose-rate table, and primary and scatter separated (PSS) data were obtained. The virtual water phantom of (201)(3) voxels (1 mm sides) was used to evaluate the calculated dose distributions. Two test case plans involving a single position of the generic HDR (192)Ir source in this phantom were prepared: (i) source centered in the phantom and (ii) source displaced 7 cm laterally from the center. Datasets were independently produced by different investigators. MC results were then compared against dose calculated using TG-43 and MBDCA methods.

RESULTS

TG-43 and PSS datasets were generated for the generic source, the PSS data for use with the ace algorithm. The dose-rate constant values obtained from seven MC simulations, performed independently using different codes, were in excellent agreement, yielding an average of 1.1109 ± 0.0004 cGy/(h U) (k = 1, Type A uncertainty). MC calculated dose-rate distributions for the two plans were also found to be in excellent agreement, with differences within type A uncertainties. Differences between commercial MBDCA and MC results were test, position, and calculation parameter dependent. On average, however, these differences were within 1% for ACUROS and 2% for ace at clinically relevant distances.

CONCLUSIONS

A hypothetical, generic HDR (192)Ir source was designed and implemented in two commercially available TPSs employing different MBDCAs. Reference dose distributions for this source were benchmarked and used for the evaluation of MBDCA calculations employing a virtual, cubic water phantom in the form of a CT DICOM image series. The implementation of a generic source of identical design in all TPSs using MBDCAs is an important step toward supporting univocal commissioning procedures and direct comparisons between TPSs.

A user-oriented procedure for the commissioning and quality assurance testing of treatment planning system dosimetry in high-dose-rate brachytherapy

PURPOSE

To develop a user-oriented procedure for testing treatment planning system (TPS) dosimetry in high-dose-rate brachytherapy, with particular focus to TPSs using model-based dose calculation algorithms (MBDCAs).

METHODS AND MATERIALS

Identical plans were prepared for three computational models using two commercially available systems and the same (192)Ir source. Reference dose distributions were obtained for each plan using the MCNP v.6.1 Monte Carlo (MC) simulation code with input files prepared via automatic parsing of plan information using a custom software tool. The same tool was used for the comparison of reference dose distributions with corresponding MBDCA exports.

RESULTS

The single source test case yielded differences due to the MBDCA spatial discretization settings. These affect points at relatively increased distance from the source, and they are abated in test cases with multiple source dwells. Differences beyond MC Type A uncertainty were also observed very close to the source(s), close to the test geometry boundaries, and within heterogeneities. Both MBDCAs studied were found equivalent to MC within 5 cm from the target volume for a clinical breast brachytherapy test case. These are in agreement with previous findings of MBDCA benchmarking in the literature.

CONCLUSIONS

The data and the tools presented in this work, that are freely available via the web, can serve as a benchmark for advanced clinical users developing their own tests, a complete commissioning procedure for new adopters of currently available TPSs using MBDCAs, a quality assurance testing tool for future updates of already installed TPSs, or as an admission prerequisite in multicentric clinical trials.

Validation of the Oncentra Brachy Advanced Collapsed cone Engine for a commercial (192)Ir source using heterogeneous geometries

PURPOSE

To validate the Advanced Collapsed cone Engine (ACE) dose calculation engine of Oncentra Brachy (OcB) treatment planning system using an (192)Ir source.

METHODS AND MATERIALS

Two levels of validation were performed, conformant to the model-based dose calculation algorithm commissioning guidelines of American Association of Physicists in Medicine TG-186 report. Level 1 uses all-water phantoms, and the validation is against TG-43 methodology. Level 2 uses real-patient cases, and the validation is against Monte Carlo (MC) simulations. For each case, the ACE and TG-43 calculations were performed in the OcB treatment planning system. ALGEBRA MC system was used to perform MC simulations.

RESULTS

In Level 1, the ray effect depends on both accuracy mode and the number of dwell positions. The volume fraction with dose error ≥2% quickly reduces from 23% (13%) for a single dwell to 3% (2%) for eight dwell positions in the standard (high) accuracy mode. In Level 2, the 10% and higher isodose lines were observed overlapping between ACE (both standard and high-resolution modes) and MC. Major clinical indices (V100, V150, V200, D90, D50, and D2cc) were investigated and validated by MC. For example, among the Level 2 cases, the maximum deviation in V100 of ACE from MC is 2.75% but up to ~10% for TG-43. Similarly, the maximum deviation in D90 is 0.14 Gy between ACE and MC but up to 0.24 Gy for TG-43.

CONCLUSION

ACE demonstrated good agreement with MC in most clinically relevant regions in the cases tested. Departure from MC is significant for specific situations but limited to low-dose (<10% isodose) regions.

Dosimetric comparison of Acuros™ BV with AAPM TG43 dose calculation formalism in breast interstitial high-dose-rate brachytherapy with the use of metal catheters

Purpose

Radiotherapy for breast cancer includes different techniques and methods. The purpose of this study is to compare dosimetric calculations using TG-43 dose formalism and Varian Acuros™ BV (GBBS) dose calculation algorithm for interstitial implant of breast using metal catheters in high-dose-rate (HDR) brachytherapy, using ¹⁹²Ir.

Material and methods

Twenty patients who were considered for breast conservative surgery (BCS), underwent lumpectomy and axillary dissection. These patients received perioperative interstitial HDR brachytherapy as upfront boost using rigid metal implants. Whole breast irradiation was delivered TG-43 after a gap of two weeks. Standard brachytherapy dose calculation was done by dosimetry. This does not take into account tissue heterogeneity, attenuation and scatter in the metal applicator, and effects of patient boundary. Acuros™ BV is a Grid Based Boltzmann Solver code (GBBS), which takes into consideration all the above, was used to compute dosimetry and the two systems were compared.

Results

Comparison of GBBS and TG-43 formalism on interstitial metal catheters shows difference in dose prescribed to CTV and other OARs. While the estimated dose to CTV was only marginally different with the two systems, there is a significant difference in estimated doses of starting from 4 to 53% in the mean value of all parameters analyzed.

Conclusions

TG-43 algorithm seems to significantly overestimate the dose to various volumes of interest; GBBS based dose calculation algorithm has impact on CTV, heart, ipsilateral lung, heart, contralateral breast, skin, and ribs of the ipsilateral breast side; the prescription changes occurred due to effect of metal catheters, inhomogeneities, and scatter conditions.

BrachyGuide: a brachytherapy-dedicated DICOM RT viewer and interface to Monte Carlo simulation software

This work presents BrachyGuide, a brachytherapy‐dedicated software tool for the automatic preparation of input files for Monte Carlo simulation from treatment plans exported in DICOM RT format, and results of calculations performed for its benchmarking. Three plans were prepared using two computational models, the image series of a water sphere and a multicatheter breast brachytherapy patient, for each of two commercially available treatment planning systems: BrachyVision and Oncentra Brachy. One plan involved a single source dwell position of an 192Ir HDR source (VS2000 or mHDR‐v2) at the center of the water sphere using the TG43 algorithm, and the other two corresponded to the TG43 and advanced dose calculation algorithm for the multicatheter breast brachytherapy patient. Monte Carlo input files were prepared using BrachyGuide and simulations were performed with MCNP v.6.1. For the TG43 patient plans, the Monte Carlo computational model was manually edited in the prepared input files to resemble TG43 dosimetry assumptions. Hence all DICOM RT dose exports were equivalent to corresponding simulation results and their comparison was used for benchmarking the use of BrachyGuide. Monte Carlo simulation results and corresponding DICOM RT dose exports agree within type A uncertainties in the majority of points in the computational models. Treatment planning system, algorithm, and source specific differences greater than type A uncertainties were also observed, but these were explained by treatment planning system‐related issues and other sources of type B uncertainty. These differences have to be taken into account in commissioning procedures of brachytherapy dosimetry algorithms. BrachyGuide is accurate and effective for use in the preparation of commissioning tests for new brachytherapy dosimetry algorithms as a user‐oriented commissioning tool and the expedition of retrospective patient cohort studies of dosimetry planning.

PACS numbers: 87.53.Bn, 87.53.Jw, 87.55.D‐, 87.55.Qr, 87.55.km, 87.55.K‐

Current state of the art brachytherapy treatment planning dosimetry algorithms

Following literature contributions delineating the deficiencies introduced by the approximations of conventional brachytherapy dosimetry, different model-based dosimetry algorithms have been incorporated into commercial systems for (192)Ir brachytherapy treatment planning. The calculation settings of these algorithms are pre-configured according to criteria established by their developers for optimizing computation speed vs accuracy. Their clinical use is hence straightforward. A basic understanding of these algorithms and their limitations is essential, however, for commissioning; detecting differences from conventional algorithms; explaining their origin; assessing their impact; and maintaining global uniformity of clinical practice.

Head and neck (192)Ir HDR-brachytherapy dosimetry using a grid-based Boltzmann solver

Purpose

To compare dosimetry for head and neck cancer patients, calculated with TG-43 formalism and a commercially available grid-based Boltzmann solver.

Material and methods

This study included 3D-dosimetry of 49 consecutive brachytherapy head and neck cancer patients, computed by a grid-based Boltzmann solver that takes into account tissue inhomogeneities as well as TG-43 formalism. 3D-treatment planning was carried out by using computed tomography.

Results

Dosimetric indices D₉₀ and V₁₀₀ for target volume were about 3% lower (median value) for the grid-based Boltzmann solver relative to TG-43-based computation (p < 0.01). The V₁₅₀ dose parameter showed 1.6% increase from grid-based Boltzmann solver to TG-43 (p < 0.01).

Conclusions

Dose differences between results of a grid-based Boltzmann solver and TG-43 formalism for high-dose-rate head and neck brachytherapy patients to the target volume were found. Distinctions in D₉₀ of CTV were low (2.63 Gy for grid-based Boltzmann solver vs. 2.71 Gy TG-43 in mean). In our clinical practice, prescription doses remain unchanged for high-dose-rate head and neck brachytherapy for the time being.

Dosimetric accuracy of a deterministic radiation transport based (192)Ir brachytherapy treatment planning system. Part III. Comparison to Monte Carlo simulation in voxelized anatomical computational models

PURPOSE

To compare TG43-based and Acuros deterministic radiation transport-based calculations of the BrachyVision treatment planning system (TPS) with corresponding Monte Carlo (MC) simulation results in heterogeneous patient geometries, in order to validate Acuros and quantify the accuracy improvement it marks relative to TG43.

METHODS

Dosimetric comparisons in the form of isodose lines, percentage dose difference maps, and dose volume histogram results were performed for two voxelized mathematical models resembling an esophageal and a breast brachytherapy patient, as well as an actual breast brachytherapy patient model. The mathematical models were converted to digital imaging and communications in medicine (DICOM) image series for input to the TPS. The MCNP5 v.1.40 general-purpose simulation code input files for each model were prepared using information derived from the corresponding DICOM RT exports from the TPS.

RESULTS

Comparisons of MC and TG43 results in all models showed significant differences, as reported previously in the literature and expected from the inability of the TG43 based algorithm to account for heterogeneities and model specific scatter conditions. A close agreement was observed between MC and Acuros results in all models except for a limited number of points that lay in the penumbra of perfectly shaped structures in the esophageal model, or at distances very close to the catheters in all models.

CONCLUSIONS

Acuros marks a significant dosimetry improvement relative to TG43. The assessment of the clinical significance of this accuracy improvement requires further work. Mathematical patient equivalent models and models prepared from actual patient CT series are useful complementary tools in the methodology outlined in this series of works for the benchmarking of any advanced dose calculation algorithm beyond TG43.

Commissioning of a grid-based Boltzmann solver for cervical cancer brachytherapy treatment planning with shielded colpostats

PURPOSE

We sought to commission a gynecologic shielded colpostat analytic model provided from a treatment planning system (TPS) library. We have reported retrospectively the dosimetric impact of this applicator model in a cohort of patients.

METHODS AND MATERIALS

A commercial TPS with a grid-based Boltzmann solver (GBBS) was commissioned for (192)Ir high-dose-rate (HDR) brachytherapy for cervical cancer with stainless steel-shielded colpostats. Verification of the colpostat analytic model was verified using a radiograph and vendor schematics. MCNPX v2.6 Monte Carlo simulations were performed to compare dose distributions around the applicator in water with the TPS GBBS dose predictions. Retrospectively, the dosimetric impact was assessed over 24 cervical cancer patients' HDR plans.

RESULTS

Applicator (TPS ID #AL13122005) shield dimensions were within 0.4 mm of the independent shield dimensions verification. GBBS profiles in planes bisecting the cap around the applicator agreed with Monte Carlo simulations within 2% at most locations; differing screw representations resulted in differences of up to 9%. For the retrospective study, the GBBS doses differed from TG-43 as follows (mean value ± standard deviation [min, max]): International Commission on Radiation units [ICRU]rectum (-8.4 ± 2.5% [-14.1, -4.1%]), ICRUbladder (-7.2 ± 3.6% [-15.7, -2.1%]), D2cc-rectum (-6.2 ± 2.6% [-11.9, -0.8%]), D2cc-sigmoid (-5.6 ± 2.6% [-9.3, -2.0%]), and D2cc-bladder (-3.4 ± 1.9% [-7.2, -1.1%]).

CONCLUSIONS

As brachytherapy TPSs implement advanced model-based dose calculations, the analytic applicator models stored in TPSs should be independently validated before clinical use. For this cohort, clinically meaningful differences (>5%) from TG-43 were observed. Accurate dosimetric modeling of shielded applicators may help to refine organ toxicity studies.

Sub-second high dose rate brachytherapy Monte Carlo dose calculations with bGPUMCD

PURPOSE

To establish the accuracy and speed of bGPUMCD, a GPU-oriented Monte Carlo code used for high dose rate brachytherapy dose calculations. The first objective is to evaluate the time required for dose calculation when full Monte Carlo generated dose distribution kernels are used for plan optimization. The second objective is to assess the accuracy and speed when recalculating pre-optimized plans, consisting of many dwell positions.

METHODS

bGPUMCD is tested with three clinical treatment plans : one prostate case, one breast case, and one rectum case with a shielded applicator. Reference distributions, generated with GEANT4, are used as a basis of comparison. Calculations of full dose distributions of pre-optimized treatment plans as well as single dwell dosimetry are performed. Single source dosimetry, based on TG-43 parameters reproduction, is also presented for the microSelectron V2 (Nucletron, Veenendaal, The Netherlands).

RESULTS

In timing experiments, the computation of single dwell position dose kernels takes between 0.25 and 0.5 s. bGPUMCD can compute full dose distributions of previously optimized plans in ∼2 s. bGPUMCD is capable of computing pre-optimized brachytherapy plans within 1% for the prostate case and 2% for the breast and shielded applicator cases, when comparing the dosimetric parameters D90 and V100 of the reference (GEANT4) and bGPUMCD distributions. For all voxels within the target, an absolute average difference of approximately 1% is found for the prostate case, less than 2% for the breast case and less than 2% for the rectum case with shielded applicator. Larger point differences (>5%) are found within bony regions in the prostate case, where bGPUMCD underdoses compared to GEANT4. Single source dosimetry results are mostly within 2% for the radial function and within 1%-4% for the anisotropic function.

CONCLUSIONS

bGPUMCD has the potential to allow for fast MC dose calculation in a clinical setting for all phases of HDR treatment planning, from dose kernel calculations for plan optimization to plan recalculation.

Physics-aspects of dose accuracy in high dose rate (HDR) brachytherapy: source dosimetry, treatment planning, equipment performance and in vivo verification techniques

This study provides a review of recent publications on the physics-aspects of dosimetric accuracy in high dose rate (HDR) brachytherapy. The discussion of accuracy is primarily concerned with uncertainties, but methods to improve dose conformation to the prescribed intended dose distribution are also noted. The main aim of the paper is to review current practical techniques and methods employed for HDR brachytherapy dosimetry. This includes work on the determination of dose rate fields around brachytherapy sources, the capability of treatment planning systems, the performance of treatment units and methods to verify dose delivery. This work highlights the determinants of accuracy in HDR dosimetry and treatment delivery and presents a selection of papers, focusing on articles from the last five years, to reflect active areas of research and development. Apart from Monte Carlo modelling of source dosimetry, there is no clear consensus on the optimum techniques to be used to assure dosimetric accuracy through all the processes involved in HDR brachytherapy treatment. With the exception of the ESTRO mailed dosimetry service, there is little dosimetric audit activity reported in the literature, when compared with external beam radiotherapy verification.

Fast patient-specific Monte Carlo brachytherapy dose calculations via the correlated sampling variance reduction technique

PURPOSE

To demonstrate potential of correlated sampling Monte Carlo (CMC) simulation to improve the calculation efficiency for permanent seed brachytherapy (PSB) implants without loss of accuracy.

METHODS

CMC was implemented within an in-house MC code family (PTRAN) and used to compute 3D dose distributions for two patient cases: a clinical PSB postimplant prostate CT imaging study and a simulated post lumpectomy breast PSB implant planned on a screening dedicated breast cone-beam CT patient exam. CMC tallies the dose difference, ΔD, between highly correlated histories in homogeneous and heterogeneous geometries. The heterogeneous geometry histories were derived from photon collisions sampled in a geometrically identical but purely homogeneous medium geometry, by altering their particle weights to correct for bias. The prostate case consisted of 78 Model-6711 (125)I seeds. The breast case consisted of 87 Model-200 (103)Pd seeds embedded around a simulated lumpectomy cavity. Systematic and random errors in CMC were unfolded using low-uncertainty uncorrelated MC (UMC) as the benchmark. CMC efficiency gains, relative to UMC, were computed for all voxels, and the mean was classified in regions that received minimum doses greater than 20%, 50%, and 90% of D(90), as well as for various anatomical regions.

RESULTS

Systematic errors in CMC relative to UMC were less than 0.6% for 99% of the voxels and 0.04% for 100% of the voxels for the prostate and breast cases, respectively. For a 1 × 1 × 1 mm(3) dose grid, efficiency gains were realized in all structures with 38.1- and 59.8-fold average gains within the prostate and breast clinical target volumes (CTVs), respectively. Greater than 99% of the voxels within the prostate and breast CTVs experienced an efficiency gain. Additionally, it was shown that efficiency losses were confined to low dose regions while the largest gains were located where little difference exists between the homogeneous and heterogeneous doses. On an AMD 1090T processor, computing times of 38 and 21 sec were required to achieve an average statistical uncertainty of 2% within the prostate (1 × 1 × 1 mm(3)) and breast (0.67 × 0.67 × 0.8 mm(3)) CTVs, respectively.

CONCLUSIONS

CMC supports an additional average 38-60 fold improvement in average efficiency relative to conventional uncorrelated MC techniques, although some voxels experience no gain or even efficiency losses. However, for the two investigated case studies, the maximum variance within clinically significant structures was always reduced (on average by a factor of 6) in the therapeutic dose range generally. CMC takes only seconds to produce an accurate, high-resolution, low-uncertainly dose distribution for the low-energy PSB implants investigated in this study.

Impact of heterogeneity-based dose calculation using a deterministic grid-based Boltzmann equation solver for intracavitary brachytherapy

PURPOSE

To investigate the dosimetric impact of the heterogeneity dose calculation Acuros (Transpire Inc., Gig Harbor, WA), a grid-based Boltzmann equation solver (GBBS), for brachytherapy in a cohort of cervical cancer patients.

METHODS AND MATERIALS

The impact of heterogeneities was retrospectively assessed in treatment plans for 26 patients who had previously received (192)Ir intracavitary brachytherapy for cervical cancer with computed tomography (CT)/magnetic resonance-compatible tandems and unshielded colpostats. The GBBS models sources, patient boundaries, applicators, and tissue heterogeneities. Multiple GBBS calculations were performed with and without solid model applicator, with and without overriding the patient contour to 1 g/cm(3) muscle, and with and without overriding contrast materials to muscle or 2.25 g/cm(3) bone. Impact of source and boundary modeling, applicator, tissue heterogeneities, and sensitivity of CT-to-material mapping of contrast were derived from the multiple calculations. American Association of Physicists in Medicine Task Group 43 (TG-43) guidelines and the GBBS were compared for the following clinical dosimetric parameters: Manchester points A and B, International Commission on Radiation Units and Measurements (ICRU) report 38 rectal and bladder points, three and nine o'clock, and (D2cm3) to the bladder, rectum, and sigmoid.

RESULTS

Points A and B, D(2) cm(3) bladder, ICRU bladder, and three and nine o'clock were within 5% of TG-43 for all GBBS calculations. The source and boundary and applicator account for most of the differences between the GBBS and TG-43 guidelines. The D(2cm3) rectum (n = 3), D(2cm3) sigmoid (n = 1), and ICRU rectum (n = 6) had differences of >5% from TG-43 for the worst case incorrect mapping of contrast to bone. Clinical dosimetric parameters were within 5% of TG-43 when rectal and balloon contrast were mapped to bone and radiopaque packing was not overridden.

CONCLUSIONS

The GBBS has minimal impact on clinical parameters for this cohort of patients with unshielded applicators. The incorrect mapping of rectal and balloon contrast does not have a significant impact on clinical parameters. Rectal parameters may be sensitive to the mapping of radiopaque packing.

The dosimetric impact of heterogeneity corrections in high-dose-rate ¹⁹²Ir brachytherapy for cervical cancer: Investigation of both conventional Point-A and volume-optimized plans

PURPOSE

To evaluate the dosimetric impact of heterogeneity corrections on both conventional and volume-optimized high-dose-rate (HDR) ¹⁹²Ir brachytherapy tandem-and-ovoid treatment plans.

METHODS AND MATERIALS

Both conventional and volume-optimized treatment plans were retrospectively created using eight unique CT data sets. In the volume-optimized plans, the clinical target volume (CTV) and organs-at-risk (rectum, bladder, and sigmoid) were contoured on the CT data sets by a single physician. For each plan, dose calculations representing homogeneous water medium were performed using the Task Group (TG-43) formalism and dose calculations with heterogeneity corrections were performed using a commercially available treatment planning system.

RESULTS

For the conventional plans, the change in dose between TG-43 and heterogeneity-corrected calculations was assessed for the following points: Point-A (left and right) and International Commission on Radiation Units and Measurements (ICRU) 38 defined rectum and bladder points. It was found that the dose to the ICRU bladder decreased the most (-2.2±0.9%), whereas ICRU rectum (-1.7±0.8%), Point-A right (-1.1±0.4%), and Point-A left (-1.0±0.3%) also showed decreases with heterogeneity-corrected calculations. For the volume-optimized plans, the change in dose between TG-43 and heterogeneity-corrected calculations was assessed for the following dose-volume histogram parameters: D(90) of the CTV and D(2cc) of the rectum, bladder, and sigmoid. It was found that D(90) of the CTV decreased by -1.9±0.7% and D(2cc) decreased by -2.6±1.4%, -1.0±0.4%, and -2.0±0.6% for the rectum, bladder and sigmoid, respectively, with heterogeneity-corrected calculations.

CONCLUSIONS

Heterogeneity corrections on high-dose rate plans were found to have only a small dosimetric impact over TG-43-based dose calculations for both conventional Point-A and volume-optimized plans.