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.
|
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.
|
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.
|
Casas, J. A., Moreno, J. M., Rius, N., Ruiz de Austri, R., & Zaldivar, B. (2011). Fair scans of the seesaw. Consequences for predictions on LFV processes. J. High Energy Phys., 03(3), 034–22pp.
Abstract: We give a straightforward procedure to scan the seesaw parameter-space, using the common “R-parametrization”, in a complete way. This includes a very simple rule to incorporate the perturbativity requirement as a condition for the entries of the R-matrix. As a relevant application, we show that the somewhat propagated belief that BR(mu -> e, gamma) in supersymmetric seesaw models depends strongly on the value of theta(13) is an “optical effect” produced by incomplete scans, and does not hold after a careful analytical and numerical study. When the complete scan is done, BR(mu -> e, gamma) gets very insensitive to theta(13). This holds even if the right-handed neutrino masses are kept constant or under control (as is required for succesful leptogenesis). In most cases the values of BR(mu -> e, gamma) are larger than the experimental upper bound. Including (unflavoured) leptogenesis does not introduce any further dependence on theta(13), although decreases the typical value of BR(mu -> e, gamma).
|