Search

PURPOSE The purpose of this work is threefold: First, to obtain the phase space of an electronic brachytherapy (eBT) system designed for surface skin treatments. Second, to explore the use of some efficiency enhancing (EFEN) strategies in the determination of the phase space. Third, to use the phase space previously obtained to perform a dosimetric characterization of the Esteya eBT system. METHODS The Monte Carlo study of the 69.5 kVp x-ray beam of the Esteya® unit (Elekta Brachytherapy, Veenendaal, The Netherlands) was performed with PENELOPE2014. The EFEN strategies included the use of variance reduction techniques and mixed Class II simulations, where transport parameters were fine-tuned. Four source models were studied varying the most relevant parameters characterizing the electron beam impinging the target: the energy spectrum (mono-energetic or Gaussian shaped), and the electron distribution over the focal spot (uniform or Gaussian shaped). Phase spaces obtained were analyzed to detect differences in the calculated data due to the EFEN strategy or the source configuration. Depth dose curves and absorbed dose profiles were obtained for each source model and compared to experimental data previously published. RESULTS In our EFEN strategy, the interaction forcing variance reduction (VRIF) technique increases efficiency by a factor ~20. Tailoring the transport parameters values (C1 and C2) does not increase the efficiency in a significant way. Applying a universal cutoff energy EABS of 10 keV saves 84% of CPU time while showing negligible impact on the calculated results. Disabling the electron transport by imposing an electron energy cutoff of 70 keV (except for the target) saves an extra 8% (losing in the process 1.2% of the photons). The Gaussian energy source (FWHM = 10%, centered at the nominal kVp, homogeneous electron distribution) shows characteristic K-lines in its energy spectrum, not observed experimentally. The average photon energy using an ideal source (mono-energetic, homogeneous electron distribution) was 36.19 ± 0.09 keV, in agreement with the published measured data of 36.2 ± 0.2 keV. The use of a Gaussian-distributed electron source (mono-energetic) increases the penumbra by 50%, which is closer to the measurement results. The maximum discrepancy of the calculated percent depth dose with the corresponding measured values is 4.5% (at the phantom surface, less than 2% beyond 1 mm depth) and 5% (for the 80% of the field) in the dose profile. Our results agree with the findings published by other authors and are consistent within the expected Type A and B uncertainties. CONCLUSIONS Our results agree with the published measurement results within the reported uncertainties. The observed differences in PDD, dose profiles, and photon spectrum come from three main sources of uncertainty: intermachine variations, measurements, and Monte Carlo calculations. It has been observed that a mono-energetic source with a Gaussian electron distribution over the focal spot is a suitable choice to reproduce the experimental data.

Purpose: High dose rate (HDR) brachytherapy for treatment of small skin lesions using the Leipzig and Valencia applicators is a widely used technique. These applicators are equipped with an attachable plastic cap to be placed during fraction delivery to ensure electronic equilibrium and to prevent secondary electrons from reaching the skin surface. The purpose of this study is to report on the dosimetric impact of the cap being absent during HDR fraction delivery, which has not been explored previously in the literature. Methods: geant4 Monte Carlo simulations (version 10.0) have been performed for the Leipzig and Valencia applicators with and without the plastic cap. In order to validate the Monte Carlo simulations, experimental measurements using radiochromic films have been done. Results: Dose absorbed within 1 mm of the skin surface increases by a factor of 1500% for the Leipzig applicators and of 180% for the Valencia applicators. Deeper than 1 mm, the overdosage flattens up to a 10% increase. Conclusions: Differences of treating with or without the plastic cap are significant. Users must check always that the plastic cap is in place before any treatment in order to avoid overdosage of the skin. Prior to skin HDR fraction delivery, the timeout checklist should include verification of the cap placement.

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 192Ir 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 penelope2008MC 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.

The purpose of this work is threefold: First, to obtain the phase space of an electronic brachytherapy (eBT) system designed for surface skin treatments. Second, to explore the use of some efficiency enhancing (EFEN) strategies in the determination of the phase space. Third, to use the phase space previously obtained to perform a dosimetric characterization of the Esteya eBT system.The Monte Carlo study of the 69.5 kVp x-ray beam of the EsteyaIn our EFEN strategy, the interaction forcing variance reduction (VRIF) technique increases efficiency by a factor ~20. Tailoring the transport parameters values (C1 and C2) does not increase the efficiency in a significant way. Applying a universal cutoff energy EABS of 10 keV saves 84% of CPU time while showing negligible impact on the calculated results. Disabling the electron transport by imposing an electron energy cutoff of 70 keV (except for the target) saves an extra 8% (losing in the process 1.2% of the photons). The Gaussian energy source (FWHM = 10%, centered at the nominal kVp, homogeneous electron distribution) shows characteristic K-lines in its energy spectrum, not observed experimentally. The average photon energy using an ideal source (mono-energetic, homogeneous electron distribution) was 36.19 ± 0.09 keV, in agreement with the published measured data of 36.2 ± 0.2 keV. The use of a Gaussian-distributed electron source (mono-energetic) increases the penumbra by 50%, which is closer to the measurement results. The maximum discrepancy of the calculated percent depth dose with the corresponding measured values is 4.5% (at the phantom surface, less than 2% beyond 1 mm depth) and 5% (for the 80% of the field) in the dose profile. Our results agree with the findings published by other authors and are consistent within the expected Type A and B uncertainties.Our results agree with the published measurement results within the reported uncertainties. The observed differences in PDD, dose profiles, and photon spectrum come from three main sources of uncertainty: intermachine variations, measurements, and Monte Carlo calculations. It has been observed that a mono-energetic source with a Gaussian electron distribution over the focal spot is a suitable choice to reproduce the experimental data.

