Dosimetric perturbations of a lead shield for surface and interstitial high-dose-rate brachytherapy.

In surface and interstitial high-dose-rate brachytherapy with either (60)Co, (192)Ir, or (169)Yb sources, some radiosensitive organs near the surface may be exposed to high absorbed doses. This may be reduced by covering the implants with a lead shield on the body surface, which results in dosimetric perturbations. Monte Carlo simulations in Geant4 were performed for the three radionuclides placed at a single dwell position. Four different shield thicknesses (0, 3, 6, and 10 mm) and three different source depths (0, 5, and 10 mm) in water were considered, with the lead shield placed at the phantom surface. Backscatter dose enhancement and transmission data were obtained for the lead shields. Results were corrected to account for a realistic clinical case with multiple dwell positions. The range of the high backscatter dose enhancement in water is 3 mm for (60)Co and 1 mm for both (192)Ir and (169)Yb. Transmission data for (60)Co and (192)Ir are smaller than those reported by Papagiannis et al (2008 Med. Phys. 35 4898-4906) for brachytherapy facility shielding; for (169)Yb, the difference is negligible. In conclusion, the backscatter overdose produced by the lead shield can be avoided by just adding a few millimetres of bolus. Transmission data provided in this work as a function of lead thickness can be used to estimate healthy organ equivalent dose saving. Use of a lead shield is justified.

Refining Treatment Planning in STereotactic Arrhythmia Radioablation: Benchmark Results and Consensus Statement From the STOPSTORM.eu Consortium.

PURPOSE: STereotactic Arrhythmia Radioablation (STAR) showed promising results in patients with refractory ventricular tachycardia. However, clinical data are scarce and heterogeneous. The STOPSTORM.eu consortium was established to investigate and harmonize STAR in Europe. The primary goal of this benchmark study was to investigate current treatment planning practice within the STOPSTORM project as a baseline for future harmonization. METHODS AND MATERIALS: Planning target volumes (PTVs) overlapping extracardiac organs-at-risk and/or cardiac substructures were generated for 3 STAR cases. Participating centers were asked to create single-fraction treatment plans with 25 Gy dose prescriptions based on in-house clinical practice. All treatment plans were reviewed by an expert panel and quantitative crowd knowledge-based analysis was performed with independent software using descriptive statistics for International Commission on Radiation Units and Measurements report 91 relevant parameters and crowd dose-volume histograms. Thereafter, treatment planning consensus statements were established using a dual-stage voting process. RESULTS: Twenty centers submitted 67 treatment plans for this study. In most plans (75%) intensity modulated arc therapy with 6 MV flattening filter free beams was used. Dose prescription was mainly based on PTV D95% (49%) or D96%-100% (19%). Many participants preferred to spare close extracardiac organs-at-risk (75%) and cardiac substructures (50%) by PTV coverage reduction. PTV D0.035cm3 ranged from 25.5 to 34.6 Gy, demonstrating a large variety of dose inhomogeneity. Estimated treatment times without motion compensation or setup ranged from 2 to 80 minutes. For the consensus statements, a strong agreement was reached for beam technique planning, dose calculation, prescription methods, and trade-offs between target and extracardiac critical structures. No agreement was reached on cardiac substructure dose limitations and on desired dose inhomogeneity in the target. CONCLUSIONS: This STOPSTORM multicenter treatment planning benchmark study not only showed strong agreement on several aspects of STAR treatment planning, but also revealed disagreement on others. To standardize and harmonize STAR in the future, consensus statements were established; however, clinical data are urgently needed for actionable guidelines for treatment planning.

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.

Dosimetry revisited for the HDR 192Ir brachytherapy source model mHDR-v2.

