RapidBrachyMCTPS: a Monte Carlo-based treatment planning system for brachytherapy applications

Penalty weight tuning in high dose rate brachytherapy using multi-objective Bayesian optimization

Objective.Treatment plan optimization in high dose rate brachytherapy often requires manual fine-tuning of penalty weights for each objective, which can be time-consuming and dependent on the planner's experience. To automate this process, this study used a multi-criteria approach called multi-objective Bayesian optimization with q-noisy expected hypervolume improvement as its acquisition function (MOBO-qNEHVI).Approach.The treatment plans of 13 prostate cancer patients were retrospectively imported to a research treatment planning system, RapidBrachyMTPS, where fast mixed integer optimization (FMIO) performs dwell time optimization given a set of penalty weights to deliver 15 Gy to the target volume. MOBO-qNEHVI was used to find patient-specific Pareto optimal penalty weight vectors that yield clinically acceptable dose volume histogram metrics. The relationship between the number of MOBO-qNEHVI iterations and the number of clinically acceptable plans per patient (acceptance rate) was investigated. The performance time was obtained for various parameter configurations.Main results.MOBO-qNEHVI found clinically acceptable treatment plans for all patients. With increasing the number of MOBO-qNEHVI iterations, the acceptance rate grew logarithmically while the performance time grew exponentially. Fixing the penalty weight of the tumour volume to maximum value, adding the target dose as a parameter, initiating MOBO-qNEHVI with 25 parallel sampling of FMIO, and running 6 MOBO-qNEHVI iterations found solutions that delivered 15 Gy to the hottest 95% of the clinical target volume while respecting the dose constraints to the organs at risk. The average acceptance rate for each patient was 89.74% ± 8.11%, and performance time was 66.6 ± 12.6 s. The initiation took 22.47 ± 7.57 s, and each iteration took 7.35 ± 2.45 s to find one Pareto solution.Significance.MOBO-qNEHVI combined with FMIO can automatically explore the trade-offs between treatment plan objectives in a patient specific manner within a minute. This approach can reduce the dependency of plan quality on planner's experience and reduce dose to the organs at risk.

RapidBrachyTG43: A Geant4-based TG-43 parameter and dose calculation module for brachytherapy dosimetry

BACKGROUND

The AAPM TG-43U1 formalism remains the clinical standard for dosimetry of low- and high-energy γ $\gamma$ -emitting brachytherapy sources. TG-43U1 and related reports provide consensus datasets of TG-43 parameters derived from various published measured data and Monte Carlo simulations. These data are used to perform standardized and fast dose calculations for brachytherapy treatment planning.

PURPOSE

Monte Carlo TG-43 dosimetry parameters are commonly derived to characterize novel brachytherapy sources. RapidBrachyTG43 is a module of RapidBrachyMCTPS, a Monte Carlo-based treatment planning system, designed to automate this process, requiring minimal user input to prepare Geant4-based Monte Carlo simulations for a source. RapidBrachyTG43 may also perform a TG-43 dose to water-in-water calculation for a plan, substantially accelerating the same calculation performed using RapidBrachyMCTPS's Monte Carlo dose calculation engine.

METHODS

TG-43 parametersS K / A $S_K/A$ , Λ $\Lambda$ ,g L ( r ) $g_L(r)$ , andF ( r , θ ) $F(r,\theta)$ were calculated using three commercial source models, one each of125 $^{125}$ I,192 $^{192}$ Ir, and60 $^{60}$ Co, and were benchmarked to published data. TG-43 dose to water was calculated for a clinical breast brachytherapy plan and was compared to a Monte Carlo dose calculation with all patient tissues, air, and catheters set to water.

RESULTS

TG-43 parameters for the three simulated sources agreed with benchmark datasets within tolerances specified by the High Energy Brachytherapy Dosimetry working group. A gamma index comparison between the TG-43 and Monte Carlo dose-to-water calculations with a dose difference and difference to agreement criterion of 1%/1 mm yielded a 98.9% pass rate, with all relevant dose volume histogram metrics for the plan agreeing within 1%. Performing a TG-43-based dose calculation provided an acceleration of dose-to-water calculation by a factor of 165.

CONCLUSIONS