Recently, a new HDR 60Co brachytherapy source, Flexisource Co-60, has been developed (Nucletron B.V.). This study aims to obtain quality dosimetric data for this source for its use in clinical practice as required by AAPM and ESTRO.Penelope2008 and GEANT4 Monte Carlo codes were used to dosimetrically characterize this source. Water composition and mass density was that recommended by AAPM. Due to the high energy of the 60Co, dose for small distances cannot be approximated by collisional kerma. Therefore, we have considered absorbed dose to water for r0.75 cm and collisional kerma from 0.75r20 cm. To provide adequate spatial resolution, cells were 0.01 cm in thickness for r2 cm from the source and a factor of 10 thicker for 2r20 cm respectively. Angular sampling was taken every 2°. Additional simulations were performed to obtain SK as recommended by AAPM. Mass-energy absorption coefficients in water and air were consistently derived and used to calculate collisional kerma. Along-away tables and TG-43 formalism parameters and functions were derived. Dosimetric data were also provided following the primary and scatter dose separation for the collapsed cone technique.TG-43 dosimetry parameters with L = 0.35 cm were obtained. Results performed with both radiation transport codes showed agreement typically within 0.2% for r0.8 cm and up to 2% closer to the source. Using Penelope2008 and GEANT4, an average of Î= 1.085±0.003 cGy/(h U) (with k = 1, Type A uncertainties) was obtained. Dose rate constant, radial dose function and anisotropy functions for the Flexisource Co-60 are compared with published data for other Co-60 sources.Dosimetric data are provided for the new Flexisource Co-60 source not studied previously in the literature. Using the data provided by this study in the treatment planning systems, it can be used in clinical practice. This project has been funded by Nucletron BV.

Purpose: Recommendations of the American Association of Physicists in Medicine (AAPM) and the European Society for Radiotherapy and Oncology (ESTRO) on dose calculations for high-energy (average energy higher than 50 keV) photon-emitting brachytherapy sources are presented, including the physical characteristics of specific192Ir, 137Cs, and 60Co source models. Methods: This report has been prepared by the High Energy Brachytherapy Source Dosimetry (HEBD) Working Group. This report includes considerations in the application of the TG-43U1 formalism to high-energy photon-emitting sources with particular attention to phantom size effects, interpolation accuracy dependence on dose calculation grid size, and dosimetry parameter dependence on source active length. Results: Consensus datasets for commercially available high-energy photon sources are provided, along with recommended methods for evaluating these datasets. Recommendations on dosimetry characterization methods, mainly using experimental procedures and Monte Carlo, are established and discussed. Also included are methodological recommendations on detector choice, detector energy response characterization and phantom materials, and measurement specification methodology. Uncertainty analyses are discussed and recommendations for high-energy sources without consensus datasets are given. Conclusions: Recommended consensus datasets for high-energy sources have been derived for sources that were commercially available as of January 2010. Data are presented according to the AAPM TG-43U1 formalism, with modified interpolation and extrapolation techniques of the AAPM TG-43U1S1 report for the 2D anisotropy function and radial dose function.

Purpose: For a given radionuclide, there are several photonspectrum choices available to dosimetry investigators for simulating the radiation emissions from brachytherapy sources. This study examines the dosimetric influence of selecting the spectra for I 192 r , I 125 , and P 103 d on the final estimations of kerma and dose. Methods: For I 192 r , I 125 , and P 103 d , 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 I 192 r , I 125 , and P 103 d 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 I 192 r , I 125 , and P 103 d , 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 photonspectra 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.

Purpose: The region of electronic disequilibrium near photon-emitting brachytherapysources of high-energy radionuclides ( C 60 o , C 137 s , I 192 r , and Y 169 b ) and contributions to total dose from emitted electrons were studied using the GEANT4 and PENELOPEMonte 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 sourcephotons, x rays, and bremsstrahlung; source β − , Auger electrons, and internal conversionelectrons; and water collisional kerma were scored. To determine if conclusions obtained for electronic equilibrium conditions and electrondose contribution to total dose for the representative spherical sources could be applied to actual sources, the I 192 r 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 C 60 o , C 137 s , I 192 r , and Y 169 b 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 C 60 o and I 192 r , respectively. Electron emissions become important (i.e., > 0.5 % ) within 3.3 mm of C 60 o and 1.7 mm of I 192 r sources, yet are negligible over all distances for C 137 s and Y 169 b . Electronic equilibrium conditions along the transversal source axis for the mHDR-v2 source are comparable to those of the spherical sources while electrondose 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.

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.

To address the limited availability of radiation shielding data for brachytherapy as well as some disparity in existing data, Monte Carlo simulation was used to generate radiation transmission data for 60Co, 137CS, 198Au, 192Ir 169Yb, 170Tm, 131Cs, 125I, and 103pd photons through concrete, stainless steel, lead, as well as lead glass and baryte concrete. Results accounting for the oblique incidence of radiation to the barrier, spectral variation with barrier thickness, and broad beam conditions in a realistic geometry are compared to corresponding data in the literature in terms of the half value layer (HVL) and tenth value layer (TVL) indices. It is also shown that radiation shielding calculations using HVL or TVL values could overestimate or underestimate the barrier thickness required to achieve a certain reduction in radiation transmission. This questions the use of HVL or TVL indices instead of the actual transmission data. Therefore, a three-parameter model is fitted to results of this work to facilitate accurate and simple radiation shielding calculations.

The two most commonly used modalities of cancer treatment in clinical brachytherapy practice today are high dose rate (HDR) and pulsed dose rate (PDR) brachytherapy. In a clinical treatment, quality dose rate distribution data sets of the brachytherapy sources are required for each source model. The purpose of this study is to obtain detailed dose rate distributions around the new BEBIG HDR and PDR Ir-192 brachytherapy sources. These distributions will then be used as input data in the treatment planning systems dedicated to brachytherapy and its calculations can be verified. The Monte Carlo method was used to obtain the dose rate distributions around the sources studied, taking into account the AAPM-ESTRO recent recommendations. A complete dosimetric data set for the BEBIG Ir-192 HDR and PDR sources, types Ir2.A85-2 and Ir2.A85-1, were obtained. This dosimetric data set is composed of the TG-43 dosimetric functions and parameters and along-away dose rate table to facilitate quality control of treatment planning systems.

Ir-192 wires have been used in low-dose-rate brachytherapy for many years. Commercially available treatment planning systems approximate the dose rate distribution of the straight or curved wires applying the superposition principle using one of the following methods: (i) The wire is modeled as a set of point sources, (ii) the wire is modeled as a set of small straight segment wires, (iii) the values of the parameters and functions of the American Association of Physicists in Medicine (AAPM) Task Group 43 protocol are obtained for wire lengths between 3 and 7 cm assuming some simplifications. The dose rate distributions obtained using these methods for linear wires of different lengths and U-shaped wires present significant deviations compared to those obtained by Monte Carlo. In the present study we propose a new method to model 192Ir wires of any length and shape, named the Two Lengths based Segmented method. This method uses the formalism stated in the AAPM Task Group 43 protocol for two straight wires only, 0.5 and 1 cm, to obtain the dose rate distribution around wires of any length (down to 0.3 cm and up to 10 cm) improving on the results of the aforementioned ones. This method can easily be applied to dose calculations around other wires, such as Pd-103 ones.