PURPOSE: Recently, the manufacturer of the HDR 192Ir mHDR-v2 brachytherapy source reported small design changes (referred to herein as mHDR-v2r) that are within the manufacturing tolerances but may alter the existing dosimetric data for this source. This study aimed to (1) check whether these changes affect the existing dosimetric data published for this source; (2) obtain new dosimetric data in close proximity to the source, including the contributions from 192Ir electrons and considering the absence of electronic equilibrium; and (3) obtain scatter dose components for collapsed cone treatment planning system implementation. METHODS: Three different Monte Carlo (MC) radiation transport codes were used: MCNP5, PENELOPE2008, and GEANT4. The source was centrally positioned in a 40 cm radius water phantom. Absorbed dose and collision kerma were obtained using 0.1 mm (0.5 mm) thick voxels to provide high-resolution dosimetry near (far from) the source. Dose-rate distributions obtained with the three MC codes were compared. RESULTS: Simulations of mHDR-v2 and mHDR-v2r designs performed with three radiation transport codes showed agreement typically within 0.2% for r > or = 0.25 cm. Dosimetric contributions from source electrons were significant for r < 0.25 cm. The dose-rate constant and radial dose function were similar to those from previous MC studies of the mHDR-v2 design. The 2D anisotropy function also coincided with that of the mHDR-v2 design for r > or = 0.25 cm. Detailed results of dose distributions and scatter components are presented for the modified source design. CONCLUSIONS: Comparison of these results to prior MC studies showed agreement typically within 0.5% for r > or = 0.25 cm. If dosimetric data for r < 0.25 cm are not needed, dosimetric results from the prior MC studies will be adequate.

Influence of photon energy spectra from brachytherapy sources on Monte Carlo simulations of kerma and dose rates in water and air.

PURPOSE: For a given radionuclide, there are several photon spectrum choices available to dosimetry investigators for simulating the radiation emissions from brachytherapy sources. This study examines the dosimetric influence of selecting the spectra for 192Ir, 125I, and 103Pd on the final estimations of kerma and dose. METHODS: For 192Ir, 125I, and 103Pd, the authors considered from two to five published spectra. Spherical sources approximating common brachytherapy sources were assessed. Kerma and dose results from GEANT4, MCNP5, and PENELOPE-2008 were compared for water and air. The dosimetric influence of 192Ir, 125I, and 103Pd spectral choice was determined. RESULTS: For the spectra considered, there were no statistically significant differences between kerma or dose results based on Monte Carlo code choice when using the same spectrum. Water-kerma differences of about 2%, 2%, and 0.7% were observed due to spectrum choice for 192Ir, 125, and 103Pd, respectively (independent of radial distance), when accounting for photon yield per Bq. Similar differences were observed for air-kerma rate. However, their ratio (as used in the dose-rate constant) did not significantly change when the various photon spectra were selected because the differences compensated each other when dividing dose rate by air-kerma strength. CONCLUSIONS: Given the standardization of radionuclide data available from the National Nuclear Data Center (NNDC) and the rigorous infrastructure for performing and maintaining the data set evaluations, NNDC spectra are suggested for brachytherapy simulations in medical physics applications.

Study of encapsulated 170Tm sources for their potential use in brachytherapy.

