|
|
Bertone, G., Cerdeño, D. G., Fornasa, M., Ruiz de Austri, R., & Trotta, R. (2010). Identification of dark matter particles with LHC and direct detection data. Phys. Rev. D, 82(5), 055008–7pp.
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 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 Omega(DM). We show that a significantly better reconstruction of the DM properties can be obtained with a combined analysis of LHC and direct detection data, by making a simple Ansatz on the weakly interacting massive particles local density rho(0)((chi) over bar1), i.e., by assuming that the local density scales with the cosmological relic abundance, (rho(0)((chi) over bar1)/rho(DM)) = (Omega(0)((chi) over bar1)/Omega(DM)). We demonstrate this method in an explicit example in the context of a 24-parameter supersymmetric model, with a neutralino lightest supersymmetric particle in the stau coannihilation region. Our results show that future ton-scale direct detection experiments will allow to break degeneracies in the supersymmetric parameter space and achieve a significantly better reconstruction of the neutralino composition and its relic density than with LHC data alone.
|
|
|
|
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.
|
|
|
|
Strege, C., Bertone, G., Besjes, G. J., Caron, S., Ruiz de Austri, R., Strubig, A., et al. (2014). Profile likelihood maps of a 15-dimensional MSSM. J. High Energy Phys., 09(9), 081–59pp.
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 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, mu), 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.5TeV 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, similar to 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 x 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.
|
|
|
|
Cabrera, M. E., Casas, J. A., Ruiz de Austri, R., & Trotta, R. (2011). Quantifying the tension between the Higgs mass and (g-2)(mu) in the constrained MSSM. Phys. Rev. D, 84(1), 015006–7pp.
Abstract: Supersymmetry has often been invoked as the new physics that might reconcile the experimental muon magnetic anomaly, a(mu), with the theoretical prediction (basing the computation of the hadronic contribution on e(+)e(-) data). However, in the context of the constrained minimal supersymmetric standard model (CMSSM), the required supersymmetric contributions (which grow with decreasing supersymmetric masses) are in potential tension with a possibly large Higgs mass (which requires large stop masses). In the limit of very large m(h) supersymmetry gets decoupled, and the CMSSM must show the same discrepancy as the standard model with a(mu). But it is much less clear for which size of m(h) does the tension start to be unbearable. In this paper, we quantify this tension with the help of Bayesian techniques. We find that for m(h) >= 125 GeV the maximum level of discrepancy given the current data (similar to 3.2 sigma) is already achieved. Requiring less than 3 sigma discrepancy, implies m(h) less than or similar to 120 GeV. For a larger Higgs mass we should give up either the CMSSM model or the computation of a(mu) based on e(+)e(-); or accept living with such an inconsistency.
|
|
|
|
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.
|
|
|
|
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., Kong, K. C., Ruiz de Austri, R., & Trotta, R. (2011). Global fits of the minimal universal extra dimensions scenario. Phys. Rev. D, 83(3), 036008–15pp.
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 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 similar to 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 fb(-1) of data should be able to probe the current best-fit UED parameters.
|
|
|
|
Roszkowski, L., Ruiz de Austri, R., Trotta, R., Tsai, Y. L. S., & Varley, T. A. (2011). Global fits of the nonuniversal Higgs model. Phys. Rev. D, 83(1), 015014–19pp.
Abstract: We carry out global fits to the nonuniversal Higgs Model (NUHM), applying all relevant present-day constraints. We present global probability maps for the NUHM parameters and observables (including collider signatures, direct, and indirect detection quantities), both in terms of posterior probabilities and in terms of profile likelihood maps. We identify regions of the parameter space where the neutralino dark matter in the model is either binolike, or else higgsinolike with mass close to 1 TeV and a spin-independent scattering cross section similar to 10(-9)-10(-8) pb. We trace the occurrence of the higgsinolike region to be a consequence of a mild focusing effect in the running of one of the Higgs masses, the existence of which in the NUHM we identify in our analysis. Although the usual binolike neutralino is more prominent, higgsinolike dark matter cannot be excluded, however its significance strongly depends on the prior and statistics used to assess it. We note that, despite experimental constraints often favoring different regions of parameter space to the constrained minimal supersymmetric standard model, most observational consequences appear fairly similar, which will make it challenging to distinguish the two models experimentally.
|
|
|
|
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.
|
|
|
|
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.
|
|