The purposes of this study are: (i) to design field flattening filters for the Leipzig applicators of 2 and 3 cm of inner diameter with the source traveling parallel to the applicator contact surface, which are accessories of the microSelectron-HDR afterloader (Nucletron, Veenendaal, The Netherlands). These filters, made of tungsten, aim to flatten the heterogeneous dose distribution obtained with the Leipzig applicators. (ii) To estimate the dose rate distributions for these Leipzig+filter applicators by means of the Monte Carlo (MC) method. (iii) To experimentally verify these distributions for prototypes of these new applicators, and (iv) to obtain the correspondence factors to measure the output of the applicators by the user using an insert into a well chamber. The MC GEANT4 code has been used to design the filters and to obtain the dose rate distributions in liquid water for the two Leipzig+filter applicators. In order to validate this specific application and to guarantee that realistic source-applicator geometry has been considered, an experimental verification procedure was implemented in this study, in accordance with the updated recommendations of the American Association of Physicists in Medicine Task Group No. 43 U1 Report. Thermoluminescent dosimeters, radiochromic film, and a pin-point ionization chamber in a plastic [polymethylmethacrylate (PMMA)] phantom were used to verify the MC results for the two applicators of a microSelectron-HDR afterloader with the mHDR-v2 source. To verify the output of the Leipzig +filter applicators, correspondence factors were deduced for the well chambers HDR100-plus (Standard Imaging, Inc., Middleton, WI) and TM33004 (PTW, Freiburg, Germany) using a specific insert for both applicators. The doses measured in the PMMA phantom agree within experimental uncertainties with the dose obtained by the MC calculations. Percentage depth dose and off-axis profiles were obtained normalized at a depth of 3 mm along the central applicator axis in a cylindrical 20 x 20 cm water phantom. A table of output factors, normalized to 1 U of source air kerma strength at this depth, is presented. Correspondence factors were obtained for the two well chambers considered. The matrix data obtained in the MC simulation with a grid separation of 0.5 mm has been used to build a data set in a convenient format to model these distributions for routine use with a brachytherapy treatment planning system.

In this work, the dose rate distribution of a new Ir-192 high dose rate source (Flexisource used in the afterloading Flexitron system, Isodose Control, Veenendaal, The Netherlands) is studied by means of Monte Carlo techniques using the GEANT4 code. The dosimetric parameters of the Task Group No. 43 Report (TG43) formalism and two-dimensional rectangular look-up tables have been obtained.

The H-type Leipzig applicators are accessories of the microSelectron-HDR system (Nucletron, Veenendaal, The Netherlands) for treatment of superficial malignancies. Recently, the dose rate distributions in liquid water for the whole set of applicators using both source models available for the microSelectron-HDR afterloaders have been obtained by means of the experimentally validated Monte Carlo (MC) code GEANT4. Also an output table (cGy/hU) at 3 mm depth on the applicator central axis was provided. The output verification of these applicators by the user, prior to their clinical use, present practical problems: small detectors such as thermoluminescent dosimeters or parallel-plate ionization chambers are not easily used for verification in a clinical environment as they require a rigid setup with the Leipzig applicator and a phantom. In contrast, well-type ionization chambers are readily available in radiotherapy departments. This study presents a technique based on the HDR1000Plus well chamber (Standar Imaging) measurements with a special insert, which allows the output verification of the H-type Leipzig applicators on a routine basis. This technique defines correspondence factors (CF) between the in water dose rate output of the Leipzig applicators (cGy/hU) obtained with MC and the reading on the well chamber with the special insert, normalized to themore » HDR calibration factor with the HDR insert and to the source strength. To commission the applicators (with the well chamber and the special insert used), the physicist should check if the CF value agrees with its tabulated values presented in this work. If the differences are within 5% the tabulated output values can be used in clinical dosimetry. This technique allows the output validation of the Leipzig applicators with a well chamber widely used for HDR Ir-192 source strength measurements. It can easily be adapted to other types of well chambers for HDR source output verification.« less

Purpose: A failure mode and effect analysis (FMEA) of skin lesions treatment process using Esteya™ device (Elekta Brachyterapy, Veenendaal, The Netherlands) was performed, with the aim of increasing the quality of the treatment and reducing the likelihood of unwanted events. Methods: A multidisciplinary team with experience in the treatment process met to establish the process map, which outlines the flow of various stages for such patients undergoing skin treatment. Potential failure modes (FM) were identified and the value of severity (S), frequency of occurrence (O), and lack of detectability (D) of the proposed FM were scored individually, each on a scale of 1 to 10 following TG-100 guidelines of the AAPM. These failure modes were ranked according to our risk priority number (RPN) and S scores. The efficiency of existing quality management tools was analyzed through a reassessment of the O and D made by consensus. Results: 149 FM were identified, 43 of which had RPN ≥ 100 and 30 had S ≥ 7. After introduction of the tools of quality management, only 3 FM had RPN ≥ 100 and 22 FM had RPN ≥ 50. These 22 FM were thoroughly analyzed and new tools for quality management were proposed. The most common cause of highest RPN FM was associated with the heavy patient workload and the continuous and accurate applicator-patient skin contact during the treatment. To overcome this second item, a regular quality control and setup review by a second individual before each treatment session was proposed. Conclusion: FMEA revealed some of the FM potentials that were not predicted during the initial implementation of the quality management tools. This exercise was useful in identifying the need of periodic update of the FMEA process as new potential failures can be identified.

Functional fits for the anisotrophy function and the radial dose function, have been studied, in a previous work, in order to characterize dose-rate distributions around some of the high-intensity 192 Ir sources. The purpose of the present work is to complete the previous one in order to include all the existing HDR and PDR 192 Ir sources. The sources addressed here are: the Buchler source from Amersham, the 12i and Plus PDR sources and the 12i and Plus HDR sources from GammaMed, and the new VariSource HDR source wire model VS2000 from Varian Oncology Systems.

The low dose rate CDC-type miniature cylindrical 137 Cs sources are available, with one or three active beads, for use in source trains in automatic and manual afterloading systems for gynecological brachytherapy. Absolute dose rate distributions in water have been calculated around these sources using the Monte CarloGEANT3 code and they are presented as conventional two-dimensional Cartesian lookup tables. The AAPM Task Group 43 formalism for dose calculation has been also applied. The dose rate constant obtained for the one bead source is Λ=1.113±0.003 cGyh −1 U −1 , and the value for the three bead source is Λ=1.103±0.003 cGyh −1 U −1 . Finally, for the treatment planning systems based on Sievert-type algorithms, the attenuation coefficients that best reproduce Monte Carlodose rate distribution are given.

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.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.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.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.