PURPOSE: High dose-rate (HDR) brachytherapy is currently performed with 192Ir sources, and 60Co has returned recently into clinical use as a source for this kind of cancer treatment. Both radionuclides have mean photon energies high enough to require specific shielded treatment rooms. In recent years, 169Yb has been explored as an alternative for HDR-brachytherapy implants. Although it has mean photon energy lower than 192Ir, it still requires extensive shielding to deliver treatment. An alternative radionuclide for brachytherapy is 170Tm (Z=69) because it has three physical properties adequate for clinical practice: (a) 128.6 day half-life, (b) high specific activity, and (c) mean photon energy of 66.39 keV. The main drawback of this radionuclide is the low photon yield (six photons per 100 electrons emitted). The purpose of this work is to study the dosimetric characteristics of this radionuclide for potential use in HDR-brachytherapy. METHODS: The authors have assumed a theoretical 170Tm cylindrical source encapsulated with stainless steel and typical dimensions taken from the currently available HDR 192Ir brachytherapy sources. The dose-rate distribution was calculated for this source using the GEANT4 Monte Carlo (MC) code considering both photon and electron 170Tm spectra. The AAPM TG-43 U1 brachytherapy dosimetry parameters were derived. To study general properties of 170Tm encapsulated sources, spherical sources encapsulated with stainless steel and platinum were also studied. Moreover, the influence of small variations in the active core and capsule dimensions on the dosimetric characteristics was assessed. Treatment times required for a 170Tm source were compared to those for 192Ir and 169Yb for the same contained activity. RESULTS: Due to the energetic beta spectrum and the large electron yield, the bremsstrahlung contribution to the dose was of the same order of magnitude as from the emitted gammas and characteristic x rays. Moreover, the electron spectrum contribution to the dose was significant up to 4 mm from the source center compared to the photon contribution. The dose-rate constant lamda of the cylindrical source was 1.23 cGy h(-1) U(-1). The behavior of the radial dose function showed promise for applications in brachytherapy. Due to the electron spectrum, the anisotropy was large for r <6 mm. Variations in manufacturing tolerances did not significantly influence the final dosimetry data when expressed in cGy h(-1) U(-1). For typical capsule dimensions, maximum reference dose rates of about 0.2, 10, and 2 Gy min(-1) would then be obtained for 170Tm, 192Ir, and 169Yb, respectively, resulting in treatment times greater than those for HDR 192Ir brachytherapy. CONCLUSIONS: The dosimetric characteristics of source designs exploiting the low photon energy of 170Tm were studied for potential application in HDR-brachytherapy. Dose-rate distributions were obtained for cylindrical and simplified spherical 170Tm source designs (stainless steel and platinum capsule materials) using MC calculations. Despite the high activity of 170Tm, calculated treatment times were much longer than for 192Ir.

Evaluation of high-energy brachytherapy source electronic disequilibrium and dose from emitted electrons.

PURPOSE: The region of electronic disequilibrium near photon-emitting brachytherapy sources of high-energy radionuclides (60Co, 137CS, 192Ir, and 169Yb) and contributions to total dose from emitted electrons were studied using the GEANT4 and PENELOPE Monte Carlo codes. METHODS: Hypothetical sources with active and capsule materials mimicking those of actual sources but with spherical shape were examined. Dose contributions due to source photons, x rays, and bremsstrahlung; source beta-, Auger electrons, and internal conversion electrons; and water collisional kerma were scored. To determine if conclusions obtained for electronic equilibrium conditions and electron dose contribution to total dose for the representative spherical sources could be applied to actual sources, the 192Ir mHDR-v2 source model (Nucletron B.V., Veenendaal, The Netherlands) was simulated for comparison to spherical source results and to published data. RESULTS: Electronic equilibrium within 1% is reached for 60Co, 137CS, 192Ir, and 169Yb at distances greater than 7, 3.5, 2, and 1 mm from the source center, respectively, in agreement with other published studies. At 1 mm from the source center, the electron contributions to total dose are 1.9% and 9.4% for 60Co and 192Ir, respectively. Electron emissions become important (i.e., > 0.5%) within 3.3 mm of 60Co and 1.7 mm of 192Ir sources, yet are negligible over all distances for 137Cs and 169Yb. Electronic equilibrium conditions along the transversal source axis for the mHDR-v2 source are comparable to those of the spherical sources while electron dose to total dose contribution are quite different. CONCLUSIONS: Electronic equilibrium conditions obtained for spherical sources could be generalized to actual sources while electron contribution to total dose depends strongly on source dimensions, material composition, and electron spectra.

An approach to using conventional brachytherapy software for clinical treatment planning of complex, Monte Carlo-based brachytherapy dose distributions.