Determination of TG-43 parameter data for novel brachytherapy sources may now be facilitated by RapidBrachyMCTPS. These parameter datasets and existing consensus or published datasets may also be used to determine the TG-43 dose for a plan in RapidBrachyMCTPS.

Deep learning for high-resolution dose prediction in high dose rate brachytherapy for breast cancer treatment

Objective.Monte Carlo (MC) simulations are the benchmark for accurate radiotherapy dose calculations, notably in patient-specific high dose rate brachytherapy (HDR BT), in cases where considering tissue heterogeneities is critical. However, the lengthy computational time limits the practical application of MC simulations. Prior research used deep learning (DL) for dose prediction as an alternative to MC simulations. While accurate dose predictions akin to MC were attained, graphics processing unit limitations constrained these predictions to large voxels of 3 mm × 3 mm × 3 mm. This study aimed to enable dose predictions as accurate as MC simulations in 1 mm × 1 mm × 1 mm voxels within a clinically acceptable timeframe.Approach.Computed tomography scans of 98 breast cancer patients treated with Iridium-192-based HDR BT were used: 70 for training, 14 for validation, and 14 for testing. A new cropping strategy based on the distance to the seed was devised to reduce the volume size, enabling efficient training of 3D DL models using 1 mm × 1 mm × 1 mm dose grids. Additionally, novel DL architecture with layer-level fusion were proposed to predict MC simulated dose to medium-in-medium (Dm,m). These architectures fuse information from TG-43 dose to water-in-water (Dw,w) with patient tissue composition at the layer-level. Different inputs describing patient body composition were investigated.Main results.The proposed approach demonstrated state-of-the-art performance, on par with the MCDm,mmaps, but 300 times faster. The mean absolute percent error for dosimetric indices between the MC and DL-predicted complete treatment plans was 0.17% ± 0.15% for the planning target volumeV100, 0.30% ± 0.32% for the skinD2cc, 0.82% ± 0.79% for the lungD2cc, 0.34% ± 0.29% for the chest wallD2ccand 1.08% ± 0.98% for the heartD2cc.Significance.Unlike the time-consuming MC simulations, the proposed novel strategy efficiently converts TG-43Dw,wmaps into preciseDm,mmaps at high resolution, enabling clinical integration.

Embedding expertise knowledge into inverse treatment planning for low-dose-rate brachytherapy of hepatic malignancies

BACKGROUND

Leveraging the precision of its radiation dose distribution and the minimization of postoperative complications, low-dose-rate (LDR) permanent seed brachytherapy is progressively adopted in addressing hepatic malignancies.

PURPOSE

The present study endeavors to devise a sophisticated treatment planning system (TPS) to optimize LDR brachytherapy for hepatic lesions.

METHODS

Our TPS encompasses four integral modules: multi-organ segmentation, seed distribution initialization, puncture pathway selection, and inverse dose planning. By amalgamating an array of deep learning models, the segmentation module proficiently labels 17 discrete abdominal targets within the images. We introduce a knowledge-based seed distribution initialization methodology that discerns the most analogous tumor shape in the reference treatment plan from the knowledge base. Subsequently, the seed distribution from the reference plan is transmuted to the current case, thus establishing seed distribution initialization. Furthermore, we parameterize the puncture needles and seeds, while concurrently constraining the puncture needle angle through the employment of a virtual puncture panel to augment planning algorithm efficiency. We also presented a user interface that includes a range of interactive features, seamlessly integrated with the treatment planning generation function.

RESULTS

The multi-organ segmentation module, which is trained by 50 cases of in-house CT scans and 694 cases of publicly available CT scans, achieved average Dice of 0.80 and Hausdorff distance of 5.2 mm in testing datasets. The results demonstrate that knowledge-based initialization exhibits a marked enhancement in expediting the convergence rate. Our TPS also demonstrates a dominant advantage in dose-volume-histogram criteria and execution time in comparison to commercial TPS.

CONCLUSION

The study proposes an innovative treatment planning system for low-dose-rate permanent seed brachytherapy for hepatic malignancies. We show that the generated treatment plans meet clinical requirement.

Benchmark of the PenRed Monte Carlo framework for HDR brachytherapy

PURPOSE

The purpose of this study is to validate the PenRed Monte Carlo framework for clinical applications in brachytherapy. PenRed is a C++ version of Penelope Monte Carlo code with additional tallies and utilities.

