|
Ballester, F., Tedgren, A. C., Granero, D., Haworth, A., Mourtada, F., Fonseca, G. P., et al. (2015). A generic high-dose rate Ir-192 brachytherapy source for evaluation of model-based dose calculations beyond the TG-43 formalism. Med. Phys., 42(6), 3048–3062.
Abstract: 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) Ir-192 source and a virtual water phantom were designed, which can be imported into a TPS. Methods: A hypothetical, generic HDR Ir-192 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 Ir-192 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 (R) Brachy with advanced collapsed-cone engine (ACE) and BrachyVision AcuRos (TM)]. 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 Ir-192 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 Ir-192 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.
|
|
|
Candela-Juan, C., Niatsetski, Y., van der Laarse, R., Granero, D., Ballester, F., Perez-Calatayud, J., et al. (2016). Design and characterization of a new high-dose-rate brachytherapy Valencia applicator for larger skin lesions. Med. Phys., 43(4), 1639–1648.
Abstract: Purpose: The aims of this study were (i) to design a new high-dose-rate (HDR) brachytherapy applicator for treating surface lesions with planning target volumes larger than 3 cm in diameter and up to 5 cm in size, using the microSelectron-HDR or Flexitron afterloader (Elekta Brachytherapy) with a Ir-192 source; (ii) to calculate by means of the Monte Carlo (MC) method the dose distribution for the new applicator when it is placed against a water phantom; and (iii) to validate experimentally the dose distributions in water. Methods: The PENELOPE2008 MC code was used to optimize dwell positions and dwell times. Next, the dose distribution in a water phantom and the leakage dose distribution around the applicator were calculated. Finally, MC data were validated experimentally for a 192Ir mHDR-v2 source by measuring (i) dose distributions with radiochromic EBT3 films (ISP); (ii) percentage depth-dose (PDD) curve with the parallel-plate ionization chamber Advanced Markus (PTW); and (iii) absolute dose rate with EBT3 films and the PinPoint T31016 (PTW) ionization chamber. Results: The new applicator is made of tungsten alloy (Densimet) and consists of a set of interchangeable collimators. Three catheters are used to allocate the source at prefixed dwell positions with preset weights to produce a homogenous dose distribution at the typical prescription depth of 3 mm in water. The same plan is used for all available collimators. PDD, absolute dose rate per unit of air kerma strength, and off-axis profiles in a cylindrical water phantom are reported. These data can be used for treatment planning. Leakage around the applicator was also scored. The dose distributions, PDD, and absolute dose rate calculated agree within experimental uncertainties with the doses measured: differences of MC data with chamber measurements are up to 0.8% and with radiochromic films are up to 3.5%. Conclusions: The new applicator and the dosimetric data provided here will be a valuable tool in clinical practice, making treatment of large skin lesions simpler, faster, and safer. Also the dose to surrounding healthy tissues is minimal.
|
|
|
Gimenez-Alventosa, V., Vijande, J., Ballester, F., & Perez-Calatayud, J. (2016). Transit dose comparisons for Co-60 and Ir-192 HDR sources. J. Radiol. Prot., 36(4), 858–864.
Abstract: The goal of this study is to evaluate the ambient dose due to the transit of high dose rate (HDR) Co-60 sources along a transfer tube as compared to Ir-192 ones in a realistic clinical scenario. This goal is accomplished by evaluating air-kerma differences with Monte Carlo calculations using PENELOPE2011. Scatter from both the afterloader and the patient was not taken into account. Two sources, mHDR-v2 and Flexisource Co-60, (Elekta Brachytherapy, Veenendaal, the Netherlands) have been considered. These sources were simulated within a standard transfer tube located in an infinite air phantom. The movement of the source was included by displacing their positions along the connecting tube from z = – 75 cm to z = + 75 cm and combining them. Since modern afterloaders like Flexitron (Elekta) or Saginova (BEBIG GmbH) are able to use equally 192Ir and 60Co sources, it was assumed that both sources are displaced with equal speed. Typical HDR source activity content values were provided by the manufacturer. 2D distributions were obtained with type-A uncertainties (k = 2) less than 0.01%. From those, the air-kerma ratio Co-60/Ir-192 was evaluated weighted by their corresponding typical activities. It was found that it varies slowly with distance (less than 10% variation at 75 cm) but strongly in time due to the shorter half-life of the 192Ir (73.83 d). The maximum ratio is located close to the tube. It reaches a value of 0.57 when the typical activity of the sources at the time when they were installed by the vendor was used. Such ratio increases up to 1.28 at the end of the recommended working life (90 d) of the Ir-192 source. Co-60/Ir-192 air-kerma ratios are almost constant (0.51-0.57) in the vicinity of the source-tube with recent installed sources. Nevertheless, air-kerma ratios increase rapidly (1.15-1.29) whenever the Ir-192 is approaching the end of its life. In case of a medical event requiring the medical staff to access the treatment room, these ratios indicate that the dosimetric impact on the medical team will be lower, with a few exceptions, in the case of Co-60-based HDR brachytherapy as compared to Ir-192-based one when typical air-kerma strength values are considered.
|
|
|
Ma, Y. Z., Vijande, J., Ballester, F., Tedgren, A. C., Granero, D., Haworth, A., et al. (2017). A generic TG-186 shielded applicator for commissioning model-based dose calculation algorithms for high-dose-rate Ir-192 brachytherapy. Med. Phys., 44(11), 5961–5976.
Abstract: PurposeA joint working group was created by the American Association of Physicists in Medicine (AAPM), the European Society for Radiotherapy and Oncology (ESTRO), and the Australasian Brachytherapy Group (ABG) with the charge, among others, to develop a set of well-defined test case plans and perform calculations and comparisons with model-based dose calculation algorithms (MBDCAs). Its main goal is to facilitate a smooth transition from the AAPM Task Group No. 43 (TG-43) dose calculation formalism, widely being used in clinical practice for brachytherapy, to the one proposed by Task Group No. 186 (TG-186) for MBDCAs. To do so, in this work a hypothetical, generic high-dose rate (HDR) Ir-192 shielded applicator has been designed and benchmarked. MethodsA generic HDR Ir-192 shielded applicator was designed based on three commercially available gynecological applicators as well as a virtual cubic water phantom that can be imported into any DICOM-RT compatible treatment planning system (TPS). The absorbed dose distribution around the applicator with the TG-186 Ir-192 source located at one dwell position at its center was computed using two commercial TPSs incorporating MBDCAs (Oncentra((R)) Brachy with Advanced Collapsed-cone Engine, ACE, and BrachyVision ACUROS) and state-of-the-art Monte Carlo (MC) codes, including ALGEBRA, BrachyDose, egs_brachy, Geant4, MCNP6, and Penelope2008. TPS-based volumetric dose distributions for the previously reported source centered in water and source displaced test cases, and the new source centered in applicator test case, were analyzed here using the MCNP6 dose distribution as a reference. Volumetric dose comparisons of TPS results against results for the other MC codes were also performed. Distributions of local and global dose difference ratios are reported. ResultsThe local dose differences among MC codes are comparable to the statistical uncertainties of the reference datasets for the source centered in water and source displaced test cases and for the clinically relevant part of the unshielded volume in the source centered in applicator case. Larger local differences appear in the shielded volume or at large distances. Considering clinically relevant regions, global dose differences are smaller than the local ones. The most disadvantageous case for the MBDCAs is the one including the shielded applicator. In this case, ACUROS agrees with MC within [-4.2%, +4.2%] for the majority of voxels (95%) while presenting dose differences within [-0.12%, +0.12%] of the dose at a clinically relevant reference point. For ACE, 95% of the total volume presents differences with respect to MC in the range [-1.7%, +0.4%] of the dose at the reference point. ConclusionsThe combination of the generic source and generic shielded applicator, together with the previously developed test cases and reference datasets (available in the Brachytherapy Source Registry), lay a solid foundation in supporting uniform commissioning procedures and direct comparisons among treatment planning systems for HDR Ir-192 brachytherapy.
|
|
|
Peppa, V., Thomson, R. M., Enger, S. A., Fonseca, G. P., Lee, C. N., Lucero, J. N. E., et al. (2023). A MC-based anthropomorphic test case for commissioning model-based dose calculation in interstitial breast 192-Ir HDR brachytherapy. Med. Phys., 50(7), 4675–4687.
Abstract: PurposeTo provide the first clinical test case for commissioning of Ir-192 brachytherapy model-based dose calculation algorithms (MBDCAs) according to the AAPM TG-186 report workflow. Acquisition and Validation MethodsA computational patient phantom model was generated from a clinical multi-catheter Ir-192 HDR breast brachytherapy case. Regions of interest (ROIs) were contoured and digitized on the patient CT images and the model was written to a series of DICOM CT images using MATLAB. The model was imported into two commercial treatment planning systems (TPSs) currently incorporating an MBDCA. Identical treatment plans were prepared using a generic Ir-192 HDR source and the TG-43-based algorithm of each TPS. This was followed by dose to medium in medium calculations using the MBDCA option of each TPS. Monte Carlo (MC) simulation was performed in the model using three different codes and information parsed from the treatment plan exported in DICOM radiation therapy (RT) format. Results were found to agree within statistical uncertainty and the dataset with the lowest uncertainty was assigned as the reference MC dose distribution. Data Format and Usage NotesThe dataset is available online at ,. Files include the treatment plan for each TPS in DICOM RT format, reference MC dose data in RT Dose format, as well as a guide for database users and all files necessary to repeat the MC simulations. Potential ApplicationsThe dataset facilitates the commissioning of brachytherapy MBDCAs using TPS embedded tools and establishes a methodology for the development of future clinical test cases. It is also useful to non-MBDCA adopters for intercomparing MBDCAs and exploring their benefits and limitations, as well as to brachytherapy researchers in need of a dosimetric and/or a DICOM RT information parsing benchmark. Limitations include specificity in terms of radionuclide, source model, clinical scenario, and MBDCA version used for its preparation.
|
|
|
Valdes-Cortez, C., Niatsetski, Y., Perez-Calatayud, J., Ballester, F., & Vijande, J. (2022). A Monte Carlo study of the relative biological effectiveness in surface brachytherapy. Med. Phys., 49, 5576–5588.
Abstract: Purpose This work aims to simulate clustered DNA damage from ionizing radiation and estimate the relative biological effectiveness (RBE) for radionuclide (rBT)- and electronic (eBT)-based surface brachytherapy through a hybrid Monte Carlo (MC) approach, using realistic models of the sources and applicators. Methods Damage from ionizing radiation has been studied using the Monte Carlo Damage Simulation algorithm using as input the primary electron fluence simulated using a state-of-the-art MC code, PENELOPE-2018. Two Ir-192 rBT applicators, Valencia and Leipzig, one Co-60 source with a Freiburg Flap applicator (reference source), and two eBT systems, Esteya and INTRABEAM, have been included in this study implementing full realizations of their geometries as disclosed by the manufacturer. The role played by filtration and tube kilovoltage has also been addressed. Results For rBT, an RBE value of about 1.01 has been found for the applicators and phantoms considered. In the case of eBT, RBE values for the Esteya system show an almost constant RBE value of about 1.06 for all depths and materials. For INTRABEAM, variations in the range of 1.12-1.06 are reported depending on phantom composition and depth. Modifications in the Esteya system, filtration, and tube kilovoltage give rise to variations in the same range. Conclusions Current clinical practice does not incorporate biological effects in surface brachytherapy. Therefore, the same absorbed dose is administered to the patients independently on the particularities of the rBT or eBT system considered. The almost constant RBE values reported for rBT support that assumption regardless of the details of the patient geometry, the presence of a flattening filter in the applicator design, or even significant modifications in the photon energy spectra above 300 keV. That is not the case for eBT, where a clear dependence on the eBT system and the characteristics of the patient geometry are reported. A complete study specific for each eBT system, including detailed applicator characteristics (size, shape, filtering, among others) and common anatomical locations, should be performed before adopting an existing RBE value.
|
|