|
Lopez-Fogliani, D. E., Perez, A. D., & Ruiz de Austri, R. (2021). Dark matter candidates in the NMSSM with RH neutrino superfields. J. Cosmol. Astropart. Phys., 04(4), 067–35pp.
Abstract: R-parity conserving supersymmetric models with right-handed (RH) neutrinos are very appealing since they could naturally explain neutrino physics and also provide a good dark matter (DM) candidate such as the lightest supersymmetric particle (LSP). In this work we consider the next-to-minimal supersymmetric standard model (NMSSM) plus RH neutrino superfields, with effective Majorana masses dynamically generated at the electroweak scale (EW). We perform a scan of the relevant parameter space and study both possible DM candidates: RH sneutrino and neutralino. Especially for the case of RH sneutrino DM we analyse the intimate relation between both candidates to obtain the correct amount of relic density. Besides the well-known resonances, annihilations through scalar quartic couplings and coannihilation mechanisms with all kind of neutralinos, are crucial. Finally, we present the impact of current and future direct and indirect detection experiments on both DM candidates.
|
|
|
Kim, J. S., Lopez-Fogliani, D. E., Perez, A. D., & Ruiz de Austri, R. (2023). Right-handed sneutrino and gravitino multicomponent dark matter in light of neutrino detectors. J. Cosmol. Astropart. Phys., 04(4), 050–32pp.
Abstract: We investigate the possibility that right-handed (RH) sneutrinos and gravitinos can coexist and explain the dark matter (DM) problem. We compare extensions of the minimal supersymmetric standard model (MSSM) and the next-to-MSSM (NMSSM) adding RH neutrinos superfields, with special emphasis on the latter. If the gravitino is the lightest supersymmetric particle (LSP) and the RH sneutrino the next-to-LSP (NLSP), the heavier particle decays to the former plus left-handed (LH) neutrinos through the mixing between the scalar partners of the LH and RH neutrinos. However, the interaction is suppressed by the Planck mass, and if the LH-RH sneutrino mixing parameter is small, << O(10-2), a long-lived RH sneutrino NLSP is possible even surpassing the age of the Universe. As a byproduct, the NLSP to LSP decay produces monochromatic neutrinos in the ballpark of current and planned neutrino telescopes like Super-Kamiokande, IceCube and Antares that we use to set constraints and show prospects of detection. In the NMSSM+RHN, assuming a gluino mass parameter M3 = 3 TeV we found the following lower limits for the gravitino mass m3/2 >= 1-600 GeV and the reheating temperature TR >= 105-3 x 107 GeV, for m nu similar to R similar to 10-800 GeV. If we take M3 = 10 TeV, then the limits on TR are relaxed by one order of magnitude.
|
|
|
Jueid, A., Kip, J., Ruiz de Austri, R., & Skands, P. (2023). Impact of QCD uncertainties on antiproton spectra from dark-matter annihilation. J. Cosmol. Astropart. Phys., 04(4), 068–15pp.
Abstract: Dark-matter particles that annihilate or decay can undergo complex sequences of processes, including strong and electromagnetic radiation, hadronisation, and hadron de-cays, 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 annihila-tion. As a proof-of-principle, we show that for a two-parameter model that depends only on the thermally-averaged annihilation cross section ((o -v)) and the dark-matter mass (Mx), QCD uncertainties can affect the best-fit mass by up to ti 14% (with large uncertainties for large DM masses), depending on the choice of Mx and the annihilation channel (bb over bar or W+W-), and (o -v) by up to ti 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 GitHub [1] repository.
|
|
|
Cabrera, M. E., Casas, J. A., & Ruiz de Austri, R. (2010). MSSM forecast for the LHC. J. High Energy Phys., 05(5), 043–48pp.
Abstract: We perform a forecast of the MSSM with universal soft terms (CMSSM) for the LHC, based on an improved Bayesian analysis. We do not incorporate ad hoc measures of the fine-tuning to penalize unnatural possibilities: such penalization arises from the Bayesian analysis itself when the experimental value of M-Z is considered. This allows to scan the whole parameter space, allowing arbitrarily large soft terms. Still the low-energy region is statistically favoured (even before including dark matter or g-2 constraints). Contrary to other studies, the results are almost unaffected by changing the upper limits taken for the soft terms. The results are also remarkable stable when using flat or logarithmic priors, a fact that arises from the larger statistical weight of the low-energy region in both cases. Then we incorporate all the important experimental constrains to the analysis, obtaining a map of the probability density of the MSSM parameter space, i.e. the forecast of the MSSM. Since not all the experimental information is equally robust, we perform separate analyses depending on the group of observables used. When only the most robust ones are used, the favoured region of the parameter space contains a significant portion outside the LHC reach. This effect gets reinforced if the Higgs mass is not close to its present experimental limit and persits when dark matter constraints are included. Only when the g-2 constraint (based on e(+)e(-) data) is considered, the preferred region (for μ> 0) is well inside the LHC scope. We also perform a Bayesian comparison of the positive- and negative-mu possibilities.
|
|
|
Caron, S., Gomez-Vargas, G. A., Hendriks, L., & Ruiz de Austri, R. (2018). Analyzing gamma rays of the Galactic Center with deep learning. J. Cosmol. Astropart. Phys., 05(5), 058–24pp.
Abstract: We present the application of convolutional neural networks to a particular problem in gamma ray astronomy. Explicitly, we use this method to investigate the origin of an excess emission of GeV gamma rays in the direction of the Galactic Center, reported by several groups by analyzing Fermi-LAT data. Interpretations of this excess include gamma rays created by the annihilation of dark matter particles and gamma rays originating from a collection of unresolved point sources, such as millisecond pulsars. We train and test convolutional neural networks with simulated Fermi-LAT images based on point and diffuse emission models of the Galactic Center tuned to measured gamma ray data. Our new method allows precise measurements of the contribution and properties of an unresolved population of gamma ray point sources in the interstellar diffuse emission model. The current model predicts the fraction of unresolved point sources with an error of up to 10% and this is expected to decrease with future work.
|
|
|
Amoroso, S., Caron, S., Jueid, A., Ruiz de Austri, R., & Skands, P. (2019). Estimating QCD uncertainties in Monte Carlo event generators for gamma-ray dark matter searches. J. Cosmol. Astropart. Phys., 05(5), 007–44pp.
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(*) -> qq-', 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 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.
|
|
|
Balazs, C. et al, Mamuzic, J., & Ruiz de Austri, R. (2021). A comparison of optimisation algorithms for high-dimensional particle and astrophysics applications. J. High Energy Phys., 05(5), 108–46pp.
Abstract: Optimisation problems are ubiquitous in particle and astrophysics, and involve locating the optimum of a complicated function of many parameters that may be computationally expensive to evaluate. We describe a number of global optimisation algorithms that are not yet widely used in particle astrophysics, benchmark them against random sampling and existing techniques, and perform a detailed comparison of their performance on a range of test functions. These include four analytic test functions of varying dimensionality, and a realistic example derived from a recent global fit of weak-scale supersymmetry. Although the best algorithm to use depends on the function being investigated, we are able to present general conclusions about the relative merits of random sampling, Differential Evolution, Particle Swarm Optimisation, the Covariance Matrix Adaptation Evolution Strategy, Bayesian Optimisation, Grey Wolf Optimisation, and the PyGMO Artificial Bee Colony, Gaussian Particle Filter and Adaptive Memory Programming for Global Optimisation algorithms.
|
|
|
Feroz, F., Cranmer, K., Hobson, M., Ruiz de Austri, R., & Trotta, R. (2011). Challenges of profile likelihood evaluation in multi-dimensional SUSY scans. J. High Energy Phys., 06(6), 042–23pp.
Abstract: Statistical inference of the fundamental parameters of supersymmetric theories is a challenging and active endeavor. Several sophisticated algorithms have been employed to this end. While Markov-Chain Monte Carlo (MCMC) and nested sampling techniques are geared towards Bayesian inference, they have also been used to estimate frequentist confidence intervals based on the profile likelihood ratio. We investigate the performance and appropriate configuration of MULTINEST, a nested sampling based algorithm, when used for profile likelihood-based analyses both on toy models and on the parameter space of the Constrained MSSM. We find that while the standard configuration previously used in the literarture is appropriate for an accurate reconstruction of the Bayesian posterior, the profile likelihood is poorly approximated. We identify a more appropriate MULTINEST configuration for profile likelihood analyses, which gives an excellent exploration of the profile likelihood (albeit at a larger computational cost), including the identification of the global maximum likelihood value. We conclude that with the appropriate configuration MULTINEST is a suitable tool for profile likelihood studies, indicating previous claims to the contrary are not well founded.
|
|
|
Casas, J. A., Gomez Vargas, G. A., Moreno, J. M., Quilis, J., & Ruiz de Austri, R. (2018). Extended Higgs-portal dark matter and the Fermi-LAT Galactic Center Excess. J. Cosmol. Astropart. Phys., 06(6), 031–16pp.
Abstract: In the present work, we show that the Galactic Center Excess (GCE) emission, as recently updated by the Fermi-LAT Collaboration, could be explained by a mixture 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, especially from DM direct detection limits. Thus we consider the most economical extension, called ESHP (for extended SHP), which consists solely in the addition of a second (more massive) scalar singlet in the dark sector. The second scalar can be integrated-out, leaving a standard SHP plus a dimension-6 operator. Mainly, 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 -> hh. We demonstrate that, despite its economy, the ESHP model provides an excellent fit to the GCE (with p-value similar to 0.6-0.7) for very reasonable values of the parameters, in particular, ms similar or equal to 130 GeV. This agreement of the DM candidate to the GCE properties does not clash with other observables and keep the S – particle relic density at the accepted value for the DM content in the universe.
|
|
|
Gomez, M. E., Lola, S., Ruiz de Austri, R., & Shafi, Q. (2018). Confronting SUSY GUT With Dark Matter, Sparticle Spectroscopy and Muon (g – 2). Front. Physics, 6, 127–9pp.
Abstract: We explore the implications of LHC and cold dark matter searches for supersymmetric particle mass spectra in two different grand unified models with left-right symmetry, SO(10) and SU(4)(c) x SU(2)(L) x SU(2)(R) (4-2-2). We identify characteristic differences between the two scenarios, which imply distinct correlations between experimental measurements and the particular structure of the GUT group. The gauge structure of 4-2-2 enhances significantly the allowed parameter space as compared to SO(10), giving rise to a variety of coannihilation scenarios compatible with the LHC data, LSP dark matter and the ongoing muon g-2 experiment.
|
|