|
Martin Lozano, V., Sanda Seoane, R. M., & Zurita, J. (2023). Z'-explorer 2.0: Reconnoitering the dark matter landscape. Comput. Phys. Commun., 288, 108729–14pp.
Abstract: We introduce version 2.0 of Z'-explorer, a software tool that provides a simple, fast, and user-friendly test of models with an extra U (1) gauge boson (Z') against experimental LHC results. The main novelty of the second version is the inclusion of missing energy searches, as the first version only included final states into SM particles. Hence Z'-explorer 2.0 is able to test dark matter models where the Z' acts as an s-channel mediator between the Standard Model and the dark sector, a widespread benchmark employed by the ATLAS and CMS experimental collaborations. To this end, we perform here the first public reinterpretation of the most recent ATLAS mono-jet search with 139 fb-1. In addition, the corresponding searches in the visible final states have also been updated. We illustrate the power of our code by re -obtaining public plots and also showing novel results. In particular, we study the cases where the Z' couples strongly to top quarks (top-philic), where dark matter couples with a mixture of vector and axial-vector couplings, and also perform a scan in the parameter space of a string inspired Stuckelberg model. Z'-explorer 2.0 is publicly available on GitHub. Program summary Program Title: Z'-explorer 2.0 CPC Library link to program files: https://doi .org /10 .17632 /k7tdp8kwgf .2 Developer's repository link: https://github .com /ro -sanda /Z--explorer-2 .0 Licensing provisions: GPLv3 Programming language: C++ and bash Nature of problem: New SM neutral gauge bosons, Z', are ubiquitously present in models of New Physics. In order to confront these models versus a large and ever-growing library of LHC searches, Z'-explorer 1.0 had already included all final states including Standard Model particles. Notably, the previous version of this tool lacked the so-called invisible final states manifested as a momentum imbalance in the transverse plane (“missing energy”). These searches help to probe mediators into a dark sector, where a dark matter candidate resides. Solution method: Z'-explorer encodes the production cross sections for Z' bosons at the LHC as a function of their mass, allowing for a fast evaluation of the exclusion limits. This version of Z'-explorer includes a careful validation of the latest search with one energetic jet (mono-jet) performed by the ATLAS collaboration. Hence one can now test if a given point in parameter space is excluded by both visible and invisible searches. The modular structure of the code has been kept, which allows for potential additions (low-energy constraints, flavor, extrapolation to future colliders).
|
|
|
Figueroa, D. G., Florio, A., Torrenti, F., & Valkenburg, W. (2023). CosmoLattice: A modern code for lattice simulations of scalar and gauge field dynamics in an expanding universe. Comput. Phys. Commun., 283, 108586–13pp.
Abstract: This paper describes CosmoGattice, a modern package for lattice simulations of the dynamics of interacting scalar and gauge fields in an expanding universe. CosmoGattice incorporates a series of features that makes it very versatile and powerful: i) it is written in C++ fully exploiting the object oriented programming paradigm, with a modular structure and a clear separation between the physics and the technical details, ii) it is MPI-based and uses a discrete Fourier transform parallelized in multiple spatial dimensions, which makes it specially appropriate for probing scenarios with well -separated scales, running very high resolution simulations, or simply very long ones, iii) it introduces its own symbolic language, defining field variables and operations over them, so that one can introduce differential equations and operators in a manner as close as possible to the continuum, iv) it includes a library of numerical algorithms, ranging from O(delta t(2)) to O(delta t(10)) methods, suitable for simulating global and gauge theories in an expanding grid, including the case of 'self-consistent' expansion sourced by the fields themselves. Relevant observables are provided for each algorithm (e.g. energy densities, field spectra, lattice snapshots) and we note that, remarkably, all our algorithms for gauge theories (Abelian or non-Abelian) always respect the Gauss constraint to machine precision. Program summary Program Title:: CosmoGattice CPC Library link to program files: https://doi .org /10 .17632 /44vr5xssc6 .1 Developer's repository link: http://github .com /cosmolattice /cosmolattice Licensing provisions: MIT Programming language: C++, MPI Nature of problem: The phenomenology of high energy physics in the early universe is typically characterized by non-linear dynamics, which cannot be captured accurately with analytical techniques. In order to fully understand the non-linearities developed in a given scenario, one needs to carry out lattice simulations. A number of public packages for lattice simulations have appeared over the years, but most of them are only capable of simulating scalar fields. However, realistic models of particle physics do contain other kind of field species, such as (Abelian or non-Abelian) gauge fields, whose non-linear dynamics can also play a relevant role in the early universe. Tensor modes representing gravitational waves are also naturally expected in many scenarios. Solution method: CosmoGattice represents a modern code for lattice simulations of scalar-gauge field theories in an expanding universe. It allows for the simulation of the evolution of interacting (singlet) scalar fields, charged scalar fields under U(1) and/or SU(2) gauge groups, and the corresponding associated Abelian and/or non-Abelian gauge fields. From version 1.1 onward, CosmoGattice also allows to simulate the production of gravitational waves. Simulations can be done either in a flat space-time background, or in a homogeneous and isotropic (spatially flat) expanding FLRW background. CosmoGattice provides symplectic integrators, with accuracy ranging from O (delta t(2)) up to O(delta t(10)), to simuate the non-linear dynamics of the appropriate fields in comoving three-dimensional lattices. The code is parallelized with MPI, and uses a discrete Fourier Transform parallelized in multiple spatial dimensions, which makes it a very powerful code for probing physical problems with well-separated scales. Moreover, the code has been designed as a `platform' to implement any system of dynamical equations suitable for discretization on a lattice.
|
|
|
Fernandez, A., Hufschmidt, D., Colaux, J. L., Valiente-Dobon, J. J., Godinho, V., Jimenez de Haro, M. C., et al. (2020). Low gas consumption fabrication of He-3 solid targets for nuclear reactions. Mater. Des., 186, 108337–10pp.
Abstract: Nanoporous solids that stabilize trapped gas nanobubbles open new possibilities to fabricate solid targets for nuclear reactions. A methodology is described based on the magnetron sputtering (MS) technique operated under quasistatic flux conditions to produce such nanocomposites films with He-3 contents of up to 16 at.% in an amorphous-silicon matrix. In addition to the characteristic low pressure (3-6 Pa) needed for the gas discharge, the method ensures almost complete reduction of the process gas flow during film fabrication. The method could produce similar materials to those obtained under classical dynamic flux conditions for MS. The drastic reduction (>99.5%) of the gas consumption is fundamental for the fabrication of targets with scarce and expensive gases. Si:He-3 and W:He-3 targets are presented together with their microstructural (scanning and transmission electron microscopy, SEM and TEM respectively) and compositional (Ion Beam Analysis, IBA) characterization. The He-3 content achieved was over 1 x 10(18) at/cm(2) for film thicknesses between 1.5 and 3 μm for both Si and W matrices. First experiments to probe the stability of the targets for nuclear reaction studies in inverse kinematics configurations are presented.
|
|
|
Gariazzo, S., de Salas, P. F., Pisanti, O., & Consiglio, R. (2022). PArthENoPE revolutions. Comput. Phys. Commun., 271, 108205–13pp.
Abstract: This paper presents the main features of a new and updated version of the program PArthENoPE, which the community has been using for many years for computing the abundances of light elements produced during Big Bang Nucleosynthesis. This is the third release of the PArthENoPE code, after the 2008 and the 2018 ones, and will be distributed from the code's website, http://parthenope.na.infn.it. Apart from minor changes, the main improvements in this new version include a revisited implementation of the nuclear rates for the most important reactions of deuterium destruction, H-2(p,gamma) He-3, H-2(d, n)He-3 and H-2(d, p)H-3, and a re-designed GUI, which extends the functionality of the previous one. The new GUI, in particular, supersedes the previous tools for running over grids of parameters with a better management of parallel runs, and it offers a brand-new set of functions for plotting the results. Program summary Program title: PArthENoPE 3.0 CPC Library link to program files: https://doi.org/10.17632/wygr7d8yt9.2 Developer's repository link: http://parthenope.na.infn.it Licensing provisions: GPLv3 Programming language: Fortran 77 and Python Nature of problem: Computation of yields of light elements synthesized in the primordial universe Solution method: Livermore Solver for Ordinary Differential Equations (LSODE) for stiff and nonstiff systems, Python GUI for running and plotting Journal reference of previous version: Comput. Phys. Commun. 233 (2018) 237-242 Does the new version supersede the previous version?: Yes Reasons for the new version: Update of the physics and improvements in the GUI Summary of revisions: Update of the physics implemented in the Fortran code and improvements in the GUI functionalities, in particular new plotting functions.
|
|
|
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.
|
|
|
Aiola, S., Amhis, Y., Billoir, P., Jashal, B. K., Henry, L., Oyanguren, A., et al. (2021). Hybrid seeding: A standalone track reconstruction algorithm for scintillating fibre tracker at LHCb. Comput. Phys. Commun., 260, 107713–5pp.
Abstract: We describe the Hybrid seeding, a stand-alone pattern recognition algorithm aiming at finding charged particle trajectories for the LHCb upgrade. A significant improvement to the charged particle reconstruction efficiency is accomplished by exploiting the knowledge of the LHCb magnetic field and the position of energy deposits in the scintillating fibre tracker detector. Moreover, we achieve a low fake rate and a small contribution to the overall timing budget of the LHCb real-time data processing.
|
|
|
KM3NeT Collaboration(Aiello, S. et al), Alves Garre, S., Calvo, D., Carretero, V., Colomer, M., Corredoira, I., et al. (2020). gSeaGen: The KM3NeT GENIE-based code for neutrino telescopes. Comput. Phys. Commun., 256, 107477–15pp.
Abstract: The gSeaGen code is a GENIE-based application developed to efficiently generate high statistics samples of events, induced by neutrino interactions, detectable in a neutrino telescope. The gSeaGen code is able to generate events induced by all neutrino flavours, considering topological differences between tracktype and shower-like events. Neutrino interactions are simulated taking into account the density and the composition of the media surrounding the detector. The main features of gSeaGen are presented together with some examples of its application within the KM3NeT project. Program summary Program Title: gSeaGen CPC Library link to program files: http://dx.doi.org/10.17632/ymgxvy2br4.1 Licensing provisions: GPLv3 Programming language: C++ External routines/libraries: GENIE [1] and its external dependencies. Linkable to MUSIC [2] and PROPOSAL [3]. Nature of problem: Development of a code to generate detectable events in neutrino telescopes, using modern and maintained neutrino interaction simulation libraries which include the state-of-the-art physics models. The default application is the simulation of neutrino interactions within KM3NeT [4]. Solution method: Neutrino interactions are simulated using GENIE, a modern framework for Monte Carlo event generators. The GENIE framework, used by nearly all modern neutrino experiments, is considered as a reference code within the neutrino community. Additional comments including restrictions and unusual features: The code was tested with GENIE version 2.12.10 and it is linkable with release series 3. Presently valid up to 5 TeV. This limitation is not intrinsic to the code but due to the present GENIE valid energy range. References: [1] C. Andreopoulos at al., Nucl. Instrum. Meth. A614 (2010) 87. [2] P. Antonioli et al., Astropart. Phys. 7 (1997) 357. [3] J. H. Koehne et al., Comput. Phys. Commun. 184 (2013) 2070. [4] S. Adrian-Martinez et al., J. Phys. G: Nucl. Part. Phys. 43 (2016) 084001.
|
|
|
KM3NeT Collaboration(Aiello, S. et al), Calvo, D., Coleiro, A., Colomer, M., Gozzini, S. R., Hernandez-Rey, J. J., et al. (2020). The Control Unit of the KM3NeT Data Acquisition System. Comput. Phys. Commun., 256, 107433–16pp.
Abstract: The KM3NeT Collaboration runs a multi-site neutrino observatory in the Mediterranean Sea. Water Cherenkov particle detectors, deep in the sea and far off the coasts of France and Italy, are already taking data while incremental construction progresses. Data Acquisition Control software is operating off-shore detectors as well as testing and qualification stations for their components. The software, named Control Unit, is highly modular. It can undergo upgrades and reconfiguration with the acquisition running. Interplay with the central database of the Collaboration is obtained in a way that allows for data taking even if Internet links fail. In order to simplify the management of computing resources in the long term, and to cope with possible hardware failures of one or more computers, the KM3NeT Control Unit software features a custom dynamic resource provisioning and failover technology, which is especially important for ensuring continuity in case of rare transient events in multi-messenger astronomy. The software architecture relies on ubiquitous tools and broadly adopted technologies and has been successfully tested on several operating systems.
|
|
|
Gonzalez-Iglesias, D., Gimeno, B., Esperante, D., Martinez-Reviriego, P., Martin-Luna, P., Fuster-Martinez, N., et al. (2024). Non-resonant ultra-fast multipactor regime in dielectric-assist accelerating structures. Results Phys., 56, 107245–12pp.
Abstract: The objective of this work is the evaluation of the risk of suffering a multipactor discharge in an S-band dielectric-assist accelerating (DAA) structure for a compact low-energy linear particle accelerator dedicated to hadrontherapy treatments. A DAA structure consists of ultra-low loss dielectric cylinders and disks with irises which are periodically arranged in a metallic enclosure, with the advantage of having an extremely high quality factor and very high shunt impedance at room temperature, and it is therefore proposed as a potential alternative to conventional disk-loaded copper structures. However, it has been observed that these structures suffer from multipactor discharges. In fact, multipactor is one of the main problems of these devices, as it limits the maximum accelerating gradient. Because of this, the analysis of multipactor risk in the early design steps of DAA cavities is crucial to ensure the correct performance of the device after fabrication. In this paper, we present a comprehensive and detailed study of multipactor in our DAA design through numerical simulations performed with an in-house developed code based on the Monte-Carlo method. The phenomenology of the multipactor (resonant electron trajectories, electron flight time between impacts, etc.) is described in detail for different values of the accelerating gradient. It has been found that in these structures an ultra-fast non-resonant multipactor appears, which is different from the types of multipactor theoretically studied in the scientific literature. In addition, the effect of several low electron emission coatings on the multipactor threshold is investigated. Furthermore, a novel design based on the modification of the DAA cell geometry for multipactor mitigation is introduced, which shows a significant increase in the accelerating gradient handling capabilities of our prototype.
|
|
|
del Rio, A., & Agullo, I. (2023). Chiral fermion anomaly as a memory effect. Phys. Rev. D, 108(10), 105025–22pp.
Abstract: We study the nonconservation of the chiral charge of Dirac fields between past and future null infinity due to the Adler-Bell-Jackiw chiral anomaly. In previous investigations [A. del Rio, Phys. Rev. D 104, 065012 (2021)], we found that this charge fails to be conserved if electromagnetic sources in the bulk emit circularly polarized radiation. In this article, we unravel yet another contribution coming from the nonzero, infrared “soft” charges of the external, electromagnetic field. This new contribution can be interpreted as another manifestation of the ordinary memory effect produced by transitions between different infrared sectors of Maxwell theory, but now on test quantum fields rather than on test classical particles. In other words, a flux of electromagnetic waves can leave a memory on quantum fermion states in the form of a permanent, net helicity. We elaborate this idea in both 1 + 1 and 3 + 1 dimensions. We also show that, in sharp contrast, gravitational infrared charges do not contribute to the fermion chiral anomaly.
|
|