The purpose of this study is to find fitted functional forms to the anisotropy function, F(r,θ), and the radial dose function, g(r), in order to characterize dose-rate distributions around all the high-intensity 192 Ir sources currently in use. Dosimetry data are at present available as tables for: the microSelectron HDR (“classic” and “new” design models), the PDR source, and the VariSource HDR source, expressed in terms of the AAPM Task Group No. 43 recommendations. There is only one paper out which introduces a functional form to fit the anisotropy function, but only for symmetric sources with respect to the transverse axis. However, dosimetric data of the HDR and PDR sources mentioned above cannot be reproduced with these functional forms. In our study F(r,θ) and g(r) published data are fitted with functional forms in such a way that appropriate limits are reached for both functions and the maximum fit error approaches the data uncertainty. The average fit error is less than 1% in all cases. These functional forms make handling data easier within the treatment planning system, avoiding the use of tabulated data.

Air-communicating well ionization chambers are commonly used to assess air kerma strength of sources used in brachytherapy. The signal produced is supposed to be proportional to the air density within the chamber and, therefore, a density-independent air kerma strength is obtained when the measurement is corrected to standard atmospheric conditions using the usual temperature and pressure correction factor. Nevertheless, when assessing low energy sources, the ionization chambers may not fulfill that condition and a residual density dependence still remains after correction. In this work, the authors examined the behavior of the PTW 34051 SourceCheck ionization chamber when measuring the air kerma strength of (125)I seeds.Four different SourceCheck chambers were analyzed. With each one of them, two series of measurements of the air kerma strength for (125)I selectSeed(TM) brachytherapy sources were performed inside a pressure chamber and varying the pressure in a range from 747 to 1040 hPa (560 to 780 mm Hg). The temperature and relative humidity were kept basically constant. An analogous experiment was performed by taking measurements at different altitudes above sea level.Contrary to other well-known ionization chambers, like the HDR1000 PLUS, in which the temperature-pressure correction factor overcorrects the measurements, in the SourceCheck ionization chamber they are undercorrected. At a typical atmospheric situation of 933 hPa (700 mm Hg) and 20 °C, this undercorrection turns out to be 1.5%. Corrected measurements show a residual linear dependence on the density and, as a consequence, an additional density dependent correction must be applied. The slope of this residual linear density dependence is different for each SourceCheck chamber investigated. The results obtained by taking measurements at different altitudes are compatible with those obtained with the pressure chamber.Variations of the altitude and changes in the weather conditions may produce significant density corrections, and that effect should be taken into account. This effect is chamber-dependent, indicating that a specific calibration is necessary for each particular chamber. To our knowledge, this correction has not been considered so far for SourceCheck ionization chambers, but its magnitude cannot be neglected in clinical practice. The atmospheric pressure and temperature at which the chamber was calibrated need to be taken into account, and they should be reported in the calibration certificate. In addition, each institution should analyze the particular response of its SourceCheck ionization chamber and compute the adequate correction factors. In the absence of a suitable pressure chamber, a possibility for this assessment is to take measurements at different altitudes, spanning a wide enough air density range.

The purpose of this study is to calculate basic dosimetry data for a CSM11 low dose rate 137Cs source in water. This source is widely used in afterloadable dome cylinders designed to homogeneously irradiate the vaginal cuff alone or additional areas of the vagina in hysterectomized patients. In this study, the Monte Carlo simulation code GEANT, incorporating in detail source geometry, is used to investigate the dosimetric characteristics of the source. The calculated data were analyzed using a fitting procedure that is described in detail. Absolute dose rate distributions in water were calculated around this source and are presented as conventional 2D Cartesian lookup tables (classically along-away tables). Also, the dose calculation formalism endorsed by the Interstitial Collaborative Working Group and the AAPM Task Group 43 have been calculated. The calculated dose rate constant for this source is lambda = 1.096 +/- 0.002 cGy h(-1) U(-1). The anisotropy function results in about 50% deviations from isotropy at positions near the long axis of the source. The radial dose function is given as a polynomial that reproduces the calculated data up to 20 cm. Best-fit values of attenuation coefficients suitable for use in Sievert integral calculations have been derived.

Monte Carlo simulations of absolute dose rate in liquid water are presented in the form of away-along tables for 1 and 5 cm 192 Ir wires of 0.5 and 0.6 mm diameter. Simulated absolute dose rate values can be used as benchmark data to verify the calculation results of treatment planning systems or directly as input data for treatment planning. Best fit value of an attenuation coefficient suitable for use in Sievert integral-type calculations has been derived based on Monte Carlo simulation results. For the treatment planning systems that are based on the TG43 formalism we have also computed the required dosimetry parameters.

The aim of this study was to obtain equivalent doses in radiosensitive organs (aside from the bladder and rectum) when applying high-dose-rate (HDR) brachytherapy to a localized prostate carcinoma using (60)Co or (192)Ir sources. These data are compared with results in a water phantom and with expected values in an infinite water medium. A comparison with reported values from proton therapy and intensity-modulated radiation therapy (IMRT) is also provided.Monte Carlo simulations in Geant4 were performed using a voxelized phantom described in International Commission on Radiological Protection (ICRP) Publication 110, which reproduces masses and shapes from an adult reference man defined in ICRP Publication 89. Point sources of (60)Co or (192)Ir with photon energy spectra corresponding to those exiting their capsules were placed in the center of the prostate, and equivalent doses per clinical absorbed dose in this target organ were obtained in several radiosensitive organs. Values were corrected to account for clinical circumstances with the source located at various positions with differing dwell times throughout the prostate. This was repeated for a homogeneous water phantom.For the nearest organs considered (bladder, rectum, testes, small intestine, and colon), equivalent doses given by (60)Co source were smaller (8%-19%) than from (192)Ir. However, as the distance increases, the more penetrating gamma rays produced by (60)Co deliver higher organ equivalent doses. The overall result is that effective dose per clinical absorbed dose from a (60)Co source (11.1 mSv/Gy) is lower than from a (192)Ir source (13.2 mSv/Gy). On the other hand, equivalent doses were the same in the tissue and the homogeneous water phantom for those soft tissues closer to the prostate than about 30 cm. As the distance increased, the differences of photoelectric effect in water and soft tissue, and appearance of other materials such as air, bone, or lungs, produced variations between both phantoms which were at most 35% in the considered organ equivalent doses. Finally, effective doses per clinical absorbed dose from IMRT and proton therapy were comparable to those from both brachytherapy sources, with brachytherapy being advantageous over external beam radiation therapy for the furthest organs.A database of organ equivalent doses when applying HDR brachytherapy to the prostate with either (60)Co or (192)Ir is provided. According to physical considerations, (192)Ir is dosimetrically advantageous over (60)Co sources at large distances, but not in the closest organs. Damage to distant healthy organs per clinical absorbed dose is lower with brachytherapy than with IMRT or protons, although the overall effective dose per Gy given to the prostate seems very similar. Given that there are several possible fractionation schemes, which result in different total amounts of therapeutic absorbed dose, advantage of a radiation treatment (according to equivalent dose to healthy organs) is treatment and facility dependent.

