|
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.
|
|
|
Asai, M., Cortes-Giraldo, M. A., Gimenez-Alventosa, V., Gimenez, V., & Salvat, F. (2021). The PENELOPE Physics Models and Transport Mechanics. Implementation into Geant4. Front. Physics, 9, 738735–20pp.
Abstract: A translation of the penelope physics subroutines to C++, designed as an extension of the Geant4 toolkit, is presented. The Fortran code system penelope performs Monte Carlo simulation of coupled electron-photon transport in arbitrary materials for a wide energy range, nominally from 50 eV up to 1 GeV. Penelope implements the most reliable interaction models that are currently available, limited only by the required generality of the code. In addition, the transport of electrons and positrons is simulated by means of an elaborate class II scheme in which hard interactions (involving deflection angles or energy transfers larger than pre-defined cutoffs) are simulated from the associated restricted differential cross sections. After a brief description of the interaction models adopted for photons and electrons/positrons, we describe the details of the class-II algorithm used for tracking electrons and positrons. The C++ classes are adapted to the specific code structure of Geant4. They provide a complete description of the interactions and transport mechanics of electrons/positrons and photons in arbitrary materials, which can be activated from the G4ProcessManager to produce simulation results equivalent to those from the original penelope programs. The combined code, named PenG4, benefits from the multi-threading capabilities and advanced geometry and statistical tools of Geant4.
|
|
|
Navarro, P., Gimeno, B., Alvarez Melcon, A., Arguedas Cuendis, S., Cogollos, C., Diaz-Morcillo, A., et al. (2022). Wide-band full-wave electromagnetic modal analysis of the coupling between dark-matter axions and photons in microwave resonators. Phys. Dark Universe, 36, 101001–14pp.
Abstract: The electromagnetic coupling axion-photon in a microwave cavity is revisited with the Boundary Integral-Resonant Mode Expansion (BI-RME) 3D technique. Such full-wave modal technique has been applied for the rigorous analysis of the excitation of a microwave cavity with an axion field. In this scenario, the electromagnetic field generated by the axion-photon coupling can be assumed to be driven by equivalent electrical charge and current densities. These densities have been inserted in the general BI-RME 3D equations, which express the RF electromagnetic field existing within a cavity as an integral involving the Dyadic Green's functions of the cavity (under Coulomb gauge) as well as such densities. This method is able to take into account any arbitrary spatial and temporal variation of both magnitude and phase of the axion field. Next, we have obtained a simple network driven by the axion current source, which represents the coupling between the axion field and the resonant modes of the cavity. With this approach, it is possible to calculate the extracted and dissipated RF power as a function of frequency along a broad band and without Cauchy-Lorentz approximations, obtaining the spectrum of the electromagnetic field generated in the cavity, and dealing with modes relatively close to the axion resonant mode. Moreover, with this technique we have a complete knowledge of the signal extracted from the cavity, not only in magnitude but also in phase. This can be an interesting issue for future analysis where the axion phase is an important parameter.
|
|
|
Garcfa-Barcelo, J. M., Melcon, A. A., Cuendis, S. A., Diaz-Morcillo, A., Gimeno, B., Kanareykin, A., et al. (2023). On the Development of New Tuning and Inter-Coupling Techniques Using Ferroelectric Materials in the Detection of Dark Matter Axions. IEEE Access, 11, 30360–30372.
Abstract: Tuning is an essential requirement for the search of dark matter axions employing haloscopes since its mass is not known yet to the scientific community. At the present day, most haloscope tuning systems are based on mechanical devices which can lead to failures due to the complexity of the environment in which they are used. However, the electronic tuning making use of ferroelectric materials can provide a path that is less vulnerable to mechanical failures and thus complements and expands current tuning systems. In this work, we present and design a novel technique for using the ferroelectric Potassium Tantalate (KTaO3 or KTO) material as a tuning element in haloscopes based on coupled microwave cavities. In this line, the structures used in the Relic Axion Detector Exploratory Setup (RADES) group are based on several cavities that are connected by metallic irises, which act as interresonator coupling elements. In this article, we also show how to use these KTaO3 films as interresonator couplings between cavities, instead of inductive or capacitive metallic windows used in the past. These two techniques represent a crucial upgrade over the current systems employed in the dark matter axions community, achieving a tuning range of 2.23% which represents a major improvement as compared to previous works (<0.1%) for the same class of tuning systems. The theoretical and simulated results shown in this work demonstrate the interest of the novel techniques proposed for the incorporation of this kind of ferroelectric media in multicavity resonant haloscopes in the search for dark matter axions.
|
|
|
Becchetti, M., Bonciani, R., Cieri, L., Coro, F., & Ripani, F. (2024). Full top-quark mass dependence in diphoton production at NNLO in QCD. Phys. Lett. B, 848, 138362–7pp.
Abstract: In this paper we consider the diphoton production in hadronic collisions at the next-to-next-to-leading order (NNLO) in perturbative QCD, taking into account for the first time the full top quark mass dependence up to two loops (full NNLO). We show selected numerical distributions, highlighting the kinematic regions where the massive corrections are more significant. We make use of the recently computed two-loop massive amplitudes for diphoton production in the quark annihilation channel. The remaining massive contributions at NNLO are also considered, and we comment on the weight of the different types of contributions to the full and complete result.
|
|