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.
|
NEXT Collaboration(Renner, J. et al), Kekic, M., Martinez-Lema, G., Alvarez, V., Benlloch-Rodriguez, J. M., Carcel, S., et al. (2019). Energy calibration of the NEXT-White detector with 1% resolution near Q(beta beta) of Xe-136. J. High Energy Phys., 10(10), 230–13pp.
Abstract: Excellent energy resolution is one of the primary advantages of electroluminescent high-pressure xenon TPCs. These detectors are promising tools in searching for rare physics events, such as neutrinoless double-beta decay (beta beta 0 nu), which require precise energy measurements. Using the NEXT-White detector, developed by the NEXT (Neutrino Experiment with a Xenon TPC) collaboration, we show for the first time that an energy resolution of 1% FWHM can be achieved at 2.6 MeV, establishing the present technology as the one with the best energy resolution of all xenon detectors for beta beta 0 nu searches.
|
Antusch, S., Figueroa, D. G., Marschall, K., & Torrenti, F. (2020). Energy distribution and equation of state of the early Universe: Matching the end of inflation and the onset of radiation domination. Phys. Lett. B, 811, 135888–7pp.
Abstract: We study the energy distribution and equation of state of the universe between the end of inflation and the onset of radiation domination (RD), considering observationally consistent single-field inflationary scenarios, with a potential 'flattening' at large field values, and a monomial shape V(phi) proportional to vertical bar phi vertical bar(p) around the origin. As a proxy for (p)reheating, we include a quadratic interaction g(2)phi X-2(2) between the inflaton phi and a light scalar 'daughter' field X, with g(2) > 0. We capture the non-perturbative and non-linear nature of the system dynamics with lattice simulations, obtaining that: i) the final energy transferred to X depends only on p, not on g(2); ii) the final transfer of energy is always negligible for 2 <= p < 4, and of order similar to 50% for p >= 4; iii) the system goes at late times to matter-domination for p = 2, and always to RD for p > 2. In the latter case we calculate the number of e-folds until RD, significantly reducing the uncertainty in the inflationary observables Tl-s and r.
|
Fanchiotti, H., Garcia Canal, C. A., & Vento, V. (2023). Energy loss of monopolium in a medium. Eur. Phys. J. Plus, 138(9), 850–11pp.
Abstract: We study the energy loss of excited monopolium in an atomic medium. We perform a classical calculation in line with a similar calculation performed for charged particles which leads in the non-relativistic limit to the Bethe-Bloch formula except for the density dependence of the medium, which we do not consider in this paper. Our result shows that for maximally deformed Rydberg states, the ionization of monopolium in a light atomic medium is similar to that of light ions.
|
Garcia-Barcelo, J. M., Diaz-Morcillo, A., & Gimeno, B. (2023). Enhancing resonant circular-section haloscopes for dark matter axion detection: approaches and limitations in volume expansion. J. High Energy Phys., 11(11), 159–30pp.
Abstract: Haloscopes, microwave resonant cavities utilized in detecting dark matter axions within powerful static magnetic fields, are pivotal in modern astrophysical research. This paper delves into the realm of cylindrical geometries, investigating techniques to augment volume and enhance compatibility with dipole or solenoid magnets. The study explores volume constraints in two categories of haloscope designs: those reliant on single cavities and those employing multicavities. In both categories, strategies to increase the expanse of elongated structures are elucidated. For multicavities, the optimization of space within magnets is explored through 1D configurations. Three subcavity stacking approaches are investigated, while the foray into 2D and 3D geometries lays the groundwork for future topological developments. The results underscore the efficacy of these methods, revealing substantial room for progress in cylindrical haloscope design. Notably, an elongated single cavity design attains a three-order magnitude increase in volume compared to a WC-109 standard waveguide-based single cavity. Diverse prototypes featuring single cavities, 1D, 2D, and 3D multicavities highlight the feasibility of leveraging these geometries to magnify the volume of tangible haloscope implementations.
|
Khatun, A., Chatterjee, S. S., Thakore, T., & Agarwalla, S. K. (2020). Enhancing sensitivity to non-standard neutrino interactions at INO combining muon and hadron information. Eur. Phys. J. C, 80(6), 533–17pp.
Abstract: In this paper, we explore the impact of flavor violating neutral current non-standard interaction (NSI) parameter epsilon(mu tau) in the oscillation of atmospheric neutrinos and antineutrinos separately using the 50 kt magnetized ICAL detector at INO. We find that due to non-zero epsilon(mu tau), nu(mu) -> nu(mu) and (nu) over bar (mu) -> (nu) over bar (mu) transition probabilities get modified substantially at higher energies and longer baselines, where vacuum oscillation dominates. We demonstrate for the first time that by adding the hadron energy information along with the muon energy and muon direction in each event, the sensitivity of ICAL to the NSI parameter epsilon(mu tau) can be enhanced significantly. The most optimistic bound on epsilon(mu tau) that we obtain is – 0.01 < epsilon(mu tau) < 0.01 at 90% C.L. using 500 kt.yr exposure and considering E-mu, cos theta(mu), and E-had' as observables in their ranges of [1, 21] GeV, [- 1, 1], and [0, 25] GeV, respectively. We discuss for the first time the importance of the charge identification capability of the ICAL detector to have better constraints on epsilon(mu t). We also study the impact of non-zero epsilon(mu tau) on mass hierarchy determination and precision measurement of oscillation parameters.
|
Chakraborty, K., Goswami, S., Gupta, C., & Thakore, T. (2019). Enhancing the hierarchy and octant sensitivity of ESS nu SB in conjunction with T2K, NO nu A and ICAL@INO. J. High Energy Phys., 05(5), 137–26pp.
Abstract: The main aim of the ESSSB proposal is the discovery of the leptonic CP phase (CP) with a high significance (5 sigma for 50% values of (CP)) by utilizing the physics at the second oscillation maxima of the P-e channel. It can achieve 3 sigma sensitivity to hierarchy for all values of (CP). In this work, we concentrate on the hierarchy and octant sensitivity of the ESSSB experiment. We show that combining the ESSSB experiment with the atmospheric neutrino data from the proposed India-based Neutrino Observatory (INO) experiment can result in an increased sensitivity to mass hierarchy. In addition, we also combine the results from the ongoing experiments T2K and NOa assuming their full run-time and present the combined sensitivity of ESSSB + ICAL@INO + T2K + NOA. We show that while by itself ESSSB can have up to 3 sigma hierarchy sensitivity, the combination of all the experiments can give up to 5 sigma sensitivity depending on the true hierarchy-octant combination. The octant sensitivity of ESSSB is low by itself. However the combined sensitivity of all the above experiments can give up to 3 sigma sensitivity depending on the choice of true hierarchy and octant. We discuss the various degeneracies and the synergies that lead to the enhanced sensitivity when combining different experimental data.
|
Pla, S., & Winstanley, E. (2023). Equivalence of the adiabatic expansion and Hadamard renormalization for a charged scalar field. Phys. Rev. D, 107(2), 025004–22pp.
Abstract: We examine the relationship between three approaches (Hadamard, DeWitt-Schwinger, and adiabatic) to the renormalization of expectation values of field operators acting on a charged quantum scalar field. First, we demonstrate that the DeWitt-Schwinger representation of the Feynman Green's function is a particular case of the Hadamard representation. Next, we restrict attention to a spatially flat Friedmann-Lemaitre-Robertson-Walker universe with time-dependent, purely electric, background electromagnetic field, considering two-, three-, and four-dimensional space-times. Working to the order required for the renormalization of the stress-energy tensor, we find the adiabatic and DeWitt-Schwinger expansions of the Green's function when the space-time points are spatially separated. In two and four dimensions, the resulting DeWitt-Schwinger and adiabatic expansions are identical. In three dimensions, the DeWittSchwinger expansion contains terms of adiabatic order 4 that are not necessary for the renormalization of the stress-energy tensor and hence absent in the adiabatic expansion. The equivalence of the DeWittSchwinger and adiabatic approaches to renormalization in the scenario considered is thereby demonstrated in even dimensions. In odd dimensions the situation is less clear and further investigation is required in order to determine whether adiabatic renormalization is a locally covariant renormalization prescription.
|
Ji, T., Dong, X. K., Albaladejo, M., Du, M. L., Guo, F. K., & Nieves, J. (2022). Establishing the heavy quark spin and light flavor molecular multiplets of the X(3872), Z(c)(3900), and X(3960) br. Phys. Rev. D, 106(9), 094002–13pp.
Abstract: Recently, the LHCb Collaboration reported a near-threshold enhancement X(3960) in the D+sD-s invariant mass distribution. We show that the data can be well described by either a bound or a virtual state below the D+sD-s threshold. The mass given by the pole position is (3928 +/- 3) MeV. Using this mass and the existing information on the X(3872) and Zc(3900) resonances, a complete spectrum of the S-wave hadronic molecules formed by a pair of ground state charmed and anticharmed mesons is established. Thus, pole positions of the partners of the X(3872) , Zc(3900) , and the newly observed D+sD-s state are predicted. Calculations have been carried out at the leading order of nonrelativistic effective field theory and considering both heavy quark spin and light flavor SU(3) symmetries, though conservative errors from the breaking of these symmetries are provided.
|
Amoroso, S., Caron, S., Jueid, A., Ruiz de Austri, R., & Skands, P. (2019). Estimating QCD uncertainties in Monte Carlo event generators for gamma-ray dark matter searches. J. Cosmol. Astropart. Phys., 05(5), 007–44pp.
Abstract: Motivated by the recent galactic center gamma-ray excess identified in the Fermi-LAT data, we perform a detailed study of QCD fragmentation uncertainties in the modeling of the energy spectra of gamma-rays from Dark-Matter (DM) annihilation. When Dark-Matter particles annihilate to coloured final states, either directly or via decays such as W(*) -> qq-', photons are produced from a complex sequence of shower, hadronisation and hadron decays. In phenomenological studies their energy spectra are typically computed using Monte Carlo event generators. These results have however intrinsic uncertainties due to the specific model used and the choice of model parameters, which are difficult to asses and which are typically neglected. We derive a new set of hadronisation parameters (tunes) for the PYTHIA 8.2 Monte Carlo generator from a fit to LEP and SLD data at the Z peak. For the first time we also derive a conservative set of uncertainties on the shower and hadronisation model parameters. Their impact on the gamma-ray energy spectra is evaluated and discussed for a range of DM masses and annihilation channels. The spectra and their uncertainties are also provided in tabulated form for future use. The fragmentation-parameter uncertainties may be useful for collider studies as well.
|