METHODS AND MATERIALS

Six benchmarking scenarios are explored to validate the use of PenRed and its improved bachytherapy-oriented capabilities for HDR brachytherapy. A new tally allowing the evaluation of collisional kerma for any material using the track length kerma estimator and the possibility to obtain the seed positions, weights and directions processing directly the DICOM file are now implemented in the PenRed distribution. The four non-clinical test cases developed by the Joint AAPM-ESTRO-ABG-ABS WG-DCAB were evaluated by comparing local and global absorbed dose differences with respect to established reference datasets. A prostate and a palliative lung cases, were also studied. For them, absorbed dose ratios, global absorbed dose differences, and cumulative dose-volume histograms were obtained and discussed.

RESULTS

The air-kerma strength and the dose rate constant corresponding to the two sources agree with the reference datatests within 0.3% (Sk) and 0.1% (Λ). With respect to the first three WG-DCAB test cases, more than 99.8% of the voxels present local (global) differences within ±1%(±0.1%) of the reference datasets. For test Case 4 reference dataset, more than 94.9%(97.5%) of voxels show an agreement within ±1%(±0.1%), better than similar benchmarking calculations in the literature. The track length kerma estimator scorer implemented increases the numerical efficiency of brachytherapy calculations two orders of magnitude, while the specific brachytherapy source allows the user to avoid the use of error-prone intermediate steps to translate the DICOM information into the simulation. In both clinical cases, only minor absorbed dose differences arise in the low-dose isodoses. 99.8% and 100% of the voxels have a global absorbed dose difference ratio within ±0.2% for the prostate and lung cases, respectively. The role played by the different segmentation and composition material in the bone structures was discussed, obtaining negligible absorbed dose differences. Dose-volume histograms were in agreement with the reference data.

CONCLUSIONS

PenRed incorporates new tallies and utilities and has been validated for its use for detailed and precise high-dose-rate brachytherapy simulations.

Applying the column generation method to the intensity modulated high dose rate brachytherapy inverse planning problem

Objective.Intensity modulated high dose rate brachytherapy (IMBT) is a rapidly developing application of brachytherapy where anisotropic dose distributions can be produced at each source dwell position. This technique is made possible by placing rotating metallic shields inside brachytherapy needles or catheters. By dynamically directing the radiation towards the tumours and away from the healthy tissues, a more conformal dose distribution can be obtained. The resulting treatment planning involves optimizing dwell position and shield angle (DPSA). The aim of this study was to investigate the column generation method for IMBT treatment plan optimization.Approach.A column generation optimization algorithm was developed to optimize the dwell times and shield angles. A retrospective study was performed on 10 prostate cases using RapidBrachyMCTPS. At every iteration, the plan was optimized with the chosen DPSA which would best improve the cost function that was added to the plan. The optimization process was stopped when the remaining DPSAs would not add value to the plan to limit the plan complexity.Main results.The average number of DPSAs and voxels were 2270 and 7997, respectively. The column generation approach yielded near-optimal treatment plans by using only 11% of available DPSAs on average in ten prostate cases. The coverage and organs at risk constraints passed in all ten cases.Significance.The column generation method produced high-quality deliverable prostate IMBT plans. The treatment plan quality reached a plateau, where adding more DPSAs had a minimal effect on dose volume histogram parameters. The iterative nature of the column generation method allows early termination of the treatment plan creation process as soon as the dosimetric indices from dose volume histogram satisfy the clinical requirements or if their values stabilize.

Validation of the TOPAS Monte Carlo toolkit for LDR brachytherapy simulations

PURPOSE

This work has the purpose of validating the Monte Carlo toolkit TOol for PArticle Simulation (TOPAS) for low-dose-rate (LDR) brachytherapy uses.

METHODS AND MATERIALS

Simulations of 12 LDR sources and 2 COMS eye plaques (10 mm and 20 mm in diameter) and comparisons with published reference data from the Carleton Laboratory for Radiotherapy Physics (CLRP), the TG-43 consensus data and the TG-129 consensus data were performed. Sources from the IROC Houston Source Registry were modeled. The OncoSeed 6711 and the SelectSeed 130.002 were also modeled for historical reasons. For each source, the dose rate constant, the radial dose function and the anisotropy functions at 0.5, 1 and 5 cm were extracted. For the eye plaques (loaded with ¹²⁵I sources), dose distribution maps, dose profiles along the central axis and transverse axis were calculated.

