Home | << 1 2 3 4 5 6 >> |
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.
|
Langer, C., Algora, A., Couture, A., Csatlos, M., Gulyas, J., Heil, M., et al. (2011). Simulations and developments of the Low Energy Neutron detector Array LENA. Nucl. Instrum. Methods Phys. Res. A, 659(1), 411–418.
Abstract: Prototypes of the Low Energy Neutron detector Array (LENA) have been tested and compared with detailed GEANT simulations. LENA will consist of plastic scintillation bars with the dimensions 1000 x 45 x 10 mm(3). The tests have been performed with gamma-ray sources and neutrons originating from the neutron-induced fission of (235)U. The simulations agreed very well with the measured response and were therefore used to simulate the response to mono-energetic neutrons with different detection thresholds. LENA will be used to detect low-energy neutrons from (p,n)-type reactions with low momentum transfer foreseen at the R(3)B and EXL setups at FAIR, Darmstadt.
|
Hornillos, M. B. G., Gorlychev, V., Caballero, R., Cortes, G., Poch, A., Pretel, C., et al. (2011). Monte Carlo Simulations for the Study of a Moderated Neutron Detector. J. Korean Phys. Soc., 59(2), 1573–1576.
Abstract: This work presents the Monte Carlo simulations performed with the MCNPX and GEANT4 codes for the design of a BEta deLayEd Neutron detector, BELEN-20. This detector will be used for the study of beta delayed neutron emission and consists of a block of polyethylene with dimensions 90 x 90 x 80 cm(3) and 20 cylindrical (3)He gas counters. The results of these simulations have been validated experimentally with a (252)Cf source in the laboratory at UPC, Barcelona. Also the first experiment with this detector has been carried out in November 2009 in JYFL, Finland. In this experiment the neutron emission probability after beta decay of the fission products (88)Br, (94,95)Rb, and (138)I has been measured; this data is still under analysis. Simulations with MCNPX and GEANT4 have been performed in order to obtain the efficiency of the BELEN-20 detector for each of the above nuclei using the neutron energy distribution corresponding to each nucleus.
|
Gimenez-Alventosa, V., Vijande, J., Ballester, F., & Perez-Calatayud, J. (2016). Transit dose comparisons for Co-60 and Ir-192 HDR sources. J. Radiol. Prot., 36(4), 858–864.
Abstract: The goal of this study is to evaluate the ambient dose due to the transit of high dose rate (HDR) Co-60 sources along a transfer tube as compared to Ir-192 ones in a realistic clinical scenario. This goal is accomplished by evaluating air-kerma differences with Monte Carlo calculations using PENELOPE2011. Scatter from both the afterloader and the patient was not taken into account. Two sources, mHDR-v2 and Flexisource Co-60, (Elekta Brachytherapy, Veenendaal, the Netherlands) have been considered. These sources were simulated within a standard transfer tube located in an infinite air phantom. The movement of the source was included by displacing their positions along the connecting tube from z = – 75 cm to z = + 75 cm and combining them. Since modern afterloaders like Flexitron (Elekta) or Saginova (BEBIG GmbH) are able to use equally 192Ir and 60Co sources, it was assumed that both sources are displaced with equal speed. Typical HDR source activity content values were provided by the manufacturer. 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 typical activities. It was found that it varies slowly with distance (less than 10% variation at 75 cm) but strongly in time due to the shorter half-life of the 192Ir (73.83 d). The maximum ratio is located close to the tube. It reaches a value of 0.57 when the typical activity of the sources at the time when they were installed by the vendor was used. Such ratio increases up to 1.28 at the end of the recommended working life (90 d) of the Ir-192 source. Co-60/Ir-192 air-kerma ratios are almost constant (0.51-0.57) in the vicinity of the source-tube with recent installed sources. Nevertheless, air-kerma ratios increase rapidly (1.15-1.29) whenever the Ir-192 is approaching the end of its life. In case of a medical event requiring the medical staff to access the treatment room, these ratios indicate that the dosimetric impact on the medical team will be lower, with a few exceptions, in the case of Co-60-based HDR brachytherapy as compared to Ir-192-based one when typical air-kerma strength values are considered.
Keywords: Monte Carlo; dosimetry; HDR brachytherapy; transit dose
|
Gimenez-Alventosa, V., Antunes, P. C. G., Vijande, J., Ballester, F., Perez-Calatayud, J., & Andreo, P. (2017). Collision-kerma conversion between dose-to-tissue and dose-to-water by photon energy-fluence corrections in low-energy brachytherapy. Phys. Med. Biol., 62(1), 146–164.
Abstract: The AAPM TG-43 brachytherapy dosimetry formalism, introduced in 1995, has become a standard for brachytherapy dosimetry worldwide; it implicitly assumes that charged-particle equilibrium (CPE) exists for the determination of absorbed dose to water at different locations, except in the vicinity of the source capsule. Subsequent dosimetry developments, based on Monte Carlo calculations or analytical solutions of transport equations, do not rely on the CPE assumption and determine directly the dose to different tissues. At the time of relating dose to tissue and dose to water, or vice versa, it is usually assumed that the photon fluence in water and in tissues are practically identical, so that the absorbed dose in the two media can be related by their ratio of mass energy-absorption coefficients. In this work, an efficient way to correlate absorbed dose to water and absorbed dose to tissue in brachytherapy calculations at clinically relevant distances for low-energy photon emitting seeds is proposed. A correction is introduced that is based on the ratio of the water-to-tissue photon energy-fluences. State-of-the art Monte Carlo calculations are used to score photon fluence differential in energy in water and in various human tissues (muscle, adipose and bone), which in all cases include a realistic modelling of low-energy brachytherapy sources in order to benchmark the formalism proposed. The energy-fluence based corrections given in this work are able to correlate absorbed dose to tissue and absorbed dose to water with an accuracy better than 0.5% in the most critical cases (e.g. bone tissue).
|
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.
Keywords: Monte Carlo; ocular brachytherapy; treatment planning
|
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.
Keywords: Monte Carlo; PenRed; Brachytherapy; DICOM; Medical physics
|
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.
|
Etxebeste, A., Dauvergne, D., Fontana, M., Letang, J. M., Llosa, G., Muñoz, E., et al. (2020). CCMod: a GATE module for Compton camera imaging simulation. Phys. Med. Biol., 65(5), 055004–17pp.
Abstract: Compton cameras are gamma-ray imaging systems which have been proposed for a wide variety of applications such as medical imaging, nuclear decommissioning or homeland security. In the design and optimization of such a system Monte Carlo simulations play an essential role. In this work, we propose a generic module to perform Monte Carlo simulations and analyses of Compton Camera imaging which is included in the open-source GATE/Geant4 platform. Several digitization stages have been implemented within the module to mimic the performance of the most commonly employed detectors (e.g. monolithic blocks, pixelated scintillator crystals, strip detectors...). Time coincidence sorter and sequence coincidence reconstruction are also available in order to aim at providing modules to facilitate the comparison and reproduction of the data taken with different prototypes. All processing steps may be performed during the simulation (on-the-fly mode) or as a post-process of the output files (offline mode). The predictions of the module have been compared with experimental data in terms of energy spectra, angular resolution, efficiency and back-projection image reconstruction. Consistent results within a 3-sigma interval were obtained for the energy spectra except for low energies where small differences arise. The angular resolution measure for incident photons of 1275 keV was also in good agreement between both data sets with a value close to 13 degrees. Moreover, with the aim of demonstrating the versatility of such a tool the performance of two different Compton camera designs was evaluated and compared.
Keywords: Monte Carlo; simulation; gamma imaging; Compton camera
|
n_TOF Collaboration(Mendoza, E. et al), Giubrone, G., & Tain, J. L. (2011). Improved Neutron Capture Cross Section Measurements with the n_TOF Total Absorption Calorimeter. J. Korean Phys. Soc., 59(2), 1813–1816.
Abstract: The n_TOF collaboration operates a Total Absorption Calorimeter (TAC) [1] for measuring neutron capture cross-sections of low-mass and/or radioactive samples. The results obtained with the TAC have led to a substantial improvement of the capture cross sections of (237)Np and (240)Pu [2]. The experience acquired during the first measurements has allowed us to optimize the performance of the TAC and to improve the capture signal to background ratio, thus opening the way to more complex and demanding measurements on rare radioactive materials. The new design has been reached by a series of detailed Monte Carlo simulations of complete experiments and dedicated test measurements. The new capture setup will be presented and the main achievements highlighted.
|