Purpose: To estimate the sensitivity of TrueBeam 2.0 Imaging System 6DoF automatic matching tool through the acquisition of cone-beam CT images in different phantoms applying submillimeter translations and rotations of tenths of a degree and registered with image simulation CT. Methods: To evaluate overall system-wide image, we consider two uncertainties source; First, the uncertainty of the manual phantom displacement (e_m). This uncertainty is calculated by a digital caliper (0.01 mm) for vertical (Vrt), lateral (Lat) and longitudinal (Lng). A digital inclinometer (0.01°) for the pitch and roll and the own phantom scale to evaluate the coordinate rotation (Rtn). The second uncertainty is the displacement detected by the algorithm system of matching (σ_d) that we obtain from the standard deviations of the different measurements. We use three different phantoms. The BrainLab Radiosurgery system for supporting masks with an anthropomorphic dummy adapted to allow displacements of 0.1 mm in Vrt, Lat and Lng dimensions and rotations of 0.1° in Pitch dimension. For the analysis of the Rtn and Roll dimensions we use two homemade phantoms (RinoRot and RinoRoll, La Fe Hospital, Valencia, Spain) that allow rotations of 0.3°. Results: In the case of manual displacement of 0.10 ± 0.03 mm in the translations, the system detect 0.10 ± 0.07 mm, 0.12 ± 0.07 mm and 0.13 ± 0.07 mm (mean ± SD) in Lat, Vrt and Lng respectively. In the case of rotational dimension, manual displacement of 0.3 ± 0.1° was detected with 0.19 ± 0.06°, 0.29 ± 0.03° and 0.27 ± 0.06° in Pitch, Roll and Rtn. Conclusion: We conclude that the sensitivity of the automatic matching system is within 0.10 mm in translations and 0.3° in rotations. These values are under the own sensitivity of the software.

Purpose: Utilization of HDR brachytherapy treatment of skin lesions using collimated applicators, such as the Valencia or Leipzig is increasing. These applicators are made of cup-shaped tungsten material in order to focalize the radiation into the lesion and to protect nearby tissues. These applicators have an attachable plastic cap that removes secondary electrons generated in the applicator and flattens the treatment surface. The purpose of this study is to examine the dosimetric impact of this cap, and the effect if the cap is not placed during the HDR fraction delivery. Methods: Monte Carlo simulations have been done using the code Geant4 for the Valencia and Leipzig applicators. Dose rate distributions have been obtained for the applicators with and without the plastic cap. An experimental study using EBT3 radiochromic film has been realized in order to verify the Monte Carlo results. Results: The Monte Carlo simulations show that absorbed dose in the first millimeter of skin can increase up to 180% for the Valencia applicator if the plastic cap is absent and up to 1500% for the Leipzig applicators. At deeper distances the increase of dose is smaller being about 10–15%. Conclusion: Important differences have been found if the plastic cap of the applicators is absent in the treatment producing an overdosage in the skin. The user should have a checklist to remind him check always before HDR fraction delivery to insure the plastic cap is placed on the applicator. This work was supported in part by Generalitat Valenciana under Project PROMETEOII/2013/010, by the Spanish Government under Project No. FIS2013-42156, and by a research agreement with Elekta Brachytherapy, Veenendaal, The Netherlands.

Purpose: The purpose of this study is to compare the transit dose due to the movement of high dose rate (HDR) Ir-192 and Co-60 sources along the transfer tube. This is performed by evaluating air-kerma differences in the vicinity of the transfer tube when both sources are moved with the same velocity from a HDR brachytherapy afterloader into a patient. Methods: Monte Carlo simulations have been performed using PENELOPE2014. mHDR-v2 and Flexisource sources have been considered. Collisional kerma has been scored. The sources were simulated within a plastic catheter located in an infinite air phantom. The movement of the seed was included by displacing their positions along the connecting catheter from z=−75 cm to z=+75 cm and combining them. Backscatter from the afterloader and the patient was not considered. Since modern afterloaders like Flexitron (Elekta) or Saginova (Bebig) are able to use equally Ir-192 and Co-60 sources it was assumed that both sources are displaced with equal speed. Typical content activity values were provided by the manufacturer (460 GBq for Ir-192 and 75 GBq for Co-60). Results: 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 activities. It was found that it varies slowly with distance (less than 10% variation) but strongly in time due to the shorter half-life of the Ir-192 (73.83 days). The maximum ratio is located close to the catheter with a value of 0.57 when both sources are installed by the manufacturer, while increasing up to 1.25 at the end of the recommended working life (90 days) of the Ir-192 source. Conclusion: Air-kerma ratios are almost constant (0.51–0.57) in the vicinity of the source. Nevertheless, air-kerma ratios increase rapidly whenever the Ir-192 is approaching the end of its life.

Purpose: Total Skin Electron Irradiation (TSEI) is a radiotherapy treatment which involves irradiating the entire body surface as homogeneously as possible. It is composed of an extensive multi-step technique in which quality management requires high consumption of resources and a fluid communication between the involved staff, necessary to improve the safety of treatment. The TG-100 proposes a new perspective of quality management in radiotherapy, presenting a systematic method of risk analysis throughout the global flow of the stages through the patient. The purpose of this work has been to apply TG-100 approach to the TSEI procedure in our institution. Methods: A multidisciplinary team specifically targeting TSEI procedure was formed, that met regularly and jointly developed the process map (PM), following TG-100 guidelines of the AAPM. This PM is a visual representation of the temporal flow of steps through the patient since start until the end of his stay in the radiotherapy service. Results: This is the first stage of the full risk analysis, which is being carried out in the center. The PM provides an overview of the process and facilitates the understanding of the team members who will participate in the subsequent analysis. Currently, the team is implementing the analysis of failure modes and effects (FMEA). The failure modes of each of the steps have been identified and assessors are assigning a value of severity (S), frequency of occurrence (O) and lack of detection (D) individually. To our knowledge, this is the first PM made for the TSEI. The developed PM can be useful for those centers that intend to implement the TSEI technique. Conclusion: The PM of TSEI technique has been established, as the first stage of full risk analysis, performed in a reference center in this treatment.