RESULTS

Dose rate constants for 11 of the 12 sources are within 4% of the consensus data and within 2% of the CLRP data. The radial dose functions and anisotropy functions are mostly within 2% of the CLRP data. In average, 92% of all voxels are within 1% of the CLRP data for the eye plaques dose distributions. The dose profiles are within 0.5% (central axis) and 1% (transverse axis) of the reference data.

CONCLUSION

The TOPAS MC toolkit was validated for LDR brachytherapy applications. Single-seed and multi-seed results agree with the published reference data. TOPAS has several benefits such as a simplified approach to MC simulations and an accessible brachytherapy package including comprehensive learning resources.

Dosimetric impact of a robust optimization approach to mitigate effects from rotational uncertainty in prostate intensity-modulated brachytherapy

BACKGROUND

Intensity-modulated brachytherapy (IMBT) is an emerging technology for cancer treatment, in which radiation sources are shielded to shape the dose distribution. The rotatable shields provide an additional degree of freedom, but also introduce an additional, directional, type of uncertainty, compared to conventional high-dose-rate brachytherapy (HDR BT).

PURPOSE

We propose and evaluate a robust optimization approach to mitigate the effects of rotational uncertainty in the shields with respect to planning criteria.

METHODS

A previously suggested prototype for platinum-shielded prostate ¹⁶⁹ Yb-based dynamic IMBT is considered. We study a retrospective patient data set (anatomical contours and catheter placement) from two clinics, consisting of six patients that had previously undergone conventional ¹⁹² Ir HDR BT treatment. The Monte Carlo-based treatment planning software RapidBrachyMCTPS is used for dose calculations. In our computational experiments, we investigate systematic rotational shield errors of ±10° and ±20°, and the same systematic error is applied to all dwell positions in each scenario. This gives us three scenarios, one nominal and two with errors. The robust optimization approach finds a compromise between the average and worst-case scenario outcomes.

RESULTS

We compare dose plans obtained from standard models and their robust counterparts. With dwell times obtained from a linear penalty model (LPM), for 10° errors, the dose to urethra ( D 0.1 c c $D_{0.1cc}$ ) and rectum ( D 0.1 c c $D_{0.1cc}$ and D 1 c c $D_{1cc}$ ) increase with up to 5% and 7%, respectively, in the worst-case scenario, while with the robust counterpart, the corresponding increases were 3% and 3%. For all patients and all evaluated criteria, the worst-case scenario outcome with the robust approach had lower deviation compared to the standard model, without compromising target coverage. We also evaluated shield errors up to 20° and while the deviations increased to a large extent with the standard models, the robust models were capable of handling even such large errors.

CONCLUSIONS

We conclude that robust optimization can be used to mitigate the effects from rotational uncertainty and to ensure the treatment plan quality of IMBT.

Emerging technologies in brachytherapy

Brachytherapy is a mature treatment modality. The literature is abundant in terms of review articles and comprehensive books on the latest established as well as evolving clinical practices. The intent of this article is to part ways and look beyond the current state-of-the-art and review emerging technologies that are noteworthy and perhaps may drive the future innovations in the field. There are plenty of candidate topics that deserve a deeper look, of course, but with practical limits in this communicative platform, we explore four topics that perhaps is worthwhile to review in detail at this time. First, intensity modulated brachytherapy (IMBT) is reviewed. The IMBT takes advantage ofanisotropicradiation profile generated through intelligent high-density shielding designs incorporated onto sources and applicators such to achieve high quality plans. Second, emerging applications of 3D printing (i.e. additive manufacturing) in brachytherapy are reviewed. With the advent of 3D printing, interest in this technology in brachytherapy has been immense and translation swift due to their potential to tailor applicators and treatments customizable to each individual patient. This is followed by, in third, innovations in treatment planning concerning catheter placement and dwell times where new modelling approaches, solution algorithms, and technological advances are reviewed. And, fourth and lastly, applications of a new machine learning technique, called deep learning, which has the potential to improve and automate all aspects of brachytherapy workflow, are reviewed. We do not expect that all ideas and innovations reviewed in this article will ultimately reach clinic but, nonetheless, this review provides a decent glimpse of what is to come. It would be exciting to monitor as IMBT, 3D printing, novel optimization algorithms, and deep learning technologies evolve over time and translate into pilot testing and sensibly phased clinical trials, and ultimately make a difference for cancer patients. Today's fancy is tomorrow's reality. The future is bright for brachytherapy.

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.

