|
|
Valdes-Cortez, C., Mansour, I., Rivard, M. J., Ballester, F., Mainegra-Hing, E., Thomson, R. M., et al. (2021). A study of Type B uncertainties associated with the photoelectric effect in low-energy Monte Carlo simulations. Phys. Med. Biol., 66(10), 105014–14pp.
Abstract: Purpose. To estimate Type B uncertainties in absorbed-dose calculations arising from the different implementations in current state-of-the-art Monte Carlo (MC) codes of low-energy photon cross-sections (<200 keV). Methods. MC simulations are carried out using three codes widely used in the low-energy domain: PENELOPE-2018, EGSnrc, and MCNP. Three dosimetry-relevant quantities are considered: mass energy-absorption coefficients for water, air, graphite, and their respective ratios; absorbed dose; and photon-fluence spectra. The absorbed dose and the photon-fluence spectra are scored in a spherical water phantom of 15 cm radius. Benchmark simulations using similar cross-sections have been performed. The differences observed between these quantities when different cross-sections are considered are taken to be a good estimator for the corresponding Type B uncertainties. Results. A conservative Type B uncertainty for the absorbed dose (k = 2) of 1.2%-1.7% (<50 keV), 0.6%-1.2% (50-100 keV), and 0.3% (100-200 keV) is estimated. The photon-fluence spectrum does not present clinically relevant differences that merit considering additional Type B uncertainties except for energies below 25 keV, where a Type B uncertainty of 0.5% is obtained. Below 30 keV, mass energy-absorption coefficients show Type B uncertainties (k = 2) of about 1.5% (water and air), and 2% (graphite), diminishing in all materials for larger energies and reaching values about 1% (40-50 keV) and 0.5% (50-75 keV). With respect to their ratios, the only significant Type B uncertainties are observed in the case of the water-to-graphite ratio for energies below 30 keV, being about 0.7% (k = 2). Conclusions. In contrast with the intermediate (about 500 keV) or high (about 1 MeV) energy domains, Type B uncertainties due to the different cross-sections implementation cannot be considered subdominant with respect to Type A uncertainties or even to other sources of Type B uncertainties (tally volume averaging, manufacturing tolerances, etc). Therefore, the values reported here should be accommodated within the uncertainty budget in low-energy photon dosimetry studies.
|
|
|
|
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.
|
|
|
|
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.
|
|
|
|
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.
|
|
|
|
Perez-Calatayud, J., Ballester, F., Tedgren, C., DeWerd, L. A., Papagiannis, P., Rivard, M. J., et al. (2022). GEC-ESTRO ACROP recommendations on calibration and traceability of HE HDR-PDR photon-emitting brachytherapy sources at the hospital level. Radiother. Oncol., 176, 108–117.
Abstract: The vast majority of radiotherapy departments in Europe using brachytherapy (BT) perform temporary implants of high-or pulsed-dose rate (HDR-PDR) sources with photon energies higher than 50 keV. Such techniques are successfully applied to diverse pathologies and clinical scenarios. These recommen-dations are the result of Working Package 21 (WP-21) initiated within the BRAchytherapy PHYsics Quality Assurance System (BRAPHYQS) GEC-ESTRO working group with a focus on HDR-PDR source cal-ibration. They provide guidance on the calibration of such sources, including practical aspects and issues not specifically accounted for in well-accepted societal recommendations, complementing the BRAPHYQS WP-18 Report dedicated to low energy BT photon emitting sources (seeds). The aim of this report is to provide a European-wide standard in HDR-PDR BT source calibration at the hospital level to maintain high quality patient treatments.
|
|
|
|
Peppa, V., Thomson, R. M., Enger, S. A., Fonseca, G. P., Lee, C. N., Lucero, J. N. E., et al. (2023). A MC-based anthropomorphic test case for commissioning model-based dose calculation in interstitial breast 192-Ir HDR brachytherapy. Med. Phys., 50(7), 4675–4687.
Abstract: PurposeTo provide the first clinical test case for commissioning of Ir-192 brachytherapy model-based dose calculation algorithms (MBDCAs) according to the AAPM TG-186 report workflow. Acquisition and Validation MethodsA computational patient phantom model was generated from a clinical multi-catheter Ir-192 HDR breast brachytherapy case. Regions of interest (ROIs) were contoured and digitized on the patient CT images and the model was written to a series of DICOM CT images using MATLAB. The model was imported into two commercial treatment planning systems (TPSs) currently incorporating an MBDCA. Identical treatment plans were prepared using a generic Ir-192 HDR source and the TG-43-based algorithm of each TPS. This was followed by dose to medium in medium calculations using the MBDCA option of each TPS. Monte Carlo (MC) simulation was performed in the model using three different codes and information parsed from the treatment plan exported in DICOM radiation therapy (RT) format. Results were found to agree within statistical uncertainty and the dataset with the lowest uncertainty was assigned as the reference MC dose distribution. Data Format and Usage NotesThe dataset is available online at ,. Files include the treatment plan for each TPS in DICOM RT format, reference MC dose data in RT Dose format, as well as a guide for database users and all files necessary to repeat the MC simulations. Potential ApplicationsThe dataset facilitates the commissioning of brachytherapy MBDCAs using TPS embedded tools and establishes a methodology for the development of future clinical test cases. It is also useful to non-MBDCA adopters for intercomparing MBDCAs and exploring their benefits and limitations, as well as to brachytherapy researchers in need of a dosimetric and/or a DICOM RT information parsing benchmark. Limitations include specificity in terms of radionuclide, source model, clinical scenario, and MBDCA version used for its preparation.
|
|
|
|
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.
|
|
|
|
Beaulieu, L., Ballester, F., Granero, D., Tedgren, A. C., Haworth, A., Lowenstein, J. R., et al. (2023). AAPM WGDCAB Report 372: A joint AAPM, ESTRO, ABG, and ABS report on commissioning of model-based dose calculation algorithms in brachytherapy. Med. Phys., 50(8), e946–e960.
Abstract: The introduction of model-based dose calculation algorithms (MBDCAs) in brachytherapy provides an opportunity for a more accurate dose calculation and opens the possibility for novel, innovative treatment modalities. The joint AAPM, ESTRO, and ABG Task Group 186 (TG-186) report provided guidance to early adopters. However, the commissioning aspect of these algorithms was described only in general terms with no quantitative goals. This report, from the Working Group on Model-Based Dose Calculation Algorithms in Brachytherapy, introduced a field-tested approach to MBDCA commissioning. It is based on a set of well-characterized test cases for which reference Monte Carlo (MC) and vendor-specific MBDCA dose distributions are available in a Digital Imaging and Communications in Medicine-Radiotherapy (DICOM-RT) format to the clinical users. The key elements of the TG-186 commissioning workflow are now described in detail, and quantitative goals are provided. This approach leverages the well-known Brachytherapy Source Registry jointly managed by the AAPM and the Imaging and Radiation Oncology Core (IROC) Houston Quality Assurance Center (with associated links at ESTRO) to provide open access to test cases as well as step-by-step user guides. While the current report is limited to the two most widely commercially available MBDCAs and only for Ir-192-based afterloading brachytherapy at this time, this report establishes a general framework that can easily be extended to other brachytherapy MBDCAs and brachytherapy sources. The AAPM, ESTRO, ABG, and ABS recommend that clinical medical physicists implement the workflow presented in this report to validate both the basic and the advanced dose calculation features of their commercial MBDCAs. Recommendations are also given to vendors to integrate advanced analysis tools into their brachytherapy treatment planning system to facilitate extensive dose comparisons. The use of the test cases for research and educational purposes is further encouraged.
|
|
|
|
Fletcher, E. M., Ballester, F., Beaulieu, L., Morrison, H., Poher, A., Rivard, M. J., et al. (2024). Generation and comparison of 3D dosimetric reference datasets for COMS eye plaque brachytherapy using model-based dose calculations. Med. Phys., 51, 694–706.
Abstract: PurposeA joint Working Group of the American Association of Physicists in Medicine (AAPM), the European Society for Radiotherapy and Oncology (ESTRO), and the Australasian Brachytherapy Group (ABG) was created to aid in the transition from the AAPM TG-43 dose calculation formalism, the current standard, to model-based dose calculations. This work establishes the first test cases for low-energy photon-emitting brachytherapy using model-based dose calculation algorithms (MBDCAs).Acquisition and Validation MethodsFive test cases are developed: (1) a single model 6711 125I brachytherapy seed in water, 13 seeds (2) individually and (3) in combination in water, (4) the full Collaborative Ocular Melanoma Study (COMS) 16-mm eye plaque in water, and (5) the full plaque in a realistic eye phantom. Calculations are done with four Monte Carlo (MC) codes and a research version of a commercial treatment planning system (TPS). For all test cases, local agreement of MC codes was within & SIM;2.5% and global agreement was & SIM;2% (4% for test case 5). MC agreement was within expected uncertainties. Local agreement of TPS with MC was within 5% for test case 1 and & SIM;20% for test cases 4 and 5, and global agreement was within 0.4% for test case 1 and 10% for test cases 4 and 5.Data Format and Usage NotesDose distributions for each set of MC and TPS calculations are available online () along with input files and all other information necessary to repeat the calculations.Potential ApplicationsThese data can be used to support commissioning of MBDCAs for low-energy brachytherapy as recommended by TGs 186 and 221 and AAPM Report 372. This work additionally lays out a sample framework for the development of test cases that can be extended to other applications beyond eye plaque brachytherapy.
|
|
|
|
Oliver, S., Gimenez-Alventosa, V., Berumen, F., Gimenez, V., Beaulieu, L., Ballester, F., et al. (2023). Benchmark of the PenRed Monte Carlo framework for HDR brachytherapy. Z. Med. Phys., 33(4), 511–528.
Abstract: Purpose: The purpose of this study is to validate the PenRed Monte Carlo framework for clinical applications in brachytherapy. PenRed is a C++ version of Penelope Monte Carlo code with additional tallies and utilities. Methods and materials: Six benchmarking scenarios are explored to validate the use of PenRed and its improved bachytherapy-oriented capabilities for HDR brachytherapy. A new tally allowing the evaluation of collisional kerma for any material using the track length kerma estimator and the possibility to obtain the seed positions, weights and directions processing directly the DICOM file are now implemented in the PenRed distribution. The four non-clinical test cases developed by the Joint AAPM-ESTRO-ABG-ABS WG-DCAB were evaluated by comparing local and global absorbed dose differences with respect to established reference datasets. A prostate and a palliative lung cases, were also studied. For them, absorbed dose ratios, global absorbed dose differences, and cumulative dose-volume histograms were obtained and discussed. Results: The air-kerma strength and the dose rate constant corresponding to the two sources agree with the reference datatests within 0.3% (Sk) and 0.1% (K). With respect to the first three WG-DCAB test cases, more than 99.8% of the voxels present local (global) differences within +/- 1%(+/- 0.1%) of the reference datasets. For test Case 4 reference dataset, more than 94.9%(97.5%) of voxels show an agreement within +/- 1%(+/- 0.1%), better than similar benchmarking calculations in the literature. The track length kerma estimator scorer implemented increases the numerical efficiency of brachytherapy calculations two orders of magnitude, while the specific brachytherapy source allows the user to avoid the use of error-prone intermediate steps to translate the DICOM information into the simulation. In both clinical cases, only minor absorbed dose differences arise in the low-dose isodoses. 99.8% and 100% of the voxels have a global absorbed dose difference ratio within +/- 0.2%for the prostate and lung cases, respectively. The role played by the different segmentation and composition material in the bone structures was discussed, obtaining negligible absorbed dose differ-ences. Dose-volume histograms were in agreement with the reference data.Conclusions: PenRed incorporates new tallies and utilities and has been validated for its use for detailed and precise high-dose-rate brachytherapy simulations.
|
|