Purpose: Tissue heterogeneities and calcifications have significant impact on the dosimetry of low energy brachytherapy (BT). RayStretch is an analytical algorithm developed in our institution to incorporate heterogeneity corrections in LDR prostate brachytherapy. The aim of this work is to study its application in clinical cases by comparing its predictions with the results obtained with TG-43 and Monte Carlo (MC) simulations. Methods: A clinical implant (71 I-125 seeds, 15 needles) from a real patient was considered. On this patient, different volumes with calcifications were considered. Its properties were evaluated in three ways by i) the Treatment planning system (TPS) (TG-43), ii) a MC study using the Penelope2009 code, and iii) RayStretch. To analyse the performance of RayStretch, calcifications located in the prostate lobules covering 11% of the total prostate volume and larger calcifications located in the lobules and underneath the urethra for a total occupied volume of 30% were considered. Three mass densities (1.05, 1.20, and 1.35 g/cm3) were explored for the calcifications. Therefore, 6 different scenarios ranging from small low density calcifications to large high density ones have been discussed. Results: DVH and D90 results given by RayStretch agree within 1% with the full MC simulations. Although no effort has been done to improve RayStretch numerical performance, its present implementation is able to evaluate a clinical implant in a few seconds to the same level of accuracy as a detailed MC calculation. Conclusion: RayStretch is a robust method for heterogeneity corrections in prostate BT supported on TG-43 data. Its compatibility with commercial TPSs and its high calculation speed makes it feasible for use in clinical settings for improving treatment quality. It will allow in a second phase of this project, its use during intraoperative ultrasound planning. This study was partly supported by a fellowship grant from the Spanish Ministry of Education, by the Generalitat Valenciana under Project PROMETEOII/2013/010, by the Spanish Government under Project No. FIS2013-42156 and by the European Commission within the SeventhFramework Program through ENTERVISION (grant agreement number 264552).

Purpose: ABS and GEC-ESTRO have recommended MRI T2 for image guided brachytherapy. Recently, a new applicator (Benidorm Template, TB) has been developed in our Department (Rodriguez et al 2015). TB is fully MRI compatible because the Titanium needles and it allows the use of intrauterine tandem. Currently, TPS applicators library are not currently available for non-rigid applicators in case of interstitial component as the TB.The purpose of this work is to present the development of a library for the TB, together with its use on a pre-planning technique. Both new goals allow a very efficient and exclusive T2 MRI based planning clinical TB implementation. Methods: The developed library has been implemented in Oncentra Brachytherapy TPS, version 4.3.0 (Elekta) and now is being implemented on Sagiplan v 2.0 TPS (Eckert&Ziegler BEBIG). To model the TB, free and open software named FreeCAD and MeshLab have been used. The reconstruction process is based on three inserted A-vitamin pellets together with the data provided by the free length. The implemented preplanning procedure is as follow: 1) A MRI T2 acquisition is performed with the template in place just with the vaginal cylinder (no uterine tube nor needles). 2) The CTV is drawn and the required needles are selected using a developed Java based application and 3) A post-implant MRI T2 is performed. Results: This library procedure has been successfully applied by now in 25 patients. In this work the use of the developed library will be illustrated with clinical examples. The preplanning procedure has been applied by now in 6 patients, having significant advantages: needle depth estimation, needle positions and number are optimized a priori, time saving, etc Conclusion: TB library and pre-plan techniques are feasible and very efficient and their use will be illustrated in this work.

