|
Athron, P., Park, J. H., Stockinger, D., & Voigt, A. (2015). FlexibleSUSY-A spectrum generator generator for supersymmetric models. Comput. Phys. Commun., 190, 139–172.
Abstract: We introduce FlexibleSUSY, a Mathematica and C++ package, which generates a fast, precise C++ spectrum generator for any SUSY model specified by the user. The generated code is designed with both speed and modularity in mind, making it easy to adapt and extend with new features. The model is specified by supplying the superpotential, gauge structure and particle content in a SARAH model file; specific boundary conditions e.g. at the GUT, weak or intermediate scales are defined in a separate FlexibleSUSY model file. From these model files, FlexibleSUSY generates C++ code for self-energies, tadpole corrections, renormalization group equations (RGEs) and electroweak symmetry breaking (EWSB) conditions and combines them with numerical routines for solving the RGEs and EWSB conditions simultaneously. The resulting spectrum generator is then able to solve for the spectrum of the model, including loop-corrected pole masses, consistent with user specified boundary conditions. The modular structure of the generated code allows for individual components to be replaced with an alternative if available. FlexibleSUSY has been carefully designed to grow as alternative solvers and calculators are added. Predefined models include the MSSM, NMSSM, E6SSM, USSM, R-symmetric models and models with right-handed neutrinos. Program Summary Program title: FlexibleSUSY Catalogue identifier: AEVIv10 Program summary URL: http://cpc.cs.qub.ac.uk/summaries/AEVIv10.html obtainable from: CPC Program Library, Queen's University, Belfast, N. Ireland Licensing provisions: GNU General Public License, version 3 No. of lines in distributed program, including test data, etc.: 129406 No. of bytes in distributed program, including test data, etc.: 854831 Distribution format: tar.gz Programming language: C++, Wolfram/Mathematica, FORTRAN, Bourne shell. Computer: Personal computer. Operating system: Tested on Linux 3.x, Mac OS X. Classification: 11.1, 11.6, 6.5. External routines: SARAH 4.0.4, Boost library, Eigen, LAPACK Nature of problem: Determining the mass spectrum and mixings for any supersymmetric model. The generated code must find simultaneous solutions to constraints which are specified at two or more different renormalization scales, which are connected by renormalization group equations forming a large set of coupled first-order differential equations. Solution method: Nested iterative algorithm and numerical minimization of the Higgs potential. Restrictions: The couplings must remain perturbative at all scales between the highest and the lowest boundary condition. FlexibleSUSY assumes that all couplings of the model are real (i.e. CP-conserving). Due to the modular nature of the generated code, adaption and extension to overcome restrictions in scope is quite straightforward. Running time: 0.06-0.2 seconds per parameter point.
|
|
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.
|
|
Bonilla, J., Brivio, I., Gavela, M. B., & Sanz, V. (2021). One-loop corrections to ALP couplings. J. High Energy Phys., 11(11), 168–57pp.
Abstract: The plethora of increasingly precise experiments which hunt for axion-like particles (ALPs), as well as their widely different energy reach, call for the theoretical understanding of ALP couplings at loop-level. We derive the one-loop contributions to ALP-SM effective couplings, including finite corrections. The complete leading-order – dimension five – effective linear Lagrangian is considered. The ALP is left off-shell, which is of particular impact on LHC and accelerator searches of ALP couplings to gamma gamma, ZZ, WW, Z gamma gluons and fermions. All results are obtained in the covariant Rg gauge. A few phenomenological consequences are also explored as illustration, with flavour diagonal channels in the case of fermions: in particular, we explore constraints on the coupling of the ALP to top quarks, that can be extracted from LHC data, from astrophysical sources and from Dark Matter direct detection experiments such as PandaX, LUX and XENONIT. Furthermore, we clarify the relation between alternative ALP bases, the role of gauge anomalous couplings and their interface with chirality-conserving and chirality-flip fermion interactions, and we briefly discuss renormalization group aspects.
|
|
Bustamante, M., Gago, A. M., & Jones Perez, J. (2011). SUSY renormalization group effects in ultra high energy neutrinos. J. High Energy Phys., 05(5), 133–26pp.
Abstract: We have explored the question of whether the renormalization group running of the neutrino mixing parameters in the Minimal Supersymmetric Standard Model is detectable with ultra-high energy neutrinos from active galactic nuclei (AGN). We use as observables the ratios of neutrino fluxes produced at the AGN, focusing on four different neutrino production models: (Phi(0)(v epsilon+(v) over bar epsilon) : Phi(0)(v mu+(v) over bar mu) : Phi(0)(v tau+(v) over bar tau)) = (1 : 2 : 0), (0 : 1 : 0), (1 : 0 : 0), and (1 : 1 : 0). The prospects for observing deviations experimentally are taken into consideration, and we find out that it is necessary to impose a cut-off on the transferred momentum of Q(2) >= 10(7) GeV(2). However, this condition, together with the expected low value of the diffuse AGN neutrino flux, yields a negligible event rate at a km-scale. Cerenkov detector such as IceCube.
|
|
Dhani, P. K., Rodrigo, G., & Sborlini, G. F. R. (2023). Triple-collinear splittings with massive particles. J. High Energy Phys., 12(12), 188–20pp.
Abstract: We analyze in detail the most singular behaviour of processes involving triple-collinear splittings with massive particles in the quasi-collinear limit, and present compact expressions for the splitting amplitudes and the corresponding splitting kernels at the squared-amplitude level. Our expressions fully agree with well-known triple-collinear splittings in the massless limit, which are used as a guide to achieve the final expressions. These results are important to quantify dominant mass effects in many observables, and constitute an essential ingredient of current high-precision computational frameworks for collider phenomenology.
|
|
Escribano, P., Reig, M., & Vicente, A. (2020). Generalizing the Scotogenic model. J. High Energy Phys., 07(7), 097–25pp.
Abstract: The Scotogenic model is an economical setup that induces Majorana neutrino masses at the 1-loop level and includes a dark matter candidate. We discuss a generalization of the original Scotogenic model with arbitrary numbers of generations of singlet fermion and inert doublet scalar fields. First, the full form of the light neutrino mass matrix is presented, with some comments on its derivation and with special attention to some particular cases. The behavior of the theory at high energies is explored by solving the Renormalization Group Equations.
|
|
Giachino, A., van Hameren, A., & Ziarko, G. (2024). A new subtraction scheme at NLO exploiting the privilege of k<sub>T</sub>-factorization. J. High Energy Phys., 06(6), 167–39pp.
Abstract: We present a subtraction method for the calculation of real-radiation integrals at NLO in hybrid k(T)-factorization. The main difference with existing methods for collinear factorization is that we subtract the momentum recoil, occurring due to the mapping from an (n + 1)-particle phase space to an n-particle phase space, from the initial-state momenta, instead of distributing it over the final-state momenta.
|
|
Mehrabankar, S., Garcia-March, M. A., Almudever, C. G., & Perez, A. (2024). Reducing the number of qubits in quantum simulations of one dimensional many-body Hamiltonians. New J. Phys., 26(8), 083023–17pp.
Abstract: We investigate the Ising and Heisenberg models using the block renormalization group method (BRGM), focusing on its behavior across different system sizes. The BRGM reduces the number of spins by a factor of 1/2 (1/3) for the Ising (Heisenberg) model, effectively preserving essential physical features of the model while using only a fraction of the spins. Through a comparative analysis, we demonstrate that as the system size increases, there is an exponential convergence between results obtained from the original and renormalized Ising Hamiltonians, provided the coupling constants are redefined accordingly. Remarkably, for a spin chain with 24 spins, all physical features, including magnetization, correlation function, and entanglement entropy, exhibit an exact correspondence with the results from the original Hamiltonian. The study of the Heisenberg model also shows this tendency, although complete convergence may appear for a size much larger than 24 spins, and is therefore beyond our computational capabilities. The success of BRGM in accurately characterizing the Ising model, even with a relatively small number of spins, underscores its robustness and utility in studying complex physical systems, and facilitates its simulation on current NISQ computers, where the available number of qubits is largely constrained.
|
|
Merle, A., Platscher, M., Rojas, N., Valle, J. W. F., & Vicente, A. (2016). Consistency of WIMP Dark Matter as radiative neutrino mass messenger. J. High Energy Phys., 07(7), 013–17pp.
Abstract: The scotogenic scenario provides an attractive approach to both Dark Matter and neutrino mass generation, in which the same symmetry that stabilises Dark Matter also ensures the radiative seesaw origin of neutrino mass. However the simplest scenario may suffer from inconsistencies arising from the spontaneous breaking of the underlying Z(2) symmetry. Here we show that the singlet-triplet extension of the simplest model naturally avoids this problem due to the presence of scalar triplets neutral under the Z(2) which affect the evolution of the couplings in the scalar sector. The scenario offers good prospects for direct WIMP Dark Matter detection through the nuclear recoil method.
|