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.
|
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.
|
Strege, C., Bertone, G., Cerdeño, D. G., Fornasa, M., Ruiz de Austri, R., & Trotta, R. (2012). Updated global fits of the cMSSM including the latest LHC SUSY and Higgs searches and XENON100 data. J. Cosmol. Astropart. Phys., 03(3), 030–22pp.
Abstract: We present new global fits of the constrained Minimal Supersymmetric Standard Model (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.
|
Strege, C., Bertone, G., Feroz, F., Fornasa, M., Ruiz de Austri, R., & Trotta, R. (2013). Global fits of the cMSSM and NUHM including the LHC Higgs discovery and new XENON100 constraints. J. Cosmol. Astropart. Phys., 04(4), 013–40pp.
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(-1) integrated luminosity null Supersymmetry searches by ATLAS, the new LHCb measurement of B R ((B) over bar (s) -> mu(+) mu(-)) 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.
|