|
Ayala, C., & Cvetic, G. (2016). anQCD: Fortran programs for couplings at complex momenta in various analytic QCD models. Comput. Phys. Commun., 199, 114–117.
Abstract: We provide three Fortran programs which evaluate the QCD analytic (holomorphic) couplings A(v)(Q(2)) for complex or real squared momenta Q(2). These couplings are holomorphic analogs of the powers a(Q(2))(v) of the underlying perturbative QCD (pQCD) coupling a(Q(2)) equivalent to alpha(s)(Q(2))/pi, in three analytic QCD models (anQCD): Fractional Analytic Perturbation Theory (FAPT), Two-delta analytic QCD (2 delta anQCD), and Massive Perturbation Theory (MPT). The index v can be noninteger. The provided programs do basically the same job as the Mathematica package anQCD.m published by us previously (Ayala and Cvetic, 2015), but are now written in Fortran. Program summary Program title: AanQCDext Catalogue identifier: AEYKv10 Program summary URL: http://cpc.cs.qub.ac.uk/summaries/AEYICv1_0.html Program obtainable from: CPC Program Library, Queen's University, Belfast, N. Ireland Licensing provisions: Standard CPC licence, http://cpc.cs.qub.ac.uk/licence/licence.html No. of lines in distributed program, including test data, etc.: 12105 No. of bytes in distributed program, including test data, etc.: 98822 Distribution format: tar.gz Programming language: Fortran. Computer: Any work-station or PC where Fortran 95/200312008 (gfortran) is running. Operating system: Operating system Linux (Ubuntu and Scientific Linux), Windows (in all cases using gfortran). Classification: 11.1, 11.5. Nature of problem: Calculation of values of the running analytic couplings A(v)(Q(2); N-f) for general complex squared momenta Q(2) equivalent to -q(2), in three analytic QCD models, where A(v)(Q(2); N-f) is the analytic (holomorphic) analog of the power (alpha(s)(Q(2); N-f)/pi)(v). Here, A(v)(Q(2); N-f) is a holomorphic function in the Q(2) complex plane, with the exception of the negative semiaxis (-infinity, -M-thr(2)], reflecting the analyticity properties of the spacelike renormalization invariant quantities D(Q(2)) in QCD. In contrast, the perturbative QCD power (alpha(s)(Q(2); N-f)/pi)(v) has singularities even outside the negative semiaxis (Landau ghosts). The three considered models are: Analytic Perturbation theory (APT); Two-delta analytic QCD (2 delta anQCD); Massive Perturbation Theory (MPT). We refer to Ref. [1] for more details and literature. Solution method: The Fortran programs for FAPT and 2 delta anQCD models contain routines and functions needed to perform two-dimensional numerical integrations involving the spectral function, in order to evaluate A(v)(Q(2)) couplings. In MPT model, one-dimensional numerical integration involving A(1)(Q(2)) is sufficient to evaluate any A(v)(Q(2)) coupling. Restrictions: For unphysical choices of the input parameters the results are meaningless. When Q(2) is close to the cut region of the couplings (Q(2) real negative), the calculations can take more time and can have less precision. Running time: For evaluation of a set of about 10 related couplings, the times vary in the range t similar to 10(1)-10(2) s. MPT requires less time, t similar to 1-10(1) s. References: [1] C. Ayala and G. Cvetic, anQCD: a Mathematica package for calculations in general analytic QCD models, Comput. Phys. Commun. 190 (2015) 182.
|
|
|
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.
|
|
|
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.
|
|
|
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.
|
|
|
Allanach, B. C., Martin, S. P., Robertson, D. G., & Ruiz de Austri, R. (2017). The inclusion of two-loop SUSYQCD corrections to gluino and squark pole masses in the minimal and next-to-minimal supersymmetric standard model: SOFTSUSY3.7. Comput. Phys. Commun., 219, 339–345.
Abstract: We describe an extension of the SOFTSUSY spectrum calculator to include two-loop supersymmetric QCD (SUSYQCD) corrections of order O(alpha(2)(s)) to gluino and squark pole masses, either in the minimal supersymmetric standard model (MSSM) or the next-to-minimal supersymmetric standard model (NMSSM). This document provides an overview of the program and acts as a manual for the new version of SOFTSUSY, which includes the increase in accuracy in squark and gluino pole mass predictions. Program summary Program title: SOFTSUSY Program Files doi: http://dx.doLorg/10.17632/sh77x9j7hs.1 Licensing provisions: GNU GPLv3 Programming language: C++, fortran, C Nature of problem: Calculating supersymmetric particle spectrum, mixing parameters and couplings in the MSSM or the NMSSM. The solution to the renormalization group equations must be consistent with theoretical boundary conditions on supersymmetry breaking parameters, as well as a weak-scale boundary condition on gauge couplings, Yukawa couplings and the Higgs potential parameters. Solution method: Nested fixed point iteration. Restrictions: SOFTSUSY will provide a solution only in the perturbative regime and it assumes that all couplings of the model are real (i.e. CP-conserving). If the parameter point under investigation is nonphysical for some reason (for example because the electroWeak potential does not have an acceptable minimum), SOFTSUSY returns an error message. The higher order corrections included are for the MSSM (R-parity conserving or violating) or the real R-parity conserving NMSSM only. Journal reference of previous version: Comput. Phys. Comm. 189 (2015) 192. Does the new version supersede the previous version?: Yes. Reasons for the new version: It is desirable to improve the accuracy of the squark and gluinos mass predictions, since they strongly affect supersymmetric particle production cross-sections at colliders. Summary of revisions: The calculation of the squark and gluino pole masses is extended to be of next-to next-to leading order in SUSYQCD, i.e. including terms up to O(g(s)(4)/(16 pi(2))(2)). Additional comments: Program obtainable from http://softsusy.hepforge.org/
|
|