|
|
Bridges, M., Cranmer, K., Feroz, F., Hobson, M., Ruiz de Austri, R., & Trotta, R. (2011). A coverage study of the CMSSM based on ATLAS sensitivity using fast neural networks techniques. J. High Energy Phys., 03(3), 012–23pp.
Abstract: We assess the coverage properties of confidence and credible intervals on the CMSSM parameter space inferred from a Bayesian posterior and the profile likelihood based on an ATLAS sensitivity study. In order to make those calculations feasible, we introduce a new method based on neural networks to approximate the mapping between CMSSM parameters and weak-scale particle masses. Our method reduces the computational effort needed to sample the CMSSM parameter space by a factor of similar to 10(4) with respect to conventional techniques. We find that both the Bayesian posterior and the profile likelihood intervals can significantly over-cover and identify the origin of this effect to physical boundaries in the parameter space. Finally, we point out that the effects intrinsic to the statistical procedure are conflated with simplifications to the likelihood functions from the experiments themselves.
|
|
|
|
Johannesson, G., Ruiz de Austri, R., Vincent, A. C., Moskalenko, I. V., Orlando, E., Porter, T. A., et al. (2016). Bayesian analysis of cosmic-ray propagation: evidence against homogeneous diffusion. Astrophys. J., 824(1), 16–19pp.
Abstract: We present the results of the most complete scan of the parameter space for cosmic ray (CR) injection and propagation. We perform a Bayesian search of the main GALPROP parameters, using the MultiNest nested sampling algorithm, augmented by the BAMBI neural network machine-learning package. This is the first study to separate out low-mass isotopes (p, (p) over bar and He) from the usual light elements (Be, B, C, N, and O). We find that the propagation parameters that best-fit p, (p) over bar, and He data are significantly different from those that fit light elements, including the B/C and Be-10/Be-9 secondary-to-primary ratios normally used to calibrate propagation parameters. This suggests that each set of species is probing a very different interstellar medium, and that the standard approach of calibrating propagation parameters using B/C can lead to incorrect results. We present posterior distributions and best-fit parameters for propagation of both sets of nuclei, as well as for the injection abundances of elements from H to Si. The input GALDEF files with these new parameters will be included in an upcoming public GALPROP update.
|
|
|
|
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.
|
|
|
|
Pato, M., Baudis, L., Bertone, G., Ruiz de Austri, R., Strigari, L. E., & Trotta, R. (2011). Complementarity of dark matter direct detection targets. Phys. Rev. D, 83(8), 083505–11pp.
Abstract: We investigate the reconstruction capabilities of the 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 similar to 4 for xenon and germanium, compared to the case when astrophysical quantities are fixed. However, the combination of argon, germanium, and xenon data increases the constraining power by a factor of similar to 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.
|
|
|
|
Bertone, G., Cerdeño, D. G., Fornasa, M., Pieri, L., Ruiz de Austri, R., & Trotta, R. (2012). Complementarity of indirect and accelerator dark matter searches. Phys. Rev. D, 85(5), 055014–10pp.
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-parameters supersymmetric model. Combining mock high-luminosity LHC data with presentday null searches for gamma rays from dwarf galaxies with the Fermi Large Area Telescope, we show that current Fermi Large Area Telescope limits already have the capability of ruling out a spurious wino-like solution which 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 the Draco dwarf galaxy with a Cherenkov-Telescope-Array-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.
|
|
|
|
Trotta, R., Johannesson, G., Moskalenko, I. V., Porter, T. A., Ruiz de Austri, R., & Strong, A. W. (2011). Constraints on Cosmic-Ray Propagation Models from a Global Bayesian Analysis. Astrophys. J., 729(2), 106–16pp.
Abstract: Research in many areas of modern physics such as, e. g., indirect searches for dark matter and particle acceleration in supernova remnant shocks rely heavily on studies of cosmic rays (CRs) and associated diffuse emissions (radio, microwave, X-rays, gamma-rays). While very detailed numerical models of CR propagation exist, a quantitative statistical analysis of such models has been so far hampered by the large computational effort that those models require. Although statistical analyses have been carried out before using semi-analytical models (where the computation is much faster), the evaluation of the results obtained from such models is difficult, as they necessarily suffer from many simplifying assumptions. The main objective of this paper is to present a working method for a full Bayesian parameter estimation for a numerical CR propagation model. For this study, we use the GALPROP code, the most advanced of its kind, which uses astrophysical information, and nuclear and particle data as inputs to self-consistently predict CRs, gamma-rays, synchrotron, and other observables. We demonstrate that a full Bayesian analysis is possible using nested sampling and Markov Chain Monte Carlo methods (implemented in the SuperBayeS code) despite the heavy computational demands of a numerical propagation code. The best-fit values of parameters found in this analysis are in agreement with previous, significantly simpler, studies also based on GALPROP.
|
|
|
|
Bertone, G., Cumberbatch, D., Ruiz de Austri, R., & Trotta, R. (2012). Dark Matter searches: the nightmare scenario. J. Cosmol. Astropart. Phys., 01(1), 004–24pp.
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 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 100 fb(-1) of integrated luminosity at 14TeV. We also demonstrate that IceCube (80 strings plus DeepCore) will be able to probe as much as similar or equal to 17% of the currently favoured parameter space after 5 years of observation.
|
|
|
|
Liem, S., Bertone, G., Calore, F., Ruiz de Austri, R., Tait, T. M. P., Trotta, R., et al. (2016). Effective field theory of dark matter: a global analysis. J. High Energy Phys., 09(9), 077–22pp.
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 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.
|
|
|
|
Roszkowski, L., Ruiz de Austri, R., & Trotta, R. (2010). Efficient reconstruction of constrained MSSM parameters from LHC data: A case study. Phys. Rev. D, 82(5), 055003–12pp.
Abstract: We present an efficient method of reconstructing the parameters of the constrained MSSM from assumed future LHC data, applied both on their own right and in combination with the cosmological determination of the relic dark matter abundance. Focusing on the ATLAS SU3 benchmark point, we demonstrate that our simple Gaussian approximation can recover the values of its parameters remarkably well. We examine two popular noninformative priors and obtain very similar results, although when we use an informative, naturalness-motivated prior, we find some sizeable differences. We show that a further strong improvement in reconstructing the SU3 parameters can by achieved by applying additional information about the relic abundance at the level of WMAP accuracy, although the expected data from Planck will have only a very limited additional impact. Further external data may be required to break some remaining degeneracies. We argue that the method presented here is applicable to a wide class of low-energy effective supersymmetric models, as it does not require one to deal with purely experimental issues, e.g., detector performance, and has the additional advantages of computational efficiency. Furthermore, our approach allows one to distinguish the effect of the model's internal structure and of the external data on the final parameters constraints.
|
|
|
|
Bertone, G., Calore, F., Caron, S., Ruiz de Austri, R., Kim, J. S., Trotta, R., et al. (2016). Global analysis of the pMSSM in light of the Fermi GeV excess: prospects for the LHC Run-II and astroparticle experiments. J. Cosmol. Astropart. Phys., 04(4), 037–20pp.
Abstract: We present a new global fit of the 19-dimensional phenomenological Minimal Supersymmetric Standard Model (pMSSM-19) that complies with all the latest experimental results from dark matter indirect, direct and accelerator dark matter searches. We show that the model provides a satisfactory explanation of the excess of gamma rays from the Galactic centre observed by the Fermi Large Area Telescope, assuming that it is produced by the annihilation of neutralinos in the Milky Way halo. We identify two regions that pass all the constraints: the first corresponds to neutralinos with a mass similar to 80 – 100 GeV annihilating into WW with a branching ratio of 95%; the second to heavier neutralinos, with mass similar to 180 – 200 GeV annihilating into (l) over barl with a branching ratio of 87%. We show that neutralinos compatible with the Galactic centre GeV excess will soon be within the reach of LHC run-II – notably through searches for charginos and neutralinos, squarks and light smuons – and of Xenon1T, thanks to its unprecedented sensitivity to spin-dependent cross-section off neutrons.
|
|
