| Home | << 1 2 3 4 5 6 7 >> |
|
Mansour, I. R., Valdes-Cortez, C., Ayala Alvarez, D. S., Berumen, F., Côte, J. S., Ndoutoume-Paquet, G., et al. (2025). Reference datasets for commissioning of model-based dose calculation algorithms for electronic brachytherapy. Med. Phys., , 11pp.
Abstract: PurposeThis work provides the first two clinical test cases for commissioning electronic brachytherapy (eBT) model-based dose calculation algorithms (MBDCAs) for skin irradiation using surface applicators.Acquisition and Validation MethodsThe test cases utilize the INTRABEAM 30 mm surface applicator. Test Case I: water phantom is used to evaluate the algorithm's performance in a uniform medium consisting of a voxelized water cube surrounded by air. Test Case II: Surface eBT represents a heterogeneous medium with four distinct layers: skin tissue, adipose tissue, cortical bone, and soft tissue. Treatment plans for both cases were created and exported into the Radiance treatment planning system (TPS). Dose-to-medium calculations were then performed using this Monte Carlo (MC)-based TPS and compared with MC simulations conducted independently by three different groups using two codes: EGSnrc and PENELOPE. The results agreed within expected Type A and B statistical uncertainties.Data Format and Usage NotesThe dataset is available online at https://doi.org/10.52519/00005. A proprietary file designed for use within Radiance containing CT images and the treatment plan for both test cases, the LINAC modeling, and the CT calibration are included, as well as reference MC and TPS dose data in RTdose format and all files required to run the MC simulations.Potential ApplicationsThis dataset serves as a valuable resource for commissioning eBT MBDCAs and lays the groundwork for developing clinical test cases for other eBT systems. It is also a helpful educational tool for exploring various eBT devices and their advantages and drawbacks. Furthermore, brachytherapy researchers seeking a benchmark for dosimetric calculations in the low-energy domain will find this dataset indispensable.
|
|
Martin-Luna, P., Esperante, D., Casaña, J. V., Fernandez Prieto, A., Fuster-Martinez, N., Rivas, I. G., et al. (2025). Effects of the passive voltage divider in a photomultiplier tube: Analytical model, simulations and experimental validation. Sens. Actuator A-Phys., 381, 116057–11pp.
Abstract: The effects of the passive resistive voltage divider network in a photomultiplier tube (PMT) have been investigated by developing an in-house Monte Carlo simulation code and compared with experimental measurements and an analytical model. The simulation code follows an iterative procedure that takes into account the transport and amplification of the electrons within the device depending on the electrostatic fields produced by the electrode voltages. The PMT gain, dynode voltages, rise time and transit time have been studied as a function of the photocathode current and supply voltage. A good agreement between the analytical model, the simulations and numerous experimental measurements using a Hamamatsu R13408-100 PMT has been obtained. The simulation results endorse the use of logistic functions within the analytical model to account for the collection efficiency in the last dynode stages. This works deepens the understanding of passive voltage dividers and develops an advanced behavioral circuit model of photomultiplier tubes. Although validated fora single PMT, the proposed methodology is applicable to any PMT model. This aids in optimizing the design of fully active voltage dividers, to be applied in extremely pulsed applications with high count rates such as prompt gamma-ray imaging during proton therapy.
|
|
Martin-Luna, P., Esperante, D., Prieto, A. F., Fuster-Martinez, N., Rivas, I. G., Gimeno, B., et al. (2024). Simulation of electron transport and secondary emission in a photomultiplier tube and validation. Sens. Actuator A-Phys., 365, 114859–10pp.
Abstract: The electron amplification and transport within a photomultiplier tube (PMT) has been investigated by developing an in-house Monte Carlo simulation code. The secondary electron emission in the dynodes is implemented via an effective electron model and the Modified Vaughan's model, whereas the transport is computed with the Boris leapfrog algorithm. The PMT gain, rise time and transit time have been studied as a function of supply voltage and external magnetostatic field. A good agreement with experimental measurements using a Hamamatsu R13408-100 PMT was obtained. The simulations have been conducted following different treatments of the underlying geometry: three-dimensional, two-dimensional and intermediate (2.5D). The validity of these approaches is compared. The developed framework will help in understanding the behavior of PMTs under highly intense and irregular illumination or varying external magnetic fields, as in the case of prompt gamma-ray measurements during pencil-beam proton therapy; and aid in optimizing the design of voltage dividers with behavioral circuit models.
|
|
Migueis, F., Casaña, J. V., Garcia-Fernandez, D., Hueso-Gonzalez, F., Llosa, G., Prieto, A. F., et al. (2025). Sensitivity of coaxial prompt gamma-ray monitoring in heterogeneous geometries: A Monte Carlo simulation study. Radiat. Phys. Chem., 232, 112639–11pp.
Abstract: Proton beams offer significant advantages over conventional radiotherapy due to their unique interaction with matter. Specifically, the ionization density caused by these beams is higher in a well-defined region (the Bragg peak) with a sharp decline in intensity beyond a specific depth. However, variations in proton range – often caused by changes in patient anatomy and morphology during treatment – can introduce uncertainties in dose distribution. To account for this, clinicians apply conservative margins, which limit the full potential of proton therapy. Efforts have been focused on developing proton range and dose distribution monitoring systems to reduce the need for large safety margins. These systems are based on detecting and analyzing the byproducts that result from the interaction between the proton beams and tissue. In this article, we focused specifically on a system that aims to detect photons called prompt gamma (PG) rays. We conducted Monte Carlo simulations of proton beams interacting with anthropomorphic phantoms of varying densities to simulate morphological changes. A single scintillation detector was positioned coaxially with the beam and behind the phantom to capture the emitted PG rays in each scenario. Our analysis focused on discrepancies in proton range that resulted from irradiating an anthropomorphic head phantom with varying brain tissue densities and detecting secondary particles resulting from these interactions. We observed potential correlations between gamma-ray signatures and variations in proton range and energy deposition, suggesting that this monitoring technique could be effective for real-world clinical applications.
Keywords: Protons; Coaxial; Proton therapy; Prompt gamma rays; Monte Carlo; Range verification
|
|
Muñoz, E., Barrio, J., Bernabeu, J., Etxebeste, A., Lacasta, C., Llosa, G., et al. (2018). Study and comparison of different sensitivity models for a two-plane Compton camera. Phys. Med. Biol., 63(13), 135004–19pp.
Abstract: Given the strong variations in the sensitivity of Compton cameras for the detection of events originating from different points in the field of view (FoV), sensitivity correction is often necessary in Compton image reconstruction. Several approaches for the calculation of the sensitivity matrix have been proposed in the literature. While most of these models are easily implemented and can be useful in many cases, they usually assume high angular coverage over the scattered photon, which is not the case for our prototype. In this work, we have derived an analytical model that allows us to calculate a detailed sensitivity matrix, which has been compared to other sensitivity models in the literature. Specifically, the proposed model describes the probability of measuring a useful event in a two-plane Compton camera, including the most relevant physical processes involved. The model has been used to obtain an expression for the system and sensitivity matrices for iterative image reconstruction. These matrices have been validated taking Monte Carlo simulations as a reference. In order to study the impact of the sensitivity, images reconstructed with our sensitivity model and with other models have been compared. Images have been reconstructed from several simulated sources, including point-like sources and extended distributions of activity, and also from experimental data measured with Na-22 sources. Results show that our sensitivity model is the best suited for our prototype. Although other models in the literature perform successfully in many scenarios, they are not applicable in all the geometrical configurations of interest for our system. In general, our model allows to effectively recover the intensity of point-like sources at different positions in the FoV and to reconstruct regions of homogeneous activity with minimal variance. Moreover, it can be employed for all Compton camera configurations, including those with low angular coverage over the scatterer.
Keywords: Compton camera imaging; MLEM; Monte Carlo simulations; image quality
|
|
n_TOF Collaboration(Alcayne, V. et al), Balibrea-Correa, J., Domingo-Pardo, C., Lerendegui-Marco, J., Babiano-Suarez, V., & Ladarescu, I. (2024). A Segmented Total Energy Detector (sTED) optimized for (n,γ) cross-section measurements at n_TOF EAR2. Radiat. Phys. Chem., 217, 111525–11pp.
Abstract: The neutron time-of-flight facility nTOF at CERN is a spallation source dedicated to measurements of neutroninduced reaction cross-sections of interest in nuclear technologies, astrophysics, and other applications. Since 2014, Experimental ARea 2 (EAR2) is operational and delivers a neutron fluence of similar to 4 center dot 10(7) neutrons per nominal proton pulse, which is similar to 50 times higher than the one of Experimental ARea 1 (EAR1) of similar to 8 center dot 10(5) neutrons per pulse. The high neutron flux at EAR2 results in high counting rates in the detectors that challenged the previously existing capture detection systems. For this reason, a Segmented Total Energy Detector (sTED) has been developed to overcome the limitations in the detector's response, by reducing the active volume per module and by using a photo-multiplier (PMT) optimized for high counting rates. This paper presents the main characteristics of the sTED, including energy and time resolution, response to gamma-rays, and provides as well details of the use of the Pulse Height Weighting Technique (PHWT) with this detector. The sTED has been validated to perform neutron-capture cross-section measurements in EAR2 in the neutron energy range from thermal up to at least 400 keV. The detector has already been successfully used in several measurements at nTOF EAR2.
Keywords: Neutron capture; PHWT; Scintillation detectors; Monte Carlo simulation
|
|
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.
|
|
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
|
|
Oliver, S., Rodriguez Bosca, S., & Gimenez-Alventosa, V. (2024). Enabling particle transport on CAD-based geometries for radiation simulations with penRed. Comput. Phys. Commun., 298, 109091–11pp.
Abstract: Geometry construction is a fundamental aspect of any radiation transport simulation, regardless of the Monte Carlo code being used. Typically, this process is tedious, time-consuming, and error-prone. The conventional approach involves defining geometries using mathematical objects or surfaces. However, this method comes with several limitations, especially when dealing with complex models, particularly those with organic shapes. Furthermore, since each code employs its own format and methodology for defining geometries, sharing and reproducing simulations among researchers becomes a challenging task. Consequently, many codes have implemented support for simulating over geometries constructed via Computer-Aided Design (CAD) tools. Unfortunately, this feature is lacking in penRed and other PENELOPE physics-based codes. Therefore, the objective of this work is to implement such support within the penRed framework. New version program summary Program Title: Parallel Engine for Radiation Energy Deposition (penRed) CPC Library link to program files: https://doi.org/10.17632/rkw6tvtngy.2 Developer's repository link: https://github.com/PenRed/PenRed Code Ocean capsule: https://codeocean.com/capsule/1041417/tree Licensing provisions: GNU Affero General Public License v3 Programming language: C++ standard 2011. Journal reference of previous version: V. Gimenez-Alventosa, V. Gimenez Gomez, S. Oliver, PenRed: An extensible and parallel Monte-Carlo framework for radiation transport based on PENELOPE, Computer Physics Communications 267 (2021) 108065. doi:https://doi.org/10.1016/j.cpc.2021.108065. Does the new version supersede the previous version?: Yes Reasons for the new version: Implements the capability to simulate on CAD constructed geometries, among many other features and fixes. Summary of revisions: All changes applied through the code versions are summarized in the file CHANGELOG.md in the repository package. Nature of problem: While Monte Carlo codes have proven valuable in simulating complex radiation scenarios, they rely heavily on accurate geometrical representations. In the same way as many other Monte Carlo codes, penRed employs simple geometric quadric surfaces like planes, spheres and cylinders to define geometries. However, since these geometric models offer a certain level of flexibility, these representations have limitations when it comes to simulating highly intricate and irregular shapes. Anatomic structures, for example, require detailed representations of organs, tissues and bones, which are difficult to achieve using basic geometric objects. Similarly, complex devices or intricate mechanical systems may have designs that cannot be accurately represented within the constraints of such geometric models. Moreover, when the complexity of the model increases, geometry construction process becomes more difficult, tedious, time-consuming and error-prone [2]. Also, as each Monte Carlo geometry library uses its own format and construction method, reproducing the same geometry among different codes is a challenging task. Solution method: To face the problems stated above, the objective of this work is to implement the capability to simulate using irregular and adaptable meshed geometries in the penRed framework. This kind of meshes can be constructed using Computer-Aided Design (CAD) tools, the use of which is very widespread and streamline the design process. This feature has been implemented in a new geometry module named “MESH_BODY” specific for this kind of geometries. This one is freely available and usable within the official penRed package1. It can be used since penRed version 1.9.3b and above.
|
|
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.
|