Patient-Specific Monte Carlo-Based Dose-Kernel Approach for Inverse Planning in Afterloading Brachytherapy

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.

DVH-Based Inverse Planning Using Monte Carlo Dosimetry for LDR Prostate Brachytherapy

PURPOSE

Inverse planning is an integral part of modern low-dose-rate brachytherapy. Current clinical planning systems do not exploit the total dose information and largely use the American Association of Physicists in Medicine TG-43 dosimetry formalism to ensure clinically acceptable planning times. Thus, suboptimal plans may be derived as a result of TG-43-related dose overestimation and nonconformity with dose distribution requirements. The purpose of this study was to propose an inverse planning approach that can improve planning quality by combining dose-volume information and precision without compromising the overall execution times.

METHODS AND MATERIALS

The dose map was generated by accumulating precomputed Monte Carlo (MC) dose kernels for each candidate source implantation site. The MC computational burden was reduced by using graphics processing unit acceleration, allowing accurate dosimetry calculations to be performed in the intraoperative environment. The proposed dose-volume histogram (DVH) fast simulated annealing optimization algorithm was evaluated using clinical plans that were delivered to 18 patients who underwent low-dose-rate prostate brachytherapy.

RESULTS

Our method generated plans in 37.5 ± 3.2 seconds with similar prostate dose coverage, improved prostate dose homogeneity of up to 6.1%, and lower dose to the urethra of up to 4.0%.

CONCLUSIONS

A DVH-based optimization algorithm using MC dosimetry was developed. The inclusion of the DVH requirements allowed for increased control over the optimization outcome. The optimal plan's quality was further improved by considering tissue heterogeneity.

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.

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.