|
|
Ferrer-Sanchez, A., Martin-Guerrero, J., Ruiz de Austri, R., Torres-Forne, A., & Font, J. A. (2024). Gradient-annihilated PINNs for solving Riemann problems: Application to relativistic hydrodynamics. Comput. Meth. Appl. Mech. Eng., 424, 116906–18pp.
Abstract: We present a novel methodology based on Physics-Informed Neural Networks (PINNs) for solving systems of partial differential equations admitting discontinuous solutions. Our method, called Gradient-Annihilated PINNs (GA-PINNs), introduces a modified loss function that forces the model to partially ignore high-gradients in the physical variables, achieved by introducing a suitable weighting function. The method relies on a set of hyperparameters that control how gradients are treated in the physical loss. The performance of our methodology is demonstrated by solving Riemann problems in special relativistic hydrodynamics, extending earlier studies with PINNs in the context of the classical Euler equations. The solutions obtained with the GA-PINN model correctly describe the propagation speeds of discontinuities and sharply capture the associated jumps. We use the relative l(2) error to compare our results with the exact solution of special relativistic Riemann problems, used as the reference ''ground truth'', and with the corresponding error obtained with a second-order, central, shock-capturing scheme. In all problems investigated, the accuracy reached by the GA-PINN model is comparable to that obtained with a shock-capturing scheme, achieving a performance superior to that of the baseline PINN algorithm in general. An additional benefit worth stressing is that our PINN-based approach sidesteps the costly recovery of the primitive variables from the state vector of conserved variables, a well-known drawback of grid-based solutions of the relativistic hydrodynamics equations. Due to its inherent generality and its ability to handle steep gradients, the GA-PINN methodology discussed in this paper could be a valuable tool to model relativistic flows in astrophysics and particle physics, characterized by the prevalence of discontinuous solutions.
|
|
|
|
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.
|
|
|
|
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.
|
|
|
|
Oberhauser, B. B., Bisio, P., Celentano, A., Depero, E., Dusaev, R. R., Kirpichnikov, D. V., et al. (2024). Development of the fully Geant4 compatible package for the simulation of Dark Matter in fixed target experiments. Comput. Phys. Commun., 300, 109199–11pp.
Abstract: The search for new comparably light (well below the electroweak scale) feebly interacting particles is an exciting possibility to explain some mysterious phenomena in physics, among them the origin of Dark Matter. The sensitivity study through detailed simulation of projected experiments is a key point in estimating their potential for discovery. Several years ago we created the DMG4 package for the simulation of DM (Dark Matter) particles in fixed target experiments. The natural approach is to integrate this simulation into the same program that performs the full simulation of particles in the experiment setup. The Geant4 toolkit framework was chosen as the most popular and versatile solution nowadays. The simulation of DM particles production by this package accommodates several possible scenarios, employing electron, muon or photon beams and involving various mediators, such as vector, axial vector, scalar, pseudoscalar, or spin 2 particles. The bremsstrahlung, annihilation or Primakoff processes can be simulated. The package DMG4 contains a subpackage DarkMatter with cross section methods weakly connected to Geant4. It can be used in different frameworks. In this paper, we present the latest developments of the package, such as extending the list of possible mediator particle types, refining formulas for the simulation and extending the mediator mass range. The user interface is also made more flexible and convenient. In this work, we also demonstrate the usage of the package, the improvements in the simulation accuracy and some cross check validations. Program summary Program title: DMG4 CPC Library link to program files: https://doi .org /10 .17632 /cmr4bcrj6j .1 Licensing provisions: GNU General Public License 3 Programming language: c++ Journal reference of previous version: Comput. Phys. Commun. 269 (2021) 108129 Does the new version supersede the previous version?: Yes Reasons for the new version: Numerous developments, addition of new features Summary of revisions: WW approximation cross sections for the muon beam are implemented and cross-checked, models with semivisible A ' (inelastic Dark Matter) and spin 2 mediators are added. The range of possible mediator masses is extended. Several important improvements for the annihilation processes are made, the number of possible annihilation processes is extended. User interface is improved. Several bugs are fixed. Nature of problem: For the simulation of Dark Matter production processes in fixed target experiments a code that can be easily integrated in programs for the full simulation of experimental setup is needed. Solution method: A fully Geant4 compatible DM simulation package DMG4 was presented in 2020. We present numerous further developments of this package.
|
|
|
|
KM3NeT Collaboration(Aiello, S. et al), Alves Garre, S., Bariego-Quintana, A., Calvo, D., Carretero, V., Cecchini, V., et al. (2024). The Power Board of the KM3NeT Digital Optical Module: Design, Upgrade, and Production. Electronics, 13(11), 2044–17pp.
Abstract: The KM3NeT Collaboration is building an underwater neutrino observatory at the bottom of the Mediterranean Sea, consisting of two neutrino telescopes, both composed of a three-dimensional array of light detectors, known as digital optical modules. Each digital optical module contains a set of 31 three-inch photomultiplier tubes distributed over the surface of a 0.44 m diameter pressure-resistant glass sphere. The module also includes calibration instruments and electronics for power, readout, and data acquisition. The power board was developed to supply power to all the elements of the digital optical module. The design of the power board began in 2013, and ten prototypes were produced and tested. After an exhaustive validation process in various laboratories within the KM3NeT Collaboration, a mass production batch began, resulting in the construction of over 1200 power boards so far. These boards were integrated in the digital optical modules that have already been produced and deployed, which total 828 as of October 2023. In 2017, an upgrade of the power board, to increase reliability and efficiency, was initiated. The validation of a pre-production series has been completed, and a production batch of 800 upgraded boards is currently underway. This paper describes the design, architecture, upgrade, validation, and production of the power board, including the reliability studies and tests conducted to ensure safe operation at the bottom of the Mediterranean Sea throughout the observatory's lifespan.
|
|
|
|
Lerendegui-Marco, J., Babiano-Suarez, V., Balibrea-Correa, J., Caballero, L., Calvo, D., Ladarescu, I., et al. (2024). Simultaneous Gamma-Neutron Vision device: a portable and versatile tool for nuclear inspections. EPJ Tech. Instrum., 11(1), 2–17pp.
Abstract: This work presents GN-Vision, a novel dual gamma-ray and neutron imaging system, which aims at simultaneously obtaining information about the spatial origin of gamma-ray and neutron sources. The proposed device is based on two position sensitive detection planes and exploits the Compton imaging technique for the imaging of gamma-rays. In addition, spatial distributions of slow- and thermal-neutron sources (<100 eV) are reconstructed by using a passive neutron pin-hole collimator attached to the first detection plane. The proposed gamma-neutron imaging device could be of prime interest for nuclear safety and security applications. The two main advantages of this imaging system are its high efficiency and portability, making it well suited for nuclear applications were compactness and real-time imaging is important. This work presents the working principle and conceptual design of the GN-Vision system and explores, on the basis of Monte Carlo simulations, its simultaneous gamma-ray and neutron detection and imaging capabilities for a realistic scenario where a Cf-252 source is hidden in a neutron moderating container.
|
|
|
|
Tagliente, G. et al, Babiano-Suarez, V., Domingo-Pardo, C., Ladarescu, I., & Tain, J. L. (2024). High-resolution cross section measurements for neutron interactions on 89Y with incident neutron energies up to 95 keV. Eur. Phys. J. A, 60(1), 21–18pp.
Abstract: The cross section of the Y-89(n,gamma) reaction has important implications in nuclear astrophysics and for advanced nuclear technology. Given its neutron magic number N = 50 and a consequent small neutron capture crosssection,89Y represents one of the key nuclides for the stellars-process. It acts as a bottleneck in the neutron capture chain between the Fe seed and the heavier elements. Moreover, it is located at the overlapping region, where both the weak and mains-process components take place.Y-89, the only stable yttrium isotope, is also used in innovative nuclear reactors. Neutron capture and transmission measurements were per-formed at the time-of-flight facilities n_TOF at CERN and GELINA at JRC-Geel. Resonance parameters of individual resonances were extracted from a resonance analysis of the experimental transmission and capture yields, up to a neutron incident energy of 95 keV. Even though a comparison with results reported in the literature shows differences in resonance parameters, the present data are consistent with the Maxwellian averaged cross section suggested by the astro-physical database KADoNiS.
|
|
|
|
n_TOF Collaboration(Wright, T. et al), Domingo-Pardo, C., & Tain, J. L. (2024). Measurement of the prompt fission γ-rays from slow neutron-induced fission of 235U with STEFF. Eur. Phys. J. A, 60(3), 70–11pp.
Abstract: The amount of energy carried by gamma-rays during the fission process is an important consideration when developing new reactor designs. Many studies of gamma-ray energy and multiplicity, from a multitude of fissioning systems, were measured during the 1970s. However the data from such experiments largely underestimates the heating effect caused by gamma-rays in the structure of a reactor. It is therefore essential to obtain more accurate measurements of the energy carried during gamma-ray emission. As such, the OECD Nuclear Energy Agency has put out a high priority request [1] for measurements of the mean gamma-ray energy and multiplicity to an accuracy better than 7.5 percent from several fissioning systems; including U-235(n(thermal)). Measurements of the rays from these fissioning nuclei were performed with the SpecTrometer for Exotic Fission Fagments (STEFF).
|
|
|
|
Abreu, L. M., Song, J., Brandao, P. C. S., & Oset, E. (2024). A note on the tensor and vector exchange contributions to K (K)over-bar → K (K)over-bar, D(D)over-bar → D(D)over-bar and π+π- → π+π- reactions. Eur. Phys. J. A, 60(3), 76–10pp.
Abstract: In this note we study the tensor and vector exchange contributions to the elastic reactions involving the pseudoscalars mesons pi(+) pi(-), K+ K- and D+D-. In the case of the tensor-exchange contributions we assume that an intermediate tensor f(2)(1270) is dynamically generated from the interaction of two virtual rho mesons, with the use of a pole approximation. The calculation of the two-loop amplitude is facilitated since the triangle loops can be factorized and computed separately. The results show very small contributions coming from the tensor-exchange mechanisms when compared with those from the vector-exchange processes. We compare our results for pi pi and K (K) over bar scattering with those obtained in other works where the f2(1270) is considered as an ordinary q (q) over bar meson. Our picture provides a smaller contribution but of similar order of magnitude for pion scattering and stabilizes the results in the case of K (K) over bar, allowing us to make estimates for D (D) over bar scattering.
|
|
|
|
Lin, J. X., Li, J. T., Liang, W. H., Chen, H. X., & Oset, E. (2024). J/ψ decays into ω(φ) f1(1285) and ω(φ) “f1(1420)”. Eur. Phys. J. C, 84(1), 52–8pp.
Abstract: We perform a theoretical study of the J/psi -> omega (Phi)K* K + c.c. -> omega(Phi)K-0 pi+ K- reactions with the assumption that the f1(1285) is dynamically generated from a single channel K*K + c.c interaction in the chiral unitary approach. Two peaks in the K-0 pi+ K- invariant mass distribution are observed, one clear peak locates at the f(1)(1285) nominal mass, the other peak locates at around 1420MeV with about 70MeV width. We conclude that the former peak is associated with the f(1)(1285) and the latter peak is not a genuine resonance but a manifestation of the kinematic effect in the higher energy region caused by the K* K + c.c. decay mode of the f(1)(1285).
|
|
