|
|
Vijande, J., Granero, D., Perez-Calatayud, J., & Ballester, F. (2013). Monte Carlo dosimetric study of the medium dose rate CSM40 source. Appl. Radiat. Isot., 82, 283–288.
Abstract: The Cs-137 medium dose rate (MDR) CSM40 source model (Eckert & Ziegler BEBIG, Germany) is in clinical use but no dosimetric dataset has been published. This study aims to obtain dosimetric data for the CSM40 source for its use in clinical practice as required by the American Association of Physicists in Medicine (AAPM) and the European Society for Radiotherapy and Oncology (ESTRO). Penelope2008 and Geant4 Monte Carlo codes were used to characterize this source dosimetrically. It was located in an unbounded water phantom with composition and mass density as recommended by AAPM and ESTRO. Due to the low photon energies of Cs-137, absorbed dose was approximated by collisional kerma. Additional simulations were performed to obtain the air-kerma strength, sic. Mass-energy absorption coefficients in water and air were consistently derived and used to calculate collisional kerma. Results performed with both radiation transport codes showed agreement typically within 0.05%. Dose rate constant, radial dose function and anisotropy function are provided for the CSM40 and compared with published data for other commercially available Cs-137 sources. An uncertainty analysis has been performed. The data provided by this study can be used as input data and verification in the treatment planning systems. (C) 2013 Elsevier Ltd. All rights reserved.
|
|
|
|
Vijande, J., Ballester, F., Ouhib, Z., Granero, D., Pujades-Claumarchirant, M. C., & Perez-Calatayud, J. (2012). Dosimetry comparison between TG-43 and Monte Carlo calculations using the Freiburg flap for skin high-dose-rate brachytherapy. Brachytherapy, 11(6), 528–535.
Abstract: PURPOSE: The purpose of this work was to evaluate whether the delivered dose to the skin surface and at the prescription depth when using a Freiburg flap applicator is in agreement with the one predicted by the treatment planning system (TPS) using the TG-43 dose-calculation formalism. METHODS AND MATERIALS: Monte Carlo (MC) simulations and radiochromic film measurements have been performed to obtain dose distributions with the source located at the center of one of the spheres and between two spheres. Primary and scatter dose contributions were evaluated to understand the role played by the scatter component. A standard treatment plan was generated using MC- and TG-43-based TPS applying the superposition principle. RESULTS: The MC model has been validated by performing additional simulations in the same conditions but transforming air and Freiburg flap materials into water to match TG-43 parameters. Both dose distributions differ less than 1%. Scatter defect compared with TG-43 data is up to 15% when the source is located at the center of the sphere and up to 25% when the source is between two spheres. Maximum deviations between TPS- and MC-based distributions are of 5%. CONCLUSIONS: The deviations in the TG-43-based dose distributions for a standard treatment plan with respect to the MC dose distribution calculated taking into account the composition and shape of the applicator and the surrounding air are lower than 5%. Therefore, this study supports the validity of the TPS used in clinical practice. (C) 2012 American Brachytherapy Society. Published by Elsevier Inc. All rights reserved.
|
|
|
|
Pujades, M. C., Granero, D., Vijande, J., Ballester, F., Perez-Calatayud, J., Papagiannis, P., et al. (2014). Air-kerma evaluation at the maze entrance of HDR brachytherapy facilities. J. Radiol. Prot., 34(4), 741–753.
Abstract: In the absence of procedures for evaluating the design of brachytherapy (BT) facilities for radiation protection purposes, the methodology used for external beam radiotherapy facilities is often adapted. The purpose of this study is to adapt the NCRP 151 methodology for estimating the air-kerma rate at the door in BT facilities. Such methodology was checked against Monte Carlo (MC) techniques using the code Geant4. Five different facility designs were studied for Ir-192 and Co-60 HDR applications to account for several different bunker layouts. For the estimation of the lead thickness needed at the door, the use of transmission data for the real spectra at the door instead of the ones emitted by Ir-192 and Co-60 will reduce the lead thickness by a factor of five for Ir-192 and ten for Co-60. This will significantly lighten the door and hence simplify construction and operating requirements for all bunkers. The adaptation proposed in this study to estimate the air-kerma rate at the door depends on the complexity of the maze: it provides good results for bunkers with a maze (i.e. similar to those used for linacs for which the NCRP 151 methodology was developed) but fails for less conventional designs. For those facilities, a specific Monte Carlo study is in order for reasons of safety and cost-effectiveness.
|
|
|
|
Garcia-Cases, F., Perez-Calatayud, J., Ballester, F., Vijande, J., & Granero, D. (2018). Peripheral dose around a mobile linac for intraoperative radiotherapy: radiation protection aspects. J. Radiol. Prot., 38(4), 1393–1411.
Abstract: The aim of this work is to analyse the scattered radiation produced by the mobile accelerator Mobetron 1000. To do so, detailed Monte Carlo simulations using two different codes, Penelope2008 and Geant4, were performed. Measurements were also done. To quantify the attenuation due to the internal structures, present in the accelerator head, on the scattered radiation produced, some of the main structural shielding in the Mobetron 1000 has been incorporated into the geometry simulation. Results are compared with measurements. Some discrepancies between the calculated and measured dose values were found. These differences can be traced back to the importance of the radiation component due to low energy scattered electrons. This encouraged us to perform additional calculations to separate the role played by this component. Ambient dose equivalent, H*(10), outside of the operating room (OR) has been evaluated using Geant4. H*(10) has been measured inside and outside the OR, being its values compatible with those reported in the literature once the low energy electron component is removed. With respect to the role played by neutrons, estimations of neutron H*(10) using Geant4 together with H*(10) measurements has been performed for the case of the 12 MeV electron beam. The values obtained agree with the experimental values existing in the literature, being much smaller than those registered in conventional accelerators. This study is a useful tool for the clinical user to investigate the radiation protection issues arising with the use of these accelerators in ORs without structural shielding. These results will also enable to better fix the maximum number of treatments that could be performed while insuring adequate radiological protection of workers and public in the hospital.
|
|
|
|
Granero, D., Vijande, J., Ballester, F., & Rivard, M. J. (2011). Dosimetry revisited for the HDR Ir-192 brachytherapy source model mHDR-v2. Med. Phys., 38(1), 487–494.
Abstract: Purpose: Recently, the manufacturer of the HDR Ir-192 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 >= 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 >= 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 >= 0.25 cm. If dosimetric data for r<0.25 cm are not needed, dosimetric results from the prior MC studies will be adequate. c 2011 American Association of Physicists in Medicine.
|
|
|
|
Ballester, F., Granero, D., Perez-Calatayud, J., Venselaar, J. L. M., & Rivard, M. J. (2010). Study of encapsulated Tm-170 sources for their potential use in brachytherapy. Med. Phys., 37(4), 1629–1637.
Abstract: Methods: The authors have assumed a theoretical Tm-170 cylindrical source encapsulated with stainless steel and typical dimensions taken from the currently available HDR Ir-192 brachytherapy sources. The dose-rate distribution was calculated for this source using the GEANT4 Monte Carlo (MC) code considering both photon and electron Tm-170 spectra. The AAPM TG-43 U1 brachytherapy dosimetry parameters were derived. To study general properties of Tm-170 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 Tm-170 source were compared to those for Ir-192 and Yb-169 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 Lambda 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 Tm-170, Ir-192, and Yb-169, respectively, resulting in treatment times greater than those for HDR Ir-192 brachytherapy. Conclusions: The dosimetric characteristics of source designs exploiting the low photon energy of Tm-170 were studied for potential application in HDR-brachytherapy. Dose-rate distributions were obtained for cylindrical and simplified spherical Tm-170 source designs (stainless steel and platinum capsule materials) using MC calculations. Despite the high activity of Tm-170, calculated treatment times were much longer than for Ir-192.
|
|
|
|
Rivard, M. J., Granero, D., Perez-Calatayud, J., & Ballester, F. (2010). Influence of photon energy spectra from brachytherapy sources on Monte Carlo simulations of kerma and dose rates in water and air. Med. Phys., 37(2), 869–876.
Abstract: Methods: For Ir-192, I-125, and Pd-103, 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 Ir-192, I-125, and Pd-103 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 Ir-192, I-125, and Pd-103, 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.
|
|
|
|
Granero, D., Perez-Calatayud, J., Vijande, J., Ballester, F., & Rivard, M. J. (2014). Limitations of the TG-43 formalism for skin high-dose-rate brachytherapy dose calculations. Med. Phys., 41(2), 021703–8pp.
Abstract: 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 x 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 Co-60 and Ir-192 sources and a hypothetical Yb-169 source were considered. The Geant4Monte 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 x 5 cm Ir-192 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 Co-60 and Ir-192, yet -3% for Yb-169. For the interstitial implant, dose differences at the skin surface were -7% for Co-60, -0.6% for Ir-192, and -2.5% for Yb-169. 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 Co-60 and Ir-192. For lower energy radionuclides like Yb-169, bolus may be needed; and (iii) for the interstitial case, at least a 0.1 cm bolus is advised for Co-60 to avoid underdosing superficial target layers. For Ir-192 and Yb-169, 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.
|
|
|
|
Ballester, F., Tedgren, A. C., Granero, D., Haworth, A., Mourtada, F., Fonseca, G. P., et al. (2015). A generic high-dose rate Ir-192 brachytherapy source for evaluation of model-based dose calculations beyond the TG-43 formalism. Med. Phys., 42(6), 3048–3062.
Abstract: Purpose: In order to facilitate a smooth transition for brachytherapy dose calculations from the American Association of Physicists in Medicine (AAPM) Task Group No. 43 (TG-43) formalism to model-based dose calculation algorithms (MBDCAs), treatment planning systems (TPSs) using a MBDCA require a set of well-defined test case plans characterized by Monte Carlo (MC) methods. This also permits direct dose comparison to TG-43 reference data. Such test case plans should be made available for use in the software commissioning process performed by clinical end users. To this end, a hypothetical, generic high-dose rate (HDR) Ir-192 source and a virtual water phantom were designed, which can be imported into a TPS. Methods: A hypothetical, generic HDR Ir-192 source was designed based on commercially available sources as well as a virtual, cubic water phantom that can be imported into any TPS in DICOM format. The dose distribution of the generic Ir-192 source when placed at the center of the cubic phantom, and away from the center under altered scatter conditions, was evaluated using two commercial MBDCAs [Oncentra (R) Brachy with advanced collapsed-cone engine (ACE) and BrachyVision AcuRos (TM)]. Dose comparisons were performed using state-of-the-art MC codes for radiation transport, including ALGEBRA, BrachyDose, GEANT4, MCNP5, MCNP6, and pENELopE2008. The methodologies adhered to recommendations in the AAPM TG-229 report on high-energy brachytherapy source dosimetry. TG-43 dosimetry parameters, an along-away dose-rate table, and primary and scatter separated (PSS) data were obtained. The virtual water phantom of (201)(3) voxels (1 mm sides) was used to evaluate the calculated dose distributions. Two test case plans involving a single position of the generic HDR Ir-192 source in this phantom were prepared: (i) source centered in the phantom and (ii) source displaced 7 cm laterally from the center. Datasets were independently produced by different investigators. MC results were then compared against dose calculated using TG-43 and MBDCA methods. Results: TG-43 and PSS datasets were generated for the generic source, the PSS data for use with the ACE algorithm. The dose-rate constant values obtained from seven MC simulations, performed independently using different codes, were in excellent agreement, yielding an average of 1.1109 +/- 0.0004 cGy/(h U) (k = 1, Type A uncertainty). MC calculated dose-rate distributions for the two plans were also found to be in excellent agreement, with differences within type A uncertainties. Differences between commercial MBDCA and MC results were test, position, and calculation parameter dependent. On average, however, these differences were within 1% for ACUROS and 2% for ACE at clinically relevant distances. Conclusions: A hypothetical, generic HDR Ir-192 source was designed and implemented in two commercially available TPSs employing different MBDCAs. Reference dose distributions for this source were benchmarked and used for the evaluation of MBDCA calculations employing a virtual, cubic water phantom in the form of a CT DICOM image series. The implementation of a generic source of identical design in all TPSs using MBDCAs is an important step toward supporting univocal commissioning procedures and direct comparisons between TPSs.
|
|
|
|
Candela-Juan, C., Niatsetski, Y., van der Laarse, R., Granero, D., Ballester, F., Perez-Calatayud, J., et al. (2016). Design and characterization of a new high-dose-rate brachytherapy Valencia applicator for larger skin lesions. Med. Phys., 43(4), 1639–1648.
Abstract: Purpose: The aims of this study were (i) to design a new high-dose-rate (HDR) brachytherapy applicator for treating surface lesions with planning target volumes larger than 3 cm in diameter and up to 5 cm in size, using the microSelectron-HDR or Flexitron afterloader (Elekta Brachytherapy) with a Ir-192 source; (ii) to calculate by means of the Monte Carlo (MC) method the dose distribution for the new applicator when it is placed against a water phantom; and (iii) to validate experimentally the dose distributions in water. Methods: The PENELOPE2008 MC code was used to optimize dwell positions and dwell times. Next, the dose distribution in a water phantom and the leakage dose distribution around the applicator were calculated. Finally, MC data were validated experimentally for a 192Ir mHDR-v2 source by measuring (i) dose distributions with radiochromic EBT3 films (ISP); (ii) percentage depth-dose (PDD) curve with the parallel-plate ionization chamber Advanced Markus (PTW); and (iii) absolute dose rate with EBT3 films and the PinPoint T31016 (PTW) ionization chamber. Results: The new applicator is made of tungsten alloy (Densimet) and consists of a set of interchangeable collimators. Three catheters are used to allocate the source at prefixed dwell positions with preset weights to produce a homogenous dose distribution at the typical prescription depth of 3 mm in water. The same plan is used for all available collimators. PDD, absolute dose rate per unit of air kerma strength, and off-axis profiles in a cylindrical water phantom are reported. These data can be used for treatment planning. Leakage around the applicator was also scored. The dose distributions, PDD, and absolute dose rate calculated agree within experimental uncertainties with the doses measured: differences of MC data with chamber measurements are up to 0.8% and with radiochromic films are up to 3.5%. Conclusions: The new applicator and the dosimetric data provided here will be a valuable tool in clinical practice, making treatment of large skin lesions simpler, faster, and safer. Also the dose to surrounding healthy tissues is minimal.
|
|
