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.
|
NEXT Collaboration(Jones, B. J. P. et al), Carcel, S., Carrion, J. V., Diaz, J., Martin-Albo, J., Martinez, A., et al. (2022). The dynamics of ions on phased radio-frequency carpets in high pressure gases and application for barium tagging in xenon gas time projection chambers. Nucl. Instrum. Methods Phys. Res. A, 1039, 167000–19pp.
Abstract: Radio-frequency (RF) carpets with ultra-fine pitches are examined for ion transport in gases at atmospheric pressures and above. We develop new analytic and computational methods for modeling RF ion transport at densities where dynamics are strongly influenced by buffer gas collisions. An analytic description of levitating and sweeping forces from phased arrays is obtained, then thermodynamic and kinetic principles are used to calculate ion loss rates in the presence of collisions. This methodology is validated against detailed microscopic SIMION simulations. We then explore a parameter space of special interest for neutrinoless double beta decay experiments: transport of barium ions in xenon at pressures from 1 to 10 bar. Our computations account for molecular ion formation and pressure dependent mobility as well as finite temperature effects. We discuss the challenges associated with achieving suitable operating conditions, which lie beyond the capabilities of existing devices, using presently available or near-future manufacturing techniques.
|
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.
|
Angles-Castillo, A., Perucho, M., Marti, J. M., & Laing, R. A. (2021). On the deceleration of Fanaroff-Riley Class I jets: mass loading of magnetized jets by stellar winds. Mon. Not. Roy. Astron. Soc., 500(1), 1512–1530.
Abstract: In this paper, we present steady-state relativistic magnetohydrodynamic simulations that include a mass-load term to study the process of jet deceleration. The mass load mimics the injection of a proton-electron plasma from stellar winds within the host galaxy into initially pair plasma jets, with mean stellar mass-losses ranging from 10(-14) to 10(-9) M-circle dot yr(-1). The spatial jet evolution covers similar to 500 pc from jet injection in the grid at 10 pc from the jet nozzle. Our simulations use a relativistic gas equation of state and a pressure profile for the ambient medium. We compare these simulations with previous dynamical simulations of relativistic, non-magnetized jets. Our results show that toroidal magnetic fields can prevent fast jet expansion and the subsequent embedding of further stars via magnetic tension. In this sense, magnetic fields avoid a runaway deceleration process. Furthermore, when the mass load is large enough to increase the jet density and produce fast, differential jet expansion, the conversion of magnetic energy flux into kinetic energy flux (i.e. magnetic acceleration), helps to delay the deceleration process with respect to non-magnetized jets. We conclude that the typical stellar population in elliptical galaxies cannot explain jet deceleration in classical Fanaroff-Riley type I radio galaxies. However, we observe a significant change in the jet composition, thermodynamical parameters, and energy dissipation along its evolution, even for moderate values of the mass load.
|
Capozzi, F., Ferreira, R. Z., Lopez-Honorez, L., & Mena, O. (2023). CMB and Lyman-alpha constraints on dark matter decays to photons. J. Cosmol. Astropart. Phys., 06(6), 060–23pp.
Abstract: Dark matter energy injection in the early universe modifies both the ionization history and the temperature of the intergalactic medium. In this work, we improve the CMB bounds on sub-keV dark matter and extend previous bounds from Lyman-& alpha; observations to the same mass range, resulting in new and competitive constraints on axion-like particles (ALPs) decaying into two photons. The limits depend on the underlying reionization history, here accounted self-consistently by our modified version of the publicly available DarkHistory and CLASS codes. Future measurements such as the ones from the CMB-S4 experiment may play a crucial, leading role in the search for this type of light dark matter candidates.
|
Watanabe, H., Watanabe, Y. X., Hirayama, Y., Andreyev, A. N., Hashimoto, T., Kondev, F. G., et al. (2021). Beta decay of the axially asymmetric ground state of Re-192. Phys. Lett. B, 814, 136088–6pp.
Abstract: The beta decay of Re-192(75)117, which lies near the boundary between the regions of predicted prolate and oblate deformations, has been investigated using the KEK Isotope Separation System (KISS) in RIKEN Nishina Center. This is the first case in which a low-energy beam of rhenium isotope has been successfully extracted from an argon gas-stopping cell using a laser-ionization technique, following production via multi-nucleon transfer between heavy ions. The ground state of Re-192 has been assigned J(pi) = (0(-)) based on the observed beta feedings and deduced logf t values towards the 0(+) and 2(+) states in Os-192, which is known as a typical gamma-soft nucleus. The shape transition from axial symmetry to axial asymmetry in the Re isotopes is discussed from the viewpoint of single-particle structure using the nuclear Skyrme-Hartree-Fock model.
|
LHCb Collaboration(Aaij, R. et al), Garcia Martin, L. M., Henry, L., Jashal, B. K., Martinez-Vidal, F., Oyanguren, A., et al. (2020). Strong constraints on the b -> s gamma photon polarisation from B-0 -> K(*0)e(+)e(-) decays. J. High Energy Phys., 12(12), 081–25pp.
Abstract: An angular analysis of the B-0 -> K*(0)e(+)e(-) decay is performed using a data sample corresponding to an integrated luminosity of 9 fb(-1) of pp collisions collected with the LHCb experiment. The analysis is conducted in the very low dielectron mass squared (q(2)) interval between 0.0008 and 0.257 GeV2, where the rate is dominated by the B-0 -> K*(0)gamma transition with a virtual photon. The fraction of longitudinal polarisation of the K*(0) meson, F-L, is measured to be F-L = (4.4 +/- 2.6 +/- 1.4)%, where the first uncertainty is statistical and the second systematic. The A(T)(Re) observable, which is related to the lepton forward-backward asymmetry, is measured to be A(T)(Re) = -0.06 +/- 0.08 +/- 0.02. The A(T)((2)) and A(T)(Im) transverse asymmetries, which are sensitive to the virtual photon polarisation, are found to be A(T)((2)) = 0.11 +/- 0.10 +/- 0.02 and A(T)(Im) = 0.02 +/- 0.10 +/- 0.01. The results are consistent with Standard Model predictions and provide the world's best constraint on the b -> s gamma photon polarisation.
|
Vijande, J., Tedgren, A. C., Ballester, F., Baltas, D., Papagiannis, P., Rivard, M. J., et al. (2021). Source strength determination in iridium-192 and cobalt-60 brachytherapy: A European survey on the level of agreement between clinical measurements and manufacturer certificates. Phys. Imag. Radiat. Oncol., 19, 108–111.
Abstract: Background and purpose: Brachytherapy treatment outcomes depend on the accuracy of the delivered dose distribution, which is proportional to the reference air-kerma rate (RAKR). Current societal recommendations require the medical physicist to compare the measured RAKR values to the manufacturer source calibration certificate. The purpose of this work was to report agreement observed in current clinical practice in the European Union. Materials and methods: A European survey was performed for high- and pulsed-dose-rate (HDR and PDR) highenergy sources (Ir-192 and Co-60), to quantify observed RAKR differences. Medical physicists at eighteen hospitals from eight European countries were contacted, providing 1,032 data points from 2001 to 2020. Results: Over the survey period, 77% of the Ir-192 measurements used a well chamber instead of the older Krieger phantom method. Mean differences with the manufacturer calibration certificate were 0.01% +/- 1.15% for Ir-192 and -0.1% +/- 1.3% for Co-60. Over 95% of RAKR measurements in the clinic were within 3% of the manufacturer calibration certificate. Conclusions: This study showed that the agreement level was generally better than that reflected in prior societal recommendations positing 5%. Future recommendations on high-energy HDR and PDR source calibrations in the clinic may consider tightened agreements levels.
|
Araujo Filho, A. A., Furtado, J., Hassanabadi, H., & Reis, J. A. A. S. (2023). Thermal analysis of photon-like particles in rainbow gravity. Phys. Dark Universe, 42, 101310–8pp.
Abstract: This work is devoted to study the thermodynamic behavior of photon-like particles within the rainbow gravity formalism. To to do this, we chose two particular ansatzs to accomplish our calculations. First, we consider a dispersion relation which avoids UV divergences, getting a positive effective cosmological constant. We provide numerical analysis for the thermodynamic functions of the system and bounds are estimated. Furthermore, a phase transition is also expected for this model. Second, we consider a dispersion relation employed in the context of Gamma Ray Bursts. Remarkably, for this latter case, the thermodynamic properties are calculated in an analytical manner and they turn out to depend on the harmonic series Hn, gamma & UGamma; (z), polygamma & psi;n(z) and zeta Riemann functions & zeta;(z).
|
Perez-Perez, J., Amare, J. C., Bandac, I. C., Bayo, A., Borjabad-Sanchez, S., Calvo-Mozota, J. M., et al. (2022). Radon Mitigation Applications at the Laboratorio Subterráneo de Canfranc (LSC). Universe, 8(2), 112–20pp.
Abstract: The Laboratorio Subterraneo de Canfranc (LSC) is the Spanish national hub for low radioactivity techniques and the associated scientific and technological applications. The concentration of the airborne radon is a major component of the radioactive budget in the neighborhood of the detectors. The LSC hosts a Radon Abatement System, which delivers a radon suppressed air with 1.1 & PLUSMN;0.2 mBq/m(3) of Rn-222. The radon content in the air is continuously monitored with an Electrostatic Radon Monitor. Measurements with the double beta decay demonstrators NEXT-NEW and CROSS and the gamma HPGe detectors show the important reduction of the radioactive background due to the purified air in the vicinity of the detectors. We also discuss the use of this facility in the LSC current program which includes NEXT-100, low background biology experiments and radiopure copper electroformation equipment placed in the radon-free clean room.
|
