Showing 1–50 of 57 results for author: de Austri, R R

Insights into Dark Matter Direct Detection Experiments: Decision Trees versus Deep Learning

Authors: Daniel E. Lopez-Fogliani, Andres D. Perez, Roberto Ruiz de Austri

Abstract: The detection of Dark Matter (DM) remains a significant challenge in particle physics. This study exploits advanced machine learning models to improve detection capabilities of liquid xenon time projection chamber experiments, utilizing state-of-the-art transformers alongside traditional methods like Multilayer Perceptrons and Convolutional Neural Networks. We evaluate various data representations… ▽ More The detection of Dark Matter (DM) remains a significant challenge in particle physics. This study exploits advanced machine learning models to improve detection capabilities of liquid xenon time projection chamber experiments, utilizing state-of-the-art transformers alongside traditional methods like Multilayer Perceptrons and Convolutional Neural Networks. We evaluate various data representations and find that simplified feature representations, particularly corrected S1 and S2 signals, retain critical information for classification. Our results show that while transformers offer promising performance, simpler models like XGBoost can achieve comparable results with optimal data representations. We also derive exclusion limits in the cross-section versus DM mass parameter space, showing minimal differences between XGBoost and the best performing deep learning models. The comparative analysis of different machine learning approaches provides a valuable reference for future experiments by guiding the choice of models and data representations to maximize detection capabilities. △ Less

Submitted 14 June, 2024; originally announced June 2024.

Comments: 26 pages, 7 figures, 2 tables

Report number: IFT-UAM/CSIC-24-89

CosmiXs: Cosmic messenger spectra for indirect dark matter searches

Authors: Chiara Arina, Mattia Di Mauro, Nicolao Fornengo, Jan Heisig, Adil Jueid, Roberto Ruiz de Austri

Abstract: The energy spectra of particles produced from dark matter (DM) annihilation or decay are one of the fundamental ingredients to calculate the predicted fluxes of cosmic rays and radiation searched for in indirect DM detection. We revisit the calculation of the source spectra for annihilating and decaying DM using the Vincia shower algorithm in Pythia to include QED and QCD final state radiation and… ▽ More The energy spectra of particles produced from dark matter (DM) annihilation or decay are one of the fundamental ingredients to calculate the predicted fluxes of cosmic rays and radiation searched for in indirect DM detection. We revisit the calculation of the source spectra for annihilating and decaying DM using the Vincia shower algorithm in Pythia to include QED and QCD final state radiation and diagrams for the Electroweak (EW) corrections with massive bosons, not present in the default Pythia shower model. We take into account the spin information of the particles during the entire EW shower and the off-shell contributions from massive gauge bosons. Furthermore, we perform a dedicated tuning of the Vincia and Pythia parameters to LEP data on the production of pions, photons, and hyperons at the $Z$ resonance and discuss the underlying uncertainties. To enable the use of our results in DM studies, we provide the tabulated source spectra for the most relevant cosmic messenger particles, namely antiprotons, positrons, $γ$ rays and the three neutrino flavors, for all the fermionic and bosonic channels and DM masses between 5 GeV and 100 TeV, on https://github.com/ajueid/CosmiXs.git. △ Less

Submitted 2 December, 2023; originally announced December 2023.

Comments: 39 pages, 14 figures, 4 tables

Report number: TTK-23-32, CTPU-PTC-23-36

AutoSourceID-Classifier. Star-Galaxy Classification using a Convolutional Neural Network with Spatial Information

Authors: F. Stoppa, S. Bhattacharyya, R. Ruiz de Austri, P. Vreeswijk, S. Caron, G. Zaharijas, S. Bloemen, G. Principe, D. Malyshev, V. Vodeb, P. J. Groot, E. Cator, G. Nelemans

Abstract: Aims. Traditional star-galaxy classification techniques often rely on feature estimation from catalogues, a process susceptible to introducing inaccuracies, thereby potentially jeopardizing the classification's reliability. Certain galaxies, especially those not manifesting as extended sources, can be misclassified when their shape parameters and flux solely drive the inference. We aim to create a… ▽ More Aims. Traditional star-galaxy classification techniques often rely on feature estimation from catalogues, a process susceptible to introducing inaccuracies, thereby potentially jeopardizing the classification's reliability. Certain galaxies, especially those not manifesting as extended sources, can be misclassified when their shape parameters and flux solely drive the inference. We aim to create a robust and accurate classification network for identifying stars and galaxies directly from astronomical images. By leveraging convolutional neural networks (CNN) and additional information about the source position, we aim to accurately classify all stars and galaxies within a survey, particularly those with a signal-to-noise ratio (S/N) near the detection limit. Methods. The AutoSourceID-Classifier (ASID-C) algorithm developed here uses 32x32 pixel single filter band source cutouts generated by the previously developed ASID-L code. ASID-C utilizes CNNs to distinguish these cutouts into stars or galaxies, leveraging their strong feature-learning capabilities. Subsequently, we employ a modified Platt Scaling calibration for the output of the CNN. This technique ensures that the derived probabilities are effectively calibrated, delivering precise and reliable results. Results. We show that ASID-C, trained on MeerLICHT telescope images and using the Dark Energy Camera Legacy Survey (DECaLS) morphological classification, outperforms similar codes like SourceExtractor. ASID-C opens up new possibilities for accurate celestial object classification, especially for sources with a S/N near the detection limit. Potential applications of ASID-C, like real-time star-galaxy classification and transient's host identification, promise significant contributions to astronomical research. △ Less

Submitted 26 July, 2023; originally announced July 2023.

Journal ref: A&A 680, A109 (2023)

AutoSourceID-FeatureExtractor. Optical image analysis using a two-step mean variance estimation network for feature estimation and uncertainty characterisation

Authors: F. Stoppa, R. Ruiz de Austri, P. Vreeswijk, S. Bhattacharyya, S. Caron, S. Bloemen, G. Zaharijas, G. Principe, V. Vodeb, P. J. Groot, E. Cator, G. Nelemans

Abstract: Aims. In astronomy, machine learning has been successful in various tasks such as source localisation, classification, anomaly detection, and segmentation. However, feature regression remains an area with room for improvement. We aim to design a network that can accurately estimate sources' features and their uncertainties from single-band image cutouts, given the approximated locations of the sou… ▽ More Aims. In astronomy, machine learning has been successful in various tasks such as source localisation, classification, anomaly detection, and segmentation. However, feature regression remains an area with room for improvement. We aim to design a network that can accurately estimate sources' features and their uncertainties from single-band image cutouts, given the approximated locations of the sources provided by the previously developed code AutoSourceID-Light (ASID-L) or other external catalogues. This work serves as a proof of concept, showing the potential of machine learning in estimating astronomical features when trained on meticulously crafted synthetic images and subsequently applied to real astronomical data. Methods. The algorithm presented here, AutoSourceID-FeatureExtractor (ASID-FE), uses single-band cutouts of 32x32 pixels around the localised sources to estimate flux, sub-pixel centre coordinates, and their uncertainties. ASID-FE employs a two-step mean variance estimation (TS-MVE) approach to first estimate the features and then their uncertainties without the need for additional information, for example the point spread function (PSF). For this proof of concept, we generated a synthetic dataset comprising only point sources directly derived from real images, ensuring a controlled yet authentic testing environment. Results.We show that ASID-FE, trained on synthetic images derived from the MeerLICHT telescope, can predict more accurate features with respect to similar codes such as SourceExtractor and that the two-step method can estimate well-calibrated uncertainties that are better behaved compared to similar methods that use deep ensembles of simple MVE networks. Finally, we evaluate the model on real images from the MeerLICHT telescope and the Zwicky Transient Facility (ZTF) to test its transfer learning abilities. △ Less

Submitted 29 October, 2023; v1 submitted 23 May, 2023; originally announced May 2023.

Journal ref: A&A 680, A108 (2023)

Gradient-Annihilated PINNs for Solving Riemann Problems: Application to Relativistic Hydrodynamics

Authors: Antonio Ferrer-Sánchez, José D. Martín-Guerrero, Roberto Ruiz de Austri, Alejandro Torres-Forné, José A. Font

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 requires the model to partially ignore high-gradients in the physical variables, achieved by introducing a suitable weighting funct… ▽ More 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 requires 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 and how the activation functions of the neural model are dynamically accounted for. The performance of our GA-PINN model 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 our 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 error obtained with a second-order, central, shock-capturing scheme. In all problems investigated, the accuracy reached by our GA-PINN model is comparable to that obtained with a shock-capturing scheme and significantly higher than that achieved by a baseline PINN algorithm. 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 ones, 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 method discussed could be a valuable tool to model relativistic flows in astrophysics and particle physics, characterized by the prevalence of discontinuous solutions. △ Less

Submitted 19 May, 2023; v1 submitted 15 May, 2023; originally announced May 2023.

Comments: 25 pages, 16 figures

The Strong Force meets the Dark Sector: a robust estimate of QCD uncertainties for anti-matter dark matter searches

