|
|
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.
|
|
|
|
KM3NeT Collaboration(Aiello, S. et al), Alves Garre, S., Calvo, D., Carretero, V., Garcia Soto, A., Gozzini, S. R., et al. (2024). Embedded software of the KM3NeT central logic board. Comput. Phys. Commun., 296, 109036–15pp.
Abstract: The KM3NeT Collaboration is building and operating two deep sea neutrino telescopes at the bottom of the Mediterranean Sea. The telescopes consist of latices of photomultiplier tubes housed in pressure-resistant glass spheres, called digital optical modules and arranged in vertical detection units. The two main scientific goals are the determination of the neutrino mass ordering and the discovery and observation of high-energy neutrino sources in the Universe. Neutrinos are detected via the Cherenkov light, which is induced by charged particles originated in neutrino interactions. The photomultiplier tubes convert the Cherenkov light into electrical signals that are acquired and timestamped by the acquisition electronics. Each optical module houses the acquisition electronics for collecting and timestamping the photomultiplier signals with one nanosecond accuracy. Once finished, the two telescopes will have installed more than six thousand optical acquisition nodes, completing one of the more complex networks in the world in terms of operation and synchronization. The embedded software running in the acquisition nodes has been designed to provide a framework that will operate with different hardware versions and functionalities. The hardware will not be accessible once in operation, which complicates the embedded software architecture. The embedded software provides a set of tools to facilitate remote manageability of the deployed hardware, including safe reconfiguration of the firmware. This paper presents the architecture and the techniques, methods and implementation of the embedded software running in the acquisition nodes of the KM3NeT neutrino telescopes. Program summary Program title: Embedded software for the KM3NeT CLB CPC Library link to program files: https://doi.org/10.17632/s847hpsns4.1 Licensing provisions: GNU General Public License 3 Programming language: C Nature of problem: The challenge for the embedded software in the KM3NeT neutrino telescope lies in orchestrating the Digital Optical Modules (DOMs) to achieve the synchronized data acquisition of the incoming optical signals. The DOMs are the crucial component responsible for capturing neutrino interactions deep underwater. The embedded software must configure and precisely time the operation of each DOM. Any deviation or timing mismatch could compromise data integrity, undermining the scientific value of the experiment. Therefore, the embedded software plays a critical role in coordinating, synchronizing, and operating these modules, ensuring they work in unison to capture and process neutrino signals accurately, ultimately advancing our understanding of fundamental particles in the Universe. Solution method: The embedded software on the DOMs provides a solution based on a C-based bare-metal application, operating without a real-time embedded OS. It is loaded into the RAM during FPGA configuration, consuming less than 256 kB of RAM. The software architecture comprises two layers: system software and application. The former offers OS-like features, including a multitasking scheduler, firmware updates, peripheral drivers, a UDP-based network stack, and error handling utilities. The application layer contains a state machine ensuring consistent program states. It is navigated via slow control events, including external inputs and autonomous responses. Subsystems within the application code control specific acquisition electronics components via the associated driver abstractions. Additional comments including restrictions and unusual features: Due to the operation conditions of the neutrino telescope, where access is restricted, the embedded software implements a fail-safe procedure to reconfigure the firmware where the embedded software runs.
|
|
|
|
Gonzalez-Iglesias, D., Gimeno, B., Esperante, D., Martinez-Reviriego, P., Martin-Luna, P., Fuster-Martinez, N., et al. (2024). Non-resonant ultra-fast multipactor regime in dielectric-assist accelerating structures. Results Phys., 56, 107245–12pp.
Abstract: The objective of this work is the evaluation of the risk of suffering a multipactor discharge in an S-band dielectric-assist accelerating (DAA) structure for a compact low-energy linear particle accelerator dedicated to hadrontherapy treatments. A DAA structure consists of ultra-low loss dielectric cylinders and disks with irises which are periodically arranged in a metallic enclosure, with the advantage of having an extremely high quality factor and very high shunt impedance at room temperature, and it is therefore proposed as a potential alternative to conventional disk-loaded copper structures. However, it has been observed that these structures suffer from multipactor discharges. In fact, multipactor is one of the main problems of these devices, as it limits the maximum accelerating gradient. Because of this, the analysis of multipactor risk in the early design steps of DAA cavities is crucial to ensure the correct performance of the device after fabrication. In this paper, we present a comprehensive and detailed study of multipactor in our DAA design through numerical simulations performed with an in-house developed code based on the Monte-Carlo method. The phenomenology of the multipactor (resonant electron trajectories, electron flight time between impacts, etc.) is described in detail for different values of the accelerating gradient. It has been found that in these structures an ultra-fast non-resonant multipactor appears, which is different from the types of multipactor theoretically studied in the scientific literature. In addition, the effect of several low electron emission coatings on the multipactor threshold is investigated. Furthermore, a novel design based on the modification of the DAA cell geometry for multipactor mitigation is introduced, which shows a significant increase in the accelerating gradient handling capabilities of our prototype.
|
|
|
|
Baran, J. et al, & Brzezinski, K. (2024). Feasibility of the J-PET to monitor the range of therapeutic proton beams. Phys. Medica, 118, 103301–9pp.
Abstract: Purpose: The aim of this work is to investigate the feasibility of the Jagiellonian Positron Emission Tomography (J -PET) scanner for intra-treatment proton beam range monitoring. Methods: The Monte Carlo simulation studies with GATE and PET image reconstruction with CASToR were performed in order to compare six J -PET scanner geometries. We simulated proton irradiation of a PMMA phantom with a Single Pencil Beam (SPB) and Spread -Out Bragg Peak (SOBP) of various ranges. The sensitivity and precision of each scanner were calculated, and considering the setup's cost-effectiveness, we indicated potentially optimal geometries for the J -PET scanner prototype dedicated to the proton beam range assessment. Results: The investigations indicate that the double -layer cylindrical and triple -layer double -head configurations are the most promising for clinical application. We found that the scanner sensitivity is of the order of 10-5 coincidences per primary proton, while the precision of the range assessment for both SPB and SOBP irradiation plans was found below 1 mm. Among the scanners with the same number of detector modules, the best results are found for the triple -layer dual -head geometry. The results indicate that the double -layer cylindrical and triple -layer double -head configurations are the most promising for the clinical application, Conclusions: We performed simulation studies demonstrating that the feasibility of the J -PET detector for PET -based proton beam therapy range monitoring is possible with reasonable sensitivity and precision enabling its pre -clinical tests in the clinical proton therapy environment. Considering the sensitivity, precision and cost-effectiveness, the double -layer cylindrical and triple -layer dual -head J -PET geometry configurations seem promising for future clinical application.
|
|
|
|
Barrientos, L., Borja-Lloret, M., Casaña, J. V., Dendooven, P., Garcia Lopez, J. G., Hueso-Gonzalez, F., et al. (2024). Gamma-ray sources imaging and test-beam results with MACACO III Compton camera. Phys. Medica, 117, 103199–10pp.
Abstract: Hadron therapy is a radiotherapy modality which offers a precise energy deposition to the tumors and a dose reduction to healthy tissue as compared to conventional methods. However, methods for real-time monitoring are required to ensure that the radiation dose is deposited on the target. The IRIS group of IFIC-Valencia developed a Compton camera prototype for this purpose, intending to image the Prompt Gammas emitted by the tissue during irradiation. The system detectors are composed of Lanthanum (III) bromide scintillator crystals coupled to silicon photomultipliers. After an initial characterization in the laboratory, in order to assess the system capabilities for future experiments in proton therapy centers, different tests were carried out in two facilities: PARTREC (Groningen, The Netherlands) and the CNA cyclotron (Sevilla, Spain). Characterization studies performed at PARTREC indicated that the detectors linearity was improved with respect to the previous version and an energy resolution of 5.2 % FWHM at 511 keV was achieved. Moreover, the imaging capabilities of the system were evaluated with a line source of 68Ge and a point-like source of 241Am-9Be. Images at 4.439 MeV were obtained from irradiation of a graphite target with an 18 MeV proton beam at CNA, to perform a study of the system potential to detect shifts at different intensities. In this sense, the system was able to distinguish 1 mm variations in the target position at different beam current intensities for measurement times of 1800 and 600 s.
|
|
