|
|
Vijande, J., Tedgren, A. C., Ballester, F., Baltas, D., Papagiannis, P., Rivard, M. J., et al. (2021). Source strength determination in iridium-192 and cobalt-60 brachytherapy: A European survey on the level of agreement between clinical measurements and manufacturer certificates. Phys. Imag. Radiat. Oncol., 19, 108–111.
Abstract: Background and purpose: Brachytherapy treatment outcomes depend on the accuracy of the delivered dose distribution, which is proportional to the reference air-kerma rate (RAKR). Current societal recommendations require the medical physicist to compare the measured RAKR values to the manufacturer source calibration certificate. The purpose of this work was to report agreement observed in current clinical practice in the European Union. Materials and methods: A European survey was performed for high- and pulsed-dose-rate (HDR and PDR) highenergy sources (Ir-192 and Co-60), to quantify observed RAKR differences. Medical physicists at eighteen hospitals from eight European countries were contacted, providing 1,032 data points from 2001 to 2020. Results: Over the survey period, 77% of the Ir-192 measurements used a well chamber instead of the older Krieger phantom method. Mean differences with the manufacturer calibration certificate were 0.01% +/- 1.15% for Ir-192 and -0.1% +/- 1.3% for Co-60. Over 95% of RAKR measurements in the clinic were within 3% of the manufacturer calibration certificate. Conclusions: This study showed that the agreement level was generally better than that reflected in prior societal recommendations positing 5%. Future recommendations on high-energy HDR and PDR source calibrations in the clinic may consider tightened agreements levels.
|
|
|
|
Ibanez-Rosello, B., Bautista-Ballesteros, J. A., Candela-Juan, C., Villaescusa, J. I., Ballester, F., Vijande, J., et al. (2017). Evaluation of the shielding in a treatment room with an electronic brachytherapy unit. J. Radiol. Prot., 37(2), N5–N12.
Abstract: Esteya (R) (Elekta Brachytherapy, Veenendaal, The Netherlands) is an electronic brachytherapy (eBT) system based on a 69.5 kVp x-ray source and a set of collimators of 1 to 3 cm in diameter, used for treating non-melanoma skin cancer lesions. This study aims to estimate room shielding requirements for this unit. The non-primary (scattered and leakage) ambient dose equivalent rates were measured with a Berthold LB-133 monitor (Berthold Technologies, Bad Wildbad, Germany). The latter ranges from 17 mSv h(-1) at 0.25 m distance from the x-ray source to 0.1 mSv h(-1) at 2.5 m. The necessary room shielding was then estimated following US and some European guidelines. The room shielding for all barriers considered was below 2 mmPb. The dose to a companion who, exceptionally, would stay with the patient during all treatment was estimated to be below 1 mSv if a leaded apron is used. In conclusion, Esteya shielding requirements are minimal.
|
|
|
|
Hueso-Gonzalez, F., Vijande, J., Ballester, F., Perez-Calatayud, J., & Siebert, F. A. (2015). A simple analytical method for heterogeneity corrections in low dose rate prostate brachytherapy. Phys. Med. Biol., 60(14), 5455–5469.
Abstract: In low energy brachytherapy, the presence of tissue heterogeneities contributes significantly to the discrepancies observed between treatment plan and delivered dose. In this work, we present a simplified analytical dose calculation algorithm for heterogeneous tissue. We compare it with Monte Carlo computations and assess its suitability for integration in clinical treatment planning systems. The algorithm, named as RayStretch, is based on the classic equivalent path length method and TG-43 reference data. Analytical and Monte Carlo dose calculations using Penelope2008 are compared for a benchmark case: a prostate patient with calcifications. The results show a remarkable agreement between simulation and algorithm, the latter having, in addition, a high calculation speed. The proposed analytical model is compatible with clinical real-time treatment planning systems based on TG-43 consensus datasets for improving dose calculation and treatment quality in heterogeneous tissue. Moreover, the algorithm is applicable for any type of heterogeneities.
|
|
|
|
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.
|
|
|
|
Oliver, S., Vijande, J., Tejedor-Aguilar, N., Miro, R., Rovira-Escutia, J. J., Ballester, F., et al. (2023). Monte Carlo flattening filter design to high energy intraoperative electron beam homogenization. Radiat. Phys. Chem., 212, 111102–6pp.
Abstract: Intraoperative radiotherapy using mobile linear accelerators is used for a wide variety of malignancies. However, when large fields are used in combination with high energies, a deterioration of the flatness dose profile is measured with respect to smaller fields and lower energies. Indeed, for the LIAC HWL of Sordina, this deterioration is observed for the 12 MeV beam combined with 10 cm (or larger) diameter applicator. Aimed to solve this problem, a flattening filter has been designed and validated evaluating the feasibility of its usage at the upper part of the applicator. The design of the filter was based on Monte Carlo simulations because of its accuracy in modeling components of clinical devices, among other purposes. The LIAC 10 cm diameter applicator was modeled and simulated independently by two different research groups using two different MC codes, reproducing the heterogeneity of the 12 MeV energy beam. Then, an iterative process of filter design was carried out. Finally, the MC designed conical filter with the optimal size and height to obtain the desired flattened beam was built in-house using a 3D printer. During the experimental validation of the applicator-filter, percentage depth dose, beam profiles, absolute and peripheral dose measurements were performed to demonstrate the effectiveness of the filter addition in the applicator. These measurements conclude that the beam has been flattened, from 5.9% with the standard configuration to 1.6% for the configuration with the filter, without significant increase of the peripheral dose. Consequently, the new filter-applicator LIAC configuration can be used also in a conventional surgery room. A reduction of 16% of the output dose and a reduction of 1.1 mm in the D50 of the percentage depth dose was measured with respect to the original configuration. This work is a proof-of-concept that demonstrates that it is possible to add a filter able to flatten the beam delivered by the Sordina LIAC HWL. Future studies will focus on more refined technical solutions fully compatible with the integrity of the applicator, including its sterilization, to be safely introduced in the clinical practice.
|
|
|
|
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.
|
|
|
|
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.
|
|
|
|
Oliver-Canamas, L., Vijande, J., Candela-Juan, C., Gimeno-Olmos, J., Pujades-Claumarchirant, M. C., Rovira-Escutia, J. J., et al. (2023). A User-Friendly System for Mailed Dosimetric Audits of Ir-192 or Co-60 HDR Brachytherapy Sources. Cancers, 15(9), 2484–14pp.
Abstract: Nowadays, the options available to perform external dosimetric audits of the high dose rate (HDR) brachytherapy treatment process are limited. In this work, we present a methodology that allows for performing dosimetric audits in this field. A phantom was designed and manufactured for this purpose. The criteria for its design, together with the in-house measurements for its characterization, are presented. The result is a user-friendly system that can be mailed to perform dosimetric audits in HDR brachytherapy on-site for systems using either Iridium-192 (Ir-192) or Cobalt-60 (Co-60) sources. Objectives: The main goal of this work is to design and characterize a user-friendly methodology to perform mailed dosimetric audits in high dose rate (HDR) brachytherapy for systems using either Iridium-192 (Ir-192) or Cobalt-60 (Co-60) sources. Methods: A solid phantom was designed and manufactured with four catheters and a central slot to place one dosimeter. Irradiations with an Elekta MicroSelectron V2 for Ir-192, and with a BEBIG Multisource for Co-60 were performed for its characterization. For the dose measurements, nanoDots, a type of optically stimulated luminescent dosimeters (OSLDs), were characterized. Monte Carlo (MC) simulations were performed to evaluate the scatter conditions of the irradiation set-up and to study differences in the photon spectra of different Ir-192 sources (Microselectron V2, Flexisource, BEBIG Ir2.A85-2 and Varisource VS2000) reaching the dosimeter in the irradiation set-up. Results: MC simulations indicate that the surface material on which the phantom is supported during the irradiations does not affect the absorbed dose in the nanoDot. Generally, differences below 5% were found in the photon spectra reaching the detector when comparing the Microselectron V2, the Flexisource and the BEBIG models. However, differences up to 20% are observed between the V2 and the Varisource VS2000 models. The calibration coefficients and the uncertainty in the dose measurement were evaluated. Conclusions: The system described here is able to perform dosimetric audits in HDR brachytherapy for systems using either Ir-192 or Co-60 sources. No significant differences are observed between the photon spectra reaching the detector for the MicroSelectron V2, the Flexisource and the BEBIG Ir-192 sources. For the Varisource VS2000, a higher uncertainty is considered in the dose measurement to allow for the nanoDot response.
|
|
|
|
Assam, I., Vijande, J., Ballester, F., Perez-Calatayud, J., Poppe, B., & Siebert, F. A. (2022). Evaluation of dosimetric effects of metallic artifact reduction and tissue assignment on Monte Carlo dose calculations for I-125 prostate implants. Med. Phys., 49, 6195–6208.
Abstract: Purpose Monte Carlo (MC) simulation studies, aimed at evaluating the magnitude of tissue heterogeneity in I-125 prostate permanent seed implant brachytherapy (BT), customarily use clinical post-implant CT images to generate a virtual representation of a realistic patient model (virtual patient model). Metallic artifact reduction (MAR) techniques and tissue assignment schemes (TAS) are implemented on the post-implant CT images to mollify metallic artifacts due to BT seeds and to assign tissue types to the voxels corresponding to the bright seed spots and streaking artifacts, respectively. The objective of this study is to assess the combined influence of MAR and TAS on MC absorbed dose calculations in post-implant CT-based phantoms. The virtual patient models used for I-125 prostate implant MC absorbed dose calculations in this study are derived from the CT images of an external radiotherapy prostate patient without BT seeds and prostatic calcifications, thus averting the need to implement MAR and TAS. Methods The geometry of the IsoSeed I25.S17plus source is validated by comparing the MC calculated results of the TG-43 parameters for the line source approximation with the TG-43U1S2 consensus data. Four MC absorbed dose calculations are performed in two virtual patient models using the egs_brachy MC code: (1) TG-43-based D-w,w-TG(43), (2) D-w,D-w-MBDC that accounts for interseed scattering and attenuation (ISA), (3) D-m,D-m that examines ISA and tissue heterogeneity by scoring absorbed dose in tissue, and (4) D-w,D-m that unlike D-m,D-m scores absorbed dose in water. The MC absorbed doses (1) and (2) are simulated in a TG-43 patient phantom derived by assigning the densities of every voxel to 1.00 g cm(-3) (water), whereas MC absorbed doses (3) and (4) are scored in the TG-186 patient phantom generated by mapping the mass density of each voxel to tissue according to a CT calibration curve. The MC absorbed doses calculated in this study are compared with VariSeed v8.0 calculated absorbed doses. To evaluate the dosimetric effect of MAR and TAS, the MC absorbed doses of this work (independent of MAR and TAS) are compared to the MC absorbed doses of different I-125 source models from previous studies that were calculated with different MC codes using post-implant CT-based phantoms generated by implementing MAR and TAS on post-implant CT images. Results The very good agreement of TG-43 parameters of this study and the published consensus data within 3% validates the geometry of the IsoSeed I25.S17plus source. For the clinical studies, the TG-43-based calculations show a D-90 overestimation of more than 4% compared to the more realistic MC methods due to ISA and tissue composition. The results of this work generally show few discrepancies with the post-implant CT-based dosimetry studies with respect to the D-90 absorbed dose metric parameter. These discrepancies are mainly Type B uncertainties due to the different I-125 source models and MC codes. Conclusions The implementation of MAR and TAS on post-implant CT images have no dosimetric effect on the I-125 prostate MC absorbed dose calculation in post-implant CT-based phantoms.
|
|
|
|
Valdes-Cortez, C., Niatsetski, Y., Perez-Calatayud, J., Ballester, F., & Vijande, J. (2022). A Monte Carlo study of the relative biological effectiveness in surface brachytherapy. Med. Phys., 49, 5576–5588.
Abstract: Purpose This work aims to simulate clustered DNA damage from ionizing radiation and estimate the relative biological effectiveness (RBE) for radionuclide (rBT)- and electronic (eBT)-based surface brachytherapy through a hybrid Monte Carlo (MC) approach, using realistic models of the sources and applicators. Methods Damage from ionizing radiation has been studied using the Monte Carlo Damage Simulation algorithm using as input the primary electron fluence simulated using a state-of-the-art MC code, PENELOPE-2018. Two Ir-192 rBT applicators, Valencia and Leipzig, one Co-60 source with a Freiburg Flap applicator (reference source), and two eBT systems, Esteya and INTRABEAM, have been included in this study implementing full realizations of their geometries as disclosed by the manufacturer. The role played by filtration and tube kilovoltage has also been addressed. Results For rBT, an RBE value of about 1.01 has been found for the applicators and phantoms considered. In the case of eBT, RBE values for the Esteya system show an almost constant RBE value of about 1.06 for all depths and materials. For INTRABEAM, variations in the range of 1.12-1.06 are reported depending on phantom composition and depth. Modifications in the Esteya system, filtration, and tube kilovoltage give rise to variations in the same range. Conclusions Current clinical practice does not incorporate biological effects in surface brachytherapy. Therefore, the same absorbed dose is administered to the patients independently on the particularities of the rBT or eBT system considered. The almost constant RBE values reported for rBT support that assumption regardless of the details of the patient geometry, the presence of a flattening filter in the applicator design, or even significant modifications in the photon energy spectra above 300 keV. That is not the case for eBT, where a clear dependence on the eBT system and the characteristics of the patient geometry are reported. A complete study specific for each eBT system, including detailed applicator characteristics (size, shape, filtering, among others) and common anatomical locations, should be performed before adopting an existing RBE value.
|
|
