|
Otten, S., Rolbiecki, K., Caron, S., Kim, J. S., Ruiz de Austri, R., & Tattersall, J. (2020). DeepXS: fast approximation of MSSM electroweak cross sections at NLO. Eur. Phys. J. C, 80(1), 12–9pp.
Abstract: We present a deep learning solution to the prediction of particle production cross sections over a complicated, high-dimensional parameter space. We demonstrate the applicability by providing state-of-the-art predictions for the production of charginos and neutralinos at the Large Hadron Collider (LHC) at the next-to-leading order in the phenomenological MSSM-19 and explicitly demonstrate the performance for pp ->(chi) over tilde (+)(1)(chi) over tilde (-)(1), (chi) over tilde (0)(2)(chi) over tilde (0)(2) and (chi) over tilde (0)(2)(chi) over tilde (+/-)(1) as a proof of concept which will be extended to all SUSY electroweak pairs. We obtain errors that are lower than the uncertainty from scale and parton distribution functions with mean absolute percentage errors of well below 0.5% allowing a safe inference at the next-to-leading order with inference times that improve the Monte Carlo integration procedures that have been available so far by a factor of O(10(7)) from O(min) to O(mu s) per evaluation.
|
|
|
Kpatcha, E., Lopez-Fogliani, D. E., Munoz, C., & Ruiz de Austri, R. (2020). Impact of Higgs physics on the parameter space of the μnu SSM. Eur. Phys. J. C, 80(4), 336–43pp.
Abstract: Given the increasing number of experimental data, together with the precise measurement of the properties of the Higgs boson at the LHC, the parameter space of supersymmetric models starts to be constrained. We carry out a detailed analysis of this issue in the framework of the μnu SSM. In this model, three families of right-handed neutrino superfields are present in order to solve the μproblem and simultaneously reproduce neutrino physics. The new couplings and sneutrino vacuum expectation values in the μnu SSM induce new mixing of states, and, in particular, the three right sneutrinos can be substantially mixed with the neutral Higgses. After diagonalization, the masses of the corresponding three singlet-like eigenstates can be smaller or larger than the mass of the Higgs, or even degenerated with it. We analyze whether these situations are still compatible with the experimental results. To address it we scan the parameter space of the Higgs sector of the model. In particular, we sample the μnu SSM using a powerful likelihood data-driven method, paying special attention to satisfy the constraints coming from Higgs sector measurements/limits (using HiggsBounds and HiggsSignals), as well as a class of flavor observables such as B and μdecays, while muon g-2 is briefly discussed. We find that large regions of the parameter space of the μnu SSM are viable, containing an interesting phenomenology that could be probed at the LHC.
|
|
|
Felea, D., Mamuzic, J., Maselek, R., Mavromatos, N. E., Mitsou, V. A., Pinfold, J. L., et al. (2020). Prospects for discovering supersymmetric long-lived particles with MoEDAL. Eur. Phys. J. C, 80(5), 431–12pp.
Abstract: We present a study on the possibility of searching for long-lived supersymmetric partners with the MoEDAL experiment at the LHC. MoEDAL is sensitive to highly ionising objects such as magnetic monopoles or massive (meta)stable electrically charged particles. We focus on prospects of directly detecting long-lived sleptons in a phenomenologically realistic model which involves an intermediate neutral long-lived particle in the decay chain. This scenario is not yet excluded by the current data from ATLAS or CMS, and is compatible with astrophysical constraints. Using Monte Carlo simulation, we compare the sensitivities of MoEDAL versus ATLAS in scenarios where MoEDAL could provide discovery reach complementary to ATLAS and CMS, thanks to looser selection criteria combined with the virtual absence of background. It is also interesting to point out that, in such scenarios, in which charged staus are the main long-lived candidates, the relevant mass range for MoEDAL is compatible with a potential role of Supersymmetry in providing an explanation for the anomalous events observed by the ANITA detector.
|
|
|
Desai, N., Domingo, F., Kim, J. S., Ruiz de Austri, R., Rolbiecki, K., Sonawane, M., et al. (2021). Constraining electroweak and strongly charged long-lived particles with CheckMATE. Eur. Phys. J. C, 81(11), 968–19pp.
Abstract: Long-lived particles have become a new frontier in the exploration of physics beyond the Standard Model. In this paper, we present the implementation of four types of long-lived particle searches, viz. displaced leptons, disappearing track, displaced vertex with either muons or with missing transverse energy, and heavy charged tracks. These four categories cover the signatures of a large range of physics models. We illustrate their potential for exclusion and discuss their mutual overlaps in mass-lifetime space for two simple phenomenological models involving either a U(1)-charged or a coloured scalar.
|
|
|
MoEDAL Collaboration(Acharya, B. et al), Musumeci, E., Mitsou, V. A., Papavassiliou, J., Ruiz de Austri, R., Santra, A., et al. (2022). Search for highly-ionizing particles in pp collisions at the LHC's Run-1 using the prototype MoEDAL detector. Eur. Phys. J. C, 82(8), 694–16pp.
Abstract: A search for highly electrically charged objects (HECOs) and magnetic monopoles is presented using 2.2 fb(-1) of p – p collision data taken at a centre of mass energy (E-CM) of 8 TeV by the MoEDAL detector during LHC's Run-1. The data were collected using MoEDAL's prototype Nuclear Track Detectord array and the Trapping Detector array. The results are interpreted in terms of Drell-Yan pair production of stable HECO and monopole pairs with three spin hypotheses (0, 1/2 and 1). The search provides constraints on the direct production of magnetic monopoles carrying one to four Dirac magnetic charges and with mass limits ranging from 590 GeV/c(2) to 1 TeV/c(2). Additionally, mass limits are placed on HECOs with charge in the range 10e to 180e, where e is the charge of an electron, for masses between 30 GeV/c(2) and 1 TeV/c(2).
|
|
|
Allanach, B. C., Bednyakov, A., & Ruiz de Austri, R. (2015). Higher order corrections and unification in the minimal supersymmetric standard model: SOFTSUSY3.5. Comput. Phys. Commun., 189, 192–206.
Abstract: We explore the effects of three-loop minimal supersymmetric standard model renormalisation group equation terms and some leading two-loop threshold corrections on gauge and Yukawa unification: each being one loop higher order than current public spectrum calculators. We also explore the effect of the higher order terms (often 2-3 GeV) on the lightest CP even Higgs mass prediction. We illustrate our results in the constrained minimal supersymmetric standard model. Neglecting threshold corrections at the grand unified scale, the discrepancy between the unification scale alpha(s) and the other two unified gauge couplings changes by 0.1% due to the higher order corrections and the difference between unification scale bottom-tau Yukawa couplings neglecting unification scale threshold corrections changes by up to 1%. The difference between unification scale bottom and top Yukawa couplings changes by a few percent. Differences due to the higher order corrections also give an estimate of the size of theoretical uncertainties in the minimal supersymmetric standard model spectrum. We use these to provide estimates of theoretical uncertainties in predictions of the dark matter relic density (which can be of order one due to its strong dependence on sparticle masses) and the LHC sparticle production cross-section (often around 30%). The additional higher order corrections have been incorporated into SOFTSUSY, and we provide details on how to compile and use the program. We also provide a summary of the approximations used in the higher order corrections. Program Summary Nature of problem: Calculating supersymmetric particle spectrum and mixing parameters in the minimal supersymmetric standard model. The solution to the renormalisation group equations must be consistent with boundary conditions on supersymmetry breaking parameters, as well as the weak-scale boundary condition on gauge couplings, Yukawa couplings and the Higgs potential parameters. Program title: SOFTSUSY Catalogue identifier: ADPMv50 Program summary URL: http://cpc.cs.qub.ac.uk/summaries/ADPMv50.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.: 240528 No. of bytes in distributed program, including test data, etc.: 2597933 Distribution format: tar.gz Programming language: C++, Fortran. Computer: Personal computer. Operating system: Tested on Linux 3.4.6. Word size: 64 bits. Classification: 11.1, 11.6. External routines: At least GiNaC1.3.5 [1] and CLN1.3.1 (both freely obtainable from http://www.ginac.de). Does the new version supersede the previous version?: Yes Catalogue identifier of previous version: ADPMv40 Journal reference of previous version: Comput. Phys. Comm. 185 (2014) 2322 Solution method: Nested iterative algorithm. Reasons for new version: Extension to include additional two and three-loop terms. Summary of revisions: All quantities in the minimal supersymmetric standard model are extended to have three-loop renormalisation group equations (including 3-family mixing) in the limit of real parameters and some leading two-loop threshold corrections are incorporated to the third family Yukawa couplings and the strong gauge coupling. 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 non-physical 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 real R-parity conserving minimal supersymmetric standard model (MSSM) only. Running time: A minute per parameter point. The tests provided with the package only take a few seconds to run.
|
|
|
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/
|
|
|
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.
|
|
|
Trotta, R., Johannesson, G., Moskalenko, I. V., Porter, T. A., Ruiz de Austri, R., & Strong, A. W. (2011). Constraints on Cosmic-Ray Propagation Models from a Global Bayesian Analysis. Astrophys. J., 729(2), 106–16pp.
Abstract: Research in many areas of modern physics such as, e. g., indirect searches for dark matter and particle acceleration in supernova remnant shocks rely heavily on studies of cosmic rays (CRs) and associated diffuse emissions (radio, microwave, X-rays, gamma-rays). While very detailed numerical models of CR propagation exist, a quantitative statistical analysis of such models has been so far hampered by the large computational effort that those models require. Although statistical analyses have been carried out before using semi-analytical models (where the computation is much faster), the evaluation of the results obtained from such models is difficult, as they necessarily suffer from many simplifying assumptions. The main objective of this paper is to present a working method for a full Bayesian parameter estimation for a numerical CR propagation model. For this study, we use the GALPROP code, the most advanced of its kind, which uses astrophysical information, and nuclear and particle data as inputs to self-consistently predict CRs, gamma-rays, synchrotron, and other observables. We demonstrate that a full Bayesian analysis is possible using nested sampling and Markov Chain Monte Carlo methods (implemented in the SuperBayeS code) despite the heavy computational demands of a numerical propagation code. The best-fit values of parameters found in this analysis are in agreement with previous, significantly simpler, studies also based on GALPROP.
|
|
|
Johannesson, G., Ruiz de Austri, R., Vincent, A. C., Moskalenko, I. V., Orlando, E., Porter, T. A., et al. (2016). Bayesian analysis of cosmic-ray propagation: evidence against homogeneous diffusion. Astrophys. J., 824(1), 16–19pp.
Abstract: We present the results of the most complete scan of the parameter space for cosmic ray (CR) injection and propagation. We perform a Bayesian search of the main GALPROP parameters, using the MultiNest nested sampling algorithm, augmented by the BAMBI neural network machine-learning package. This is the first study to separate out low-mass isotopes (p, (p) over bar and He) from the usual light elements (Be, B, C, N, and O). We find that the propagation parameters that best-fit p, (p) over bar, and He data are significantly different from those that fit light elements, including the B/C and Be-10/Be-9 secondary-to-primary ratios normally used to calibrate propagation parameters. This suggests that each set of species is probing a very different interstellar medium, and that the standard approach of calibrating propagation parameters using B/C can lead to incorrect results. We present posterior distributions and best-fit parameters for propagation of both sets of nuclei, as well as for the injection abundances of elements from H to Si. The input GALDEF files with these new parameters will be included in an upcoming public GALPROP update.
|
|