Balazs, C. et al, Mamuzic, J., & Ruiz de Austri, R. (2021). A comparison of optimisation algorithms for high-dimensional particle and astrophysics applications. J. High Energy Phys., 05(5), 108–46pp.
Abstract: Optimisation problems are ubiquitous in particle and astrophysics, and involve locating the optimum of a complicated function of many parameters that may be computationally expensive to evaluate. We describe a number of global optimisation algorithms that are not yet widely used in particle astrophysics, benchmark them against random sampling and existing techniques, and perform a detailed comparison of their performance on a range of test functions. These include four analytic test functions of varying dimensionality, and a realistic example derived from a recent global fit of weak-scale supersymmetry. Although the best algorithm to use depends on the function being investigated, we are able to present general conclusions about the relative merits of random sampling, Differential Evolution, Particle Swarm Optimisation, the Covariance Matrix Adaptation Evolution Strategy, Bayesian Optimisation, Grey Wolf Optimisation, and the PyGMO Artificial Bee Colony, Gaussian Particle Filter and Adaptive Memory Programming for Global Optimisation algorithms.
|
Caron, S., Ruiz de Austri, R., & Zhang, Z. Y. (2023). Mixture-of-Theories training: can we find new physics and anomalies better by mixing physical theories? J. High Energy Phys., 03(3), 004–37pp.
Abstract: Model-independent search strategies have been increasingly proposed in recent years because on the one hand there has been no clear signal for new physics and on the other hand there is a lack of a highly probable and parameter-free extension of the standard model. For these reasons, there is no simple search target so far. In this work, we try to take a new direction and ask the question: bearing in mind that we have a large number of new physics theories that go beyond the Standard Model and may contain a grain of truth, can we improve our search strategy for unknown signals by using them “in combination”? In particular, we show that a signal hypothesis based on a large, intermingled set of many different theoretical signal models can be a superior approach to find an unknown BSM signal. Applied to a recent data challenge, we show that “mixture-of-theories training” outperforms strategies that optimize signal regions with a single BSM model as well as most unsupervised strategies. Applications of this work include anomaly detection and the definition of signal regions in the search for signals of new physics.
|
Jueid, A., Kip, J., Ruiz de Austri, R., & Skands, P. (2024). The Strong Force meets the Dark Sector: a robust estimate of QCD uncertainties for anti-matter dark matter searches. J. High Energy Phys., 02(2), 119–48pp.
Abstract: In dark-matter annihilation channels to hadronic final states, stable particles – such as positrons, photons, antiprotons, and antineutrinos – are produced via complex sequences of phenomena including QED/QCD radiation, hadronisation, and hadron decays. These processes are normally modelled by Monte Carlo (MC) event generators whose limited accuracy imply intrinsic QCD uncertainties on the predictions for indirect-detection experiments like Fermi-LAT, Pamela, IceCube or Ams-02. In this article, we perform a comprehensive analysis of QCD uncertainties, meaning both perturbative and nonperturbative sources of uncertainty are included – estimated via variations of MC renormalization-scale and fragmentation-function parameters, respectively – in antimatter spectra from dark-matter annihilation, based on parametric variations of the Pythia 8 event generator. After performing several retunings of light-quark fragmentation functions, we define a set of variations that span a conservative estimate of the QCD uncertainties. We estimate the effects on antimatter spectra for various annihilation channels and final-state particle species, and discuss their impact on fitted values for the dark-matter mass and thermally-averaged annihilation cross section. We find dramatic impacts which can go up to O(10%) for the annihilation cross section. We provide the spectra in tabulated form including QCD uncertainties and code snippets to perform fast dark-matter fits, in this github repository.
|
Kim, J. S., Lopez-Fogliani, D. E., Perez, A. D., & Ruiz de Austri, R. (2022). The new (g-2)(mu) and right-handed sneutrino dark matter. Nucl. Phys. B, 974, 115637–23pp.
Abstract: In this paper we investigate the (g – 2)(mu) discrepancy in the context of the R-parity conserving next-to minimal supersymmetric Standard Model plus right-handed neutrinos superfields. The model has the ability to reproduce neutrino physics data and includes the interesting possibility to have the right-handed sneutrino as the lightest supersymmetric particle and a viable dark matter candidate. Since right-handed sneutrinos are singlets, no new contributions for delta a(mu) with respect to the MSSM and NMSSM are present. However, the possibility to have the right-handed sneutrino as the lightest supersymmetric particle opens new ways to escape Large Hadron Collider and direct detection constraints. In particular, we find that dark matter masses within 10 less than or similar to m((upsilon) over tildeR) less than or similar to 600 GeV are fully compatible with current experimental constraints. Remarkably, not only spectra with light sleptons are needed, but we obtain solutions with m((mu) over tilde) greater than or similar to 600 GeV in the entire dark matter mass range that could be probed by new (g – 2)(mu) data in the near future. In addition, dark matter direct detection experiments will be able to explore a sizable portion of the allowed parameter space with mvR < 300 GeV, while indirect detection experiments will be able to probe a much smaller fraction within 200 less than or similar to m((nu)over tilde>R) less than or similar to 350 GeV.
|
Kim, J. S., Reuter, J., Rolbiecki, K., & Ruiz de Austri, R. (2016). A resonance without resonance: Scrutinizing the diphoton excess at 750 GeV. Phys. Lett. B, 755, 403–408.
Abstract: Motivated by the recent diphoton excesses reported by both ATLAS and CMS collaborations, we suggest that a new heavy spinless particle is produced in gluon fusion at the LHC and decays to a couple of lighter pseudoscalars which then decay to photons. The new resonances could arise from a new strongly interacting sector and couple to Standard Model gauge bosons only via the corresponding Wess-Zumino-Witten anomaly. We present a detailed recast of the newest 13 TeV data from ATLAS and CMS together with the 8 TeV data to scan the consistency of the parameter space for those resonances.
|
MoEDAL Collaboration(Acharya, B. et al), Bernabeu, J., Garcia, C., Mamuzic, J., Mitsou, V. A., Ruiz de Austri, R., et al. (2018). Search for magnetic monopoles with the MoEDAL forward trapping detector in 2.11 fb(-1) of 13 TeV proton-proton collisions at the LHC. Phys. Lett. B, 782, 510–516.
Abstract: We update our previous search for trapped magnetic monopoles in LHC Run 2 using nearly six times more integrated luminosity and including additional models for the interpretation of the data. The MoEDAL forward trapping detector, comprising 222 kg of aluminium samples, was exposed to 2.11 fb(-1) of 13 TeV proton-proton collisions near the LHCb interaction point and analysed by searching for induced persistent currents after passage through a superconducting magnetometer. Magnetic charges equal to the Dirac charge or above are excluded in all samples. The results are interpreted in Drell-Yan production models for monopoles with spins 0, 1/2 and 1: in addition to standard point-like couplings, we also consider couplings with momentum-dependent form factors. The search provides the best current laboratory constraints for monopoles with magnetic charges ranging from two to five times the Dirac charge.
|
Ghosh, P., Lara, I., Lopez-Fogliani, D. E., Muñoz, C., & Ruiz de Austri, R. (2018). Searching for left sneutrino LSP at the LHC. Int. J. Mod. Phys. A, 33(18-19), 1850110–62pp.
Abstract: We analyze relevant signals expected at the LHC for a left sneutrino as the lightest supersymmetric particle (LSP). The discussion is carried out in the “mu from nu” supersymmetric standard model (mu nu SSM), where the presence of R-parity breaking couplings involving right-handed neutrinos solves the μproblem and reproduces neutrino data. The sneutrinos are pair produced via a virtual W, Z or gamma in the s channel. From the prompt decay of a pair of left sneutrinos LSPs of any family, a significant diphoton signal plus missing transverse energy (MET) from neutrinos can be present in the mass range 118-132 GeV, with 13 TeV center-of-mass energy and an integrated luminosity of 100 fb(-1). In addition, in the case of a pair of tau left sneutrinos LSPs, given the large value of the tau Yukawa coupling diphoton plus leptons and/or multileptons can appear. We find that the number of expected events for the multilepton signal, together with properly adopted search strategies, is sufficient to give a significant evidence for a sneutrino of mass in the range 130-310 GeV, even with the integrated luminosity of 20 fb(-1). In the case of the signal producing diphoton plus leptons, an integrated luminosity of 100 fb(-1) is needed to give a significant evidence in the mass range 95-145 GeV. Finally, we discuss briefly the presence of displaced vertices and the associated range of masses.
|
Flores, M. M., Kim, J. S., Rolbiecki, K., & Ruiz de Austri, R. (2023). Updated LHC bounds on MUED after run 2. Int. J. Mod. Phys. A, 38(1), 2350002–14pp.
Abstract: We present updated LHC limits on the minimal universal extra dimensions (MUEDs) model from the Run 2 searches. We scan the parameter space against a number of searches implemented in the public code CheckMATE and derive up-to-date limits on the MUED parameter space from 13TeV searches. The strongest constraints come from a search dedicated to squarks and gluinos with one isolated lepton, jets and missing transverse energy. In the procedure, we take into account initial state radiation and stress its importance in the MUED searches, which is not always appreciated.
|
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.
|
MoEDAL Collaboration(Acharya, B. et al), Bernabeu, J., Garcia, C., Mamuzic, J., Mitsou, V. A., Ruiz de Austri, R., et al. (2017). Search for Magnetic Monopoles with the MoEDAL Forward Trapping Detector in 13 TeV Proton-Proton Collisions at the LHC. Phys. Rev. Lett., 118(6), 061801–6pp.
Abstract: MoEDAL is designed to identify new physics in the form of long-lived highly ionizing particles produced in high-energy LHC collisions. Its arrays of plastic nuclear-track detectors and aluminium trapping volumes provide two independent passive detection techniques. We present here the results of a first search for magnetic monopole production in 13 TeV proton-proton collisions using the trapping technique, extending a previous publication with 8 TeV data during LHC Run 1. A total of 222 kg of MoEDAL trapping detector samples was exposed in the forward region and analyzed by searching for induced persistent currents after passage through a superconducting magnetometer. Magnetic charges exceeding half the Dirac charge are excluded in all samples and limits are placed for the first time on the production of magnetic monopoles in 13 TeV pp collisions. The search probes mass ranges previously inaccessible to collider experiments for up to five times the Dirac charge.
|