Certain brachytherapy dose distributions, such as those for LDR prostate implants, are readily modeled by treatment planning systems (TPS) that use the superposition principle of individual seed dose distributions to calculate the total dose distribution. However, dose distributions for brachytherapy treatments using high-Z shields or having significant material heterogeneities are not currently well modeled using conventional TPS. The purpose of this study is to establish a new treatment planning technique (Tufts technique) that could be applied in some clinical situations where the conventional approach is not acceptable and dose distributions present cylindrical symmetry. Dose distributions from complex brachytherapy source configurations determined with Monte Carlo methods were used as input data. These source distributions included the 2 and 3 cm diameter Valencia skin applicators from Nucletron, 4-8 cm diameter AccuBoost peripheral breast brachytherapy applicators from Advanced Radiation Therapy, and a 16 mm COMS-based eye plaque using 103Pd, 125I, and 131Cs seeds. Radial dose functions and 2D anisotropy functions were obtained by positioning the coordinate system origin along the dose distribution cylindrical axis of symmetry. Origin:tissue distance and active length were chosen to minimize TPS interpolation errors. Dosimetry parameters were entered into the PINNACLE TPS, and dose distributions were subsequently calculated and compared to the original Monte Carlo-derived dose distributions. The new planning technique was able to reproduce brachytherapy dose distributions for all three applicator types, producing dosimetric agreement typically within 2% when compared with Monte Carlo-derived dose distributions. Agreement between Monte Carlo-derived and planned dose distributions improved as the spatial resolution of the fitted dosimetry parameters improved. For agreement within 5% throughout the clinical volume, spatial resolution of dosimetry parameter data < or = 0.1 cm was required, and the virtual brachytherapy source data set included over 5000 data points. On the other hand, the lack of consideration for applicator heterogeneity effect caused conventional dose overestimates exceeding an order of magnitude in regions of clinical interest. This approach is rationalized by the improved dose estimates. In conclusion, a new technique was developed to incorporate complex Monte Carlo-based brachytherapy dose distributions into conventional TPS. These results are generalizable to other brachytherapy source types and other TPS.

TG-43 U1 based dosimetric characterization of model 67-6520 Cs-137 brachytherapy source.

PURPOSE: Brachytherapy treatment has been a cornerstone for management of various cancer sites, particularly for the treatment of gynecological malignancies. In low dose rate brachytherapy treatments, 137Cs sources have been used for several decades. A new 137Cs source design has been introduced (model 67-6520, source B3-561) by Isotope Products Laboratories (IPL) for clinical application. The goal of the present work is to implement the TG-43 U1 protocol in the characterization of the aforementioned 137Cs source. METHODS: The dosimetric characteristics of the IPL 137Cs source are measured using LiF thermoluminescent dosimeters in a Solid Water phantom material and calculated using Monte Carlo simulations with the GEANT4 code in Solid Water and liquid water. The dose rate constant, radial dose function, and two-dimensional anisotropy function of this source model were obtained following the TG-43 U1 recommendations. In addition, the primary and scatter dose separation (PSS) formalism that could be used in convolution/superposition methods to calculate dose distributions around brachytherapy sources in heterogeneous media was studied. RESULTS: The measured and calculated dose rate constants of the IPL 137Cs source in Solid Water were found to be 0.930 (+/-7.3%) and 0.928 (+/-2.6%) cGy h(-1) U(-1), respectively. The agreement between these two methods was within our experimental uncertainties. The Monte Carlo calculated value in liquid water of the dose rate constant was null set=0.948 (+/-2.6%) cGy h(-1) U(-1). Similarly, the agreement between measured and calculated radial dose functions and the anisotropy functions was found to be within +/-5%. In addition, the tabulated data that are required to characterize the source using the PSS formalism were derived. CONCLUSIONS: In this article the complete dosimetry of the newly designed 137Cs IPL source following the AAPM TG-43 U1 dosimetric protocol and the PSS formalism is provided.

Equivalent phantom sizes and shapes for brachytherapy dosimetric studies of 192Ir and 137Cs.