How to design, fabricate, and validate a customized COMS-style eye plaque: Illustrated with a narrow-slotted plaque example

PURPOSE

A customized Collaborative Ocular Melanoma Study (COMS)-style eye plaque may provide superior dosimetric coverage compared with standard models for certain intraocular tumor locations and shapes. This work provides a recipe for developing and validating such customized plaques.

METHODS AND MATERIALS

The concept-into-clinical treatment process for a customized COMS-style eye plaque begins with a CAD model design that meets the specifications of the radiation oncologist and surgeon based on magnetic resonance, ultrasound, and clinical measurements, as well as a TG-43 hybrid heterogeneity-corrected dose prediction to model the dose distribution. Next, a 3D printed plastic prototype is created and reviewed. After design approval, a Modulay plaque is commercially fabricated. Quality assurance (QA) is subsequently performed to verify the physical measurements of the Modulay and Silastic and also includes dosimetric measurement of the calibration, depth dose, and dose profiles. Sterilization instructions are provided by the commercial fabricator. This customization procedure and QA methodology is demonstrated with a narrow-slotted plaque that was recently constructed for the treatment of a circumpapillary (e.g., surrounding the optic disk) ocular tumor.

RESULTS

The production of a customized COMS-style eye plaque is a multistep process. Dosimetric modeling is recommended to ensure that the design will meet the patient's needs, and QA is essential to confirm that the plaque has the proper dimensions and dose distribution. The customized narrow-slotted plaque presented herein was successfully implemented in the clinic, and provided superior dose coverage of juxtapapillary and circumpapillary tumors compared with standard or notched COMS-style plaques. Plaque development required approximately 30 h of physicist time and a fabrication cost of $1500.

CONCLUSION

Customized eye plaques may be used to treat intraocular tumors that cannot be adequately managed with standard models. The procedure by which a customized COMS-style plaque may be designed, fabricated, and validated was presented along with a clinical example.

Monte Carlo study of TG-43 dosimetry parameters of GammaMed Plus high dose rate 192 Ir brachytherapy source using TOPAS

PURPOSE

To develop a simulation model for GammaMed Plus high dose rate ¹⁹² Ir brachytherapy source in TOPAS Monte Carlo software and validate it by calculating the TG-43 dosimetry parameters and comparing them with published data.

METHODS

We built a model for GammaMed Plus high dose rate brachytherapy source in TOPAS. The TG-43 dosimetry parameters including air-kerma strength S_K , dose-rate constant Λ, radial dose function g_L (r), and 2D anisotropy function F(r,θ) were calculated using Monte Carlo simulation with Geant4 physics models and NNDC ¹⁹² Ir spectrum. Calculations using an old ¹⁹² Ir spectrum were also carried out to evaluate the impact of incident spectrum and cross sections. The results were compared with published data.

RESULTS

For calculations using the NNDC spectrum, the air-kerma strength per unit source activity S_K /A and Λ were 1.0139 × 10^-7 U/Bq and 1.1101 cGy.h^-1 .U^-1 , which were 3.56% higher and 0.62% lower than the reference values, respectively. The g_L (r) agreed with reference values within 1% for radial distances from 2 mm to 20 cm. For radial distances of 1, 3, 5, and 10 cm, the agreements between F(r,θ) from this work and the reference data were within 1.5% for 15° < θ < 165°, and within 4% for all θ values. The discrepancies were attributed to the updated source spectrum and cross sections. They caused deviations of the S_K /A of 2.90% and 0.64%, respectively. As for g_L (r), they caused average deviations of -0.22% and 0.48%, respectively. Their impact on F(r,θ) was not quantified for the relatively high statistical uncertainties, but basically they did not result in significant discrepancies.

CONCLUSION

A model for GammaMed Plus high dose rate ¹⁹² Ir brachytherapy source was developed in TOPAS and validated following TG-43 protocols, which can be used for future studies. The impact of updated incident spectrum and cross sections on the dosimetry parameters was quantified.

