|
|
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.
|
|
|
|
Monerris-Belda, O., Cervera Marin, R., Rodriguez Jodar, M., Diaz-Caballero, E., Alcaide Guillen, C., Petit, J., et al. (2021). High Power RF Discharge Detection Technique Based on the In-Phase and Quadrature Signals. IEEE Trans. Microw. Theory Tech., 69(12), 5429–5438.
Abstract: High power radio frequency (RF) breakdown testing is a subject of great relevance in the space industry, due to the increasing need of higher transmission power and smaller devices. This work presents a novel RF breakdown detection system, which monitors the same parameters as the microwave nulling system but with several advantages. Where microwave nulling-a de facto standard in RF breakdown testing-is narrowband and requires continuous tuning to keep its sensitivity, the proposed technique is broadband and maintains its performance for any RF signal. On top of that, defining the detection threshold is cumbersome due to the lack of an international standardized criterion. Small responses may appear in the detection system during the test and, sometimes, it is not possible to determine if these are an actual RF breakdown or random noise. This new detection system uses a larger analysis bandwidth, thus reducing the cases in which a small response is difficult to be classified. The proposed detection method represents a major step forward in high power testing as it runs without human intervention, warning the operator or decreasing the RF power automatically much faster than any human operator.
|
|
|
|
Gonzalez-Iglesias, D., Esperante, D., Gimeno, B., Boronat, M., Blanch, C., Fuster-Martinez, N., et al. (2021). Analytical RF Pulse Heating Analysis for High Gradient Accelerating Structures. IEEE Trans. Nucl. Sci., 68(2), 78–91.
Abstract: The main aim of this work is to present a simple method, based on analytical expressions, for obtaining the temperature increase due to the Joule effect inside the metallic walls of an RF accelerating component. This technique relies on solving the 1-D heat-transfer equation for a thick wall, considering that the heat sources inside the wall are the ohmic losses produced by the RF electromagnetic fields penetrating the metal with finite electrical conductivity. Furthermore, it is discussed how the theoretical expressions of this method can be applied to obtain an approximation to the temperature increase in realistic 3-D RF accelerating structures, taking as an example the cavity of an RF electron photoinjector and a traveling wave linac cavity. These theoretical results have been benchmarked with numerical simulations carried out with commercial finite-element method (FEM) software, finding good agreement among them. Besides, the advantage of the analytical method with respect to the numerical simulations is evidenced. In particular, the model could be very useful during the design and optimization phase of RF accelerating structures, where many different combinations of parameters must be analyzed in order to obtain the proper working point of the device, allowing to save time and speed up the process. However, it must be mentioned that the method described in this article is intended to provide a quick approximation to the temperature increase in the device, which of course is not as accurate as the proper 3-D numerical simulations of the component.
|
|
|
|
Wilson, J. N. et al, & Algora, A. (2021). Angular momentum generation in nuclear fission. Nature, 590(7847), 566–570.
Abstract: When a heavy atomic nucleus splits (fission), the resulting fragments are observed to emerge spinning(1); this phenomenon has been a mystery in nuclear physics for over 40 years(2,3). The internal generation of typically six or seven units of angular momentum in each fragment is particularly puzzling for systems that start with zero, or almost zero, spin. There are currently no experimental observations that enable decisive discrimination between the many competing theories for the mechanism that generates the angular momentum(4-12). Nevertheless, the consensus is that excitation of collective vibrational modes generates the intrinsic spin before the nucleus splits (pre-scission). Here we show that there is no significant correlation between the spins of the fragment partners, which leads us to conclude that angular momentum in fission is actually generated after the nucleus splits (post-scission). We present comprehensive data showing that the average spin is strongly mass-dependent, varying in saw-tooth distributions. We observe no notable dependence of fragment spin on the mass or charge of the partner nucleus, confirming the uncorrelated post-scission nature of the spin mechanism. To explain these observations, we propose that the collective motion of nucleons in the ruptured neck of the fissioning system generates two independent torques, analogous to the snapping of an elastic band. A parameterization based on occupation of angular momentum states according to statistical theory describes the full range of experimental data well. This insight into the role of spin in nuclear fission is not only important for the fundamental understanding and theoretical description of fission, but also has consequences for the gamma-ray heating problem in nuclear reactors(13,14), for the study of the structure of neutron-rich isotopes(15,16), and for the synthesis and stability of super-heavy elements(17,18). gamma-ray spectroscopy experiments on the origin of spin in the products of nuclear fission of spin-zero nuclei suggest that the fission fragments acquire their spin after scission, rather than before.
|
|
|
|
Khachatryan, M. et al, Coloma, P. (2021). Electron-beam energy reconstruction for neutrino oscillation measurements. Nature, 599(7886), 565–570.
Abstract: Neutrinos exist in one of three types or 'flavours'-electron, muon and tau neutrinos-and oscillate from one flavour to another when propagating through space. This phenomena is one of the few that cannot be described using the standard model of particle physics (reviewed in ref. (1)), and so its experimental study can provide new insight into the nature of our Universe (reviewed in ref. (2)). Neutrinos oscillate as a function of their propagation distance (L) divided by their energy (E). Therefore, experiments extract oscillation parameters by measuring their energy distribution at different locations. As accelerator-based oscillation experiments cannot directly measure E, the interpretation of these experiments relies heavily on phenomenological models of neutrino-nucleus interactions to infer E. Here we exploit the similarity of electron-nucleus and neutrino-nucleus interactions, and use electron scattering data with known beam energies to test energy reconstruction methods and interaction models. We find that even in simple interactions where no pions are detected, only a small fraction of events reconstruct to the correct incident energy. More importantly, widely used interaction models reproduce the reconstructed energy distribution only qualitatively and the quality of the reproduction varies strongly with beam energy. This shows both the need and the pathway to improve current models to meet the requirements of next-generation, high-precision experiments such as Hyper-Kamiokande (Japan)(3) and DUNE (USA)(4). Electron scattering measurements are shown to reproduce only qualitatively state-of-the-art lepton-nucleus energy reconstruction models, indicating that improvements to these particle-interaction models are required to ensure the accuracy of future high-precision neutrino oscillation experiments.
|
|