The impact of phantom size and shape in brachytherapy dosimetry was assessed using Monte Carlo methods in liquid water for 192Ir and 137Cs point sources. This is needed since differences in published dosimetry data, both measurements and simulations, employ a variety of phantom sizes and shapes which can cause dose differences exceeding 30% near the phantom periphery. Spheres of radius, Rsph, 10-40 cm were examined to determine the equivalent spherical phantom size to a variety of cylinder and cube sizes, Rcyl and Rcube, respectively. These sizes ranged from 10 to 30 cm. The equivalent Rsph for a given size cylinder or cube was determined using a figure of merit (FOM) function to minimize differences between radial dose functions, g(r). Using the FOM approach, a linear fit (R2 > 0.99) was obtained for the equivalent Rsph for a given size cylinder or cube. The equivalent phantom for a cylinder, of 40 cm diameter and length 40 cm, is a sphere of 21 cm in radius and the equivalent phantom for a cube of 30 cm on each side is a sphere of 17.5 in radius. When normalizing all results to r=1 cm for g(r) comparisons of phantom shape, the absolute dose rates were equivalent within 0.1% for Rsph > or =10 cm for both 192Ir and 137Cs. Correlation factors to permit comparisons of unbounded g(r) data for r < or =15 cm in 20 published datasets resulted in agreement generally within 2%. Residual differences with four datasets were attributed to methodological uncertainties in the published references.

Development and clinical characterization of a novel 2041 liquid-filled ionization chambers array for high-resolution verification of radiotherapy treatments.

PURPOSE: The aim of this study was to present a novel 2041 liquid-filled ionization chamber array for high-resolution verification of radiotherapy treatments. MATERIALS AND METHODS: The prototype has 2041 ionization chambers of 2.5 × 2.5 mm2 area filled with isooctane. The detection elements are arranged in a central square grid of 43 × 43, totally covering an area of 107.5 × 107.5 mm2 . The central inline and cross-line are extended to 227 mm and the diagonals to 321 mm to be able to perform profile measurements of large fields. We have studied stability, pixel response uniformity, dose rate dependence, depth and field size dependence and anisotropy. We present results for output factors, tongue-and-groove, garden fence, small field profiles, irregular fields, and verification of dose planes of patient treatments. RESULTS: Comparison with other detectors used for small field dosimetry (SFD, CC13, microDiamond) has shown good agreement. Output factors measured with the device for square fields ranging from 10 × 10 to 100 × 100 mm2 showed relative differences within 1%. The response of the detector shows a strong dependence on the angle of incident radiation that needs to be corrected for. On the other hand, inter-pixel relative response variations in the 0.95-1.08 range have been found and corrected for. The application of the device for the verification of dose planes of several treatments has shown gamma passing rates above 97% for tolerances of 2% and 2 mm. The verification of other clinical fields, like small fields and irregular fields used in the commissioning of the TPS, also showed large passing rates. The verification of garden fence and tongue-and-groove fields was affected by volume-averaging effects. CONCLUSIONS: The results show that the liquid filled ionization chamber prototype here presented is appropriate for the verification of radiotherapy treatments with high spatial resolution. Recombination effects do not affect very much the verification of relative dose distributions. However, verification of absolute dose distributions may require normalization to a radiation field which is representative of the dose rate of the treatment delivered.

A generic TG-186 shielded applicator for commissioning model-based dose calculation algorithms for high-dose-rate 192 Ir brachytherapy.