On the impact of absorbed dose specification, tissue heterogeneities, and applicator heterogeneities on Monte Carlo-based dosimetry of Ir-192, Se-75, and Yb-169 in conventional and intensity-modulated brachytherapy for the treatment of cervical cancer

PURPOSE

The purpose of this study was to evaluate the impact of dose reporting schemes and tissue/applicator heterogeneities for ¹⁹² Ir-, ⁷⁵ Se-, and ¹⁶⁹ Yb-based MRI-guided conventional and intensity-modulated brachytherapy.

METHODS AND MATERIALS

Treatment plans using a variety of dose reporting and tissue/applicator segmentation schemes were generated for a cohort (n = 10) of cervical cancer patients treated with ¹⁹² Ir-based Venezia brachytherapy. Dose calculations were performed using RapidBrachyMCTPS, a Geant4-based research Monte Carlo treatment planning system. Ultimately, five dose calculation scenarios were evaluated: (a) dose to water in water (D_w,w ); (b) D_w,w taking the applicator material into consideration (D_w,wApp ); (c) dose to water in medium (D_w,m ); (d and e) dose to medium in medium with mass densities assigned either nominally per structure (D_{m,m (Nom)} ) or voxel-by-voxel (D_m,m ).

RESULTS

Ignoring the plastic Venezia applicator (D_w,wApp ) overestimates D_m,m by up to 1% (average) with high energy source (¹⁹² Ir and ⁷⁵ Se) and up to 2% with ¹⁶⁹ Yb. Scoring dose to water (D_w,wApp or D_w,m ) generally overestimates dose and this effect increases with decreasing photon energy. Reporting dose other than D_m,m (or D_{m,m Nom} ) for ¹⁶⁹ Yb-based conventional and intensity-modulated brachytherapy leads to a simultaneous overestimation (up to 4%) of CTV_HR D₉₀ and underestimation (up to 2%) of bladder D_2cc due to a significant dip in the mass-energy absorption ratios at the depths of nearby targets and OARs. Using a nominal mass-density assignment per structure, rather than a CT-derived voxel-by-voxel assignment for MRI-guided brachytherapy, amounts to a dose error up to 1% for all radionuclides considered.

CONCLUSIONS

The effects of the considered dose reporting schemes trend correspondingly between conventional and intensity-modulated brachytherapy. In the absence of CT-derived mass densities, MRI-only-based dosimetry can adequately approximate D_m,m by assigning nominal mass densities to structures. Tissue and applicator heterogeneities do not significantly impact dosimetry for ¹⁹² Ir and ⁷⁵ Se, but do for ¹⁶⁹ Yb; dose reporting must be explicitly defined since D_w,m and D_w,w may overstate the dosimetric benefits.

A novel and fast methodology to calculate doses in LDR brachytherapy

We present the concept of a new methodology for faster simulation of the doses in brachytherapy with permanent implants, based on the knowledge of the seeds arrangement, adding previously simulated doses in an equivalent medium in terms of the atomic composition of the organ in question. To perform the doses calculations we use Monte Carlo simulations. We simulated a cylindrical I-125 seed and compared our results against published data. Our proposal is to have the doses simulated previously in different arrangement of seed-absorbents, and then, considering the spacial positions of the seeds after the implants, these doses can be directly added, obtaining a very fast computation of the total dose. Two phantoms of prostates with permanent implant seeds in 2D and 3D arrangements were simulated. The results of the proposed methodology were compared with two complete Monte Carlo simulations in 2D and 3D designs. Differences in doses were analysed, obtaining statistical discrepancies of less than 1% and reducing the simulation time by more than 4 orders of magnitude. With the proposed methodology, it is possible to perform rapid dose calculations in brachytherapy, using laptop or desktop computers.

A novel minimally invasive dynamic-shield, intensity-modulated brachytherapy system for the treatment of cervical cancer

PURPOSE

To present a novel, MRI-compatible dynamicshield intensity modulated brachytherapy (IMBT) applicator and delivery system using ¹⁹² Ir, ⁷⁵ Se, and ¹⁶⁹ Yb radioisotopes for the treatment of locally advanced cervical cancer. Needle-free IMBT is a promising technique for improving target coverage and organs at risk (OAR) sparing.

METHODS AND MATERIALS

