|
|
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.
|
|
|
|
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.
|
|
|
|
Gimenez-Alventosa, V., Gimenez, V., & Oliver, S. (2021). PenRed: An extensible and parallel Monte-Carlo framework for radiation transport based on PENELOPE. Comput. Phys. Commun., 267, 108065–12pp.
Abstract: Monte Carlo methods provide detailed and accurate results for radiation transport simulations. Unfortunately, the high computational cost of these methods limits its usage in real-time applications. Moreover, existing computer codes do not provide a methodology for adapting these kinds of simulations to specific problems without advanced knowledge of the corresponding code system, and this restricts their applicability. To help solve these current limitations, we present PenRed, a general-purpose, standalone, extensible and modular framework code based on PENELOPE for parallel Monte Carlo simulations of electron-photon transport through matter. It has been implemented in C++ programming language and takes advantage of modern object-oriented technologies. In addition, PenRed offers the capability to read and process DICOM images as well as to construct and simulate image-based voxelized geometries, so as to facilitate its usage in medical applications. Our framework has been successfully verified against the original PENELOPE Fortran code. Furthermore, the implemented parallelism has been tested showing a significant improvement in the simulation time without any loss in precision of results. Program summary Program title: PenRed: Parallel Engine for Radiation Energy Deposition. CPC Library link to program files: https://doi .org /10 .17632/rkw6tvtngy.1 Licensing provision: GNU Affero General Public License (AGPL). Programming language: C++ standard 2011. Nature of problem: Monte Carlo simulations usually require a huge amount of computation time to achieve low statistical uncertainties. In addition, many applications necessitate particular characteristics or the extraction of specific quantities from the simulation. However, most available Monte Carlo codes do not provide an efficient parallel and truly modular structure which allows users to easily customise their code to suit their needs without an in-depth knowledge of the code system. Solution method: PenRed is a fully parallel, modular and customizable framework for Monte Carlo simulations of the passage of radiation through matter. It is based on the PENELOPE [1] code system, from which inherits its unique physics models and tracking algorithms for charged particles. PenRed has been coded in C++ following an object-oriented programming paradigm restricted to the C++11 standard. Our engine implements parallelism via a double approach: on the one hand, by using standard C++ threads for shared memory, improving the access and usage of the memory, and, on the other hand, via the MPI standard for distributed memory infrastructures. Notice that both kinds of parallelism can be combined together in the same simulation. Moreover, both threads and MPI processes, can be balanced using the builtin load balance system (RUPER-LB [30]) to maximise the performance on heterogeneous infrastructures. In addition, PenRed provides a modular structure with methods designed to easily extend its functionality. Thus, users can create their own independent modules to adapt our engine to their needs without changing the original modules. Furthermore, user extensions will take advantage of the builtin parallelism without any extra effort or knowledge of parallel programming. Additional comments including restrictions and unusual features: PenRed has been compiled in linux systems withg++ of GCC versions 4.8.5, 7.3.1, 8.3.1 and 9; clang version 3.4.2 and intel C++ compiler (icc) version 19.0.5.281. Since it is a C++11-standard compliant code, PenRed should be able to compile with any compiler with C++11 support. In addition, if the code is compiled without MPI support, it does not require any non standard library. To enable MPI capabilities, the user needs to install whatever available MPI implementation, such as openMPI [24] or mpich [25], which can be found in the repositories of any linux distribution. Finally, to provide DICOM processing support, PenRed can be optionally compiled using the dicom toolkit (dcmtk) [32] library. Thus, PenRed has only two optional dependencies, an MPI implementation and the dcmtk library.
|
|
|
|
Hueso-Gonzalez, F., Ballester, F., Perez-Calatayud, J., Siebert, F. A., & Vijande, J. (2017). Towards clinical application of RayStretch for heterogeneity corrections in LDR permanent I-125 prostate brachytherapy. Brachytherapy, 16(3), 616–623.
Abstract: PURPOSE: RayStretch is a simple algorithm proposed for heterogeneity corrections in low-dose-rate brachytherapy. It is built on top of TG-43 consensus data, and it has been validated with Monte Carlo (MC) simulations. In this study, we take a real clinical prostate implant with 71 1251 seeds as reference and we apply RayStretch to analyze its performance in worst-case scenarios. METHODS AND MATERIALS: To do so, we design two cases where large calcifications are located in the prostate lobules. RayStretch resilience under various calcification density values is also explored. Comparisons against MC calculations are performed. RESULTS: Dose volume histogram related parameters like prostate D-90, rectum D-2cc, or urethra D-10 obtained with RayStretch agree within a few percent with the detailed MC results for all cases considered. CONCLUSIONS: The robustness and compatibility of RayStretch with commercial treatment planning systems indicate its applicability in clinical practice for dosimetric corrections in prostate calcifications. Its use during intraoperative ultrasound planning is foreseen.
|
|
|
|
Plaza, J., Martinez, T., Becares, V., Cano-Ott, D., Villamarin, D., de Rada, A. P., et al. (2023). Thermal neutron background at Laboratorio Subterraneo de Canfranc (LSC). Astropart Phys., 146, 102793–9pp.
Abstract: The thermal neutron background at Laboratorio Subterraneo de Canfranc (LSC) has been determined using several He-3 proportional counter detectors. Bare and Cd shielded counters were used in a series of long measurements. Pulse shape discrimination techniques were applied to discriminate between neutron and gamma signals as well as other intrinsic contributions. Montecarlo simulations allowed us to estimate the sensitivity of the detectors and calculate values for the background flux of thermal neutrons inside Hall-A of LSC. The obtained value is (3.5 +/- 0.8)x10(-6) n/cm(2)s, and is within an order of magnitude compared to similar facilities.
|
|
|
|
Molina, F., Aguilera, P., Romero-Barrientos, J., Arellano, H. F., Agramunt, J., Medel, J., et al. (2017). Energy distribution of the neutron flux measurements at the Chilean Reactor RECH-1 using multi-foil neutron activation and the Expectation Maximization unfolding algorithm. Appl. Radiat. Isot., 129, 28–34.
Abstract: We present a methodology to obtain the energy distribution of the neutron flux of an experimental nuclear reactor, using multi-foil activation measurements and the Expectation Maximization unfolding algorithm, which is presented as an alternative to well known unfolding methods such as GRAVEL. Self-shielding flux corrections for energy bin groups were obtained using MCNP6 Monte Carlo simulations. We have made studies at the at the Dry Tube of RECH-1 obtaining fluxes of 1.5(4) x 10(13) cm(-2) s(-1) for the thermal neutron energy region, 1.9(5) x 10(12) cm(-2) s(-1) for the epithermal neutron energy region, and 4.3(11) x 10(11) cm(-2) s(-1) for the fast neutron energy region.
|
|
|
|
Gomez-Cadenas, J. J., Martin-Albo, J., Menendez, J., Mezzetto, M., Monrabal, F., & Sorel, M. (2024). The search for neutrinoless double-beta decay. Riv. Nuovo Cimento, 46, 619–692.
Abstract: Neutrinos are the only particles in the Standard Model that could be Majorana fermions, that is, completely neutral fermions that are their own antiparticles. The most sensitive known experimental method to verify whether neutrinos are Majorana particles is the search for neutrinoless double-beta decay. The last 2 decades have witnessed the development of a vigorous program of neutrinoless double-beta decay experiments, spanning several isotopes and developing different strategies to handle the backgrounds masking a possible signal. In addition, remarkable progress has been made in the understanding of the nuclear matrix elements of neutrinoless double-beta decay, thus reducing a substantial part of the theoretical uncertainties affecting the particle-physics interpretation of the process. On the other hand, the negative results by several experiments, combined with the hints that the neutrino mass ordering could be normal, may imply very long lifetimes for the neutrinoless double-beta decay process. In this report, we review the main aspects of such process, the recent progress on theoretical ideas and the experimental state of the art. We then consider the experimental challenges to be addressed to increase the sensitivity to detect the process in the likely case that lifetimes are much longer than currently explored, and discuss a selection of the most promising experimental efforts.
|
|
|
|
Garcia Canal, C. A., Tarutina, T., & Vento, V. (2023). Analysis of Nuclear Effects in Structure Functions and Their Connection with the Binding Energy of Nuclei. Braz. J. Phys., 53(6), 161–8pp.
Abstract: We describe nuclear effects in structure functions of nuclei in DIS by means of a multiplicative factor beta(A)(x) which differentiates the structure function of the bound nucleons from that of the free nucleons. Our analysis determines that beta(A)(x) establishes a relation between the quark-gluon dynamics expressed by the bound nucleon structure functions and the nuclear dynamics as described by the well-known semi-empirical Bethe-Weizsacker mass formula. This relation corroborates a connection between the underlying quark-gluon dynamics and the phenomenological nuclear dynamics.
|
|
|
|
Fanchiotti, H., Garcia Canal, C. A., Mayosky, M., Veiga, A., & Vento, V. (2023). The Geometric Phase in Classical Systems and in the Equivalent Quantum Hermitian and Non-Hermitian PT-Symmetric Systems. Braz. J. Phys., 53(6), 143–11pp.
Abstract: The decomplexification procedure allows one to show mathematically (stricto sensu) the equivalence (isomorphism) between the quantum dynamics of a system with a finite number of basis states and a classical dynamics system. This unique way of connecting different dynamics was used in the past to analyze the relationship between the well-known geometric phase present in the quantum evolution discovered by Berry and its generalizations, with their analogs, the Hannay phases, in the classical domain. In here, this analysis is carried out for several quantum hermitian and non-hermitian PT-symmetric Hamiltonians and compared with the Hannay phase analysis in their classical isomorphic equivalent systems. As the equivalence ends in the classical domain with oscillator dynamics, we exploit the analogy to propose resonant electric circuits coupled with a gyrator, to reproduce the geometric phase coming from the theoretical solutions, in simulated laboratory experiments.
|
|
|
|
Albiol, F., Corbi, A., & Albiol, A. (2017). Evaluation of modern camera calibration techniques for conventional diagnostic X-ray imaging settings. Radiol. Phys. Technol., 10(1), 68–81.
Abstract: We explore three different alternatives for obtaining intrinsic and extrinsic parameters in conventional diagnostic X-ray frameworks: the direct linear transform (DLT), the Zhang method, and the Tsai approach. We analyze and describe the computational, operational, and mathematical background differences for these algorithms when they are applied to ordinary radiograph acquisition. For our study, we developed an initial 3D calibration frame with tin cross-shaped fiducials at specific locations. The three studied methods enable the derivation of projection matrices from 3D to 2D point correlations. We propose a set of metrics to compare the efficiency of each technique. One of these metrics consists of the calculation of the detector pixel density, which can be also included as part of the quality control sequence in general X-ray settings. The results show a clear superiority of the DLT approach, both in accuracy and operational suitability. We paid special attention to the Zhang calibration method. Although this technique has been extensively implemented in the field of computer vision, it has rarely been tested in depth in common radiograph production scenarios. Zhang's approach can operate on much simpler and more affordable 2D calibration frames, which were also tested in our research. We experimentally confirm that even three or four plane-image correspondences achieve accurate focal lengths.
|
|