Recommendations of the American Association of Physicists in Medicine (AAPM) and the European Society for Radiotherapy and Oncology (ESTRO) on dose calculations for high-energy (average energy higher than 50 keV) photon-emitting brachytherapy sources are presented, including the physical characteristics of specific (192)Ir, (137)Cs, and (60)Co source models.This report has been prepared by the High Energy Brachytherapy Source Dosimetry (HEBD) Working Group. This report includes considerations in the application of the TG-43U1 formalism to high-energy photon-emitting sources with particular attention to phantom size effects, interpolation accuracy dependence on dose calculation grid size, and dosimetry parameter dependence on source active length.Consensus datasets for commercially available high-energy photon sources are provided, along with recommended methods for evaluating these datasets. Recommendations on dosimetry characterization methods, mainly using experimental procedures and Monte Carlo, are established and discussed. Also included are methodological recommendations on detector choice, detector energy response characterization and phantom materials, and measurement specification methodology. Uncertainty analyses are discussed and recommendations for high-energy sources without consensus datasets are given.Recommended consensus datasets for high-energy sources have been derived for sources that were commercially available as of January 2010. Data are presented according to the AAPM TG-43U1 formalism, with modified interpolation and extrapolation techniques of the AAPM TG-43U1S1 report for the 2D anisotropy function and radial dose function.

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.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.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 r6 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.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.

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.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.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.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.

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.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.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.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.

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 dataor = 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.

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 (R20.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 Rsphor =10 cm for both 192Ir and 137Cs. Correlation factors to permit comparisons of unbounded g(r) data for ror =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.

Purpose: The Esteya electronic brachytherapy system (Elekta Brachytherapy, The Netherlands) has been recently developed for HDR brachytherapy treatment of skin lesions, using a 69.5 kVp x-ray source and a set of applicators that generate beams of 1 to 3 cm diameter. The purpose of this study was to evaluate the surface dose rate of the Esteya (and its estimated uncertainty) by using different detectors and calibration protocols. Methods: The surface dose rates for the different applicators were measured with the plane parallel ionization chamber T34013 (PTW-Freiburg, Germany) and the parallel plate ionization chamber Exradin A20 (Standard Imaging Inc., USA). The first one was calibrated in both air and water, whereas the second one was provided with calibration in air. Two different chambers of each type were used to evaluate detector reproducibility. Measurements in air were performed according to dosimetric protocol established by the TG-61 of the AAPM, whereas measurements in water were performed in a CIRS Plastic Water LR phantom, using the methodology proposed by the TRS-398 of IAEA. The correction factors of each protocol and ionization chamber were carefully evaluated. Results: Differences between surface dose rates measured in-air and in-water using the T34013 chamber range from 1.6% to 3.3%. Differences are below 3.7% when measurements with the A20 and the T34013 chambers calibrated in air are compared. However, absorbed dose rates measured with the A20 chamber in air and the T34013 chamber in water differ up to 6%. Estimated uncertainties are of the order of 2.5% (with coverage factor k=1). The output factors measured with the PTW chambers differ by less than 1.1% when compared to the values measured with the A20 chamber. Conclusion: The two ionization chambers and the two dosimetric protocols provide consistent results considering the combined uncertainties. Both chambers have negligible stem effect. This study was supported in part by Generalitat Valenciana under Project PROMETEOII/2013/010 and by Spanish Government under Project No. FIS2013-42156. Funding was also received by Elekta and PTW.

Purpose: Esteya™ (Elekta) is a brachytherapy electronic system used for treating skin cancer lesions. This unit is based on a 69.5 kV X-ray source and a surface applicator which produces a circular beam of a specific diameter, which varies from 1 to 3 cm. This study aims to establish the radiation protection measures for the system. Methods: A characterization of the scattered and leakage radiation of the system was implemented by means of experimental measurements. The scattered radiation measurements were performed with a Berthold LB133 ionization chamber, evaluating the dose rate for different distances from the applicator surface. The patient was simulated through water equivalent solid material of 9 cm thickness. Leakage was evaluated by setting 2 cm thick lead at the output of the radiation beam.From the obtained measurements, the necessary shielding has been estimated following the Report 147 by ICRP, and particularizing those values to a hypothetical installation. The dose received by a person who exceptionally would have to eventually stay inside the room with the patient was also analyzed. Results: The scattered dose rate ranges from 17.3 mSv/h at 0.25 m distance to 0.1 mSv/h at 2.5 m. Leakage represents 0.24% of the scattered dose rate. With these measurements and some assumptions on the workload and dimensions of a hypothetical treatment room it was estimated that the overall lead shielding for primary and secondary barriers was well below 1 mm, although it depends on each particular case. The overall absorbed dose to an eventual companion placed 0.5 m away from the patient is below 2 mSv after the complete treatment. Conclusion: Esteya shielding requirements are minimal, which facilitate its integration into conventional rooms. Furthermore, the presence of a companion during treatment is possible if absolutely necessary, recommending the use of a leaded apron and thyroid shielding.

Purpose: To describe how we moved to VMAT in the craniospinal radiotherapy planning process, the actual procedure details, and the results for the patients treated. Methods: Twelve patients underwent craniospinal irradiation with the new procedure, based on the paper by Lee et al. (IJROBP 82, 2012), with some additional modifications. Patients were treated in supine position in Varian Clinac iX linacs with 6 MV RapidArc; prescription doses ranged from 23.4 to 40 Gy (13 to 20 fractions); depending on the PTV length, 2 or 3 isocenters were used, all coordinates being equal except the longitudinal one, setting a few centimeter-long overlapping region; 2 arcs (RA) sharing isocentre for the cranial region, RA1 encompassing cranium and superior spinal region, and RA2 intended to improve conformity, only for cranium; for spine, 1 or 2 isocenters were employed; optimization was performed with Eclipse (V 13.0) using AAA algorithm, establishing sets of optimization parameters to give high conformity while sparing OAR. In pediatric patients, homogeneous irradiation of the vertebrae was also required.Conformity (CI) and heterogeneity (HI) indices (same as Lee et al.), and mean and maximum doses for OAR were calculated. Several pre-treatment verification methods were used: Octavius4D (PTW) for each isocentre, point dose at the junction region, Portal Dosimetry (when possible), and independent MU verification software (Diamond, PTW). Results: CI median value was 1.02 (0.99–1.07) and HI, 1.07 (1.06–1.09); a great reduction was observed for CI and OAR mean doses with respect to Lee et al. data; median maximum eye lens dose was 7.3 Gy (4.0–12.0); mean LungV20Gy was 1.9%; in children, vertebrae were homogeneously irradiated (D95%=20.8 Gy, Dmean= 23.2 Gy).All pre-treatment verifications were found within our action levels except for Portal Dosimetry. Conclusion: A RapidArc planning process for craniospinal axis irradiation has been implemented with significant advantages on conformity, homogeneity, feasibility and efficiency. and increase brain tissue sparing. Variations in volume decrease may be related to shape or location of the tumor.

Purpose: The aim of this study was to obtain equivalent doses to radiosensitive organs when applying high‐dose‐rate (HDR) brachytherapy to the prostate using60 Co or 192 Ir sources, and in comparison to external‐beam radiotherapy (EBRT) modalities. Methods: Monte Carlo simulations in Geant4 were performed using a voxelized adult reference man described in Publication 110 by the International Commission on Radiological Protection (ICRP). Point sources of 60Co or 192Ir with photon energy spectra corresponding to those exiting their capsules were placed in the center of the prostate. Equivalent doses per therapeutic absorbed dose to the prostate were obtained in several radiosensitive organs. Values were corrected to account for realistic source dwell times and positions throughout the prostate. This was repeated for a homogeneous water phantom to assess the adequacy of a homogeneous phantom for HDR brachytherapy dose calculations. Brachytherapy doses were compared to reported values from intensity‐modulated radiation therapy (IMRT) and proton therapy. Results: For the nearest organs considered (bladder, rectum, testes, small intestine, and colon), equivalent doses given by 60 Co source were smaller than from 192Ir. As the distance increased, 60Co delivered higher equivalent doses. Relative to 192Ir, the dose ratio increased up to a factor 5 for the brain. Differences between the heterogeneous and the water phantom became apparent at distances > 30 cm. For the closest organs, equivalent doses from the three EBRT modalities were within an order of magnitude to equivalent dosesdeposited by both brachytherapy sources. However, equivalent doses from brachytherapy to the farthest organs were up to two orders of magnitude smaller than EBRT doses. Conclusion: According to physical considerations, “ Co appears dosimetrically advantageous over 192 Ir sources at close distances, but not for the farthest organs. Equivalent doses to distant healthy organs were lower for brachytherapy than the studied EBRT modalities.

In this study a complete set of dosimetric data for the Plus and 12i Gammamed PDR (pulsed dose rate) 192 Ir sources is presented. These data have been calculated using the Monte Carlo simulation code GEANT3. Absolute dose rate distributions in water around these sources were calculated and are presented in form of conventional two dimensional (2D) Cartesian look-up tables. All dosimetric quantities recommended by the AAPM Task Group 43 report have been also calculated. These quantities are dose rate constant, radial dose function, anisotropy function and anisotropy factor. The dose rate distribution of the 12i source was compared with the corresponding data for the microselectron PDR source showing large differences between both sources.

In this study a complete set of dosimetric data for the GammaMed high dose rate (HDR) 12i and Plus 192 Ir sources are presented. These data have been calculated by means of the Monte Carlo simulation code GEANT3. Absolute dose rate distributions in water are presented as conventional two dimensional (2D) Cartesian look-up tables, and in the TG43 formalism.

Purpose: In the commissioning and QA of surface isotope-based applicators, source-indexer distance (SID) has a great influence in the flatness, symmetry and output. To these purposes, methods described in the literature are the use of a special insert at the entrance of dwell chamber or radiochromic films. Here we present the experience with a micro-chamber array to perform the commissioning and QA of Valencia applicators. Methods: Valencia applicators have been used, the classic and the new extra-shielded version. A micro-chamber array has been employed, 1000 SRS (PTW), with 977 liquid filled, 2.3×2.3×0.5 mm3 sized ion chambers covering 11×11 cm2, which spacing is 2.5 mm in the central 5.5×5.5 cm2, dedicated mainly in principle, in conjunction with Octavius 4D (PTW), to IMRT, VMAT, SBRT verifications. Verisoft software that allows for 3D and planar analysis has been used to evaluate the results. Applicators were located on the surface of the array. To verify the SID, measurements corresponding to the reference value, SID ± 1 mm and SID ± 2 mm were acquired (integration time was fixed in order to discard the influence of the source entrance/exit). Once SID was determined, standard protocol treatments corresponding to 3 Gy and 7 Gy were acquired in order to establish typical patient dose distribution. Results: The method is fast and sensitive. The SID obtained was 1321 mm which is the nominal value included in the applicator manual. For example at 1319 mm an asymmetry of ±8% with respect to the central value was measured, along with a central deviation of −4% referred to 1321 mm. Conclusion: A practical method for the commissioning and QA of Valencia applicators has been described. It has been shown that it is an efficient and accurate tool for these purposes as well as for the verification of the absolute output constancy.

Purpose: The aim of this work is to characterize dosimetrically a clinical intraoperative electron beam accelerator, Mobetron (IntraOp Medical, Inc.) in clinical use in our Hospital. Once this first step is completed our purpose is to evaluate shielding requirements for such a device by preparing adequate phase space files. Methods: It is known that electron beam simulation parameters required for state-of-the-art Monte Carlo codes to obtain a good match with measured data, like the mean energy or the FWHM, may not be code-independent due to the different set of process simulated and formalisms involved. Then, to cross-check our results against any issue in the simulation we have compared experimental data (PDD and profiles for electrons in the range 4 to 12 MeV) with simulations performed independently using both Penelope2011 and Geant4 codes. To do so, the geometry and materials of the head of the accelerator have been fully characterized following information provided by the manufacturer. Results: Both simulations agree with experimental data within experimental uncertainties (±1 mm displacement), although small variations (less than 10%) in the mean energy and FWHM are required to match measured values depending on the code used. Conclusion: Independent Monte Carlo simulations were used to obtain an excellent match to measured electron dose distributions. This opens the road to use such data for evaluating shielding requirements which is the main objective of this project.

Purpose: To study the performance of Dosimetry Check (DC), an EPID-based dosimetry software, which allows performing transit dosimetry, in low density medium, by comparing calculations in-phantom, and analysing results for 15 lung patients. Methods: DC software (v.3.8, pencil beam-based algorithm) has been tested, for plans (Eclipse v.10.0 TPS) delivered in two Varian Clinac iX equipped with aS1000 EPIDs.In the CIRS lung phantom, comparisons between DC and Eclipse (Acuros) were performed for several plans: (1) four field box; (2) square field delivered in arc mode; (3) RapidArc lung patient plan medially centred; (4) RapidArc lung patient plan centred in one lung. Reference points analysed: P1 (medial point, plans 1–3) and P2 (located inside one lung, plan 4).For fifteen lung patients treated with RapidArc, the isocentre and 9 additional points inside the PTV as well as the gamma passing rate (3%/3mm) for the PTV and at the main planes were studied. Results: In-phantom:P1: Per-field differences in plan 1: good agreement for AP-PA fields; discrepancy of 7% for the lateral fields. Global differences (plans 1–3): about 4%, showing a compensating effect of the individual differences.P2: Global difference (plan 4): 15 %. This represents the worst case situation as it is a point surrounded by lung tissue, where the DC pencil beam algorithm is expected to give the greater difference against Acuros.Lung patients: Mean point difference inside the PTV:(5.4±4.2) %. Gamma passing rate inside the PTV:(45±12) %. Conclusion: The performance of DC in heterogeneous lung medium was studied with a special phantom and the results for 15 patients were analysed. The found deviations show that even though DC is a highly promising in vivo dosimetry tool, there is a need of incorporating a more accurate algorithm mainly for plans with low density regions involved.

Purpose: With the establishment of the IGRT as a standard technique, the extra dose that is given to the patients should be taken into account. Furthermore, it has been a recent decrease of the dose threshold in the lens, reduced to 0.5 Gy (ICRP ref 4825-3093-1464 on 21st April, 2011).The purpose of this work was to evaluate the extra dose that the lens is receive due to the Cone-Beam (CBCT) location systems in Head-and-Neck treatments. Methods: The On-Board Imaging (OBI) v 1.5 of the two Varian accelerators, one Clinac iX and one True Beam, were used to obtain the dose that this OBI version give to the lens in the Head-and-Neck location treatments. All CBCT scans were acquired with the Standard Dose Head protocol (100 kVp, 80 mA, 8 ms and 200 degree of rotation).The measurements were taken with thermoluminescence (TLD) EXTRAD (Harshaw) dosimeters placed in an anthropomorphic phantom over the eye and under 3 mm of bolus material to mimic the lens position. The center of the head was placed at the isocenter. To reduce TLD energy dependence, they were calibrated at the used beam quality. Results: The average lens dose at the lens in the OBI v 1.5 systems of the Clinac iX and the True Beam is 0.071 and 0.076 cGy/CBCT, respectively. Conclusions: The extra absorbed doses that receive the eye lenses due to one CBCT acquisition with the studied protocol is far below the new ICRP recommended threshold for the lens. However, the addition effect of several CBCT acquisition during the whole treatment should be taken into account.