PURPOSE: A 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) 192 Ir shielded applicator has been designed and benchmarked. METHODS: A generic HDR 192 Ir 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 192 Ir source located at one dwell position at its center was computed using two commercial TPSs incorporating MBDCAs (Oncentra® 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. RESULTS: The 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. CONCLUSIONS: The 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 192 Ir brachytherapy.

Broad-beam transmission data for new brachytherapy sources, Tm-170 and Yb-169.

The characteristics of the radionuclides (170)Tm and (169)Yb are highly interesting for their use as high dose-rate brachytherapy sources. The introduction of brachytherapy equipment containing these sources will lead to smaller required thicknesses of the materials used in radiation protection barriers compared with the use of conventional sources such as (192)Ir and (137)Cs. The purpose of this study is to determine the required thicknesses of protection material for the design of the protecting walls. Using the Monte Carlo method, transmission data were derived for broad-beam geometries through lead and concrete barriers, from which the first half value layer and tenth value layer are obtained. In addition, the dose reduction in a simulated patient was studied to determine whether transmission in the patient is a relevant factor in radiation protection calculations.

Commissioning and quality assurance procedures for the HDR Valencia skin applicators.

The Valencia applicators (Nucletron, an Elekta company, Elekta AB, Stockholm, Sweden) are cup-shaped tungsten applicators with a flattening filter used to collimate the radiation produced by a high-dose-rate (HDR) 192Ir source, and provide a homogeneous absorbed dose at a given depth. This beam quality provides a good option for the treatment of skin lesions at shallow depth (3-4 mm). The user must perform commissioning and periodic testing of these applicators to guarantee the proper and safe delivery of the intended absorbed dose, as recommended in the standards in radiation oncology. In this study, based on AAPM and GEC-ESTRO guidelines for brachytherapy units and our experience, a set of tests for the commissioning and periodic testing of the Valencia applicators is proposed. These include general considerations, verification of the manufacturer documentation and physical integrity, evaluation of the source-to-indexer distance and reproducibility, setting the library plan in the treatment planning system, evaluation of flatness and symmetry, absolute output and percentage depth dose verification, independent calculation of the treatment time, and visual inspection of the applicator before each treatment. For each test, the proposed methodology, equipment, frequency, expected results, and tolerance levels (when applicable) are provided.

A dosimetric study of Leipzig applicators.

PURPOSE: To obtain the absolute dose-rate distribution in liquid water for all six cup-shaped Leipzig applicators by means of an experimentally validated Monte Carlo (MC) code. These six applicators were used in high-dose-rate (HDR) afterloaders with the "classic" and v2 (192)Ir sources. The applicators have an inner diameter of 1, 2, and 3 cm, with the source traveling parallel or perpendicular to the contact surface. METHODS AND MATERIALS: The MC GEANT4 code was used to obtain the dose-rate distribution in liquid water for the six applicators and the two HDR source models. To normalize the applicator output factors, a MC simulation for the "classic" and v2 sources in air was performed to estimate the air-kerma strength. To validate this specific application and to guarantee that realistic source-applicator geometry was considered, an experimental verification procedure was implemented in this study, in accordance with the TG43U1 recommendations. Thermolumniscent dosimeter chips and a parallel plate ionization chamber in a polymethyl methacrylate (PMMA) phantom were used to verify the MC results for the six applicators in a microSelectronHDR afterloader with the "classic" source. Dose-rate distributions dependence on phantom size has been evaluated using two different phantom sizes. RESULTS: Percentage depth dose and off-axis profiles were obtained normalized at a depth of 3 mm along the central axis for both phantom sizes. A table of output factors, normalized to 1 U of source kerma strength at this depth, is presented. The dose measured in the PMMA phantom agrees within experimental uncertainties with the dose obtained by the MC GEANT4 code calculations. The phantom size influence on dose-rate distributions becomes significant at depths greater than 5 cm. CONCLUSIONS: MC-detailed simulation was performed for the Nucletron Leipzig HDR applicators. The matrix data obtained, with a grid separation of 0.5 mm, can be used to build a dataset in a convenient format to model these distributions for routine use with a brachytherapy treatment planning system.

Monte Carlo calculation of the TG-43 dosimetric parameters of a new BEBIG Ir-192 HDR source.

BACKGROUND AND PURPOSE: High dose rate (HDR) brachytherapy is a highly extended practice in clinical brachytherapy today. Quality dose rate distribution datasets of the HDR sources used in a clinical treatment are required. Because of the different source designs, a specific dosimetry dataset is required for each source model. In the recently published BRAPHYQS-ESTRO Report, an overview of available dosimetric data for all HDR Ir-192 sources is given, pointing out the lack of data for one of the sources that is used by the BEBIG MultiSource afterloading system (BEBIG GmbH, Germany). The purpose of this study is to obtain detailed dose rate distributions in liquid water media around this source. MATERIAL AND METHODS: The Monte Carlo code GEANT4 was used to estimate dose rate in water and air-kerma strength around the Ir-192 source. All the details of the stainless steel encapsulated BEBIG HDR 1.1mm in external diameter has been included in the simulation. RESULTS: A complete dosimetric dataset for the BEBIG Ir-192 HDR source is presented. TG43 dosimetric functions and parameters have been obtained as well as a 2D rectangular dose rate table, consistent with the TG43 dose calculation formalism. The dosimetric parameters and functions obtained for the BEBIG HDR source have been compared with that obtained in the literature for others HDR sources, showing that the use of specific datasets for this new source is justified. CONCLUSIONS: This dataset can be used as input in the TPS and to validate its calculations. As policy of BRAPHYQS-ESTRO task group, this dataset will be incorporated to the website: available to users in excel format.

A Monte Carlo study of intersource effects in dome-type applicators loaded with LDR Cs-137 sources.

In this study, the dose rate distributions produced by low dose rate Cs-137 sources loaded in afterloadable dome applicators are studied using the Monte Carlo method. Dose differences between Monte Carlo results and calculations done using the superposition principle are within 1-3% in front of the applicator and between 3 and 10% near and along the longitudinal source axis. Consequently, the real doses to lateral vaginal wall, rectum and bladder are very close to the doses estimated applying the superposition principle, while the dose to the vaginal cuff has been overestimated by up to 10%.

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.

Limitations of the TG-43 formalism for skin high-dose-rate brachytherapy dose calculations.

PURPOSE: In skin high-dose-rate (HDR) brachytherapy, sources are located outside, in contact with, or implanted at some depth below the skin surface. Most treatment planning systems use the TG-43 formalism, which is based on single-source dose superposition within an infinite water medium without accounting for the true geometry in which conditions for scattered radiation are altered by the presence of air. The purpose of this study is to evaluate the dosimetric limitations of the TG-43 formalism in HDR skin brachytherapy and the potential clinical impact. METHODS: Dose rate distributions of typical configurations used in skin brachytherapy were obtained: a 5 cm × 5 cm superficial mould; a source inside a catheter located at the skin surface with and without backscatter bolus; and a typical interstitial implant consisting of an HDR source in a catheter located at a depth of 0.5 cm. Commercially available HDR(60)Co and (192)Ir sources and a hypothetical (169)Yb source were considered. The Geant4 Monte Carlo radiation transport code was used to estimate dose rate distributions for the configurations considered. These results were then compared to those obtained with the TG-43 dose calculation formalism. In particular, the influence of adding bolus material over the implant was studied. RESULTS: For a 5 cm × 5 cm(192)Ir superficial mould and 0.5 cm prescription depth, dose differences in comparison to the TG-43 method were about -3%. When the source was positioned at the skin surface, dose differences were smaller than -1% for (60)Co and (192)Ir, yet -3% for (169)Yb. For the interstitial implant, dose differences at the skin surface were -7% for (60)Co, -0.6% for (192)Ir, and -2.5% for (169)Yb. CONCLUSIONS: This study indicates the following: (i) for the superficial mould, no bolus is needed; (ii) when the source is in contact with the skin surface, no bolus is needed for either (60)Co and (192)Ir. For lower energy radionuclides like (169)Yb, bolus may be needed; and (iii) for the interstitial case, at least a 0.1 cm bolus is advised for (60)Co to avoid underdosing superficial target layers. For (192)Ir and (169)Yb, no bolus is needed. For those cases where no bolus is needed, its use might be detrimental as the lack of radiation scatter may be beneficial to the patient, although the 2% tolerance for dose calculation accuracy recommended in the AAPM TG-56 report is not fulfilled.