The IMBT delivery system dynamically controls the rotation of a novel tungsten shield placed inside an MRI-compatible, 6-mm wide intrauterine tandem. Using 36 cervical cancer cases, conventional intracavitary brachytherapy (IC-BT) and intracavitary/interstitial brachytherapy (IC/IS-BT) (10Ci ¹⁹² Ir) plans were compared to IMBT (10Ci ¹⁹² Ir; 11.5Ci ⁷⁵ Se; 44Ci ¹⁶⁹ Yb). All plans were generated using the Geant4-based Monte Carlo dose calculation engine, RapidBrachyMC. Treatment plans were optimized then normalized to the same high-risk clinical target volume (HR-CTV) D₉₀ and the D_2cc for bladder, rectum, and sigmoid in the research brachytherapy planning system, RapidBrachyMCTPS. Plans were renormalized until either of the three OAR reached dose limits to calculate the maximum achievable HR-CTV D₉₀ and D₉₈ .

RESULTS

Compared to IC-BT, IMBT with either of the three radionuclides significantly improves the HR-CTV D₉₀ and D₉₈ by up to 5.2% ± 0.3% (P < 0.001) and 6.7% ± 0.5% (P < 0.001), respectively, with the largest dosimetric enhancement when using ¹⁶⁹ Yb followed by ⁷⁵ Se and then ¹⁹² Ir. Similarly, D_2cc for all OAR improved with IMBT by up to 7.7% ± 0.6% (P < 0.001). For IC/IS-BT cases, needle-free IMBT achieved clinically acceptable plans with ¹⁶⁹ Yb-based IMBT further improving HR-CTV D₉₈ by 1.5% ± 0.2% (P = 0.034) and decreasing sigmoid D_2cc by 1.9% ± 0.4% (P = 0.048). Delivery times for IMBT are increased by a factor of 1.7, 3.3, and 2.3 for ¹⁹² Ir, ⁷⁵ Se, and ¹⁶⁹ Yb, respectively, relative to conventional ¹⁹² Ir BT.

CONCLUSIONS

Dynamic shield IMBT provides a promising alternative to conventional IC- and IC/IS-BT techniques with significant dosimetric enhancements and even greater improvements with intermediate energy radionuclides. The ability to deliver a highly conformal, OAR-sparing dose without IS needles provides a simplified method for improving the therapeutic ratio less invasively and in a less resource intensive manner.

Monte Carlo dosimetric characterization of a new high dose rate 169 Yb brachytherapy source and independent verification using a multipoint plastic scintillator detector

PURPOSE

A prototype 169 Yb source was developed in combination with a dynamic rotating platinum shield system (AIM-Brachy) to deliver intensity modulated brachytherapy (IMBT). The purpose of this study was to evaluate the dosimetric characteristics of the bare/shielded 169 Yb source using Monte Carlo (MC) simulations and perform an independent dose verification using a dosimetry platform based on a multipoint plastic scintillator detector (mPSD).

METHODS

The TG-43U1 dosimetric parameters were calculated for the source model using RapidBrachyMCTPS. Real-time dose rate measurements were performed in a water tank for both the bare/shielded source using a custom remote afterloader. For each dwell position, the dose rate was independently measured by the three scintillators (BCF-10, BCF-12, and BCF-60). For the bare source, dose rate was measured at distances up to 3 cm away from the source over a range of 7 cm along the catheter. For the shielded source, measurements were performed with the mPSD placed at 1 cm from the source at four different azimuthal angles ( 0 ∘ , 9 0 ∘ , 18 0 ∘ , and 27 0 ∘ ).

RESULTS

The dosimetric parameters were tabulated for the source model. For the bare source, differences between measured and calculated along-away dose rates were generally below 5-10%. Along the transverse axis, deviations were, on average (range), 3.3% (0.6-6.2%) for BCF-10, 1.7% (0.9-2.9%) for BCF-12, and 2.2% (0.3-4.4%) for BCF-60. The maximum dose rate reduction due to shielding at a radial distance of 1 cm was 88.8 ± 1.2%, compared to 83.5 ± 0.5% as calculated by MC.

CONCLUSIONS

The dose distribution for the bare/shielded 169 Yb source was independently verified using mPSD with good agreement in regions close to the source. The 169 Yb source coupled with the partial-shielding system is an effective technique to deliver IMBT.

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.