Authors: Adil Jueid, Jochem Kip, Roberto Ruiz de Austri, Peter Skands

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 event generators whose limited accuracy imply intrinsic QCD uncertainties on the prediction… ▽ More 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 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 complete analysis of QCD uncertainties 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 $\mathcal{O}(40)$ GeV for uncertainties on the dark-matter mass and up to $\mathcal{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 https://github.com/ajueid/qcd-dm.github.io.git repository. △ Less

Submitted 20 March, 2023; originally announced March 2023.

Comments: v1: 52 pages, 21 figures, and 9 tables

Report number: CTPU-PTC-23-08

Mind the gap: The discrepancy between simulation and reality drives interpretations of the Galactic Center Excess

Authors: Sascha Caron, Christopher Eckner, Luc Hendriks, Guðlaugur Jóhannesson, Roberto Ruiz de Austri, Gabrijela Zaharijas

Abstract: The Galactic Center Excess (GCE) in GeV gamma rays has been debated for over a decade, with the possibility that it might be due to dark matter annihilation or undetected point sources such as millisecond pulsars (MSPs). This study investigates how the gamma-ray emission model ($γ$EM) used in Galactic center analyses affects the interpretation of the GCE's nature. To address this issue, we constru… ▽ More The Galactic Center Excess (GCE) in GeV gamma rays has been debated for over a decade, with the possibility that it might be due to dark matter annihilation or undetected point sources such as millisecond pulsars (MSPs). This study investigates how the gamma-ray emission model ($γ$EM) used in Galactic center analyses affects the interpretation of the GCE's nature. To address this issue, we construct an ultra-fast and powerful inference pipeline based on convolutional Deep Ensemble Networks. We explore the two main competing hypotheses for the GCE using a set of $γ$EMs with increasing parametric freedom. We calculate the fractional contribution ($f_{\mathrm{src}}$) of a dim population of MSPs to the total luminosity of the GCE and analyze its dependence on the complexity of the $γ$EM. For the simplest $γ$EM, we obtain $f_{\mathrm{src}} = 0.10 \pm 0.07$, while the most complex model yields $f_{\mathrm{src}} = 0.79 \pm 0.24.$ In conclusion, we find that the statement about the nature of the GCE (dark matter or not) strongly depends on the assumed $γ$EM. The quoted results for $f_{\mathrm{src}}$ do not account for the additional uncertainty arising from the fact that the observed gamma-ray sky is out-of-distribution concerning the investigated $γ$EM iterations. We quantify the reality gap between our $γ$EMs using deep-learning-based One-Class Deep Support Vector Data Description networks, revealing that all employed $γ$EMs have gaps to reality. Our study casts doubt on the validity of previous conclusions regarding the GCE and dark matter, and underscores the urgent need to account for the reality gap and consider previously overlooked ''out of domain'' uncertainties in future interpretations. △ Less

Submitted 23 May, 2023; v1 submitted 17 November, 2022; originally announced November 2022.

Comments: 56 pages, 25 figures; comments welcome! Accepted for submission to JCAP; text coincides with the published version

Report number: LAPTH-068/22

Estimating QCD uncertainties on antiproton spectra from dark-matter annihilation

Authors: Adil Jueid, Jochem Kip, Roberto Ruiz de Austri, Peter Skands

Abstract: In this talk, we discuss the physics modeling of antiproton spectra arising from dark matter (DM) annihilation or decay in a model-independent manner. The modeling of antiproton spectra contains some intrinsic uncertainties related to QCD parton showers and hadronisation of baryons. We briefly assess the sources of these uncertainties and their impact on antiproton energy spectra for a few selecte… ▽ More In this talk, we discuss the physics modeling of antiproton spectra arising from dark matter (DM) annihilation or decay in a model-independent manner. The modeling of antiproton spectra contains some intrinsic uncertainties related to QCD parton showers and hadronisation of baryons. We briefly assess the sources of these uncertainties and their impact on antiproton energy spectra for a few selected DM scenarios. The results are provided in tabulated form for future analyses. △ Less

Submitted 11 March, 2022; originally announced March 2022.

Comments: 16 pages, 4 figures. Contribution for the talk presented at the Workshop on Computational Tools for High Energy Physics and Cosmology - CompTools2021, 22 - 26 November, 2021, Institut de Physique des 2 Infinis (IP2I), Lyon, France. Based on arXiv:2202.11546. Data is tabulated and can be found at Github: https://github.com/ajueid/qcd-dm.github.io.git

Report number: KIAS-Q22005

Impact of QCD uncertainties on antiproton spectra from dark-matter annihilation

Authors: Adil Jueid, Jochem Kip, Roberto Ruiz de Austri, Peter Skands

Abstract: Dark-matter particles that annihilate or decay can undergo complex sequences of processes, including strong and electromagnetic radiation, hadronisation, and hadron decays, before particles that are stable on astrophysical time scales are produced. Antiprotons produced in this way may leave footprints in experiments such as AMS--02. Several groups have reported an excess of events in the antiproto… ▽ More Dark-matter particles that annihilate or decay can undergo complex sequences of processes, including strong and electromagnetic radiation, hadronisation, and hadron decays, before particles that are stable on astrophysical time scales are produced. Antiprotons produced in this way may leave footprints in experiments such as AMS--02. Several groups have reported an excess of events in the antiproton flux in the rigidity range of $10$--$20$ GV. However, the theoretical modeling of baryon production is not straightforward and relies in part on phenomenological models in Monte Carlo event generators. In this work, we assess the impact of QCD uncertainties on the spectra of antiprotons from dark-matter annihilation. As a proof-of-principle, we show that for a two-parameter model that depends only on the thermally-averaged annihilation cross section ($\langle σv \rangle$) and the dark-matter mass ($M_χ$), QCD uncertainties can affect the best-fit mass by up to $\sim 14 \%$ (with large uncertainties for large DM masses), depending on the choice of $M_χ$ and the annihilation channel ($b\bar{b}$ or $W^+ W^-$), and $\langle σv \rangle$ by up to $\sim 10\%$. For comparison, changes to the underlying diffusion parameters are found to be within $1\%$--$5\%$, and the results are also quite resilient to the choice of cosmic-ray propagation model. These findings indicate that QCD uncertainties need to be included in future DM analyses. To facilitate full-fledged analyses, we provide the spectra in tabulated form including QCD uncertainties and code snippets to perform mass interpolations and quick DM fits. The code can be found in this \href{https://github.com/ajueid/qcd-dm.github.io.git}{github} repository. △ Less

Submitted 17 February, 2023; v1 submitted 23 February, 2022; originally announced February 2022.

Comments: v2: 8 pages + references, 3 figures, 2 tables. Spectra with QCD uncertainties can be found in https://github.com/ajueid/qcd-dm.github.io.git

Report number: CTPU-PTC-22-22

AutoSourceID-Light. Fast Optical Source Localization via U-Net and Laplacian of Gaussian

Authors: Fiorenzo Stoppa, Paul Vreeswijk, Steven Bloemen, Saptashwa Bhattacharyya, Sascha Caron, Guðlaugur Jóhannesson, Roberto Ruiz de Austri, Chris van den Oetelaar, Gabrijela Zaharijas, Paul. J. Groot, Eric Cator, Gijs Nelemans

Abstract: $\textbf{Aims}… ▽ More $\textbf{Aims}$. With the ever-increasing survey speed of optical wide-field telescopes and the importance of discovering transients when they are still young, rapid and reliable source localization is paramount. We present AutoSourceID-Light (ASID-L), an innovative framework that uses computer vision techniques that can naturally deal with large amounts of data and rapidly localize sources in optical images. $\textbf{Methods}$. We show that the AutoSourceID-Light algorithm based on U-shaped networks and enhanced with a Laplacian of Gaussian filter (Chen et al. 1987) enables outstanding performances in the localization of sources. A U-Net (Ronneberger et al. 2015) network discerns the sources in the images from many different artifacts and passes the result to a Laplacian of Gaussian filter that then estimates the exact location. $\textbf{Results}$. Application on optical images of the MeerLICHT telescope demonstrates the great speed and localization power of the method. We compare the results with the widely used SExtractor (Bertin & Arnouts 1996) and show the out-performances of our method. AutoSourceID-Light rapidly detects more sources not only in low and mid crowded fields, but particularly in areas with more than 150 sources per square arcminute. △ Less

Submitted 15 July, 2022; v1 submitted 1 February, 2022; originally announced February 2022.

Journal ref: A&A 662, A109 (2022)

Identification of point sources in gamma rays using U-shaped convolutional neural networks and a data challenge

Authors: Boris Panes, Christopher Eckner, Luc Hendriks, Sascha Caron, Klaas Dijkstra, Guðlaugur Jóhannesson, Roberto Ruiz de Austri, Gabrijela Zaharijas

Abstract: At GeV energies, the sky is dominated by the interstellar emission from the Galaxy. With limited statistics and spatial resolution, accurately separating point sources is therefore challenging. Here we present the first application of deep learning based algorithms to automatically detect and classify point sources from gamma-ray data. For concreteness we refer to this approach as AutoSourceID. To… ▽ More At GeV energies, the sky is dominated by the interstellar emission from the Galaxy. With limited statistics and spatial resolution, accurately separating point sources is therefore challenging. Here we present the first application of deep learning based algorithms to automatically detect and classify point sources from gamma-ray data. For concreteness we refer to this approach as AutoSourceID. To detect point sources, we utilized U-shaped convolutional networks for image segmentation and {\it k}-means for source clustering and localization. We also explored the Centroid-Net algorithm, which is designed to find and count objects. The training data are based on 9.5 years of exposure from The Fermi Large Area Telescope (Fermi-LAT) and we used source properties of active galactic nuclei (AGNs) and pulsars (PSRs) from the fourth Fermi-LAT source catalog (4FGL) in addition to several models of background interstellar emission. The results of the localization algorithm are fed into a classification neural network that is trained to separate the three general source classes (AGNs, PSRs, and FAKE sources). We compared our localization algorithms qualitatively with traditional methods and find them to have similar detection thresholds. We also demonstrate the robustness of our source localization algorithms to modifications in the interstellar emission models, which presents a clear advantage over traditional methods. The classification network is able to discriminate between the three classes with typical accuracy of $\sim$ 70%, as long as balanced data sets are used in classification training. In https://github.com/bapanes/AutoSourceID, we publish our training data sets and analysis scripts and invite the community to join the data challenge aimed to improve the localization and classification of gamma-ray point sources. △ Less

Submitted 25 October, 2021; v1 submitted 19 March, 2021; originally announced March 2021.

Comments: 19 pages (two-column format), 14 figures, accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 656, A62 (2021)

Simple and statistically sound recommendations for analysing physical theories

Authors: Shehu S. AbdusSalam, Fruzsina J. Agocs, Benjamin C. Allanach, Peter Athron, Csaba Balázs, Emanuele Bagnaschi, Philip Bechtle, Oliver Buchmueller, Ankit Beniwal, Jihyun Bhom, Sanjay Bloor, Torsten Bringmann, Andy Buckley, Anja Butter, José Eliel Camargo-Molina, Marcin Chrzaszcz, Jan Conrad, Jonathan M. Cornell, Matthias Danninger, Jorge de Blas, Albert De Roeck, Klaus Desch, Matthew Dolan, Herbert Dreiner, Otto Eberhardt , et al. (50 additional authors not shown)

Abstract: Physical theories that depend on many parameters or are tested against data from many different experiments pose unique challenges to statistical inference. Many models in particle physics, astrophysics and cosmology fall into one or both of these categories. These issues are often sidestepped with statistically unsound ad hoc methods, involving intersection of parameter intervals estimated by mul… ▽ More Physical theories that depend on many parameters or are tested against data from many different experiments pose unique challenges to statistical inference. Many models in particle physics, astrophysics and cosmology fall into one or both of these categories. These issues are often sidestepped with statistically unsound ad hoc methods, involving intersection of parameter intervals estimated by multiple experiments, and random or grid sampling of model parameters. Whilst these methods are easy to apply, they exhibit pathologies even in low-dimensional parameter spaces, and quickly become problematic to use and interpret in higher dimensions. In this article we give clear guidance for going beyond these procedures, suggesting where possible simple methods for performing statistically sound inference, and recommendations of readily-available software tools and standards that can assist in doing so. Our aim is to provide any physicists lacking comprehensive statistical training with recommendations for reaching correct scientific conclusions, with only a modest increase in analysis burden. Our examples can be reproduced with the code publicly available at https://doi.org/10.5281/zenodo.4322283. △ Less

Submitted 11 April, 2022; v1 submitted 17 December, 2020; originally announced December 2020.

Comments: 15 pages, 4 figures. extended discussions. closely matches version accepted for publication

Report number: PSI-PR-20-23, BONN-TH-2020-11, CP3-20-59, KCL-PH-TH/2020-75, P3H-20-080, TTP20-044, TUM-HEP-1310/20, IFT-UAM/CSIC-20-180, TTK-20-47, CERN-TH-2020-215, FTPI-MINN-20-36, UMN-TH-4005/20, HU-EP-20/37, DESY 20-222, ADP-20-33/T1143, Imperial/TP/2020/RT/04, UCI-TR-2020-19, gambit-review-2020

Journal ref: Rep. Prog. Phys. 85 052201 (2022)

Particle spectra from dark matter annihilation: physics modeling and QCD uncertainties

Authors: Simone Amoroso, Sascha Caron, Adil Jueid, Roberto Ruiz de Austri, Peter Skands

Abstract: In this talk, we discuss the physics modelling of particle spectra arising from dark matter (DM) annihilation or decay. In the context of the indirect searches of DM, the final state products will, in general, undergo a set of complicated processes such as resonance decays, QED/QCD radiation, hadronisation and hadron decays. This set of processes lead to stable particles (photons, positrons, anti-… ▽ More In this talk, we discuss the physics modelling of particle spectra arising from dark matter (DM) annihilation or decay. In the context of the indirect searches of DM, the final state products will, in general, undergo a set of complicated processes such as resonance decays, QED/QCD radiation, hadronisation and hadron decays. This set of processes lead to stable particles (photons, positrons, anti-protons, and neutrinos among others) which travel for very long distances before reaching the detectors. The modelling of their spectra contains some uncertainties which are often neglected in the relevant analyses. We discuss the sources of these uncertainties and estimate their impact on photon energy spectra for benchmark DM scenarios with $m_χ\in [10, 1000]\,$GeV. Instructions for how to retrieve complete tables from Zenodo are also provided. △ Less

Submitted 16 December, 2020; originally announced December 2020.

Comments: 14 pages, 5 figures, Contribution for the talk presented in Tools for High Energy Physics and Cosmology workshop (TOOLS2020), 2-6 November 2020, IP2I Lyon, France. Based on arXiv:1812.07424. Relevant data can be found on Zenodo (http://doi.org/10.5281/zenodo.3764809)

Prospects for discovering supersymmetric long-lived particles with MoEDAL

Authors: D. Felea, J. Mamuzic, R. Masełek, N. E. Mavromatos, V. A. Mitsou, J. L. Pinfold, R. Ruiz de Austri, K. Sakurai, A. Santra, O. Vives

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 neutra… ▽ More 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. △ Less

Submitted 24 May, 2020; v1 submitted 16 January, 2020; originally announced January 2020.

Comments: 12 pages, 6 figures; preliminary results presented in arXiv:1903.11022; matches published version in EPJC

Report number: IFIC/19-58, FTUV-19-1216, KCL-PH-TH/2019-93

Journal ref: Eur.Phys.J.C 80 (2020) 5, 431

Estimating QCD uncertainties in Monte Carlo event generators for gamma-ray dark matter searches

Authors: Simone Amoroso, Sascha Caron, Adil Jueid, Roberto Ruiz de Austri, Peter Skands

Abstract: Motivated by the recent galactic center gamma-ray excess identified in the Fermi-LAT data, we perform a detailed study of QCD fragmentation uncertainties in the modeling of the energy spectra of gamma-rays from Dark-Matter (DM) annihilation. When Dark-Matter particles annihilate to coloured final states, either directly or via decays such as $W^{(*)}\to q\bar{q}'$, photons are produced from a comp… ▽ More Motivated by the recent galactic center gamma-ray excess identified in the Fermi-LAT data, we perform a detailed study of QCD fragmentation uncertainties in the modeling of the energy spectra of gamma-rays from Dark-Matter (DM) annihilation. When Dark-Matter particles annihilate to coloured final states, either directly or via decays such as $W^{(*)}\to q\bar{q}'$, photons are produced from a complex sequence of shower, hadronisation and hadron decays. In phenomenological studies, their energy spectra are typically computed using Monte Carlo event generators. These results have however intrinsic uncertainties due to the specific model used and the choice of model parameters, which are difficult to asses and which are typically neglected. We derive a new set of hadronisation parameters (tunes) for the \textsc{Pythia~8.2} Monte Carlo generator from a fit to LEP and SLD data at the $Z$ peak. For the first time, we also derive a conservative set of uncertainties on the shower and hadronisation model parameters. Their impact on the gamma-ray energy spectra is evaluated and discussed for a range of DM masses and annihilation channels. The spectra and their uncertainties are also provided in tabulated form for future use. The fragmentation-parameter uncertainties may be useful for collider studies as well. △ Less

Submitted 1 May, 2020; v1 submitted 18 December, 2018; originally announced December 2018.

Comments: Tables for the particle spectra with their uncertainties can be found in https://zenodo.org/record/3764809#.XqvK0y-cZQJ; Discussion slightly improved, and some references are added. Comments are welcome!

Journal ref: JCAP 1905 (2019) 007

Statistical challenges in the search for dark matter

Authors: Sara Algeri, Melissa van Beekveld, Nassim Bozorgnia, Alyson Brooks, J. Alberto Casas, Jessi Cisewski-Kehe, Francis-Yan Cyr-Racine, Thomas D. P. Edwards, Fabio Iocco, Bradley J. Kavanagh, Judita Mamužić, Siddharth Mishra-Sharma, Wolfgang Rau, Roberto Ruiz de Austri, Benjamin R. Safdi, Pat Scott, Tracy R. Slatyer, Yue-Lin Sming Tsai, Aaron C. Vincent, Christoph Weniger, Jennifer Rittenhouse West, Robert L. Wolpert

Abstract: The search for the particle nature of dark matter has given rise to a number of experimental, theoretical and statistical challenges. Here, we report on a number of these statistical challenges and new techniques to address them, as discussed in the DMStat workshop held Feb 26 - Mar 3 2018 at the Banff International Research Station for Mathematical Innovation and Discovery (BIRS) in Banff, Albert… ▽ More The search for the particle nature of dark matter has given rise to a number of experimental, theoretical and statistical challenges. Here, we report on a number of these statistical challenges and new techniques to address them, as discussed in the DMStat workshop held Feb 26 - Mar 3 2018 at the Banff International Research Station for Mathematical Innovation and Discovery (BIRS) in Banff, Alberta. △ Less

Submitted 24 July, 2018; originally announced July 2018.

Comments: 32 pages, 8 figures, 331 references

Report number: IPPP/18/60

Identifying WIMP dark matter from particle and astroparticle data

Authors: Gianfranco Bertone, Nassim Bozorgnia, Jong Soo Kim, Sebastian Liem, Christopher McCabe, Sydney Otten, Roberto Ruiz de Austri

Abstract: One of the most promising strategies to identify the nature of dark matter consists in the search for new particles at accelerators and with so-called direct detection experiments. Working within the framework of simplified models, and making use of machine learning tools to speed up statistical inference, we address the question of what we can learn about dark matter from a detection at the LHC a… ▽ More One of the most promising strategies to identify the nature of dark matter consists in the search for new particles at accelerators and with so-called direct detection experiments. Working within the framework of simplified models, and making use of machine learning tools to speed up statistical inference, we address the question of what we can learn about dark matter from a detection at the LHC and a forthcoming direct detection experiment. We show that with a combination of accelerator and direct detection data, it is possible to identify newly discovered particles as dark matter, by reconstructing their relic density assuming they are weakly interacting massive particles (WIMPs) thermally produced in the early Universe, and demonstrating that it is consistent with the measured dark matter abundance. An inconsistency between these two quantities would instead point either towards additional physics in the dark sector, or towards a non-standard cosmology, with a thermal history substantially different from that of the standard cosmological model. △ Less

Submitted 3 May, 2018; v1 submitted 13 December, 2017; originally announced December 2017.

Comments: 24 pages (+21 pages of appendices and references) and 14 figures. v2: Updated to match JCAP version; includes minor clarifications in text and updated references

Report number: IPPP/17/94, KCL-PH-TH/2017-55, TTK-17-45

Journal ref: JCAP 1803 (2018) no.03, 026

Extended Higgs-portal dark matter and the Fermi-LAT Galactic Center Excess

Authors: J. A. Casas, G. A. Gomez-Vargas, J. M. Moreno, J. Quilis, R. Ruiz de Austri

Abstract: We show that the Galactic Center Excess (GCE) emission, as recently updated by the Fermi-LAT Collaboration, could be explained by the sum of Fermi-bubbles-like emission plus dark matter (DM) annihilation, in the context of a scalar-singlet Higgs portal scenario (SHP). In fact, the standard SHP, where the DM particle, $S$, only has renormalizable interactions with the Higgs, is non-operational due… ▽ More We show that the Galactic Center Excess (GCE) emission, as recently updated by the Fermi-LAT Collaboration, could be explained by the sum of Fermi-bubbles-like emission plus dark matter (DM) annihilation, in the context of a scalar-singlet Higgs portal scenario (SHP). In fact, the standard SHP, where the DM particle, $S$, only has renormalizable interactions with the Higgs, is non-operational due to strong constraints, specially from DM direct detection limits. Thus we consider a most economical extension, called ESHP (for extended SHP), which simply consists in the addition of a second (heavier) scalar singlet in the dark sector. The second scalar can be integrated-out, leaving a standard SHP plus a dimension-6 operator. Essentially, this model has only two relevant parameters (the DM mass and the coupling of the dim-6 operator). DM annihilation occurs mainly into two Higgs bosons, $SS\rightarrow hh$. We demonstrate that, despite its economy, the ESHP model provides excellent fits to the GCE (with p-value $\sim 0.6-0.7$) for very reasonable values of the parameters, in particular $m_S \simeq 130$ GeV. This is achieved without conflicting with other observables and keeping the $S-$particle relic density at the accepted value for the DM content in the universe. △ Less

Submitted 25 July, 2018; v1 submitted 29 November, 2017; originally announced November 2017.

Comments: 20 pages, 7 figures. Version published in JCAP

Report number: IFT-UAM/CSIC-17-113

Implications of the Fermi-LAT Pass 8 Galactic Center Excess on Supersymmetric Dark Matter

Authors: Abraham Achterberg, Melissa van Beekveld, Sascha Caron, Germán A. Gómez-Vargas, Luc Hendriks, Roberto Ruiz de Austri

Abstract: The Fermi Collaboration has recently updated their analysis of gamma rays from the center of the Galaxy. They reconfirm the presence of an unexplained emission feature which is most prominent in the region of $1-10$ GeV, known as the Galactic Center GeV excess (GCE). Although the GCE is now firmly detected, an interpretation of this emission as a signal of self-annihilating dark matter (DM) partic… ▽ More The Fermi Collaboration has recently updated their analysis of gamma rays from the center of the Galaxy. They reconfirm the presence of an unexplained emission feature which is most prominent in the region of $1-10$ GeV, known as the Galactic Center GeV excess (GCE). Although the GCE is now firmly detected, an interpretation of this emission as a signal of self-annihilating dark matter (DM) particles is not unambiguously possible due to systematic effects in the gamma-ray modeling estimated in the Galactic Plane. In this paper we build a covariance matrix, collecting different systematic uncertainties investigated in the Fermi Collaboration's paper that affect the GCE spectrum. We show that models where part of the GCE is due to annihilating DM can still be consistent with the new data. We also re-evaluate the parameter space regions of the minimal supersymmetric Standard Model (MSSM) that can contribute dominantly to the GCE via neutralino DM annihilation. All recent constraints from DM direct detection experiments such as PICO, LUX, PandaX and Xenon1T, limits on the annihilation cross section from dwarf spheroidal galaxies and Large Hadron Collider limits are considered in this analysis. Due to a slight shift in the energy spectrum of the GC excess with respect to the previous Fermi analysis, and the recent limits from direct detection experiments, we find a slightly shifted parameter region of the MSSM compared to our previous analysis that is consistent with the GCE. Neutralinos with a mass between $85-220$ GeV can describe the excess via annihilation into a pair of $W$-bosons or top quarks. Remarkably, there are low fine-tuning models among the regions that we have found. The complete set of solutions will be probed by upcoming direct detection experiments and with dedicated searches in the upcoming data of the Large Hadron Collider. △ Less

Submitted 10 January, 2018; v1 submitted 29 September, 2017; originally announced September 2017.

Comments: Replaced with the published version in JCAP. Auxiliary material on the construction of the covariance matrix and the GCE spectrum is included

Journal ref: JCAP12(2017)040

Analyzing γ-rays of the Galactic Center with Deep Learning

Authors: Sascha Caron, Germán A. Gómez-Vargas, Luc Hendriks, Roberto Ruiz de Austri

Abstract: We present a new method to interpret the $γ$-ray data of our inner Galaxy as measured by the Fermi Large Area Telescope (Fermi LAT). We train and test convolutional neural networks with simulated Fermi-LAT images based on models tuned to real data. We use this method to investigate the origin of an excess emission of GeV $γ$-rays seen in previous studies. Interpretations of this excess include… ▽ More We present a new method to interpret the $γ$-ray data of our inner Galaxy as measured by the Fermi Large Area Telescope (Fermi LAT). We train and test convolutional neural networks with simulated Fermi-LAT images based on models tuned to real data. We use this method to investigate the origin of an excess emission of GeV $γ$-rays seen in previous studies. Interpretations of this excess include $γ$ rays created by the annihilation of dark matter particles and $γ$ rays originating from a collection of unresolved point sources, such as millisecond pulsars. Our new method allows precise measurements of the contribution and properties of an unresolved population of $γ$-ray point sources in the interstellar diffuse emission model. △ Less

Submitted 26 May, 2018; v1 submitted 22 August, 2017; originally announced August 2017.

Comments: 24 pages, 11 figures

Journal ref: Journal of Cosmology and Astroparticle Physics, Volume 2018, May 2018

Global fits of GUT-scale SUSY models with GAMBIT

Authors: The GAMBIT Collaboration, Peter Athron, Csaba Balázs, Torsten Bringmann, Andy Buckley, Marcin Chrząszcz, Jan Conrad, Jonathan M. Cornell, Lars A. Dal, Joakim Edsjö, Ben Farmer, Paul Jackson, Abram Krislock, Anders Kvellestad, Farvah Mahmoudi, Gregory D. Martinez, Antje Putze, Are Raklev, Christopher Rogan, Roberto Ruiz de Austri, Aldo Saavedra, Christopher Savage, Pat Scott, Nicola Serra, Christoph Weniger , et al. (1 additional authors not shown)

Abstract: We present the most comprehensive global fits to date of three supersymmetric models motivated by grand unification: the Constrained Minimal Supersymmetric Standard Model (CMSSM), and its Non-Universal Higgs Mass generalisations NUHM1 and NUHM2. We include likelihoods from a number of direct and indirect dark matter searches, a large collection of electroweak precision and flavour observables, dir… ▽ More We present the most comprehensive global fits to date of three supersymmetric models motivated by grand unification: the Constrained Minimal Supersymmetric Standard Model (CMSSM), and its Non-Universal Higgs Mass generalisations NUHM1 and NUHM2. We include likelihoods from a number of direct and indirect dark matter searches, a large collection of electroweak precision and flavour observables, direct searches for supersymmetry at LEP and Runs I and II of the LHC, and constraints from Higgs observables. Our analysis improves on existing results not only in terms of the number of included observables, but also in the level of detail with which we treat them, our sampling techniques for scanning the parameter space, and our treatment of nuisance parameters. We show that stau co-annihilation is now ruled out in the CMSSM at more than 95\% confidence. Stop co-annihilation turns out to be one of the most promising mechanisms for achieving an appropriate relic density of dark matter in all three models, whilst avoiding all other constraints. We find high-likelihood regions of parameter space featuring light stops and charginos, making them potentially detectable in the near future at the LHC. We also show that tonne-scale direct detection will play a largely complementary role, probing large parts of the remaining viable parameter space, including essentially all models with multi-TeV neutralinos. △ Less

Submitted 23 October, 2018; v1 submitted 22 May, 2017; originally announced May 2017.

Comments: 50 pages, 21 figures, 7 tables, v2 accepted for publication in EPJC, v3 update Zenodo link

Report number: CoEPP-MN-17-9, gambit-physics-2017

Search for sharp and smooth spectral signatures of $μν$SSM gravitino dark matter with Fermi-LAT

Authors: German A. Gomez-Vargas, Daniel E. Lopez-Fogliani, Carlos Munoz, Andres D. Perez, Roberto Ruiz de Austri

Abstract: The $μν$SSM solves the $μ$ problem of supersymmetric models and reproduces neutrino data, simply using couplings with right-handed neutrinos $ν$'s. Given that these couplings break explicitly $R$ parity, the gravitino is a natural candidate for decaying dark matter in the $μν$SSM. In this work we carry out a complete analysis of the detection of $μν$SSM gravitino dark matter through $γ$-ray observ… ▽ More The $μν$SSM solves the $μ$ problem of supersymmetric models and reproduces neutrino data, simply using couplings with right-handed neutrinos $ν$'s. Given that these couplings break explicitly $R$ parity, the gravitino is a natural candidate for decaying dark matter in the $μν$SSM. In this work we carry out a complete analysis of the detection of $μν$SSM gravitino dark matter through $γ$-ray observations. In addition to the two-body decay producing a sharp line, we include in the analysis the three-body decays producing a smooth spectral signature. We perform first a deep exploration of the low-energy parameter space of the $μν$SSM taking into account that neutrino data must be reproduced. Then, we compare the $γ$-ray fluxes predicted by the model with Fermi-LAT observations. In particular, with the 95$\%$ CL upper limits on the total diffuse extragalactic $γ$-ray background using 50 months of data, together with the upper limits on line emission from an updated analysis using 69.9 months of data. For standard values of bino and wino masses, gravitinos with masses larger than about 4 GeV, or lifetimes smaller than $10^{28}$ s, produce too large fluxes and are excluded as dark matter candidates. However, when limiting scenarios with large and close values of the gaugino masses are considered, the constraints turn out to be less stringent, excluding masses larger than 17 GeV and lifetimes smaller than $4\times 10^{25}$ s. △ Less

Submitted 17 March, 2017; v1 submitted 30 August, 2016; originally announced August 2016.

Comments: Minor changes, references added, version published in JCAP. 23 pages, 7 figures, 3 tables

Effective Field Theory of Dark Matter: a Global Analysis

Authors: Sebastian Liem, Gianfranco Bertone, Francesca Calore, Roberto Ruiz de Austri, Tim M. P. Tait, Roberto Trotta, Christoph Weniger

Abstract: We present global fits of an effective field theory description of real, and complex scalar dark matter candidates. We simultaneously take into account all possible dimension 6 operators consisting of dark matter bilinears and gauge invariant combinations of quark and gluon fields. We derive constraints on the free model parameters for both the real (five parameters) and complex (seven) scalar dar… ▽ More We present global fits of an effective field theory description of real, and complex scalar dark matter candidates. We simultaneously take into account all possible dimension 6 operators consisting of dark matter bilinears and gauge invariant combinations of quark and gluon fields. We derive constraints on the free model parameters for both the real (five parameters) and complex (seven) scalar dark matter models obtained by combining Planck data on the cosmic microwave background, direct detection limits from LUX, and indirect detection limits from the Fermi Large Area Telescope. We find that for real scalars indirect dark matter searches disfavour a dark matter particle mass below 100 GeV. For the complex scalar dark matter particle current data have a limited impact due to the presence of operators that lead to p-wave annihilation, and also do not contribute to the spin-independent scattering cross- section. Although current data are not informative enough to strongly constrain the theory parameter space, we demonstrate the power of our formalism to reconstruct the theoretical parameters compatible with an actual dark matter detection, by assuming that the excess of gamma rays observed by the Fermi Large Area Telescope towards the Galactic centre is entirely due to dark matter annihilations. Please note that the excess can very well be due to astrophysical sources such as millisecond pulsars. We find that scalar dark matter interacting via effective field theory operators can in principle explain the Galactic centre excess, but that such interpretation is in strong tension with the non-detection of gamma rays from dwarf galaxies in the real scalar case. In the complex scalar case there is enough freedom to relieve the tension. △ Less

Submitted 18 March, 2016; originally announced March 2016.

Comments: 20 pages, 5 figures

Bayesian analysis of cosmic-ray propagation: evidence against homogeneous diffusion

Authors: G. Jóhannesson, R. Ruiz de Austri, A. C. Vincent, I. V. Moskalenko, E. Orlando, T. A. Porter, A. W. Strong, R. Trotta, F. Feroz, P. Graff, M. P. Hobson

Abstract: We present the results of the most complete ever 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 such study to separate out low-mass isotopes ($p$, $\bar p$ and He) from the usual lig… ▽ More We present the results of the most complete ever 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 such study to separate out low-mass isotopes ($p$, $\bar p$ and He) from the usual light elements (Be, B, C, N, O). We find that the propagation parameters that best fit $p$, $\bar p$, He data are significantly different from those that fit light elements, including the B/C and $^{10}$Be/$^9$Be secondary-to-primary ratios normally used to calibrate propagation parameters. This suggests 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. △ Less

Submitted 12 April, 2016; v1 submitted 6 February, 2016; originally announced February 2016.

Comments: 20 pages, 10 figures, 5 tables. v2 accepted for publication in ApJ

Report number: IPPP/16/10

The case for 100 GeV bino dark matter: A dedicated LHC tri-lepton search

Authors: Melissa van Beekveld, Wim Beenakker, Sascha Caron, Roberto Ruiz de Austri

Abstract: Global fit studies performed in the pMSSM and the photon excess signal originating from the Galactic Center seem to suggest compressed electroweak supersymmetric spectra with a $\sim$100 GeV bino-like dark matter particle. We find that these scenarios are not probed by traditional electroweak supersymmetry searches at the LHC. We propose to extend the ATLAS and CMS electroweak supersymmetry search… ▽ More Global fit studies performed in the pMSSM and the photon excess signal originating from the Galactic Center seem to suggest compressed electroweak supersymmetric spectra with a $\sim$100 GeV bino-like dark matter particle. We find that these scenarios are not probed by traditional electroweak supersymmetry searches at the LHC. We propose to extend the ATLAS and CMS electroweak supersymmetry searches with an improved strategy for bino-like dark matter, focusing on chargino plus next-to-lightest neutralino production, with a subsequent decay into a tri-lepton final state. We explore the sensitivity for pMSSM scenarios with $Δm = m_{\rm NLSP} - m_{\rm LSP} \sim (5 - 50)$ GeV in the $\sqrt{s} = 14$ TeV run of the LHC. Counterintuitively, we find that the requirement of low missing transverse energy increases the sensitivity compared to the current ATLAS and CMS searches. With 300 fb$^{-1}$ of data we expect the LHC experiments to be able to discover these supersymmetric spectra with mass gaps down to $Δm \sim 9$ GeV for DM masses between 40 and 140 GeV. We stress the importance of a dedicated search strategy that targets precisely these favored pMSSM spectra. △ Less

Submitted 2 May, 2016; v1 submitted 1 February, 2016; originally announced February 2016.

Comments: Published in JHEP

A description of the Galactic Center excess in the Minimal Supersymmetric Standard Model

Authors: Abraham Achterberg, Simone Amoroso, Sascha Caron, Luc Hendriks, Roberto Ruiz de Austri, Christoph Weniger

Abstract: Observations with the Fermi Large Area Telescope (LAT) indicate an excess in gamma rays originating from the center of our Galaxy. A possible explanation for this excess is the annihilation of Dark Matter particles. We have investigated the annihilation of neutralinos as Dark Matter candidates within the phenomenological Minimal Supersymmetric Standard Model (pMSSM). An iterative particle filter a… ▽ More Observations with the Fermi Large Area Telescope (LAT) indicate an excess in gamma rays originating from the center of our Galaxy. A possible explanation for this excess is the annihilation of Dark Matter particles. We have investigated the annihilation of neutralinos as Dark Matter candidates within the phenomenological Minimal Supersymmetric Standard Model (pMSSM). An iterative particle filter approach was used to search for solutions within the pMSSM. We found solutions that are consistent with astroparticle physics and collider experiments, and provide a fit to the energy spectrum of the excess. The neutralino is a Bino/Higgsino or Bino/Wino/Higgsino mixture with a mass in the range $84-92$~GeV or $87-97$~GeV annihilating into W bosons. A third solutions is found for a neutralino of mass $174-187$~GeV annihilating into top quarks. The best solutions yield a Dark Matter relic density $0.06 < Ωh^2 <0.13$. These pMSSM solutions make clear forecasts for LHC, direct and indirect DM detection experiments. If the MSSM explanation of the excess seen by Fermi-LAT is correct, a DM signal might be discovered soon. △ Less

Submitted 31 March, 2015; v1 submitted 19 February, 2015; originally announced February 2015.

Comments: Large extension of previous paper: 2 more solutions found in the MSSM (Bino-Higgsino, Bino-Wino-Higgsino into WW and Bino into ttbar), added description on extra fit uncertainties, added description on flavor observables, added discussion on dwarf limits

Profile likelihood maps of a 15-dimensional MSSM

Authors: C. Strege, G. Bertone, G. J. Besjes, S. Caron, R. Ruiz de Austri, A. Strubig, R. Trotta

Abstract: We present statistically convergent profile likelihood maps obtained via global fits of a phenomenological Minimal Supersymmetric Standard Model with 15 free parameters (the MSSM-15), based on over 250M points. We derive constraints on the model parameters from direct detection limits on dark matter, the Planck relic density measurement and data from accelerator searches. We provide a detailed ana… ▽ More We present statistically convergent profile likelihood maps obtained via global fits of a phenomenological Minimal Supersymmetric Standard Model with 15 free parameters (the MSSM-15), based on over 250M points. We derive constraints on the model parameters from direct detection limits on dark matter, the Planck relic density measurement and data from accelerator searches. We provide a detailed analysis of the rich phenomenology of this model, and determine the SUSY mass spectrum and dark matter properties that are preferred by current experimental constraints. We evaluate the impact of the measurement of the anomalous magnetic moment of the muon ($g-2$) on our results, and provide an analysis of scenarios in which the lightest neutralino is a subdominant component of the dark matter. The MSSM-15 parameters are relatively weakly constrained by current data sets, with the exception of the parameters related to dark matter phenomenology ($M_1$, $M_2$, $μ$), which are restricted to the sub-TeV regime, mainly due to the relic density constraint. The mass of the lightest neutralino is found to be < 1.5 TeV at 99% C.L., but can extend up to 3 TeV when excluding the $g - 2$ constraint from the analysis. Low-mass bino-like neutralinos are strongly favoured, with spin-independent scattering cross-sections extending to very small values, $\sim 10^{-20}$ pb. ATLAS SUSY null searches strongly impact on this mass range, and thus rule out a region of parameter space that is outside the reach of any current or future direct detection experiment. The best-fit point obtained after inclusion of all data corresponds to a squark mass of 2.3 TeV, a gluino mass of 2.1 TeV and a 130 GeV neutralino with a spin-independent cross-section of $2.4 \times 10^{-10}$ pb, which is within the reach of future multi-ton scale direct detection experiments and of the upcoming LHC run at increased centre-of-mass energy. △ Less

Submitted 3 May, 2014; originally announced May 2014.

Comments: 52 pages, 21 figures, submitted to JHEP

Impact of nucleon matrix element uncertainties on the interpretation of direct and indirect dark matter search results

Authors: R. Ruiz de Austri, C. Pérez de los Heros

Abstract: We study in detail the impact of the current uncertainty in nucleon matrix elements on the sensitivity of direct and indirect experimental techniques for dark matter detection. We perform two scans in the framework of the cMSSM: one using recent values of the pion-sigma term obtained from Lattice QCD, and the other using values derived from experimental measurements. The two choices correspond to… ▽ More We study in detail the impact of the current uncertainty in nucleon matrix elements on the sensitivity of direct and indirect experimental techniques for dark matter detection. We perform two scans in the framework of the cMSSM: one using recent values of the pion-sigma term obtained from Lattice QCD, and the other using values derived from experimental measurements. The two choices correspond to extreme values quoted in the literature and reflect the current tension between different ways of obtaining information about the structure of the nucleon. All other inputs in the scans, astrophysical and from particle physics, are kept unchanged. We use two experiments, XENON100 and IceCube, as benchmark cases to illustrate our case. We find that the interpretation of dark matter search results from direct detection experiments is more sensitive to the choice of the central values of the hadronic inputs than the results of indirect search experiments. The allowed regions of cMSSM parameter space after including XENON100 constrains strongly differ depending on the assumptions on the hadronic matrix elements used. On the other hand, the constraining potential of IceCube is almost independent of the choice of these values. △ Less

Submitted 22 October, 2013; v1 submitted 25 July, 2013; originally announced July 2013.

Comments: Version accepted in JCAP

Journal ref: JCAP11(2013)049

Global Fits of the cMSSM and NUHM including the LHC Higgs discovery and new XENON100 constraints

Authors: C. Strege, G. Bertone, F. Feroz, M. Fornasa, R. Ruiz de Austri, R. Trotta

Abstract: We present global fits of the constrained Minimal Supersymmetric Standard Model (cMSSM) and the Non-Universal Higgs Model (NUHM), including the most recent CMS constraint on the Higgs boson mass, 5.8/fb integrated luminosity null Supersymmetry searches by ATLAS, the new LHCb measurement of the Bs to mu+mu- branching ratio and the 7-year WMAP dark matter relic abundance determination. We include th… ▽ More We present global fits of the constrained Minimal Supersymmetric Standard Model (cMSSM) and the Non-Universal Higgs Model (NUHM), including the most recent CMS constraint on the Higgs boson mass, 5.8/fb integrated luminosity null Supersymmetry searches by ATLAS, the new LHCb measurement of the Bs to mu+mu- branching ratio and the 7-year WMAP dark matter relic abundance determination. We include the latest dark matter constraints from the XENON100 experiment, marginalising over astrophysical and particle physics uncertainties. We present Bayesian posterior and profile likelihood maps of the highest resolution available today, obtained from up to 350M points. We find that the new constraint on the Higgs boson mass has a dramatic impact, ruling out large regions of previously favoured cMSSM and NUHM parameter space. In the cMSSM, light sparticles and predominantly gaugino-like dark matter with a mass of a few hundred GeV are favoured. The NUHM exhibits a strong preference for heavier sparticle masses and a Higgsino-like neutralino with a mass of 1 TeV. The future ton-scale XENON1T direct detection experiment will probe large portions of the currently favoured cMSSM and NUHM parameter space. The LHC operating at 14 TeV collision energy will explore the favoured regions in the cMSSM, while most of the regions favoured in the NUHM will remain inaccessible. Our best-fit points achieve a satisfactory quality-of-fit, with p-values ranging from 0.21 to 0.35, so that none of the two models studied can be presently excluded at any meaningful significance level. △ Less

Submitted 11 December, 2012; originally announced December 2012.

Comments: 37 pages, 11 figures, submitted to JCAP

Journal ref: JCAP04(2013)013

Updated global fits of the cMSSM including the latest LHC SUSY and Higgs searches and XENON100 data

Authors: C. Strege, G. Bertone, D. G. Cerdeno, M. Fornasa, R. Ruiz de Austri, R. Trotta

Abstract: We present new global fits of the cMSSM, including LHC 1/fb integrated luminosity SUSY exclusion limits, recent LHC 5/fb constraints on the mass of the Higgs boson and XENON100 direct detection data. Our analysis fully takes into account astrophysical and hadronic uncertainties that enter the analysis when translating direct detection limits into constraints on the cMSSM parameter space. We provid… ▽ More We present new global fits of the cMSSM, including LHC 1/fb integrated luminosity SUSY exclusion limits, recent LHC 5/fb constraints on the mass of the Higgs boson and XENON100 direct detection data. Our analysis fully takes into account astrophysical and hadronic uncertainties that enter the analysis when translating direct detection limits into constraints on the cMSSM parameter space. We provide results for both a Bayesian and a Frequentist statistical analysis. We find that LHC 2011 constraints in combination with XENON100 data can rule out a significant portion of the cMSSM parameter space. Our results further emphasise the complementarity of collider experiments and direct detection searches in constraining extensions of Standard Model physics. The LHC 2011 exclusion limit strongly impacts on low-mass regions of cMSSM parameter space, such as the stau co-annihilation region, while direct detection data can rule out regions of high SUSY masses, such as the Focus-Point region, which is unreachable for the LHC in the near future. We show that, in addition to XENON100 data, the experimental constraint on the anomalous magnetic moment of the muon plays a dominant role in disfavouring large scalar and gaugino masses. We find that, should the LHC 2011 excess hinting towards a Higgs boson at 126 GeV be confirmed, currently favoured regions of the cMSSM parameter space will be robustly ruled out from both a Bayesian and a profile likelihood statistical perspective. △ Less

Submitted 22 February, 2012; v1 submitted 18 December, 2011; originally announced December 2011.

Comments: 18 pages, 5 figures, 2 tables; matches JCAP accepted version, profile likelihood contours slightly modified from v1, main conclusions unchanged

Report number: IFT-UAM/CSIC-11-99; FTUAM-11-66

Complementarity of Indirect and Accelerator Dark Matter Searches

Authors: G. Bertone, D. G. Cerdeno, M. Fornasa, L. Pieri, R. Ruiz de Austri, R. Trotta

Abstract: Even if Supersymmetric particles are found at the Large Hadron Collider (LHC), it will be difficult to prove that they constitute the bulk of the Dark Matter (DM) in the Universe using LHC data alone. We study the complementarity of LHC and DM indirect searches, working out explicitly the reconstruction of the DM properties for a specific benchmark model in the coannihilation region of a 24-parame… ▽ More Even if Supersymmetric particles are found at the Large Hadron Collider (LHC), it will be difficult to prove that they constitute the bulk of the Dark Matter (DM) in the Universe using LHC data alone. We study the complementarity of LHC and DM indirect searches, working out explicitly the reconstruction of the DM properties for a specific benchmark model in the coannihilation region of a 24-parameters supersymmetric model. Combining mock high-luminosity LHC data with present-day null searches for gamma-rays from dwarf galaxies with the Fermi LAT, we show that current Fermi LAT limits already have the capability of ruling out a spurious Wino-like solution that would survive using LHC data only, thus leading to the correct identification of the cosmological solution. We also demonstrate that upcoming Planck constraints on the reionization history will have a similar constraining power, and discuss the impact of a possible detection of gamma-rays from DM annihilation in Draco with a CTA-like experiment. Our results indicate that indirect searches can be strongly complementary to the LHC in identifying the DM particles, even when astrophysical uncertainties are taken into account. △ Less

Submitted 10 November, 2011; originally announced November 2011.

Comments: 9 pages, 3 figures

Report number: FTUAM-11-60, IFT-UAM/CSIC-11-79

Dark Matter Searches: The Nightmare Scenario

Authors: Gianfranco Bertone, Daniel Cumberbatch, Roberto Ruiz de Austri, Roberto Trotta

Abstract: The unfortunate case where the Large Hadron Collider (LHC) fails to discover physics Beyond the Standard Model (BSM) is sometimes referred to as the "Nightmare scenario" of particle physics. We study the consequences of this hypothetical scenario for Dark Matter (DM), in the framework of the constrained Minimal Supersymmetric Standard Model (cMSSM). We evaluate the surviving regions of the cMSSM p… ▽ More The unfortunate case where the Large Hadron Collider (LHC) fails to discover physics Beyond the Standard Model (BSM) is sometimes referred to as the "Nightmare scenario" of particle physics. We study the consequences of this hypothetical scenario for Dark Matter (DM), in the framework of the constrained Minimal Supersymmetric Standard Model (cMSSM). We evaluate the surviving regions of the cMSSM parameter space after null searches at the LHC, using several different LHC configurations, and study the consequences for DM searches with ton-scale direct detectors and the IceCube neutrino telescope. We demonstrate that ton-scale direct detection experiments will be able to conclusively probe the cMSSM parameter space that would survive null searches at the LHC with 100fb$^{-1}$ of integrated luminosity at 14TeV. We also demonstrate that IceCube (80 strings plus DeepCore) will be able to probe as much as 17% of the currently favoured parameter space after 5 years of observation. △ Less

Submitted 4 January, 2012; v1 submitted 28 July, 2011; originally announced July 2011.

Comments: V2: 24 pages, 6 figures, 4 tables. Replaced to match version published in JCAP. Minor revisions made to address referee's comments

Journal ref: JCAP 01 (2012) 004

Global fits of the cMSSM including the first LHC and XENON100 data

Authors: G. Bertone, D. G. Cerdeno, M. Fornasa, R. Ruiz de Austri, C. Strege, R. Trotta

Abstract: We present updated global fits of the constrained Minimal Supersymmetric Standard Model (cMSSM), including the most recent constraints from the ATLAS and CMS detectors at the LHC, as well as the most recent results of the XENON100 experiment. Our robust analysis takes into account both astrophysical and hadronic uncertainties that enter in the calculation of the rate of WIMP-induced recoils in dir… ▽ More We present updated global fits of the constrained Minimal Supersymmetric Standard Model (cMSSM), including the most recent constraints from the ATLAS and CMS detectors at the LHC, as well as the most recent results of the XENON100 experiment. Our robust analysis takes into account both astrophysical and hadronic uncertainties that enter in the calculation of the rate of WIMP-induced recoils in direct detection experiment. We study the consequences for neutralino Dark Matter, and show that current direct detection data already allow to robustly rule out the so-called Focus Point region, therefore demonstrating the importance of particle astrophysics experiments in constraining extensions of the Standard Model of Particle Physics. We also observe an increased compatibility between results obtained from a Bayesian and a Frequentist statistical perspective. We find that upcoming ton-scale direct detection experiments will probe essentially the entire currently favoured region (at the 99% level), almost independently of the statistical approach used. Prospects for indirect detection of the cMSSM are further reduced. △ Less

Submitted 19 January, 2012; v1 submitted 8 July, 2011; originally announced July 2011.

Comments: 14 pages, 4 figures, minor changes following referee report, matches published version

Report number: IFT-UAM/CSIC-11-52; FTUAM-11-50

Journal ref: JCAP01(2012)015

Is light neutralino as dark matter still viable?

Authors: Daniel T. Cumberbatch, Daniel E. López-Fogliani, Leszek Roszkowski, Roberto Ruiz de Austri, Yue-Lin S. Tsai

Abstract: Motivated by the recent re-confirmation by CoGENT of the low-energy excess of events observed last year, and the recent improved limits from the XENON-100 experiment that are in contention with the CoGENT data, we re-examine the low mass neutralino region of the Minimal Supersymmetric Standard Model and of the Next-to-Minimal Supersymmetric Standard Model, both without assuming gaugino mass unific… ▽ More Motivated by the recent re-confirmation by CoGENT of the low-energy excess of events observed last year, and the recent improved limits from the XENON-100 experiment that are in contention with the CoGENT data, we re-examine the low mass neutralino region of the Minimal Supersymmetric Standard Model and of the Next-to-Minimal Supersymmetric Standard Model, both without assuming gaugino mass unification. We make several focused scans for each model, determining conservative constraints on input parameters. We then determine how these constraints are made increasingly stringent as we re-invoke our experimental constraints involving the dark matter relic abundance, collider constraints from LEP and the Tevatron, and then from flavour physics, as a series of successive 2 sigma hard cuts. We find that for both models, when all relevant constraints are applied in this fashion, we do not generate neutralino LSPs that possess a spin-independent scattering cross section in excess of 10^-5 pb and a mass ~7 GeV < m_chi < ~9 GeV that is necessary in order to explain the CoGENT observations. △ Less

Submitted 8 July, 2011; originally announced July 2011.

Comments: 16 pages, 7 figures, 3 tables

Complementarity of Dark Matter Direct Detection Targets

Authors: Miguel Pato, Laura Baudis, Gianfranco Bertone, Roberto Ruiz de Austri, Louis E. Strigari, Roberto Trotta

Abstract: We investigate the reconstruction capabilities of Dark Matter mass and spin-independent cross-section from future ton-scale direct detection experiments using germanium, xenon or argon as targets. Adopting realistic values for the exposure, energy threshold and resolution of Dark Matter experiments which will come online within 5 to 10 years, the degree of complementarity between different targets… ▽ More We investigate the reconstruction capabilities of Dark Matter mass and spin-independent cross-section from future ton-scale direct detection experiments using germanium, xenon or argon as targets. Adopting realistic values for the exposure, energy threshold and resolution of Dark Matter experiments which will come online within 5 to 10 years, the degree of complementarity between different targets is quantified. We investigate how the uncertainty in the astrophysical parameters controlling the local Dark Matter density and velocity distribution affects the reconstruction. For a 50 GeV WIMP, astrophysical uncertainties degrade the accuracy in the mass reconstruction by up to a factor of $\sim 4$ for xenon and germanium, compared to the case when astrophysical quantities are fixed. However, combination of argon, germanium and xenon data increases the constraining power by a factor of $\sim 2$ compared to germanium or xenon alone. We show that future direct detection experiments can achieve self-calibration of some astrophysical parameters, and they will be able to constrain the WIMP mass with only very weak external astrophysical constraints. △ Less

Submitted 20 June, 2011; v1 submitted 15 December, 2010; originally announced December 2010.

Comments: 11 pages, 5 figures, 4 tables, matches published version

Journal ref: Phys.Rev.D83:083505,2011

Constraints on cosmic-ray propagation models from a global Bayesian analysis

Authors: R. Trotta, G. Johannesson, I. V. Moskalenko, T. A. Porter, R. Ruiz de Austri, A. W. Strong

Abstract: Research in many areas of modern physics such as, e.g., indirect searches for dark matter and particle acceleration in SNR 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… ▽ More Research in many areas of modern physics such as, e.g., indirect searches for dark matter and particle acceleration in SNR 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, that uses astrophysical information, nuclear and particle data as input 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. △ Less

Submitted 6 January, 2011; v1 submitted 29 October, 2010; originally announced November 2010.

Comments: 19 figures, 3 tables, emulateapj.sty. A typo is fixed. To be published in the Astrophysical Journal v.728 (February 10, 2011 issue). Supplementary material can be found at http://www.g-vo.org/pub/GALPROP/GalpropBayesPaper/

Journal ref: ApJ, 729, 106 (2011)

Global fits of the Minimal Universal Extra Dimensions scenario

Authors: Gianfranco Bertone, Kyoungchul Kong, Roberto Ruiz de Austri, Roberto Trotta

Abstract: In theories with Universal Extra-Dimensions (UED), the gamma_1 particle, first excited state of the hypercharge gauge boson, provides an excellent Dark Matter (DM) candidate. Here we use a modified version of the SuperBayeS code to perform a Bayesian analysis of the minimal UED scenario, in order to assess its detectability at accelerators and with DM experiments. We derive in particular the most… ▽ More In theories with Universal Extra-Dimensions (UED), the gamma_1 particle, first excited state of the hypercharge gauge boson, provides an excellent Dark Matter (DM) candidate. Here we use a modified version of the SuperBayeS code to perform a Bayesian analysis of the minimal UED scenario, in order to assess its detectability at accelerators and with DM experiments. We derive in particular the most probable range of mass and scattering cross sections off nucleons, keeping into account cosmological and electroweak precision constraints. The consequences for the detectability of the gamma_1 with direct and indirect experiments are dramatic. The spin-independent cross section probability distribution peaks at ~ 10^{-11} pb, i.e. below the sensitivity of ton-scale experiments. The spin-dependent cross-section drives the predicted neutrino flux from the center of the Sun below the reach of present and upcoming experiments. The only strategy that remains open appears to be direct detection with ton-scale experiments sensitive to spin-dependent cross-sections. On the other hand, the LHC with 1 1/fb of data should be able to probe the current best-fit UED parameters. △ Less

Submitted 16 December, 2010; v1 submitted 11 October, 2010; originally announced October 2010.

Comments: 15 pages, 11 figures; minor quantitative differences wrt v1 due to the correction of a numerical problem, main results unchanged

Journal ref: Phys.Rev.D83:036008,2011

Identification of Dark Matter particles with LHC and direct detection data

Authors: Gianfranco Bertone, David G. Cerdeno, Mattia Fornasa, Roberto Ruiz de Austri, Roberto Trotta

Abstract: Dark matter (DM) is currently searched for with a variety of detection strategies. Accelerator searches are particularly promising, but even if Weakly Interacting Massive Particles (WIMPs) are found at the Large Hadron Collider (LHC), it will be difficult to prove that they constitute the bulk of the DM in the Universe. We show that a significantly better reconstruction of the DM properties can be… ▽ More Dark matter (DM) is currently searched for with a variety of detection strategies. Accelerator searches are particularly promising, but even if Weakly Interacting Massive Particles (WIMPs) are found at the Large Hadron Collider (LHC), it will be difficult to prove that they constitute the bulk of the DM in the Universe. We show that a significantly better reconstruction of the DM properties can be obtained with a combined analysis of LHC and direct detection (DD) data, by making a simple Ansatz on the WIMP local density, i.e. by assuming that the local density scales with the cosmological relic abundance. We demonstrate this method in an explicit example in the context of a 24-parameter supersymmetric model, with a neutralino LSP in the stau co-annihilation region. Our results show that future ton-scale DD experiments will allow to break degeneracies in the SUSY parameter space and achieve a significantly better reconstruction of the neutralino composition and its relic density than with LHC data alone. △ Less

Submitted 12 October, 2010; v1 submitted 24 May, 2010; originally announced May 2010.

Comments: 7 pages, 2 figures, comments: replaced with the journal version

Journal ref: Phys.Rev.D82:055008,2010

The degenerate gravitino scenario

Authors: Lotfi Boubekeur, Ki Young Choi, Roberto Ruiz de Austri, Oscar Vives

Abstract: In this work, we explore the "degenerate gravitino" scenario where the mass difference between the gravitino and the lightest MSSM particle is much smaller than the gravitino mass itself. In this case, the energy released in the decay of the next to lightest sypersymmetric particle (NLSP) is reduced. Consequently the cosmological and astrophysical constraints on the gravitino abundance, and henc… ▽ More In this work, we explore the "degenerate gravitino" scenario where the mass difference between the gravitino and the lightest MSSM particle is much smaller than the gravitino mass itself. In this case, the energy released in the decay of the next to lightest sypersymmetric particle (NLSP) is reduced. Consequently the cosmological and astrophysical constraints on the gravitino abundance, and hence on the reheating temperature, become softer than in the usual case. On the other hand, such small mass splittings generically imply a much longer lifetime for the NLSP. We find that, in the constrained MSSM (CMSSM), for neutralino LSP or NLSP, reheating temperatures compatible with thermal leptogenesis are reached for small splittings of order 10^{-2} GeV. While for stau NLSP, temperatures of 4x10^9 GeV can be obtained even for splittings of order of tens of GeVs. This "degenerate gravitino" scenario offers a possible way out to the gravitino problem for thermal leptogenesis in supersymmetric theories. △ Less

Submitted 7 April, 2010; v1 submitted 2 February, 2010; originally announced February 2010.

Comments: 27 pages, 10 figures and 1 table. Minor typos and references fixed. Matches published version in JCAP.

Report number: FTUV-10-0126, IFIC/10-03, PNUTP-10-A03.

Journal ref: JCAP 1004:005,2010

Gamma-ray detection from gravitino dark matter decay in the $μν$SSM

Authors: Ki-Young Choi, Daniel E. Lopez-Fogliani, Carlos Munoz, Roberto Ruiz de Austri

Abstract: The $μν$SSM provides a solution to the $μ$-problem of the MSSM and explains the origin of neutrino masses by simply using right-handed neutrino superfields. Given that R-parity is broken in this model, the gravitino is a natural candidate for dark matter since its lifetime becomes much longer than the age of the Universe. We consider the implications of gravitino dark matter in the $μν$SSM, anal… ▽ More The $μν$SSM provides a solution to the $μ$-problem of the MSSM and explains the origin of neutrino masses by simply using right-handed neutrino superfields. Given that R-parity is broken in this model, the gravitino is a natural candidate for dark matter since its lifetime becomes much longer than the age of the Universe. We consider the implications of gravitino dark matter in the $μν$SSM, analyzing in particular the prospects for detecting gamma rays from decaying gravitinos. If the gravitino explains the whole dark matter component, a gravitino mass larger than 20 GeV is disfavored by the isotropic diffuse photon background measurements. On the other hand, a gravitino with a mass range between 0.1-20 GeV gives rise to a signal that might be observed by the FERMI satellite. In this way important regions of the parameter space of the $μν$SSM can be checked. △ Less

Submitted 25 February, 2010; v1 submitted 19 June, 2009; originally announced June 2009.

Comments: Final version to appear in JCAP, 13 pages, 3 figures

Report number: FTUAM 09/15, IFT-UAM/CSIC-09-29

Prospects for dark matter detection with IceCube in the context of the CMSSM

Authors: R. Trotta, R. Ruiz de Austri, C. Pérez de los Heros

Abstract: We study in detail the ability of the nominal configuration of the IceCube neutrino telescope (with 80 strings) to probe the parameter space of the Constrained MSSM (CMSSM) favoured by current collider and cosmological data. Adopting conservative assumptions about the galactic halo model and the expected experiment performance, we find that IceCube has a probability between 2% and 12% of achievi… ▽ More We study in detail the ability of the nominal configuration of the IceCube neutrino telescope (with 80 strings) to probe the parameter space of the Constrained MSSM (CMSSM) favoured by current collider and cosmological data. Adopting conservative assumptions about the galactic halo model and the expected experiment performance, we find that IceCube has a probability between 2% and 12% of achieving a 5sigma detection of dark matter annihilation in the Sun, depending on the choice of priors for the scalar and gaugino masses and on the astrophysical assumptions. We identify the most important annihilation channels in the CMSSM parameter space favoured by current constraints, and we demonstrate that assuming that the signal is dominated by a single annihilation channel canlead to large systematic errors in the inferred WIMP annihilation cross section. We demonstrate that ~ 66% of the CMSSM parameter space violates the equilibrium condition between capture and annihilation in the center of the Sun. By cross-correlating our predictions with direct detection methods, we conclude that if IceCube does detect a neutrino flux from the Sun at high significance while direct detection experiments do not find a signal above a spin-independent cross section sigma_SI^p larger than 5x10^{-9} pb, the CMSSM will be strongly disfavoured, given standard astrophysical assumptions for the WIMP distribution. This result is robust with respect to a change of priors. We argue that the proposed low-energy DeepCore extension of IceCube will be an ideal instrument to focus on relevant CMSSM areas of parameter space. △ Less

Submitted 27 July, 2009; v1 submitted 2 June, 2009; originally announced June 2009.

Comments: 32 pages, 12 figures. Updated discussion of comparison with direct detection. References added. Main results unchanged. Matches version accepted by JCAP

Journal ref: JCAP 0908:034,2009

Flat Tree-level Inflationary Potentials in Light of CMB and LSS Data

Authors: G. Ballesteros, J. A. Casas, J. R. Espinosa, R. Ruiz de Austri, R. Trotta

Abstract: We use cosmic microwave background and large scale structure data to test a broad and physically well-motivated class of inflationary models: those with flat tree-level potentials (typical in supersymmetry). The non-trivial features of the potential arise from radiative corrections which give a simple logarithmic dependence on the inflaton field, making the models very predictive. We also consid… ▽ More We use cosmic microwave background and large scale structure data to test a broad and physically well-motivated class of inflationary models: those with flat tree-level potentials (typical in supersymmetry). The non-trivial features of the potential arise from radiative corrections which give a simple logarithmic dependence on the inflaton field, making the models very predictive. We also consider a modified scenario with new physics beyond a certain high-energy cut-off showing up as non-renormalizable operators (NRO) in the inflaton field. We find that both kinds of models fit remarkably well CMB and LSS data, with very few free parameters. Besides, a large part of these models naturally predict a reasonable number of e-folds. A robust feature of these scenarios is the smallness of tensor perturbations (r < 10^{-3}). The NRO case can give a sizeable running of the spectral index while achieving a sufficient number of e-folds. We use Bayesian model comparison tools to assess the relative performance of the models. We believe that these scenarios can be considered as a standard physical class of inflationary models, on a similar footing with monomial potentials. △ Less

Submitted 21 November, 2007; originally announced November 2007.

Comments: 42 LaTeX pages, 8 figures

Report number: IFT-UAM/CSIC-07-60, CERN-PH-TH/2007-225

Journal ref: JCAP 0803:018,2008

Determining Reheating Temperature at Colliders with Axino or Gravitino Dark Matter

Authors: Ki-Young Choi, Leszek Roszkowski, Roberto Ruiz de Austri

Abstract: After a period of inflationary expansion, the Universe reheated and reached full thermal equilibrium at the reheating temperature T_R. In this work we point out that, in the context of effective low-energy supersymmetric models, LHC measurements may allow one to determine T_R as a function of the mass of the dark matter particle assumed to be either an axino or a gravitino. An upper bound on t… ▽ More After a period of inflationary expansion, the Universe reheated and reached full thermal equilibrium at the reheating temperature T_R. In this work we point out that, in the context of effective low-energy supersymmetric models, LHC measurements may allow one to determine T_R as a function of the mass of the dark matter particle assumed to be either an axino or a gravitino. An upper bound on their mass may also be derived. △ Less

Submitted 25 March, 2008; v1 submitted 17 October, 2007; originally announced October 2007.

Comments: 19 pages, some improvements, JHEP version

Journal ref: JHEP0804:016,2008

On prospects for dark matter indirect detection in the Constrained MSSM

Authors: Leszek Roszkowski, Roberto Ruiz de Austri, Joe Silk, Roberto Trotta

Abstract: In the framework of the Constrained MSSM we derive the most probable ranges of the diffuse gamma radiation flux from the direction of the Galactic center and of the positron flux from the Galactic halo due to neutralino dark matter annihilation. We find that, for a given halo model, and assuming flat priors, the 68% probability range of the integrated gamma-ray flux spans about one order of magn… ▽ More In the framework of the Constrained MSSM we derive the most probable ranges of the diffuse gamma radiation flux from the direction of the Galactic center and of the positron flux from the Galactic halo due to neutralino dark matter annihilation. We find that, for a given halo model, and assuming flat priors, the 68% probability range of the integrated gamma-ray flux spans about one order of magnitude, while the 95% probability range can be much larger and extend over four orders of magnitude (even exceeding five for a tiny region at small neutralino mass). The detectability of the signal by GLAST depending primarily on the cuspiness of the halo profile. The positron flux, on the other hand, appears to be too small to be detectable by PAMELA, unless the boost factor is at least of order ten and/or the halo profile is extremely cuspy. We also briefly discuss the sensitivity of our results to the choice of priors. △ Less

Submitted 14 November, 2008; v1 submitted 4 July, 2007; originally announced July 2007.

Comments: final PLB version, discussion of prior dependence added, paper based on a new software tool SuperBayeS (available from http://superbayes.org)

Journal ref: Phys.Lett.B671:10-14,2009

Implications for the Constrained MSSM from a new prediction for b to s gamma

Authors: Leszek Roszkowski, Roberto Ruiz de Austri, Roberto Trotta

Abstract: We re-examine the properties of the Constrained MSSM in light of updated constraints, paying particular attention to the impact of the recent substantial shift in the Standard Model prediction for BR(B to X_s gamma). With the help of a Markov Chain Monte Carlo scanning technique, we vary all relevant parameters simultaneously and derive Bayesian posterior probability maps. We find that the case… ▽ More We re-examine the properties of the Constrained MSSM in light of updated constraints, paying particular attention to the impact of the recent substantial shift in the Standard Model prediction for BR(B to X_s gamma). With the help of a Markov Chain Monte Carlo scanning technique, we vary all relevant parameters simultaneously and derive Bayesian posterior probability maps. We find that the case of μ>0 remains favored, and that for μ<0 it is considerably more difficult to find a good global fit to current constraints. In both cases we find a strong preference for a focus point region. This leads to improved prospects for detecting neutralino dark matter in direct searches, while superpartner searches at the LHC become more problematic, especially when μ<0. In contrast, prospects for exploring the whole mass range of the lightest Higgs boson at the Tevatron and the LHC remain very good, which should, along with dark matter searches, allow one to gain access to the otherwise experimentally challenging focus point region. An alternative measure of the mean quality-of-fit which we also employ implies that present data are not yet constraining enough to draw more definite conclusions. We also comment on the dependence of our results on the choice of priors and on some other assumptions. △ Less

Submitted 13 December, 2007; v1 submitted 14 May, 2007; originally announced May 2007.

Comments: JHEP version

Journal ref: JHEP 0707:075,2007

Prospects for direct dark matter detection in the Constrained MSSM

Authors: Roberto Trotta, Roberto Ruiz de Austri, Leszek Roszkowski

Abstract: We outline the WIMP dark matter parameter space in the Constrained MSSM by performing a comprehensive statistical analysis that compares with experimental data predicted superpartner masses and other collider observables as well as a cold dark matter abundance. We include uncertainties arising from theoretical approximations as well as from residual experimental errors on relevant Standard Model… ▽ More We outline the WIMP dark matter parameter space in the Constrained MSSM by performing a comprehensive statistical analysis that compares with experimental data predicted superpartner masses and other collider observables as well as a cold dark matter abundance. We include uncertainties arising from theoretical approximations as well as from residual experimental errors on relevant Standard Model parameters. We present high-probability regions for neutralino dark matter direct detection cross section, and we find that 10^-10 pb < sigma_p^SI < 10^-8 for direct WIMP detection (with details slightly dependent on the assumptions made). We highlight a complementarity between LHC and WIMP dark matter searches in exploring the CMSSM parameter space. We conclude that most of the 95% probability region for the cross section will be explored by future one-tonne detectors, that will therefore cover most of the currently favoured region of parameter space. △ Less

Submitted 5 September, 2006; originally announced September 2006.

Comments: 5 pages, 3 figures. Contribution to the proceedings of the Francesco Melchiorri Memorial conference, Rome, April 2006

Journal ref: NewAstron.Rev.51:316-320,2007

A Markov Chain Monte Carlo Analysis of the CMSSM

Authors: Roberto Ruiz de Austri, Roberto Trotta, Leszek Roszkowski

Abstract: We perform a comprehensive exploration of the Constrained MSSM parameters using a Markov Chain Monte Carlo technique and a Bayesian analysis. We compute superpartner masses and other collider observables as well as a cold dark matter abundance, and compare them with experimental data. We include uncertainties arising from theoretical approximations as well as from residual experimental errors on… ▽ More We perform a comprehensive exploration of the Constrained MSSM parameters using a Markov Chain Monte Carlo technique and a Bayesian analysis. We compute superpartner masses and other collider observables as well as a cold dark matter abundance, and compare them with experimental data. We include uncertainties arising from theoretical approximations as well as from residual experimental errors on relevant SM parameters. We delineate probability distributions of the CMSSM parameters, the collider and cosmological observables as well as a dark matter direct detection cross section. The 68% probability intervals of the CMSSM parameters are: 0.52 TeV < m_{1/2} < 1.26 TeV, m_0 <2.10 TeV, -0.34 TeV < A_0 < 2.41 TeV and 38.5< tan(beta) <54.6. Generally, large fractions of high probability ranges of the superpartner masses will be probed at the LHC. We highlight a complementarity between LHC and WIMP dark matter searches in exploring the CMSSM parameter space. We further expose a number of correlations among the observables, in particular between BR(B_s \to μ^+ μ^-) and BR({\bar B}\to X_sγ) or sigma_p^{SI}. Once SUSY is discovered, this and other correlations may prove helpful in distinguishing the CMSSM from other supersymmetric models. The robustness of our results is investigated in terms of the assumed ranges of CMSSM parameters and the effect of the (g-2)_mu anomaly which shows some tension with the other observables. We find that the results for m_0, and the observables which strongly depend on it, are sensitive to our assumptions, while our conclusions for the other variables are robust. △ Less

Submitted 28 September, 2006; v1 submitted 2 February, 2006; originally announced February 2006.

Comments: 52 pages, 18 PS figures, published version, minor corrections

Journal ref: JHEP 0605:002,2006

Solving the Cosmic Lithium Problems with Gravitino Dark Matter in the CMSSM

Authors: Karsten Jedamzik, Ki-Young Choi, Leszek Roszkowski, Roberto Ruiz de Austri

Abstract: Standard Big Bang nucleosynthesis at baryonic density as inferred by WMAP implies a primordial Li7 abundance factor of two to three larger than that inferred by observations of low--metallicity halo stars. Recent observations of Li6 in halo stars suggest a pre--galactic origin of this isotope, which is exceedingly difficult to explain by putative high redshift cosmic ray populations. We investig… ▽ More Standard Big Bang nucleosynthesis at baryonic density as inferred by WMAP implies a primordial Li7 abundance factor of two to three larger than that inferred by observations of low--metallicity halo stars. Recent observations of Li6 in halo stars suggest a pre--galactic origin of this isotope, which is exceedingly difficult to explain by putative high redshift cosmic ray populations. We investigate if one or both of these lithium problems may be solved by late decaying relic particles in the Constrained Minimal Supersymmetric Standard Model (CMSSM) coupled to gravity. Assuming that the gravitino is the lightest supersymmetric particle (LSP) we find that in large parts of the CMSSM parameter space either of these problems can be potentially solved. In particular, we also find solutions where both lithium problems may be solved simultaneously. These entail the hadronic decay of relic 1 TeV staus into 100 GeV gravitinos at 1000 sec after the Big Bang, as proposed by one of us before. Moreover, the produced gravitinos naturally contribute a large fraction, or all, to the by WMAP required dark matter density. A prediction of this model is the dark matter to be lukewarm with free--streaming lengths of the order of a 4 keV early freezing--out relic particle. Such lukewarm gravitino dark matter may have beneficial effects to the formation of galactic structures and may be constrained by the reionisation history of the Universe. The current version of the paper presents results for both cases (a) when catalytic nuclear reations are included and (b) when they are neglected. △ Less

Submitted 30 October, 2008; v1 submitted 2 December, 2005; originally announced December 2005.

Comments: catalytic effects are fully included in revised version, results are shown with and without catalytic effects, conclusions hardly changed

Journal ref: JCAP 0607 (2006) 007

Gravitino Dark Matter in the CMSSM With Improved Constraints from BBN

Authors: David G. Cerdeno, Ki-Young Choi, Karsten Jedamzik, Leszek Roszkowski, Roberto Ruiz de Austri

Abstract: In the framework of the Constrained MSSM we re--examine the gravitino as the lightest superpartner and a candidate for cold dark matter in the Universe. Unlike in other recent studies, we include both a thermal contribution to its relic population from scatterings in the plasma and a non--thermal one from neutralino or stau decays after freeze--out. Relative to a previous analysis [1] we update,… ▽ More In the framework of the Constrained MSSM we re--examine the gravitino as the lightest superpartner and a candidate for cold dark matter in the Universe. Unlike in other recent studies, we include both a thermal contribution to its relic population from scatterings in the plasma and a non--thermal one from neutralino or stau decays after freeze--out. Relative to a previous analysis [1] we update, extend and considerably improve our treatment of constraints from observed light element abundances on additional energy released during BBN in association with late gravitino production. Assuming the gravitino mass in the GeV to TeV range, and for natural ranges of other supersymmetric parameters, the neutralino region is excluded, while for smaller values of the gravitino mass it becomes allowed again. The gravitino relic abundance is consistent with observational constraints on cold dark matter from BBN and CMB in some well defined domains of the stau region but, in most cases, only due to a dominant contribution of the thermal population. This implies, depending on the gravitino mass, a large enough reheating temperature. If $\mgravitino>1$ GeV then $T_R>10^7$ GeV, if allowed by BBN and other constraints but, for light gravitinos, if $\mgravitino>100$ keV then $T_R>3\times 10^3$ GeV. On the other hand, constraints mostly from BBN imply an upper bound $T_R \lsim {a few}x 10^8\times10^9$ GeV which appears inconsistent with thermal leptogenesis. Finally, most of the preferred stau region corresponds to the physical vacuum being a false vacuum. The scenario can be partially probed at the LHC. △ Less

Submitted 6 July, 2007; v1 submitted 25 September, 2005; originally announced September 2005.

Comments: Version with Erratum. Numerical bug fixed. An upper bound on the reheating temperature strengthened by about an order of magnitude

Journal ref: JCAP 0606:005,2006

Constraints on a mixed inflaton and curvaton scenario for the generation of the curvature perturbation

Authors: G. Lazarides, R. Ruiz de Austri, R. Trotta

Abstract: We consider a supersymmetric grand unified model which naturally solves the strong CP and mu problems via a Peccei-Quinn symmetry and leads to the standard realization of hybrid inflation. We show that the Peccei-Quinn field of this model can act as curvaton. In contrast to the standard curvaton hypothesis, both the inflaton and the curvaton contribute to the total curvature perturbation. The mo… ▽ More We consider a supersymmetric grand unified model which naturally solves the strong CP and mu problems via a Peccei-Quinn symmetry and leads to the standard realization of hybrid inflation. We show that the Peccei-Quinn field of this model can act as curvaton. In contrast to the standard curvaton hypothesis, both the inflaton and the curvaton contribute to the total curvature perturbation. The model predicts an isocurvature perturbation too which has mixed correlation with the adiabatic one. The cold dark matter of the universe is mostly constituted by axions plus a small amount of lightest sparticles. The predictions of the model are confronted with the Wilkinson microwave anisotropy probe and other cosmic microwave background radiation data. We analyze two representative choices of parameters and derive bounds on the curvaton contribution to the adiabatic perturbation. We find that, for the choice which provides the best fitting of the data, the curvaton contribution to the adiabatic amplitude must be smaller than about 67% (at 95% confidence level). The best-fit power spectra are dominated by the adiabatic part of the inflaton contribution. We use Bayesian model comparison to show that this choice of parameters is disfavored with respect to the pure inflaton scale-invariant case with odds of 50 to 1. For the second choice of parameters, the adiabatic mode is dominated by the curvaton, but this choice is strongly disfavored relative to the pure inflaton scale-invariant case (with odds of 10^7 to 1). We conclude that in the present framework the perturbations must be dominated by the adiabatic component from the inflaton. △ Less

Submitted 29 September, 2004; originally announced September 2004.

Comments: 27 pages including 16 figures, uses Revtex

Report number: UT-STPD-3/04

Journal ref: Phys.Rev. D70 (2004) 123527